Search

While generative models have shown great success in generating high-dimensional samples conditional on low-dimensional descriptors (learning e.g. stroke thickness in MNIST, hair color in CelebA, or speaker identity in Wavenet), their generation out-of-sample poses fundamental problems. The conditional variational autoencoder (CVAE) as a simple conditional generative model does not explicitly relate conditions during training and, hence, has no incentive of learning a compact joint distribution across conditions. We overcome this limitation by matching their distributions using maximum mean discrepancy (MMD) in the decoder layer that follows the bottleneck. This introduces a strong regularization both for reconstructing samples within the same condition and for transforming samples across conditions, resulting in much improved generalization. We refer to the architecture as \emph{transformer} VAE (trVAE). Benchmarking trVAE on high-dimensional image and tabular data, we demonstrate higher robustness and higher accuracy than existing approaches. In particular, we show qualitatively improved predictions for cellular perturbation response to treatment and disease based on high-dimensional single-cell gene expression data, by tackling previously problematic minority classes and multiple conditions. For generic tasks, we improve Pearson correlations of high-dimensional estimated means and variances with their ground truths from 0.89 to 0.97 and 0.75 to 0.87, respectively., Comment: Added reference to Johansson et al. (2016) and removed sentences from Lopez et al. (2018) in the background section (see acknowledgements)

We present a new impurity solver for dynamical mean-field theory based on imaginary-time evolution of matrix product states. This converges the self-consistency loop on the imaginary-frequency axis and obtains real-frequency information in a final real-time evolution. Relative to computations on the real-frequency axis, required bath sizes are much smaller and less entanglement is generated, so much larger systems can be studied. The power of the method is demonstrated by solutions of a three band model in the single and two-site dynamical mean-field approximation. Technical issues are discussed, including details of the method, efficiency as compared to other matrix product state based impurity solvers, bath construction and its relation to real-frequency computations and the analytic continuation problem of quantum Monte Carlo, the choice of basis in dynamical cluster approximation, and perspectives for off-diagonal hybridization functions., Comment: 8 pages + 4 pages appendix, 9 figures

Many numerical techniques for the description of quantum systems that are coupled to a continuous bath require the discretization of the latter. To this end, a wealth of methods has been developed in the literature, which we classify as (i) direct discretization, (ii) orthogonal polynomial, and (iii) numerical optimization strategies. We recapitulate strategies (i) and (ii) to clarify their relation. For quadratic Hamiltonians, we show that (ii) is the best strategy in the sense that it gives the numerically exact time evolution up to a maximum time $t_\text{max}$, for which we give a simple expression. For non-quadratic Hamiltonians, we show that no such best strategy exists. We present numerical examples relevant to open quantum systems and strongly correlated systems, as treated by dynamical mean-field theory (DMFT).

We study the unitary time evolution of antiferromagnetic order in the Hubbard model after a quench starting from the perfect N\'eel state. In this setup, which is well suited for experiments with cold atoms, one can distinguish fundamentally different pathways for melting of long-range order at weak and strong interaction. In the Mott insulating regime, melting of long-range order occurs due to the ultra-fast transfer of energy from charge excitations to the spin background, while local magnetic moments and their exchange coupling persist during the process. The latter can be demonstrated by a local spin-precession experiment. At weak interaction, local moments decay along with the long-range order. The dynamics is governed by residual quasiparticles, which are reflected in oscillations of the off-diagonal components of the momentum distribution. Such oscillations provide an alternative route to study the prethermalization phenomenon and its influence on the dynamics away from the integrable (noninteracting) limit. The Hubbard model is solved within nonequilibrium dynamical mean-field theory, using the density matrix-renormalization group as an impurity solver., Comment: 11 pages

We link linear prediction of Chebyshev and Fourier expansions to analytic continuation. We push the resolution in the Chebyshev-based computation of $T=0$ many-body spectral functions to a much higher precision by deriving a modified Chebyshev series expansion that allows to reduce the expansion order by a factor $\sim\frac{1}{6}$. We show that in a certain limit the Chebyshev technique becomes equivalent to computing spectral functions via time evolution and subsequent Fourier transform. This introduces a novel recursive time evolution algorithm that instead of the group operator $e^{-iHt}$ only involves the action of the generator $H$. For quantum impurity problems, we introduce an adapted discretization scheme for the bath spectral function. We discuss the relevance of these results for matrix product state (MPS) based DMRG-type algorithms, and their use within dynamical mean-field theory (DMFT). We present strong evidence that the Chebyshev recursion extracts less spectral information from $H$ than time evolution algorithms when fixing a given amount of created entanglement., Comment: 12 pages + 6 pages appendix, 11 figures

We introduce a strictly single-site DMRG algorithm based on the subspace expansion of the Alternating Minimal Energy (AMEn) method. The proposed new MPS basis enrichment method is sufficient to avoid local minima during the optimisation, similarly to the density matrix perturbation method, but computationally cheaper. Each application of $\hat H$ to $|\Psi\rangle$ in the central eigensolver is reduced in cost for a speed-up of $\approx (d + 1)/2$, with $d$ the physical site dimension. Further speed-ups result from cheaper auxiliary calculations and an often greatly improved convergence behaviour. Runtime to convergence improves by up to a factor of 2.5 on the Fermi-Hubbard model compared to the previous single-site method and by up to a factor of 3.9 compared to two-site DMRG. The method is compatible with real-space parallelisation and non-abelian symmetries., Comment: 9 pages, 6 figures; added comparison with two-site DMRG

We solve the nonequilibrium dynamical mean-field theory (DMFT) using matrix product states (MPS). This allows us to treat much larger bath sizes and by that reach substantially longer times (factor $\sim$ 2 -- 3) than with exact diagonalization. We show that the star geometry of the underlying impurity problem can have substantially better entanglement properties than the previously favoured chain geometry. This has immense consequences for the efficiency of an MPS-based description of general impurity problems: in the case of equilibrium DMFT, it leads to an orders-of-magnitude speedup. We introduce an approximation for the two-time hybridization function that uses time-translational invariance, which can be observed after a certain relaxation time after a quench to a time-independent Hamiltonian., Comment: 11 pages + 3 pages appendix, 14 figures

We compute the spectral functions for the two-site dynamical cluster theory and for the two-orbital dynamical mean-field theory in the density-matrix renormalization group (DMRG) framework using Chebyshev expansions represented with matrix product states (MPS). We obtain quantitatively precise results at modest computational effort through technical improvements regarding the truncation scheme and the Chebyshev rescaling procedure. We furthermore establish the relation of the Chebyshev iteration to real-time evolution, and discuss technical aspects as computation time and implementation in detail., Comment: 20 pages, 12 figures

We study the quantum Rabi model within the framework of the analytical solution developed in Phys. Rev. Lett. 107,100401 (2011). In particular, through time-dependent correlation functions, we give a quantitative criterion for classifying two regions of the quantum Rabi model, involving the Jaynes-Cummings, the ultrastrong, and deep strong coupling regimes. In addition, we find a stationary qubit-field entangled basis that governs the whole dynamics as the coupling strength overcomes the mode frequency., Comment: 8 pages, 8 figures. Revised version, accepted for publication in Physical Review A

We use the analytical solution of the quantum Rabi model to obtain absolutely convergent series expressions of the exact eigenstates and their scalar products with Fock states. This enables us to calculate the numerically exact time evolution of <\sigma_x(t)> and <\sigma_z(t)> for all regimes of the coupling strength, without truncation of the Hilbert space. We find a qualitatively different behavior of both observables which can be related to their representations in the invariant parity subspaces., Comment: 8 pages, 7 figures, published version

We study the expansion of bosonic Mott insulators in the presence of an optical lattice after switching off a confining potential. We use the Gutzwiller mean-field approximation and consider two different setups. In the first one, the expansion is restricted to one direction. We show that this leads to the emergence of two condensates with well defined momenta, and argue that such a construct can be used to create atom lasers in optical lattices. In the second setup, we study Mott insulators that are allowed to expand in all directions in the lattice. In this case, a simple condensate is seen to develop within the mean-field approximation. However, its constituent bosons are found to populate many nonzero momentum modes. An analytic understanding of both phenomena in terms of the exact dispersion relation in the hard-core limit is presented., Comment: 11 pages, 9 figures. Figures 2,3,4 corrected

A quantum many-body system which is prepared in the ground state of an integrable Hamiltonian does not directly thermalize after a sudden small parameter quench away from integrability. Rather, it will be trapped in a prethermalized state and can thermalize only at a later stage. We discuss several examples for which this prethermalized state shares some properties with the nonthermal steady state that emerges in the corresponding integrable system. These examples support the notion that nonthermal steady states in integrable systems may be viewed as prethermalized states that never decay further. Furthermore we show that prethermalization plateaus are under certain conditions correctly predicted by generalized Gibbs ensembles, which are the appropriate extension of standard statistical mechanics in the presence of many constants of motion. This establishes that the relaxation behaviors of integrable and nearly integrable systems are continuously connected and described by the same statistical theory., Comment: 11 pages, 2 figures

We present a theoretical study of the quantum corrections to the revival time due to finite tunneling in the collapse and revival of matter wave interference after a quantum quench. We study hard-core bosons in a superlattice potential and the Bose-Hubbard model by means of exact numerical approaches and mean-field theory. We consider systems without and with a trapping potential present. We show that the quantum corrections to the revival time can be used to accurately determine the value of the hopping parameter in experiments with ultracold bosons in optical lattices., Comment: 10 pages, 12 figures, typos in section 3A corrected