$J/\psi$ Photoproduction from Threshold to HERA: Leading-Twist Convolution, Small-$x$ Pathology, and Eikonal Unitarization

Die Arbeit zeigt, dass eine eikonale Unitarisierung der führenden-Twist-Faltung moderner, bei kleinen $x$ singulärer PDFs die $J/\psi$-Photoproduktion sowohl im HERA-Bereich als auch nahe der Schwelle konsistent beschreibt, wobei der dispersive Realteil und OPE-Subtraktionskonstanten die Schwellendynamik dominieren.

Das Wesentliche

Die Geschichte vom „J/ψ-Ball" und dem unsichtbaren „Gluon-Wetter"

Stellen Sie sich vor, Sie werfen einen kleinen, extrem schweren und dichten Ball (das ist das J/ψ-Meson, ein Teilchen aus einem schweren Quark und seinem Antiquark) gegen eine große, weiche Wand (das Proton oder Nukleon).

Physiker wollen genau verstehen: Was passiert, wenn dieser Ball die Wand trifft? Wie stark prallt er ab? Und was sagt uns das über die Wand selbst?

In diesem Papier untersucht der Autor Arkadiy Syamtomov genau diesen Vorgang, aber er nutzt dabei moderne Werkzeuge, um zu sehen, ob seine Berechnungen mit der Realität übereinstimmen. Er stößt dabei auf zwei große Probleme und findet eine clevere Lösung.

1. Das Problem mit den „Wetterkarten" (Die PDFs)

Um zu berechnen, wie die Wand reagiert, braucht man eine Art „Wetterkarte" für das Innere des Protons. Diese Karte zeigt, wie viele unsichtbare Teilchen namens Gluonen (die Kleber der Materie) sich wo befinden.

In den 1990er Jahren hatten die Physiker eine sehr einfache, glatte Wetterkarte. Damals passte ihre Berechnung gut zu den Daten.
Heute haben wir moderne, hochauflösende Wetterkarten (die sogenannten NNLO-PDFs wie ABMP16, MSHT20 usw.). Diese Karten zeigen etwas Neues: Bei sehr kleinen Entfernungen (oder sehr hohen Energien) gibt es dort eine riesige Explosion an Gluonen. Es ist, als würde das Wetter an einer bestimmten Stelle extrem stürmisch werden.

2. Der erste Fehler: Der „Zauberstab", der verrückt spielt

Der Autor hat versucht, die alten Berechnungsmethoden mit diesen neuen, stürmischen Karten zu nutzen.

• Die Metapher: Stellen Sie sich vor, Sie versuchen, die Form eines Kuchens zu erraten, indem Sie nur die Gesamtmenge an Mehl und Zucker abwiegen (das nennt man „Momente"). Bei der alten, glatten Karte funktionierte das gut.
• Das Problem: Bei der neuen Karte mit dem „Gluon-Sturm" am Rand (bei kleinem $x$) funktioniert dieser Trick nicht mehr. Die Rechnung sagt plötzlich: „Der Kuchen muss am Anfang extrem steil ansteigen und dann sofort wieder abfallen!"
• Das Ergebnis: Die Berechnung sagt voraus, dass der Ball fast sofort nach dem Aufprall explodiert und dann verschwindet. Das ist physikalisch unsinnig. Die modernen Karten haben die alte Rechenmethode „geblendet".

3. Der zweite Fehler: Der „Überschuss" bei hohen Geschwindigkeiten

Also probierte der Autor einen anderen Weg: Er ließ den Ball direkt durch die Wetterkarte fliegen (die „Direkte Faltung").

• Das Ergebnis am Start: Super! Bei niedrigen Geschwindigkeiten (nahe der Schwelle, wo die Experimente wie GlueX gemessen haben) passte die Rechnung perfekt. Der Ball prallte genau so ab, wie die Messungen zeigten.
• Das Problem bei hoher Geschwindigkeit: Als der Ball aber sehr schnell wurde (wie in den alten Experimenten am HERA-Beschleuniger), sagte die Rechnung: „Der Ball prallt 7- bis 12-mal stärker ab als in der Realität!"
• Die Metapher: Es ist, als würde Ihr Wettermodell sagen: „Bei Sturm wird der Ball so schnell fliegen, dass er durch die Wand bricht!" Aber in der Realität prallt er nur ab. Die Rechnung ist zu optimistisch.

4. Die Lösung: Der „Sicherheitsgurt" (Eikonal-Unitarisierung)

Warum ist die Rechnung bei hohen Geschwindigkeiten zu stark? Weil sie annimmt, dass die Gluonen einfach nur addiert werden. Aber in der Quantenwelt gibt es eine Regel: Wenn zu viele Gluonen aufeinandertreffen, fangen sie an, sich gegenseitig zu blockieren. Sie werden „gesättigt".

Der Autor schlägt vor, einen Sicherheitsgurt einzubauen.

• Die Metapher: Stellen Sie sich vor, Sie fahren mit dem Auto. Bei niedriger Geschwindigkeit ist alles okay. Aber wenn Sie zu schnell werden, greift automatisch ein Bremssystem ein, damit Sie nicht gegen eine Wand fliegen.
• Die Anwendung: Der Autor fügt eine Formel hinzu, die diesen „Bremseffekt" (Unitarisierung) simuliert. Wenn die Geschwindigkeit (die Energie $W$) zu hoch wird, dämpft diese Formel die vorherige, zu starke Vorhersage.
• Das Ergebnis: Plötzlich passt die Rechnung perfekt zu allen Daten – sowohl bei langsamen als auch bei schnellen Bällen. Der „Sicherheitsgurt" hat den Überschuss genau richtig weggebremst.

5. Das Geheimnis am Startpunkt

Ein besonders interessanter Fund am Anfang der Geschichte:
Nahe dem Startpunkt (der „Schwelle"), wo der Ball gerade erst die Wand berührt, ist das, was man sieht, fast nur das Gegenteil von dem, was man erwartet.

• Normalerweise denkt man, die Stärke des Aufpralls kommt von der Menge der Gluonen (dem Imaginärteil).
• Aber hier zeigt sich: Das Wichtigste ist eine Art „Trägheit" oder „Gedächtnis" der Wand (der Realteil), das durch die OPE-Methode (Operator-Produkt-Entwicklung) berechnet wird.
• Vergleich: Es ist, als würde der Ball nicht durch den Aufprall selbst abprallen, sondern weil die Wand selbst eine gewisse „Steifigkeit" hat, die man aus der Struktur der Materie ableiten kann. Diese Steifigkeit ist der Schlüssel, um die ersten Experimente zu verstehen.

Fazit für den Alltag

Diese Arbeit sagt uns im Grunde:

1. Alte Methoden scheitern an neuen Daten: Wenn wir moderne, detaillierte Karten des Universums nutzen, funktionieren alte Rechen-Tricks nicht mehr (sie werden verrückt).
2. Einfache Modelle sind zu stark: Wenn wir einfach nur die Gluonen zusammenzählen, denken wir, die Welt sei gewaltiger, als sie ist.
3. Die Natur bremst sich selbst: Bei extremen Bedingungen (hohe Energie) gibt es einen natürlichen „Sicherheitsmechanismus" (Sättigung), der verhindert, dass die Kräfte ins Unendliche wachsen.
4. Der Anfang ist anders: Ganz am Anfang des Prozesses wird das Ergebnis nicht durch den Aufprall bestimmt, sondern durch eine fundamentale Eigenschaft der Materie selbst.

Der Autor hat also gezeigt, wie man die alten Theorien mit den neuen Daten in Einklang bringt, indem er einen cleveren „Sicherheitsgurt" für die Formeln einführt. Das hilft uns, die Bausteine der Materie besser zu verstehen.

Technische Zusammenfassung

Titel: J/ψ-Photoproduktion von der Schwelle bis HERA: Leading-Twist-Faltung, kleine-x-Pathologien und eikonale Unitarisierung

Autor: Arkadiy I. Syamtomov (Bogolyubov-Institut für Theoretische Physik, Ukraine)

---

1. Problemstellung und Motivation

Die Arbeit untersucht die Photoproduktion von $J/\psi$-Mesonen an Nukleonen ($\gamma N \to J/\psi N$) über einen weiten Energiebereich, von der kinematischen Schwelle ($W \approx 4$ GeV) bis zu den hohen Energien des HERA-Beschleunigers ($W \approx 300$ GeV).

Das zentrale Ziel ist die Überprüfung der Gültigkeit des Leading-Twist-Formalismus (basierend auf der Operatorproduktentwicklung, OPE) unter Verwendung moderner NNLO-Gluon-Partonverteilungsfunktionen (PDFs) (ABMP16, MSHT20, CT18, NNPDF4.0).

• Hintergrund: Die $J/\psi$-Photoproduktion dient als Sonde für die gluonische Struktur des Nukleons. Da das $J/\psi$ klein und stark gebunden ist, kann die Wechselwirkung im Rahmen der Störungstheorie (OPE) behandelt werden.
• Das Dilemma: Frühere Analysen (1999) nutzten skalierende PDFs ohne starke Singularitäten bei kleinen $x$. Moderne PDFs weisen jedoch ein starkes Wachstum bei kleinen $x$ auf ($g(x) \sim x^{\alpha_g}$ mit $\alpha_g < 0$). Es ist unklar, wie sich diese moderne Eingabe auf die Rekonstruktion des Wirkungsquerschnitts auswirkt, insbesondere in Kombination mit Zielmasse-Korrekturen (Target-Mass Corrections, TMC) und Dispersionsrelationen.

2. Methodik

Der Autor kombiniert mehrere theoretische Werkzeuge:

1. Vektor-Meson-Dominanz (VMD): Die Photoproduktion wird auf die Vorwärtsstreuung $J/\psi$-Nukleon zurückgeführt. Der Wirkungsquerschnitt hängt vom Imaginärteil (optischer Theorem) und dem Realteil (Dispersionsrelation) der Amplitude ab.
2. Operatorproduktentwicklung (OPE) und Summenregeln:
   • Die Amplitude wird als Summe von Mellin-Momenten der Gluonverteilung $A_n(Q) = \int dx \, x^{n-2} g(x, Q)$ dargestellt.
   • Zielmasse-Korrekturen (TMC): Diese werden explizit berücksichtigt, um nicht-parametrisch kleine Terme der Ordnung $m_N^2/\epsilon_0^2$ zu erfassen. Sie führen zu einem gewichteten Moment $\tilde{A}_n$.
3. Zwei Rekonstruktionsansätze:
   • Momentenbasierte Rekonstruktion: Ein Ansatz mit zwei Parametern ($C, a$) für den Wirkungsquerschnitt $\sigma \sim (s-s_{th})^a / s^{a+1}$, der durch Anpassung an die Summenregeln für $n=4$ und $n=6$ bestimmt wird.
   • Direkte Faltung (Direct Convolution): Eine partonische Faltung (erste Iteration), die den Wirkungsquerschnitt direkt aus dem Integral über die Gluonverteilung berechnet, wobei die kinematische Grenze $x > \epsilon_0/\lambda$ die Integration bei der Schwelle automatisch einschränkt.
4. Dispersionsrelation: Der Realteil der Amplitude wird aus dem Imaginärteil über eine einmal subtrahierte Dispersionsrelation berechnet. Die Subtraktionskonstante $M_{\psi N}(0)$ wird durch die OPE bei Null-Energie bestimmt.
5. Eikonale Unitarisierung: Um Diskrepanzen bei hohen Energien zu beheben, wird eine eikonale Screening-Funktion eingeführt, die Mehrfachstreuungseffekte (Saturation) bei kleinen $x$ modelliert.

3. Wichtige Ergebnisse und Pathologien

Die Analyse identifiziert zwei komplementäre "Pathologien" beim Vergleich von Momenten-basierten und Faltungs-Ansätzen mit modernen PDFs:

A. Die Momenten-basierte Pathologie (Schwellenbereich)

• Phänomen: Bei Verwendung moderner PDFs mit kleinem-$x$-Singularität führt die momentenbasierte Rekonstruktion zu einem unphysikalischen Schwellenexponenten $a \simeq 17\text{--}20$.
• Ursache: Die kleine-$x$-Singularität verzerrt die Hierarchie der Mellin-Momente. Das Verhältnis $\tilde{A}_6 / \tilde{A}_4$ wird um einen Faktor $\simeq 1.7$ im Vergleich zu skalierenden PDFs reduziert. Da der Exponent $a$ durch dieses Verhältnis festgelegt wird, muss er drastisch ansteigen, um die Summenregeln zu erfüllen.
• Folge: Der resultierende Wirkungsquerschnitt steigt extrem steil an der Schwelle an und fällt bei hohen Energien zu schnell ab ($\sim W^{-4}$), was die HERA-Daten um Größenordnungen unterschätzt, wenn er an die Schwellendaten normiert wird.

B. Die Faltungs-Pathologie (Hohe Energien)

• Phänomen: Der direkte Faltungsansatz vermeidet die Schwellen-Pathologie und beschreibt die Daten von GlueX, Cornell und Fermilab (Schwellenbereich) hervorragend für alle vier modernen PDF-Sets.
• Problem: Bei hohen Energien ($W \gtrsim 90$ GeV) überschätzt die Faltung die HERA-Messungen um einen Faktor von 7 bis 12.
• Ursache: Dies ist eine inhärente Eigenschaft der Leading-Twist-Faltung mit PDFs, die bei kleinem $x$ singulär wachsen ($g(x) \sim x^{\alpha_g}$). Der kinematische untere Integrationsgrenze $x_{min} \to 0$ bei hohen Energien erlaubt es dem starken kleinen-$x$-Wachstum, direkt in den Wirkungsquerschnitt einzugehen.
• Vergleich: Die alte Analyse von 1999 mit einer nicht-singulären skalierenden PDF zeigte keine solche Diskrepanz, da dort das kleine-$x$-Wachstum fehlte.

C. Dominanz des Realteils an der Schwelle

• Nahe der kinematischen Schwelle ($W \approx 4.1$ GeV) dominiert der Realteil der Amplitude den Wirkungsquerschnitt, da der Imaginärteil (der mit der inelastischen Streuung verknüpft ist) gegen Null geht.
• Der Realteil wird durch die OPE-Subtraktionskonstante $M_{\psi N}(0) \simeq 36\text{--}39 \, \text{GeV}^{-2}$ verankert und durch das Dispersionsintegral weiter verstärkt (Gesamtwert $\approx 150 \, \text{GeV}^{-2}$).
• Das Verhältnis $\rho = \text{Re}M / \text{Im}M$ ist an der Schwelle sehr groß ($\approx 30\text{--}40$) und fällt mit steigender Energie ab.

4. Lösung: Eikonale Unitarisierung

Um die Diskrepanz bei hohen Energien zu lösen, ohne die erfolgreiche Beschreibung des Schwellenbereichs zu zerstören, wird eine minimale eikonale Unitarisierung vorgeschlagen:
$$\text{Im } M_{\text{unit}} = M_{\text{sat}}(W) \left[ 1 - \exp\left( -\frac{\text{Im } M_{\text{conv}}}{M_{\text{sat}}(W)} \right) \right]$$

• Mechanismus: Dies modelliert die Absättigung (Saturation) der gluonischen Dichte bei kleinen $x$, wenn der Wirkungsquerschnitt die geometrische Größe des Dipols erreicht.
• Parameter: Die Sättigungsamplitude $M_{\text{sat}}(W)$ wächst mit der Energie ($W^\delta$). Die Parameter werden an die HERA-Daten angepasst.
• Ergebnis:
  • Der Schwellenbereich bleibt unverändert, da $\text{Im } M_{\text{conv}} \to 0$ und der Screening-Faktor gegen 1 geht.
  • Bei hohen Energien ($W > 30$ GeV) wird der Leading-Twist-Wirkungsquerschnitt um den Faktor 7–12 unterdrückt.
  • Die angepasste Kurve beschreibt den gesamten $W$-Bereich (Schwelle bis HERA) mit einem $\chi^2/N_{\text{data}} \approx 3/21$.

5. Bedeutung und Schlussfolgerungen

1. Validierung der OPE: Der Leading-Twist-Ansatz ist im Schwellenbereich gültig, aber die Rekonstruktion über Mellin-Momente ist bei modernen, kleinen-$x$-singulären PDFs problematisch. Der direkte Faltungsansatz ist hier überlegen.
2. Physik jenseits von Leading-Twist: Die Diskrepanz bei hohen Energien ist kein Fehler in der PDF-Fit, sondern ein Hinweis darauf, dass Leading-Twist-Effekte bei kleinen $x$ durch Unitarisierung (Saturation) unterdrückt werden müssen. Es handelt sich nicht um einen fehlenden Austauschmechanismus (wie einen zusätzlichen Pomeron-Term), sondern um eine Unterdrückung des Leading-Twist-Beitrags.
3. Schwellenphysik: Die Arbeit bestätigt, dass die $J/\psi$-Photoproduktion nahe der Schwelle fast ausschließlich durch den dispersive Realteil der Amplitude getrieben wird, der durch die OPE-Subtraktionskonstante festgelegt ist. Dies bietet eine Verbindung zu den gravitativen Formfaktoren des Nukleons.
4. Zukunftsausblick: Die Arbeit liefert eine robuste Basis für zukünftige Studien, einschließlich der Anwendung auf Bottomonium (bessere Kurzstreckenkontrolle) und Kernziele (Medium-Effekte), sowie die quantitative Untersuchung der Verbindung zwischen $M_{\psi N}(0)$ und dem Protonenradius.

Zusammenfassend demonstriert das Paper, wie moderne QCD-Daten (NNLO PDFs) in Kombination mit präzisen theoretischen Rahmenwerken (OPE, TMC, Dispersion) neue Einsichten in die Struktur der Nukleonen-Gluonenverteilung liefern, gleichzeitig aber die Grenzen perturbativer Ansätze bei hohen Energien aufzeigen, die durch nicht-perturbative Unitarisierungseffekte korrigiert werden müssen.