Pole positions and residues from pion photoproduction using the Laurent+Pietarinen expansion method

Abstract: We have applied a new approach to determine the pole positions and residues from pion photoproduction multipoles. The method is based on a Laurent expansion of the partial wave T-matrices, with a Pietarinen series representing the regular part of energy-dependent and single-energy photoproduction solutions. The method has been applied to multipole fits generated by the MAID and GWU/SAID groups. We show that the number and properties of poles extracted from photoproduction data correspond very well to results from $π$N elastic data and values cited by Particle Data Group (PDG). The photoproduction residues provide new information for the electromagnetic current at the pole position, which are independent of background parameterizations, as opposed to the Breit-Wigner representation. Finally, we present the photo-decay amplitudes from the current MAID and SAID solutions at the pole, for all four-star nucleon resonances below W=2 GeV.

• The paper introduces the Laurent+Pietarinen expansion as a systematic method to separate pole contributions from background in pion photoproduction data.
• It demonstrates high precision in extracting resonance parameters, achieving sub-MeV accuracy for the Δ(1232) multipole and robust photo-decay amplitudes.
• The study highlights ambiguities in distinguishing genuine resonances from complex branch effects, underscoring the need for multi-channel analyses.

Pole Extraction in Pion Photoproduction via the Laurent+Pietarinen Expansion

Motivation and Context

The characterization of baryon resonances in hadronic and electromagnetic reactions demands precise determination of pole parameters—positions and residues—of the scattering matrix $T(W)$ in the complex energy plane. Traditional experimental analyses have relied on Breit-Wigner (BW) fits with ad hoc backgrounds, but such representations are parametrization-dependent and do not provide unique definitions of resonance properties. Recent consensus in hadron spectroscopy, reflected, for example, in PDG recommendations, is that pole parameters, not BW quantities, faithfully encode resonance information connecting phenomenology with QCD.

While analytic continuation of partial-wave amplitudes to the complex plane is possible in elaborate coupled-channel models, robust and efficient single-channel methods applicable to both model-generated and experimental data remain under active development. This paper introduces and applies the Laurent+Pietarinen (L+P) expansion—a controlled analytic expansion which systematically and efficiently separates pole and regular background contributions in single-channel processes. The method is applied to pion photoproduction multipoles from state-of-the-art MAID and GWU/SAID analyses, enabling the direct extraction of pole positions, complex residues, and subsequently photo-decay amplitudes for all four-star nucleon resonances under $W=2$ GeV (1404.1544).

Analytical Framework: Laurent+Pietarinen Method

The Laurent+Pietarinen expansion is a hybrid analytic technique based on the Mittag-Leffler/Laurent decomposition of $T(W)$ about its simple pole singularities,

$$T(W) = \sum_{i=1}^k \frac{a_{-1}^{(i)}}{W - W_i} + B^L(W),$$

where $a_{-1}^{(i)}$ and $W_i$ are complex residues and pole positions, and $B^L(W)$ is a regular piece (analytic in the domain of interest). The novelty lies in the modeling of $B^L(W)$: instead of polynomial backgrounds of arbitrary or unphysical structure, it is parametrized by rapidly convergent Pietarinen series constructed from conformal mapping around branch points,

$$F(W) = \sum_{n=0}^N c_n X(W)^n, \quad X(W) = \frac{\alpha - \sqrt{x_P - W}}{\alpha + \sqrt{x_P - W}}$$

with $x_P$ the branch point of each channel/cut (physical or unphysical), and $\alpha$, $c_n$ free real parameters. Multiple cuts—elastic and prominent inelastic—are each represented by their own Pietarinen series with individual mapping parameters, allowing for a systematic and data-driven background description.

This construction analytically encodes known physical thresholds and their associated cuts in the $T$-matrix—embedding correct analytic properties at all physical singularities—and, crucially, factors the pole contributions cleanly from the non-singular background. Minimization of model dependence is achieved by selecting the number of background expansions (series and order) only as necessitated by the fit quality.

Fitting Strategy, Error Analysis, and Physical Interpretation

The method is applied both to energy-dependent (ED) and single-energy (SE) partial-wave amplitudes, generated by MAID and SAID, for all major multipoles of pion photoproduction. For each partial wave, optimal numbers and types of Pietarinen series are chosen to accommodate the analytic structure: a left-hand cut (effectively representing all negative-energy processes, e.g., nucleon exchange, $t$/$u$-channel processes), the physical $\pi N$ elastic threshold, and principal inelastic thresholds (e.g., $\eta N$, $\pi\pi N$). The status of each threshold—real or complex—is also tested, with the possibility to represent quasi-two-body channels with complex branch points when intermediate resonances contribute in three-body final states.

Fitting parameters include all pole positions and residues, Pietarinen expansion coefficients, and possibly the locations of the branch points themselves. Optimization is performed via $\chi^2$ minimization (SE) or a discrepancy parameter $D_{dp}$ for model-generated data with no experimental uncertainties.

Error analysis is twofold: statistical errors are determined by fit uncertainties, while systematic errors are evaluated by varying the modeling of branch points (fixing to physical values or treating as free parameters) and the expansion order, probing model dependence. Significant tests are performed to diagnose ambiguities in distinguishing genuine resonance poles from complex branch point effects resulting from intermediate channel coupling—e.g., the indistinguishability of genuine $N(1710)$ pole contributions from effective $\rho N$ branch cuts in the $P_{11}$ wave.

Numerical Results: Pole Spectrum and Comparisons

The L+P analysis extracts the complete set of pole parameters for all dominant nucleon and $\Delta$ resonances up to $W=2$ GeV with high accuracy, using both MAID and SAID ED/SE datasets. The extracted pole positions show excellent agreement (within uncertainties) between different datasets and with values from $\pi N$ elastic scattering. Particularly, the poles of $\Delta(1232)$ are reproduced in the $P_{33}$ multipole with sub-MeV precision in mass and few MeV in width. Consistent information is obtained for $N(1520)$ and other four-star resonances.

Strong discrepancies arise only for less-established states (e.g., $N(1710)$, $N(2000)$), where the fits demonstrate the possible interchangeability between introducing a new pole and allowing a complex branch point to absorb strength—a critical ambiguity given only single-channel data. For example, the $N(1710)$ resonance in the $P_{11}$ multipole is either manifest as a physical pole (with real branch points; needed to describe the amplitude) or, if a complex $\rho N$ branch point is introduced, its distinct pole disappears from the best-fit model. This dichotomy underscores the intrinsic limitation of single-channel analyses in resolving resonance dynamics when significant channel coupling is present.

Residues of the multipole amplitudes are robustly extracted, enabling computation of photo-decay amplitudes $A_{1/2}$, $A_{3/2}$ at pole positions. For the best-established resonances, values are consistent between MAID and SAID with uncertainties at the $<10\%$ level. Direct comparison with other recent coupled-channel extractions (BnGa, Jülich, ANL-Osaka) is provided, showing general consistency for dominant amplitudes, though significant phase differences remain for certain subdominant contributions.

Theoretical and Practical Implications

The L+P method supplies a disciplined, parsimonious, and analytically controlled tool for pole parameter extraction in photoproduction, bridging the gap between phenomenological fits and formal scattering theory. Unlike the fixed-form backgrounds of phenomenological BW fits, the Pietarinen expansion responds adaptively to the analytic structure as dictated by the underlying physics and the actual data, minimizing ad hoc contributions.

Notable technical claims of the work include:

• The number and nature of resonances extracted from single-pion photoproduction with L+P expansion reproduce the established spectrum from $\pi N$ elastic analyses, confirming the method's reliability for leading states.
• The residues extracted at the pole are independent of background parameterization, unlike BW residues, providing direct access to electromagnetic couplings cleanly separated from non-resonant effects.
• Ambiguities remain for less-established resonances in the presence of significant multi-channel coupling, reflected as alternative fits of similar quality with either additional poles or substituted complex branch points.

Practical implications include the capacity to standardize resonance properties for forthcoming PDG listings and the provision of residues for global coupled-channel analyses and reaction modeling, facilitating further connection to QCD-inspired calculations (e.g., Lattice QCD, dynamical coupled-channel models).

Open Questions and Prospects

The inability of single-channel partial-wave analyses to incontrovertibly distinguish between higher-lying resonance poles and complex branch points highlights the necessity for multidimensional analysis, including channels with explicit three-body final states and polarization observables. Improved experimental coverage, including double-polarization observables and complete experiments, is vital to resolving these ambiguities.

The L+P method is readily extensible to other production reactions (e.g., kaon photoproduction, Compton scattering), providing an analytic and systematic means to extract resonance pole information from disparate datasets for global analyses.

Conclusion

The application of the Laurent+Pietarinen expansion to pion photoproduction yields a coherent, controlled extraction of nucleon and $\Delta$ resonance pole positions and electromagnetic residues from MAID and SAID multipole amplitudes. The method achieves high fidelity for leading four-star states and exposes the model sensitivity and physical ambiguities underlying less-certain resonances, especially when only single-channel inputs are available. As PDG and the broader hadronic physics community continue to prioritize pole parameters for resonance listings, the L+P approach supplies an essential tool for this transition, while simultaneously indicating domains where future experimental and theoretical work is necessary—most notably, in disentangling genuine resonant and multichannel-coupling effects.

---

Reference: "Pole positions and residues from pion photoproduction using the Laurent+Pietarinen expansion method" (1404.1544).