Search

We revisit the dynamics of a Bose-Einstein condensate in a double-well potential, from the regime of Josephson plasma oscillations to the self-trapping regime, by means of the Bogoliubov quasiparticle projection method. For very small imbalance between left and right wells only the lowest Bogoliubov mode is significantly occupied. In this regime the system performs plasma oscillations at the corresponding frequency, and the evolution of the condensate is characterized by a periodic transfer of population between the ground and the first excited state. As the initial imbalance is increased, more excited modes -- though initially not macroscopically occupied -- get coupled during the evolution of the system. Since their population also varies with time, the frequency spectrum of the imbalance turns out to be still peaked around a single frequency, which is continuously shifted towards lower values. The nonlinear mixing between Bogoliubov modes eventually drives the system into the the self-trapping regime, when the population of the ground state can be transferred completely to the excited states at some time during the evolution. For simplicity, here we consider a one-dimensional setup, but the results are expected to hold also in higher dimensions., Comment: 7 pages, 5 figures

We investigate the correspondence between the tight-binding Floquet Hamiltonian of a periodically modulated honeycomb lattice and the Haldane model. We show that - though the two systems share the same topological phase diagram, as reported in a breakthrough experiment with ultracold atoms in a stretched honeycomb lattice [Jotzu et al., Nature 515, 237 (2014)] - the corresponding Hamiltonians are not equivalent, the one of the shaken lattice presenting a much richer structure., Comment: 8 pages, 7 figures

Tight-binding models for ultracold atoms in optical lattices can be properly defined by using the concept of maximally localized Wannier functions for composite bands. The basic principles of this approach are reviewed here, along with different applications to lattice potentials with two minima per unit cell, in one and two spatial dimensions. Two independent methods for computing the tight-binding coefficients - one ab initio, based on the maximally localized Wannier functions, the other through analytic expressions in terms of the energy spectrum - are considered. In the one dimensional case, where the tight-binding coefficients can be obtained by designing a specific gauge transformation, we consider both the case of quasi resonance between the two lowest bands, and that between s and p orbitals. In the latter case, the role of the Wannier functions in the derivation of an effective Dirac equation is also reviewed. Then, we consider the case of a two dimensional honeycomb potential, with particular emphasis on the Haldane model, its phase diagram, and the breakdown of the Peierls substitution. Tunable honeycomb lattices, characterized by movable Dirac points, are also considered. Finally, general considerations for dealing with the interaction terms are presented., Comment: minor corrections, references added; 24 pages, 32 figures; to appear on SCIENCE CHINA Physics, Mechanics & Astronomy

We present an ab initio analysis of a continuous Hamiltonian that maps into the celebrated Haldane model. The tunnelling coefficients of the tight-binding model are computed by means of two independent methods - one based on the maximally localized Wannier functions, the other through analytic expressions in terms of gauge-invariant properties of the spectrum - that provide a remarkable agreement and allow to accurately reproduce the exact spectrum of the continuous Hamiltonian. By combining these results with the numerical calculation of the Chern number, we are able to draw the phase diagram in terms of the physical parameters of the microscopic model. Remarkably, we find that only a small fraction of the original phase diagram of the Haldane model can be accessed, and that the topological insulator phase is suppressed in the deep tight-binding regime., Comment: 11 pages, 9 figures

The physics of one dimensional optical superlattices with resonant $s$-$p$ orbitals is reexamined in the language of appropriate Wannier functions. It is shown that details of the tight binding model realized in different optical potentials crucially depend on the proper determination of Wannier functions. We discuss the properties of a superlattice model which quasi resonantly couples $s$ and $p$ orbitals and show its relation with different tight binding models used in other works., Comment: 9pp, 10 figures, updated references, comments to kuba@if.uj.edu.pl

We present two independent calculations of the tight-binding parameters for a specific realization of the Haldane model with ultracold atoms. The tunneling coefficients up to next-to-nearest neighbors are computed ab-initio by using the maximally localized Wannier functions, and compared to analytical expressions written in terms of gauge invariant, measurable properties of the spectrum. The two approaches present a remarkable agreement and evidence the breakdown of the Peierls substitution: (i) the phase acquired by the next-to-nearest tunneling amplitude $t_{1}$ presents quantitative and qualitative differences with respect to that obtained by the integral of the vector field A, as considered in the Peierls substitution, even in the regime of low amplitudes of A; (ii) for larger values, also $|t_{1}|$ and the nearest-neighbor tunneling $t_{0}$ have a marked dependence on A. The origin of this behavior and its implications are discussed., Comment: 5 pages, 5 figures

We review the derivation of the effective Dirac equation for ultracold atoms in one-dimensional bichromatic optical lattices, following the proposal by Witthaut et al. Phys. Rev. A 84, 033601 (2011). We discuss how such a derivation - based on a suitable rotation of the Bloch basis and on a coarse graining approximation - is affected by the choice of the Wannier functions entering the coarsening procedure. We show that in general the Wannier functions obtained by rotating the maximally localized Wannier functions for the original Bloch bands can be sufficiently localized for justifying the coarse graining approximation. We also comment on the relation between the rotation needed to achieve the Dirac form and the standard Foldy-Wouthuysen transformation. Our results provide a solid ground for the interpretation of the experimental results by Salger et al. Phys. Rev. Lett. 107, 240401 (2011) in terms of an effective Dirac dynamics., Comment: 6 pages, 5 figures

We present an accurate ab initio tight-binding model, capable of describing the dynamics of Dirac points in tunable honeycomb optical lattices following a recent experimental realization [L. Tarruell et al., Nature 483, 302 (2012)]. Our scheme is based on first-principle maximally localized Wannier functions for composite bands. The tunneling coefficients are calculated for different lattice configurations, and the spectrum properties are well reproduced with high accuracy. In particular, we show which tight binding description is needed in order to accurately reproduce the position of Dirac points and the dispersion law close to their merging, for different laser intensities., Comment: 11 pages, 16 figures

We discuss how to construct tight-binding models for ultra cold atoms in honeycomb potentials, by means of the maximally localized Wannier functions (MLWFs) for composite bands introduced by Marzari and Vanderbilt [1]. In particular, we work out the model with up to third-nearest neighbors, and provide explicit calculations of the MLWFs and of the tunneling coefficients for the graphene-lyke potential with two degenerate minima per unit cell. Finally, we discuss the degree of accuracy in reproducing the exact Bloch spectrum of different tight-binding approximations, in a range of typical experimental parameters., Comment: 5 pages, 5figures

We discuss a method for constructing generalized Wannier functions that are maximally localized at the minima of a one-dimensional periodic potential with a double-well per unit cell. By following the approach of (Marzari M and Vanderbilt D 1997 Phys. Rev. B 56, 12847), we consider a set of band-mixing Wannier functions with minimal spread, and design a specific two-step gauge transformation of the Bloch functions for a composite two band system. This method is suited to efficiently compute the tight-binding coefficients needed for mapping the continuous system to a discrete lattice model. Their behaviour is analyzed here as a function of the symmetry properties of the double-well (including the possibility of parity-breaking), in a range of feasible experimental parameters., Comment: 21 pages, 10 figures; revised version; corrected Figs. 2 and 8; added references

We investigate the dynamics of a wave packet in a parity-breaking one-dimensional periodic potential slowly varied in time and perturbed by a linear potential. Parity is broken by considering an asymmetric double well per unit cell. By comparing the prediction of the semiclassical dynamics with the full Schr\"odinger solution, we show that Bloch oscillations are strongly affected by anomalous velocity corrections related to Berry's phase. We characterize how these effects depend on the degree of parity breaking of the potential and on the modulation parameters. We also discuss how to measure the effects of the anomalous velocity in current experiments with non-interacting Bose-Einstein condensates in bichromatic optical lattices, under the effect of gravity., Comment: 8 pages, 10 figures

It is shown that the nonlinear wave equation $\partial_t^2\phi - \partial^2_x \phi -\mu_0\partial_x(\partial_x\phi)^3 =0$, which is the continuum limit of the Fermi-Pasta-Ulam (FPU) beta model, has a positive Lyapunov exponent lambda_1, whose analytic energy dependence is given. The result (a first example for field equations) is achieved by evaluating the lattice-spacing dependence of lambda_1 for the FPU model within the framework of a Riemannian description of Hamiltonian chaos. We also discuss a difficulty of the statistical mechanical treatment of this classical field system, which is absent in the dynamical description., Comment: 4 pages, 1 figure