Source author record

Jakub Zakrzewski

Jakub Zakrzewski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

54works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

54 published item(s)

preprint2026arXiv

Facets of Many Body Localization

Many-body localization (MBL) appears to be a robust example of ergodicity breaking in many-body interacting systems. Here, we review different aspects of MBL, concentrating on various ways the disorder may be introduced into the system studied. In particular, we consider both the random and quasiperiodic diagonal (i.e., on-site) disorder as well as bond disorder as realized in randomly distributed atoms interacting via long-range interactions. We also review the quantum sun model, which seems to be the ideal, albeit artificial, model exhibiting MBL.

preprint2022arXiv

A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation

A precise understanding of mechanisms governing the dynamics of electrons in atoms and molecules subjected to intense laser fields has a key importance for the description of attosecond processes such as the high-harmonic generation and ionization. From the theoretical point of view, this is still a challenging task, as new approaches to solve the time-dependent Schrödinger equation with both good accuracy and efficiency are still emerging. Until recently, the purely numerical methods of real-time propagation of the wavefunction using finite grids have been frequently and successfully used to capture the electron dynamics in small one- or two-electron systems. However, as the main focus of attoscience shifts toward many-electron systems, such techniques are no longer effective and need to be replaced by more approximate but computationally efficient ones. In this paper, we explore the increasingly popular method of expanding the wavefunction of the examined system into a linear combination of atomic orbitals and present a novel systematic scheme for constructing an optimal Gaussian basis set suitable for the description of excited and continuum atomic or molecular states. We analyze the performance of the proposed basis sets by carrying out a series of time-dependent configuration interaction calculations for the hydrogen atom in fields of intensity varying from $5 \times 10^{13}\:\rm W/cm^2$ to $5 \times 10^{14}\:\rm W/cm^2$. We also compare the results with the data obtained using Gaussian basis sets proposed previously by other authors.

preprint2022arXiv

Bond order via cavity-mediated interactions

We numerically study the phase diagram of bosons tightly trapped in the lowest band of an optical lattice and dispersively coupled to a single-mode cavity field. The dynamics is encompassed by an extended Bose-Hubbard model. Here, the cavity-mediated interactions are described by a two-body potential term with a global range and by a correlated tunneling term where the hopping amplitude depends on a global observable. We determine the ground state properties in one dimension by means of the density matrix renormalization group algorithm, focusing on the effects due to the correlated tunneling. The latter is responsible for the onset of bond orders, manifesting in one insulating and two gapless bond ordered phases. We discuss the resulting phases for different geometries that correspond to different relative strengths of the correlated tunneling coefficient. We finally analyze the scaling of the entanglement entropy in the gapless bond ordered phases that appear entirely due to global interactions and determine the corresponding central charges.

preprint2022arXiv

Challenges to observation of many-body localization

We study time dynamics of 1D disordered Heisenberg spin-1/2 chain focusing on a regime of large system sizes and a long time evolution. This regime is relevant for observation of many-body localization (MBL), a phenomenon that is expected to freeze the dynamics of the system and prevent it from reaching thermal equilibrium. Performing extensive numerical simulations of the imbalance, a quantity often employed in the experimental studies of MBL, we show that the regime of a slow power-law decay of imbalance persists to disorder strengths exceeding by at least a factor of 2 the current estimates of the critical disorder strength for MBL. Even though we investigate time evolution up to few thousands tunneling times, we observe no signs of the saturation of imbalance that would suggest freezing of system dynamics and provide a smoking gun evidence of MBL. We demonstrate that the situation is qualitatively different when the disorder is replaced by a quasiperiodic potential. In this case, we observe an emergence of a pattern of oscillations of the imbalance that is stable with respect to changes in the system size. This suggests that the dynamics of quasiperiodic systems remain fully local at the longest time scales we reach provided that the quasiperiodic potential is sufficiently strong. Our study identifies challenges in an unequivocal experimental observation of the phenomenon of MBL.

preprint2022arXiv

Devil's staircase of topological Peierls insulators and Peierls supersolids

We consider a mixture of ultracold bosonic atoms on a one-dimensional lattice described by the XXZ-Bose-Hubbard model, where the tunneling of one species depends on the spin state of a second deeply trapped species. We show how the inclusion of antiferromagnetic interactions among the spin degrees of freedom generates a Devil's staircase of symmetry-protected topological phases for a wide parameter regime via a bosonic analog of the Peierls mechanism in electron-phonon systems. These topological Peierls insulators are examples of symmetry-breaking topological phases, where long-range order due to spontaneous symmetry breaking coexists with topological properties such as fractionalized edge states. Moreover, we identify a region of supersolid phases that do not require long-range interactions. They appear instead due to a Peierls incommensurability mechanism, where competing orders modify the underlying crystalline structure of Peierls insulators, becoming superfluid. Our work show the possibilities that ultracold atomic systems offer to investigate strongly-correlated topological phenomena beyond those found in natural materials.

preprint2022arXiv

Many-body localization regime for cavity induced long-range interacting models

Many-body localization (MBL) features are studied here for a large spin chain model with long range interactions. The model corresponds to cold atoms placed inside a cavity and driven by an external laser field with long range interactions coming from rescattering of cavity photons. Earlier studies were limited to small sizes amenable to exact diagonalization. It is shown that nonergodic features and MBL may exist in this model for random disorder as well as in the presence of tilted potential on experimental time scales also for experimentally relevant system sizes using tensor networks algorithms.

preprint2022arXiv

Phase diagram of 1+1D Abelian-Higgs model and its critical point

We determine the phase diagram of the Abelian-Higgs model in one spatial dimension and time (1+1D) on a lattice. We identify a line of first order phase transitions separating the Higgs region from the confined one. This line terminates in a quantum critical point above which the two regions are connected by a smooth crossover. We analyze the critical point and find compelling evidences for its description as the product of two non-interacting systems, a massless free fermion and a massless free boson. However, we find also some surprizing results that cannot be explained by our simple picture, suggesting this newly discovered critical point to be an unusual one.

preprint2022arXiv

Quantum boomerang effect in systems without time reversal symmetry

In an Anderson localized system, a quantum particle with a nonzero initial velocity returns, on average, to its origin. This recently discovered behavior is known as the quantum boomerang effect. Time reversal invariance was initially thought to be a necessary condition for the existence of this phenomenon. We theoretically analyze the impact of the symmetry breaking on the phenomenon using a one-dimensional system with a spin-orbit coupling and show that the time reversal invariance is not necessary for the boomerang effect to occur.

preprint2022arXiv

Quantum phases of dipolar bosons in one-dimensional optical lattices

We theoretically analyze the phase diagram of a quantum gas of bosons that interact via repulsive dipolar interactions. The bosons are tightly confined by an optical lattice in a quasi one-dimensional geometry. In the single-band approximation, their dynamics is described by an extended Bose-Hubbard model where the relevant contributions of the dipolar interactions consist of density-density repulsion and correlated tunneling terms. We evaluate the phase diagram for unit density using numerical techniques based on the density-matrix renormalization group algorithm. Our results predict that correlated tunneling can significantly modify the parameter range of the topological insulator phase. At vanishing values of the onsite interactions, moreover, correlated tunneling promotes the onset of a phase with a large number of low energy metastable configurations.

preprint2022arXiv

Scar States in Deconfined $\mathbb{Z}_2$ Lattice Gauge Theories

The weak ergodicity breaking induced by quantum many-body scars (QMBS) represents an intriguing concept that has received great attention in recent years due to its relation to unusual non-equilibrium behaviour. Here we reveal that this phenomenon can occur in a previously unexplored regime of a lattice gauge theory, where QMBS emerge due to the presence of an extensive number of local constraints. In particular, by analyzing the gauged Kitaev model, we provide an example where QMBS appear in a regime where charges are deconfined. By means of both numerical and analytical approaches, we find a variety of scarred states far away from the regime where the model is integrable. The presence of these states is revealed both by tracing them directly from the analytically reachable limit, as well as by quantum quenches showing persistent oscillations for specific initial states.

preprint2022arXiv

Unsupervised detection of decoupled subspaces: many-body scars and beyond

Highly excited eigenstates of quantum many-body systems are typically featureless thermal states. Some systems, however, possess a small number of special, low-entanglement eigenstates known as quantum scars. We introduce a quantum-inspired machine learning platform based on a Quantum Variational Autoencoder (QVAE) that detects families of scar states in spectra of many-body systems. Unlike a classical autoencoder, QVAE performs a parametrized unitary operation, allowing us to compress a single eigenstate into a smaller number of qubits. We demonstrate that the autoencoder trained on a scar state is able to detect the whole family of scar states sharing common features with the input state. We identify families of quantum many-body scars in the PXP model beyond the $\mathbb{Z}_2$ and $\mathbb{Z}_3$ families and find dynamically decoupled subspaces in the Hilbert space of disordered, interacting spin ladder model. The possibility of an automatic detection of subspaces of scar states opens new pathways in studies of models with a weak breakdown of ergodicity and fragmented Hilbert spaces.

preprint2021arXiv

Cold atoms meet lattice gauge theory

The central idea of this review is to consider quantum field theory models relevant for particle physics and replace the fermionic matter in these models by a bosonic one. This is mostly motivated by the fact that bosons are more ``accessible'' and easier to manipulate for experimentalists, but this ``substitution'' also leads to new physics and novel phenomena. It allows us to gain new information about among other things confinement and the dynamics of the deconfinement transition. We will thus consider bosons in dynamical lattices corresponding to the bosonic Schwinger or Z$_2$ Bose-Hubbard models. Another central idea of this review concerns atomic simulators of paradigmatic models of particle physics theory such as the Creutz-Hubbard ladder, or Gross-Neveu-Wilson and Wilson-Hubbard models. Finally, we will briefly describe our efforts to design experimentally friendly simulators of these and other models relevant for particle physics.

preprint2021arXiv

Many-body localization with synthetic gauge fields in disordered Hubbard chains

We analyze the localization properties of the disordered Hubbard model in the presence of a synthetic magnetic field. An analysis of level spacing ratio shows a clear transition from ergodic to many-body localized phase. The transition shifts to larger disorder strengths with increasing magnetic flux. Study of dynamics of local correlations and entanglement entropy indicates that charge excitations remain localized whereas spin degree of freedom gets delocalized in the presence of the synthetic flux. This residual ergodicity is enhanced by the presence of the magnetic field with dynamical observables suggesting incomplete localization at large disorder strengths. Furthermore, we examine the effect of quantum statistics on the local correlations and show that the long-time spin oscillations of a hard-core boson system are destroyed as opposed to the fermionic case.

preprint2021arXiv

Nonergodic dynamics in disorder-free potentials

We review the dynamics of interacting particles in disorder-free potentials concentrating on a combination of a harmonic binding with a constant tilt. We show that a simple picture of an effective local tilt describes a variety of cases. Our examples include spinless fermions (as modeled by Heisenberg spin chain in a magnetic field), spinful fermions as well as bosons that enjoy a larger local on-site Hilbert space. We also discuss the domain-wall dynamics that reveals nonergodic features even for a relatively weak tilt as suggested by Doggen et. al. [arXiv:2012.13722]. By adding a harmonic potential on top of the static field we confirm that the surprizing regular dynamics is not entirely due to Hilbert space shuttering. It seems better explained by the inhibited transport within the domains of identically oriented spins. Once the spin-1/2 restrictions are lifted as, e.g., for bosons, the dynamics involve stronger entanglement generation. Again for domain wall melting, the effect of the harmonic potential is shown to lead mainly to an effective local tilt.

preprint2020arXiv

Ab initio study of time-dependent dynamics in strong-field triple ionization

An ab initio analysis of strong-field three-electron ionization in a restricted-dimensionality model reveals the dynamics of the ionization process and the dominant channels for double (DI) and triple ionization (TI). Simulations using wave functions that respect the Pauli principle show that the most likely channel is a sequence of single ionization (SI) and DI, while direct TI has a much lower probability. The dominant DI process has the highest probability for a singlet of up- and down-spin electrons. The results demonstrate the significance of the Pauli principle for the selection of dominant pathways in ionization and possibly other many-electron processes in strong fields.

preprint2020arXiv

Fidelity susceptibility in Gaussian Random Ensembles

The fidelity susceptibility measures sensitivity of eigenstates to a change of an external parameter. It has been fruitfully used to pin down quantum phase transitions when applied to ground states (with extensions to thermal states). Here we propose to use the fidelity susceptibility as a useful dimensionless measure for complex quantum systems. We find analytically the fidelity susceptibility distributions for Gaussian orthogonal and unitary universality classes for arbitrary system size. The results are verified by a comparison with numerical data.

preprint2020arXiv

Many-body localization in Bose-Hubbard model: evidence for the mobility edge

Motivated by recent experiments on interacting bosons in quasi-one-dimensional optical lattice [Nature {\bf 573}, 385 (2019)] we analyse theoretically properties of the system in the crossover between delocalized and localized regimes. Comparison of time dynamics for uniform and density wave like initial states enables demonstration of the existence of the mobility edge. To this end we define a new observable, the mean speed of transport at long times. It gives us an efficient estimate of the critical disorder for the crossover. We also show that the mean velocity growth of occupation fluctuations close to the edges of the system carries the similar information. Using the quantum quench procedure we show that it is possible to probe the mobility edge for different energies.

preprint2020arXiv

Quantum boomerang effect for interacting particles

When a quantum particle is launched with a finite velocity in a disordered potential, it may surprisingly come back to its initial position at long times and remain there forever. This phenomenon, dubbed ``quantum boomerang effect'', was introduced in [Phys. Rev. A 99, 023629 (2019)]. Interactions between particles, treated within the mean-field approximation, are shown to partially destroy the boomerang effect: the center of mass of the wave packet makes a U-turn, but does not completely come back to its initial position. We show that this phenomenon can be quantitatively interpreted using a single parameter, the average interaction energy.

preprint2020arXiv

Strong-field ionization of atoms with $p^3$ valence shell: Two versus three active electrons

For a model atom with the $p^3$ valence shell we construct consistent three- and two-active electrons models enabling their direct comparison. Within these models, we study the influence of the third active electron on the double ionization yield in strong femtosecond laser fields. We reveal proportionality between double ionization signals obtained with both models in the field intensity region where non-sequential ionization dominates. We derive analytically a correspondence rule connecting the double ionization yields obtained within the three- and two-active electrons models.

preprint2020arXiv

Superfluid phases induced by the dipolar interactions

We determine the quantum ground state of dipolar bosons in a quasi-one-dimensional optical lattice and interacting via $s$-wave scattering. The Hamiltonian is an extended Bose-Hubbard model which includes hopping terms due to the interactions. We identify the parameter regime for which the coefficients of the interaction-induced hopping terms become negative. For these parameters we numerically determine the phase diagram for a canonical ensemble and by means of density matrix renormalization group. We show that at sufficiently large values of the dipolar strength there is a quantum interference between the tunneling due to single-particle effects and the one due to the interactions. Because of this phenomenon, incompressible phases appear at relatively large values of the single-particle tunneling rates. This quantum interference cuts the phase diagram into two different, disconnected superfluid phases. In particular, at vanishing kinetic energy, the phase is always superfluid with a staggered superfluid order parameter. These dynamics emerge from quantum interference phenomena between quantum fluctuations and interactions and shed light into their role in determining the thermodynamic properties of quantum matter.

preprint2020arXiv

Thouless time analysis of Anderson and many-body localization transitions

Spectral statistics of disordered systems encode Thouless and Heisenberg time scales whose ratio determines whether the system is chaotic or localized. Identifying similarities between system size and disorder strength scaling of Thouless time for disordered quantum many-body systems with results for 3D and 5D Anderson models, we argue that the two-parameter scaling breaks down in the vicinity of the transition to the localized phase signalling subdiffusive dynamics.

preprint2020arXiv

Time dynamics with matrix product states: Many-body localization transition of large systems revisited

We compare accuracy of two prime time evolution algorithms involving Matrix Product States - tDMRG (time-dependent density matrix renormalization group) and TDVP (time-dependent variational principle). The latter is supposed to be superior within a limited and fixed auxiliary space dimension. Surprisingly, we find that the performance of algorithms depends on the model considered. In particular, many-body localized systems as well as the crossover regions between localized and delocalized phases are better described by tDMRG, contrary to the delocalized regime where TDVP indeed outperforms tDMRG in terms of accuracy and reliability. As an example, we study many-body localization transition in a large size Heisenberg chain. We discuss drawbacks of previous estimates [Phys. Rev. B 98, 174202 (2018)] of the critical disorder strength for large systems.

preprint2016arXiv

Different lattice geometries with synthetic dimension

The possibility of creating different geometries with the help of an extra synthetic dimension in optical lattices is studied. Additional linear potential and Raman assisted tunnelings are used to engineer well controlled tunnelings between available states. The great flexibility of the system allows us to obtain different geometries of synthetic lattices with possibility of adding synthetic gauge fields.

preprint2016arXiv

Locating the quantum critical point of the Bose-Hubbard model through singularities of simple observables

We show that the critical point of the two-dimensional Bose-Hubbard model can be easily found through studies of either on-site atom number fluctuations or the nearest-neighbor two-point correlation function (the expectation value of the tunnelling operator). Our strategy to locate the critical point is based on the observation that the derivatives of these observables with respect to the parameter that drives the superfluid-Mott insulator transition are singular at the critical point in the thermodynamic limit. Performing the quantum Monte Carlo simulations of the two-dimensional Bose-Hubbard model, we show that this technique leads to the accurate determination of the position of its critical point. Our results can be easily extended to the three-dimensional Bose-Hubbard model and different Hubbard-like models. They provide a simple experimentally-relevant way of locating critical points in various cold atomic lattice systems.

preprint2016arXiv

Topological Rice-Mele model in an emergent lattice: Exact diagonalization approach

Using exact diagonalization methods we study possible phases in a one dimensional model of two differently populated fermionic species in a periodically driven optical lattice. The shaking amplitude and frequency are chosen to resonantly drive $s-p$ transition while minimizing the standard intraband tunnelings. We numerically verify that in the vicinity of vanishing intraband tunnelings the system, for an appropriate filling, shows an emergent density wave configuration of composites. The majority fermions moving in such a lattice mimic the celebrated Rice-Mele model. Far away from that region structure changes to clustered phase, with intermediate phase abundantly populated by defects of the density wave. These defects lead to loclaized modes carrying fractional particle charge. The results obtained are compared with earlier mean field approximation predictions.

preprint2015arXiv

Controlling disorder with periodically modulated interactions

We investigate a celebrated problem of one dimensional tight binding model in the presence of disorder leading to Anderson localization from a novel perspective. A binary disorder is assumed to be created by immobile heavy particles for the motion of the lighter, mobile species in the limit of no interaction between mobile particles. Fast periodic modulations of interspecies interactions allow us to produce an effective model with small diagonal and large off-diagonal disorder unexplored in cold atoms experiments. We present an expression for an approximate Anderson localization length and verify the existence of the well known extended resonant mode and analyze the influence of nonzero next-nearest neighbor hopping terms. We point out that periodic modulation of interaction allow disorder to work as a tunable band-pass filter for momenta.

preprint2015arXiv

Dynamics of heat and mass transport in a quantum insulator

The real time evolution of two pieces of quantum insulators, initially at different temperatures, is studied when they are glued together. Specifically, each subsystem is taken as a Bose-Hubbard model in a Mott insulator state. The process of temperature equilibration via heat transfer is simulated in real time using the Minimally Entangled Typical Thermal States algorithm. The analytic theory based on quasiparticles transport is also given.

preprint2015arXiv

Properties of the one-dimensional Bose-Hubbard model from a high-order perturbative expansion

We employ a high-order perturbative expansion to characterize the ground state of the Mott phase of the one-dimensional Bose-Hubbard model. We compute for different integer filling factors the energy per lattice site, the two-point and density-density correlations, and expectation values of powers of the on-site number operator determining the local atom number fluctuations (variance, skewness, kurtosis). We compare these expansions to numerical simulations of the infinite-size system to determine their range of applicability. We also discuss a new sum rule for the density-density correlations that can be used in both equilibrium and non-equilibrium systems.

preprint2015arXiv

Role of correlations and off-diagonal terms in binary disordered one dimensional systems

We investigate one dimensional tight binding model in the presence of a correlated binary disorder. The disorder is due to the interaction of particles with heavy immobile other species. Off-diagonal disorder is created by means of a fast periodic modulation of interspecies interaction. The method based on transfer matrix techniques allows us to calculate the energies of extended modes in the correlated binary disorder. We focus on $N$-mer correlations and regain known results for the case of purely diagonal disorder. For off-diagonal disorder we find resonant energies. We discuss ambiguous properties of those states and compare analytical results with numerical calculations. Separately we describe a special case of the dual random dimer model.

preprint2015arXiv

Spontaneous magnetization and anomalous Hall effect in an emergent Dice lattice

Ultracold atoms in optical lattices serve as a tool to model different physical phenomena appearing originally in condensed matter. To study magnetic phenomena one needs to engineer synthetic fields as atoms are neutral. Appropriately shaped optical potentials force atoms to mimic charged particles moving in a given field. We present the realization of artificial gauge fields for the observation of anomalous Hall effect. Two species of attractively interacting ultracold fermions are considered to be trapped in a shaken two dimensional triangular lattice. A combination of interaction induced tunneling and shaking can result in an emergent Dice lattice. In such a lattice the staggered synthetic magnetic flux appears and it can be controlled with external parameters. The obtained synthetic fields are non-Abelian. Depending on the tuning of the staggered flux we can obtain either anomalous Hall effect or its quantized version. Our results are reminiscent of Anomalous Hall conductivity in spin-orbit coupled ferromagnets.

preprint2015arXiv

Time reversals of irreversible quantum maps

We introduce the notion of time reversal in open quantum systems as represented by linear quantum operations, and a related generalization of classical entropy production in the environment. This functional is the ratio of the probability to observe a transition between two states under the forward and the time reversed dynamics, and leads, as in the classical case, to fluctuation relations as tautological identities. As in classical dynamics in contact with a heat bath, time reversal is not unique, and we discuss several possibilities. For any bistochastic map its dual map preserves the trace and describes a legitimate dynamics reversed in time, in that case the entropy production in the environment vanishes. For a generic stochastic map we construct a simple quantum operation which can be interpreted as a time reversal. For instance, the decaying channel, which sends the excited state into the ground state with a certain probability, can be reversed into the channel transforming the ground state into the excited state with the same probability.

preprint2014arXiv

Bose-Hubbard model with random impurities: Multiband and nonlinear hopping effects

We investigate the phase diagrams of theoretical models describing bosonic atoms in a lattice in the presence of randomly localized impurities. By including multiband and nonlinear hopping effects we enrich the standard model containing only the chemical-potential disorder with the site-dependent hopping term. We compare the extension of the MI and the BG phase in both models using a combination of the local mean-field method and a Hartree-Fock-like procedure, as well as, the Gutzwiller-ansatz approach. We show analytical argument for the presence of triple points in the phase diagram of the model with chemical-potential disorder. These triple points however, cease to exists after the addition of the hopping disorder.

preprint2014arXiv

Non-standard Hubbard models in optical lattices: a review

Originally, the Hubbard model has been derived for describing the behaviour of strongly-correlated electrons in solids. However, since over a decade now, variations of it are also routinely being implemented with ultracold atoms in optical lattices. We review some of the rich literature on this subject, with a focus on more recent non-standard forms of the Hubbard model. After an introduction to standard (fermionic and bosonic) Hubbard models, we discuss briefly common models for mixtures, as well as the so called extended Bose-Hubbard models, that include interactions between neighboring sites, next-neighboring sites, and so on. The main part of the review discusses the importance of additional terms appearing when refining the tight-binding approximation on the original physical Hamiltonian. Even when restricting the models to the lowest Bloch band is justified, the standard approach neglects the density-induced tunneling (which has the same origin as the usual on-site interaction). The importance of these contributions is discussed for both contact and dipolar interactions. For sufficiently strong interactions, also the effects related to higher Bloch bands become important even for deep optical lattices. Different approaches that aim at incorporating these effects, mainly via dressing the basis Wannier functions with interactions, leading to effective, density-dependent Hubbard-type models, are reviewed. We discuss also examples of Hubbard-like models that explicitly involve higher $p$-orbitals, as well as models that couple dynamically spin and orbital degrees of freedom. Finally, we review mean-field nonlinear-Schrödinger models of the Salerno type that share with the non-standard Hubbard models the nonlinear coupling between the adjacent sites. In that part, discrete solitons are the main subject of the consideration. We conclude by listing some future open problems.

preprint2014arXiv

One-dimensional s-p superlattice

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.

preprint2014arXiv

Rice-Mele model with topological solitons in an optical lattice

Attractive ultra-cold fermions trapped in a one-dimensional periodically shaken opticla lattices are considered. For an appropriate resonant shaking the system realizes paradigmatic dimes physics described by Rice-Mele model. The important feature of our system is the possible presence of controlled defects. They result in the creation of topologically protected loclaized modes carrying fractional particle number. Their possible experimental signatures are discussed.

preprint2013arXiv

Density dependent tunneling in the extended Bose-Hubbard model

Recently, it has become apparent that, when the interactions between polar molecules in optical lattices becomes strong, the conventional description using the extended Hubbard model has to be modified by additional terms, in particular a density-dependent tunneling term. We investigate here the influence of this term on the ground-state phase diagrams of the two dimensional extended Bose-Hubbard model. Using Quantum Monte Carlo simulations, we investigate the changes of the superfluid, supersolid, and phase-separated parameter regions in the phase diagram of the system. By studying the interplay of the density-dependent hopping with the usual on-site interaction U and nearest-neighbor repulsion V, we show that the ground-state phase diagrams differ significantly from the ones that are expected from the standard extended Bose-Hubbard model. However we find no indication of pair-superfluid behavior in this two dimensional square lattice study in contrast to the one-dimensional case.

preprint2013arXiv

Fast dynamics for atoms in optical lattices

Cold atoms in optical lattices allow for accurate studies of many body dynamics. Rapid time-dependent modifications of optical lattice potentials may result in significant excitations in atomic systems. The dynamics in such a case is frequently quite incompletely described by standard applications of tight-binding models (such as e.g. Bose-Hubbard model or its extensions) that typically neglect the effect of the dynamics on the transformation between the real space and the tight-binding basis. We illustrate the importance of a proper quantum mechanical description using a multi-band extended Bose-Hubbard model with time-dependent Wannier functions. We apply it to situations, directly related to experiments.

preprint2013arXiv

Many-body Anderson localization in one dimensional systems

We show, using quasi-exact numerical simulations, that Anderson localization of one-dimensional particles in a disordered potential survives in the presence of attractive interaction between particles. The localization length of the composite particle can be computed analytically for weak disorder and is in good agreement with the quasi-exact numerical observations using Time Evolving Block Decimation. Our approach allows for simulation of the entire experiment including the final measurement of all atom positions.

preprint2013arXiv

Tunneling-Induced Restoration of the Degeneracy and the Time-Reversal Symmetry Breaking in Optical Lattices

We study the ground-state properties of bosons loaded into the $p$-band of a one dimensional optical lattice. We show that the phase diagram of the system is substantially affected by the anharmonicity of the lattice potential. In particular, for a certain range of tunneling strength, the full many-body ground state of the system becomes degenerate. In this region, an additional symmetry of the system, namely the parity of the occupation number of the chosen orbital, is spontaneously broken. The state with nonvanishing staggered angular momentum, which breaks the time-reversal symmetry, becomes the true ground state of the system.

preprint2013arXiv

Variational Bose-Hubbard model revisited

For strongly interacting bosons in optical lattices the standard description using Bose-Hubbard model becomes questionable. The role of excited bands becomes important. In such a situation we compare results of simulations using multiband Bose-Hubbard model with a recent proposition based on a time dependent variational approach. It is shown that the latter, in its original formulation, uses too small variational space leading often to spurious effects. Possible expansion of variational approach is discussed.

preprint2012arXiv

Dynamics of cold bosons in optical lattices: Effects of higher Bloch bands

The extended effective multiorbital Bose-Hubbard-type Hamiltonian which takes into account higher Bloch bands, is discussed for boson systems in optical lattices, with emphasis on dynamical properties, in relation with current experiments. It is shown that the renormalization of Hamiltonian parameters depends on the dimension of the problem studied. Therefore, mean field phase diagrams do not scale with the coordination number of the lattice. The effect of Hamiltonian parameters renormalization on the dynamics in reduced one-dimensional optical lattice potential is analyzed. We study both the quasi-adiabatic quench through the superfluid-Mott insulator transition and the absorption spectroscopy, that is energy absorption rate when the lattice depth is periodically modulated.

preprint2012arXiv

Frustration and time reversal symmetry breaking for Fermi and Bose-Fermi systems

The modulation of an optical lattice potential that breaks time-reversal symmetry enables the realization of complex tunneling amplitudes in the corresponding tight-binding model. For a superfluid Fermi gas in a triangular lattice potential with complex tunnelings the pairing function acquires a complex phase, so the frustrated magnetism of fermions can be realized.Bose-Fermi mixture of bosonic molecules and unbound fermions in the lattice shows also an interesting behavior. Due to boson-fermion coupling, the fermions become slaved by the bosons and the corresponding pairing function takes the complex phase determined by bosons. In the presence of bosons the Fermi system can reveal both gap and gapless superfluidity.

preprint2012arXiv

Numerical Computation of Dynamically Important Excited States of Many-Body Systems

We present an extension of the time-dependent Density Matrix Renormalization Group (t-DMRG), also known as Time Evolving Block Decimation algorithm (TEBD), allowing for the computation of dynamically important excited states of one-dimensional many-body systems. We show its practical use for analyzing the dynamical properties and excitations of the Bose-Hubbard model describing ultracold atoms loaded in an optical lattice from a Bose-Einstein condensate. This allows for a deeper understanding of nonadiabaticity in experimental realizations of insulating phases.

preprint2012arXiv

Quantum spin models with long-range interactions and tunnelings: A quantum Monte Carlo study

We use a quantum Monte Carlo method to investigate various classes of 2D spin models with long-range interactions at low temperatures. In particular, we study a dipolar XXZ model with U(1) symmetry that appears as a hard-core boson limit of an extended Hubbard model describing polarized dipolar atoms or molecules in an optical lattice. Tunneling, in such a model, is short-range, whereas density-density couplings decay with distance following a cubic power law. We investigate also an XXZ model with long-range couplings of all three spin components - such a model describes a system of ultracold ions in a lattice of microtraps. We describe an approximate phase diagram for such systems at zero and at finite temperature, and compare their properties. In particular, we compare the extent of crystalline, super?uid, and supersolid phases. Our predictions apply directly to current experiments with mesoscopic numbers of polar molecules and trapped ions.

preprint2011arXiv

Dipolar bosons on an optical lattice ring

We consider an ultra-small system of polarized bosons on an optical lattice with a ring topology interacting via long range dipole-dipole interactions. Dipoles polarized perpendicular to the plane of the ring reveal sharp transitions between different density wave phases. As the strength of the dipolar interactions is varied the behavior of the transitions is first-order like. For dipoles polarized in the plane of the ring the transitions between possible phases show pronounced sensitivity to the lattice depth. The abundance of possible configurations may be useful for quantum information applications.

preprint2011arXiv

Disordered spinor Bose-Hubbard model

We study the zero temperature phase diagram of the disordered spin-1 Bose-Hubbard model in a 2-dimensional square lattice. To this aim, we use a mean field Gutzwiller ansatz and a probabilistic mean field perturbation theory. The spin interaction induces two different regimes corresponding to a ferromagnetic and antiferromagnetic order. In the ferromagnetic case, the introduction of disorder reproduces analogous features of the disordered scalar Bose-Hubbard model, consisting in the formation of a Bose glass phase between Mott insulator lobes. In the antiferromagnetic regime the phase diagram differs more from the scalar case. Disorder in the chemical potential can lead to the disappearance of Mott insulator lobes with odd integer filling factor and, for sufficiently strong spin coupling, to Bose glass of singlets between even filling Mott insulator lobes. Disorder in the spinor coupling parameter results in the appearance of a Bose glass phase only between the n and n+1 lobes for n odd. Disorder in the scalar Hubbard interaction inhibits Mott insulator regions for occupation larger than a critical value.

preprint2011arXiv

Spin effects in Bose-Glass phases

We study the mechanism of formation of Bose glass (BG) phases in the spin-1 Bose Hubbard model when diagonal disorder is introduced. To this aim, we analyze first the phase diagram in the zero-hopping limit, there disorder induces superposition between Mott insulator (MI) phases with different filling numbers. Then BG appears as a compressible but still insulating phase. The phase diagram for finite hopping is also calculated with the Gutzwiller approximation. The bosons' spin degree of freedom introduces another scattering channel in the two-body interaction modifying the stability of MI regions with respect to the action of disorder. This leads to some peculiar phenomena such as the creation of BG of singlets, for very strong spin correlation, or the disappearance of BG phase in some particular cases where fluctuations are not able to mix different MI regions.

preprint2011arXiv

Tunable dipolar resonances and Einstein-de Haas effect in a Rb-87 atoms condensate

We study a spinor condensate of Rb-87 atoms in F = 1 hyperfine state confined in an optical dipole trap. Putting initially all atoms in mF = 1 component we observe a significant transfer of atoms to other, initially empty Zeeman states exclusively due to dipolar forces. Because of conservation of a total angular momentum the atoms going to other Zeeman components acquire an orbital angular momentum and circulate around the center of the trap. This is a realization of Einstein-de Haas effect in a system of cold gases. We show that the transfer of atoms via dipolar interaction is possible only when the energies of the initial and the final sates are equal. This condition can be fulfilled utilizing a resonant external magnetic field, which tunes energies of involved states via the linear Zeeman effect. We found that there are many final states of different spatial density which can be tuned selectively to the initial state. We show a simple model explaining high selectivity and controllability of weak dipolar interactions in the condensate of Rb-87 atoms.

preprint2009arXiv

Anderson localization of solitons

At low temperature, a quasi-one-dimensional ensemble of atoms with attractive interaction forms a bright soliton. When exposed to a weak and smooth external potential, the shape of the soliton is hardly modified, but its center-of-mass motion is affected. We show that in a spatially correlated disordered potential, the quantum motion of a bright soliton displays Anderson localization. The localization length can be much larger than the soliton size and could be observed experimentally.

preprint2009arXiv

Breakdown of adiabaticity when loading ultra-cold atoms in optical lattices

Realistic simulations of current ultra-cold atoms experiments in optical lattices show that the ramping up of the optical lattice is significantly nonadiabatic, implying that experimentally prepared Mott insulators are not really in the ground state of the atomic system. The nonadiabaticity is even larger in the presence of a secondary quasi-periodic lattice simulating "disorder". Alternative ramping schemes are suggested that improve the adiabaticity when the disorder is not too large.

preprint2004arXiv

Non-sequential double ionization of molecules

Double ionization of diatomic molecules by short linearly polarized laser pulses is analyzed. We consider the final stage of the ionization process, that is the decay of a highly excited two electron molecule, which is formed after re-scattering. The saddles of the effective adiabatic potential energy close to which simultaneous escape of electrons takes place are identified. Numerical simulations of the ionization of molecules show that the process can be dominated by either sequential or non-sequential events. In order to increase the ratio of non-sequential to sequential ionizations very short laser pulses should be applied.

preprint2002arXiv

Non-dispersive wave packets in periodically driven quantum systems

With the exception of the harmonic oscillator, quantum wave-packets usually spread as time evolves. We show here that, using the nonlinear resonance between an internal frequency of a system and an external periodic driving, it is possible to overcome this spreading and build non-dispersive (or non-spreading) wave-packets which are well localized and follow a classical periodic orbit without spreading. From the quantum mechanical point of view, the non-dispersive wave-packets are time periodic eigenstates of the Floquet Hamiltonian, localized in the nonlinear resonance island. We discuss the general mechanism which produces the non-dispersive wave-packets, with emphasis on simple realization in the electronic motion of a Rydberg electron driven by a microwave field. We show the robustness of such wavepackets for a model one-dimensional as well as for realistic three dimensional atoms. We consider their essential properties such as the stability versus ionization, the characteristic energy spectrum and long lifetimes. The requirements for experiments aimed at observing such non-dispersive wave-packets are also considered. We also discuss other related phenomena in atomic and molecular physics as well as possible further extensions of the theory.

preprint1996arXiv

Nonuniversality in level dynamics

Statistical properties of parametric motion in ensembles of Hermitian banded random matrices are studied. We analyze the distribution of level velocities and level curvatures as well as their correlation functions in the crossover regime between three universality classes. It is shown that the statistical properties of level dynamics are in general non-universal and strongly depend on the way in which the parametric dynamics is introduced.