The coupling between the spin degrees of freedom and the orbital angular momentum has a profound effect on the properties of nuclei, atoms, and condensed-matter systems. Recently, synthetic gauge fields have been realized experimentally in neutral cold-atom systems, giving rise to a spin-orbit coupling term with ``strength''
kso. This paper investigates the interplay between the single-particle spin-orbit coupling term of Rashba type and the short-range two-body
s-wave interaction for cold atoms under external confinement. Specifically, we consider two different harmonically trapped two-atom systems. The first system consists of an atom with spin-orbit coupling that interacts with a structureless particle through a short-range two-body potential. The second system consists of two atoms that both feel the spin-orbit coupling term and that interact through a short-range two-body potential. Treating the spin-orbit term perturbatively, we determine the correction to the ground-state energy for various generic parameter combinations. Selected excited states are also treated. An important aspect of our study is that the perturbative treatment is not limited to small
s-wave scattering lengths but provides insights into the system behavior over a wide range of scattering lengths, including the strongly interacting unitary regime. We find that the interplay between the spin-orbit coupling term and the
s-wave interaction generically enters, depending on the exact parameter combinations of the
s-wave scattering lengths, at order
kso2 or
kso4 for the ground state and leads to a shift of the energy of either sign. While the absence of a term proportional to
kso follows straightforwardly from the functional form of the spin-orbit coupling term, the absence of a term proportional to
kso2 for certain parameter combinations is unexpected. The well-known fact that the spin-orbit coupling term couples the relative and center-of-mass degrees of freedom has interesting consequences for the trapped two-particle systems. For example, we find that, for certain parameter combinations, the spin-orbit coupling term turns sharp crossings into avoided crossings with an energy splitting proportional to
kso. Our perturbative results are confirmed by numerical calculations that expand the eigenfunctions of the two-particle Hamiltonian in terms of basis functions that contain explicitly correlated Gaussians.
We provide a detailed study of the properties of a few interacting spin 1 / 2 fermions trapped in a one-dimensional harmonic oscillator potential. The interaction is assumed to be well represented by a contact delta potential. Numerical results obtained by means of direct diagonalization techniques are combined with analytical expressions for both the non-interacting and strongly interacting regime. The N = 2 case is used to benchmark our numerical techniques with the known exact solution of the problem. After a detailed description of the numerical methods, in a tutorial-like manner, we present the static properties of the system for N = 2 , 3 , 4 and 5 particles, e.g., low-energy spectrum, one-body density matrix, ground-state densities. Then, we consider dynamical properties of the system exploring first the excitation of the breathing mode, using the dynamical structure function and corresponding sum-rules, and then a sudden quench of the interaction strength.
We study spin-1/2 fermions in spin-dependent potentials under the spin-model approximation, in which interatomic collisions that change the total occupation of single-particle modes are ignored. The spin-model approximation maps the interacting fermion problem to an ensemble of lattice spin models in energy space, where spin-spin interactions are long range and spin anisotropic. We consider a Ramsey-type sequence in which quantum degenerate, initially spin-polarized fermions are rotated to a spin superposition, and the differences in potential for the two spin states induce
s-wave interactions. Applied to such a setup, we show that the spin-model approximation qualitatively captures the collective spin dynamics for weak interactions compared to the harmonic-oscillator frequency, and often captures quantitative features at timescales much longer than would be expected from perturbation theory. We explore corrections to the spin model, and the relative importance of corrections for both collective spin and fidelity metrics when realistic anharmonic potential corrections are taken into account, demonstrating that anharmonicity can significantly improve the fidelity of the spin-model approximation. Additionally, we present numerical techniques that are useful for analysis of spin models on an energy lattice, including enacting a change of single-particle basis on a many-body state as an effective time evolution, and fitting of spatially inhomogeneous long-range interactions with exponentials. This latter technique is useful for constructing matrix product operators for use in density-matrix renormalization group analyses, and may have broader applicability within the tensor network community.
An accurate description of the nonequilibrium dynamics of systems with coupled spin and bosonic degrees of freedom remains theoretically challenging, especially for large system sizes and in higher than one dimension. Phase-space methods such as the truncated Wigner approximation (TWA) have the advantage of being easily scalable and applicable to arbitrary dimensions. In this work we adapt the TWA to generic spin-boson models by making use of recently developed algorithms for discrete phase spaces [J. Schachenmayer, A. Pikovski, and A. M. Rey, Phys. Rev. X 5, 011022 (2015)]. Furthermore we go beyond the standard TWA approximation by applying a scheme based on the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy of equations to our coupled spin-boson model. This allows us, in principle, to study how systematically adding higher-order corrections improves the convergence of the method. To test various levels of approximation we study an exactly solvable spin-boson model, which is particularly relevant for trapped-ion arrays. Using TWA and its BBGKY extension we accurately reproduce the time evolution of a number of one- and two-point correlation functions in several dimensions and for an arbitrary number of bosonic modes.
We employ ab initio methods of quantum chemistry to investigate spin-1/2 fermions interacting via a two-body contact potential in a one-dimensional harmonic trap. The convergence of the total energy with the size of the one-particle basis set is analytically investigated for the two-body problem and the same form of the convergence formula is numerically confirmed to be valid for the many-body case. Benchmark calculations for two to six fermions with the full configuration interaction method equivalent to the exact diagonalization approach, and the coupled cluster (CC) method including single, double, triple, and quadruple excitations are presented. The convergence of the correlation energy with the level of excitations included in the CC model is analyzed. The range of the interaction strength for which single-reference CC methods work is examined. Next, the CC method restricted to single, double, and noniterative triple excitations, CCSD(T), is employed to study a two-component Fermi gas composed of 6–80 atoms in a one-dimensional harmonic trap. The density profiles of trapped atomic clouds are also reported. Finally, a comparison with experimental results for few-fermion systems is presented. Upcoming possible applications and extensions of the presented approach are discussed.
We theoretically investigate the properties of two interacting ultracold highly magnetic atoms trapped in a one-dimensional harmonic potential. The atoms interact via an anisotropic long-range dipole-dipole interaction, which in one dimension effectively can be modeled by the contact interaction. We investigate the interplay of the external magnetic field, the spin-spin interaction, and the trapping potential and how they affect the magnetization of the system. We show the role of indistinguishability and symmetries in the dynamics by studying the time evolution of the observables that could be measured experimentally. The presented model may depict the on-site interaction of the extended Hubbard models, therefore giving a better understanding of the fundamental building block of the respective many-body quantum simulators.