Following up on a recent analysis of two cold atoms in a time-dependent harmonic trap in one dimension, we explore the entanglement entropy of two and three fermions in the same situation when driven through a parametric resonance. We find that the presence of such a resonance in the two-particle system leaves a clear imprint on the entanglement entropy. We show how the signal is modified by attractive and repulsive contact interactions, and how it remains present for the three-particle system. Additionaly, we extend the work of recent experiments to demonstrate how restricting observation to a limited subsystem gives rise to locally thermal behavior.
Following up on a recent analysis of two cold atoms in a time-dependent harmonic trap in one dimension, we explore the entanglement entropy of two and three fermions in the same situation when driven through a parametric resonance. We find that the presence of such a resonance in the two-particle system leaves a clear imprint on the entanglement entropy. We show how the signal is modified by attractive and repulsive contact interactions, and how it remains present for the three-particle system. Additionaly, we extend the work of recent experiments to demonstrate how restricting observation to a limited subsystem gives rise to locally thermal behavior.
We investigate quantum correlations in the ground state of noninteracting Fermi gases of N particles trapped by an external space-dependent harmonic potential, in any dimension. For this purpose we compute one-particle correlations, particle fluctuations, and bipartite entanglement entropies of extended space regions, and study their large-N scaling behaviors. The half-space von Neumann entanglement entropy is computed for any dimension, obtaining ${S}_{\mathrm{HS}}\ensuremath{\approx}{c}_{l}{N}^{(d\ensuremath{-}1)/d}\mathrm{ln}N$ analogously to homogenous systems, with cl=1/6,1/(62),1/(66) in one, two, and three dimensions, respectively. We show that the asymptotic large-N relation ${S}_{A}\ensuremath{\approx}{\ensuremath{\pi}}^{2}{V}_{A}/3$, between the von Neumann entanglement entropy SA and particle variance VA of an extended space region A, holds for any subsystem A and in any dimension, analogously to homogeneous noninteracting Fermi gases.
We investigate the entanglement properties of the nonequilibrium dynamics of one-dimensional noninteracting Fermi gases released from a trap. The gas of N particles is initially in the ground state within a hard-wall or a harmonic trap, and then it expands after dropping the trap. We compute the time dependence of the von Neumann and R\'enyi entanglement entropies and the particle fluctuations of spatial intervals around the original trap in the limit of a large number N of particles. The results for these observables apply to one-dimensional gases of impenetrable bosons as well. We identify different dynamical regimes at small and large times, depending also on the initial condition, i.e., whether it is that of a hard-wall or a harmonic trap. In particular, we show analytically that the expansion from hard-wall traps is characterized by the asymptotic small-time behavior $S\ensuremath{\approx}(1/3)\mathrm{ln}(1/t)$ of the von Neumann entanglement entropy and the relation $S\ensuremath{\approx}{\ensuremath{\pi}}^{2}V/3$ (where V is the particle variance), which are analogous to the equilibrium behaviors whose leading logarithms are essentially determined by the corresponding conformal field theory with central charge c=1. The time dependence of the entanglement entropy of extended regions during the expansion from harmonic traps has a remarkable property in that it can be expressed as a global time-dependent rescaling of the space dependence of the initial equilibrium entanglement entropy.
We investigate finite-number effects in collisions between two states of an initially well-known number of identical bosons with contact interactions, oscillating in the presence of harmonic confinement in one dimension. We investigate two N/2 (interacting) ground states, which are initially displaced from the trap center, and the effects of varying interaction strength. The numerics focus on the simplest case of N=4. In the noninteracting case, such a system would display periodic oscillation with a half harmonic oscillator period (due to the left-right symmetry). With the addition of contact interactions between the bosons, collisions generate entanglement between each of the states and distribute energy into other modes of the oscillator. We study the system numerically via an exact diagonalization of the Hamiltonian with a finite basis, investigating left-right number uncertainty as our primary measure of entanglement. Additionally, we study the time evolution and equilibration of the single-body von Neumann entropy for both the attractive and repulsive cases. We identify parameter regimes for which attractive interactions create behavior qualitatively different from that of repulsive interactions, due to the presence of bound states (quantum solitons), and explain the processes behind this.
Motivated by recent development of the concept of the disorder operator and its relation with entanglement entropy in bosonic systems, here we show the disorder operator successfully probes many aspects of quantum entanglement in fermionic many-body systems. From both analytical and numerical computations in free and interacting fermion systems in 1D and 2D, we find the disorder operator and the entanglement entropy exhibit similar universal scaling behavior, as a function of the boundary length of the subsystem, but with subtle yet important differences. In 1D they both follow the \log{L} log L scaling behavior with the coefficient determined by the Luttinger parameter for disorder operator, and the conformal central charge for entanglement entropy. In 2D they both show the universal L\log L L log L scaling behavior in free and interacting Fermi liquid states, with the coefficients depending on the geometry of the Fermi surfaces. However at a 2D quantum critical point with non-Fermi-liquid state, extra symmetry information is needed in the design of the disorder operator, so as to reveal the critical fluctuations as does the entanglement entropy. Our results demonstrate the fermion disorder operator can be used to probe quantum many-body entanglement related to global symmetry, and provide new tools to explore the still largely unknown territory of highly entangled fermion quantum matter in 2 or higher dimensions.
Investigation of the emergence of entanglement in a classical two-particle Brownian system through nonequilibrium driving is carried out. The Brownian particles are in contact with independent thermal bath, where one particle is driven by a time-dependent potential and the other remains in a static harmonic trap. The particles are harmonically coupled, enabling thermal correlation and information exchange. Using stochastic simulations, we compute an entanglement witness based on osmotic velocities and position correlations. We show that stronger driving enhances entanglement generation, with protocol-specific features. Thermodynamic quantities such as entropy production and the Jarzynski equality are also examined, revealing a close link between nonequilibrium control and classical entanglement in mesoscopic systems.
The thermodynamic and magnetic properties of a spin-mixture of 40K atoms in a 1D harmonic trap, in the presence of an external magnetic field H, are studied using the static fluctuation approximation (SFA). The two-body interaction used is the contact potential. The thermodynamic properties include the chemical potential, energy per particle, specific heat capacity, and entropy per particle. The magnetic properties comprise the magnetization per particle and susceptibility per particle. A phase transition is observed for the repulsive interactions around 0.35TF and 0.5TF in the channels {–9/2, –7/2} and {–9/2, –5/2}, respectively. A peak appears in the entropy and susceptibility at effective scattering length |a1d| ≈ 0.16 μm. For the attractive interactions, three different behaviors depending on H are found in the low-temperature limit. The first behavior is when H < Hc, where Hc is the critical magnetic field, and typical peaks emerge. The second behavior is at H = Hc, and a magnetic phase transition occurs. The third behavior is when H > Hc and the system becomes a partially polarized phase. Our results are consistent with the reported results obtained by the local density approximation, thermal Bethe ansatz method, mean-field approximation, and quantum transfer matrix method.