# Compton scattering from the proton: An analysis using the delta expansion up to N3LO

###### Abstract

We report on a chiral effective field theory calculation of Compton scattering from the proton. Our calculation includes pions, nucleons, and the Delta(1232) as explicit degrees of freedom. It uses the “delta expansion”, and so implements the hierarchy of scales , through a counting . In this expansion the power counting in the vicinity of the Delta peak changes, and resummation of the loop graphs associated with the Delta width is indicated. This is designed to extend the region of PT applicability to momenta .

We have computed the nucleon Compton amplitude in the delta expansion up to NLO——for photon energies . This is the first order at which the proton Compton scattering amplitudes receive contributions from contact operators which encode contributions to the spin-independent polarisabilities from states with energies . We fit the coefficients of these two operators to the experimental proton Compton data that has been taken in the relevant photon-energy domain, and are in a position to extract new results for and .

Compton scattering from the proton: An analysis using the delta expansion up to N3LO

Deepshikha Shukla

George Washington University^{†}^{†}thanks: Present affiliation: The University of North Carolina at Chapel Hill

E-mail: choudhur@gwu.edu

\abstract@cs

## 1 Introduction to Compton scattering and Chiral Dynamics

Chiral dynamics in the baryonic sector is typically thought of as the study of the interactions of pions and nucleons. However the dictates of electromagnetic gauge invariance mean that chiral symmetry also strongly constrains the interactions of both with photons, and so Compton scattering from the nucleon is as fundamental a probe of chiral dynamics as pion-nucleon or nucleon-nucleon scattering. The lowest-order term in the Compton scattering amplitude (the long-wavelength limit) is the Thomson term which is reproduced by PT but which, depending as it does only on the nucleon charge and mass, is independent of chiral dynamics. However, at shorter wavelengths the probing proton starts to be sensitive to the structure of the target. At NLO in heavy baryon chiral perturbation theory (HBPT) the dominant new contribution comes from a single pion loop with photons coupling to the pion or to the N vertex (see Fig. 1), and hence a prediction can be made for these structure effects. This includes, but is not limited to, the numbers known as the polarisabilities of the nucleon; the latter are the first terms in an expansion of the six scattering amplitudes in powers of the photon energy.

The application of chiral dynamics to Compton scattering dates right back to the dawn of baryon PT, and most famously the lowest-order predictions for the electric and magnetic polarisability of the proton in HBPT, fm and fm [1] are in extremely good agreement with experimental determinations. Higher-order calculations [2, 3], and calculations of spin polarisabilities [4, 5, 6, 7, 8] followed. The first studies to systematically compare the full fourth-order predicted cross-section to a compendium of experimental data were published by some of the current authors [9, 10], first with values of and taken from the Particle Data Group, then using them as fit parameters (at fourth order, they are no longer predicted but have contributions from LECs). These studies obtained an excellent fit at low energies, and clearly demonstrated the pion-production cusp. Furthermore the fourth-order results are in markedly better agreement with the data than the third-order ones (even when the freedom to fit the polarisabilities is not exploited). However they also showed that the cross-section completely failed to keep pace with the sharp rise in the experimental data above the cusp, particularly at backward angles. This is not surprising: from around 200 MeV the experimental cross section is dominated by the Delta resonance. While contributions from the Delta enter the Delta-less theory via the LECs, these clearly cannot describe the resonance. Nonetheless, by restricting the fit to appropriate regions of energy and angle, the following results were obtained [10]:

Work has been done on including the Delta in chiral effective theories from the earliest days of baryon PT, motivated partly by its role as the large- degenerate partner of the nucleon, and partly by its experimental relevance. However there is no universal agreement as to the correct power-counting; it is not obvious how to treat the scale . The most commonly adopted solution is to treat them as proportional to the same scale (the “small scale expansion”), a pragmatic treatment which however is not strictly in the spirit of PT as does not vanish in the limit. With this treatment loops come in at the same order as loops. In the particular case of Compton scattering the lowest-order prediction for and is destroyed with this counting. More recently Pascalutsa and Phillips proposed a more sophisticated treatment in which the power counting differs in different energy regions, depending on which scales are enhanced. At low energies they introduce the ratio , with so that . With this “delta-expansion”, the loops enter at one order higher than loops.

Compton scattering cross sections calculated at lowest non-trivial order in the small scale expansion were compared to data by Hildebrandt et al. [11]. Because of the very large contribution of the Delta to and , they had to introduce counterterms for these, which they then fit to data in a similar fashion to Ref. [10]. Interestingly, in spite of the rather different input and data set, the results were very similar:

However the Delta-full results were able to describe data up to around 200 MeV lab energy, unlike the Delta-less results.

Clearly the next step is to include both the Delta and higher-order effects, and that is the aim of the current program. Very recently Lensky and Pascalutsa have also taken steps in this direction, using a relativistic formulation. However they have not reached the order at which they must fit and , and they have so far only compared to a subset of world data [12].

## 2 Comments on the proton Compton database

Modern experiments measure the differential cross section for tagged photons over a range of energies and angles. Between 1956 and 1995, twelve Compton scattering experiments took data at energies below 200 MeV. In a useful paper Baranov et al [13] examined the world data on proton Compton scattering, and demonstrated that the data from the 50s and 60s was compatible with the more modern data from 1974 onwards, and was useful in reducing errors. In addition, the most recent and extensive experiment used the TAPS detector and MAMI tagged photon facility at Mainz [14].

As part of the current project we have evaluated the consistency of the various data sets in as model-independent a fashion as we can. We find that it is crucial to allow the normalizations to float (as done in Ref. [13, 10] but not in [11]). For the Mainz data we increase the statistical error by including a “random systematic” error of 5% in quadrature, as described in Ref. [14]. Some of the very early data cannot be accommodated, nor—as is well known—can the Bonn data of Genzel et al. [15]. With these caveats, we find that the world data is reasonably consistent in the range 0–240 MeV. There is enough noise in the data sets that a of one per degree of freedom is not obtainable though.

Such data has traditionally been analyzed in one of two ways. Early experiments typically took very low energy data (up to around 100 MeV) and fitted to a low-energy expansion of the cross section which included only Born contributions and the electric and magnetic polarisabilities, with these as fit parameters. However at these energies the Born cross section dominates, and the sensitivity to and is limited. A desire to use the full range of data up to the region of the pion-production threshold meant that later authors turned to dispersion relations to obtain a prediction for the variation of the cross section with energy, leaving only certain combinations of the polarisabilities as free parameters to be fit. Olmos de León et al. [14], for instance, fit , and . Unfortunately the dispersion relation approach is not entirely free of assumptions; for instance the asymptotic form of the backward spin-independent amplitude is assumed to be dominated by exchange. There is only one attempt, by Baranov et al. [13] to fit the whole data set (up to 150 MeV) which was unfortunately done before the Mainz data was released.

The use of HBPT with or without the Delta provides an alternative, model-independent approach to the problem, and the history of such attempts is given above.

## 3 An analysis of proton Compton data using the delta-expansion

In figure 1 we show the type of diagrams which contribute at each order.

Most of the parameters which enter are well known. The values of the second order scattering LECs which we use must of course be ones appropriate for a Delta-full theory. The coupling constant (variously or in the literature) is not well known; though it has been fit to resonant electroproduction, the power counting in that regime differs from that appropriate to the low-energy regime, and so it is not clear that the same value should be used. We take the same approach as Hildebrandt et al. [11] and include it as a fit parameter.

There is a problem with the fit when including both loops and higher-order terms. Both individually correct the tendency of the third-order cross section to be too low in the vicinity of the cusp at backward angles. Together they over-correct. It is easy to deduce that the reason is related to the spin polarisabilities: there is exactly the same over-correction in the value of . We therefore promote a counterterm from higher order to allow us to set this to the best estimates of its experimental value [16, 11] or to fit it to data.

There is an obvious tension in deciding how high to take our energy cut-off in fitting data. On the one hand the higher we take it, the more data we can include and the lower the statistical error. On the other hand we know that we are not using the correct power counting for the Delta resonance region and we would not expect a good description of the data there. There are also issues associated with the choice of frame for the calculation which become more acute at higher energies. We have looked at cut-offs up to 240 MeV lab energy.

At this point our fits are preliminary, and we will not quote results. Figure 2 however shows a fit demonstrating that we can obtain a good description of data up to at least 200 MeV. The two curves illustrate the effect of differences in the treatment of the Delta: in one the treatment strictly follows the low-energy power counting, in the other we have included the Delta width and corrections to the vertex. The parameters used are consistent with those obtained by Beane et al. [10] and Hildebrandt et al. [11]. In particular a largish value of is preferred, a result consistent with the recent analysis of Lensky and Pascalutsa [12]

A main goal of this study is to produce reliable single-nucleon amplitudes for incorporation in a code for scattering from the deuteron [17], which will allow the analysis the data that is currently being taken at MAX-lab [18]. We believe we have essentially reached that goal. However the prospects for reducing the error bars on and are not so good. More high-quality lower-energy data (below pion production) is urgently needed.

## Acknowledgments

This research was supported by the US Department of Energy under grants DE-FG02-93ER40756 (DRP) and DE-FG02-95ER-40907(HWG), by the Science and Technology Facilities Council of the UK under grants R105132 (DRP) and R103273 (JMcG), and by the National Science Foundation under career grant PHY-0645498 (HWG and DS). JMcG acknowledges the hospitality of George Washington University and The Institute for Nuclear Theory, University of Washington, where some of the work was carried out.

## References

- [1] V. Bernard, N. Kaiser, J. Kambor and U.-G. Meißner, Nucl. Phys. B 388 (1992) 315.
- [2] V. Bernard, N. Kaiser, A. Schmidt and U.-G. Meißner, Phys. Lett. B 319 (1993) 269.
- [3] V. Bernard, N. Kaiser, U.-G. Meißner and A. Schmidt, Z. Phys. A348 (1993) 317.
- [4] V. Bernard, N. Kaiser and U.-G. Meißner, Int. J. Mod. Phys. E 4 (1995) 193.
- [5] X. Ji, C-W. Kao and J. Osborne, Phys. Rev. D 61 (2000) 074003.
- [6] K. B. V. Kumar, J. A. McGovern and M. C. Birse, hep-ph/9909442.
- [7] G. C. Gellas, T. R. Hemmert and U.-G. Meißner, Phys. Rev. Lett. 85 (2000) 14.
- [8] K. B. V. Kumar, J. A. McGovern, M. C. Birse, Phys. Lett. B479 (2000) 167.
- [9] J.A. McGovern, Phys. Rev. C 63 (2001) 064608.
- [10] S. R. Beane, M. Malheiro, J. A. McGovern, D. R. Phillips and U. van Kolck, Phys. Lett. B567 (2003) 200; Erratum-ibid. B607 (2005) 320; Nucl. Phys. A747 (2005) 311.
- [11] R. P. Hildebrandt, H. W. Griesshammer, T. R. Hemmert and B. Pasquini, Eur. Phys. J. A 20 (2004) 293.
- [12] V. Lensky and V. Pascalutsa, JETP Lett. 89 (2009) 108. See also the talk by V. Pascalutsa in the proceedings of this conference PoS(CD09)095.
- [13] P.S. Baranov, A.I. L’vov, V.A. Petrun’kin, and L.N. Shtarkov, Phys. Part. Nucl. 32 (2001) 376.
- [14] V. Olmos de León et al., Eur. Phys. J. A 10 (2001) 207.
- [15] H. Genzel et al, Z. Phys A 279 (1976) 399.
- [16] D. Drechsel, M. Gorchtein, B. Pasquini and M. Vanderhaeghen, Phys. Rev. C 61 (1999) 015204.
- [17] H. W. Grießhammer, in the proceedings of this conference PoS(CD09)060.
- [18] G. Feldman, talk at INT Program 08-39W ‘‘Soft Photons and Light Nuclei’’, June 2008.