Real-time dynamics in 2+1D compact QED using complex periodic Gaussian states

We investigate tensor-train approaches to the solution of the time-dependent Schrödinger equation for chain-like quantum systems with on-site and nearest-neighbor interactions only. Using efficient low-rank tensor train representations, we aim at reducing memory consumption and computational costs. As an example, coupled excitons and phonons modeled in terms of Fröhlich-Holstein type Hamiltonians are studied here. By comparing our tensor-train-based results with semi-analytical results, we demonstrate the key role of the ranks of the quantum state vectors. Typically, an excellent quality of solutions is found only when the maximum number of ranks exceeds a certain value. One class of propagation schemes builds on splitting the Hamiltonian into two groups of interleaved nearest-neighbor interactions commutating within each of the groups. In particular, the fourth-order Yoshida-Neri and the eighth-order Kahan-Li symplectic composition yield results close to machine precision. Similar results are found for fourth and eighth order global Krylov scheme. However, the computational effort currently restricts the use of these four propagators to rather short chains, which also applies to propagators based on the time-dependent variational principle, typically used for matrix product states. Yet, another class of propagators involves explicit, time-symmetrized Euler integrators. Especially, the fourth-order variant is recommended for quantum simulations of longer chains, even though the high precision of the splitting schemes cannot be reached. Moreover, the scaling of the computational effort with the dimensions of the local Hilbert spaces is much more favorable for the differencing than for splitting or variational schemes.

We construct a random unitary Gaussian circuit for continuous-variable (CV) systems subject to Gaussian measurements. We show that when the measurement rate is nonzero, the steady state entanglement entropy saturates to an area-law scaling. This is different from a many-body qubit system, where a generic entanglement transition is widely expected. Due to the unbounded local Hilbert space, the time scale to destroy entanglement is always much shorter than the one to build it, while a balance could be achieved for a finite local Hilbert space. By the same reasoning, the absence of transition should also hold for other non-unitary Gaussian CV dynamics.

We investigate the ground state of interacting spin-$\half$ fermions (3D) at a finite density ($\rho \sim \kf^3$) in the presence of a uniform non-Abelian gauge field. The gauge field configuration (GFC) described by a vector $\blam \equiv (\lambda_x, \lambda_y, \lambda_z)$, whose magnitude $\lambda$ determines the gauge coupling strength, generates a generalized Rashba spin-orbit interaction. For a weak attractive interaction in the singlet channel described by a small negative scattering length $(\kf |\as| \lesssim 1)$, the ground state in the absence of the gauge field ($\lambda=0$) is a BCS (Bardeen-Cooper-Schrieffer) superfluid with large overlapping pairs. With increasing gauge coupling strength, a non-Abelian gauge field engenders a crossover of this BCS ground state to a BEC (Bose-Einstein condensate) ground state of bosons even with a weak attractive interaction that fails to produce a two-body bound state in free vacuum. For large gauge couplings $(\lambda/\kf \gg 1)$, the BEC attained is a condensate of bosons whose properties are solely determined by the gauge field (and not by the scattering length so long as it is non-zero) -- we call these bosons "rashbons". In the absence of interactions ($\as = 0^-$), the shape of the Fermi surface of the system undergoes a topological transition at a critical gauge coupling $\lambda_T$. For high symmetry gauge field configurations we show that the crossover from the BCS superfluid to the rashbon BEC occurs in the regime of $\lambda$ near $\lambda_T$. In the context of cold atomic systems, this work makes an interesting suggestion of obtaining BCS-BEC crossover through a route other than tuning the interaction between the fermions.

Characterization of the ground states of a simplified CuO2 model in strong- and weak-coupling regions

Exact classical solutions of the massless σ-model with gauge fields in minkowski space

Confinement is an ubiquitous phenomenon when matter couples to gauge fields, which manifests itself in a linear string potential between two static charges. Although gauge fields can be integrated out in one dimension, they can mediate nonlocal interactions which in turn influence the paradigmatic Luttinger liquid properties. However, when the charges become dynamical and their densities finite, understanding confinement becomes challenging. Here we show that confinement in 1D Z_{2} lattice gauge theories, with dynamical matter fields and arbitrary densities, is related to translational symmetry breaking in a nonlocal basis. The exact transformation to this string-length basis leads us to an exact mapping of Luttinger parameters reminiscent of a Luther-Emery rescaling. We include the effects of local, but beyond contact, interactions between the matter particles, and show that confined mesons can form a Mott-insulating state when the deconfined charges cannot. While the transition to the Mott state cannot be detected in the Green's function of the charges, we show that the metallic state is characterized by hidden off-diagonal quasi-long-range order. Our predictions provide new insights to the physics of confinement of dynamical charges, and can be experimentally addressed in Rydberg-dressed quantum gases in optical lattices.

We present the exact Bethe ansatz solution for the two-dimensional BCS pairing Hamiltonian with ${p}_{x}+i{p}_{y}$ symmetry. Using both mean-field theory and the exact solution we obtain the ground-state phase diagram parametrized by the filling fraction and the coupling constant. It consists of three phases that are denoted weak-coupling BCS, weak pairing, and strong pairing. The first two phases are separated by a topologically protected line where the exact ground state is given by the Moore-Read pfaffian state. In the thermodynamic limit the ground-state energy is discontinuous on this line. The other two phases are separated by the critical line, also topologically protected, previously found by Read and Green. We establish a duality relation between the weak and strong pairing phases, whereby ground states of the weak phase are ``dressed'' versions of the ground states of the strong phase by zero energy (Moore-Read) pairs and characterized by a topological order parameter.

It is well known that symmetry-protected topological (SPT) phases host nontrivial boundaries that cannot be mimicked in a lower-dimensional system with a conventional realization of symmetry. However, for SPT phases of bosons (fermions) within the cohomology (supercohomology) classification the boundary can be recreated without the bulk at the cost of a non-onsite-symmetry action. This raises the following question: Can one also mimic the boundaries of SPT phases which lie outside the (super)cohomology classification? In this paper, we study this question in the context of $(2+1)$-dimensional $[(2+1)\mathrm{D}]$ fermion SPTs. We focus on the root SPT phase for the symmetry group $G={Z}_{2}\ifmmode\times\else\texttimes\fi{}{Z}_{2}^{f}$. Starting with an exactly solvable model for the bulk of this phase constructed by Tarantino and Fidkowski, we derive an effective one-dimensional (1D) lattice model for the boundary. Crucially, the Hilbert space of this 1D model does not have a local tensor product structure, but rather is obtained by placing a local constraint on a local tensor product Hilbert space. We derive the action of the ${Z}_{2}$ symmetry on this Hilbert space and find a simple three-site Hamiltonian that respects this symmetry. We study this Hamiltonian numerically using exact diagonalization and DMRG and find strong evidence that it realizes an Ising conformal field theory where the ${Z}_{2}$ symmetry acts as the Kramers-Wannier duality; this is the expected stable gapless boundary state of the present SPT. A simple modification of our construction realizes the boundary of the $(2+1)\mathrm{D}$ topological superconductor protected by time-reversal symmetry $\mathcal{T}$ with ${\mathcal{T}}^{2}={(\ensuremath{-}1)}^{\mathcal{F}}$.

Using a simple ansatz, Yang’s K gauge equation for SU(2) gauge fields is reduced to a system of nonlinear ordinary differential equations. Exact solutions for the equations are obtained together with corresponding gauge potentials.

Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born-Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For the intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.

Entanglement is the fundamental difference between classical and quantum\nsystems and has become one of the guiding principles in the exploration of\nhigh- and low-energy physics. The calculation of entanglement entropies in\ninteracting quantum field theories, however, remains challenging. Here, we\npresent the first method for the explicit computation of reduced density\nmatrices of interacting quantum field theories using truncated Hamiltonian\nmethods. The method is based on constructing an isomorphism between the Hilbert\nspace of the full system and the tensor product of Hilbert spaces of\nsub-intervals. This naturally enables the computation of the von Neumann and\narbitrary R\\'enyi entanglement entropies as well as mutual information,\nlogarithmic negativity and other measures of entanglement. Our method is\napplicable to equilibrium states and non-equilibrium evolution in real time. It\nis model independent and can be applied to any Hamiltonian truncation method\nthat uses a free basis expansion. We benchmark the method on the free\nKlein-Gordon theory finding excellent agreement with the analytic results. We\nfurther demonstrate its potential on the interacting sine-Gordon model,\nstudying the scaling of von Neumann entropy in ground states and real time\ndynamics following quenches of the model.\n

We reformulate the Thirring model in $D$ $(2 \le D < 4)$ dimensions as a gauge theory by introducing $U(1)$ hidden local symmetry (HLS) and study the dynamical mass generation of the fermion through the Schwinger-Dyson (SD) equation. By virtue of such a gauge symmetry we can greatly simplify the analysis of the SD equation by taking the most appropriate gauge (``nonlocal gauge'') for the HLS. In the case of even-number of (2-component) fermions, we find the dynamical fermion mass generation as the second order phase transition at certain fermion number, which breaks the chiral symmetry but preserves the parity in (2+1) dimensions ($D=3$). In the infinite four-fermion coupling (massless gauge boson) limit in (2+1) dimensions, the result coincides with that of the (2+1)-dimensional QED, with the critical number of the 4-component fermion being $N_{\rm cr} = \frac{128}{3\pi^{2}}$. As to the case of odd-number (2-component) fermion in (2+1) dimensions, the regularization ambiguity on the induced Chern-Simons term may be resolved by specifying the regularization so as to preserve the HLS. Our method also applies to the (1+1) dimensions, the result being consistent with the exact solution. The bosonization mechanism in (1+1) dimensional Thirring model is also reproduced in the context of dual-transformed theory for the HLS.

The critical value of the atom-field coupling strength for a finite number of\natoms is deter- mined by means of both, semiclassical and exact solutions. In\nthe semiclassical approach we use a variational procedure with coherent and\nsymmetry-adapted states, while for the exact quantum solution the concept of\nfidelity is employed. These procedures allow for the determination of the phase\ntransitions in the model, and coincide in the thermodynamic limit. For the\nthree cases men- tioned above, universal parametric curves are obtained for the\nexpectation values of both the first quadrature of the electromagnetic field,\nand the atomic relative population, as implicit functions of the atom-field\ncoupling parameter, valid for the ground- and first-excited states.\n

In this paper we extend our previously discovered exact solution for an SU(2) Yang-Mills-Higgs theory, to the general group SU(N+1). Using the first-order formalism of Bogomolny, an exact, spherically symmetric solution for the gauge and scalar fields is found. This solution is similar to the Schwarzschild solution of general relativity, in that the gauge and scalar fields become infinite on a spherical shell of radiusr 0=K. However in the Schwarzschild case the singularity at the event horizon is a coordinate singularity while for the present solution the singularity is a true singularity. It is speculated that this solution may give a confinement mechanism for non-Abelian gauge theories, since any particle which carries the SU(N+1) charge would become permanently trapped inside the regionr<r 0.

In this paper we extend our previously discovered exact solution for an SU(2) Yang-Mills-Higgs theory, to the general group SU(N+1). Using the first-order formalism of Bogomolny, an exact, spherically symmetric solution for the gauge and scalar fields is found. This solution is similar to the Schwarzschild solution of general relativity, in that the gauge and scalar fields become infinite on a spherical shell of radiusr 0=K. However in the Schwarzschild case the singularity at the event horizon is a coordinate singularity while for the present solution the singularity is a true singularity. It is speculated that this solution may give a confinement mechanism for non-Abelian gauge theories, since any particle which carries the SU(N+1) charge would become permanently trapped inside the regionr<r 0.