The dynamically generated $N(1535)$ state in the $\Lambda_c^+ \to p\bar{K}^0 \pi^0$ decay

This paper presents a theoretical analysis of the $\Lambda_c^+ \to p\bar{K}^0 \pi^0$ decay using the chiral unitary approach, demonstrating that the inclusion of a dynamically generated $N(1535)$ resonance alongside other intermediate states successfully reproduces recent Belle experimental data and supports the interpretation of $N(1535)$ as a meson-baryon coupled-channel state.

The Gist

Imagine the subatomic world as a bustling, chaotic dance floor. In this paper, physicists are trying to figure out exactly how the dancers move when a specific partner, the $\Lambda_c^+$ baryon (a heavy particle containing a charm quark), decides to break up and split into three new dancers: a proton, a neutral kaon, and a neutral pion.

Here is the story of what they found, explained simply:

1. The Mystery of the "Ghost" Dancer

When the heavy $\Lambda_c^+$ particle decays, it doesn't just fall apart randomly. It often passes through a "middleman" stage where temporary, short-lived particles (resonances) form before disappearing.

One of the most famous middlemen is called N(1535). For decades, scientists have argued about what this particle actually is.

• The Old Theory: It's just a standard "three-quark" particle, like a normal proton but excited.
• The New Theory (The Paper's Focus): It's a "molecular" structure. Think of it not as a single solid ball, but as a temporary handshake between two other particles (a meson and a baryon) that stick together just long enough to be noticed.

The authors wanted to prove that N(1535) is indeed this "molecular" dancer, created dynamically by the interactions of the other particles on the dance floor.

2. The Detective Work: Comparing Models to Reality

The researchers used a powerful mathematical toolkit called the Chiral Unitary Approach. You can think of this as a high-tech simulation engine. They built a model of the decay process and asked: "If N(1535) is a dynamically generated molecule, will our simulation look like the real data?"

They compared their simulation against real-world data collected by the Belle Collaboration (a giant particle physics experiment in Japan).

The Cast of Characters in their Simulation:
To get the dance right, they couldn't just look at N(1535). They had to include other "guests" who also show up at the party:

• N(1650): Another excited proton-like dancer.
• K*(892) & K*0(1430): Excited versions of the kaon.
• N(1440) & Σ(1750): Later additions to the model to fix small errors.

3. The "Before and After" Photos

The paper presents two versions of their simulation:

• Model A (The First Draft): They included the main dancers (N(1535), N(1650), and the kaon resonances).

  • Result: It looked pretty good! The simulation successfully recreated the big "hump" (peak) in the data around 1535 MeV, which is the signature of the N(1535). This was a huge win for the "molecular" theory.
  • The Glitch: However, there were still some bumps in the data the model couldn't explain, specifically around 1750 MeV and some wiggles in the middle of the energy spectrum.

• Model B (The Final Cut): The scientists realized they missed a few dancers. They added N(1440) and Σ(1750) to the mix.

  • Result: The fit became almost perfect. The simulation now matched the real experimental data (the "Belle data") across the entire board, including the tricky 1750 MeV region.

4. The "Dalitz Plot": The Dance Floor Map

To visualize this, the authors used something called a Dalitz Plot. Imagine a map of the dance floor where every dot represents one single decay event.

• The authors' simulation (the red line in their graphs) traced the exact same pattern as the real experimental dots.
• This map showed that the "N(1535)" isn't just a random bump; it's a distinct, organized structure formed by the way the particles interact.

5. The Big Conclusion

What does this all mean?

1. N(1535) is a "Molecule": The fact that the model works so well when it treats N(1535) as a dynamically generated state (a temporary bond between particles) rather than a rigid three-quark object strongly supports the idea that it is a hadronic molecule. It's like a snowflake: it exists because of how the water molecules stick together, not because it's a single, solid block of ice.
2. The Power of Collaboration: The decay of the $\Lambda_c^+$ is a complex party. You can't understand the music (the data) by listening to just one instrument. You need to hear the whole orchestra (all the intermediate resonances) to get the full picture.
3. Future Steps: While this model explains the current data beautifully, the authors note that future experiments will need to do even more detailed "partial-wave analyses" to confirm the roles of even more exotic dancers (like $\Delta$ and $\Lambda$ resonances) that might be hiding in the background.

In a nutshell: The paper is a successful detective story. By using a sophisticated mathematical model to simulate a particle decay, the authors proved that a mysterious particle called N(1535) is likely a "molecular" structure formed by the dynamic interplay of other particles, and they showed exactly how to account for all the other guests at the party to make the math match reality.

Technical Summary

1. Problem Statement

The internal structure of the $N(1535)$ resonance ($J^P = 1/2^-$) remains a subject of debate in hadron physics. While constituent quark models struggle to explain its mass inversion relative to the $N(1440)$ and its strong coupling to strangeness channels (suggesting a significant $s\bar{s}$ component), alternative theories propose it is a dynamically generated state (a molecular state) arising from coupled-channel meson-baryon interactions (specifically $\pi N$, $\eta N$, $K\Lambda$, and $K\Sigma$).

Recent experimental data from the Belle Collaboration on the decay $\Lambda^+_c \to p K^0_S \pi^0$ revealed distinct peak structures in the invariant mass distributions, particularly near the $p\eta$ threshold, which could be associated with the $N(1535)$ and $N(1650)$ resonances. However, Belle did not perform a full amplitude analysis to disentangle the contributions of these intermediate states. The authors aim to theoretically analyze this decay to:

1. Confirm the dynamical generation mechanism of the $N(1535)$.
2. Quantify the contributions of other intermediate resonances ($N(1650)$, $K^*(892)$, $K^*_0(1430)$, etc.).
3. Provide a comprehensive description of the decay dynamics consistent with Belle's data.

2. Methodology

The authors employ the Chiral Unitary Approach, a framework that combines chiral symmetry with unitarity to describe meson-baryon interactions non-perturbatively.

A. Quark-Level Mechanism

The decay is modeled via an internal emission mechanism:

• The charm quark in the $\Lambda^+_c$ decays weakly ($c \to s + W^+$).
• The resulting $uud$ cluster and the $W^+$-induced $q\bar{q}$ pair hadronize into meson-baryon pairs.
• The primary weak vertex produces a $\bar{K}^0$ and a superposition of meson-baryon states ($\pi N$, $\eta N$, $K\Lambda$).

B. Hadron-Level Dynamics

The final state interactions (FSI) are treated using the Bethe-Salpeter equation in the on-shell factorization scheme:
$$T = [1 - V G]^{-1} V$$

• $V$ (Interaction Kernel): Derived from the leading-order chiral Lagrangian, respecting $SU(3)$ flavor symmetry.
• $G$ (Loop Function): Represents the meson-baryon propagation. The authors use a cutoff regularization ($q_{max} = 1150$ MeV) for most channels but adopt a specific form for the $\pi N$ channel due to its threshold being far from the $N(1535)$ pole.
• Dynamical Generation: The $N(1535)$ is not inserted as a pre-existing pole but emerges naturally as a pole in the scattering amplitude $T$ generated by the coupled-channel interactions.

C. Model Construction

The total decay amplitude ($T_{Model}$) is constructed as a coherent sum of several contributions:

1. Dynamically Generated $N(1535)$: Calculated via the coupled-channel FSI loop.
2. Explicit Resonances:
   • $N(1650)$: S-wave contribution.
   • $K^*(892)$: P-wave contribution in the $\bar{K}^0 \pi^0$ channel.
   • $K^*_0(1430)$: S-wave contribution in the $\bar{K}^0 \pi^0$ channel (incorporating Flatté effects due to its proximity to thresholds).
   • Extended Model (Model B): Includes $N(1440)$ (P-wave) and $\Sigma(1750)$ (S-wave) to address discrepancies in specific mass regions.

The total amplitude includes interference phases ($\phi_i$) between these components. The differential decay widths are calculated and fitted to the Belle invariant mass distributions and Dalitz plots.

3. Key Contributions

• First Full Amplitude Analysis: This work provides the first theoretical amplitude analysis of the $\Lambda^+_c \to p \bar{K}^0 \pi^0$ decay, going beyond simple peak fitting to include interference effects and coupled-channel dynamics.
• Validation of Molecular Nature: By successfully reproducing the $N(1535)$ peak structure solely through coupled-channel interactions without inserting a bare $N(1535)$ seed, the study offers strong support for the molecular interpretation of the $N(1535)$.
• Comprehensive Resonance Catalog: The study systematically evaluates the roles of $N(1535)$, $N(1650)$, $K^*(892)$, $K^*_0(1430)$, $N(1440)$, and $\Sigma(1750)$, clarifying their specific contributions to different invariant mass regions.
• Isospin Symmetry Insight: The paper highlights that the suppression of the $\Delta$ contribution in the neutral channel ($\Lambda^+_c \to p \bar{K}^0 \pi^0$) compared to the charged channel ($\Lambda^+_c \to p K^- \pi^+$) makes the former an ideal laboratory for studying nucleon excited states.

4. Results

The authors fitted two models to the Belle data:

• Model A (Baseline): Includes $N(1535)$, $N(1650)$, $K^*(892)$, and $K^*_0(1430)$.

  • Results: Successfully reproduces the peak at $\sim 1510$ MeV in the $M_{p\pi^0}$ distribution (attributed to $N(1535)$) and the peak at $\sim 1650$ MeV ($N(1650)$). It also captures the $K^*(892)$ peak near 900 MeV in $M_{\bar{K}^0\pi^0}$.
  • Limitations: Discrepancies remained in the $p\bar{K}^0$ distribution around 1750 MeV and the $\bar{K}^0\pi^0$ distribution between 1000–1200 MeV.

• Model B (Extended): Adds $N(1440)$ and $\Sigma(1750)$.

  • Results: The inclusion of $\Sigma(1750)$ perfectly reproduces the structure near 1750 MeV in the $M_{p\bar{K}^0}$ distribution. The $N(1440)$ significantly improves the description of the $\bar{K}^0\pi^0$ spectrum in the 1000–1200 MeV range.
  • Statistical Fit: Model B yields a significantly better fit quality ($\chi^2/\text{d.o.f} = 2.65$) compared to Model A ($\chi^2/\text{d.o.f} = 7.99$).
  • Dalitz Plot: The calculated Dalitz plot for Model B shows excellent agreement with the experimental data, confirming the interference patterns between the various resonances.

5. Significance

• Theoretical Confirmation: The study provides robust theoretical evidence supporting the molecular nature of the $N(1535)$, demonstrating that it can be dynamically generated from meson-baryon interactions in a realistic decay process.
• Experimental Guidance: The analysis clarifies the complex interference patterns in the $\Lambda^+_c \to p K^0_S \pi^0$ decay, offering a template for future partial-wave analyses by experimental collaborations (Belle, LHCb, BESIII).
• Understanding Charmed Baryon Decays: It advances the understanding of non-factorizable processes in charmed baryon decays, showing how final-state interactions can generate or enhance specific resonance signals.
• Isospin Symmetry: The work reinforces the utility of comparing charged and neutral decay modes to isolate specific isospin contributions (e.g., suppressing $\Delta$ states to focus on nucleon excitations).

In conclusion, the paper successfully bridges the gap between recent experimental observations and theoretical models, confirming that the $N(1535)$ observed in $\Lambda^+_c$ decays is consistent with a dynamically generated state arising from coupled-channel meson-baryon interactions.