CONAN - The cruncher of local exchange coefficients for strongly interacting confined systems in one dimension

Abstract We consider a one-dimensional system of particles with strong zero-range interactions. This system can be mapped onto a spin chain of the Heisenberg type with exchange coefficients that depend on the external trap. In this paper, we present an algorithm that can be used to compute these exchange coefficients. We introduce an open source code CONAN (Coefficients of One-dimensional N -Atom Networks) which is based on this algorithm. CONAN works with arbitrary external potentials and we have tested its reliability for system sizes up to around 35 particles. As illustrative examples, we consider a harmonic trap and a box trap with a superimposed asymmetric tilted potential. For these examples, the computation time typically scales with the number of particles as O ( N 3.5 ± 0.4 ) . Computation times are around 10 s for N = 10 particles and less than 10 min for N = 20 particles. Program summary Program title: CONAN Program Files doi: http://dx.doi.org/10.17632/tw87vdy68b.1 Licensing provisions: GNU General Public License 3 (GPL) Programming language: C Nature of problem: A system of N atoms (fermions or bosons) with a strong zero-range interaction confined in a one-dimensional potential V ( x ) can be described using a spin chain Heisenberg type Hamiltonian. This effective spin chain Hamiltonian is defined through N − 1 exchange coefficients (also called geometric coefficients, α k ). The exchange coefficients depend only on the integer N and the function V ( x ) , but each coefficient is formally given as an ( N − 1 )-dimensional integral. Given a number of particles N and a confining potential V ( x ) , we wish to compute the exchange coefficients, but carrying out the ( N − 1 )-dimensional integral numerically is not, to say the least, a method that scales well with the system size. Solution method: We wish to compute the exchange coefficients for a given system, but to do this we need to express them in a way that is more well-suited for a numerical implementation, i.e. in a way that does not involve an ( N − 1 )-dimensional integral. In the submitted manuscript, we derive such an expression for the exchange coefficients. Our program, CONAN, is the numerical implementation of this formula for the exchange coefficients. Thus, CONAN takes as physical inputs the system size N and a smooth potential V ( x ) , and returns the corresponding N − 1 exchange coefficients appearing in the spin chain Hamiltonian.

[1]  U. Schollwoeck The density-matrix renormalization group in the age of matrix product states , 2010, 1008.3477.

[2]  D. Blume,et al.  Correlations of the upper branch of 1D harmonically trapped two-component fermi gases. , 2013, Physical review letters.

[3]  C. Forss'en,et al.  Fermionization of two-component few-fermion systems in a one-dimensional harmonic trap , 2013, 1304.2992.

[4]  S. Jochim,et al.  Deterministic Preparation of a Tunable Few-Fermion System , 2011, Science.

[5]  D. Weiss,et al.  A quantum Newton's cradle , 2006, Nature.

[6]  S. Reimann,et al.  Quantum magnetism without lattices in strongly interacting one-dimensional spinor gases , 2013, 1310.3705.

[7]  U. Schollwoeck The density-matrix renormalization group , 2004, cond-mat/0409292.

[8]  F. Altomare,et al.  One-Dimensional Superconductivity in Nanowires , 2013 .

[9]  Thierry Giamarchi,et al.  Quantum physics in one dimension , 2004 .

[10]  D. Petrosyan,et al.  Engineering the dynamics of effective spin-chain models for strongly interacting atomic gases , 2014, 1408.3414.

[11]  Mattias Gustavsson,et al.  Pinning quantum phase transition for a Luttinger liquid of strongly interacting bosons , 2010, Nature.

[12]  Maxim Olshanii Atomic Scattering in the Presence of an External Confinement and a Gas of Impenetrable Bosons , 1998 .

[13]  P. Massignan,et al.  Magnetism in Strongly Interacting One-Dimensional Quantum Mixtures. , 2015, Physical review letters.

[14]  S. Reimann,et al.  Antiferromagnetic Heisenberg Spin Chain of a Few Cold Atoms in a One-Dimensional Trap. , 2015, Physical review letters.

[15]  A. Jensen,et al.  Comparing and contrasting nuclei and cold atomic gases , 2013, 1303.7351.

[16]  J. Dalibard,et al.  Many-Body Physics with Ultracold Gases , 2007, 0704.3011.

[17]  M. Girardeau,et al.  Relationship between Systems of Impenetrable Bosons and Fermions in One Dimension , 1960 .

[18]  S. Jochim,et al.  Fermionization of two distinguishable fermions. , 2011, Physical review letters.

[19]  Immanuel Bloch,et al.  Tonks–Girardeau gas of ultracold atoms in an optical lattice , 2004, Nature.

[20]  V. Saligrama,et al.  Experimental realization of long-distance entanglement between spins in antiferromagnetic quantum spin chains , 2015, Nature Physics.

[21]  F. Schäfer,et al.  A one-dimensional liquid of fermions with tunable spin , 2014, Nature Physics.

[22]  Li Yang,et al.  Strongly interacting quantum gases in one-dimensional traps , 2015, 1502.01706.

[23]  Ericka Stricklin-Parker,et al.  Ann , 2005 .

[24]  Minoru Takahashi Thermodynamics of One-Dimensional Solvable Models , 1999 .

[25]  Xiaoling Cui,et al.  Effective spin-chain model for strongly interacting one-dimensional atomic gases with an arbitrary spin , 2015, 1510.06087.

[26]  T. Ho,et al.  Spin-orbit-coupled one-dimensional Fermi gases with infinite repulsion , 2013, 1310.1993.

[27]  F. Deuretzbacher,et al.  Exact solution of strongly interacting quasi-one-dimensional spinor Bose gases. , 2007, Physical review letters.

[28]  S. Blundell,et al.  Organic and molecular magnets , 2004 .

[29]  S. Jochim,et al.  From Few to Many: Observing the Formation of a Fermi Sea One Atom at a Time , 2013, Science.

[30]  A. Jensen,et al.  Strongly interacting confined quantum systems in one dimension , 2013, Nature Communications.

[31]  Toshiya Kinoshita,et al.  Observation of a One-Dimensional Tonks-Girardeau Gas , 2004, Science.

[32]  Zhong-Qi Ma,et al.  Exact solution for infinitely strongly interacting Fermi gases in tight waveguides. , 2008, Physical review letters.

[33]  Emil H. Eriksen,et al.  Computation of local exchange coefficients in strongly interacting one-dimensional few-body systems: local density approximation and exact results , 2015, 1508.07164.

[34]  L. Tonks The Complete Equation of State of One, Two and Three-Dimensional Gases of Hard Elastic Spheres , 1936 .

[35]  H. Rønnow,et al.  Fractional spinon excitations in the quantum Heisenberg antiferromagnetic chain , 2013, Nature Physics.

[36]  A. Yacoby,et al.  Electron liquids and solids in one dimension , 2010, Nature.

[37]  A. E. Thomsen,et al.  Comparing models for the ground state energy of a trapped one-dimensional Fermi gas with a single impurity , 2015, 1508.05917.

[38]  P. Zoller,et al.  Condensed matter theory of dipolar quantum gases. , 2012, Chemical reviews.

[39]  J. Danzl,et al.  Realization of an Excited, Strongly Correlated Quantum Gas Phase , 2009, Science.

[40]  D. Petrosyan,et al.  Tunable self-assembled spin chains of strongly interacting cold atoms for demonstration of reliable quantum state transfer , 2015, 1512.06722.

[41]  Ogata,et al.  Bethe-ansatz wave function, momentum distribution, and spin correlation in the one-dimensional strongly correlated Hubbard model. , 1990, Physical review. B, Condensed matter.

[42]  T. Esslinger Fermi-Hubbard Physics with Atoms in an Optical Lattice , 2010, 1007.0012.

[43]  Vincent M. Klinkhamer,et al.  Two fermions in a double well: exploring a fundamental building block of the Hubbard model. , 2014, Physical review letters.

[44]  P. Massignan,et al.  Strong-coupling ansatz for the one-dimensional Fermi gas in a harmonic potential , 2014, Science Advances.