Large distance asymptotic behavior of the emptiness formation probability of the XXZ spin- 1 2 Heisenberg chain

Similar documents
Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

Non-interacting Spin-1/2 Particles in Non-commuting External Magnetic Fields

Convexity preserving interpolation by splines of arbitrary degree

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning

1 Matrix representations of canonical matrices

Tensor Smooth Length for SPH Modelling of High Speed Impact

Appendix B. The Finite Difference Scheme

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

APPENDIX A Some Linear Algebra

PHYS 705: Classical Mechanics. Canonical Transformation II

V.C The Niemeijer van Leeuwen Cumulant Approximation

THE WEIGHTED WEAK TYPE INEQUALITY FOR THE STRONG MAXIMAL FUNCTION

PHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions

Advanced Quantum Mechanics

A how to guide to second quantization method.

The Order Relation and Trace Inequalities for. Hermitian Operators

A Hybrid Variational Iteration Method for Blasius Equation

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lecture 12: Discrete Laplacian

Module 3: Element Properties Lecture 1: Natural Coordinates

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

International Conference on Advanced Computer Science and Electronics Information (ICACSEI 2013) equation. E. M. E. Zayed and S. A.

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Algebraic Solutions to Multidimensional Minimax Location Problems with Chebyshev Distance

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Binomial transforms of the modified k-fibonacci-like sequence

Formulas for the Determinant

Integrals and Invariants of Euler-Lagrange Equations

More metrics on cartesian products

Math 217 Fall 2013 Homework 2 Solutions

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Report on Image warping

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

Supplemental document

Difference Equations

12. The Hamilton-Jacobi Equation Michael Fowler

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Implicit Integration Henyey Method

Layer Potentials. Chapter Fundamental and Singular Solutions

Randić Energy and Randić Estrada Index of a Graph

Composite Hypotheses testing

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2

Finite Element Modelling of truss/cable structures

A Radon-Nikodym Theorem for Completely Positive Maps

CHAPTER 14 GENERAL PERTURBATION THEORY

PHYS 705: Classical Mechanics. Calculus of Variations II

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Dynamics of a Superconducting Qubit Coupled to an LC Resonator

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

Chapter 13: Multiple Regression

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

The Feynman path integral

Numerical Heat and Mass Transfer

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

The Second Anti-Mathima on Game Theory

EEE 241: Linear Systems

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

SUPPLEMENTARY INFORMATION

Numerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation

An Introduction to Morita Theory

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

14 The Postulates of Quantum mechanics

Lecture Notes 7: The Unruh Effect

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Indeterminate pin-jointed frames (trusses)

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Group Analysis of Ordinary Differential Equations of the Order n>2

Homework & Solution. Contributors. Prof. Lee, Hyun Min. Particle Physics Winter School. Park, Ye

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

The Geometry of Logit and Probit

LECTURE 9 CANONICAL CORRELATION ANALYSIS

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Appendix B. Criterion of Riemann-Stieltjes Integrability

Transcription:

LPENSL-TH-0/0 Large dstance asymptotc behavor of the emptness formaton probablty of the XXZ spn- Hesenberg chan N. Ktanne, J. M. Mallet, N. A. Slavnov 3, V. Terras 4 Abstract Usng ts multple ntegral representaton, we compute the large dstance asymptotc behavor of the emptness formaton probablty of the XXZ spn- Hesenberg chan n the massless regme. Graduate School of Mathematcal Scences, Unversty of Tokyo, Japan, ktanne@ms.u-tokyo.ac.jp On leave of absence from Steklov Insttute at St. Petersburg, Russa Laboratore de Physque, UMR 567 du CNRS, ENS Lyon, France, mallet@ens-lyon.fr 3 Steklov Mathematcal Insttute, Moscow, Russa, nslavnov@m.ras.ru 4 Department of Physcs and Astronomy, Rutgers Unversty, USA, vterras@physcs.rutgers.edu On leave of absence from LPMT, UMR 585 du CNRS, Montpeller, France

Emptness formaton probablty at large dstance The Hamltonan of the XXZ spn- Hesenberg chan s gven by H = M σ x m σm+ x + σy m σy m+ + σz m σz m+.. m= Here s the ansotropy parameter, σm x,y,z denote the usual Paul matrces actng on the quantum space at ste m of the chan. The emptness formaton probablty τm the probablty to fnd n the ground state a ferromagnetc strng of length m s defned as the followng expectaton value m σk z τm = ψ g ψ g,. k= where ψ g denotes the normalzed ground state. In the thermodynamc t M, ths quantty can be expressed as a multple ntegral wth m ntegratons [,, 3, 4, 5]. Recently, n the artcle [6], a new multple ntegral representaton for τm was obtaned. It leads n a drect way to the known answer at the free fermon pont = 0 [9], n partcular usng a saddle pont method, and to ts frst exact determnaton outsde the free fermon pont, namely at = [0]. The purpose of ths letter s to present the evaluaton of the asymptotc behavor of τm at large dstance m, n the massless regme < <, va the saddle pont method. We fnd log τm m m =log π + R 0 dω ω snh ω ω π cosh,.3 snh πω ω snh cosh ω where cos =,0<<π. If s commensurate wth π n other words f e s a root of unty, then the ntegral n.3 can be taken explctly n terms of ψ-functon logarthmc dervatve of Γ-functon. In partcular for = π and = π 3 respectvely = 0 and = / we obtan from.3 log τm m m = log, =0, log τm m m = 3 log 3 3log, =, whch concdes wth the known results obtaned respectvely n [7, 8, 9] and n [, 0]. For the partcular case of the XXX chan =, = 0 the asymptotc behavor can be evaluated also by the saddle pont method and t s gven by log τm m m =log Γ 3 4 Γ Γ 4.4 log0.599,.5

whch s n good agreement wth the known numercal result log0.598, obtaned n []. Below, we explan the man features of our method. A more detaled account of the proofs and technques nvolved wll be publshed later. The saddle pont method The multple ntegral representaton for τm obtaned n [6] can be wrtten n the form τm = sn m m m π D d m λ F {λ},m wth m a>b F {λ},m = snh π λ a λ b snhλ a λ b snhλ a λ b + ξ,...ξ m a>b m snhλa snhλ a + m,. a= cosh π λ a sn det m m.. snhλ snhξ a ξ b j ξ k snhλ j ξ k Here the ntegraton doman D s <λ <λ < <λ m <. Followng the standard arguments of the saddle pont method we estmate the ntegral. by the maxmal value of the ntegrand. Let {λ } be the set of parameters correspondng to ths maxmum. They satsfy the saddle pont equatons and for large m we assume that ther dstrbuton can be descrbed by a densty functon ρλ : ρλ j = m mλ j+.3 λ j. Thus for large m, one can replace sums over the set {λ } by ntegrals. ntegrable on the real axs, then Namely, f fλ s m m m fλ j j= m j= j k fλ j λ j λ k fλρλ dλ, V.P. fλ λ λ ρλ dλ, k m..4 Due to.4 t s easy to see that n the pont λ,...,λ m the products n the second lne of. behave as expc m. 3

Our goal s now to estmate the behavor of the term F {λ },m. To do ths we factorze the determnant n. as follows for large m: sn det m snhλ j ξ k snhλ j ξ k wth Indeed, for m one has = π m det m δ jk Kλ j λ k πmρλ k det m Kλ = det m δ jk Kλ j λ k πmρλ k det m snh π λ j ξ k snh π λ j ξ k,.5 sn snhλ snhλ +..6 m =det m snh π λ j ξ k Kλ j λ l πmρλ l snh π λ l ξ k l= det m snh π λ j ξ k Kλ j µ π dµ snh π µ ξ k m sn = det m π snhλ j ξ ksnhλ j ξ. k.7 Here we have used the fact that the functon / snh π λ j ξ solves the Leb ntegral equaton for the densty of the ground state of the XXZ magnet [3] and we have used the notatons of [6]. The second determnant n the r.h.s. of.5 s a Cauchy determnant, hence, F {λ },m= m π m +m m snh π λ a λ b a>b m cosh m π λ a a= det m δ jk Kλ j λ k πmρλ k..8 The behavor of the determnant n.8 can be estmated va Hadamard nequalty det m a jk max a jk m m m..9 appled to the above determnant and to the determnant of the nverse matrx, whch shows that m m log det m δ jk Kλ j λ k πmρλ k =0..0 4

The last equaton means that det m δ jk Kλ j λ k /πmρλ k does not contrbute to the leadng term of the asymptotcs. Hence, t can be excluded from our consderatons. Thus, up to subleadng correctons of the exponental type the emptness formaton probablty behaves as π m τm e m S 0, m,. wth S 0 S{λ }= m snh π m log λ a λ b snhλ a λ b snhλ a λ b + + m a= a>b m snhλ a / snhλ a + / log cosh π.. λ a Here the parameters {λ } are the solutons of the saddle pont equatons In our case the system.3 has the form π tanh πλ j = m m k= k j S 0 λ j cothλ j / cothλ j + / =0..3 π coth π λ j λ k cothλ j λ k cothλ j λ k +..4 Usng.4 we transform.4 nto the ntegral equaton for the densty ρλ π tanh πλ cothλ / cothλ + / π = V.P. coth π λ µ cothλ µ cothλ µ + ρµ dµ..5 Respectvely the acton S 0 takes the form snhλ / snhλ + / S 0 = dλρλlog cosh π λ + snh π λ µ dµdλρλρµlog..6 snhλ µ snhλ µ + 5

Snce the kernel of the ntegral operator n.5 depends on the dfference of the arguments, ths equaton can be solved va Fourer transform. Then ˆρω = Makng the nverse Fourer transform we fnd e ωλ ρλ dλ = cosh ω cosh ω..7 ρλ = πλ cosh cosh πλ,.8 whch obvously satsfes the needed normalsaton condton for densty ntegral on the real axs equals one. It remans to substtute.7,.8 nto.6, and after straghtforward calculatons we arrve at S 0 = R 0 dω ω snh ω ω π cosh..9 snh πω ω snh cosh ω Thus, we have obtaned.3. In the case of the XXX chan = one should rescale λ j λ j, ξ j ξ j n the orgnal multple ntegral representaton. for τm and then proceed to the t 0. The remanng computatons are then very smlar to the ones descrbed above, therefore we present here only the man results. The behavor of τm snowgvenby The acton S 0 n the saddle pont has the form λ /λ + / S 0 = log cosh ρλ dλ πλ π tanh πλ λ λ + 4 + τm π m e m S 0, m..0 = V.P. dµdλρλρµlog The analog of the ntegral equaton.5 n the XXX case s π coth πλ µ snh πλ µ λ µ λ µ + λ µ λ µ +.. ρµ dµ.. The soluton of ths equaton s ρλ = πλ cosh..3 coshπλ Substtutng.3 nto. we fnally arrve at.5. 6

Acknowledgments N. K. s supported by JSPS grant P077. N. K. would lke to thank M. Jmbo for help. N. S. s supported by the grants RFBR-0-0-00484, Foundaton of the Support of Russan Scence, Leadng Scentfc Schools 00-5-96046, the Program Nonlnear Dynamcs and Soltons and by CNRS. J.M. M. s supported by CNRS. V. T s supported by DOE grant DE-FG0-96ER40959 and by CNRS. N. K, N. S. and V. T. would lke to thank the Theoretcal Physcs group of the Laboratory of Physcs at ENS Lyon for hosptalty, whch makes ths collaboraton possble. We also would lke to thank the organzors of the 6th Internatonal Workshop Conformal Feld Theory and Integrable Models held n Chernogolovka, September 5-, 00, for the nce and stmulatng scentfc and extra-scentfc atmosphere they succeeded to generate. References [] M. Jmbo, K. Mk, T. Mwa and A. Nakayashk, Phys. Lett. A 68 99 56. [] M. Jmbo and T. Mwa, Journ. Phys. A: Math. Gen., 9 996 93. [3] M. Jmbo and T. Mwa, Algebrac analyss of solvable lattce models AMS, 995. [4] N. Ktanne, J. M. Mallet and V. Terras, Nucl. Phys. B, 554 [FS] 999 647, mathph/980700. [5] N. Ktanne, J. M. Mallet and V. Terras, Nucl. Phys. B, 567 [FS] 000 554; mathph/990709. [6] N. Ktanne, J. M. Mallet, N. A. Slavnov, V. Terras, Nucl. Phys. B 64 [FS] 00 487; hep-th/00045. [7] A. R. Its, A. G. Izergn, V. E. Korepn and N. A. Slavnov, Phys. Rev. Lett., 70 993 704. [8] M. Shrosh, M. Takahash and Y. Nshyama, Emptness Formaton Probablty for the One-Dmensonal Isotropc XY Model, cond-mat/00606. [9] N. Ktanne, J. M. Mallet, N. A. Slavnov, V. Terras, Correlaton functons of the XXZ- Hesenberg chan at the free fermon pont from ther multple ntegral representatons, hep-th/00369, to appear n Nuclear Physcs B, 00. [0] N. Ktanne, J. M. Mallet, N. A. Slavnov, V. Terras, J. Phys. A: Math. Gen. 35 00 L385-L388; hep-th/0034. [] A. V. Razumov and Yu. G. Stroganov, J. Phys. A: Math. Gen. 34 00 385-90; condmat/004. 7

[] H. E. Boos, V. E. Korepn, Y. Nshyama, M. Shrosh, Quantum Correlatons and Number Theory, cond-mat/0034. [3] E. Leb, T. Shultz and D. Matts, Ann. Phys., 6 96 407. 8