Dynamical correlations of the spin and the energy density are investigated in the critical region of the random transverse-field Ising chain by numerically exact calculations in large finite systems (L ≤ 128). The spin-spin autocorrelation function is found to decay proportional to (log t) −2xm and (log t) −2xs m in the bulk and on the surface, respectively, with xm and x s m the bulk and surface magnetization exponents, respectively. On the other hand the critical energy-energy autocorrelation functions have a power law decay, which are characterized by novel critical exponents ηe ≈ 2.2 in the bulk and η s e ≈ 2.5 at the surface, respectively. The numerical results are compared with the predictions of a scaling theory. 05.50.+q, 64.60.Ak, 68.35.Rh Typeset using REVTEX 1 The asymptotic behavior of the time-dependent correlation functions for interacting many-body systems turned out a very difficult subject of theoretical research. Exact results in this field are scarce, one can mention the one-dimensional spin 1/2 XY-model [1] and the Ising chain in a transverse-field [2]. Both models can be mapped onto a system of non-interacting fermions, where the equal-position correlation functions are calculated by the Pfaffian method utilizing the theory of Töplitz determinants. In this Letter we consider- at the first time- the critical dynamical correlations of an interacting quantum system in the presence of quenched (i.e. time-dependent) disorder. It has recently become clear that quenched disorder has rather different effects on phase transitions in quantum systems [3] than on those thermally driven phase transitions. For example, in the Griffiths phase, which is situated at the disordered side of the critical point, the susceptibility has an essential singularity in classical systems, whereas in a quantum system the corresponding singularity is stronger, it is in a power law form. Here we consider the prototype of random quantum systems the one-dimensional random transverse-field Ising model defined by the Hamiltonian:Dynamical correlations of the spin and the energy density are investigated in the critical region of the random transverse-field Ising chain by numerically exact calculations in large finite systems (L ≤ 128). The spin-spin autocorrelation function is found to decay proportional to (log t) −2xm and (log t) −2xs m in the bulk and on the surface, respectively, with xm and x s m the bulk and surface magnetization exponents, respectively. On the other hand the critical energy-energy autocorrelation functions have a power law decay, which are characterized by novel critical exponents ηe ≈ 2.2 in the bulk and η s e ≈ 2.5 at the surface, respectively. The numerical results are compared with the predictions of a scaling theory. 05.50.+q, 64.60.Ak, 68.35.Rh Typeset using REVTEX 1 The asymptotic behavior of the time-dependent correlation functions for interacting many-body systems turned out a very difficult subject of theoretical research. Exact results in this field are scarce, one can mention the one-dimensional spin 1/2 XY-model [1] and the Ising chain in a transverse-field [2]. Both models can be mapped onto a system of non-interacting fermions, where the equal-position correlation functions are calculated by the Pfaffian method utilizing the theory of Töplitz determinants. In this Letter we consider- at the first time- the critical dynamical correlations of an interacting quantum system in the presence of quenched (i.e. time-dependent) disorder. It has recently become clear that quenched disorder has rather different effects on phase transitions in quantum systems [3] than on those thermally driven phase transitions. For example, in the Griffiths phase, which is situated at the disordered side of the critical point, the susceptibility has an essential singularity in classical systems, whereas in a quantum system the corresponding singularity is stronger, it is in a power law form. Here we consider the prototype of random quantum systems the one-dimensional random transverse-field Ising model defined by the Hamiltonian:

Jan 1, 19978 citations

Dynamical critical properties of the random transverse-field Ising spin chain

Jens Kisker, A. P. Young

We study the dynamical properties of the random transverse-field Ising chain at criticality using a mapping to free fermions with which we can obtain numerically exact results for system sizes L as large as 256. The probability distribution of the local imaginary time correlation function $S(\ensuremath{\tau})$ is investigated and found to be simply a function of $\ensuremath{\alpha}\ensuremath{\equiv}\ensuremath{-}\mathrm{ln}S(\ensuremath{\tau})/\mathrm{ln}\ensuremath{\tau}.$ This scaling behavior implies that the typical correlation function decays algebraically, ${S}_{\mathrm{typ}}(\ensuremath{\tau})\ensuremath{\sim}{\ensuremath{\tau}}^{\ensuremath{-}{\ensuremath{\alpha}}_{\mathrm{typ}}},$ where the exponent ${\ensuremath{\alpha}}_{\mathrm{typ}}$ is determined from $P(\ensuremath{\alpha}),$ the distribution of $\ensuremath{\alpha}.$ The precise value for ${\ensuremath{\alpha}}_{\mathrm{typ}}$ depends on exactly how the ``typical'' correlation function is defined. The form of $P(\ensuremath{\alpha})$ for small $\ensuremath{\alpha}$ gives a contribution to the average correlation function ${S}_{\mathrm{av}}(\ensuremath{\tau})$ namely, ${S}_{\mathrm{av}}(\ensuremath{\tau})\ensuremath{\sim}(\mathrm{ln}\ensuremath{\tau}{)}^{\ensuremath{-}{2x}_{m}},$ where ${x}_{m}$ is the bulk magnetization exponent, which was reported recently in Europhys. Lett. 39, 135 (1997). These results represent a type of ``multiscaling'' different from the well-known ``multifractal'' behavior.We study the dynamical properties of the random transverse-field Ising chain at criticality using a mapping to free fermions with which we can obtain numerically exact results for system sizes L as large as 256. The probability distribution of the local imaginary time correlation function $S(\ensuremath{\tau})$ is investigated and found to be simply a function of $\ensuremath{\alpha}\ensuremath{\equiv}\ensuremath{-}\mathrm{ln}S(\ensuremath{\tau})/\mathrm{ln}\ensuremath{\tau}.$ This scaling behavior implies that the typical correlation function decays algebraically, ${S}_{\mathrm{typ}}(\ensuremath{\tau})\ensuremath{\sim}{\ensuremath{\tau}}^{\ensuremath{-}{\ensuremath{\alpha}}_{\mathrm{typ}}},$ where the exponent ${\ensuremath{\alpha}}_{\mathrm{typ}}$ is determined from $P(\ensuremath{\alpha}),$ the distribution of $\ensuremath{\alpha}.$ The precise value for ${\ensuremath{\alpha}}_{\mathrm{typ}}$ depends on exactly how the ``typical'' correlation function is defined. The form of $P(\ensuremath{\alpha})$ for small $\ensuremath{\alpha}$ gives a contribution to the average correlation function ${S}_{\mathrm{av}}(\ensuremath{\tau})$ namely, ${S}_{\mathrm{av}}(\ensuremath{\tau})\ensuremath{\sim}(\mathrm{ln}\ensuremath{\tau}{)}^{\ensuremath{-}{2x}_{m}},$ where ${x}_{m}$ is the bulk magnetization exponent, which was reported recently in Europhys. Lett. 39, 135 (1997). These results represent a type of ``multiscaling'' different from the well-known ``multifractal'' behavior.

arXiv:cond-mat/9807025Physical review. B, Condensed matterDec 1, 199816 citations

Griffiths-McCoy singularities in the random transverse-field Ising spin chain

Ferenc Iglói, Róbert Juhász, Heiko Rieger

We consider the paramagnetic phase of the random transverse-field Ising spin chain and study the dynamical properties by numerical methods and scaling considerations. We extend our previous work [Phys. Rev. B 57, 11 404 (1998)] to new quantities, such as the nonlinear susceptibility, higher excitations, and the energy-density autocorrelation function. We show that in the Griffiths phase all the above quantities exhibit power-law singularities and the corresponding critical exponents, which vary with the distance from the critical point, can be related to the dynamical exponent z, the latter being the positive root of $[{(J/h)}^{1/z}{]}_{\mathrm{av}}=1$ . Particularly, whereas the average spin autocorrelation function in imaginary time decays as $[G{]}_{\mathrm{av}}(\ensuremath{\tau})\ensuremath{\sim}{\ensuremath{\tau}}^{\ensuremath{-}1/z}$, the average energy-density autocorrelations decay with another exponent as $[{G}^{e}{]}_{\mathrm{av}}(\ensuremath{\tau})\ensuremath{\sim}{\ensuremath{\tau}}^{\ensuremath{-}2\ensuremath{-}1/z}$.We consider the paramagnetic phase of the random transverse-field Ising spin chain and study the dynamical properties by numerical methods and scaling considerations. We extend our previous work [Phys. Rev. B 57, 11 404 (1998)] to new quantities, such as the nonlinear susceptibility, higher excitations, and the energy-density autocorrelation function. We show that in the Griffiths phase all the above quantities exhibit power-law singularities and the corresponding critical exponents, which vary with the distance from the critical point, can be related to the dynamical exponent z, the latter being the positive root of $[{(J/h)}^{1/z}{]}_{\mathrm{av}}=1$. Particularly, whereas the average spin autocorrelation function in imaginary time decays as $[G{]}_{\mathrm{av}}(\ensuremath{\tau})\ensuremath{\sim}{\ensuremath{\tau}}^{\ensuremath{-}1/z}$, the average energy-density autocorrelations decay with another exponent as $[{G}^{e}{]}_{\mathrm{av}}(\ensuremath{\tau})\ensuremath{\sim}{\ensuremath{\tau}}^{\ensuremath{-}2\ensuremath{-}1/z}$.

arXiv:cond-mat/9811369Physical review. B, Condensed matterMay 1, 199919 citations

Transverse-field Ising spin chain with inhomogeneous disorder

Dragi Karevski, Róbert Juhász, Loïc Turban +1 more

We consider the critical and off-critical properties at the boundary of the random transverse-field Ising spin chain when the distribution of the couplings and/or transverse fields, at a distance l from the surface, deviates from its uniform bulk value by terms of order ${l}^{\ensuremath{-}\ensuremath{\kappa}}$ with an amplitude A. Exact results are obtained using a correspondence between the surface magnetization of the model and the surviving probability of a random walk with time-dependent absorbing boundary conditions. For slow enough decay, $\ensuremath{\kappa}<1/2,$ the inhomogeneity is relevant: Either the surface stays ordered at the bulk critical point or the average surface magnetization displays an essential singularity, depending on the sign of A. In the marginal situation $\ensuremath{\kappa}=1/2,$ the average surface magnetization decays as a power law with a continuously varying, A-dependent, critical exponent which is obtained analytically. The behavior of the critical and off-critical autocorrelation functions as well as the scaling form of the probability distributions for the surface magnetization and the first gaps are determined through a phenomenological scaling theory. In the Griffiths phase, the properties of the Griffiths-McCoy singularities are not affected by the inhomogeneity. The various results are checked using numerical methods based on a mapping to free fermions.We consider the critical and off-critical properties at the boundary of the random transverse-field Ising spin chain when the distribution of the couplings and/or transverse fields, at a distance l from the surface, deviates from its uniform bulk value by terms of order ${l}^{\ensuremath{-}\ensuremath{\kappa}}$ with an amplitude A. Exact results are obtained using a correspondence between the surface magnetization of the model and the surviving probability of a random walk with time-dependent absorbing boundary conditions. For slow enough decay, $\ensuremath{\kappa}<1/2,$ the inhomogeneity is relevant: Either the surface stays ordered at the bulk critical point or the average surface magnetization displays an essential singularity, depending on the sign of A. In the marginal situation $\ensuremath{\kappa}=1/2,$ the average surface magnetization decays as a power law with a continuously varying, A-dependent, critical exponent which is obtained analytically. The behavior of the critical and off-critical autocorrelation functions as well as the scaling form of the probability distributions for the surface magnetization and the first gaps are determined through a phenomenological scaling theory. In the Griffiths phase, the properties of the Griffiths-McCoy singularities are not affected by the inhomogeneity. The various results are checked using numerical methods based on a mapping to free fermions.

arXiv:cond-mat/9903388Physical review. B, Condensed matterAug 1, 19998 citations

Finite-temperature and dynamical properties of the random transverse-field Ising spin chain

A. P. Young

We study numerically the paramagnetic phase of the spin- $1/2$ random transverse-field Ising chain, using a mapping to noninteracting fermions. We extend our earlier work, Phys. Rev. B 53, 8486 (1996), to finite temperatures and to dynamical properties. Our results are consistent with the idea that there are Griffiths-McCoy singularities in the paramagnetic phase described by a continuously varying exponent $z(\ensuremath{\delta}),$ where \ensuremath{\delta} measures the deviation from criticality. There are some discrepancies between the values of $z(\ensuremath{\delta})$ obtained from different quantities, but this may be due to corrections to scaling. The average on-site time dependent correlation function decays with a power law in the paramagnetic phase, namely, ${\ensuremath{\tau}}^{\ensuremath{-}1/z(\ensuremath{\delta})},$ where \ensuremath{\tau} is imaginary time. However, the typical value decays with a stretched exponential behavior, $\mathrm{exp}(\ensuremath{-}c{\ensuremath{\tau}}^{1/\ensuremath{\mu}}),$ where \ensuremath{\mu} may be related to $z(\ensuremath{\delta}).$ We also obtain results for the full probability distribution of time dependent correlation functions at different points in the paramagnetic phase.We study numerically the paramagnetic phase of the spin-$1/2$ random transverse-field Ising chain, using a mapping to noninteracting fermions. We extend our earlier work, Phys. Rev. B 53, 8486 (1996), to finite temperatures and to dynamical properties. Our results are consistent with the idea that there are Griffiths-McCoy singularities in the paramagnetic phase described by a continuously varying exponent $z(\ensuremath{\delta}),$ where \ensuremath{\delta} measures the deviation from criticality. There are some discrepancies between the values of $z(\ensuremath{\delta})$ obtained from different quantities, but this may be due to corrections to scaling. The average on-site time dependent correlation function decays with a power law in the paramagnetic phase, namely, ${\ensuremath{\tau}}^{\ensuremath{-}1/z(\ensuremath{\delta})},$ where \ensuremath{\tau} is imaginary time. However, the typical value decays with a stretched exponential behavior, $\mathrm{exp}(\ensuremath{-}c{\ensuremath{\tau}}^{1/\ensuremath{\mu}}),$ where \ensuremath{\mu} may be related to $z(\ensuremath{\delta}).$ We also obtain results for the full probability distribution of time dependent correlation functions at different points in the paramagnetic phase.

arXiv:cond-mat/9707060Physical review. B, Condensed matterNov 1, 199758 citations

Random transverse Ising spin chain and random walks

Ferenc Iglói, Heiko Rieger

We study the critical and off-critical (Griffiths-McCoy) regions of the random transverse-field Ising spin chain by analytical and numerical methods and by phenomenological scaling considerations. Here we extend previous investigations to surface quantities and to the ferromagnetic phase. The surface magnetization of the model is shown to be related to the surviving probability of an adsorbing walk and several critical exponents are exactly calculated. Analyzing the structure of low-energy excitations we present a phenomenological theory which explains both the scaling behavior at the critical point and the nature of Griffiths-McCoy singularities in the off-critical regions. In the numerical part of the work we used the free-fermion representation of the model and calculated the critical magnetization profiles, which are found to follow very accurately the conformal predictions for different boundary conditions. In the off-critical regions we demonstrated that the Griffiths-McCoy singularities are characterized by a single, varying exponent, the value of which is related through duality in the paramagnetic and ferromagnetic phases.We study the critical and off-critical (Griffiths-McCoy) regions of the random transverse-field Ising spin chain by analytical and numerical methods and by phenomenological scaling considerations. Here we extend previous investigations to surface quantities and to the ferromagnetic phase. The surface magnetization of the model is shown to be related to the surviving probability of an adsorbing walk and several critical exponents are exactly calculated. Analyzing the structure of low-energy excitations we present a phenomenological theory which explains both the scaling behavior at the critical point and the nature of Griffiths-McCoy singularities in the off-critical regions. In the numerical part of the work we used the free-fermion representation of the model and calculated the critical magnetization profiles, which are found to follow very accurately the conformal predictions for different boundary conditions. In the off-critical regions we demonstrated that the Griffiths-McCoy singularities are characterized by a single, varying exponent, the value of which is related through duality in the paramagnetic and ferromagnetic phases.

arXiv:cond-mat/9709260Physical review. B, Condensed matterMay 1, 1998105 citations

{x}_{m}

is the bulk magnetization exponent, which was reported recently in Europhys. Lett. 39, 135 (1997). These results represent a type of ``multiscaling'' different from the well-known ``multifractal'' behavior.

[{(J/h)}^{1/z}{]}_{\mathrm{av}}=1

. Particularly, whereas the average spin autocorrelation function in imaginary time decays as $[G{]}_{\mathrm{av}}(\ensuremath{\tau})\ensuremath{\sim}{\ensuremath{\tau}}^{\ensuremath{-}1/z}$, the average energy-density autocorrelations decay with another exponent as $[{G}^{e}{]}_{\mathrm{av}}(\ensuremath{\tau})\ensuremath{\sim}{\ensuremath{\tau}}^{\ensuremath{-}2\ensuremath{-}1/z}$.

We study numerically the paramagnetic phase of the spin-

1/2

random transverse-field Ising chain, using a mapping to noninteracting fermions. We extend our earlier work, Phys. Rev. B 53, 8486 (1996), to finite temperatures and to dynamical properties. Our results are consistent with the idea that there are Griffiths-McCoy singularities in the paramagnetic phase described by a continuously varying exponent $z(\ensuremath{\delta}),$ where \ensuremath{\delta} measures the deviation from criticality. There are some discrepancies between the values of $z(\ensuremath{\delta})$ obtained from different quantities, but this may be due to corrections to scaling. The average on-site time dependent correlation function decays with a power law in the paramagnetic phase, namely, ${\ensuremath{\tau}}^{\ensuremath{-}1/z(\ensuremath{\delta})},$ where \ensuremath{\tau} is imaginary time. However, the typical value decays with a stretched exponential behavior, $\mathrm{exp}(\ensuremath{-}c{\ensuremath{\tau}}^{1/\ensuremath{\mu}}),$ where \ensuremath{\mu} may be related to $z(\ensuremath{\delta}).$ We also obtain results for the full probability distribution of time dependent correlation functions at different points in the paramagnetic phase.