BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Martin Ganahl

Martin Ganahl appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elquant-phcond-mat.mes-hallcond-mat.quant-gascond-mat.dis-nnphysics.comp-phcond-mat.stat-mechDistributed, Parallel, and Cluster Computing

Catalog footprint

What is connected

13works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Simulation of quantum many-body dynamics with Tensor Processing Units: Floquet prethermalization

Tensor Processing Units (TPUs) are specialized hardware accelerators developed by Google to support large-scale machine-learning tasks, but they can also be leveraged to accelerate and scale other linear-algebra-intensive computations. In this paper we demonstrate the usage of TPUs for massively parallel, classical simulations of quantum many-body dynamics on long timescales. We apply our methods to study the phenomenon of Floquet prethermalization, i.e., exponentially slow heating in quantum spin chains subject to high-frequency periodic driving. We simulate the dynamics of L=34 qubits for over $10^5$ Floquet periods, corresponding to circuits with millions of two-qubit gates. The circuits simulated have no additional symmetries and represent a pure-state evolution in the full $2^L$-dimensional Hilbert space. This is achieved by distributing the computation over 128 TPU cores. On that size TPU cluster, we find speedups in wall-clock runtime of 230x and 15x when compared to reference CPU and single-GPU simulations, respectively, for shorter 30-qubit simulations that can be handled by all three platforms. We study the computational cost of the simulations, as a function of both the number of qubits and the number of TPU cores used, up to our maximum capacity of L=40 qubits, which requires a ``full pod" of 2048 TPU cores with tens of terabytes of memory in total. For these simulations, an 8-TPU-core machine is comparable to a single A100 GPU, and thus the full TPU pod is comparable to a machine with hundreds of GPUs. However, the TPU pod is more energy and cost efficient, and readily accessible (via Google Cloud), unlike such large many-GPU configurations. We also study the accumulation of numerical error as a function of circuit depth in very deep circuits. Our work demonstrates that TPUs can offer significant advantages for state-of-the-art simulations of quantum many-body dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large Scale Distributed Linear Algebra With Tensor Processing Units

We have repurposed Google Tensor Processing Units (TPUs), application-specific chips developed for machine learning, into large-scale dense linear algebra supercomputers. The TPUs' fast inter-core interconnects (ICI)s, physically two-dimensional network topology, and high-bandwidth memory (HBM) permit distributed matrix multiplication algorithms to rapidly become computationally bound. In this regime, the matrix-multiply units (MXU)s dominate the runtime, yielding impressive scaling, performance, and raw size: operating in float32 precision, a full 2048-core pod of third generation TPUs can multiply two matrices with linear size $N= 220= 1 048 576$ in about 2 minutes. Via curated algorithms emphasizing large, single-core matrix multiplications, other tasks in dense linear algebra can similarly scale. As examples, we present (i) QR decomposition; (ii) resolution of linear systems; and (iii) the computation of matrix functions by polynomial iteration, demonstrated by the matrix polar factorization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Recurrent Neural Network Wave Functions

A core technology that has emerged from the artificial intelligence revolution is the recurrent neural network (RNN). Its unique sequence-based architecture provides a tractable likelihood estimate with stable training paradigms, a combination that has precipitated many spectacular advances in natural language processing and neural machine translation. This architecture also makes a good candidate for a variational wave function, where the RNN parameters are tuned to learn the approximate ground state of a quantum Hamiltonian. In this paper, we demonstrate the ability of RNNs to represent several many-body wave functions, optimizing the variational parameters using a stochastic approach. Among other attractive features of these variational wave functions, their autoregressive nature allows for the efficient calculation of physical estimators by providing independent samples. We demonstrate the effectiveness of RNN wave functions by calculating ground state energies, correlation functions, and entanglement entropies for several quantum spin models of interest to condensed matter physicists in one and two spatial dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Auxiliary master equation approach within matrix product states: Spectral properties of the nonequilibrium Anderson impurity model

Within the recently introduced auxiliary master equation approach it is possible to address steady state properties of strongly correlated impurity models, small molecules or clusters efficiently and with high accuracy. It is particularly suited for dynamical mean field theory in the nonequilibrium as well as in the equilibrium case. The method is based on the solution of an auxiliary open quantum system, which can be made quickly equivalent to the original impurity problem. In its first implementation a Krylov space method was employed. Here, we aim at extending the capabilities of the approach by adopting matrix product states for the solution of the corresponding auxiliary quantum master equation. This allows for a drastic increase in accuracy and permits us to access the Kondo regime for large values of the interaction. In particular, we investigate the nonequilibrium steady state of a single impurity Anderson model and focus on the spectral properties for temperatures T below the Kondo temperature T_K and for small bias voltages phi. For the two cases considered, with T=T_K/4 and T=T_K/10 we find a clear splitting of the Kondo resonance into a two-peak structure for phi close above T_K. In the equilibrium case (phi=0) and for T=T_K/4, the obtained spectral function essentially coincides with the one from numerical renormalization group.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient DMFT impurity solver using real-time dynamics with Matrix Product States

We propose to calculate spectral functions of quantum impurity models using the Time Evolving Block Decimation (TEBD) for Matrix Product States. The resolution of the spectral function is improved by a so-called linear prediction approach. We apply the method as an impurity solver within the Dynamical Mean Field Theory (DMFT) for the single- and two-band Hubbard model on the Bethe lattice. For the single-band model we observe sharp features at the inner edges of the Hubbard bands. A finite size scaling shows that they remain present in the thermodynamic limit. We analyze the real time-dependence of the double occupation after adding a single electron and observe oscillations at the same energy as the sharp feature in the Hubbard band, indicating a long-lived coherent superposition of states that correspond to the Kondo peak and the side peaks. For a two-band Hubbard model we observe an even richer structure in the Hubbard bands, which cannot be related to a multiplet structure of the impurity, in addition to sharp excitations at the band edges of a type similar to the single-band case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lieb-Liniger model with exponentially-decaying interactions: a continuous matrix product state study

The Lieb-Liniger model describes one-dimensional bosons interacting through a repulsive contact potential. In this work, we introduce an extended version of this model by replacing the contact potential with a decaying exponential. Using the recently developed continuous matrix product states techniques, we explore the ground state phase diagram of this model by examining the superfluid and density correlation functions. At weak coupling superfluidity governs the ground state, in a similar way as in the Lieb-Liniger model. However, at strong coupling quasi-crystal and super-Tonks-Girardeau regimes are also found, which are not present in the original Lieb-Liniger case. Therefore the presence of the exponentially-decaying potential leads to a superfluid/super-Tonks-Girardeau/quasi-crystal crossover, when tuning the coupling strength from weak to strong interactions. This corresponds to a Luttinger liquid parameter in the range $K \in (0, \infty)$; in contrast with the Lieb-Liniger model, where $K \in [1, \infty)$, and the screened long-range potential, where $K \in (0, 1]$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonequilibrium, spatio-temporal formation of the Kondo screening-cloud on a lattice

We study the nonequilibrium formation of a spin screening cloud that accompanies the quenching of a local magnetic moment immersed in a Fermi sea at zero temperature. Based on high precision density matrix renormalization group results for the interacting single impurity Anderson model we discuss the real time evolution after a quantum quench in the impurity-reservoir hybridization using time evolving block decimation. We report emergent length and time scales in the spatiotemporal structure of non-local correlation functions in the spin- and the charge density channel. At equilibrium, our data for the correlation functions and the extracted length-scales show good agreement with existing results, as do local time dependent observables at the impurity. In the time-dependent data, we identify a major signal which defines a "light cone" moving at the Fermi velocity and a ferromagnetic component in its wake. Inside the light cone we find that the structure of the nonequilibrium correlation functions emerges on two time scales. Initially, the qualitative structure of the correlation functions develops rapidly at the lattice Fermi velocity. Subsequently the spin correlations converge to the equilibrium results on a much larger time scale. This process sets a dynamic energy scale, which we identify to be proportional to the Kondo temperature. Outside the light cone we observe two different power law decays of the correlation functions in space, with time and interaction strength independent exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Chebyshev expansion for Impurity Models using Matrix Product States

We improve a recently developed expansion technique for calculating real frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost and, in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than lowering it. We validate the method for the exactly solvable resonating level model and the single impurity Anderson model. It is capable of resolving sharp Kondo resonances, as well as peaks within the Hubbard bands when employed as an impurity solver for dynamical mean-field theory (DMFT). Our method works at zero temperature and allows for arbitrary discretization of the bath spectrum. It achieves similar precision as the dynamical density matrix renormalization group (DDMRG), at lower cost. We also propose an alternative expansion, of 1-exp(-tau H) instead of the usual H, which opens the possibility of using established methods for the time evolution of matrix product states to calculate spectral functions directly.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum Flutter: Signatures and Robustness

We investigate the motion of an impurity particle injected with finite velocity into an interacting one-dimensional quantum gas. Using large-scale numerical simulations based on matrix product states, we observe and quantitatively analyze long-lived oscillations of the impurity momentum around a non-zero saturation value, called quantum flutter. We show that the quantum flutter frequency is equal to the energy difference between two branches of collective excitations of the model. We propose an explanation of the finite saturation momentum of the impurity based on the properties of the edge of the excitation spectrum. Our results indicate that quantum flutter exists away from integrability, and provide parameter regions in which it could be observed in experiments with ultracold atoms using currently available technology.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Steady-state and quench dependent relaxation of a quantum dot coupled to one-dimensional leads

We study the time evolution and steady state of the charge current in a single-impurity Anderson model, using matrix product states techniques. A nonequilibrium situation is imposed by applying a bias voltage across one-dimensional tight-binding leads. Focusing on particle-hole symmetry, we extract current-voltage characteristics from universal low-bias up to high-bias regimes, where band effects start to play a dominant role. We discuss three quenches, which after strongly quench-dependent transients yield the same steady-state current. Among these quenches we identify those favorable for extracting steady-state observables. The period of short-time oscillations is shown to compare well to real-time renormalization group results for a simpler model of spinless fermions. We find indications that many-body effects play an important role at high-bias voltage and finite bandwidth of the metallic leads. The growth of entanglement entropy after a certain time scale (proportional to the inverse of Delta) is the major limiting factor for calculating the time evolution. We show that the magnitude of the steady-state current positively correlates with entanglement entropy. The role of high-energy states for the steady-state current is explored by considering a damping term in the time evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Clustered Wigner crystal phases of cold polar molecules in arrays of one-dimensional tubes

We analyze theoretically polar molecules confined in planar arrays of one dimensional tubes. In the classical limit, if the number of tubes is finite, new types of "clustered Wigner crystals" with increasingly many molecules per unit cell can be stabilized by tuning the in-plane angle between the dipolar moments and the tube direction. Quantum mechanically, these phases melt into distinct "clustered Luttinger liquids." We calculate the phase diagram of the system and study the quantum melting of the clustered phases. We find that the requirements for exploring these phases are reachable in current experiments and discuss possible experimental signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Steady-state spectra, current and stability diagram of a quantum dot: a non-equilibrium Variational Cluster Approach

We calculate steady-state properties of a strongly correlated quantum dot under voltage bias by means of non-equilibrium Cluster Perturbation Theory and the non-equilibrium Variational Cluster Approach, respectively. Results for the steady-state current are benchmarked against data from accurate Matrix Product State based time evolution. We show that for low to medium interaction strength, non-equilibrium Cluster Perturbation Theory already yields good results, while for higher interaction strength the self-consistent feedback of the non-equilibrium Variational Cluster Approach significantly enhances the accuracy. We report the current-voltage characteristics for different interaction strengths. Furthermore we investigate the non-equilibrium local density of states of the quantum dot and illustrate that within the variational approach a linear splitting and broadening of the Kondo resonance is predicted which depends on interaction strength. Calculations with applied gate voltage, away from particle hole symmetry, reveal that the maximum current is reached at the crossover from the Kondo regime to the doubly-occupied or empty quantum dot. Obtained stability diagrams compare very well to recent experimental data [Phys. Rev. B, 84, 245316 (2011)].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Observing complex bound states in the spin-1/2 Heisenberg XXZ chain using local quantum quenches

We consider the non-equilibrium evolution in the spin-1/2 XXZ Heisenberg chain for fixed magnetization after a local quantum quench. This model is equivalent to interacting spinless fermions. Initially an infinite magnetic field is applied to n consecutive sites and the ground state is calculated. At time t=0 the field is switched off and the time evolution of observables such as the z-component of spin is computed using the Time Evolving Block Decimation (TEBD) algorithm. We find that the observables exhibit strong signatures of linearly propagating spinon and bound state excitations. These persist even when integrability-breaking perturbations are included. Since bound states ("strings") are notoriously difficult to observe using conventional probes such as inelastic neutron scattering, we conclude that local quantum quenches are an ideal setting for studying their properties. We comment on implications of our results for cold atom experiments.

0 citations0 reviewsNo code linkedOpen access