We investigate the quantum many-body dynamics of bosonic atoms hopping in a two-leg ladder with strong on-site contact interactions. We observe that when the atoms are prepared in a staggered pattern with pairs of atoms on every other rung, singlon defects, i.e. rungs with only one atom, can localize due to an emergent topological model, even though the underlying model in the absence of interactions admits only topologically trivial states. This emergent topological localization results from the formation of a zero-energy edge mode in an effective lattice formed by two adjacent chains with alternating strong and weak hoping links (Su-Schrieffer-Heeger chains) and opposite staggering which interface at the defect position. Our findings open the opportunity to dynamically generate non-trivial topological behaviors without the need for complex Hamiltonian engineering.

Random defects do not constitute the unique source of electron localization in two dimensions. Lattice quasidisorder generated from two inplane superimposed rotated, main and secondary, square lattices, namely monolayers where moiré patterns are formed, leads to a sharp localized to delocalized single‐particle transition. This is demonstrated here for both, discrete and continuum models of moiré patterns that arise as the twisting angle between the main and the secondary lattices is varied in the interval . Localized to delocalized transition is recognized as the moiré patterns depart from being perfect square crystals to non‐crystalline structures. Extended single‐particle states are found for rotation angles associated with Pythagorean triples that produce perfectly periodic structures. Conversely, angles not arising from such Pythagorean triples lead to non‐commensurate or quasidisordered structures, thus originating localized states. These conclusions are drawn from a stationary analysis where the standard inverse participation ratio (IPR) parameter measuring localization allowed to detect the transition. While both, ground state and excited states are analyzed for the discrete model, where the secondary lattice is considered as a perturbation of the main one, the sharp transition is tracked back for the fundamental state in the continuous scenario where the secondary lattice is not a perturbation any more.

Ultracold dipolar hard-core bosons in optical ladders provide exciting possibilities for the quantum simulation of anisotropic spin ladders. We show that introducing a tilt along the rungs results in a rich phase diagram at unit rung filling. In particular, for a sufficiently strong dipolar strength, the interplay between the long-range tail of the dipolar interactions and the tilting leads to the emergence of a quantum floating phase, a critical phase with incommensurate density-density correlations. Interestingly, the floating phase is topological, constituting an intermediate gapless stage in the melting of a crystal into a gapped topological Haldane phase. This novel scenario for topological floating phases in dipolar spin ladders can be investigated in ongoing experiments.

Ultracold dipolar particles pinned in optical lattices or tweezers provide an excellent platform for studying out-of-equilibrium quantum magnetism with dipole-mediated couplings. Starting with an initial state with spins of opposite orientation in each of the legs of a ladder lattice, we show that spin relaxation displays an unexpected dependence on inter-leg distance and dipole orientation. This intricate dependence, stemming from the interplay between intra- and inter-leg interactions, results in three distinct dynamical regimes: (i) ergodic, characterized by the fast relaxation towards equilibrium of correlated pairs of excitations generated at exponentially fast rates from the initial state; (ii) metastable, in which the state is quasi-localized in the initial state and only decays at exceedingly long timescales, resembling false vacuum decay; and, surprisingly, (iii) partially-relaxed, with coexisting fast partial relaxation and very long-lived partial quasi-localization. Realizing these intriguing dynamics is within reach of current state-of-the-art experiments in dipolar gases.

Quasiparticles and their interactions are a key part of our understanding of quantum many-body systems. Quantum simulation experiments with cold atoms have in recent years advanced our understanding of isolated quasiparticles, but so far they have provided limited information regarding their interactions and possible bound states. Here, we show how exploring mobile impurities immersed in a Bose-Einstein condensate (BEC) in a two-dimensional lattice can address this problem. First, the spectral properties of individual impurities are examined, and in addition to the attractive and repulsive polarons known from continuum gases, we identify a new kind of quasiparticle stable for repulsive boson-impurity interactions. The spatial properties of polarons are calculated showing that there is an increased density of bosons at the site of the impurity both for repulsive and attractive interactions. We then derive an effective Schrödinger equation describing two polarons interacting via the exchange of density oscillations in the BEC, which takes into account strong impurity-boson two-body correlations. Using this, we show that the attractive nature of the effective interaction between two polarons combined with the two-dimensionality of the lattice leads to the formation of bound states -- i.e. bipolarons. The wave functions of the bipolarons are examined showing that the ground state is symmetric under particle exchange and therefore relevant for bosonic impurities, whereas the first excited state is doubly degenerate and odd under particle exchange making it relevant for fermionic impurities. Our results show that quantum gas microscopy in optical lattices is a promising platform to explore the spatial properties of polarons as well as to finally observe the elusive bipolarons.

The use of machine learning algorithms to address classification problems in several scientific branches has increased over the past years. In particular, the supervised learning technique with artificial neural networks has been successfully employed in classifying phases of matter. In this article, we use a fully connected feed-forward neural network to classify extended and localized single-particle states that arise from quasiperiodic one-dimensional lattices. We demonstrate that our neural network achieves to correctly uncover the nature of the single-particle states even when the wave functions come from a more complex Hamiltonian than the one used to train the network.

Polar lattice gases are usually assumed to have an intersite interaction that decays with the interparticle distance r as 1/r3. However, a loose-enough transversal confinement may strongly modify the dipolar decay in one-dimensional lattices. We show that this modification alters significantly the ground-state properties of hard-core dipolar bosons. For repulsive intersite interactions, the corrected decay alters the conditions for devil's staircase insulators, affecting significantly the particle distribution in the presence of an overall harmonic confinement. For attractive interactions, it results in a reduction of the critical dipole interaction for the formation of self-bound clusters, and for a marked enhancement of the region of liquefied lattice droplets.

We show that the long-distance behavior of the two-body density correlation functions and the Cooper-pair probability density of a balanced mixture of a two-component Fermi gas at T = 0, is universal along the BEC-BCS crossover. Our result is demonstrated by numerically solving the mean-field BCS model for different finite short-range atomic interaction potentials. We find an analytic expression for the correlation length in terms of the chemical potential and the energy gap at zero momentum.

We study the temporal growth and spatial propagation of quantum correlations in a two-dimensional bilayer realising a spin-1/2 quantum XXZ model with couplings mediated by long-range and anisotropic dipolar interactions. Starting with an initial state consisting of spins with opposite magnetization in each of the layers, we predict the emergence of a momentum-dependent dynamic instability in the spin structure factor that results, at short times, in the creation of pairs of excitations at exponentially fast rates. The created pairs present a characteristic momentum distribution that can be tuned by controlling the dipolar orientation, the layer separation or the dipolar couplings. The predicted behavior remains observable at very low filling fractions, making it accessible in state-of-the-art experiments with Rydberg atoms, magnetic atoms, and polar molecule arrays.

Pair localization in one-dimensional quasicrystals with nearest-neighbor hopping is independent of whether short-range interactions are repulsive or attractive. We numerically demonstrate that this symmetry is broken when the hopping follows a power law 1 /r α . In particular, for repulsively bound states, we find that the critical quasiperiodicity that signals the transition to localization is always bounded by the standard Aubry-Andr´e critical point, whereas attractively bound dimers get localized at larger quasiperiodic modulations when the range of the hopping increases. Extensive numerical calculations establish the contrasting nature of the pair energy gap for repulsive and attractive interactions, as well as the behavior of the algebraic localization of the pairs as a function of quasiperiodicity, interaction strength, and power-law hops. The results here discussed are of direct relevance to the study of the quantum dynamics of systems with power-law couplings.

We study the effects of interparticle interactions and power-law tunneling couplings on quantum walks executed by both a single one and a pair of hard-core bosons moving in clean and disordered one-dimensional lattices. For this purpose, we perform exact diagonalization to explicitly evaluate the short and long time probabilities of finding the walkers within a surveillance area. Our main conclusions, summarized in phase diagrams in the disorder-power-law and interaction-disorder spaces, allowed us to discern two different scenarios for the single and two quantum walkers dynamics. While in the single particle case the transition to localized and extended regimes is identified for well defined values of the disorder amplitude and power law hopping, those frontiers are replaced by diffuse contours in the interacting two particle case. In fact, counterintuitive transport regimes as diffusion enhanced by disorder, and space constrained dynamics assisted by both interactions and short tunneling range are found. Our results are of direct relevance for quantum systems with long-range interactions that are currently realized in the laboratory.

In this manuscript, we discuss the emergence of p-wave superfluidity in a dipolar Fermi gas confined in a double layer array of parallel optical lattices in two dimensions. The dipole moments of the molecules placed at the sites of the optical lattices, separated a distance L and pointing in opposite directions produce an effective attractive interaction among them, except between those dipoles situated one on top of the other. Such interaction between dipoles is precisely the origin of the non-conventional superfluid state. We present the analysis for the ground state of the many-body system within the mean-field scheme. In particular, we study the stable regions, as a function of the system parameters, namely the effective interaction between dipoles and the filling factor n, for which the superfluid state can exist. Following the BKT scheme, we estimate the critical temperature of the superfluid state.

We investigate the persistence, in time and space, of ferromagnetic domains in two dimensions subjected to the influence of both the static disorder of variable strength and weak interactions. The domains are represented by a two-species bosonic mixture of Rb-87 ultracold atoms in different hyperfine states, such that initially one lies on the left half and the other on the right half of a square lattice. The dynamics of the double domain is followed by describing the two-component superfluid through the time-dependent Gross-Pitaevskii coupled equations, with values of the intra- and interspecies interaction that guarantee miscibility of the components. A robust analysis of the magnetization dynamics for several values of the interspecies interaction, reachable in current experimental setups, and the investigation of the density-weighted magnetization correlator lead us to conclude that the presence of structural disorder yields a slowdown the process of destruction of the initial ferromagnetic order. As shown by our numerical calculations, magnetization is maintained up to 50% of its initial value for the largest disorder amplitude considered.

We investigate the relaxation process of ferromagnetic domains in 2D subjected to the influence of both, static disorder of variable strength and weak interactions. The domains are represented by a two species bosonic mixture of $^{87}$Rb ultracold atoms, such that initially each specie lies on left and right halves of a square lattice. The dynamics of the double domain is followed by describing the two-component superfluid, at mean field level, through the time dependent Gross-Pitaevskii coupled equations, considering values of the intra and inter-species interaction, reachable in current experimental setups, that guaranty miscibility of the components. A robust analysis for several values inter-species interaction leads us to conclude that the presence of structural disorder leads to slowdown the relaxation process of the initial ferromagnetic order. As shown by our numerical experiments, magnetization is maintained up to 60 percent of its initial value for the largest disorder amplitude.

We present a thorough pedagogical analysis of the single particle localization phenomenon in a quasiperiodic lattice in one dimension. Description of disorder in the lattice is represented by the Aubry-Andr\'e model. Characterization of localization is performed through the analysis of both, stationary and dynamical properties. The stationary properties investigated are the inverse participation ratio (IPR), the normalized participation ratio (NPR) and the energy spectrum as a function of the disorder strength. As expected, the distinctive Hofstadter pattern is found. Two dynamical quantities allow discerning the localization phenomenon, being the spreading of an initially localized state and the evolution of population imbalance in even and odd sites across the lattice.

The superfluid p=px+ipy phases in an ultracold gas of dipolar Fermi molecules lying in two parallel square lattices in 2D are investigated. As shown by a two-body study, dipole moments oriented in opposite directions in each layer are the key ingredients in our mean-field analysis from which unconventional superfluidity is predicted. The T=0 phase diagram summarizes our findings: stable and metastable superfluid phases appear as a function of both, the dipole-dipole interaction coupling parameter and filling factor. A first-order phase transition, and thus a mixture of superfluid phases at different densities, is revealed from the coexistence curves in the metastable region. The model predicts that these superfluid phases can be observed experimentally at 10 nK in molecules of NaK confined in optical lattices of size a=532 nm. Other routes to reach higher temperatures require the use of subwavelength confinement technique.

By exact numerical solutions of the Gross-Pitaevskii (GP) equation in 3D, we assess the validity of 1D and 2D approximations in the study of Bose-Einstein condensates confined in harmonic trap potentials. Typically, these approximations are performed when one or more of the harmonic frequencies are much greater than the remaining ones, using arguments based on the adiabatic evolution of the initial approximated state. Deviations from the 3D solution are evaluated as a function of both the effective interaction strength and the ratio between the trap frequencies that define the reduced dimension where the condensate is confined. The observables analyzed are both stationary and dynamical character, namely, the chemical potential, the wave function profiles, and the time evolution of the approximated 1D and 2D stationary states, considered as initial states in the 3D GP equation. Our study, besides setting quantitative limits on approximations previously developed, should be useful in actual experimental studies where quasi-1D and quasi-2D conditions are assumed. From a qualitative perspective, 1D and 2D approximations certainly become valid when the anisotropy is large, but in addition, the interaction strength needs to be above a certain threshold.

We investigate the formation of Cooper pairs, bound dimers and the dimer‐dimer elastic scattering of ultracold dipolar Fermi molecules confined in a 2D optical lattice bilayer configuration. While the energy and their associated bound states are determined in a variational way, the correlated two‐molecule pair is addressed as in the original Cooper formulation. We demonstrate that the 2D lattice confinement favors the formation of zero center mass momentum bound states. Regarding the Cooper pairs binding energy, this depends on the molecule populations in each layer. Maximum binding energies occur for non‐zero (zero) pair momentum when the Fermi system is polarized (unpolarized). We find an analytic expression for the dimer‐dimer effective interaction in the deep BEC regime. The present analysis represents a route for addressing the BCS‐BEC crossover in dipolar Fermi gases confined in 2D optical lattices within the current experimental panorama.