Reduced density matrix approach to one-dimensional ultracold bosonic systems

This paper presents a variational method using the two-boson reduced density matrix to accurately calculate the ground-state energies, densities, and correlation functions of one-dimensional harmonically trapped bosons across a wide range of particle numbers ($N=2$ to $10^4$) and interaction strengths, effectively bridging the gap between few-body and mean-field regimes.

The Gist

The "Universal Translator" for Quantum Crowds: A Simple Guide

Imagine you are trying to describe a crowd of people.

If you are looking at two people dancing in a room, you can describe every single detail: where their hands are, how they move, and exactly how they interact. This is easy, but it’s "microscopic" detail.

If you are looking at a million people in a stadium, you don't care about individual hand movements. Instead, you look at the "big picture": the density of the crowd, the waves moving through the stands, and the general flow of the crowd. This is "macroscopic" detail.

The Problem in Physics:
For a long time, physicists have had two different sets of tools. One tool is like a high-powered microscope—it’s amazing for studying two or three particles, but if you try to use it on a million, the computer explodes because there is too much data. The other tool is like a satellite photo—it’s great for seeing the whole crowd, but it’s "blurry" and misses the subtle, important interactions between individuals.

In the world of ultracold bosons (tiny particles cooled to near absolute zero), scientists have struggled to find a single tool that works for both a tiny group of two particles and a massive "cloud" of thousands.

---

The Solution: The "Buddy System" (The 2-RDM Method)

The researchers in this paper used a clever mathematical shortcut called the Two-Particle Reduced Density Matrix (2-RDM).

Instead of trying to track every single person in the stadium (which is impossible), or just looking at the blurry satellite photo (which is inaccurate), they decided to focus on pairs.

The Analogy: The "Buddy System"
Imagine you want to understand the behavior of a massive music festival without tracking every individual. Instead of a list of 100,000 names, you simply keep a massive ledger of how any two people interact.

• Do they bump into each other?
• Do they stay apart?
• Do they dance in sync?

By understanding the "rules of the pair," you can mathematically reconstruct the behavior of the entire crowd. This paper proves that if you know how the "buddies" behave, you can accurately predict the energy, the density, and the "vibe" (the correlations) of the whole system, whether it’s a tiny duo or a massive army.

---

What did they actually find?

The researchers tested this "Buddy System" across a massive range, from 2 particles to 10,000 particles. Here is what they discovered:

1. It’s a Universal Tool: It didn't matter if the particles were barely touching or pushing against each other with extreme force; the math held up. It bridged the gap between the "micro" and the "macro."
2. The "Traffic Jam" Effect: They observed a phenomenon called "fermionisation." Imagine a crowd of people who usually walk through each other like ghosts. Suddenly, you turn up the "interaction strength," and it’s as if everyone suddenly becomes a solid object. They can no longer pass one another, creating a "quantum traffic jam." The 2-RDM method captured this transition perfectly.
3. Smooth Transitions: They showed that as you add more and more particles, the system doesn't suddenly "break" or change its fundamental nature in a weird way. It evolves smoothly, like adding more water to a pool—the physics stays consistent.

Why does this matter?

This is like finding a single mathematical language that can describe both a conversation between two people and the roar of a stadium.

By mastering this "Buddy System" approach, scientists can now study complex quantum materials and gases more efficiently. It gives them a way to simulate massive, real-world quantum systems on computers that would otherwise be too slow to handle the task. It’s a bridge between the tiny world of individual atoms and the massive world of quantum matter.

Technical Summary

Technical Summary: Reduced Density Matrix Approach to One-Dimensional Ultracold Bosonic Systems

1. Problem Statement

The study of ultracold bosonic systems in one dimension (1D) presents a significant theoretical challenge due to the lack of a single, unified methodology capable of describing the entire parameter space. Existing approaches are generally bifurcated into two categories:

• Size-extensive methods: (e.g., direct Hamiltonian diagonalization, Bethe ansatz, or exact solutions for $N=2$) which are highly accurate for few-body systems but scale poorly or are restricted to integrable, translationally invariant models.
• Non-size-extensive methods: (e.g., Mean-field theories like the 1D Gross-Pitaevskii Equation (GPE) or the Non-Polynomial Schrödinger Equation (NPSE)) which are efficient for many-body systems but fail in the "few-body" regime or when correlations are strong (e.g., the Tonks-Girardeau limit).

There is a distinct gap in theoretical tools to describe the crossover regime—the transition from few-body to many-body physics—across varying interaction strengths (from weakly interacting BECs to strongly interacting "fermionised" gases).

2. Methodology

The authors employ a variational Two-Particle Reduced Density Matrix (2-RDM) approach. Instead of attempting to solve for the complex $N$-particle wavefunction $\Psi(x_1, \dots, x_N)$, the method focuses on the 2-RDM, which contains all necessary information to calculate the energy of a system with pairwise interactions.

Key technical components include:

• Variational Scheme: The ground-state energy is expressed as a linear functional of the 1-RDM and 2-RDM.
• N-Representability Constraints: To ensure the trial 2-RDM corresponds to a physical $N$-boson state, the authors enforce specific semidefinite constraints:
  • D-condition: Ensures the 1-RDM and 2-RDM are positive semidefinite.
  • G-condition: Constrains the distribution of one hole and one boson to be non-negative.
  • Note: The authors demonstrate that for bosonic systems, the more complex $Q$, $T_1$, and $T_2$ conditions (often required for fermions) are not strictly necessary to capture accurate structural properties.
• Computational Implementation: The problem is solved using Semidefinite Programming (SDP) via the SDPA and SDPNAL+ algorithms, utilizing a basis of 20 harmonic oscillator eigenfunctions.
• System Model: $N$ bosons in a harmonic trap interacting via a Dirac $\delta$-function contact interaction.

3. Key Contributions

• Demonstration of Universality: The paper proves that the 2-RDM method is "size-extensive" in a practical sense, successfully bridging the gap between $N=2$ and $N=10^4$ particles.
• Validation of the Crossover Regime: The authors provide evidence that the method maintains continuity and physical accuracy in the intermediate regime ($N \sim 100$) where mean-field and few-body methods typically fail.
• Efficiency via Minimal Constraints: The study shows that the relatively simple D and G conditions are sufficient to extract accurate ground-state energies, densities, and correlation functions for bosonic systems.

4. Results

• Ground-State Energies: The 2-RDM results show excellent agreement with the analytic Busch solution for $N=2$ and align with the highly accurate 1D NPSE for large $N$ and strong interactions ($\beta$). It successfully captures the energy plateau in the Tonks-Girardeau (TG) limit.
• Densities: For large $N$ ($10^3$), the 2-RDM accurately reproduces the density profiles of the 1D GPE (weak interaction) and the 1D NPSE (strong interaction), where the standard GPE fails.
• Correlations: The method accurately captures "fermionisation." In the strongly interacting limit, the pair correlation function $\sigma(z,0)$ is suppressed at the origin ($z \to 0$), a hallmark of the Tonks-Girardeau gas.
• Continuity: A plot of the correlation function at the origin $\sigma(0,0)$ vs. interaction strength $\beta$ shows a smooth, continuous evolution across all $N$ values, confirming the method's robustness in the crossover regime.

5. Significance

This work establishes the 2-RDM method as a powerful, versatile tool for quantum gas research. Its ability to handle both few-body and many-body limits within a single framework makes it uniquely suited for studying the transition between different quantum phases. The authors suggest this methodology could be extended to study more complex phenomena, such as solitons, quantum droplets, and providing corrections to the 3D GPE by accounting for quantum fluctuations.