On operator product expansion in the spin-orbit coupled bosonic system

Sintesi Tecnica: Espansione del Prodotto di Operatori in Sistemi Bosonici con Accoppiamento Spin-Orbita

Definizione del Problema
I sistemi di bosoni ultra-freddi, in particolare quelli che esibiscono un accoppiamento spin-orbita (SOC) indotto da laser Raman contro-propaganti, mostrano fasi quantistiche ricche come stati a momento zero, a onda piana e supersolidi. Sebbene l'Espansione del Prodotto di Operatori (OPE) sia stata ampiamente utilizzata nei sistemi fermionici per derivare relazioni universali (ad esempio, le relazioni di Tan) e comprendere la fisica a breve distanza, la sua applicazione a sistemi bosonici con accoppiamento spin-orbita e accoppiamento coerente rimane poco esplorata. L'obiettivo primario di questo lavoro è derivare l'OPE per due operatori di campo locali, $\psi^\dagger_\sigma(\vec{R} - \frac{1}{2}\vec{r})$ e $\psi_{\sigma'}(\vec{R} + \frac{1}{2}\vec{r})$, in questi sistemi bosonici. Nello specifico, gli autori mirano a isolare il termine di "densità di contatto" — il principale contributo non analitico nella distanza di separazione $r$ — che governa il decadimento asintotico della distribuzione di momento per grandi numeri d'onda ($k \to \infty$).

Metodologia
Lo studio impiega una combinazione di tecniche di teoria quantistica dei campi e analisi diagrammatica:

1. Hamiltoniana e Trasformazione della Base: Gli autori partono dall'azione di un sistema bosonico con accoppiamento spin-orbita che coinvolge due stati iperfini ($\psi_1, \psi_2$) accoppiati dalla frequenza Rabi $\Omega$ e dal trasferimento di momento $\kappa_0$. Eseguono una trasformazione unitaria per diagonalizzare l'Hamiltoniana a singola particella, introducendo operatori di particelle "vestite" ($\alpha_\pm$). Tuttavia, riconoscendo che le interazioni diventano dipendenti dal momento in questa base vestita, i calcoli di scattering vengono esegtti nella base originale dei campi ($\psi_{1,2}$).
2. Calcolo dell'Ampiezza di Scattering: Utilizzando la formula di riduzione LSZ, gli autori calcolano le ampiezze di scattering $2 \to 2'$. Calcolano le funzioni di Green nello spazio dei momenti, includendo gli effetti della frequenza di accoppiamento Raman $\Omega$. L'ampiezza di scattering è derivata risolvendo un'equazione matriciale che coinvolge integrali a un loop ($Z_i(E)$) che tengono conto dell'autenergia e delle correzioni al vertice.
3. Derivazione dell'OPE tramite Valori di Aspettazione: L'OPE viene derivata valutando il valore di aspettazione del prodotto di operatori bilocali in stati a due particelle. Questo è rappresentato diagrammaticamente come una somma di contributi da:
   • Nessuno scattering (analitico in $\vec{r}$).
   • Singolo scattering (analitico in $\vec{r}$).
   • Doppio scattering (potenzialmente non analitico in $\vec{r}$).
     Gli autori si concentrano sui diagrammi di doppio scattering, poiché l'integrazione sul momento duale a $\vec{r}$ genera termini non analitici in $\vec{r}$ (specificamente proporzionali a $r$).
4. Corrispondenza dei Coefficienti: I termini non analitici derivati dai diagrammi di doppio scattering vengono confrontati con i valori di aspettazione dei termini locali di densità di contatto a quattro operatori (ad esempio, $(\psi^\dagger \psi)^2$). Ciò consente la determinazione dei coefficienti di Wilson per l'OPE.
5. Analisi delle Fasi: Le densità di contatto derivate vengono valutate attraverso le tre fasi distinte del sistema SOC: la fase a Onda Piana, la fase a Momento Zero e la fase a Strisce (Stripe), utilizzando un ansatz di funzione d'onda coerente variazionale.

Contributi Chiave e Risultati

• Derivazione delle OPE: Il documento fornisce espressioni esplicite per l'OPE degli operatori di campo sia nei regimi a accoppiamento coerente ($\kappa_0 = 0$) che a accoppiamento spin-orbita ($\kappa_0 \neq 0$). Il principale termine non analitico (ordine $r$) per l'operatore diagonale $\psi^\dagger_1(\vec{R} - \frac{1}{2}\vec{r})\psi_1(\vec{R} + \frac{1}{2}\vec{r})$ è trovato come:
  $$\psi^\dagger_1 \psi_1(\vec{R}) + \vec{r} \cdot \psi^\dagger_1 \overset{\leftrightarrow}{\nabla} \psi_1(\vec{R}) - \frac{r}{8\pi} \left[ g^2(\psi^\dagger_1 \psi_1)^2 + \lambda^2 \psi^\dagger_1 \psi_1 \psi^\dagger_2 \psi_2 \right](\vec{R}) + \dots$$
  Simili espressioni sono derivate per gli operatori off-diagonal e per la componente $\psi_2$.
• Robustezza rispetto all'Accoppiamento Spin-Orbita: Un risultato significativo è che la forma dei termini OPE non analitici principali rimane invariata in presenza del momento di rinculo $\kappa_0$. Gli autori sostengono che, sebbene $\kappa_0$ rompa l'invarianza di Galilei e introduca una dipendenza dal momento nelle ampiezze di scattering, gli integrali specifici responsabili della dipendenza non analitica in $r$ producono risultati strutturali identici al caso a accoppiamento coerente.
• Densità di Contatto nelle Diverse Fasi: Gli autori calcolano i valori di aspettazione degli operatori di densità di contatto nelle tre fasi:
  • Fase a Momento Zero: Le densità di contatto sono costanti e continue attraverso il confine di fase con la fase a Onda Piana.
  • Fase a Onda Piana: Le densità dipendono dai parametri di accoppiamento e si connettono continuamente alla fase a Momento Zero.
  • Fase a Strisce (Stripe): Le densità esibiscono una modulazione spaziale ($\cos(2\kappa_1 x - \delta)$). Fondamentalmente, gli autori riscontrano che al confine di fase tra la fase a Strisce e quella a Onda Piana, le densità di contatto sono discontinue. Questa discontinuità funge da firma della natura del primo ordine di questa specifica transizione di fase.

Significatività e Rivendicazioni
Il documento sostiene di colmare una lacuna nella letteratura derivando esplicitamente l'OPE e le densità di contatto per sistemi bosonici con accoppiamento spin-orbita e accoppiamento coerente, estendendo concetti precedentemente stabiliti principalmente per i fermioni. Il lavoro dimostra che il comportamento universale del grande-$k$ della distribuzione di momento in questi sistemi bosonici è controllato dalle densità di contatto derivate.

Gli autori inquadrano modestamente i loro risultati come un passo fondamentale. Essi notano esplicitamente di aver escluso la fisica a tre corpi (l'effetto Efimov), che è nota per introdurre ulteriori densità di contatto nei sistemi bosonici, e suggeriscono questo come una direzione necessaria per investigazioni future. Inoltre, identificano l'OPE delle correnti conservate come una questione aperta per lavori futuri per esplorare la potenziale fisica universale nelle correlazioni di corrente. Lo studio stabilisce che la natura delle transizioni di fase (continua vs primo ordine) può essere osservata direttamente attraverso la continuità o la discontinuità delle densità di contatto.