High order non-unitary split-step decomposition of unitary operators

We propose a high order numerical decomposition of exponentials of hermitean operators in terms of a product of exponentials of simple terms, following an idea which has been pioneered by M. Suzuki, however implementing it for complex coefficients. We outline a convenient fourth order formula which can be written compactly for arbitrary number of noncommuting terms in the Hamiltonian and which is superiour to the optimal formula with real coefficients, both in complexity and accuracy. We show asymptotic stability of our method for sufficiently small time step and demonstrate its efficiency and accuracy in different numerical models.Comment: 10 pages, 4 figures (5 eps files) Submitted to J. of Phys. A: Math. Ge

Berry-Robnik level statistics in a smooth billiard system

Berry-Robnik level spacing distribution is demonstrated clearly in a generic quantized plane billiard for the first time. However, this ultimate semi-classical distribution is found to be valid only for extremely small semi-classical parameter (effective Planck's constant) where the assumption of statistical independence of regular and irregular levels is achieved. For sufficiently larger semiclassical parameter we find (fractional power-law) level repulsion with phenomenological Brody distribution providing an adequate global fit.Comment: 10 pages in LaTeX with 4 eps figures include

Exact solution of Markovian master equations for quadratic fermi systems: thermal baths, open XY spin chains, and non-equilibrium phase transition

We generalize the method of third quantization to a unified exact treatment of Redfield and Lindblad master equations for open quadratic systems of n fermions in terms of diagonalization of 4n x 4n matrix. Non-equilibrium thermal driving in terms of the Redfield equation is analyzed in detail. We explain how to compute all physically relevant quantities, such as non-equilibrium expectation values of local observables, various entropies or information measures, or time evolution and properties of relaxation. We also discuss how to exactly treat explicitly time dependent problems. The general formalism is then applied to study a thermally driven open XY spin 1/2 chain. We find that recently proposed non-equilibrium quantum phase transition in the open XY chain survives the thermal driving within the Redfield model. In particular, the phase of long-range magnetic correlations can be characterized by hypersensitivity of the non-equilibrium-steady state to external (bath or bulk) parameters. Studying the heat transport we find negative thermal conductance for sufficiently strong thermal driving, as well as non-monotonic dependence of the heat current on the strength of the bath coupling.Comment: 24 pages, 12 figures, submitted to New Journal of Physics, Focus issue "Quantum Information and Many-Body Theory

The Effects of Symmetries on Quantum Fidelity Decay

We explore the effect of a system's symmetries on fidelity decay behavior. Chaos-like exponential fidelity decay behavior occurs in non-chaotic systems when the system possesses symmetries and the applied perturbation is not tied to a classical parameter. Similar systems without symmetries exhibit faster-than-exponential decay under the same type of perturbation. This counter-intuitive result, that extra symmetries cause the system to behave in a chaotic fashion, may have important ramifications for quantum error correction.Comment: 5 pages, 3 figures, to be published Phys. Rev. E Rapid Communicatio

A map from 1d Quantum Field Theory to Quantum Chaos on a 2d Torus

Dynamics of a class of quantum field models on 1d lattice in Heisenberg picture is mapped into a class of `quantum chaotic' one-body systems on configurational 2d torus (or 2d lattice) in Schr\" odinger picture. Continuum field limit of the former corresponds to quasi-classical limit of the latter.Comment: 4 pages in REVTeX, 1 eps-figure include

Engineering fidelity echoes in Bose-Hubbard Hamiltonians

We analyze the fidelity decay for a system of interacting bosons described by a Bose-Hubbard Hamiltonian. We find echoes associated with "non-universal" structures that dominate the energy landscape of the perturbation operator. Despite their classical origin, these echoes persist deep into the quantum (perturbative) regime and can be described by an improved random matrix modeling. In the opposite limit of strong perturbations (and high enough energies), classical considerations reveal the importance of self-trapping phenomena in the echo efficiency.Comment: 6 pages, use epl2.cls class, 5 figures Cross reference with nlin, quant-phy

Regular and Irregular States in Generic Systems

In this work we present the results of a numerical and semiclassical analysis of high lying states in a Hamiltonian system, whose classical mechanics is of a generic, mixed type, where the energy surface is split into regions of regular and chaotic motion. As predicted by the principle of uniform semiclassical condensation (PUSC), when the effective $\hbar$ tends to 0, each state can be classified as regular or irregular. We were able to semiclassically reproduce individual regular states by the EBK torus quantization, for which we devise a new approach, while for the irregular ones we found the semiclassical prediction of their autocorrelation function, in a good agreement with numerics. We also looked at the low lying states to better understand the onset of semiclassical behaviour.Comment: 25 pages, 14 figures (as .GIF files), high quality figures available upon reques

Separating the regular and irregular energy levels and their statistics in Hamiltonian system with mixed classical dynamics

We look at the high-lying eigenstates (from the 10,001st to the 13,000th) in the Robnik billiard (defined as a quadratic conformal map of the unit disk) with the shape parameter $\lambda=0.15$. All the 3,000 eigenstates have been numerically calculated and examined in the configuration space and in the phase space which - in comparison with the classical phase space - enabled a clear cut classification of energy levels into regular and irregular. This is the first successful separation of energy levels based on purely dynamical rather than special geometrical symmetry properties. We calculate the fractional measure of regular levels as $\rho_1=0.365\pm 0.01$ which is in remarkable agreement with the classical estimate $\rho_1=0.360\pm 0.001$. This finding confirms the Percival's (1973) classification scheme, the assumption in Berry-Robnik (1984) theory and the rigorous result by Lazutkin (1981,1991). The regular levels obey the Poissonian statistics quite well whereas the irregular sequence exhibits the fractional power law level repulsion and globally Brody-like statistics with $\beta = 0.286\pm0.001$. This is due to the strong localization of irregular eigenstates in the classically chaotic regions. Therefore in the entire spectrum we see that the Berry-Robnik regime is not yet fully established so that the level spacing distribution is correctly captured by the Berry-Robnik-Brody distribution (Prosen and Robnik 1994).Comment: 20 pages, file in plain LaTeX, 7 figures upon request submitted to J. Phys. A. Math. Gen. in December 199

Dephasing-induced diffusive transport in anisotropic Heisenberg model

We study transport properties of anisotropic Heisenberg model in a disordered magnetic field experiencing dephasing due to external degrees of freedom. In the absence of dephasing the model can display, depending on parameter values, the whole range of possible transport regimes: ideal ballistic conduction, diffusive, or ideal insulating behavior. We show that the presence of dephasing induces normal diffusive transport in a wide range of parameters. We also analyze the dependence of spin conductivity on the dephasing strength. In addition, by analyzing the decay of spin-spin correlation function we discover a presence of long-range order for finite chain sizes. All our results for a one-dimensional spin chain at infinite temperature can be equivalently rephrased for strongly-interacting disordered spinless fermions.Comment: 15 pages, 9 PS figure

Driven-dissipative many-body pairing states for cold fermionic atoms in an optical lattice

We discuss the preparation of many-body states of cold fermionic atoms in an optical lattice via controlled dissipative processes induced by coupling the system to a reservoir. Based on a mechanism combining Pauli blocking and phase locking between adjacent sites, we construct complete sets of jump operators describing coupling to a reservoir that leads to dissipative preparation of pairing states for fermions with various symmetries in the absence of direct inter-particle interactions. We discuss the uniqueness of these states, and demonstrate it with small-scale numerical simulations. In the late time dissipative dynamics, we identify a "dissipative gap" that persists in the thermodynamic limit. This gap implies exponential convergence of all many-body observables to their steady state values. We then investigate how these pairing states can be used as a starting point for the preparation of the ground state of Fermi-Hubbard Hamiltonian via an adiabatic state preparation process also involving the parent Hamiltonian of the pairing state. We also provide a proof-of-principle example for implementing these dissipative processes and the parent Hamiltonians of the pairing states, based on Yb171 atoms in optical lattice potentials.Comment: as published in the focus issue on Bose Condensation Phenomena in Atomic and Solid State Physics, New J. Phys. 14 (2012