Momentum distribution of a freely expanding Lieb-Liniger gas

We numerically study free expansion of a few Lieb-Liniger bosons, which are initially in the ground state of an infinitely deep hard-wall trap. Numerical calculation is carried out by employing a standard Fourier transform, as follows from the Fermi-Bose transformation for a time-dependent Lieb-Liniger gas. We study the evolution of the momentum distribution, the real-space single-particle density, and the occupancies of natural orbitals. Our numerical calculation allows us to explore the behavior of these observables in the transient regime of the expansion, where they are non-trivially affected by the particle interactions. We derive analytically (by using the stationary phase approximation) the formula which connects the asymptotic shape of the momentum distribution and the initial state. For sufficiently large times the momentum distribution coincides (up to a simple scaling transformation) with the shape of the real-space single-particle density (the expansion is asymptotically ballistic). Our analytical and numerical results are in good agreement.Comment: small changes; references correcte

Geometry of quantum observables and thermodynamics of small systems

The concept of ergodicity---the convergence of the temporal averages of observables to their ensemble averages---is the cornerstone of thermodynamics. The transition from a predictable, integrable behavior to ergodicity is one of the most difficult physical phenomena to treat; the celebrated KAM theorem is the prime example. This Letter is founded on the observation that for many classical and quantum observables, the sum of the ensemble variance of the temporal average and the ensemble average of temporal variance remains constant across the integrability-ergodicity transition. We show that this property induces a particular geometry of quantum observables---Frobenius (also known as Hilbert-Schmidt) one---that naturally encodes all the phenomena associated with the emergence of ergodicity: the Eigenstate Thermalization effect, the decrease in the inverse participation ratio, and the disappearance of the integrals of motion. As an application, we use this geometry to solve a known problem of optimization of the set of conserved quantities---regardless of whether it comes from symmetries or from finite-size effects---to be incorporated in an extended thermodynamical theory of integrable, near-integrable, or mesoscopic systems

Polarization Suppression and Nonmonotonic Local Two-Body Correlations in the Two-Component Bose Gas in One Dimension

We study the interplay of quantum statistics, strong interactions and finite temperatures in the two-component (spinor) Bose gas with repulsive delta-function interactions in one dimension. Using the Thermodynamic Bethe Ansatz, we obtain the equation of state, population densities and local density correlation numerically as a function of all physical parameters (interaction, temperature and chemical potentials), quantifying the full crossover between low-temperature ferromagnetic and high-temperature unpolarized regimes. In contrast to the single-component, Lieb-Liniger gas, nonmonotonic behaviour of the local density correlation as a function of temperature is observed.Comment: 4 pages, 6 figure

Invariants of Collective Neutrino Oscillations

We consider the flavor evolution of a dense neutrino gas by taking into account both vacuum oscillations and self interactions of neutrinos. We examine the system from a many-body perspective as well as from the point of view of an effective one-body description formulated in terms of the neutrino polarization vectors. We show that, in the single angle approximation, both the many-body picture and the effective one-particle picture possess several constants of motion. We write down these constants of motion explicitly in terms of the neutrino isospin operators for the many-body case and in terms of the polarization vectors for the effective one-body case. The existence of these constants of motion is a direct consequence of the fact that the collective neutrino oscillation Hamiltonian belongs to the class of Gaudin Hamiltonians. This class of Hamiltonians also includes the (reduced) BCS pairing Hamiltonian describing superconductivity. We point out the similarity between the collective neutrino oscillation Hamiltonian and the BCS pairing Hamiltonian. The constants of motion manifest the exact solvability of the system. Borrowing the well established techniques of calculating the exact BCS spectrum, we present exact eigenstates and eigenvalues of both the many-body and the effective one-particle Hamiltonians describing the collective neutrino oscillations. For the effective one-body case, we show that spectral splits of neutrinos can be understood in terms of the adiabatic evolution of some quasi-particle degrees of freedom from a high density region where they coincide with flavor eigenstates to the vacuum where they coincide with mass eigenstates. We write down the most general consistency equations which should be satisfied by the effective one-body eigenstates and show that they reduce to the spectral split consistency equations for the appropriate initial conditions.Comment: 26 pages with one figure. Published versio

Exact Results for Three-Body Correlations in a Degenerate One-Dimensional Bose Gas

Motivated by recent experiments we derive an exact expression for the correlation function entering the three-body recombination rate for a one-dimensional gas of interacting bosons. The answer, given in terms of two thermodynamic parameters of the Lieb-Liniger model, is valid for all values of the dimensionless coupling $\gamma$ and contains the previously known results for the Bogoliubov and Tonks-Girardeau regimes as limiting cases. We also investigate finite-size effects by calculating the correlation function for small systems of 3, 4, 5 and 6 particles.Comment: 4 pages, 2 figure

More on the exact solution of the O(n) model on a random lattice and an investigation of the case |n|>2

For $n\in [-2,2]$ the $O(n)$ model on a random lattice has critical points to which a scaling behaviour characteristic of 2D gravity interacting with conformal matter fields with $c\in [-\infty,1]$ can be associated. Previously we have written down an exact solution of this model valid at any point in the coupling constant space and for any $n$. The solution was parametrized in terms of an auxiliary function. Here we determine the auxiliary function explicitly as a combination of $\theta$-functions, thereby completing the solution of the model. Using our solution we investigate, for the simplest version of the model, hitherto unexplored regions of the parameter space. For example we determine in a closed form the eigenvalue density without any assumption of being close to or at a critical point. This gives a generalization of the Wigner semi-circle law to $n\neq 0$. We also study the model for $|n|>2$. Both for $n2$ we find that the model is well defined in a certain region of the coupling constant space. For $n<-2$ we find no new critical points while for $n>2$ we find new critical points at which the string susceptibility exponent $\gamma_{str}$ takes the value $+\frac{1}{2}$.Comment: 27 pages, LaTeX file (uses epsf) + 3 eps figures, formulas involving the string susceptibility corrrected, no change in conclusion

Spectra and Symmetry in Nuclear Pairing

We apply the algebraic Bethe ansatz technique to the nuclear pairing problem with orbit dependent coupling constants and degenerate single particle energy levels. We find the exact energies and eigenstates. We show that for a given shell, there are degeneracies between the states corresponding to less and more than half full shell. We also provide a technique to solve the equations of Bethe ansatz.Comment: 15 pages of REVTEX with 2 eps figure

New Formula for the Eigenvectors of the Gaudin Model in the sl(3) Case

We propose new formulas for eigenvectors of the Gaudin model in the \sl(3) case. The central point of the construction is the explicit form of some operator P, which is used for derivation of eigenvalues given by the formula $| w_1, w_2) = \sum_{n=0}^\infty P^n/n! | w_1, w_2,0>$, where $w_1$, $w_2$ fulfil the standard well-know Bethe Ansatz equations

Exact Diagonalisation of The XY-Hamiltonian of Open Linear Chains with Periodic Coupling Constants and Its Application to Dynamics of One-Dimensional Spin Systems

A new method of diagonalisation of the XY-Hamiltonian of inhomogeneous open linear chains with periodic (in space) changing Larmor frequencies and coupling constants is developed. As an example of application, multiple quantum dynamics of an inhomogeneous chain consisting of 1001 spins is investigated. Intensities of multiple quantum coherences are calculated for arbitrary inhomogeneous chains in the approximation of the next nearest interactions. {\it Key words:} linear spin chain, nearest--neighbour approximation, three--diagonal matrices, diagonalisation, fermions, multiple--quantum NMR, multiple--quantum coherence, intensities of multiple--quantum coherences. {\it PACS numbers:} 05.30.-d; 76.20.+qComment: 21 pages + 1 figure (to download separately via ps-format

Exactly Solvable Pairing Model Using an Extension of Richardson-Gaudin Approach

We introduce a new class of exactly solvable boson pairing models using the technique of Richardson and Gaudin. Analytical expressions for all energy eigenvalues and first few energy eigenstates are given. In addition, another solution to Gaudin's equation is also mentioned. A relation with the Calogero-Sutherland model is suggested.Comment: 9 pages of Latex. In the proceedings of Blueprints for the Nucleus: From First Principles to Collective Motion: A Festschrift in Honor of Professor Bruce Barrett, Istanbul, Turkey, 17-23 May 200