Time-evolving Block Decimation Research Articles

Single-excitation correlation dynamics of polar molecules on 1D lattices: driving by external electric fields

Abstract We study a LiCs strongly-interacting 1D molecular lattice system in the Mott insulator regime, allowing one molecule per site prepared in the lowest electronic and vibrational state in presence of a perpendicular electric field. At large intermolecular distance and low filling, the dipole-dipole interaction in the nearest-neighbor approximation governs the dynamics exchange of the two allowed rotational levels, we analyze it for several field strengths. The generated von Neumann entanglement entropy and its monotonically growth in short-term evolutions is shown. The potential of these molecular systems to be used in quantum information protocols is enhanced by the early emergence of a long-range order in the single-particle density matrix correlator throughout the lattice as here we present. Contrary of the widely used two-level approach, our results clearly show that the full four set of internal rotational projections must be taken into account. The numerical simulations are performed by means of the Time-Evolving Block Decimation algorithm based on the Matrix Product State formalism and the Susuki-Trotter decomposition.

Topological transition between competing orders in quantum spin chains

We study quantum phase transitions between competing orders in one-dimensional spin systems. We focus on systems that can be mapped to a dual-field double sine-Gordon model as a bosonized effective field theory. This model contains two pinning potential terms of dual fields that stabilize competing orders and allows different types of quantum phase transition to happen between two ordered phases. At the transition point, elementary excitations change from the topological soliton of one of the dual fields to that of the other, thus it can be characterized as a topological transition. We compute the dynamical susceptibilities and the entanglement entropy, which gives us access to the central charge, of the system using a numerical technique of infinite time-evolving block decimation and characterize the universality class of the transition as well as the nature of the order in each phase. The possible realizations of such transitions in experimental systems both for condensed matter and cold atomic gases are also discussed.

Energy transport in a disordered spin chain with broken U(1) symmetry: Diffusion, subdiffusion, and many-body localization

We explore the physics of the disordered XYZ spin chain using two complementary numerical techniques: exact diagonalization (ED) on chains of up to 17 spins, and time-evolving block decimation (TEBD) on chains of up to 400 spins. Our principal findings are as follows. First, the clean XYZ spin chain shows ballistic energy transport for all parameter values that we investigated. Second, for weak disorder there is a stable diffusive region that persists up to a critical disorder strength that depends on the XY anisotropy. Third, for disorder strengths above this critical value energy transport becomes increasingly subdiffusive. Fourth, the many-body localization transition moves to significantly higher disorder strengths as the XY anisotropy is increased. We discuss these results, and their relation to our current physical picture of subdiffusion in the approach to many-body localization.

For high-precision bosonic Josephson junctions, many-body effects matter

Typical treatments of superconducting or superfluid Josephson junctions rely on mean-field or two-mode models; we explore many-body dynamics of an isolated, ultracold, Bose-gas long Josephson junction using time-evolving block decimation simulations. We demonstrate that with increasing repulsive interaction strength, localized dynamics emerge that influence macroscopic condensate behavior and can lead to the formation of solitons that directly oppose the symmetry of the junction. Initial state population and phase yield insight into dynamic tunneling regimes of a quasi one-dimensional double well potential, from Josephson oscillations to macroscopic-self-trapping. Population imbalance simulations reveal substantial deviation of many-body dynamics from mean-field Gross–Pitaevskii predictions, particularly as the barrier height and interaction strength increase. In addition, the sudden approximation supports localized particle–hole formation after a diabatic quench, and correlation measures unveil a new dynamic regime: the Fock flashlight.

Sharp entanglement thresholds in the logarithmic negativity of disjoint blocks in the transverse-field Ising chain

Entanglement has developed into an essential concept for the characterization of phases and phase transitions in ground states of quantum many-body systems. In this work we use the logarithmic negativity to study the spatial entanglement structure in the transverse-field Ising chain both in the ground state and at nonzero temperatures. Specifically, we investigate the entanglement between two disjoint blocks as a function of their separation, which can be viewed as the entanglement analog of a spatial correlation function. We find sharp entanglement thresholds at a critical distance beyond which the logarithmic negativity vanishes exactly and thus the two blocks become unentangled, which holds even in the presence of long-ranged quantum correlations, i.e., at the system’s quantum critical point. Using time-evolving block decimation, we explore this feature as a function of temperature and size of the two blocks and present a simple model to describe our numerical observations.

Quantum phase transition modulation in an atomtronic Mott switch

Mott insulators provide stable quantum states and long coherence times due to small number fluctuations, making them good candidates for quantum memory and atomic circuits. We propose a proof-of-principle for a 1D Mott switch using an ultracold Bose gas and optical lattice. With time-evolving block decimation simulations—efficient matrix product state methods—we design a means for transient parameter characterization via a local excitation for ease of engineering into more complex atomtronics. We perform the switch operation by tuning the intensity of the optical lattice, and thus the interaction strength through a conductance transition due to the confined modifications of the ‘wedding cake’ Mott structure. We demonstrate the time-dependence of Fock state transmission and fidelity of the excitation as a means of tuning up the device in a double well and as a measure of noise performance. Two-point correlations via the g(2) measure provide additional information regarding superfluid fragments on the Mott insulating background due to the confinement of the potential.

Variational optimization algorithms for uniform matrix product states

We combine the Density Matrix Renormalization Group (DMRG) with Matrix Product State tangent space concepts to construct a variational algorithm for finding ground states of one dimensional quantum lattices in the thermodynamic limit. A careful comparison of this variational uniform Matrix Product State algorithm (VUMPS) with infinite Density Matrix Renormalization Group (IDMRG) and with infinite Time Evolving Block Decimation (ITEBD) reveals substantial gains in convergence speed and precision. We also demonstrate that VUMPS works very efficiently for Hamiltonians with long range interactions and also for the simulation of two dimensional models on infinite cylinders. The new algorithm can be conveniently implemented as an extension of an already existing DMRG implementation.

Ultra-fast relaxation, decoherence, and localization of photoexcited states in π-conjugated polymers.

The exciton relaxation dynamics of photoexcited electronic states in poly(p-phenylenevinylene) are theoretically investigated within a coarse-grained model, in which both the exciton and nuclear degrees of freedom are treated quantum mechanically. The Frenkel-Holstein Hamiltonian is used to describe the strong exciton-phonon coupling present in the system, while external damping of the internal nuclear degrees of freedom is accounted for by a Lindblad master equation. Numerically, the dynamics are computed using the time evolving block decimation and quantum jump trajectory techniques. The values of the model parameters physically relevant to polymer systems naturally lead to a separation of time scales, with the ultra-fast dynamics corresponding to energy transfer from the exciton to the internal phonon modes (i.e., the C-C bond oscillations), while the longer time dynamics correspond to damping of these phonon modes by the external dissipation. Associated with these time scales, we investigate the following processes that are indicative of the system relaxing onto the emissive chromophores of the polymer: (1) Exciton-polaron formation occurs on an ultra-fast time scale, with the associated exciton-phonon correlations present within half a vibrational time period of the C-C bond oscillations. (2) Exciton decoherence is driven by the decay in the vibrational overlaps associated with exciton-polaron formation, occurring on the same time scale. (3) Exciton density localization is driven by the external dissipation, arising from "wavefunction collapse" occurring as a result of the system-environment interactions. Finally, we show how fluorescence anisotropy measurements can be used to investigate the exciton decoherence process during the relaxation dynamics.

Probabilistic low-rank factorization accelerates tensor network simulations of critical quantum many-body ground states.

We provide evidence that randomized low-rank factorization is a powerful tool for the determination of the ground-state properties of low-dimensional lattice Hamiltonians through tensor network techniques. In particular, we show that randomized matrix factorization outperforms truncated singular value decomposition based on state-of-the-art deterministic routines in time-evolving block decimation (TEBD)- and density matrix renormalization group (DMRG)-style simulations, even when the system under study gets close to a phase transition: We report linear speedups in the bond or local dimension of up to 24 times in quasi-two-dimensional cylindrical systems.

Topological orders and quantum phase transitions in a one-dimensional extended quantum compass model

By using the infinite time evolving block decimation in the presentation of infinite matrix product states, we study an extended quantum compass model (EQCM). This model does not only include extremely rich phase diagrams due to competitions of orbital degrees of freedom and anisotropic couplings between pseudospin-1/2 operators but also have the capacity to describe property of protected qubits for quantum computation which leads to lots of attentions paid to the phase boundaries of the EQCM. However, few attentions are paid to long-range topological string correlation order parameters of the EQCM. To study order parameters, one should understand spontaneous symmetry breaking which relates to Landau quantum phase transitions theory. Once spontaneous symmetry breaking happens, there should exist some local order which can be described by a local order parameter. This order parameter can be used to distinguish the phase from others. For continuous quantum phase transitions, in the critical regime, critical exponents can be extracted. Unfortunately, the long-range topological string correlation orders are beyond Landau quantum phase transitions theory, one can not directly use two paradigms of Landau-Ginzburg-Wilson. Usually, one can define a local order parameter by local magnetization. Naturally, one can also refer to this way to define the long-range topological string correlation order parameters by long-range topological string correlations on the following conditions, i.e. the quantum system undergoes a hidden spontaneous symmetry breaking; the long-range topological string correlation order parameter can be used to distinguish the phase from others; for continuous quantum phase transitions, the long-range topological string correlation order parameter satisfies scaling law when control parameter getting close to critical point. Based on above idea, in order to characterize the topological ordered phases and quantum phase transitions in the EQCM, even/odd long-range topological string correlations are introduced based on even/odd bonds. Hereafter, fidelity per lattice site, even/odd long-range topological string correlations, the saturation behavior of odd long-range topological string correlations and order parameters are calculated. The long-range topological string correlations show three distinguished behaviors which include decaying to zero, monotonic saturation and oscillatory saturation. By the above characterizations, oscillatory/monotonic odd long-range topological string correlation order parameter is derived. Then ground-state phase diagram of order parameters is computed which includes oscillatory/monotonic odd long-range topological string correlation order phase and antiferromagnetic phase. In the critical regime, critical exponent β=1/8 extracted from monotonic odd long-range topological string correlation order parameter and local magnetization shows the phase transition belongs to Ising universality. In addition, the phase transition points, the order of the phase transitions of fidelity show consistent with the results of order parameters.

Macroscopic quantum tunneling escape of Bose-Einstein condensates

Recent experiments on macroscopic quantum tunneling reveal a non-exponential decay of the number of atoms trapped in a quasibound state behind a potential barrier. Through both experiment and theory, we demonstrate this non-exponential decay results from interactions between atoms. Quantum tunneling of tens of thousands of 87 Rb atoms in a Bose-Einstein condensate is modeled by a modified Jeffreys-Wentzel-Kramers-Brillouin model, taking into account the effective time-dependent barrier induced by the mean-field. Three-dimensional Gross-Pitaevskii simulations corroborate a mean-field result when compared with experiments. However, with one-dimensional modeling using time-evolving block decimation, we present an effective renormalized mean-field theory that suggests many-body dynamics for which a bare mean-field theory may not apply.

Quantum phase transition in the one-dimensional quantum Heisenberg XYZ model with Dzyaloshinskii–Moriya interaction

We investigated the quantum phase transition occurred in one-dimensional quantum Heisenberg XYZ model with Dzyaloshinskii–Moriya interaction via the infinite matrix product state representation with the infinite time evolving block decimation method. Entanglement entropy and local order parameter in and near the transition point are given. Scaling relation plays crucial roles on identifying a quantum system with a physically different phase. The scaling relation of the entanglement entropy, local order parameter and finite correlation length with the truncation dimension are also obtained. All the interesting results give a theoretical justification for the high accuracy of infinite time evolved block decimation algorithm which works in the thermodynamical limit.

Hybrid Semiclassical Theory of Quantum Quenches in One-Dimensional Systems.

We develop a hybrid semiclassical method to study the time evolution of one-dimensional quantum systems in and out of equilibrium. Our method handles internal degrees of freedom completely quantum mechanically by a modified time-evolving block decimation method while treating orbital quasiparticle motion classically. We can follow dynamics up to time scales well beyond the reach of standard numerical methods to observe the crossover between preequilibrated and locally phase equilibrated states. As an application, we investigate the quench dynamics and phase fluctuations of a pair of tunnel-coupled one-dimensional Bose condensates. We demonstrate the emergence of soliton-collision-induced phase propagation, soliton-entropy production, and multistep thermalization. Our method can be applied to a wide range of gapped one-dimensional systems.

Photoinduced charge-order melting dynamics in a one-dimensional interacting Holstein model

Transient quantum dynamics in an interacting fermion-phonon system are investigated. In particular, a charge order (CO) melting after a short optical-pulse irradiation and roles of the quantum phonons on the transient dynamics are focused on. A spinless-fermion model in a one-dimensional chain coupled with local phonons is analyzed numerically. The infinite time-evolving block decimation algorithm is adopted as a reliable numerical method for one-dimensional quantum many-body systems. Numerical results for the photoinduced CO melting dynamics without phonons are well interpreted by the soliton picture for the CO domains. This interpretation is confirmed by the numerical simulation for an artificial local excitation and the classical soliton model. In the case of the large phonon frequency corresponding to the antiadiabatic condition, the CO melting is induced by propagations of the polaronic solitons with the renormalized soliton velocity. On the other hand, in the case of the small phonon frequency corresponding to the adiabatic condition, the first stage of the CO melting dynamics occurs due to the energy transfer from the fermionic to phononic systems, and the second stage is brought about by the soliton motions around the bottom of the soliton band. Present analyses provide a standard reference for the photoinduced CO melting dynamics in low-dimensional many-body quantum systems.

Impurity and soliton dynamics in a Fermi gas with nearest-neighbor interactions

We study spinless fermions with repulsive nearest-neighbor interactions perturbed by an impurity particle or a local potential quench. Using the numerical time-evolving block decimation method and a simplified analytic model, we show that the pertubations create a soliton-antisoliton pair. If solitons are already present in the bath, the two excitations have a drastically different dynamics: The antisoliton does not annihilate with the solitons and is therefore confined close to its origin while the soliton excitation propagates. We discuss the consequences for experiments with ultracold gases.

Entangled Dynamics in Macroscopic Quantum Tunneling of Bose-Einstein Condensates.

Tunneling of a quasibound state is a nonsmooth process in the entangled many-body case. Using time-evolving block decimation, we show that repulsive (attractive) interactions speed up (slow down) tunneling. While the escape time scales exponentially with small interactions, the maximization time of the von Neumann entanglement entropy between the remaining quasibound and escaped atoms scales quadratically. Stronger interactions require higher-order corrections. Entanglement entropy is maximized when about half the atoms have escaped.

Driven Bose-Hubbard model with a parametrically modulated harmonic trap

We investigate a one-dimensional Bose-Hubbard model in a parametrically driven global harmonic trap. The delicate interplay of both the local atom interaction and the global driving allows to control the dynamical stability of the trapped quantum many-body state. The mechanism is illustrated for weak interaction by a discretized Gross-Pitaevskii equation within a Gaussian variational ansatz, yielding to a Mathieu equation for the condensate width. The parametric resonance condition can be tuned by the atom interaction strength. For stronger interaction, this mechanism is confirmed by results of the numerically exact time-evolving block decimation scheme. The global modulation also induces an effective time-independent inhomogeneous hopping strength for the atoms.

Entanglement dynamics in critical random quantum Ising chain with perturbations

We simulate the entanglement dynamics in a critical random quantum Ising chain with generic perturbations using the time-evolving block decimation algorithm. Starting from a product state, we observe super-logarithmic growth of entanglement entropy with time. The numerical result is consistent with the analytical prediction of Vosk and Altman using a real-space renormalization group technique.

The Approximation of Bosonic System by Fermion in Quantum Cellular Automaton

In one-dimensional multiparticle Quantum Cellular Automaton (QCA), the approximation of the bosonic system by fermion (boson-fermion correspondence) can be derived in a rather simple and intriguing way, where the principle to impose zero-derivative boundary conditions of one-particle QCA is also analogously used in particle-exchange boundary conditions. As a clear cut demonstration of this approximation, we calculate the ground state of few-particle systems in a box using imaginary time evolution simulation in 2nd quantization form as well as in 1st quantization form. Moreover in this 2nd quantized form of QCA calculation, we use Time Evolving Block Decimation (TEBD) algorithm. We present this demonstration to emphasize that the TEBD is most natu-rally regarded as an approximation method to the 2nd quantized form of QCA.

Quantum entanglement and critical exponents in one-dimensional spin-1 bond-alternating XXZ chains

The characterization of the quantum phase transition in a lowdimensional system has attracted a considerable amount of attention in quantum manybody systems. As one of the fundamental models in quantum magnetism, spin-1 models have richer phase diagrams and show more complex physical phenomena. In the spin-1 antiferromagnetic XXZ model, the Haldane phase and the Nel phase are the gapped topologic phases which cannot be characterized by the local order parameters. To characterize the nature in such phases, one has to calculate the non-local long range order parameters. Normally, the non-local order parameter in the topological phase is obtained from the extrapolation of finite-sized system in numerical study. However, it is difficult to extract the critical exponents with such an extrapolated non-local order parameter due to the numerical accuracy. In a recently developed tensor network representation, i.e., the infinite matrix product state (iMPS) algorithm, it was shown that the non-local order can be directly calculated from a very large lattice distance in an infinite-sized system rather than an extrapolated order parameter in a finite-sized system. Therefore, it is worthwhile using this convenient technique to study the non-local orders in the topological phases and characterize the quantum criticalities in the topological quantum phase transitions. In this paper, by utilizing the infinite matrix product state algorithm based on the tensor network representation and infinite time evolving block decimation method, the quantum entanglement, fidelity, and critical exponents of the topological phase transition are investigated in the one-dimensional infinite spin-1 bond-alternating XXZ Heisenberg model. It is found that there is always a local dimerization order existing in the whole parameter range when the bond-alternative strength parameter changes from 0 to 1. Also, due to the effect of the bond-alternating, there appears a quantum phase transition from the long-rang ordering topological Nel phase to the local ordering dimerization phase. The von Neumann entropy, fidelity per lattice site, and order parameters all give the same phase transition point at c = 0.638. To identify the type of quantum phase transition, the central charge c 0.5 is extracted from the ground state von Neumann entropy and the finite correlation length, which indicates that the phase transition belongs to the two-dimensional Ising universality class. Furthermore, it is found that the Nel order and the susceptibility of Nel order have power-law relations to |-c|. From the numerical fitting of the Nel order and its susceptibility, we obtain the characteristic critical exponents ' = 0.236 and ' = 0.838. It indicates that such critical exponents from our method characterize the nature of the quantum phase transition. Our critical exponents from the iMPS method can provide guidance for studying the properties of the phase transition in quantum spin systems.