Random Transverse-field Ising Chain Research Articles

Critical behavior of dirty free parafermionic chains

Abstract A family of $\mathbb Z_n$-symmetric non-Hermitian models of Baxter was shown by Fendley to be exactly solvable via a parafermionic generalization of the Clifford algebra. We study these models with spatially random couplings, and obtain several exact results on thermodynamic singularities as the distributions of couplings are varied. We find that these singularities, independent of $n$, are identical to those in the random transverse-field Ising chain; correspondingly the models host infinite-randomness critical points. Similarities in structure to exact methods for random Ising models, a strong-disorder renormalization group, and generalizations to other models with free spectra, are discussed.

Spin- and flux-gap renormalization in the random Kitaev spin ladder

We study the Kitaev spin ladder with random couplings by using the real-space renormalization group technique. This model is the minimum model in Kitaev systems that has conserved plaquette fluxes, and its quasi-one-dimensional geometry makes it possible to study the strong-disorder fixed points for both spin- and flux- excitation gaps. In the Ising limit, the behavior of the spin gap is consistent with the familiar random transverse-field Ising chain, but the flux gap is dominated by the y-coupling. In the XX limit, while the x- and y-couplings are renormalized simultaneously, the z-couplings are not renormalized drastically and lead to non-universal disorder criticality at low-energy scales.

Critical quench dynamics of random quantum spin chains: ultra-slow relaxation from initial order and delayed ordering from initial disorder

By means of free fermionic techniques combined with multiple precision arithmetic we study the time evolution of the average magnetization, , of the random transverse-field Ising chain after global quenches. We observe different relaxation behaviors for quenches starting from different initial states to the critical point. Starting from a fully ordered initial state, the relaxation is logarithmically slow described by , and in a finite sample of length L the average magnetization saturates at a size-dependent plateau here the two exponents satisfy the relation . Starting from a fully disordered initial state, the magnetization stays at zero for a period of time until with and then starts to increase until it saturates to an asymptotic value , with . For both quenching protocols, finite-size scaling is satisfied in terms of the scaled variable . Furthermore, the distribution of long-time limiting values of the magnetization shows that the typical and the average values scale differently and the average is governed by rare events. The non-equilibrium dynamical behavior of the magnetization is explained through semi-classical theory.

Strong-disorder renormalization-group study of the one-dimensional tight-binding model

We formulate a strong-disorder renormalization-group (SDRG) approach to study the beta function of the tight-binding model in one dimension with both diagonal and off-diagonal disorder for states at the band center. We show that the SDRG method, when used to compute transport properties, yields exact results since it is identical to the transfer matrix method. The beta function is shown to be universal when only off-diagonal disorder is present even though single-parameter scaling is known to be violated. A different single-parameter scaling theory is formulated for this particular (particle-hole symmetric) case. Upon breaking particle-hole symmetry (by adding diagonal disorder), the beta function is shown to crossover from the universal behavior of the particle-hole symmetric case to the conventional non-universal one in agreement with the two-parameter scaling theory. We finally draw an analogy with the random transverse-field Ising chain in the paramagnetic phase. The particle-hole symmetric case corresponds to the critical point of the quantum Ising model while the generic case corresponds to the Griffiths paramagnetic phase.

Infinite-disorder critical points of models with stretched exponential interactions

We show that an interaction decaying as a stretched exponential function of distance, , is able to alter the universality class of short-range systems having an infinite-disorder critical point. To do so, we study the low-energy properties of the random transverse-field Ising chain with the above form of interaction by a strong-disorder renormalization group (SDRG) approach. We find that the critical behavior of the model is controlled by infinite-disorder fixed points different from those of the short-range model if 0 < a < 1/2. In this range, the critical exponents calculated analytically by a simplified SDRG scheme are found to vary with a, while, for a > 1/2, the model belongs to the same universality class as its short-range variant. The entanglement entropy of a block of size L increases logarithmically with L at the critical point but, unlike the short-range model, the prefactor is dependent on disorder in the range 0 < a < 1/2. Numerical results obtained by an improved SDRG scheme are found to be in agreement with the analytical predictions. The same fixed points are expected to describe the critical behavior of, among others, the random contact process with stretched exponentially decaying activation rates.

Random transverse-field Ising chain with long-range interactions

We study the low-energy properties of the long-range random transverse-field Ising chain with ferromagnetic interactions decaying as a power of the distance. Using variants of the strong-disorder renormalization group method, the critical behavior is found to be controlled by a strong-disorder fixed point with a finite dynamical exponent . Approaching the critical point, the correlation length diverges exponentially. In the critical point, the magnetization shows an -independent logarithmic finite-size scaling and the entanglement entropy satisfies the area law. These observations are argued to hold for other systems with long-range interactions, even in higher dimensions.

Distribution of dynamical quantities in the contact process, random walks, and quantum spin chains in random environments.

We study the distribution of dynamical quantities in various one-dimensional disordered models, the critical behavior of which is described by an infinite randomness fixed point. In the disordered contact process, the survival probability P(t) is found to show multiscaling in the critical point, meaning that P(t)=t-δ, where the (environment and time-dependent) exponent δ has a universal limit distribution when t→∞. The limit distribution is determined by the strong disorder renormalization group method analytically in the end point of a semi-infinite lattice, where it is found to be exponential, while, in the infinite system, conjectures on its limiting behaviors for small and large δ, which are based on numerical results, are formulated. By the same method, the survival probability in the problem of random walks in random environments is also shown to exhibit multiscaling with an exponential limit distribution. In addition to this, the (imaginary-time) spin-spin autocorrelation function of the random transverse-field Ising chain is found to have a form similar to that of survival probability of the contact process at the level of the renormalization approach. Consequently, a relationship between the corresponding limit distributions in the two problems can be established. Finally, the distribution of the spontaneous magnetization in this model is also discussed.

Rounding of a first-order quantum phase transition to a strong-coupling critical point

We investigate the effects of quenched disorder on first-order quantum phase transitions on the example of the $N$-color quantum Ashkin-Teller model. By means of a strong-disorder renormalization group, we demonstrate that quenched disorder rounds the first-order quantum phase transition to a continuous one for both weak and strong coupling between the colors. In the strong-coupling case, we find a distinct type of infinite-randomness critical point characterized by additional internal degrees of freedom. We investigate its critical properties in detail and find stronger thermodynamic singularities than in the random transverse field Ising chain. We also discuss the implications for higher spatial dimensions as well as unusual aspects of our renormalization-group scheme.

Dissipation effects in random transverse-field Ising chains

We study the effects of Ohmic, super-Ohmic, and sub-Ohmic dissipation on the zero-temperature quantum phase transition in the random transverse-field Ising chain by means of an (asymptotically exact) analytical strong-disorder renormalization-group approach. We find that Ohmic damping destabilizes the infinite-randomness critical point and the associated quantum Griffiths singularities of the dissipationless system. The quantum dynamics of large magnetic clusters freezes completely, which destroys the sharp phase transition by smearing. The effects of sub-Ohmic dissipation are similar and also lead to a smeared transition. In contrast, super-Ohmic damping is an irrelevant perturbation; the critical behavior is thus identical to that of the dissipationless system. We discuss the resulting phase diagrams, the behavior of various observables, and the implications to higher dimensions and experiments.

Entanglement entropy dynamics of disordered quantum spin chains

By means of free fermionic techniques we study the time evolution of the entanglement entropy, S(t), of a block of spins in the random transverse-field Ising chain after a sudden change of the parameters of the Hamiltonian. We consider global quenches, when the parameters are modified uniformly in space, as well as local quenches, when two disconnected blocks are suddenly joined together. For a non-critical final state, the dynamical entanglement entropy is found to approach a finite limiting value for both types of quenches. If the quench is performed to the critical state, the entropy grows for an infinite block as S(t) \sim ln ln t. This type of ultraslow increase is explained through the strong disorder renormalization group method.

Density of critical clusters in strips of strongly disordered systems

We consider two models with disorder-dominated critical points and study the distribution of clusters that are confined in strips and touch one or both boundaries. For the classical random bond Potts model in the large- q limit, we study optimal Fortuin-Kasteleyn clusters using a combinatorial optimization algorithm. For the random transverse-field Ising chain, clusters are defined and calculated through the strong-disorder renormalization group method. The numerically calculated density profiles close to the boundaries are shown to follow scaling predictions. For the random bond Potts model, we have obtained accurate numerical estimates for the critical exponents and demonstrated that the density profiles are well described by conformal formulas.

Theory of Smeared Quantum Phase Transitions

We present an analytical strong-disorder renormalization group theory of the quantum phase transition in the dissipative random transverse-field Ising chain. For Ohmic dissipation, we solve the renormalization flow equations analytically, yielding asymptotically exact results for the low-temperature properties of the system. We find that the interplay between quantum fluctuations and Ohmic dissipation destroys the quantum critical point by smearing. We also determine the phase diagram and the behavior of observables in the vicinity of the smeared quantum phase transition.

Finite-size scaling properties of random transverse-field Ising chains: Comparison between canonical and microcanonical ensembles for the disorder

The random transverse-field Ising chain is the simplest disordered model presenting a quantum phase transition at $T=0.$ We compare analytically its finite-size scaling properties in two different ensembles for the disorder: (i) the canonical ensemble, where the disorder variables are independent and (ii) the microcanonical ensemble, where there exists a global constraint on the disorder variables. The observables under study are the surface magnetization, the correlation of the two surface magnetizations, the gap and the end-to-end spin-spin correlation $C(L)$ for a chain of length L. At criticality, each observable decays typically as ${e}^{\ensuremath{-}w\sqrt{L}}$ in both ensembles, but the probability distributions of the rescaled variable w are different in the two ensembles, in particular, in their asymptotic behaviors. As a consequence, the dependence on L of averaged observables differs in the two ensembles. For instance, the correlation $C(L)$ decays algebraically as $1/L$ in the canonical ensemble, but subexponentially as ${e}^{\ensuremath{-}{\mathrm{cL}}^{1/3}}$ in the microcanonical ensemble. Off criticality, probability distributions of rescaled variables are governed by the critical exponent $\ensuremath{\nu}=2$ in both ensembles, but the following observables are governed by the exponent $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\nu}}=1$ in the microcanonical ensemble, instead of the exponent $\ensuremath{\nu}=2$ in the canonical ensemble: (a) In the disordered phase, the averaged surface magnetization, the averaged correlation of the two surface magnetizations and the averaged end-to-end spin-spin correlation; (b) in the ordered phase, the averaged gap. In conclusion, the measure of the rare events that dominate various averaged observables can be very sensitive to the microcanonical constraint.

Ensemble dependence in the random transverse-field Ising chain

In a disordered system one can either consider a microcanonical ensemble, where there is a precise constraint on the random variables, or a canonical ensemble where the variables are chosen according to a distribution without constraints. We address the question as to whether critical exponents in these two cases can differ through a detailed study of the random transverse-field Ising chain. We find that the exponents are the same in both ensembles, though some critical amplitudes vanish in the microcanonical ensemble for correlations which span the whole system and are particularly sensitive to the constraint. This can \textit{appear} as a different exponent. We expect that this apparent dependence of exponents on ensemble is related to the integrability of the model, and would not occur in non-integrable models.

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

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.

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

We study numerically the paramagnetic phase of the spin-1/2 random transverse-field Ising chain, using a mapping to non-interacting fermions. We extend our earlier work, Phys. Rev. 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(\delta)$, where $\delta$ measures the deviation from criticality. There are some discrepancies between the values of $z(\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 $\tau^{-1/z(\delta)}$, where $\tau$ is imaginary time. However, the typical value decays with a stretched exponential behavior, $\exp(-c\tau^{1/\mu})$, where $\mu$ may be related to $z(\delta)$. We also obtain results for the full probability distribution of time dependent correlation functions at different points in the paramagnetic phase.

Numerical study of the random transverse-field Ising spin chain.

We study numerically the critical region and the disordered phase of the random transverse-field Ising chain. By using a mapping of Lieb, Schultz, and Mattis to noninteracting fermions, we can obtain a numerically exact solution for rather large system sizes, L\ensuremath{\le}128. Our results confirm the striking predictions of earlier analytical work and, in addition, give results for some probability distributions and scaling functions. \textcopyright{} 1996 The American Physical Society.