Quasi-one-dimensional (Q1D) systems, i.e., three- and two-dimensional (3D/2D) arrays composed of weakly coupled one-dimensional lattices of interacting quantum particles, exhibit rich and fascinating physics. They are studied across various areas of condensed matter and ultracold atomic lattice-gas physics, and are often marked by dimensional crossover as the coupling between one-dimensional systems is increased or temperature decreased, i.e., the Q1D system goes from appearing largely 1D to largely 3D. Phase transitions occurring along the crossover can strongly enhance this effect. Understanding these crossovers and associated phase transitions can be challenging due to the very different elementary excitations of 1D systems compared to higher-dimensional ones. In the present work, we combine numerical matrix product state (MPS) methods with mean-field (MF) theory to study paradigmatic cases of dimensional crossovers and the associated phase transitions in systems of both hard-core and soft-core lattice bosons, with relevance to both condensed matter physics and ultracold atomic gases. We show that the superfluid-to-insulator transition is a first order one, as opposed to the isotropic cases, and calculate transition temperatures for the superfluid states, finding excellent agreement with analytical theory. At the same time, our $\mathrm{MPS}+\mathrm{MF}$ approach keeps functioning well where the current analytical framework cannot be applied. We further confirm the qualitative and quantitative reliability of our approach by comparison to exact quantum Monte Carlo calculations for the full 3D arrays.
Quasi-one-dimensional (Q1D) systems, i.e., three- and two-dimensional (3D/2D) arrays composed of weakly coupled one-dimensional lattices of interacting quantum particles, exhibit rich and fascinating physics. They are studied across various areas of condensed matter and ultracold atomic lattice-gas physics, and are often marked by dimensional crossover as the coupling between one-dimensional systems is increased or temperature decreased, i.e., the Q1D system goes from appearing largely 1D to largely 3D. Phase transitions occurring along the crossover can strongly enhance this effect. Understanding these crossovers and associated phase transitions can be challenging due to the very different elementary excitations of 1D systems compared to higher-dimensional ones. In the present work, we combine numerical matrix product state (MPS) methods with mean-field (MF) theory to study paradigmatic cases of dimensional crossovers and the associated phase transitions in systems of both hard-core and soft-core lattice bosons, with relevance to both condensed matter physics and ultracold atomic gases. We show that the superfluid-to-insulator transition is a first order one, as opposed to the isotropic cases, and calculate transition temperatures for the superfluid states, finding excellent agreement with analytical theory. At the same time, our MPS+MF approach keeps functioning well where the current analytical framework cannot be applied. We further confirm the qualitative and quantitative reliability of our approach by comparison to exact quantum Monte Carlo calculations for the full 3D arrays.
Correlated electron states are at the root of many important phenomena including unconventional superconductivity (USC), where electron pairing arises from repulsive interactions.Computing the properties of correlated electrons, such as the critical temperature T c for the onset of USC, efficiently and reliably from the microscopic physics with quantitative methods remains a major challenge for almost all models and materials.In this theoretical work, we combine matrix product states (MPS) with static mean field (MF) to provide a solution to this challenge for quasi-one-dimensional (Q1D) systems: two-and threedimensional materials comprised of weakly coupled correlated 1D fermions.This MPS þ MF framework for the ground state and thermal equilibrium properties of Q1D fermions is developed and validated for attractive Hubbard systems first, and further enhanced via analytical field theory.We then deploy it to compute T c for superconductivity in 3D arrays of weakly coupled, doped, and repulsive Hubbard ladders.The MPS þ MF framework thus enables the quantitative study of USC and high-T c superconductivity-and potentially many more correlated phases-in fermionic Q1D systems based directly on their microscopic parameters, in ways inaccessible to previous methods.This approach further allows one to treat competing macroscopic orders, such as superconducting and insulating ones, on an equal footing.Benchmarks of the framework using auxiliary-field quantum Monte Carlo techniques show that the overestimation of, e.g., T c due to its mean-field component, is near constant in microscopic parameters.These features of the MPS þ MF approach to correlated fermions open up the possibility of designing deliberately optimized Q1D superconductors, from experiments in ultracold gases to synthesizing new materials.
Pairing and topological phases in cold atoms with long range interactions Chapter 1 Introduction The systematic study and classification of phase transitions has become a popular research area in physics over the past decades. The most familiar phenomena associated to the concept of phase transition are the ones giving rise to macroscopic changes in the properties of strongly correlated many-body systems due to thermal fluctuations. The most prominent examples of this kind of behaviour are the ferromagnetic-paramagnetic transition in lattice spin models and the liquid-gas transition. On the other hand, when dealing with quantum-mechanical systems changes of matter from one state to another can arise as the result of the variation of parameters other than temperature. In particular, when a given many-body quantum problem is studied at zero temperature, thereby suppressing thermal fluctuations, it is the interplay between kinetic energy and interactions that can cause the system to jump into different phases as the model parameters are varied. Phase transitions occurring at zero temperature are widely known as quantum phase transitions. In the framework of the physics of phase transitions, the role of dimensionality is a noticeable one. The major challenge in this perspective is the treatment of low dimensional systems, where standard mean-field treatments and ordinary perturbation theory are known to fail due to the enhancement of fluctuations. The goal of describing collective behaviours in reduced dimensionality is usually pursued in one-dimensional (1D) systems, which turn out to be simple enough to be studied in an effective way through both numerical and analytical techniques, while exhibiting a rich and deeply fascinating phenomenology. As a matter of fact, since particles confined to one dimension cannot avoid interactions with each other, the typical low energy excitations of the system are represented by collective density waves, which substitute nearly free quasiparticle excitations characterizing the behaviour of three-dimensional Fermi liquids in the description of the effects of the interplay between kinetic fluctuations and interactions. On the other hand, the tools employed in the activity of tackling 1D manybody problems are efficient and well developed. The analytical techniques aremainly based on the bosonization technique, which allows for the reformulation of strongly correlated 1D systems onto free bosonic theories, which can then be analyzed by means of standard path integral techniques. Meanwhile, the density matrix renormalization group (DMRG) algorithm represents the state-of-theart numerical machinery in the study of 1D physics, since it allows to efficiently extract the low energy properties of the model Hamiltonian and consequently characterize the corresponding phase diagram by means of the behaviour of properly chosen order parameters. As a final remark pointing out the relevance of the research effort in the direction of a better understanding of 1D many-body quantum systems, it is worth noticing that 1D systems do not represent only a playground for theorists, but have been experimentally realized by means of setups with cold atoms trapped in optical lattices, where the experimentalist is able to tune the model parameters by varying the characteristic features of the trapping laser beams and meanwhile has a comfortable access to the most relevant observables.
Abstract The Bose–Einstein condensation (BEC) in three-dimensional (3D) anisotropic lattices is studied. We present theoretical results for the critical temperature for BEC, chemical potential, condensate fraction and relevant thermodynamic quantities like: internal energy, entropy, specific heat and compressibility as a function of anisotropy parameter being the ratio of the nearest-neighbor in-plane ( t∥ t ‖ ) and out-of-plane ( t⊥ t ⊥ ) hopping amplitudes. In particular, considered scenarios include weakly coupled two-dimensional (2D) planes ( t⊥/t∥≪1 t ⊥ / t ‖ ≪ 1 , relevant for layered structures) as well as a rod-like geometry of interacting one-dimensional (1D) chains ( t∥/t⊥≪1 t ‖ / t ⊥ ≪ 1 ). The impact of the dimensional crossover as the system is tuned away from a set of disconnected 2D layers, or traverses from a set of separate 1D chains to a regime where a fully isotropic 3D structure emerges is elucidated. Both numerical and analytic approaches are employed, (the latter in a form of series expansions involving t∥,t⊥ t ‖ , t ⊥ amplitudes) for internal energy, entropy, specific heat and isothermal compressibility. The theoretical outcome of the present study may be of interest to a number of scenarios in solid-state physics, where the relevant quasi-particles are bosonic-like, as well as might be applicable to the physics of cold bosons loaded in artificially engineered 3D optical lattices.
We investigate three different types of quantum phase transition occurring in quasi one-dimensional systems theoretically and numerically. First, we study the band-insulator to Mott-insulator transition occurring in charge-transfer complexes, for which the half-filled one-dimensional ionic Hubbard model is considered to be the prototype model. The study is carried out by first deriving an effective spin-one model, and then studying the model numerically using the density matrix renormalization group. We perform a careful finite-size scaling analysis of the mass gaps, order-parameters, and relative susceptibility. We confirm the existence of two quantum critical points. Analysis of the critical exponents confirms that the band-insulator-to-spontaneously-dimerized phase transition belong to the 2D Ising class. The spontaneously dimerized phase undergoes a phase transition to the Mott-insulator which is an infinite-order. Second, we investigate the Mott metal-insulator transition for the half-filled Hubbard model with both nearest-neighbor t and next-nearest-neighbor t′ hopping terms. We study the model using the bosonization approach and density matrix renormalization group simulations. An effective low-energy Hamiltonian that describes the insulator-metal transition is derived. We present results of density matrix renormalization group calculations of spin and charge distribution in various sectors of the phase diagram. The numerical results support the picture derived from the effective theory and give evidence for the complete separation of the transitions involving the spin and the charge degrees of freedom. Finally, we investigate quantum phase transitions phases in low-dimensional fermionic and spin models that go from uniform to spatially inhomogeneous, i.e., dimerized, trimerized, or incommensurate, phases. We propose a new approach based on studying the length dependence of the von Neumann entropy and its corresponding Fourier spectrum for finite segments in the ground state of finite chains. Peaks at a nonzero wave vector are indicators of oscillatory behavior in decaying correlation functions and also provide significant information about certain relevant features of the excitation spectrum; in particular, they can identify the wave vector of soft modes in critical models.
The density matrix renormalization group (DMRG) is a numerical method for studying low dimensional strongly correlated models. The study of these models is one of the most active areas of condensed matter physics today and has been spurred not only by the discovery of high temperature superconductivity but also by experiments on new compounds. In this thesis we study fundamentals of the density matrix renormalization group method as well as applications of the method to Heisenberg spin chains. The introductory text is meant as a brief overview of the DMRG as well as Heisenberg spin systems and provides an introduction to the papers. In papers I and II we investigate the density matrix renormalization group with the aim to better understand the underlying structure of the method. We find that if the renormalization converges to a fixed point in the thermodynamic limit, the resulting basis states can be written as a matrix product form. We also show an explicit relation of the DMRG to variational calculations using the fact that the matrix product states can be rederived through a variational ansatz without any reference to the DMRG calculation. It is then shown how the matrix product states can be used to define a set of Bloch waves for describing elementary excitations of the Heisenberg spin-1 chain. The low lying energy spectrum is calculated. In papers III and IV we investigate the thermodynamics of impurities in a spin-1/2 chain. For this purpose we have shown that a recently developed transfer matrix DMRG method can be used to calculate impurity properties. Both impurity corrections to thermodynamic properties and local properties at finite temperature can be obtained. The method gives quicker and more accurate results than quantum Monte Carlo. The numerical results show good agreement with field theory calculations for a number of impurity configurations and the cross-over function for the magnetic susceptibility from high- to low-temperature behavior shows the expected data collapse. Our results also give quantitative predictions for experiments on doped spin-1/2 chains.
Superfluid and superconducting phase transitions of interacting bosons are observed in diverse physical systems such as helium, superconducting materials, ultra-cold atomic gases, etc. For all these systems one-, two-, and three-dimensional (1D, 2D, and 3D) realizations are experimentally accessible. However, the current experimental focus is shifting to the low dimensional systems (e.g., films and wires) due to facilitated measurement techniques and a constant demand of technological applications. From a theoretical standpoint large fluctuations in 1D and 2D systems make them a challenging problem, especially when disorder is present. Here we study peculiarities associated with disorder in low-dimensional cases. Specifically, we are interested in the class of the systems where superfluid or superconducting transitions can be mapped to the finite temperature 2D superfluid-normal Kosterlitz-Thouless (KT) transition: (a) 2D bosonic system with a strong diagonal disorder, (b) 1D system of lattice bosons with off-diagonal disorder in the limit of large integer filling factor at zero temperature, and (c) 3D layered superconductors with no tunnelling between layers in zero external magnetic field.