Researcher profile

Jakub Zakrzewski

Jakub Zakrzewski contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
14topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

22 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.