Lieb-Liniger model

The quantum many-body problem—understanding how vast numbers of interacting particles behave collectively—stands as one of the most formidable challenges in physics. Most such systems are too complex to solve exactly, forcing physicists to rely on approximations. However, a select few "integrable" models can be solved exactly, offering a pristine window into the intricate dance of quantum mechanics. The Lieb-Liniger model, describing bosons confined to a single dimension, is a cornerstone of this exclusive club. It addresses the fundamental question of how particle interactions and quantum statistics intertwine in a low-dimensional world, a question that has moved from a theoretical curiosity to an experimental reality.

This article provides a comprehensive exploration of this remarkable model. We will first delve into its core theoretical foundations in the chapter on ​​Principles and Mechanisms​​, dissecting the nature of contact interactions and uncovering the genius of the Bethe ansatz solution. We will see how this leads to profound phenomena like fermionization, where repulsive bosons magically mimic fermions. Following this, in the chapter on ​​Applications and Interdisciplinary Connections​​, we will journey beyond the theory to witness its power in action. We will explore how the model provides the definitive language for cutting-edge experiments with ultracold atoms, serves as a theoretical laboratory for the frontier of non-equilibrium physics, and reveals astonishing connections to seemingly unrelated fields like random growth and polymer physics.

Now that we have been introduced to the curious world of one-dimensional bosons, let's roll up our sleeves and explore the machinery that makes it tick. The beauty of the Lieb-Liniger model is that it all boils down to a single, simple idea: how two particles interact. From this one seed, a whole forest of complex, beautiful, and often surprising many-body physics grows.

Imagine a line of tiny quantum particles. What happens when they feel each other's presence? In the simplified world of the Lieb-Liniger model, this interaction is as local as it gets: a "contact" potential, described by the Dirac delta function, $g \delta(x_i - x_j)$. You can think of it as an infinitely sharp spike of potential energy that particles experience only when they are at the exact same position. The nature of this interaction depends entirely on the sign of the coupling constant, $g$.

Let's first consider the case of ​​attraction​​, where $g$ is negative ($g = -c$, with $c > 0$). Attraction, as you might guess, pulls particles together. If you have just two bosons, this attraction can be strong enough to bind them into a stable pair, a sort of quantum "diatomic molecule." This bound state is a new entity, with a lower energy than the two particles would have if they were far apart and at rest. The energy required to break them apart is called the ​​binding energy​​, $E_B$. By solving the simple two-body Schrödinger equation, one finds this binding energy has a beautifully compact form:

This result is quite intuitive: the stronger the attraction (larger $c$), the more tightly bound the pair is, and the more energy it takes to separate them.

But now, let's flip the sign. Let's make the interaction ​​repulsive​​ ($g > 0$). This is where the story truly takes a fascinating turn. Instead of pulling together, the particles now push each other away. They don't want to be in the same place at the same time. How does quantum mechanics handle this game of mutual avoidance? The answer lies not in binding, but in scattering.

When two classical billiard balls collide, they bounce off each other. In the quantum world, particles are waves, so when they meet, they "scatter." Their wavefunctions are altered by the encounter. For the delta-function repulsion, this alteration takes the form of a ​​scattering phase shift​​.

Think of two particles with momenta $k_1$ and $k_2$. Their combined wavefunction can be thought of as a superposition of a wave where particle 1 has momentum $k_1$ and particle 2 has $k_2$, and another wave where they've swapped momenta. The interaction forces a specific relationship between these two parts of the wavefunction right at the point of collision, $x_1 = x_2$. The solution to this reveals that the interaction's sole effect is to introduce a momentum-dependent phase shift, $\theta(k)$, into the wavefunction. For the Lieb-Liniger model, this crucial phase shift is given by:

where $k = k_1 - k_2$ is the relative momentum between the two particles. This equation is the heart of the entire model. It tells us precisely how two particles "talk" to each other.

Now, what about three particles? Or a thousand? This is where the genius of Hans Bethe comes in. In 1931, he proposed a breathtakingly elegant solution, the ​​Bethe ansatz​​ (which is just a fancy German word for "guess"). He wagered that the wavefunction for any number of particles could be constructed using only this two-body phase shift. The central idea is that any many-body collision can be broken down into a sequence of independent two-body collisions. Amazingly, for delta-function interactions in one dimension, this works perfectly.

This "guess" leads to a set of conditions that the particle momenta must satisfy, the famous ​​Bethe Ansatz Equations​​. In their logarithmic form, for $N$ particles on a ring of length $L$, they look like this:

where the $I_j$ are a set of integer quantum numbers that label the state. Look closely! The second term on the right is just a sum of the two-body phase shifts we found earlier. What this equation tells us is profound. The allowed momentum for particle $j$, $k_j$, is no longer a simple, independent value. It depends on the momenta of every other particle in the system. It's a deeply collective, quantum-mechanical democracy, where the state of each particle is intimately linked to all others through these pairwise "conversations."

These equations are notoriously difficult to solve in general, but we can explore them in two opposite limits of interaction strength $g$, revealing two completely different physical pictures.

The Weakly-Interacting Gas

Let's start with a very weak repulsion, a gentle nudge ($g \to 0$). Here, we can treat the interaction as a small perturbation. For non-interacting bosons ($g=0$), the ground state is a Bose-Einstein condensate (BEC), where all $N$ particles happily pile into the zero-momentum state. What happens when we turn on a tiny $g$? The particles are now slightly repelled by each other, and this costs energy. Using standard perturbation theory, we can calculate the energy shift. To first order, the ground state energy increases by:

This makes perfect sense. The energy cost is proportional to the interaction strength $g$ and the total number of interacting pairs, which is $\frac{N(N-1)}{2}$. The system is still very much a gas of bosons, but it's a bit less comfortable with all the particles being in the same spot. (Nature, in its subtlety, also adds a negative second-order correction, a hint that even this simple limit has hidden depths!)

The Strongly-Interacting Gas: Fermionization

Now for the opposite extreme: let's crank up the repulsion to infinity ($g \to \infty$). The particles now despise each other, acting like impenetrable points. We're in the ​​Tonks-Girardeau limit​​. What does our magic phase shift formula tell us now? As $g \to \infty$, the argument of the arctangent function also goes to infinity, and $\arctan(z) \to \pi/2$. So, the phase shift becomes:

A phase shift of $\pi$ or $-\pi$ is equivalent to multiplying the wavefunction by $e^{\pm i\pi} = -1$. This means that if you swap the positions of any two particles, the total wavefunction flips its sign: $\Psi(\dots, x_i, \dots, x_j, \dots) = -\Psi(\dots, x_j, \dots, x_i, \dots)$.

But wait! This property—the wavefunction changing sign upon particle exchange—is the defining characteristic of ​​fermions​​! This is a spectacular result. In one dimension, bosons with infinitely strong repulsion behave exactly as if they were non-interacting fermions. This phenomenon is called ​​fermionization​​. The bosons haven't turned into fermions; rather, the intense repulsion forces them to arrange themselves in a way that mimics the Pauli exclusion principle, which forbids two fermions from occupying the same quantum state.

We can see this directly in the energy. In the $g \to \infty$ limit, the Bethe equations simplify dramatically, and we can solve them. For two particles on a ring of length $L$, the ground state energy is found to be $E_0 = \frac{\pi^2\hbar^2}{mL^2}$. For three particles, it's $E_0 = \frac{4\pi^2\hbar^2}{mL^2}$ (for $L=1, \hbar=1, m=1/2$, this is $8\pi^2$). In both cases, these are precisely the ground state energies you would calculate for two or three non-interacting, spinless fermions on a ring.

This astounding correspondence holds even for a macroscopic number of particles. In the thermodynamic limit (huge $N$ and $L$, constant density $\rho = N/L$), the ground state energy per particle becomes:

This is none other than the famous ground state energy of a one-dimensional Fermi gas. The bosonic system has completely disguised itself as its fermionic counterpart.

This is a beautiful theoretical story, but could an experimentalist ever see it? You can't directly measure a wavefunction's sign flip. You need a more tangible signature. This is where the ​​two-particle correlation function​​, $g^{(2)}(z)$, comes in. It answers a simple question: "If I find a particle at some position, what is the relative probability of finding another one a distance $z$ away?" The value at zero separation, $g^{(2)}(0)$, tells us how much the particles like to cluster.

For typical, non-interacting bosons, they love to bunch up (a phenomenon called "bunching"), and $g^{(2)}(0) = 2$. For fermions, the exclusion principle keeps them apart, so $g^{(2)}(0) = 0$. What about our Lieb-Liniger bosons?

Here we can use a wonderfully elegant piece of physics known as the ​​Hellmann-Feynman theorem​​. It provides a clever link between the system's energy and its microscopic structure. The theorem states that the derivative of the system's energy with respect to a parameter in the Hamiltonian (like our interaction strength $g$) is equal to the expectation value of the derivative of the Hamiltonian itself. For our model, this works out to a simple and powerful relation: the rate of change of the energy density with respect to the interaction strength is directly proportional to $g^{(2)}(0)$.

By using known results for how the energy behaves in the strong-coupling regime, we can turn this relationship around and calculate the correlation function. The result for large, but not infinite, interaction strength $\gamma$ (a dimensionless version of $g$) is striking:

As the interaction strength $\gamma$ becomes very large, $g^{(2)}(0)$ rushes towards zero! This is the smoking gun of fermionization. The stronger the repulsion, the lower the probability of finding two bosons at the same spot, until in the limit, they perfectly avoid each other, just like fermions. This very behavior—the suppression of $g^{(2)}(0)$ with increasing repulsion—has been exquisitely measured in experiments with ultracold atoms trapped in narrow, one-dimensional "tubes," providing a stunning confirmation of this half-century-old theory. From a simple delta-function potential, we have uncovered a rich tapestry of quantum behavior, where particles can magically transform their statistical nature, a testament to the profound and unified beauty of the laws of physics.

Now that we have painstakingly assembled the beautiful machinery of the Lieb-Liniger model and admired the elegance of its exact Bethe ansatz solution, a natural and pressing question arises: What is it good for? Is it merely a pristine exhibition piece in the museum of mathematical physics, or is it a working tool, a key that unlocks doors to real-world phenomena? The answer, you will be delighted to find, is a resounding "yes" to the latter. The model is not an isolated island; it is a continental hub, a crossroads where paths from seemingly distant fields of science meet in the most surprising and illuminating ways.

In this chapter, we will embark on a journey to explore this rich landscape of applications. We will see how this "simple" model of bosons on a line provides the definitive language for describing cutting-edge experiments with ultracold atoms, how it serves as a theoretical laboratory for the wild frontier of non-equilibrium physics, and, most astonishingly, how it secretly underpins the description of phenomena as diverse as the random growth of a bacterial colony and the winding path of a polymer in a disordered medium. Let us begin.

The most direct and spectacular validation of the Lieb-Liniger model has come from the realm of ultracold atomic physics. Experimentalists have become masters of confinement, using lasers and magnetic fields to create "atomic waveguides"—cigar-shaped traps so tight in two dimensions that quantum particles are forced to move effectively only in the third, creating a near-perfect one-dimensional world. When bosonic atoms are cooled to nanokelvin temperatures in these traps, they constitute a near-perfect realization of the Lieb-Liniger gas.

Taming Quantum Correlations

What does it mean to be a boson? A key characteristic is their gregarious nature. Unlike their antisocial fermion cousins, which obey the Pauli exclusion principle, bosons prefer to be in the same state, to bunch together. For a gas of non-interacting bosons, the probability of finding two particles at the same location is twice that of finding them far apart. We can quantify this with the local pair correlation function, $g^{(2)}(0)$, which is 2 for an ideal Bose gas.

But what happens when we switch on the repulsion described by the Lieb-Liniger model? The particles, while still bosons, are now told to keep their distance. In the extreme limit of infinite repulsion, the so-called Tonks-Girardeau gas, the particles become impenetrable. The probability of finding two at the same point plummets to zero, $g^{(2)}(0) = 0$. The bosons have been "fermionized" by the strong interactions; they behave as if an exclusion principle has been imposed upon them.

This is where the power of the Lieb-Liniger model shines. It doesn't just describe these two extremes; it describes the entire crossover with exquisite precision. Using the model, one can derive a direct relationship between the system's ground state energy and this local correlation. The result is a beautiful and testable prediction: in the regime of strong, but not infinite, repulsion (large interaction parameter $\gamma$), the tendency for bosons to avoid each other is nearly perfect, with a small residual probability of overlap that vanishes as $g^{(2)}(0) \propto 1/\gamma^2$. This is not just a theoretical nicety; it is a quantitative law that has been confirmed in experiments, providing a powerful demonstration of how interactions can fundamentally alter the quantum statistics of a many-body system.

The Symphony of Excitations

A quantum system is more than its ground state; its true character is revealed in how it responds to being disturbed. If you "pluck" the Lieb-Liniger gas, what kind of ripples, or excitations, propagate through it? The Bethe ansatz solution provides a complete answer. It tells us that the elementary excitations are not simply individual atoms moving around, but collective entities called "quasiparticles."

The model allows us to calculate the full energy-momentum relationship, or dispersion relation, for these quasiparticles. For instance, in the strongly interacting regime, one type of fundamental excitation is a "hole," created by removing a particle from the fermion-like ground state. The model predicts that this hole moves through the gas as if it were a particle itself, with a well-defined energy and momentum. Most interestingly, the interactions with the surrounding medium modify its properties. While a hole in a non-interacting fermion gas would have a simple quadratic dispersion $\varepsilon \propto p^2$, the Lieb-Liniger solution reveals that interactions "dress" this excitation, changing its effective mass. In the strong coupling limit, this interaction-driven "dressing" of the hole leads to a renormalized effective mass, $m^*$, which is different from the bare particle mass. Understanding these excitation spectra is crucial for predicting thermodynamic properties, transport coefficients, and the dynamic response of the gas, all of which are measurable in cold-atom experiments.

For much of its history, many-body physics was primarily concerned with systems in thermal equilibrium. But the universe is rarely so placid. What happens when a system is violently shaken and thrown far from equilibrium? For most systems, this is an impossibly complex question. But for the Lieb-Liniger model, its "integrability"—a special mathematical property related to an infinite number of conservation laws—makes it a unique and invaluable theoretical laboratory for exploring the wild territory of non-equilibrium dynamics.

After the Quench: A System That Never Forgets

A standard way to explore non-equilibrium physics in cold-atom experiments is the "quantum quench." One prepares a system in the ground state of a given Hamiltonian and then suddenly changes a parameter, like the interaction strength $g_{1D}$. For instance, one can prepare a non-interacting Bose-Einstein condensate ($g_{1D}=0$), where all atoms happily occupy the zero-momentum state, and then abruptly switch on a strong repulsion.

What happens next? A generic, non-integrable system would eventually thermalize. It would chaotically exchange energy among its degrees of freedom until it settled into a boring thermal state described by a single number, its temperature. It would completely forget the details of its initial state. But an integrable system like the Lieb-Liniger gas behaves differently. Because of its vast number of conserved quantities, it cannot fully thermalize. Instead, it relaxes to a novel type of steady state known as a Generalized Gibbs Ensemble (GGE), which retains a detailed memory of its initial conditions. The Lieb-Liniger model has been central to developing and testing this GGE hypothesis. In the quench from a non-interacting gas, the model makes a startling prediction: in the final GGE steady state, the kinetic and potential energies are perfectly equipartitioned, a non-trivial consequence of its hidden symmetries that has no analogue in standard thermal ensembles.

The Laws of Quantum Flow

Armed with the knowledge of integrability, physicists have recently developed a powerful new theory called Generalized Hydrodynamics (GHD). Think of it as an upgraded version of the familiar fluid dynamics of water or air, but tailor-made for integrable quantum systems. Where classical hydrodynamics tracks the local conservation of mass, momentum, and energy, GHD tracks the infinite family of conserved quantities present in models like Lieb-Liniger.

This framework allows us to solve problems that were previously intractable. Consider, for example, the quantum equivalent of a dam breaking: a dense cloud of Lieb-Liniger gas initially confined to one half of a line, which is then allowed to expand into the vacuum on the other half. GHD provides an exact solution for the subsequent evolution of the density and current profiles. It shows that the system evolves in a "self-similar" fashion, depending only on the ratio $x/t$, and provides an exact formula for the particle current that flows out of the initial cloud. GHD can also be used to probe profound properties like superfluidity at finite temperatures, providing expressions for the "normal" and "superfluid" components of the quantum fluid. The Lieb-Liniger model serves as the archetypal platform for developing and benchmarking this exciting new theory of quantum transport.

Perhaps the most profound lesson the Lieb-Liniger model teaches us is one of universality. The same mathematical structures that describe cold atoms on a line appear, as if by magic, in completely different corners of the scientific world. These connections are not mere coincidences; they reveal a deep unity in the mathematical description of nature.

Random Walks and Growing Piles: The KPZ Universe

Imagine a line of fire spreading across a sheet of paper, the edge of a growing bacterial colony, or the deposition of atoms onto a surface. The height of these fluctuating, randomly growing interfaces is described by a famous stochastic equation known as the Kardar-Parisi-Zhang (KPZ) equation. For decades, the statistical properties of KPZ growth were a notoriously difficult problem.

The breakthrough came from an astonishing connection: the problem of calculating the statistical moments of the KPZ height field can be mapped exactly onto a quantum mechanics problem—specifically, the ground state of the attractive Lieb-Liniger gas!. This incredible duality means that the Bethe ansatz solution, originally developed for repulsive bosons, could be deployed to solve long-standing problems in the theory of random growth. A related problem, that of a directed polymer finding its optimal path through a random energy landscape, also maps onto the very same quantum model. The quenched free energy of the polymer, a key quantity in statistical mechanics, turns out to be nothing but the ground state energy of the attractive Lieb-Liniger system. This startling correspondence between quantum integrability and classical stochastic processes is one of the most beautiful discoveries in modern mathematical physics.

The Sound of One Dimension: Luttinger Liquids

At low energies, the world often simplifies. It has long been known that many different one-dimensional quantum systems—interacting electrons, spin chains, and Bose gases—exhibit a universal behavior at low energies, described by an effective field theory called a Luttinger liquid. This theory successfully predicts, for example, that correlation functions decay as power laws. However, the Luttinger liquid theory itself is a description of "what," not "why." It contains phenomenological parameters, like the Luttinger parameter $K$, which it cannot compute from first principles.

Here, the Lieb-Liniger model acts as a "Rosetta Stone." Because it is an exactly solvable microscopic model that becomes a Luttinger liquid at low energies, it allows us to compute the universal parameters of the effective theory directly from the microscopic interaction strength $g_{1D}$. It provides the crucial bridge between the microscopic world of individual particles and the universal, collective phenomena of the low-energy world, turning a phenomenological description into a predictive, quantitative theory.

A Quantum Compass: Berry Phases and Topology

Finally, the Lieb-Liniger model provides insights into the subtle topological and geometric properties of quantum many-body states. Imagine our ring of interacting bosons is threaded by a magnetic flux, via the Aharonov-Bohm effect. If we slowly increase the flux from zero to one "flux quantum" ($h/e$), the system's ground state wavefunction must evolve. While there is a dynamical phase associated with the energy, there is also a purely geometric component known as the Berry phase, which depends only on the path taken in the space of Hamiltonians, not on how fast the path was traversed.

For the Lieb-Liniger ground state, a clever gauge transformation reveals that winding the magnetic flux by one quantum is equivalent to twisting the boundary conditions on the many-body wavefunction by $2\pi$. The model predicts that the resulting many-body Berry phase is not just some complicated number, but is robustly quantized to be exactly $2\pi$ times the number of particles, $\gamma_g = 2\pi N$. This integer multiple hints at an underlying topological structure in the many-body state. It is a collective, quantized response of the entire system, a beautiful manifestation of quantum geometry in an interacting system.

From cold atoms to burning paper, from quantum hydrodynamics to topology, the Lieb-Liniger model continues to surprise and inspire. It is a testament to the power of a simple, elegant idea to cast light on a vast and interconnected scientific landscape, reminding us that sometimes, the deepest truths about our universe can be found by looking very carefully at a line.