Handling the Cornell potential within the Lagrange-mesh method in momentum space

Titolo: Gestione del potenziale di Cornell nel metodo della mesh di Lagrange nello spazio degli impulsi

1. Il Problema

Le equazioni di Schrödinger indipendenti dal tempo sono fondamentali nella fisica degli adroni (ad esempio, nei modelli a quark costituenti) per descrivere sistemi a due corpi. Sebbene il Metodo della Mesh di Lagrange (LMM) sia stato ampiamente sviluppato e dimostrato estremamente preciso nello spazio delle configurazioni (coordinate spaziali $r$), la sua applicazione nello spazio degli impulsi ($p$) presenta delle limitazioni significative.

In particolare, la metodologia precedente per il LMM nello spazio degli impulsi (sviluppata da Lacroix et al., Ref. [11]) si basa sul calcolo della trasformata di Fourier del potenziale. Questo approccio funziona bene per potenziali a corto raggio (come Gaussiani o Yukawa), ma fallisce numericamente per potenziali a lungo raggio o singolari, come:

• Il potenziale Coulombiano ($1/r$).
• Il potenziale Lineare ($r$), essenziale per il confinamento dei quark.
• Il potenziale di Cornell ($V(r) = -\kappa/r + ar + C$), che combina entrambi ed è cruciale nella fenomenologia degli adroni.

Nell'approccio precedente, il calcolo degli elementi di matrice per $i=j$ (diagonali) porta a divergenze matematiche (funzioni di Legendre di seconda specie che divergono quando l'argomento è 1), rendendo impossibile trattare questi potenziali fondamentali.

2. Metodologia Proposta

Gli autori propongono una metodologia alternativa per calcolare gli elementi di matrice del potenziale nello spazio degli impulsi, ispirandosi a come il LMM nello spazio delle configurazioni gestisce la cinematica semi-relativistica. Invece di trasformare il potenziale nello spazio degli impulsi, il metodo lavora direttamente con l'operatore posizione $r^2$ nello spazio degli impulsi.

La procedura si articola in quattro passaggi chiave:

1. Calcolo dell'operatore $r^2$: Si calcolano gli elementi di matrice dell'operatore $r^2$ sulla base delle funzioni di Lagrange nello spazio degli impulsi. Poiché $r^2$ agisce come un operatore differenziale nello spazio degli impulsi (analogamente a come $p^2$ agisce nello spazio delle configurazioni), si utilizzano le stesse formule derivate per la cinematica nello spazio delle configurazioni.
   $$R^2_{ij} = \langle f_i | r^2 | f_j \rangle \approx \frac{1}{h^2} \left( r^2_{ij} + \frac{l(l+1)}{x_i}\delta_{ij} \right)$$
   dove $x_i$ sono le radici dei polinomi di Laguerre e $h$ è un parametro di scala.
2. Diagonalizzazione: La matrice finita $R^2$ viene diagonalizzata per ottenere una matrice diagonale $D^2$ (contenente gli autovalori di $r^2$) e la matrice di transizione $S$.
   $$R^2 = S D^2 S^{-1}$$
3. Applicazione del Potenziale: Si applica la funzione del potenziale $V$ agli autovalori della matrice diagonale $D^2$. Poiché $D^2$ è diagonale, l'operatore $V(r^2)$ diventa semplicemente $V(\text{autovalori})$.
   $$V_{D^2} = V(D^2)$$
4. Trasformazione di ritorno: Si ritorna alla base originale delle funzioni di Lagrange applicando la trasformazione inversa:
   $$V_{R^2} = S V_{D^2} S^{-1}$$

Questo approccio evita il calcolo diretto della trasformata di Fourier del potenziale e le relative divergenze, permettendo di trattare qualsiasi forma di potenziale, inclusi quelli singolari.

3. Contributi Chiave

• Estensione della gamma di potenziali: Il metodo risolve il problema della divergenza per i potenziali Coulombiani e lineari nello spazio degli impulsi, rendendo possibile l'uso del LMM per il potenziale di Cornell.
• Validazione analitica: Il metodo è stato testato su un sistema Coulombiano non relativistico, confrontando energie e densità di probabilità con soluzioni analitiche esatte.
• Gestione delle densità di probabilità: Sono state sviluppate formule per trasformare le funzioni d'onda dallo spazio degli impulsi a quello delle configurazioni, permettendo il calcolo delle densità di probabilità radiali $R(r)$ e degli osservabili dipendenti dalla posizione.
• Efficienza computazionale: Nonostante le approssimazioni intrinseche del metodo (uso di quadrature di Gauss-Laguerre), la precisione raggiunta è eccezionale con un numero relativamente basso di punti mesh.

4. Risultati

• Sistema Coulombiano:
  • Le energie degli stati legati ($1S, 2S, 1P$) sono state calcolate con alta precisione. Con $N=300$ punti mesh, si ottengono più di 5 cifre significative corrette rispetto ai valori esatti.
  • Le densità di probabilità nello spazio degli impulsi $P(p)$ e nello spazio delle configurazioni $R(r)$ coincidono quasi perfettamente con le soluzioni analitiche per mesh sufficientemente grandi ($N \ge 135$).
  • È stato osservato che per mesh piccoli si possono verificare oscillazioni non fisiche nelle code delle distribuzioni spaziali, mitigabili regolando il parametro di scala $h$.
• Modello di Mesoni (Potenziale di Cornell):
  • Il metodo è stato applicato a un modello semi-relativistico di mesoni leggeri con potenziale di Cornell (parametri presi dalla letteratura [3]).
  • I risultati energetici ottenuti nello spazio degli impulsi sono in ottimo accordo (fino a 4 cifre decimali) con quelli ottenuti in precedenza nello spazio delle configurazioni [10] e con i dati originali [3].
  • Sebbene la convergenza sia leggermente più lenta rispetto al LMM nello spazio delle configurazioni (a causa dell'uso di un parametro $h$ fisso per semplicità), il metodo dimostra di essere robusto e accurato.
• Confronto con il metodo precedente: In un'applicazione a un potenziale Gaussiano, il nuovo metodo mostra una precisione e una convergenza comparabili al metodo precedente [11], confermando che l'approccio alternativo non degrada le prestazioni per i potenziali "facili", ma amplia drasticamente il campo di applicazione.

5. Significato e Impatto

Questo lavoro è significativo per la fisica degli adroni e la cromodinamica quantistica (QCD) fenomenologica perché:

1. Abilita l'uso dello spazio degli impulsi per potenziali realistici: Molti modelli moderni di adroni e descrizioni basate sulla teoria quantistica dei campi sono più naturali nello spazio degli impulsi. Questo metodo rimuove l'ostacolo principale che impediva l'uso del potente LMM in questo dominio per potenziali di confinamento.
2. Flessibilità: Permette ai ricercatori di scegliere la rappresentazione (spazio delle configurazioni o degli impulsi) più adatta al problema specifico senza sacrificare l'accuratezza o la capacità di trattare potenziali complessi come quello di Cornell.
3. Validazione numerica: Dimostra che il LMM, pur non essendo un metodo variazionale rigoroso (a causa dell'uso di quadrature approssimate), mantiene un livello di accuratezza impressionante, rendendolo uno strumento competitivo per lo studio di sistemi a due corpi relativistici e semi-relativistici.

In sintesi, gli autori hanno colmato una lacuna metodologica critica, permettendo l'applicazione efficiente e precisa del metodo della mesh di Lagrange nello spazio degli impulsi a una classe di potenziali fondamentali per la fisica delle particelle.