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

Resumo Técnico: Expansão de Produto de Operadores em Sistemas Bosônicos com Acoplamento Spin-Órbita

Enunciado do Problema
Sistemas bosônicos ultra-resfriados, particularmente aqueles que exibem acoplamento spin-órbita (SOC) induzido por lasers Raman contra-propagantes, apresentam fases quânticas ricas, como estados de momento zero, de onda plana e supersólidos. Embora a Expansão de Produto de Operadores (OPE) tenha sido extensivamente utilizada em sistemas fermiônicos para derivar relações universais (ex: relações de Tan) e compreender a física de curta distância, sua aplicação a sistemas bosônicos com acoplamento spin-órbita e acoplamento coerente permanece subexplorada. O objetivo primário deste trabalho é derivar a OPE para dois operadores de campo locais, $\psi^\dagger_\sigma(\vec{R} - \frac{1}{2}\vec{r})$ e $\psi_{\sigma'}(\vec{R} + \frac{1}{2}\vec{r})$, nestes sistemas bosônicos. Especificamente, os autores visam isolar o termo da "densidade de contato" — a principal contribuição não analítica na distância de separação $r$ — que governa o decaimento assintótico da distribuição de momento para grandes números de onda ($k \to \infty$).

Metodologia
O estudo emprega uma combinação de técnicas de teoria de campo quântico e análise diagramática:

1. Hamiltoniana e Transformação de Base: Os autores partem da ação de um sistema bosônico com acoplamento spin-órbita envolvendo dois estados hiperfinos ($\psi_1, \psi_2$) acoplados pela frequência de Rabi $\Omega$ e transferência de momento $\kappa_0$. Eles realizam uma transformação unitária para diagonalizar o Hamiltoniano de partícula única, introduzindo operadores de "partículas vestidas" ($\alpha_\pm$). No entanto, reconhecendo que as interações tornam-se dependentes do momento nesta base vestida, os cálculos de espalhamento são realizados na base de campo original ($\psi_{1,2}$).
2. Cálculo da Amplitude de Espalhamento: Utilizando a fórmula de redução LSZ, os autores calculam as amplitudes de espalhamento $2 \to 2'$. Eles calculam as funções de Green no espaço de momentos, incluindo os efeitos do acoplamento Raman $\Omega$. A amplitude de espalhamento é derivada resolvendo uma equação matricial envolvendo integrais de um loop ($Z_i(E)$) que contabilizam a autoenergia e as correções de vértice.
3. Derivação da OPE via Valores Esperados: A OPE é derivada avaliando o valor esperado do produto de operador bilocal em estados de duas partículas. Isso é representado diagramaticamente como uma soma de contribuições de:
   • Sem espalhamento (analítico em $\vec{r}$).
   • Espalhamento simples (analítico em $\vec{r}$).
   • Espalhamento duplo (potencialmente não analítico em $\vec{r}$).
     Os autores focam nos diagramas de espalhamento duplo, pois a integração sobre o momento dual a $\vec{r}$ gera termos não analíticos em $\vec{r}$ (especificamente proporcionais a $r$).
4. Correspondência de Coeficientes: Os termos não analíticos derivados dos diagramas de espalhamento duplo são confrontados com os valores esperados de termos locais de densidade de contato de quatro operadores (ex: $(\psi^\dagger \psi)^2$). Isso permite a determinação dos coeficientes de Wilson para a OPE.
5. Análise de Fase: As densidades de contato derivadas são avaliadas através das três fases distintas do sistema SOC, utilizando um ansatz de função de onda coerente variacional.

Principais Contribuições e Resultados

• Derivação de OPEs: O artigo fornece expressões explícitas para a OPE de operadores de campo tanto no regime de acoplamento coerente ($\kappa_0 = 0$) quanto no regime de acoplamento spin-órbita ($\kappa_0 \neq 0$). O termo não analítico principal (ordem $r$) para o operador diagonal $\psi^\dagger_1(\vec{R} - \frac{1}{2}\vec{r})\psi_1(\vec{R} + \frac{1}{2}\vec{r})$ é encontrado como:
  $$\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$$
  Expressões semelhantes são derivadas para operadores fora da diagonal e para o componente $\psi_2$.
• Robustez contra o Acoplamento Spin-Órbita: Um achado significativo é que a forma dos termos de OPE não analíticos principais permanece inalterada na presença do momento de recuo $\kappa_0$. Os autores argumentam que, embora $\kappa_0$ quebre a invariância de Galileu e introduza dependência de momento nas amplitudes de espalhamento, as integrais específicas responsáveis pela dependência não analítica em $r$ produzem os mesmos resultados estruturais que o caso de acoplamento coerente.
• Densidade de Contato em Diferentes Fases: Os autores calculam os valores esperados dos operadores de densidade de contato nas três fases:
  • Fase de Momento Zero: As densidades de contato são constantes e contínuas através da fronteira de fase com a fase de Onda Plana.
  • Fase de Onda Plana: As densidades dependem dos parâmetros de acoplamento e conectam-se continuamente à fase de Momento Zero.
  • Fase de Listras (Stripe): As densidades exibem modulação espacial ($\cos(2\kappa_1 x - \delta)$). Crucialmente, os autores descobrem que, na fronteira de fase entre as fases de Listras e Onda Plana, as densidades de contato são descontínuas. Esta descontinuidade serve como uma assinatura da natureza de primeira ordem desta transição de fase específica.

Significância e Alegações
O artigo afirma preencher uma lacuna na literatura ao derivar explicitamente a OPE e as densidades de contato para sistemas bosônicos com acoplamento spin-órbita e acoplamento coerente, estendendo conceitos previamente estabelecidos principalmente para férmions. O trabalho demonstra que o comportamento universal de grande-$k$ da distribuição de momento nestes sistemas é controlado pelas densidades de contato derivadas.

Os autores enquadram seus resultados de forma modesta como um passo fundamental. Eles explicitamente notam que excluíram a física de três corpos (o efeito Efimov), que é conhecida por introduzir densidades de contato adicionais em sistemas bosônicos, e sugerem isso como uma direção necessária para investigações futuras. Além disso, identificam a OPE de correntes conservadas como uma questão aberta para trabalhos futuros para explorar potencial física universal em correlações de corrente. O estudo estabelece que a natureza das transições de fase (contínua vs. primeira ordem) pode ser observada diretamente através da continuidade ou descontinuidade das densidades de contato.