Publications

Order by: publication date number of citations

• Nevanlinna.jl: A Julia implementation of Nevanlinna analytic continuation

  Kosuke Nogaki, Jiani Fei, Emanuel Gull, Hiroshi Shinaoka

  SciPost Phys. Codebases 19 (2023) · published 13 November 2023 | Toggle abstract · pdf

  We introduce a Julia implementation of the recently proposed Nevanlinna analytic continuation method. The method is based on Nevanlinna interpolants and inherently preserves the causality of a response function due to its construction. For theoretical calculations without statistical noise, this continuation method is a powerful tool to extract real-frequency information from numerical input data on the Matsubara axis. This method has been applied to first-principles calculations of correlated materials. This paper presents its efficient and full-featured open-source implementation of the method including the Hamburger moment problem and smoothing.

  2 citations

• A quantum register using collective excitations in a Bose-Einstein condensate

  Elisha Haber, Zekai Chen, Nicholas P. Bigelow

  SciPost Phys. 15, 188 (2023) · published 8 November 2023 | Toggle abstract · pdf

  A qubit made up of an ensemble of atoms is attractive due to its resistance to atom losses. In this work, we consider an experimentally feasible protocol to coherently load a spin-dependent optical lattice from a spatially overlapping Bose-Einstein condensate. Identifying each lattice site as a qubit, with an empty or filled site as the qubit basis, we discuss how high-fidelity single-qubit operations, two-qubit gates between arbitrary pairs of qubits, and nondestructive measurements could be performed. In this setup, the effect of atom losses has been mitigated, the atoms never need to be removed from the ground state manifold, and separate storage and computational bases for the qubits are not required, all of which can be significant sources of decoherence in many other types of atomic qubits.

• Out of equilibrium many-body expansion dynamics of strongly interacting bosons

  Rhombik Roy, Barnali Chakrabarti, Arnaldo Gammal

  SciPost Phys. Core 6, 073 (2023) · published 3 November 2023 | Toggle abstract · pdf

  We solve the Schrödinger equation from first principles to investigate the many-body effects in the expansion dynamics of one-dimensional repulsively interacting bosons released from a harmonic trap. We utilize the multiconfigurational time-dependent Hartree method for bosons (MCTDHB) to solve the many-body Schrödinger equation at high level of accuracy. The MCTDHB basis sets are explicitly time-dependent and optimised by variational principle. We probe the expansion dynamics by three key measures; time evolution of one-, two- and three-body densities. We observe when the mean-field theory results to unimodal expansion, the many-body calculation exhibits trimodal expansion dynamics. The many-body features how the initially fragmented bosons independently spreads out with time whereas the mean-field pictures the expansion of the whole cloud. We also present the three different time scale of dynamics of the inner core, outer core and the cloud as a whole. We analyze the key role played by the dynamical fragmentation during expansion. A Strong evidence of the many-body effects is presented in the dynamics of two- and three-body densities which exhibit correlation hole and pronounced delocalization effect.

• Interference dynamics of matter-waves of SU($N$) fermions

  Wayne Jordan Chetcuti, Andreas Osterloh, Luigi Amico, Juan Polo

  SciPost Phys. 15, 181 (2023) · published 30 October 2023 | Toggle abstract · pdf

  Interacting N-component fermions spatially confined in ring-shaped potentials display specific coherence properties. The orbital angular momentum per particle of such systems can be quantized to fractional values specifically depending on the particle-particle interaction. Here we demonstrate how to monitor the state of the system through homodyne (momentum distribution) and self-heterodyne system's expansion. For homodyne protocols, the momentum distribution is affected by the particle statistics in two distinctive ways. The first effect is a purely statistical one: at zero interactions, the characteristic hole in the momentum distribution around the momentum k=0 opens up once half of the SU(N) Fermi sphere is displaced. The second effect originates from the interaction: The fractionalization in the interacting system manifests itself by an additional "delay" in the flux for the occurrence of the hole, that now becomes a characteristic minimum at k=0. We demonstrate that the angular momentum fractional quantization is reflected in the self-heterodyne interference as specific dislocations in interferograms. Our analysis demonstrates how the study of the interference fringes grants us access to both number of particles and number of components of SU(N) fermions.

  1 citation

• Resolving phase transitions with discontinuous Galerkin methods

  Eduardo Grossi, Nicolas Wink

  SciPost Phys. Core 6, 071 (2023) · published 25 October 2023 | Toggle abstract · pdf

  We demonstrate the applicability and advantages of Discontinuous Galerkin (DG) schemes in the context of the Functional Renormalization Group (fRG). We investigate the $O(N)$-model in the large $N$ limit. It is shown that the flow equation for the effective potential can be cast into a conservative form. We discuss results for the Riemann problem, as well as initial conditions leading to a first and second order phase transition. In particular, we unravel the mechanism underlying first order phase transitions, based on the formation of a shock in the derivative of the effective potential.

  1 citation

• A time-dependent momentum-resolved scattering approach to core-level spectroscopies

  Krissia Zawadzki, Alberto Nocera, Adrian E. Feiguin

  SciPost Phys. 15, 166 (2023) · published 17 October 2023 | Toggle abstract · pdf

  While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

• Interfacially adsorbed bubbles determine the shape of droplets

  Alessio Squarcini, Antonio Tinti

  SciPost Phys. 15, 164 (2023) · published 16 October 2023 | Toggle abstract · pdf

  The characterization of density correlations in the presence of strongly fluctuating interfaces has always been considered a difficult problem in statistical mechanics. Here we study - by using recently developed exact field-theoretical techniques - density correlations for an interface with endpoints on a wall forming a droplet in 2D. Our framework applies to interfaces entropically repelled by a hard wall as well as to wetting transitions. In the former case bubbles adsorbed on the interface are taken into account by the theory which yields a systematic treatment of finite-size corrections to one- and two-point functions and show how these are related to Brownian excursions. Our analytical predictions are confirmed by Monte Carlo simulations without free parameters. We also determine one- and two-point functions at wetting by using integrable boundary field theory. We show that correlations are long ranged for entropic repulsion and at wetting. For both regimes we investigate correlations in momentum space by generalizing the notion of interface structure factor to semi-confined systems. Distinctive signatures of the two regimes manifest in the structure factor through a term that we identify on top of the capillary-wave one.

• Skyrmion crystals in the triangular Kondo lattice model

  Zhentao Wang, Cristian D. Batista

  SciPost Phys. 15, 161 (2023) · published 13 October 2023 | Toggle abstract · pdf

  We present a systematic study of the formation of skyrmion crystals in a triangular Kondo Lattice model for generic electron filling fractions. Our results indicate that the four-sublattice chiral antiferromagnetic ordering that was reported more than one decade ago can be understood as the dense limit of a sequence of skyrmion crystals whose lattice parameter is dictated by the Fermi wave-vector. This observation has important implications for the ongoing search of skyrmion crystals in metallic materials with localized magnetic moments.

  4 citations

• Tensor network variational optimizations for real-time dynamics: Application to the time-evolution of spin liquids

  Ravi Teja Ponnaganti, Matthieu Mambrini, Didier Poilblanc

  SciPost Phys. 15, 158 (2023) · published 12 October 2023 | Toggle abstract · pdf

  Within the Projected Entangled Pair State (PEPS) tensor network formalism, a simple update (SU) method has been used to investigate the time evolution of a two-dimensional U(1) critical spin-1/2 spin liquid under Hamiltonian quench [Phys. Rev. B 106, 195132 (2022)]. Here we introduce two different variational frameworks to describe the time dynamics of SU(2)-symmetric translationally-invariant PEPS, aiming to improve the accuracy. In one approach, after using a Trotter-Suzuki decomposition of the time evolution operator in term of two-site elementary gates, one considers a single bond embedded in an environment approximated by a Corner Transfer Matrix Renormalization Group (CTMRG). A variational update of the two tensors on the bond is performed under the application of the elementary gate and then, after symmetrization of the site tensors, the environment is updated. In the second approach, a cluster optimization is performed on a finite (periodic) cluster, maximizing the overlap of the exact time-evolved state with a symmetric finite-size PEPS ansatz. Observables are then computed on the infinite lattice contracting the infinite-PEPS (iPEPS) by CTMRG. We show that the variational schemes outperform the SU method and remain accurate over a significant time interval before hitting the entanglement barrier. Studying the spectrum of the transfer matrix, we find that the asymptotic correlations are very well preserved under time evolution, including the critical nature of the singlet correlations, as expected from the Lieb-Robinson (LR) bound theorem. Consistently, the system (asymptotic) boundary is found to be described by the same Conformal Field Theory of central charge $c = 1$ during time evolution. We also compute the time-evolution of the short distance spin-spin correlations and estimate the LR velocity.

• Monte Carlo matrix-product-state approach to the false vacuum decay in the monitored quantum Ising chain

  Jeffrey Allan Maki, Anna Berti, Iacopo Carusotto, Alberto Biella

  SciPost Phys. 15, 152 (2023) · published 11 October 2023 | Toggle abstract · pdf

  In this work we characterize the false vacuum decay in the ferromagnetic quantum Ising chain with a weak longitudinal field subject to continuous monitoring of the local magnetization. Initializing the system in a metastable state, the false vacuum, we study the competition between coherent dynamics, which tends to create resonant bubbles of the true vacuum, and measurements which induce heating and reduce the amount of quantum correlations. To this end we exploit a numerical approach based on the combination of matrix product states with stochastic quantum trajectories which allows for the simulation of the trajectory-resolved non-equilibrium dynamics of interacting many-body systems in the presence of continuous measurements. We show how the presence of measurements affects the false vacuum decay: At short times the departure from the local minimum is accelerated while at long times the system thermalizes to an infinite-temperature incoherent mixture. For large measurement rates the system enters a quantum Zeno regime. The false vacuum decay and the thermalization physics are characterized in terms of the magnetization, connected correlation function, and the trajectory-resolved entanglement entropy.

  1 citation