Fermionization is what happens to the state of strongly interacting repulsive bosons interacting with contact interactions in one spatial dimension. Crystallization is what happens for sufficiently strongly interacting repulsive bosons with dipolar interactions in one spatial dimension. Crystallization and fermionization resemble each other: in both cases - due to their repulsion - the bosons try to minimize their spatial overlap. We trace these two hallmark phases of strongly correlated one-dimensional bosonic systems by exploring their ground state properties using the one- and two-body density matrix. We solve the N-body Schrödinger equation accurately and from first principles using the multiconfigurational time-dependent Hartree for bosons (MCTDHB) and for fermions (MCTDHF) methods. Using the one- and two-body density, fermionization can be distinguished from crystallization in position space. For N interacting bosons, a splitting into an N-fold pattern in the one-body and two-body density is a unique feature of both, fermionization and crystallization. We demonstrate that this splitting is incomplete for fermionized bosons and restricted by the confinement potential. This incomplete splitting is a consequence of the convergence of the energy in the limit of infinite repulsion and is in agreement with complementary results that we obtain for fermions using MCTDHF. For crystalline bosons, in contrast, the splitting is complete: the interaction energy is capable of overcoming the confinement potential. Our results suggest that the spreading of the density as a function of the dipolar interaction strength diverges as a power law. We describe how to distinguish fermionization from crystallization experimentally from measurements of the one- and two-body density.
Fermionization is what happens to the state of strongly interacting repulsive bosons interacting with contact interactions in one spatial dimension. Crystallization is what happens for sufficiently strongly interacting repulsive bosons with dipolar interactions in one spatial dimension. Crystallization and fermionization resemble each other: in both cases - due to their repulsion - the bosons try to minimize their spatial overlap. We trace these two hallmark phases of strongly correlated one-dimensional bosonic systems by exploring their ground state properties using the one- and two-body density matrix. We solve the N-body Schrödinger equation accurately and from first principles using the multiconfigurational time-dependent Hartree for bosons (MCTDHB) and for fermions (MCTDHF) methods. Using the one- and two-body density, fermionization can be distinguished from crystallization in position space. For N interacting bosons, a splitting into an N-fold pattern in the one-body and two-body density is a unique feature of both, fermionization and crystallization. We demonstrate that this splitting is incomplete for fermionized bosons and restricted by the confinement potential. This incomplete splitting is a consequence of the convergence of the energy in the limit of infinite repulsion and is in agreement with complementary results that we obtain for fermions using MCTDHF. For crystalline bosons, in contrast, the splitting is complete: the interaction energy is capable of overcoming the confinement potential. Our results suggest that the spreading of the density as a function of the dipolar interaction strength diverges as a power law. We describe how to distinguish fermionization from crystallization experimentally from measurements of the one- and two-body density.
Fermionization is what happens when repulsive interacting bosons in one spatial dimension feel infinite interaction. Due to their strong repulsion the bosons try to minimize their spatial overlap like the fermions. Fermionized bosons and noninteracting fermions resemble each other and we trace these similarities of strongly correlated one-dimensional bosonic systems and fermionic system by exploring their ground state properties using the one- and two-body density matrix in coordinate space. We solve the N-body Schrödinger equation accurately and from first principles using the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method. Using the one- and two-body density in momentum space, fermionized bosons can be distinguished from noninteracting fermions.
In this thesis, we study the physical properties of several ultracold few-boson systems depending on the interactions between their constituents. Nowadays, experimentally, it is possible to have great control with high precision over the geometry and the interactions between the particles, making them an excellent setup to test directly the principles of quantum mechanics. A very interesting point is to study the evolution of their properties with the number of particles. The theoretical study of these systems pretends to microscopically understand the current experimental results and give support to new experimental developments.\n\n\t\t\t\t The method that will be used is the exact diagonalization of the Hamiltonian of the system. As we will see, in spite of the attempts to improve it, the method is limited by the fact that, in practice, it is only useful to study few-particle systems. The method has several advantages. First of all, one has access to both the ground and the excited states. In second place, the method is variational and converges to the exact solution as long as the Hilbert space in which we diagonalize is enlarged.\n\t\t\t\t Moreover, since we have access to the states of the system, it is possible to calculate any observable quantity of interest.\n\n\t\t\t\t First, we will study a system of spinless bosons trapped in a two-dimensional harmonic potential. The effect of the trap is to keep the system bound. It will be seen how the presence of a repulsive interaction changes the energy spectrum and other properties of the system. For instance, the density profile, which is usually measurable, and also the two-body distribution function, which is intimately related to the existence of correlations.\n\n\t\t\t\t Afterwards, the focus will be on the particular case of having only two bosons in the system interacting through a strong repulsive force. Inspired by the one-dimensional case where the fermionization phenomenon takes place in the strongly-interacting limit, we will study whether in two dimensions there is a resembling reminiscent effect. In other words, we will analyze if there are properties of the two strongly-interacting bosons in two dimensions that are like the ones of two noninteracting fermions.\n\n\t\t\t\t After that, we will tackle the localization phenomenon in a one-dimensional system that is caused by an external speckle potential that introduces disorder in the system. We will show that the localization is a robust phenomenon against repulsive contact interactions.\n\n\t\t\t\t Finally, we will study the influence of the spin-orbit coupling in a system of bosons with two possible pseudospin components, associated, for instance, to two hyperfine levels, confined in a two-dimensional harmonic trap. We will present an exhaustive analysis of the combined effects of the interaction and the spin-orbit coupling in the spectrum and the properties of the system. In particular we show the existence of a crossover in the ground state of the system susceptible to be experimentally identified.
Quasiperiodic potentials and dipolar interactions each impose long-range order in quantum systems, but their interplay unlocks a rich landscape of unexplored quantum phases. In this work, we investigate how dipolar bosonic crystals respond to correlated disorder in the form of quasiperiodic potentials. Using exact numerical simulations and a suite of observables—including order parameters, energy, density distributions, and two-body coherence measures—we explore one-dimensional dipolar bosons in quasiperiodic lattices at both commensurate and incommensurate fillings. Our results reveal a complex competition between superfluid, Mott-insulator, density-wave, and crystalline phases, governed by the intricate balance of dipolar interactions, kinetic energy, and disorder strength. Crucially, we identify mechanisms that influence dipolar crystals, showing their surprising robustness even in the presence of strong quasiperiodic disorder. Strikingly, we challenge previous claims by demonstrating that a crystal phase—expected to precede full crystallization—does not emerge in the ground state. Instead, its traits appear only under moderate disorder, but never fully develop, giving way to a direct transition from a charge density wave to a crystal state. These findings provide new insights into the resilience of many-body quantum phases in complex environments and pave the way for engineering exotic quantum states in ultracold atomic systems.
Many-body densities and correlation functions are of paramount importance for understanding quantum many-body physics. Here, we present a method to compute them; our approach is general and based on the action of bosonic or fermionic annihilation field operators on the many-body wavefunction. We analyze N=6 quasi-one-dimensional harmonically-trapped bosons with weak to strong contact interaction strength up to the Tonks-Girardeau limit with infinite repulsion using the MultiConfigurational Time-Dependent Hartree method for indistinguishable particles (MCTDH-X). We compare our MCTDH-X solutions to the analytical ones in the infinite repulsion regime as well as to the so-called correlated pair wavefunction approach and find a good agreement. Since numerical approximations are not bound to the cases where analytical solutions are known, we thus demonstrate a general method to investigate high-order reduced density matrices and correlation functions in systems for which analytical solutions are unknown. We trace the build-up of correlation features in the crossover from weak interactions to the Tonks-Girardeau limit and find that the higher-order correlation functions and densities resemble those in the Tonks-Girardeau limit for way smaller interactions than anticipated from just the one-body density.
We explore the consequences of disorder on phase coherence in the Mott insulator phases in an optical lattice. Few bosons with contact interaction in small optical lattice can feature varieties of insulating phases---weakly interacting Mott in deep lattice, maximally fragmented and strongly interacting Mott in intermediate lattice, weak Mott with double filling and intrawell coherence, fermionized Mott with strong intrawell coherence. Utilization of the multiconfigurational time-dependent Hartree method for bosons to solve the many-boson Schr\"odinger equation, facilitates to understand the microscopic effect of disorder on the different kinds of Mott phases in the primary lattice. The many-body properties are analyzed by distinct measures of the reduced one-body density in real and momentum space, fragmentation, order parameter, variance of spatial single-shot measurements, compressibility, and the Glauber normalized correlation functions. We find very complex competition of localization due to disorder and Mott correlation. We observe distinct response of four different Mott phases in the disordered lattice. When the weakly interacting Mott exhibits the Bose glass phase, the strongly correlated and fully fragmented Mott exhibits leakage, melting and central localization. In contrast, weak Mott with double filling which is a spatially separated dimer in each site, exhibits only some dissociation of intrawell coherence and assist in the development of interwell coherence in the presence of strong disorder. For the fermionized Mott, when the density in each well is fragmented, strong disorder interferes with the intrawell correlation and the characteristic dip in each site starts to disappear leading to simple Mott localization with a pair of bosons.
Strongly interacting finite ensembles of dipolar bosons in commensurately filled one-dimensional optical lattices exhibit diverse quantum phases that are rich in physics. As the strength of the long-range boson–boson interaction increases, the system transitions across different phases: from a superfluid, through a Mott-insulator and a Tonks–Girardeau gas to a crystal state. The signature of these phases and their transitions can be unequivocally identified by an experimentally detectable order parameter, recently described in Phys. Rev. A 98 235301 (2018 [ 33 ]). Herein, we calculate the momentum distributions and the normalized Glauber correlation functions of dipolar bosons in a one-dimensional optical lattice in order to characterize all their phases. To understand the behavior of the correlations across the phase transitions, we first investigate the eigenfunctions and eigenvalues of the one-body reduced density matrix as a function of the dipolar interaction strength. We then analyze the real- and momentum-space Glauber correlation functions, thereby gaining a spatially and momentum-resolved insight into the coherence properties of these quantum phases. We find an intriguing structure of non-local correlations that, independently of other observables, reveal the phase transitions of the system. In particular, spatial localization and momentum delocalization accompany the formation of correlated islands in the density as interactions become stronger. Our study showcases that precise control of intersite correlations is possible through the manipulation of the depth of the lattice, while intrasite correlations can be influenced by changing the dipolar interaction strength.