DOCSLIB.ORG

Arxiv:1906.04243V2 [Cond-Mat.Str-El]

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Arxiv:1906.04243V2 [Cond-Mat.Str-El] Applying the tempered Lefschetz thimble method to the Hubbard model away from half-filling 1, 1, 2, Masafumi Fukuma, ∗ Nobuyuki Matsumoto, † and Naoya Umeda ‡ 1Department of Physics, Kyoto University, Kyoto 606-8502, Japan 2PricewaterhouseCoopers Aarata LLC, Otemachi Park Building, 1-1-1 Otemachi, Chiyoda-ku, Tokyo 100-0004, Japan The tempered Lefschetz thimble method is a parallel-tempering algorithm towards solving the numerical sign problem. It uses the flow time of the gradient flow as a tempering parameter and is expected to tame both the sign and multimodal problems simultaneously. In this paper, we further develop the algorithm so that the expectation values can be estimated precisely with a criterion ensuring global equilibrium and the sufficiency of the sample size. To demonstrate that this algorithm works well, we apply it to the quantum Monte Carlo simulation of the Hubbard model away from half-filling on a two-dimensional lattice of small size, and show that the numerical results agree nicely with exact values. I. INTRODUCTION gradient flow (see also [15] for a similar idea). In this paper, we further develop the TLTM, propos- The sign problem is one of the major obstacles when ing an algorithm which allows the precise estimation of performing numerical calculations in various fields of expectation values with a criterion ensuring global equi- physics. Typical examples include finite density QCD [1], librium and the sufficiency of the sample size. The key quantum Monte Carlo (QMC) calculations of quantum is the use of the fact that the expectation values should statistical systems [2–4], and the numerical simulations be the same for all flow times. To demonstrate that this of real-time quantum field theories. algorithm works well, we apply it to the QMC simulation Among a variety of approaches, two algorithms have of the Hubbard model away from half-filling. taken attention as potential candidates to generically The application of Lefschetz thimble methods to the solve the sign problem for systems with complex action; Hubbard model has already been considered by several one is the complex Langevin method [5], and the other groups [20–22] (see also [23, 24] for recent study), and is a class of algorithms utilizing the Lefschetz thimbles the relevance of the contributions from multiple thimbles [6–15]. Although both the algorithms make use of com- has been reported. In this paper, we consider a two- plexification of variables and analytic continuation of in- dimensional periodic square lattice of size N = 2 2 s × tegrands, their methodologies are fairly different; the for- with the inverse temperature decomposed to Nτ = 5 mer algorithm attempts to replace the complex Boltz- pieces, and numerically evaluate the expectation values mann weight by a real positive weight defined in the of observables as functions of the chemical potential with whole complex space, while the latter deforms the in- other parameters fixed to some values. We show that the tegration region in the complex space so as to reduce the TLTM (the implementation of tempering combined with phase oscillation. At this stage, each algorithm has its the above algorithm for precise estimation) give results own advantage and disadvantage. The former is advan- that agree nicely with exact values, simultaneously re- tageous in that it is relatively fast with computational solving the sign and multimodal problems. cost O(N) (N: the degrees of freedom), but it suffers We comment that the extent of seriousness of the from the so-called wrong convergence problem [16–19]. sign problem in the QMC simulation of the Hubbard The latter is generally free from the wrong convergence model depends heavily on the choice of the Hubbard- problem if only a single thimble is relevant in evaluating Stratonovich variables. In this paper, in order to apply the expectation values of physical observables of inter- the Lefschetz thimble method, we exclusively consider a est. The disadvantage is its expensive numerical cost, arXiv:1906.04243v2 [cond-mat.str-el] 27 Nov 2019 3 Gaussian Hubbard-Stratonovich variable that leads to a which is O(N ) because of the need to calculate the Ja- complex action. There the sign problem is actually se- cobian determinant. When multiple thimbles are rele- vere as we will see below, and one needs to seriously con- vant, one needs to take care of the multimodality of the sider a dilemma between the sign and multimodal prob- distribution. The tempered Lefschetz thimble method lems, which can be solved by the TLTM as stated above. (TLTM) was thus proposed in [14] to tame both the sign However, the temporal size considered here is still small and multimodal problems simultaneously, where the sys- (Nτ = 5), and for such a high temperature regime one tem is tempered by the flow time of the antiholomorphic can resort to other methods than the Lefschetz thimble methods with a different type of Hubbard-Stratonovich variables (see discussions in section III). ∗ [email protected] This paper is organized as follows. In section II after † [email protected] briefly reviewing the TLTM [14], we give a new algorithm ‡ [email protected] which allows the precise estimation of expectation values 2 with a criterion ensuring global equilibrium and the suf- with H(z) (∂i∂j S(z)). Under the flow (2), the ac- ficiency of the sample size. This algorithm is applied to ≡ 2 tion changes as (d/dt)S(zt(x)) = ∂iS(zt(x)) 0, and the Hubbard model in section III, and we discuss about ≥ thus Re S(zt(x)) increases except at the critical points the obtained numerical results. We there also make a z (∂iS(z ) = 0), while Im S(zt(x)) is kept constant. In comment on the sign averages obtained by other meth- particular,∗ ∗ in the limit t , the deformed region will ods. Section IV is devoted to conclusion and outlook. approach a union of N-dimensional→ ∞ submanifolds (Lef- Five appendices are given for more detailed discussions schetz thimbles) on each of which Im S(z) is constant, on various topics. and thus the sign problem is expected to disappear there (except for a possible residual sign problem arising from the phase of the complex measure dz and a possible global II. TEMPERED LEFSCHETZ THIMBLE sign problem caused by phase cancellations among differ- METHOD ent thimbles). However, in the Monte Carlo calculation one cannot take the t limit na¨ıvely, because the Let x = (xi) RN be a real N-dimensional dynamical potential barriers between→ ∞ different thimbles become in- ∈ variable with action S(x) which may take complex values. finitely high so that the whole configuration space cannot Our main concern is to estimate the expectation values be explored sufficiently. This multimodality of distribu- tion makes the Monte Carlo calculation impractical, es- S(x) N dx e− (x) R O pecially when contributions from more than one thimble (x) S S(x) . (1) hO i ≡ R RN dx e− are relevant to estimating expectation values. A key pro- R posal in [12] is to use a finite value of flow time that S(z) S(z) We assume that e− and e− (z) are entire func- is large enough to avoid the sign problem but simulta- tions over CN when x is complexifiedO to z = (zi) CN . neously is not too large so that the exploration in the Then, due to Cauchy’s theorem for higher dimensions,∈ configuration space is still possible. However, it is a dif- the right-hand side does not change under continuous de- ficult task to find such value of flow time in a systematic formations of the integration region as long as the bound- way, as we will discuss at the end of section III and in ary at infinity is kept fixed so that the integrals converge. Appendix E. The sign problem will get reduced if Im S(z) is almost The TLTM [14] is a tempering algorithm that uses the constant on the new integration region. flow time as a tempering parameter. There, the global In [11–15] such a deformation x zt(x) (t 0) is relaxation of the multimodal distribution is prompted by made according to the antiholomorphic→ gradient≥ flow: enabling configurations around different modes to easily communicate through transitions in ensembles at smaller i i i z˙t = [∂iS(zt)]∗, zt=0 = x . (2) flow times. Among other possible tempering algorithms, the parallel tempering algorithm [25, 26] (also known as Equation (1) can then be rewritten as the replica exchange MCMC method; see [27] for a re- view) is adopted in the TLTM [14] because it does not S(z) Σ dz e− (z) need to introduce the probability weight factors of ensem- (x) = t O (Σ z (RN )), (3) hO iS R dz e S(z) t ≡ t bles at various flow times and because most of relevant Σt − R steps can be done in parallel processes. which can be further rewritten as a ratio of reweighted In the TLTM (see Appendix A for the summary of N integrals over R by using the Jacobian matrix Jt(x) the algorithm), we first fix the maximum flow time T i j ≡ ∂zt(x)/∂x [11]: which should be sufficiently large such that the sign prob- lem is reduced there. A possible criterion is that the S(z (x)) iθT (x) t sign average e eff is O(1) in the absence of tem- RN dx detJt(x) e− (zt(x)) ST (x) S = O |h i | S(zt(x)) pering. This process can be carried out by a test run hO i R RN dx detJt(x) e− with small statistics.
Load more

Recommended publications
  • Arxiv:2007.10198V2 [Math.NA] 5 Nov 2020
    ON THE VALIDITY OF COMPLEX LANGEVIN METHOD FOR PATH INTEGRAL COMPUTATIONS∗ ZHENNING CAIy , XIAOYU DONGz , AND YANG KUANGx Abstract. The complex Langevin (CL) method is a classical numerical strategy to alleviate the numerical sign problem in the computation of lattice field theories. Mathematically, it is a simple numerical tool to compute a wide class of high- dimensional and oscillatory integrals. However, it is often observed that the CL method converges but the limiting result is incorrect. The literature has several unclear or even conflicting statements, making the method look mysterious. By an in-depth analysis of a model problem, we reveal the mechanism of how the CL result turns biased as the parameter changes, and it is demonstrated that such a transition is difficult to capture. Our analysis also shows that the method works for any observables only if the probability density function generated by the CL process is localized. To generalize such observations to lattice field theories, we formulate the CL method on general groups using rigorous mathematical languages for the first time, and we demonstrate that such localized probability density function does not exist in the simulation of lattice field theories for general compact groups, which explains the unstable behavior of the CL method. Fortunately, we also find that the gauge cooling technique creates additional velocity that helps confine the samples, so that we can still see localized probability density functions in certain cases, as significantly broadens the application of the CL method. The limitations of gauge cooling are also discussed. In particular, we prove that gauge cooling has no effect for Abelian groups, and we provide an example showing that biased results still exist when gauge cooling is insufficient to confine the probability density function.
    [Show full text]
  • Quantum Mechanics Quantum Chromodynamics (QCD)
    Quantum Mechanics_quantum chromodynamics (QCD) In theoretical physics, quantum chromodynamics (QCD) is a theory ofstrong interactions, a fundamental forcedescribing the interactions between quarksand gluons which make up hadrons such as the proton, neutron and pion. QCD is a type of Quantum field theory called a non- abelian gauge theory with symmetry group SU(3). The QCD analog of electric charge is a property called 'color'. Gluons are the force carrier of the theory, like photons are for the electromagnetic force in quantum electrodynamics. The theory is an important part of the Standard Model of Particle physics. A huge body of experimental evidence for QCD has been gathered over the years. QCD enjoys two peculiar properties: Confinement, which means that the force between quarks does not diminish as they are separated. Because of this, when you do split the quark the energy is enough to create another quark thus creating another quark pair; they are forever bound into hadrons such as theproton and the neutron or the pion and kaon. Although analytically unproven, confinement is widely believed to be true because it explains the consistent failure of free quark searches, and it is easy to demonstrate in lattice QCD. Asymptotic freedom, which means that in very high-energy reactions, quarks and gluons interact very weakly creating a quark–gluon plasma. This prediction of QCD was first discovered in the early 1970s by David Politzer and by Frank Wilczek and David Gross. For this work they were awarded the 2004 Nobel Prize in Physics. There is no known phase-transition line separating these two properties; confinement is dominant in low-energy scales but, as energy increases, asymptotic freedom becomes dominant.
    [Show full text]
  • Field Theory Research Team
    Chapter 5 Field Theory Research Team 5.1 Members Yoshinobu Kuramashi (Team Leader) Yoshifumi Nakamura (Research Scientist) Hiroya Suno (Research Scientist, Joint Position with the Nishina Center for Accelerator-based Research) Eigo Shintani (Research Scientist) Yuya Shimizu (Postdoctoral Researcher) Yusuke Yoshimura (Postdoctoral Researcher) Ken-Ichi Ishikawa (Visiting Scientist, Hiroshima University) Takeshi Yamazaki (Visiting Scientist, University of Tsukuba) Shinji Takeda (Visiting Scientist, Kanazawa University) 5.2 Research Activities Our research field is physics of elementary particles and nuclei, which tries to answer questions in history of mankind: What is the smallest component of matter and what is the most fundamental interactions? This research subject is related to the early universe and the nucleosynthesis through Big Bang cosmology. Another important aspect is quantum properties, which play an essential role in the world of elementary particles and nuclei as well as in the material physics at the atomic or molecular level. We investigate nonperturbative properties of elementary particles and nuclei through numerical simulations with the use of lattice QCD (Quantum ChromoDynamics). The research is performed in collaboration with applied mathematicians, who are experts in developing and improving algorithms, and computer scientists responsible for research and development of software and hardware systems. Lattice QCD is one of the most advanced case in quantum sciences: Interactions between quarks, which are elementary particles known to date, are described by QCD formulated with the quantum field theory. We currently focus on two research subjects: (1) QCD at finite temperature and finite density. We try to understand the early universe and the inside of neutron star by investigating the phase structure and the equation of state.
    [Show full text]
  • New Insights Into the Sign Problem of Qcd and Heavy Tetraquarks
    ABSTRACT Title of dissertation: NEW INSIGHTS INTO THE SIGN PROBLEM OF QCD AND HEAVY TETRAQUARKS Yiming Cai Doctor of Philosophy, 2020 Dissertation directed by: Thomas Cohen Department of Physics Quantum Chromodynamics (QCD) describes strong interactions among the fundamental particles known as quarks and gluons. In principle, QCD can be used to explain complicated phenomena in the strong sector. However, at the energy scale of hadronic physics, the strong coupling constant is so large that the traditional per- turbative method is not applicable. There are two powerful alternative approaches utilized in the non-pertubative regime: (1) lattice QCD, which discretizes space-time and utilizes Monte Carlo computer simulations, and (2) finding new systematic ex- pansion regimes to obtain physical insights in certain limits. In this dissertation, three problems are studied in the context of these two approaches. In Chapter 2, a notorious numeric problem in lattice QCD known as the sign problem is explored. A subtle phenomenon caused by the interplay between the sign problem and the infinite volume limit is discussed and explained using the saddle point approximation. This work provides insight into the sign problem and the physics of the QCD θ-vacuum. Chapter 3 and Chapter 4 focus on tetraquarks, which are unconventional hadrons containing four valence quarks. Despite numerous tetraquark candidates seen in experiments, there is no unified and well-accepted theoretical descriptions of the tetraquark state yet. This dissertation examines the existence of tetraquarks in the heavy quark mass limit. A powerful systematic expansion regime can be built when the heavy quark mass is extremely large.
    [Show full text]
  • Bulk Nuclear Properties from Dynamical Description of Heavy-Ion Collisions
    BULK NUCLEAR PROPERTIES FROM DYNAMICAL DESCRIPTION OF HEAVY-ION COLLISIONS By Jun Hong A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of Physics-Doctor of Philosophy 2016 ABSTRACT BULK NUCLEAR PROPERTIES FROM DYNAMICAL DESCRIPTION OF HEAVY-ION COLLISIONS By Jun Hong Mapping out the equation of state (EOS) of nuclear matter is a long standing problem in nuclear physics. Both experimentalists and theoretical physicists spare no effort in improving understand- ing of the EOS. In this thesis, we examine observables sensitive to the EOS within the pBUU transport model based on the Boltzmann equation. By comparing theoretical predictions with ex- perimental data, we arrive at new constraints for the EOS. Further we propose novel promising observables for analysis of future experimental data. One set of observables that we examine within the pBUU model are pion yields. First, we find that net pion yields in central heavy-ion collisions (HIC) are strongly sensitive to the momentum dependence of the isoscalar nuclear mean field. We reexamine the momentum dependence that is assumed in the Boltzmann equation model for the collisions and optimize that dependence to describe the FOPI measurements of pion yields from the Au+Au collisions at different beam ener- gies. Alas such optimized dependence yields a somewhat weaker baryonic elliptic flow than seen in measurements. Subsequently, we use the same pBUU model to generate predictions for baryonic elliptic flow observable in HIC, while varying the incompressibility of nuclear matter. In parallel, we test the sensitivity of pion multiplicity to the density dependence of EOS, and in particular to incompress- ibility, and optimize that dependence to describe both the elliptic flow and pion yields.
    [Show full text]
  • Arxiv:2002.01217V1 [Hep-Ph] 4 Feb 2020 Diagram Is Reached Yet
    Nuclear Physics A Nuclear Physics A 00 (2020) 1–4 www.elsevier.com/locate/procedia XXVIIIth International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2019) QCD equation of state at vanishing and high baryon density: Chiral Mean Field model Anton Motornenkoa,b, Jan Steinheimerb, Volodymyr Vovchenkoa,b, Stefan Schrammyb, Horst Stoeckera,b,c aInstitut f¨urTheoretische Physik, Goethe Universit¨atFrankfurt, D-60438 Frankfurt am Main, Germany bFrankfurt Institute for Advanced Studies, Giersch Science Center, D-60438 Frankfurt am Main, Germany cGSI Helmholtzzentrum f¨urSchwerionenforschung GmbH, D-64291 Darmstadt, Germany Abstract The thermodynamic properties of high temperature and high density QCD-matter are studied using the Chiral SU(3)- flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The CMF model provides a proper descrip- tion of lattice QCD data, heavy-ions physics, and static neutron stars. The behavior of lines of constant pressure with increase of baryon density is discussed. The rapid change of pressure behavior at µB=T ≈ 3 suggests a strong contri- bution of baryons to thermodynamic properties at this region. The position of this region is very close to the radius of convergence for a Taylor expansion of the QCD pressure. The role of mesons and unstable hadrons in the hydrodynamic expansion of strongly interacting matter is also discussed. Keywords: QCD equation of state, parity-doublet model, QCD phase diagram 1. Introduction The well established theory of strong interactions – Quantum Chromodynamics (QCD) – suffers calcula- tion problems in the non-perutrbative regime: large coupling constant disfavors perturbative methods while the numerical sign problem prohibits lattice (LQCD) calculations at finite densities.
    [Show full text]
  • Fermi Surface Reconstruction Without Symmetry Breaking
    Fermi surface reconstruction without symmetry breaking Snir Gazit Racah Institute of Physics and The Fritz Haber Research Center for Molecular Dynamics, The Hebrew University of Jerusalem, Jerusalem 91904, Israel Fakher F. Assaad Institut f¨urTheoretische Physik und Astrophysik, Universit¨atW¨urzburg, Am Hubland, D-97074 W¨urzburg, Germany and W¨urzburg-Dresden Cluster of Excellence ct.qmat, Universit¨atW¨urzburg, Am Hubland, D-97074 W¨urzburg, Germany Subir Sachdev Department of Physics, Harvard University, Cambridge MA 02138, USA Abstract We present a sign-problem free quantum Monte Carlo study of a model that exhibits quantum phase transitions without symmetry breaking and associated changes in the size of the Fermi surface. The model is an Ising gauge theory on the square lattice coupled to an Ising matter field and spinful `orthogonal' fermions at half-filling, both carrying Ising gauge charges. In contrast to previous studies, our model hosts an electron-like, gauge-neutral fermion excitation providing access to Fermi liquid phases. One of the phases of the model is a previously studied orthogonal semi-metal, which has Z2 topological order, and Luttinger-volume violating Fermi points with gapless orthogonal fermion excitations. We elucidate the global phase diagram of the model: along with a conventional Fermi liquid phase with a large Luttinger-volume Fermi surface, we also find a ‘deconfined’ Fermi liquid in which the large Fermi surface co-exists with fractionalized excitations. We present results for the electron spectral function, showing its evolution from the orthogonal arXiv:1906.11250v3 [cond-mat.str-el] 26 Jan 2021 semi-metal with a spectral weight near momenta {±π=2; ±π=2g, to a large Fermi surface.
    [Show full text]
  • Van Der Waals Interactions in the Hadron Resonance Gas Model
    Van der Waals Interactions in the Hadron Resonance Gas Model by Aaron Boggs A dissertation submitted to the Department of Physics, College of Natural Sciences and Mathematics in partial fulfillment of the requirements for the degree of Bachelor of Science in Physics Chair of Committee: Claudia Ratti Committee Member: Rene Bellwied Committee Member: Ralph Metcalfe University of Houston August 2019 Copyright 2019, Aaron Boggs DEDICATION/EPIGRAPH I dedicate this thesis to TikTok for getting me through my many struggles and hardships while working on this thesis. iii ACKNOWLEDGMENTS Firstly, I would like to thank Dr. Claudia Ratti for supporting me, teaching me and giving me the opportunity to work with her for the past three years. I would also like to thank Paolo Parotto for spending many hours helping me with the work in this thesis and giving me much guidance. They have both taught me so much and I will continue to carry what I have learned from both of them for a long time to come. Secondly, I would like to thank my family and friends who encouraged me and gave me motiva- tion throughout my time at the University of Houston to continue following my dreams. Without their support, I would not have gotten through the past four years. iv ABSTRACT The Quark-Gluon Plasma (QGP) and its phase transition on the Quantum Chromodynamics (QCD) phase diagram have been at the forefront of high energy physics research for the past few decades. In order to study the QGP and its thermodynamic behavior, many experiments have been undertaken to recreate this state of matter at particle colliders like the Large Hadron Col- lider (LHC) at CERN and the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven National Laboratory.
    [Show full text]
  • Universal Four-Fermion Interaction in Lattice Effective Field Theory
    Mississippi State University Scholars Junction Theses and Dissertations Theses and Dissertations 1-1-2015 Universal Four-Fermion Interaction in Lattice Effective Field Theory Kristopher Katterjohn Follow this and additional works at: https://scholarsjunction.msstate.edu/td Recommended Citation Katterjohn, Kristopher, "Universal Four-Fermion Interaction in Lattice Effective Field Theory" (2015). Theses and Dissertations. 4849. https://scholarsjunction.msstate.edu/td/4849 This Graduate Thesis - Open Access is brought to you for free and open access by the Theses and Dissertations at Scholars Junction. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of Scholars Junction. For more information, please contact [email protected]. Universal four-fermion interaction in lattice effective feld theory By Kristopher J. Katterjohn A Thesis Submitted to the Faculty of Mississippi State University in Partial Fulfllment of the Requirements for the Degree of Master of Science in Physics in the Department of Physics and Astronomy Mississippi State, Mississippi August 2015 Copyright by Kristopher J. Katterjohn 2015 Universal four-fermion interaction in lattice effective feld theory By Kristopher J. Katterjohn Approved: Gautam Rupak Lan Tai Moong (Major Professor) Dipangkar Dutta (Committee Member) R. Torsten Clay (Committee Member) Hendrik F. Arnoldus (Graduate Coordinator) R. Gregory Dunaway Professor and Dean College of Arts & Sciences Name: Kristopher J. Katterjohn Date of Degree: August 14, 2015 Institution: Mississippi State University Major Field: Physics Major Professor: Dr. Gautam Rupak Lan Tai Moong Title of Study: Universal four-fermion interaction in lattice effective feld theory Pages of Study: 42 Candidate for Degree of Master of Science In this thesis we study non-relativistic, low-energy, s-wave scattering in a four-body spin-1/2 fermion system.
    [Show full text]
  • Arxiv:2008.11206V2 [Cond-Mat.Str-El] 30 Apr 2021 Topological Paramagnets, Protected by Z2 and Z2 Symmetries, 47, 49, 52–55]
    <latexit sha1_base64="9wYEiU3lOHr7PEMg5nq4gTnAkMc=">AAAB83icbVDLSsNAFL2pr1pfVZduBovgqiSloO6KblxWsK3YhDKZTtqhk0mYh1BCf8ONC0Xc+jPu/BsnbRbaemDgcM693DMnTDlT2nW/ndLa+sbmVnm7srO7t39QPTzqqsRIQjsk4Yl8CLGinAna0Uxz+pBKiuOQ0144ucn93hOViiXiXk9TGsR4JFjECNZW8v0Y63EYZo+zQWNQrbl1dw60SryC1KBAe1D98ocJMTEVmnCsVN9zUx1kWGpGOJ1VfKNoiskEj2jfUoFjqoJsnnmGzqwyRFEi7RMazdXfGxmOlZrGoZ3MM6plLxf/8/pGR5dBxkRqNBVkcSgyHOkE5QWgIZOUaD61BBPJbFZExlhiom1NFVuCt/zlVdJt1L1m/equWWtdF3WU4QRO4Rw8uIAW3EIbOkAghWd4hTfHOC/Ou/OxGC05xc4x/IHz+QP7v5Gr</latexit> <latexit sha1_base64="UNTpOsNPwaEFlt8H9lupFz9e49I=">AAAB9XicbVDLSsNAFL2pr1pfVZduBovgqiQ1oMuiG5cV7APbtEymk3boZBJmJkoJ/Q83LhRx67+482+ctFlo64GBwzn3cs8cP+ZMadv+tgpr6xubW8Xt0s7u3v5B+fCopaJEEtokEY9kx8eKciZoUzPNaSeWFIc+p21/cpP57UcqFYvEvZ7G1AvxSLCAEayN1O+FWI99P32YDWr9i0G5YlftOdAqcXJSgRyNQfmrN4xIElKhCcdKdR071l6KpWaE01mplygaYzLBI9o1VOCQKi+dp56hM6MMURBJ84RGc/X3RopDpaahbyazlGrZy8T/vG6igysvZSJONBVkcShIONIRyipAQyYp0XxqCCaSmayIjLHERJuiSqYEZ/nLq6RVqzpu1b1zK/XrvI4inMApnIMDl1CHW2hAEwhIeIZXeLOerBfr3fpYjBasfOcY/sD6/AEptZJL</latexit> <latexit sha1_base64="MmiC13il+j/rEQ0mZF3fuzKucTA=">AAAB7XicbVDLSgNBEOz1GeMr6tHLYBA8hV0J6DHoxWME84BkCb2TSTJmdmaZmRXCkn/w4kERr/6PN//GSbIHTSxoKKq66e6KEsGN9f1vb219Y3Nru7BT3N3bPzgsHR03jUo1ZQ2qhNLtCA0TXLKG5VawdqIZxpFgrWh8O/NbT0wbruSDnSQsjHEo+YBTtE5qdlEkI+yVyn7Fn4OskiAnZchR75W+un1F05hJSwUa0wn8xIYZasupYNNiNzUsQTrGIes4KjFmJszm107JuVP6ZKC0K2nJXP09kWFszCSOXGeMdmSWvZn4n9dJ7eA6zLhMUsskXSwapIJYRWavkz7XjFoxcQSp5u5WQkeokVoXUNGFECy/vEqal5WgWqneV8u1mzyOApzCGVxAAFdQgzuoQwMoPMIzvMKbp7wX7937WLSuefnMCfyB9/kDjcePHw==</latexit> <latexit sha1_base64="MmiC13il+j/rEQ0mZF3fuzKucTA=">AAAB7XicbVDLSgNBEOz1GeMr6tHLYBA8hV0J6DHoxWME84BkCb2TSTJmdmaZmRXCkn/w4kERr/6PN//GSbIHTSxoKKq66e6KEsGN9f1vb219Y3Nru7BT3N3bPzgsHR03jUo1ZQ2qhNLtCA0TXLKG5VawdqIZxpFgrWh8O/NbT0wbruSDnSQsjHEo+YBTtE5qdlEkI+yVyn7Fn4OskiAnZchR75W+un1F05hJSwUa0wn8xIYZasupYNNiNzUsQTrGIes4KjFmJszm107JuVP6ZKC0K2nJXP09kWFszCSOXGeMdmSWvZn4n9dJ7eA6zLhMUsskXSwapIJYRWavkz7XjFoxcQSp5u5WQkeokVoXUNGFECy/vEqal5WgWqneV8u1mzyOApzCGVxAAFdQgzuoQwMoPMIzvMKbp7wX7937WLSuefnMCfyB9/kDjcePHw==</latexit>
    [Show full text]
  • Quantum Monte Carlo Studies of Fermi Systems in Lattice Effective Field Theory
    QUANTUM MONTE CARLO STUDIES OF FERMI SYSTEMS IN LATTICE EFFECTIVE FIELD THEORY By Rongzheng He A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of Physics { Doctor of Philosophy 2019 ABSTRACT QUANTUM MONTE CARLO STUDIES OF FERMI SYSTEMS IN LATTICE EFFECTIVE FIELD THEORY By Rongzheng He With the rapid growth of the computational resources, quantum Monte Carlo (QMC) methodology has become a powerful tool for numerical simulations, especially for lattice effective field theory. Those QMC simulations have successfully described the physics of the few- and many-body systems. In this thesis, we investigate the Fermi systems with balanced and unbalanced populations of up and down spins using QMC with ab initio techniques. We also present a newly developed method called eigenvector continuation (EC) and its promising applications to some numerically unavoidable problems, like the sign problem in Monte Carlo simulations. For the eigenvector continuation method, we demonstrate that although Hamiltonian is usually represented as a matrix in a linear space with enormous dimensions, the eigenvector trajectory generated by a smoothly changed Hamitonian matrix is well approximated by a low-dimensional space. We use analytic continuation theory to prove this statement and propose an algorithm to implement our method. In the simulation with strong numerical sign oscillations, we first \learn" the subspace where the trajectory is approximately spanned by a finite number of accurately computable eigenvectors and then apply eigenvetor continuation to solve the physics system where there is a severe sign problem. Our results converge rapidly as we include more eigenvectors. The results show that for the same computational cost, the EC method reduces errors by an order of magnitude compared to the direct calculations in cases when we have a strong sign problem.
    [Show full text]
  • Arxiv:2006.07654V1 [Math.NA] 13 Jun 2020 C278 USA
    Numerical analysis for inchworm Monte Carlo method: Sign problem and error growth Zhenning Cai1, Jianfeng Lu2, Siyao Yang1 Abstract We consider the numerical analysis of the inchworm Monte Carlo method, which is proposed recently to tackle the numerical sign problem for open quantum systems. We focus on the growth of the numerical error with respect to the simulation time, for which the inchworm Monte Carlo method shows a flatter curve than the direct application of Monte Carlo method to the classical Dyson series. To better understand the underlying mechanism of the inchworm Monte Carlo method, we distinguish two types of exponential error growth, which are known as the numerical sign problem and the error amplification. The former is due to the fast growth of variance in the stochastic method, which can be observed from the Dyson series, and the latter comes from the evolution of the numerical solution. Our analysis demonstrates that the technique of partial resummation can be considered as a tool to balance these two types of error, and the inchworm Monte Carlo method is a successful case where the numerical sign problem is effectively suppressed by such means. We first demonstrate our idea in the context of ordinary differential equations, and then provide complete analysis for the inchworm Monte Carlo method. Several numerical experiments are carried out to verify our theoretical results. Keywords: Open quantum system, inchworm Monte Carlo method, numerical sign problem, error growth 1. Introduction In quantum mechanics, an open quantum system refers to a quantum system interacting with the envi- ronment. In reality, no quantum system is absolutely isolated, and therefore the theory of open quantum systems has wide applications including quantum thermodynamics [14], quantum information science [37], and quantum biology [1].
    [Show full text]