Search

The lecture notes on "Many-body Quantum Dynamics with MCTDH-X," adapted from the 2023 Heidelberg MCTDH Summer School, provide an in-depth exploration of the Multiconfigurational Time-Dependent Hartree approach for indistinguishable particles. They serve as a comprehensive guide for understanding and utilizing the MCTDH-X software for both bosonic and fermionic systems. The tutorial begins with an introduction to the MCTDH-X software, highlighting its capability to handle various quantum systems, including those with internal degrees of freedom and long-range interactions. The theoretical foundation is then laid out on how to solve the time-dependent and time-independent Schr\"odinger equations for many-body systems. The workflow section provides practical instructions on setting up and executing simulations using MCTDH-X. Detailed benchmarks against exact solutions are presented, showcasing the accuracy and reliability of the software in ground-state and dynamic simulations. The notes then delve into the dynamics of quantum systems, covering relaxation processes, time evolution, and the analysis of propagation for both bosonic and fermionic particles. The discussion includes the interpretation of various physical quantities such as energy, density distributions, and orbital occupations. Advanced features of MCTDH-X are also explored in the last section, including the calculation of correlation functions and the creation of visualizations through video tutorials. The notes conclude with a Linux/UNIX command cheat sheet, facilitating ease of use for users operating the software on different systems. Overall, these lecture notes provide a valuable resource for researchers and students in the field of quantum dynamics, offering both theoretical insights and practical guidance on the use of MCTDH-X for studying complex many-body systems., Comment: Adapted from a series of lectures given during the 2023 MCTDH Summer School in Heidelberg (Germany); 48 pages, 30 figures

Despite acute interest in the dynamics of non-Hermitian systems, there is a lack of consensus in the mathematical formulation of non-Hermitian quantum mechanics in the community. Different methodologies are used in the literature to study non-Hermitian dynamics. This ranges from consistent frameworks like biorthogonal quantum mechanics and metric approach characterized by modified inner products, to normalization by time-dependent norms inspired by open quantum systems. In this work, we systematically explore the similarities and differences among these various methods. Utilizing illustrative models with exact solutions, we demonstrate that these methods produce not only quantitatively different results but also distinct physical interpretations. For dissipative systems where non-Hermiticity arises as an approximation, we find that the normalization method in the $\mathcal{PT}$-broken regime closely aligns with the full master equation solutions. In contrast, for quantum systems where non-Hermiticity can be engineered exactly, incorporating metric dynamics is crucial for the probabilistic interpretation of quantum mechanics, necessitating the generalization of unitary symmetry to non-Hermitian systems. This study lays the groundwork for further exploration of non-Hermitian Hamiltonians, potentially leveraging generalized symmetries for novel physical phenomena., Comment: 12 pages, 2 figures

Quasiperiodic potentials can be used to interpolate between localization and delocalization in one dimension. With the rise of optical platforms engineering dipolar interactions, a key question is the stability of quasicrystalline phases under these long-range interactions. In this work, we study repulsive ultracold dipolar fermions in a quasiperiodic optical lattice to characterize the behavior of interacting quasicrystals. We simulate the full time evolution of the typical experimental protocols used to probe quasicrystalline order and localization properties. We extract experimentally measurable dynamical observables and correlation functions to characterize the three phases observed in the noninteracting setting: localized, intermediate, and extended. We then study the stability of such phases to repulsive dipolar interactions. We find that dipolar interactions can completely alter the shape of the phase diagram by stabilizing the intermediate phase, mostly at the expense of the extended phase. Moreover, in the strongly interacting regime, a resonance-like behavior characterized by density oscillations appears. Remarkably, strong dipolar repulsions can also localize particles even in the absence of quasiperiodicity if the primary lattice is sufficiently deep. Our work shows that dipolar interactions in a quasiperiodic potential can give rise to a complex, tuneable coexistence of localized and extended quantum states., Comment: 7 pages, 4 figures; supplementary material available with 11 pages, 10 figures. Updated version with new discussion and explanation of results

We provide a quantitative blueprint for realizing two-dimensional quantum simulators employing ultracold dipolar molecules or magnetic atoms by studying their accuracy in predicting ground state properties of lattice models with long-range interactions. For experimentally relevant ranges of potential depths, interaction strengths, particle fillings, and geometric configurations, we map out the agreement between the state prepared in the quantum simulator and the target lattice state. We do so by separately calculating numerically exact many-body wave functions in the continuum and single- or multi-band lattice representations, and building their many-body state overlaps. While the agreement between quantum simulator and single-band models is good for deep optical lattices with weaker interactions and low particle densities, the higher band population rapidly increases for shallow lattices, stronger interactions, and in particular above half filling. This induces drastic changes to the properties of the simulated ground state, potentially leading to false predictions. Furthermore, we show that the interplay between commensurability and interactions can lead to quasidegeneracies, rendering a faithful ground state preparation even more challenging., Comment: 8 pages, 4 figures; supplementary material available with 8 pages, 4 figures

A super-Tonks-Giradeau gas is a highly excited yet stable quantum state of strongly attractive bosons confined to one dimension. This state can be obtained by quenching the interparticle interactions from the ground state of a strongly repulsive Tonks-Girardeau gas to the strongly attractive regime. While the super-Tonks-Girardeau quench with contact interactions has been thoroughly studied, less is known about the stability of such a procedure when long-range interactions come into play. This is a particularly important question in light of recent advances in controlling ultracold atoms with dipole-dipole interactions. In this study, we thus simulate a super-Tonks-Girardeau quench on dipolar bosons in a one-dimensional optical lattice and investigate their dynamics for many different initial states and fillings. By calculating particle density, correlations, entropy measures, and natural occupations, we establish the regimes of stability as a function of dipolar interaction strength. For an initial unit-filled Mott state, stability is retained at weak dipolar interactions. For cluster states and doubly-filled Mott states, instead, dipolar interactions eventually lead to complete evaporation of the initial state and thermalization consistent with predictions from random matrix theory. Remarkably, though, dipolar interactions can be tuned to achieve longer-lived prethermal states before the eventual thermalization. Our study highlights the potential of long-range interactions to explore new mechanisms to steer and stabilize excited quantum states of matter., Comment: 16 pages, 17 figures, 4 appendices

We devise a generic and experimentally accessible recipe to prepare boundary states of topological or nontopological quantum systems through an interplay between coherent Hamiltonian dynamics and local dissipation. Intuitively, our recipe harnesses the spatial structure of boundary states which vanish on sublattices where losses are suitably engineered. This yields unique nontrivial steady states that populate the targeted boundary states with infinite lifetimes while all other states are exponentially damped in time. Remarkably, applying loss only at one boundary can yield a unique steady state localized at the very same boundary. We detail our construction and rigorously derive full Liouvillian spectra and dissipative gaps in the presence of a spectral mirror symmetry for a one-dimensional Su-Schrieffer-Heeger model and a two-dimensional Chern insulator. We outline how our recipe extends to generic noninteracting systems.

We explore anomalous skin effects at non-Hermitian impurities by studying their interplay with potential disorder and by exactly solving a minimal lattice model. A striking feature of the solvable single-impurity model is that the presence of anisotropic hopping terms can induce a scale-free accumulation of all eigenstates opposite to the bulk hopping direction, although the nonmonotonic behavior is fine tuned and further increasing such hopping weakens and eventually reverses the effect. The interplay with bulk potential disorder, however, qualitatively enriches this phenomenology leading to a robust nonmonotonic localization behavior as directional hopping strengths are tuned. Nonmonotonicity persists even in the limit of an entirely Hermitian bulk with a single non-Hermitian impurity., Comment: Accepted version

Non-Hermiticity in quantum Hamiltonians leads to nonunitary time evolution and possibly complex energy eigenvalues, which can lead to a rich phenomenology with no Hermitian counterpart. In this work, we study the dynamics of an exactly solvable non-Hermitian system, hosting both $\mathcal{PT}$-symmetric and $\mathcal{PT}$-broken modes subject to a linear quench. Employing a fully consistent framework, in which the Hilbert space is endowed with a nontrivial dynamical metric, we analyze the dynamics of the generated defects. In contrast to Hermitian systems, our study reveals that PT -broken time evolution leads to defect freezing and hence the violation of adiabaticity. This physics necessitates the so-called metric framework, as it is missed by the oft used approach of normalizing quantities by the time-dependent norm of the state. Our results are relevant for a wide class of experimental systems., Comment: Main text: 7 pages and 3 figures

With rapid progress in control and manipulation of ultracold magnetic atoms and dipolar molecules, the quantum simulation of lattice models with strongly interacting dipole-dipole interactions (DDI) and high densities is now within experimental reach. This rapid development raises the issue about the validity of quantum simulation in such regimes. In this study, we address this question by performing a full quantitative comparison between the continuum description of a one-dimensional gas of dipolar bosons in an optical lattice, and the single-band Bose-Hubbard lattice model that it quantum simulates. By comparing energies and density distributions, and by calculating direct overlaps between the continuum and lattice many-body wavefunctions, we demonstrate that in regimes of strong DDI and high densities the continuum system fails to recreate the desired lattice model. Two-band Hubbard models become necessary to reduce the discrepancy observed between continuum and lattice descriptions, but appreciable deviations in the density profile still remain. Our study elucidates the role of strong DDI in generating physics beyond lowest-band descriptions and should offer a guideline for the calibration of near-term dipolar quantum simulators., Comment: Fixed typos, added more references

The ground states of noninteracting fermions in one-dimension with chiral symmetry form a class of topological band insulators, described by a topological invariant that can be related to the Zak phase. Recently, a generalization of this quantity to mixed states - known as the ensemble geometric phase (EGP) - has emerged as a robust way to describe topology at non-zero temperature. By using this quantity, we explore the nature of topology allowed for dissipation beyond a Lindblad description, to allow for coupling to external baths at finite temperatures. We introduce two main aspects to the theory of mixed state topology. First, we discover topological phase transitions as a function of the temperature T, manifesting as changes in winding number of the EGP accumulated over a closed loop in parameter space. We characterize the nature of these transitions and reveal that the corresponding non-equilibrium steady state at the transition can exhibit a nontrivial structure - contrary to previous studies where it was found to be in a fully mixed state. Second, we demonstrate that the EGP itself becomes quantized when key symmetries are present, allowing it to be viewed as a topological marker which can undergo equilibrium topological transitions at non-zero temperatures., Comment: 13 pages, 6 figures; comments welcome

We employ quench dynamics as an effective tool to probe different universality classes of topological phase transitions. Specifically, we study a model encompassing both Dirac-like and nodal loop criticalities. Examining the Kibble-Zurek scaling of topological defect density, we discover that the scaling exponent is reduced in the presence of extended nodal loop gap closures. For a quench through a multicritical point, we also unveil a path-dependent crossover between two sets of critical exponents. Bloch state tomography finally reveals additional differences in the defect trajectories for sudden quenches. While the Dirac transition permits a static trajectory under specific initial conditions, we find that the underlying nodal loop leads to complex time-dependent trajectories in general. In the presence of a nodal loop, we find, generically, a mismatch between the momentum modes where topological defects are generated and where dynamical quantum phase transitions occur. We also find notable exceptions where this correspondence breaks down completely., Comment: 8 pages, 7 figures; references added

In 2D semiconductors and insulators, the Chern number of the valence band Bloch state is an important quantity that has been linked to various material properties, such as the topological order. We elaborate that the opacity of 2D materials to circularly polarized light over a wide range of frequencies, measured in units of the fine structure constant, can be used to extract a spectral function that frequency-integrates to the Chern number, offering a simple optical experiment to measure it. This method is subsequently generalized to finite temperature and locally on every lattice site by a linear response theory, which helps to extract the Chern marker that maps the Chern number to lattice sites. The long range response in our theory corresponds to a Chern correlator that acts like the internal fluctuation of the Chern marker, and is found to be enhanced in the topologically nontrivial phase. Finally, from the Fourier transform of the valence band Berry curvature, a nonlocal Chern marker is further introduced, whose decay length diverges at topological phase transitions and therefore serves as a faithful indicator of the transitions, and moreover can be interpreted as a Wannier state correlation function. The concepts discussed in this work explore multi-faceted aspects of topology and should help address the impact of system inhomogeneities., Comment: Reworked version

Pauli crystals are ordered geometric structures that emerge in trapped noninteracting fermionic systems due to their underlying Pauli repulsion. The deformation of Pauli crystals - often called melting - has been recently observed in experiments, but the mechanism that leads to it remains unclear. We address this question by studying the melting dynamics of N=6 fermions as a function of periodic driving and experimental imperfections in the trap (anisotropy and anharmonicity) by employing a combination of numerical simulations and Floquet theory. Surprisingly, we reveal that the melting of Pauli crystals is not simply a direct consequence of an increase in system energy, but is instead related to the trap geometry and the population of the Floquet modes. We show that the melting is absent in traps without imperfections and triggered only by a sufficiently large shaking amplitude in traps with imperfections., Comment: 14 pages, 9 figures

We study a one-dimensional array of bosons with infinite-range interactions mediated by a laser-driven dissipative optical cavity. The cavity-mediated infinite-range interactions open up a new pathway to fermionization, hitherto only known for dipolar bosons due to their long-range interactions. In parameter ranges attainable in state-of-the-art experiments, we systematically compare observables for bosons and fermions with infinite-range interactions. At large enough laser pump power, many observables, including density distributions in real and momentum space, correlation functions, eigenvalues of the one-body density matrix, and superradiance order parameter, become identical for bosons and fermions. We map out the emergence of this cavity-induced fermionization as a function of pump power and contact interactions. We discover that cavity-mediated interactions can compensate a reduction by several orders of magnitude in the strength of the contact interactions needed to trigger fermionization., Comment: main text: 7 pages, 3 figures - supplement: 7 pages, 4 figures

Topological order in solid state systems is often calculated from the integration of an appropriate curvature function over the entire Brillouin zone. At topological phase transitions where the single particle spectral gap closes, the curvature function diverges and changes sign at certain high symmetry points in the Brillouin zone. These generic properties suggest the introduction of a supervised machine learning scheme that uses only the curvature function at the high symmetry points as input data. We apply this scheme to a variety of interacting topological insulators in different dimensions and symmetry classes, and demonstrate that an artificial neural network trained with the noninteracting data can accurately predict all topological phases in the interacting cases with very little numerical effort. Intriguingly, the method uncovers a ubiquitous interaction-induced topological quantum multicriticality in the examples studied., Comment: 8 pages, 3 figures - corrected typos in Fig. 1

The competition between short-range and cavity-mediated infinite-range interactions in a cavity-boson system leads to the existence of a superfluid phase and a Mott-insulator phase within the self-organized regime. In this work, we quantitatively compare the steady-state phase boundaries of this transition measured in experiments and simulated using the Multiconfigurational Time-Dependent Hartree Method for Indistinguishable Particles. To make the problem computationally feasible, we represent the full system by the exact many-body wave function of a two-dimensional four-well potential. We argue that the validity of this representation comes from the nature of both the cavity-atomic system and the Bose-Hubbard physics. Additionally we show that the chosen representation only induces small systematic errors, and that the experimentally measured and theoretically predicted phase boundaries agree reasonably. We thus demonstrate a new approach for the quantitative numerical determination of the superfluid--Mott-insulator phase boundary., Comment: 30 pages, 9 figures, 1 table. Submission to SciPost

The Kitaev model on the honeycomb lattice is a paradigmatic system known to host a wealth of nontrivial topological phases and Majorana edge modes. In the static case, the Majorana edge modes are nondispersive. When the system is periodically driven in time, such edge modes can disperse and become chiral. We obtain the full phase diagram of the driven model as a function of the coupling and the driving period. We characterize the quantum criticality of the different topological phase transitions in both the static and driven model via the notions of Majorana-Wannier state correlation functions and momentum-dependent fidelity susceptibilities. We show that the system hosts cross-dimensional universality classes: although the static Kitaev model is defined on a 2D honeycomb lattice, its criticality falls into the universality class of 1D linear Dirac models. For the periodically driven Kitaev model, besides the universality class of prototype 2D linear Dirac models, an additional 1D nodal loop type of criticality exists owing to emergent time-reversal and mirror symmetries, indicating the possibility of engineering multiple universality classes by periodic driving. The manipulation of time-reversal symmetry allows the periodic driving to control the chirality of the Majorana edge states., Comment: 14 pages, 11 figures

Single-shot images are the standard readout of experiments with ultracold atoms -- the tarnished looking glass into their many-body physics. The efficient extraction of observables from single-shot images is thus crucial. Here, we demonstrate how artificial neural networks can optimize this extraction. In contrast to standard averaging approaches, machine learning allows both one- and two-particle densities to be accurately obtained from a drastically reduced number of single-shot images. Quantum fluctuations and correlations are directly harnessed to obtain physical observables for bosons in a tilted double-well potential at an unprecedented accuracy. Strikingly, machine learning also enables a reliable extraction of momentum-space observables from real-space single-shot images and vice versa. This obviates the need for a reconfiguration of the experimental setup between in-situ and time-of-flight imaging, thus potentially granting an outstanding reduction in resources., Comment: 7+8 pages, 3+8 figures, software available at http://ultracold.org

We explore the impact of commensurate multifrequency driving protocols on the stability of topological edge modes in topological 1D systems of spinless fermions. Using Floquet theory, we show that all the topological phase transitions can be mapped to their single-frequency counterparts in terms of effective driving parameters. This drastic simplification can be explained by considering the gap closures in the quasienergy dispersion. These gap closures are pinned to high-symmetry points, where all the effective Floquet Hamiltonians collapse to the same form. While for all protocols all topological phase transitions coincide, the gap size and the number of edge modes in the quasienergy spectra vary considerably depending on the chosen driving pattern. This gives rise to a full range of edge states with different degrees of localization, from highly localized to completely delocalized. Switching between different driving protocols then suggests a dynamical control of the localization and stability of topological edge modes, with possible applications in quantum computation. We illustrate our findings on three paradigmatic fermionic systems -- namely the Kitaev chain, the Su-Schrieffer-Heeger model, and the Creutz ladder, and demonstrate how to control the localization length of the corresponding edge states., Comment: 16 pages, 15 figures

Topological phase transitions track changes in topological properties of a system and occur in real materials as well as quantum engineered systems, all of which differ greatly in terms of dimensionality, symmetries, interactions, and driving, and hence require a variety of techniques and concepts to describe their topological properties. For instance, depending on the system, topology may be accessed from single-particle Bloch wave functions, Green's functions, or many-body wave functions. We demonstrate that despite this diversity, all topological phase transitions display a universal feature: namely, a divergence of the curvature function that composes the topological invariant at the critical point. This feature can be exploited via a renormalization-group-like methodology to describe topological phase transitions. This approach serves to extend notions of correlation function, critical exponents, scaling laws and universality classes used in Landau theory to characterize topological phase transitions in a unified manner., Comment: 7 pages, 3 figures, to appear in Europhysics Letters

We introduce and describe the multiconfigurational time-depenent Hartree for indistinguishable particles (MCTDH-X) software. This powerful tool allows the investigation of ground state properties and dynamics of interacting quantum many-body systems in different spatial dimensions. The MCTDH-X software is a set of programs and scripts to compute, analyze, and visualize solutions for the time-dependent and time-independent many-body Schr\"{o}dinger equation for indistinguishable quantum particles. As the MCTDH-X software represents a general solver for the Schr\"{o}dinger equation, it is applicable to a wide range of problems in the fields of atomic, optical, molecular physics as well as condensed matter systems. In particular, it can be used to study light-matter interactions, correlated dynamics of electrons, as well as some aspects related to quantum information and computing. The MCTDH-X software solves a set of non-linear coupled working equations based on the application of the variational principle to the Schr\"{o}dinger equation. These equations are obtained by using an ansatz for the many-body wavefunction that is a time-dependent expansion in a set of time-dependent many-body basis states. The time-dependence of the basis set enables MCTDH-X to deal with quantum dynamics at a superior accuracy as compared to, for instance, exact diagonalization approaches. Herein, we give an introduction to the MCTDH-X software via an easy-to-follow tutorial with a focus on accessibility. We use the double well to illustrate the fermionization of bosonic particles, the crystallization of fermionic particles, characteristics of the superfluid and Mott-insulator quantum phases in Hubbard models, and even dynamical quantum phase transitions. Our tutorial guides the potential user to apply the MCTDH-X software also to more complex systems., Comment: 38 pages, 15 figures (22 pages, 9 figures in the supplementary material). Software is available on http://ultracold.org, input files for the tutorial on http://ultracold.org/data/tutorial_input_files.zip, tutorial videos on https://www.youtube.com/playlist?list=PLJIFUqmSeGBKxmLcCuk6dpILnni_uIFGu, and a supplementary video on https://www.youtube.com/watch?v=l2UsTPmJ6po

We explore the role of atomic correlations in a harmonically trapped Bose-Einstein condensate coupled to a dissipative cavity, where both the atoms and the cavity are blue detuned from the external pumping laser. Using a genuine many-body approach that goes beyond mean-field, we extract density distributions and many-body correlations to unveil a pathway to chaos at large pump power through a hierarchical self-organization of the atoms, where the atoms transition from a single-well optical lattice to a double-well optical lattice. Correlated states of the atoms emerge and are characterized by local superfluid correlations in phases which are globally superfluid or Mott insulating. Local superfluid-Mott transitions are precluded by a dynamical instability to chaos which occurs via quasiperiodic attractors. Our results explain the mechanism behind the dynamical instabilities observed in experiments., Comment: 14 pages, 9 figures (including 7 pages, 5 figures in the supplementary material)

We demonstrate that the prototypical two-dimensional Chern insulator hosts exotic quantum multi-criticality in the presence of an appropriate periodic driving: a linear Dirac-like transition coexists with a nodal loop-like transition caused by emerging symmetries. The existence of multiple universality classes and scaling laws can be unambiguously captured by a single renormalization group approach based on the stroboscopic Floquet Hamiltonian, regardless of whether the topological transition is associated with the anomalous edge modes or not., Comment: 10 pages, 7 figures, includes supplementary material

We investigate harmonically-trapped, laser-pumped bosons with infinite-range interactions induced by a dissipative high-finesse red-detuned optical cavity with numerical and analytical methods. We obtain multiple cavity and atomic observables as well as the full phase diagram of the system using the multiconfigurational time-dependent Hartree method for indistinguishable particles (MCTDH-X) approach. Besides the transition from an unorganized normal phase to a superradiant phase where atoms self-organize, we focus on an in-depth investigation of the self-organized superfluid to self-organized Mott insulator phase transition in the superradiant phase as a function of the cavity-atom coupling. The numerical results are substantiated by an analytical study of an effective Bose-Hubbard model. We numerically analyze cavity fluctuations and emergent strong correlations between atoms in the many-body state across the Mott transition via the atomic density distributions and Glauber correlation functions. Unexpectedly, the weak harmonic trap leads to features like a lattice switching between the two symmetry-broken $\mathbb{Z}_2$ configurations of the untrapped system and a reentrance of superfluidity in the Mott insulating phase. Our analytical considerations quantitatively explain the numerically observed correlation features., Comment: 17 pages, 9 figures

The topological phase transitions in static and periodically driven Kitaev chains are investigated by means of a renormalization group (RG) approach. These transitions, across which the numbers of static or Floquet Majorana edge modes change, are accompanied by divergences of the relevant Berry connections. These divergences at certain high symmetry points in momentum space form the basis of the RG approach, through which topological phase boundaries are identified as a function of system parameters. We also introduce several aspects to characterize the quantum criticality of the topological phase transitions in both static and Floquet systems: a correlation function that measures the overlap of Majorana-Wannier functions, the decay length of the Majorana edge mode and a scaling law relating the critical exponents. These indicate a common universal critical behavior for topological phase transitions, in both static and periodically driven chains. For the latter, the RG flows additionally display intriguing features related to gap closures at non-high symmetry points due to momentarily frozen dynamics., Comment: 14 pages, 12 figures

We consider laser-pumped one-dimensional two-component bosons in a parabolic trap embedded in a high-finesse optical cavity. Above a threshold pump power, the photons that populate the cavity modify the effective atom trap and mediate a coupling between the two components of the Bose-Einstein condensate. We calculate the ground state of the laser-pumped system and find different stages of self-organization depending on the power of the laser. The modified potential and the laser-mediated coupling between the atomic components give rise to rich many-body physics: an increase of the pump power triggers a self-organization of the atoms while an even larger pump power causes correlations between the self-organized atoms -- the BEC becomes fragmented and the reduced density matrix acquires multiple macroscopic eigenvalues. In this fragmented superradiant state, the atoms can no longer be described as two-level systems and the mapping of the system to the Dicke model breaks down., Comment: 8 pages, 3 figures, software available at http://ultracold.org

We numerically obtain the full time-evolution of a parametrically-driven dissipative Bose-Einstein condensate in an optical cavity and investigate the implications of driving for the phase diagram. Beyond the normal and superradiant phases, a third nonequilibrium phase emerges as a manybody parametric resonance. This dynamical normal phase switches between two symmetry-broken superradiant configurations. The switching implies a breakdown of the system's mapping to the Dicke model. Unlike the other phases, the dynamical normal phase shows features of nonintegrability and thermalization., Comment: 5 pages, 3 figures

Time periodic modulations of the transverse field in the closed XY spin-1/2 chain generate a very rich dynamical phase diagram, with a hierarchy of topological phases characterized by differing numbers of Floquet-Majorana modes. We show that this rich phase diagram survives when the system is coupled to dissipative end reservoirs. Circumventing the obstacle of preparing and measuring quasi energy configurations endemic to Floquet-Majorana detection schemes, we show that stroboscopic heat transport and spin density are robust observables to detect both the dynamical phase transitions and Majorana modes. In particular, we find that the derivative of the heat current, with respect to a control parameter, changes sign at the boundaries separating topological phases with different numbers of Floquet Majorana modes. We present a simple scheme to directly count the number of Floquet-Majorana modes in a phase from the Fourier transform of the local spin density profile. Our results are valid provided the anisotropies are not strong and can be easily implemented in quantum engineered systems., Comment: 5 pages, 3 figures