Search

We present an open-source tensor network Python library for quantum many-body simulations. At its core is an abelian-symmetric tensor, implemented as a sparse block structure managed by logical layer on top of dense multi-dimensional array backend. This serves as the basis for higher-level tensor networks algorithms, operating on matrix product states and projected entangled pair states, implemented here. Using appropriate backend, such as PyTorch, gives direct access to automatic differentiation (AD) for cost-function gradient calculations and execution on GPUs or other supported accelerators. We show the library performance in simulations with infinite projected entangled-pair states, such as finding the ground states with AD, or simulating thermal states of the Hubbard model via imaginary time evolution. We quantify sources of performance gains in those challenging examples allowed by utilizing symmetries., Comment: 10 pages, 6 figures, source code at https://github.com/yastn/yastn

The Minimally Entangled Typical Thermal States (METTS) are an ensemble of pure states, equivalent to the Gibbs thermal state, that can be efficiently represented by tensor networks. In this article, we use the Projected Entangled Pair States (PEPS) ansatz as to represent METTS on a two-dimensional (2D) lattice. While Matrix Product States (MPS) are less efficient for 2D systems due to their complexity growing exponentially with the lattice size, PEPS provide a more tractable approach. To substantiate the prowess of PEPS in modeling METTS (dubbed as PEPS-METTS), we benchmark it against the purification method for the 2D quantum Ising model at its critical temperature. Our analysis reveals that PEPS-METTS achieves accurate long-range correlations with significantly lower bond dimensions. We further corroborate this finding in the 2D Fermi Hubbard model at half-filling. At a technical level, we introduce an efficient \textit{zipper} method to obtain PEPS boundary matrix product states needed to compute expectation values. The imaginary time evolution is performed with the neighbourhood tensor update., Comment: 14 pages, 14 figures

The Hubbard model is a longstanding problem in the theory of strongly correlated electrons and a very active one in the experiments with ultracold fermionic atoms. Motivated by current and prospective quantum simulations, we apply a two-dimensional tensor network, an infinite projected entangled pair state, evolved in imaginary time by the neighborhood tensor update algorithm working directly in the thermodynamic limit. With U(1)xU(1) symmetry and the bond dimensions up to 29, we generate thermal states down to the temperature of 0.17 times the hopping rate. We obtain results for spin and charge correlators, unaffected by boundary effects. The spin correlators, measurable in prospective ultracold atoms experiments attempting to approach the thermodynamic limit, provide evidence of disruption of the antiferromagnetic background with mobile holes in a slightly doped Hubbard model. The charge correlators reveal the presence of hole-doublon pairs near half filling and signatures of hole-hole repulsion on doping. We also obtain specific heat in the slightly doped regime., Comment: 18 pages; minor changes; accepted in PRB

Metastability is a quintessential feature of first order quantum phase transitions, which is lost either by dynamical instability or by nucleating bubbles of a true vacuum through quantum tunneling. By considering a drive across the first order quantum phase transition in the quantum Ising chain in the presence of both transverse and longitudinal fields, we reveal multiple regions in the parameter space where the initial metastable state loses its metastability in successive stages. The mechanism responsible is found to be semi-degenerate resonant tunnelings to states with specific bubble sizes. We show that such dynamics of quantized bubble nucleations can be understood in terms of Landau-Zener transitions, which provide quantitative predictions of nucleation probabilities for different bubble sizes., Comment: 7+4 pages, 5+3 figures. Close to published version

Driving a homogeneous system across a quantum phase transition in a quench-time $\tau_Q$ generates excitations on wavelengths longer than the Kibble-Zurek (KZ) length $\hat\xi\propto\tau_Q^{\nu/(1+z\nu)}$ within the KZ time window $\hat t\propto\tau_Q^{z\nu/(1+z\nu)}$, where $z$ and $\nu$ are the critical exponents. Quenches designed with local time-dependent inhomogeneity can introduce a gap in the spectrum. For a variety of setups with short-range interactions, they have been shown to suppress excitations if the spatial velocity of the inhomogenous front is below the characteristic KZ velocity $\hat v \propto \hat\xi/\hat t$. Ising-like models with long-range interactions can have no sonic horizon, spreading information instantaneously across the system. Usually, this should imply that inhomogenous transitions will render the dynamics adiabatic regardless of the front velocity. However, we show that we get an adiabatic transition with no defects only when the inhomogeneous front moves slower than a characteristic crossover velocity $\tilde v \propto \theta^{(z-1)\nu/(1+\nu)}$, where $\theta$ is the slope of the inhomogeneous front at the critical point. The existence of this crossover velocity and adiabaticity of the model is a consequence of the energy gap in the quasiparticle spectrum that is opened by the inhomogeneity. This effect can be employed for efficient adiabatic quantum state preparation in systems with long-range interactions., Comment: 12 pages, 10 figures (closed to published version)

A system gradually driven through a symmetry-breaking phase transition is subject to the Kibble-Zurek mechanism (KZM). As a consequence of the critical slowing down, its state cannot follow local equilibrium, and its evolution becomes non-adiabatic near the critical point. In the simplest approximation, that stage can be regarded as "impulse" where the state of the system remains unchanged. It leads to the correct KZM scaling laws. However, such "freeze-out" might suggest that the coherence length of the nascent order parameter remains unchanged as the critical region is traversed. By contrast, the original causality-based discussion emphasized the role of the {\it sonic horizon}: domains of the broken symmetry phase can expand with a velocity limited by the speed of the relevant sound. This effect was demonstrated in the quantum Ising chain where the dynamical exponent $z=1$ and quasiparticles excited by the transition have a fixed speed of sound. To elucidate the role of the sonic horizon, in this paper we study two systems with $z>1$ where the speed of sound is no longer fixed, and the fastest excited quasiparticles set the size of the sonic horizon. Their effective speed decays with the increasing transition time. In the extreme case, the dynamical exponent $z$ can diverge such as in the Griffiths region of the random Ising chain where localization of excited quasiparticles freezes the growth of the correlation range when the critical region is traversed. Of particular interest is an example with $z<1$ -- the long-range extended Ising chain, where there is no upper limit to the velocity of excited quasiparticles with small momenta. Initially, the power-law tail of the correlation function grows adiabatically, but in the non-adiabatic stage it lags behind the adiabatic evolution -- in accord with a generalized Lieb-Robinson bound., Comment: minor corrections, accepted in PRB

The Aubry-Andr\'e 1D lattice model describes a particle hopping in a pseudo-random potential. Depending on its strength $\lambda$, all eigenstates are either localized ($\lambda>1$) or delocalized ($\lambda<1$). Near the transition, the localization length diverges like $\xi\sim(\lambda-1)^{-\nu}$ with $\nu=1$. We show that when the particle is initially prepared in a localized ground state and the potential strength is slowly ramped down across the transition, then -- in analogy with the Kibble-Zurek mechanism -- it enters the delocalized phase having finite localization length $\hat\xi\sim\tau_Q^{\nu/(1+z\nu)}$. Here $\tau_Q$ is ramp/quench time and $z$ is a dynamical exponent. At $\lambda=1$ we determine $z\simeq2.37$ from the power law scaling of energy gap with lattice size $L$. Even though for infinite $L$ the model is gapless, we show that the gap relevant for excitation during the ramp remains finite. Close to the critical point it scales like $\xi^{-z}$ with the value of $z$ determined by the finite size scaling. It is the gap between the ground state and the lowest of those excited states that overlap with the ground state enough to be accessible for excitation. We propose an experiment with a non-interacting BEC to test our prediction. Our hypothesis is further supported by considering a generalized version of Aubry-Andr\'{e} model possessing an energy-dependent mobility edge., Comment: 7 pages, 9 figures; close to the published version; new section on generalized Aubry-Andr\'e model

A short quasi-monochromatic wave packet incident on a semi-infinite disordered medium gives rise to a reflected wave. The intensity of the latter decays as a power law $1/t^{\alpha}$ in the long-time limit. Using the one-dimensional Aubry-Andr\'{e} model, we show that in the vicinity of the critical point of Anderson localization transition, the decay slows down and the power-law exponent $\alpha$ becomes smaller than both $\alpha = 2$ found in the Anderson localization regime and $\alpha = 3/2$ expected for a one-dimensional random walk of classical particles., Comment: 9 pages, 6 figures. Revised text

I have first discussed how averaging theory can be an effective tool in solving weakly non-linear oscillators. Then I have applied this technique for a Van der Pol oscillator and extended the stability criterion of a Van der Pol oscillator for any integer n(odd or even)., Comment: results are not as generic as claimed. moreover it is not useful to the community