Emptiness formation in the Lieb-Liniger gas: hydrodynamic instantons and a conjectured rate function

기술 요약: Lieb–Liniger 가스에서의 공백 형성: 유체역학적 인스턴톤과 추정된 속도 함수

문제 정의
본 논문은 반발하는 1차원 Lieb–Liniger 보스 가스의 바닥 상태에서 나타나는 공백 형성 확률(Emptiness Formation Probability, EFP)을 조사한다. EFP는 $P(R)$로 표기되며, 이는 길이가 $2R$인 공간 구간 내에 입자가 존재하지 않을 확률을 의미한다. 거시적 구간의 경우, 주 점근적 거동은 속도 함수 $f(\gamma_0)$에 의해 지배된다. 여기서 $-\log P(R) \sim (\rho_0 R)^2 f(\gamma_0)$이며, $\rho_0$는 평균 밀도이고 $\gamma_0 = c/\rho_0$는 무차원 상호작용 강도이다. 자유 페르미온 및 특정 극한(Tonks–Girardeau 및 약한 결합)에 대해서는 정확한 결과가 존재하지만, 임의의 반발 결합에 대한 일반적인 해석적 표현 $f(\gamma_0)$를 찾는 것은 오랫동안 미해결 과제로 남아 있었다. Korepin, Its, Waldron(KIW)의 선행 연구는 정확한 이중 필드 프레드홀름 행렬식(dual-field Fredholm-determinant) 표현을 제공했으나, 최근의 유체역학적 연구들은 그들이 제안한 거대 거리 점근식의 추측이 약한 결합 극한에서 파라미터적으로 부정확할 수 있음을 시사했다.

방법론
저자들은 가적분 연산자 이론과 유체역학적 변분법을 결합한 이중 접근 방식을 채택한다:

1. 가적분 연산자 분석: 저자들은 정확한 KIW 이중 필드 프레드홀름 행렬식 표현으로부터 출발하여 거대 거리 점근성을 분석한다. 이들은 기존 KIW 분석에서의 비균일성(non-uniformity)을 식별하였다: 고정된 이중 필드에 대한 행렬식 전개에서 하위 차수의 항들이, 관련 이중 필드가 $O(R)$로 스케일링되기 때문에 이중 필드 진공 평균을 수행한 후에는 $O(R^2)$의 주 차수 항이 된다. 이를 균일한 이중 스케일링(uniform double-scaling) 문제로 취급함으로써(필드 $\phi \to x\Phi$로 재스케일링), 저자들은 결합된 사델-포인트(saddle-point) 조건을 도출한다. 이는 매개변수가 없는 선형 적분 방정식의 해로 표현되는 새로운 속도 함수에 대한 추측으로 이어진다.

2. 유체역학적 변분 접근법: EFP는 유클리드 유체역학 작용량(Euclidean hydrodynamic action)에 대한 변분 문제로 정식화된다. 저자들은 $\tau=0$에서 빈 구간이라는 경계 조건을 만족하도록 이 작용량을 수치적으로 최소화한다. '공백 인스턴톤(emptiness instanton)'이라 불리는 이 최소화 요소는 최적의 희귀 요동 경로를 기술한다. 수치적 구현에는 빈 영역 경계의 아스트로이드(astroid) 형태를 해상하기 위한 경계 적합 좌표계가 사용되었으며, 광범위한 결합 강도에 걸쳐 밀도 및 속도장을 구하기 위해 Adam 및 L-BFGS와 같은 고급 최적화 알고리즘이 사용되었다.

주요 기여

• 속도 함수에 대한 새로운 추측: 저자들은 $g(t)$를 결정하는 속도 함수 $f(\gamma_0)$에 대한 특정 적분 방정식을 제안한다:
  $$\int_{-1}^{1} g(s) \log\left(1 + \frac{r^2}{(t-s)^2}\right) ds = 4t$$
  여기서 $r = c/q$는 Lieb 방정식에 의해 결정된다. 속도 함수는 다음과 같이 주어진다: $f(\gamma_0) = \frac{1}{2}(\gamma_0/r)^2 \int_{-1}^{1} t g(t) dt$. 이 공식은 매개변수가 없으며 Lieb–Liniger 모델의 이체 산란 데이터(two-body scattering data)를 직접 인코딩한다.
• KIW 불일치의 해결: 저자들은 기존 KIW의 추측이 이중 필드의 2차 기여를 간과함으로써 약한 결합 거동을 제대로 포착하지 못함을 보여준다. 새로운 추측은 결합된 사델-포인트 분석을 통해 필요한 2차 항을 포함함으로써 이를 교정한다.
• 수치적 검증: 저자들은 결합 강도 $\gamma_0$에 대해 4개 이상의 차수에 걸친 광범위한 결합 범위에서 Lieb–Liniger 유체역학 작용량을 수치적으로 최소화하는 과정을 제시한다.

결과

• 극한과의 일치: 제안된 속도 함수는 약한 결합($\gamma_0 \to 0$) 및 강한 결합($\gamma_0 \to \infty$) 극한 모두에서 알려진 1차 결과들을 정확하게 재현한다. 또한, 이 모델은 이전에 유체역학적 섭동 이론을 통해 독립적으로 도출된 양쪽 극한에서의 첫 번째 비자명한 하위 차수 보정항들도 정확히 재현한다.
• 수치적 비교: 수치적으로 최소화된 유체역학적 속도 함수 $f_{\text{num}}(\gamma_0)$는 전체 상호작용 교차 구간에서 추측된 곡선 $f_{\text{conj}}(\gamma_0)$를 밀접하게 따른다. 체계적인 상대적 차이는 약 3%–5%이며, 수치 결과가 추측보다 약간 높게 나타난다.
• 수렴 분석: 격자 세분화(grid refinement, $200\times200$에서 $300\times300$으로 증가) 및 확장된 최적화 사이클을 포함한 수렴 테스트 결과, 격자 해상도가 개선됨에 따라 불일치가 감소함을 보여준다. 이는 잔여 3%–5%의 차이가 근본적인 추측의 결함이라기보다는 유한한 격자 해상도 및 불완전한 최적화 수렴에 기인함을 시사한다.
• 인스턴톤 기하학: 수치적 최소화 요소들은 자유 페르미온 및 폴리트로픽 유체(polytropic fluids)의 결과와 일치하는 아스트로이드 형태의 진공 영역을 가진 공백 인스턴톤을 드러낸다. 경계 형상은 약한 결합과 강한 결합 극한 사이를 부드럽게 진화한다.

의의 및 주장
본 논문은 모든 상호작용 강도에 걸쳐 반발하는 Lieb–Liniger 가스의 주한 leading EFP 지수를 설명하는 통합된, 매개변수가 없는 설명을 제공한다고 주장한다. 그 의의는 다음과 같다:

1. 형식적 유도: 정확한 프레드홀름 행렬식 표현으로부터 추측된 적분 방정식으로 이어지는 형식적 유도 경로를 제공하며, 엄밀한 증명이 필요한 특정 점근 단계(uniform double-scaling)를 식별한다.
2. 경험적 견고성: 극한 사례에서의 독립적인 해석적 검증과 고정밀 수치 유체역학을 통해 추측을 뒷받침한다.
3. 미시적 연결: 적분 방정식의 로그 커널이 Lieb–Liniger 산란 커널의 적분임을 강조하며, 이는 거시적인 대편차 속도(large-deviation rate)에 미시적 산란 데이터가 직접 인코딩되어 있음을 시사한다.

저자들은 이 추측의 상태에 대해 겸손한 태도를 유지하며, 수치적 및 점근적 증거에 의해 강력하게 뒷받침되기는 하지만, 균일한 이중 스케일링 점근성에 대한 완전하고 엄밀한 수학적 유도는 여전히 열려 있는 문제임을 밝히고 있다. 또한, 하위 차수 항(exponent $R^2$에 대한 선형 및 로그 보정항)의 구조 역시 미해결 과제로 남아 있다고 언급한다.