Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1805.03464

Universal scaling relationship between classical and quantum correlations in critical quantum spin chains

Yan-Wei Dai, Xi-Hao Chen, Sam Young Cho, Huan-Qiang Zhou, Dao‐Xin Yao

Published May 9, 2018• 2 citations
arXivHTMLPDF

Abstract

We numerically investigate classical and quantum correlations in one-dimensional quantum critical systems. The infinite matrix product state (iMPS) representation is employed in order to consider an infinite-size spin chain. By using the infinite time-evolving block decimation algorithm, iMPS ground state wave functions are obtained at critical points for the transverse-field spin-$1/2$ XY model. From the ground state wave functions, we calculate classical and quantum correlations and mutual information. All of the correlations are found to exhibit a power-law decay with the increments of the lattice distance for both the transition lines of the Ising universality class and the Gaussian universality class. Such power-law scaling behaviors of the correlations manifest the existence of diversing correlation lengths, which means scale invariance. The critical features of the correlations can be characterized by introducing a critical exponent of the power-law decaying correlations. Similar to the critical exponent $η$ of the spin-spin correlation for the universality classes in the transverse-field XY model, we calculate the critical exponents of the two-spin classical and quantum correlations as well as that of the corresponding mutual information. All of the correlations have the same critical exponents, i.e., $η^{I}=η^{C}=η^{D}$ at a critical point, where the superscripts $I$, $C$, and $D$ stand for mutual information, classical correlation, and quantum correlation, respectively. Furthermore, the critical exponent $η$ of the spin-spin correlation is shown to relate to $η= η^α/2$ with $α\in \{ I, C, D\}$.

DOI
10.48550/arxiv.1805.03464

Related papers

Critical exponents of block-block mutual information in one-dimensional infinite lattice systems

Yan-Wei Dai, Xi Hao Chen, Sam Young Cho +1 more

We study the mutual information between two lattice blocks in terms of von Neumann entropies for one-dimensional infinite lattice systems. Quantum q-state Potts model and transverse-field spin-1/2 XY model are considered numerically by using the infinite matrix product state approach. As a system parameter varies, block-block mutual information exhibit singular behaviors that enable us to identify the critical points for the quantum phase transitions. As happens with von Neumann entanglement entropy of single block, at critical points, block-block mutual information for two adjacent blocks show a logarithmic leading behavior with increasing the size of the blocks, which yields the central charge c of the underlying conformal field theory, as it should be. It seems that two disjoint blocks show a similar logarithmic growth of the mutual information as a characteristic property of critical systems but the proportional coefficients of the logarithmic term are very different from the central charges. As the separation between the two lattice blocks increases, the mutual information reveals a consistent power-law decaying behavior for various truncation dimensions and lattice-block sizes. The critical exponent of block-block mutual information in the thermodynamic limit is estimated by extrapolating the exponents of power-law decaying regions for finite truncation dimensions. For a given lattice-block size ℓ, the critical exponents for the same universality classes seem to have very close values each other. Whereas the critical exponents have different values to a degree of distinction for the different universality classes. As the lattice-block size becomes bigger, the critical exponent becomes smaller. We discuss a relation between the exponents of block-block mutual information and correlation with the Shatten one-norm of block-block correlation.

arXiv:2005.07924•Physical review. E•Oct 29, 2021•5 citations

Scaling of entanglement support for matrix product states

Luca Tagliacozzo, Thiago Rodrigues de Oliveira, S. Iblisdir +1 more

The power of matrix product states to describe infinite-size translational-invariant critical spin chains is investigated. At criticality, the accuracy with which they describe ground-state properties of a system is limited by the size $\ensuremath{\chi}$ of the matrices that form the approximation. This limitation is quantified in terms of the scaling of the half-chain entanglement entropy. In the case of the quantum Ising model, we find $S\ensuremath{\sim}\frac{1}{6}\text{log}\text{ }\ensuremath{\chi}$ with high precision. This result can be understood as the emergence of an effective finite correlation length ${\ensuremath{\xi}}_{\ensuremath{\chi}}$ ruling all the scaling properties in the system. We produce six extra pieces of evidence for this finite-$\ensuremath{\chi}$ scaling, namely, the scaling of the correlation length, the scaling of magnetization, the shift of the critical point, the scaling of the entanglement entropy for a finite block of spins, the existence of scaling functions, and the agreement with analogous classical results. All our computations are consistent with a scaling relation of the form ${\ensuremath{\xi}}_{\ensuremath{\chi}}\ensuremath{\sim}{\ensuremath{\chi}}^{\ensuremath{\kappa}}$, with $\ensuremath{\kappa}=2$ for the Ising model. In the case of the Heisenberg model, we find similar results with the value $\ensuremath{\kappa}\ensuremath{\sim}1.37$. We also show how finite-$\ensuremath{\chi}$ scaling allows us to extract critical exponents. These results are obtained using the infinite time evolved block decimation algorithm which works in the thermodynamical limit and are verified to agree with density-matrix renormalization-group results and their classical analog obtained with the corner transfer-matrix renormalization group.

arXiv:0712.1976•Physical Review B•Jul 14, 2008•314 citations

Non-local Correlations in the Haldane Phase for an XXZ Spin-1 Chain: A Perspective from Infinite Matrix Product State Representation

Yao Heng Su, Sam Young Cho, Bo Li +2 more

String correlations are investigated in an infinite-size XXZ spin-1 chain. By using the infinite matrix product state representation, we calculate a long-range string order. In the XY phase, the string correlations decay within a relatively very large lattice distance, which makes a finite-size study difficult to verify the non-existence of the string order. Thus, in the Haldane phase, the non-vanishing string correlations in the limit of a very large distance allow to characterize the phase boundaries to the XY phase and the Neel phase, which implies that the transverse long-range string order is the order parameter for the Haldane phase. In addition, the singular behaviors of the von Neumann entropy and the fidelity per lattice site are shown to capture clearly the phase transition points that are consistent with the results from the string order. The estimated critical points including a BKT transition from the XY phase to the Haldane phase agree well with the previous results: $Δ_{c2} = 0$ for the XY-Haldane phase transition and $Δ_{c3} = 1.185$ for the Haldane-N'eel phase transition from the density renormalization group. From a finite-entanglement scaling of the von Neumann entropy with respect to the truncation dimension, the central charges are found to be $c \simeq 1.0$ at $Δ_{c2} = 0$ and $c \simeq 0.5$ at $Δ_{c3} = 1.185$, respectively, which shows that the XY-Haldane phase transition at $Δ_{c2} = 0$ belongs to the Heisenberg universality class, while the Haldane-Neel phase transition at $Δ_{c2} = 1.185$ belongs to the two-dimensional classical Ising universality class. It is also shown that, the long-range order parameters and the von Neumann entropy, as well as the fidelity per site approach, can be applied to characterize quantum phase transitions as a universal phase transition indicator for one-dimensional lattice many-body systems.

arXiv:1202.6147•Journal of the Physical Society of Japan•Jun 20, 2012•36 citations

Global geometric entanglement in transverse-field XY spin chains: finite and infinite systems

Tzu-Chieh Wei, Smitha Vishveshwara, Paul M. Goldbart

The entanglement in quantum XY spin chains of arbitrary length is investigated via the geometric measure of entanglement. The emergence of entanglement is explained intuitively from the perspective of perturbations. The model is solved exactly and the energy spectrum is determined and analyzed in particular for the lowest two levels for both finite and infinite systems. The overlaps for these two levels are calculated analytically for arbitrary number of spins. The entanglement is hence obtained by maximizing over a single parameter. The corresponding ground-state entanglement surface is then determined over the entire phase diagram, and its behavior can be used to delineate the boundaries in the phase diagram. For example, the field-derivative of the entanglement becomes singular along the critical line. The form of the divergence is derived analytically and it turns out to be dictated by the universality class controlling the quantum phase transition. The behavior of the entanglement near criticality can be understood via a scaling hypothesis, analogous to that for free energies. The entanglement density vanishes along the so-called disorder line in the phase diagram, the ground space is doubly degenerate and spanned by two product states. The entanglement for the superposition of the lowest two states is also calculated. The exact value of the entanglement depends on the specific form of superposition. However, in the thermodynamic limit the entanglement density turns out to be independent of the superposition. This proves that the entanglement density is insensitive to whether the ground state is chosen to be the spontaneously $Z_2$ symmetry broken one or not. The finite-size scaling of entanglement at critical points is also investigated from two different view points. First, the maximum in the field-derivative of the entanglement density is computed and fitted to a logarithmic dependence of the system size, thereby deducing the correlation length exponent for the Ising class using only the behavior of entanglement. Second, the entanglement density itself is shown to possess a correction term inversely proportional to the system size, with the coefficient being universal (but with different values for the ground state and the first excited state, respectively).

arXiv:1012.4114•Quantum Information and Computation•Mar 1, 2011•10 citations

Quantum entanglement and critical exponents in one-dimensional spin-1 bond-alternating XXZ chains

Yaoheng Su, Chen Ai-Min, Honglei Wang +1 more

The characterization of the quantum phase transition in a lowdimensional system has attracted a considerable amount of attention in quantum manybody systems. As one of the fundamental models in quantum magnetism, spin-1 models have richer phase diagrams and show more complex physical phenomena. In the spin-1 antiferromagnetic XXZ model, the Haldane phase and the Nel phase are the gapped topologic phases which cannot be characterized by the local order parameters. To characterize the nature in such phases, one has to calculate the non-local long range order parameters. Normally, the non-local order parameter in the topological phase is obtained from the extrapolation of finite-sized system in numerical study. However, it is difficult to extract the critical exponents with such an extrapolated non-local order parameter due to the numerical accuracy. In a recently developed tensor network representation, i.e., the infinite matrix product state (iMPS) algorithm, it was shown that the non-local order can be directly calculated from a very large lattice distance in an infinite-sized system rather than an extrapolated order parameter in a finite-sized system. Therefore, it is worthwhile using this convenient technique to study the non-local orders in the topological phases and characterize the quantum criticalities in the topological quantum phase transitions. In this paper, by utilizing the infinite matrix product state algorithm based on the tensor network representation and infinite time evolving block decimation method, the quantum entanglement, fidelity, and critical exponents of the topological phase transition are investigated in the one-dimensional infinite spin-1 bond-alternating XXZ Heisenberg model. It is found that there is always a local dimerization order existing in the whole parameter range when the bond-alternative strength parameter changes from 0 to 1. Also, due to the effect of the bond-alternating, there appears a quantum phase transition from the long-rang ordering topological Nel phase to the local ordering dimerization phase. The von Neumann entropy, fidelity per lattice site, and order parameters all give the same phase transition point at c = 0.638. To identify the type of quantum phase transition, the central charge c 0.5 is extracted from the ground state von Neumann entropy and the finite correlation length, which indicates that the phase transition belongs to the two-dimensional Ising universality class. Furthermore, it is found that the Nel order and the susceptibility of Nel order have power-law relations to |-c|. From the numerical fitting of the Nel order and its susceptibility, we obtain the characteristic critical exponents ' = 0.236 and ' = 0.838. It indicates that such critical exponents from our method characterize the nature of the quantum phase transition. Our critical exponents from the iMPS method can provide guidance for studying the properties of the phase transition in quantum spin systems.

Acta Physica Sinica•Jan 1, 2017•2 citations

Universal short-time quantum critical dynamics of finite-size systems

Yu-Rong Shu, Shuai Yin, Dao‐Xin Yao

We investigate the short-time quantum critical dynamics in the imaginary-time relaxation processes of finite-size systems. Universal scaling behaviors exist in the imaginary-time evolution. In particular, the system undergoes a critical initial slip stage characterized by an exponent $\ensuremath{\theta}$, in which an initial power-law increase emerges in the imaginary-time correlation function when the initial state has a zero order parameter and a vanishing correlation length. Under different initial conditions, the quantum critical point and critical exponents can be determined from the universal scaling behaviors. We apply the method to the one- and two-dimensional transverse field Ising models using quantum Monte Carlo (QMC) simulations. In the one-dimensional case, we locate the quantum critical point at ${(h/J)}_{c}=1.000\phantom{\rule{0.16em}{0ex}}03(8)$ in the thermodynamic limit, and we estimate the critical initial slip exponent $\ensuremath{\theta}=0.3734(2)$ and the static exponent $\ensuremath{\beta}/\ensuremath{\nu}=0.1251(2)$ by analyzing data on chains of length $L=32--256$ and 48--256, respectively. For the two-dimensional square-lattice system, the critical coupling ratio is given by $3.044\phantom{\rule{0.16em}{0ex}}51(7)$ in the thermodynamic limit, while the critical exponents are $\ensuremath{\theta}=0.209(4)$ and $\ensuremath{\beta}/\ensuremath{\nu}=0.518(1)$ estimated by data on systems of size $L=24$--64 and 32--64, respectively. Remarkably, the critical initial slip exponents obtained in both models are notably distinct from their classical counterparts due to the essential differences between classical and quantum dynamics. The short-time critical dynamics and the imaginary-time relaxation QMC approach can be readily adapted to various models.

arXiv:1705.05931•Physical review. B./Physical review. B•Sep 12, 2017•17 citations
We numerically investigate classical and quantum correlations in one-dimensional quantum critical systems. The infinite matrix product state (iMPS) representation is employed in order to consider an infinite-size spin chain. By using the infinite time-evolving block decimation algorithm, iMPS ground state wave functions are obtained at critical points for the transverse-field spin-1/21/21/2 XY model. From the ground state wave functions, we calculate classical and quantum correlations and mutual information. All of the correlations are found to exhibit a power-law decay with the increments of the lattice distance for both the transition lines of the Ising universality class and the Gaussian universality class. Such power-law scaling behaviors of the correlations manifest the existence of diversing correlation lengths, which means scale invariance. The critical features of the correlations can be characterized by introducing a critical exponent of the power-law decaying correlations. Similar to the critical exponent ηηη of the spin-spin correlation for the universality classes in the transverse-field XY model, we calculate the critical exponents of the two-spin classical and quantum correlations as well as that of the corresponding mutual information. All of the correlations have the same critical exponents, i.e., ηI=ηC=ηDη^{I}=η^{C}=η^{D}ηI=ηC=ηD at a critical point, where the superscripts III, CCC, and DDD stand for mutual information, classical correlation, and quantum correlation, respectively. Furthermore, the critical exponent ηηη of the spin-spin correlation is shown to relate to η=ηα/2η= η^α/2η=ηα/2 with α∈{I,C,D}α\in \{ I, C, D\}α∈{I,C,D}.
We study the mutual information between two lattice blocks in terms of von Neumann entropies for one-dimensional infinite lattice systems. Quantum q-state Potts model and transverse-field spin-1/2 XY model are considered numerically by using the infinite matrix product state approach. As a system parameter varies, block-block mutual information exhibit singular behaviors that enable us to identify the critical points for the quantum phase transitions. As happens with von Neumann entanglement entropy of single block, at critical points, block-block mutual information for two adjacent blocks show a logarithmic leading behavior with increasing the size of the blocks, which yields the central charge c of the underlying conformal field theory, as it should be. It seems that two disjoint blocks show a similar logarithmic growth of the mutual information as a characteristic property of critical systems but the proportional coefficients of the logarithmic term are very different from the central charges. As the separation between the two lattice blocks increases, the mutual information reveals a consistent power-law decaying behavior for various truncation dimensions and lattice-block sizes. The critical exponent of block-block mutual information in the thermodynamic limit is estimated by extrapolating the exponents of power-law decaying regions for finite truncation dimensions. For a given lattice-block size ℓ, the critical exponents for the same universality classes seem to have very close values each other. Whereas the critical exponents have different values to a degree of distinction for the different universality classes. As the lattice-block size becomes bigger, the critical exponent becomes smaller. We discuss a relation between the exponents of block-block mutual information and correlation with the Shatten one-norm of block-block correlation.
The power of matrix product states to describe infinite-size translational-invariant critical spin chains is investigated. At criticality, the accuracy with which they describe ground-state properties of a system is limited by the size $\ensuremath{\chi}$ of the matrices that form the approximation. This limitation is quantified in terms of the scaling of the half-chain entanglement entropy. In the case of the quantum Ising model, we find $S\ensuremath{\sim}\frac{1}{6}\text{log}\text{ }\ensuremath{\chi}$ with high precision. This result can be understood as the emergence of an effective finite correlation length ${\ensuremath{\xi}}_{\ensuremath{\chi}}$ ruling all the scaling properties in the system. We produce six extra pieces of evidence for this finite-$\ensuremath{\chi}$ scaling, namely, the scaling of the correlation length, the scaling of magnetization, the shift of the critical point, the scaling of the entanglement entropy for a finite block of spins, the existence of scaling functions, and the agreement with analogous classical results. All our computations are consistent with a scaling relation of the form ${\ensuremath{\xi}}_{\ensuremath{\chi}}\ensuremath{\sim}{\ensuremath{\chi}}^{\ensuremath{\kappa}}$, with $\ensuremath{\kappa}=2$ for the Ising model. In the case of the Heisenberg model, we find similar results with the value $\ensuremath{\kappa}\ensuremath{\sim}1.37$. We also show how finite-$\ensuremath{\chi}$ scaling allows us to extract critical exponents. These results are obtained using the infinite time evolved block decimation algorithm which works in the thermodynamical limit and are verified to agree with density-matrix renormalization-group results and their classical analog obtained with the corner transfer-matrix renormalization group.
String correlations are investigated in an infinite-size XXZ spin-1 chain. By using the infinite matrix product state representation, we calculate a long-range string order. In the XY phase, the string correlations decay within a relatively very large lattice distance, which makes a finite-size study difficult to verify the non-existence of the string order. Thus, in the Haldane phase, the non-vanishing string correlations in the limit of a very large distance allow to characterize the phase boundaries to the XY phase and the Neel phase, which implies that the transverse long-range string order is the order parameter for the Haldane phase. In addition, the singular behaviors of the von Neumann entropy and the fidelity per lattice site are shown to capture clearly the phase transition points that are consistent with the results from the string order. The estimated critical points including a BKT transition from the XY phase to the Haldane phase agree well with the previous results: Δc2=0Δ_{c2} = 0Δc2​=0 for the XY-Haldane phase transition and Δc3=1.185Δ_{c3} = 1.185Δc3​=1.185 for the Haldane-N'eel phase transition from the density renormalization group. From a finite-entanglement scaling of the von Neumann entropy with respect to the truncation dimension, the central charges are found to be c≃1.0c \simeq 1.0c≃1.0 at Δc2=0Δ_{c2} = 0Δc2​=0 and c≃0.5c \simeq 0.5c≃0.5 at Δc3=1.185Δ_{c3} = 1.185Δc3​=1.185, respectively, which shows that the XY-Haldane phase transition at Δc2=0Δ_{c2} = 0Δc2​=0 belongs to the Heisenberg universality class, while the Haldane-Neel phase transition at Δc2=1.185Δ_{c2} = 1.185Δc2​=1.185 belongs to the two-dimensional classical Ising universality class. It is also shown that, the long-range order parameters and the von Neumann entropy, as well as the fidelity per site approach, can be applied to characterize quantum phase transitions as a universal phase transition indicator for one-dimensional lattice many-body systems.
The entanglement in quantum XY spin chains of arbitrary length is investigated via the geometric measure of entanglement. The emergence of entanglement is explained intuitively from the perspective of perturbations. The model is solved exactly and the energy spectrum is determined and analyzed in particular for the lowest two levels for both finite and infinite systems. The overlaps for these two levels are calculated analytically for arbitrary number of spins. The entanglement is hence obtained by maximizing over a single parameter. The corresponding ground-state entanglement surface is then determined over the entire phase diagram, and its behavior can be used to delineate the boundaries in the phase diagram. For example, the field-derivative of the entanglement becomes singular along the critical line. The form of the divergence is derived analytically and it turns out to be dictated by the universality class controlling the quantum phase transition. The behavior of the entanglement near criticality can be understood via a scaling hypothesis, analogous to that for free energies. The entanglement density vanishes along the so-called disorder line in the phase diagram, the ground space is doubly degenerate and spanned by two product states. The entanglement for the superposition of the lowest two states is also calculated. The exact value of the entanglement depends on the specific form of superposition. However, in the thermodynamic limit the entanglement density turns out to be independent of the superposition. This proves that the entanglement density is insensitive to whether the ground state is chosen to be the spontaneously Z2Z_2Z2​ symmetry broken one or not. The finite-size scaling of entanglement at critical points is also investigated from two different view points. First, the maximum in the field-derivative of the entanglement density is computed and fitted to a logarithmic dependence of the system size, thereby deducing the correlation length exponent for the Ising class using only the behavior of entanglement. Second, the entanglement density itself is shown to possess a correction term inversely proportional to the system size, with the coefficient being universal (but with different values for the ground state and the first excited state, respectively).
The characterization of the quantum phase transition in a lowdimensional system has attracted a considerable amount of attention in quantum manybody systems. As one of the fundamental models in quantum magnetism, spin-1 models have richer phase diagrams and show more complex physical phenomena. In the spin-1 antiferromagnetic XXZ model, the Haldane phase and the Nel phase are the gapped topologic phases which cannot be characterized by the local order parameters. To characterize the nature in such phases, one has to calculate the non-local long range order parameters. Normally, the non-local order parameter in the topological phase is obtained from the extrapolation of finite-sized system in numerical study. However, it is difficult to extract the critical exponents with such an extrapolated non-local order parameter due to the numerical accuracy. In a recently developed tensor network representation, i.e., the infinite matrix product state (iMPS) algorithm, it was shown that the non-local order can be directly calculated from a very large lattice distance in an infinite-sized system rather than an extrapolated order parameter in a finite-sized system. Therefore, it is worthwhile using this convenient technique to study the non-local orders in the topological phases and characterize the quantum criticalities in the topological quantum phase transitions. In this paper, by utilizing the infinite matrix product state algorithm based on the tensor network representation and infinite time evolving block decimation method, the quantum entanglement, fidelity, and critical exponents of the topological phase transition are investigated in the one-dimensional infinite spin-1 bond-alternating XXZ Heisenberg model. It is found that there is always a local dimerization order existing in the whole parameter range when the bond-alternative strength parameter changes from 0 to 1. Also, due to the effect of the bond-alternating, there appears a quantum phase transition from the long-rang ordering topological Nel phase to the local ordering dimerization phase. The von Neumann entropy, fidelity per lattice site, and order parameters all give the same phase transition point at c = 0.638. To identify the type of quantum phase transition, the central charge c 0.5 is extracted from the ground state von Neumann entropy and the finite correlation length, which indicates that the phase transition belongs to the two-dimensional Ising universality class. Furthermore, it is found that the Nel order and the susceptibility of Nel order have power-law relations to |-c|. From the numerical fitting of the Nel order and its susceptibility, we obtain the characteristic critical exponents ' = 0.236 and ' = 0.838. It indicates that such critical exponents from our method characterize the nature of the quantum phase transition. Our critical exponents from the iMPS method can provide guidance for studying the properties of the phase transition in quantum spin systems.
We investigate the short-time quantum critical dynamics in the imaginary-time relaxation processes of finite-size systems. Universal scaling behaviors exist in the imaginary-time evolution. In particular, the system undergoes a critical initial slip stage characterized by an exponent $\ensuremath{\theta}$, in which an initial power-law increase emerges in the imaginary-time correlation function when the initial state has a zero order parameter and a vanishing correlation length. Under different initial conditions, the quantum critical point and critical exponents can be determined from the universal scaling behaviors. We apply the method to the one- and two-dimensional transverse field Ising models using quantum Monte Carlo (QMC) simulations. In the one-dimensional case, we locate the quantum critical point at (h/J)c=1.00003(8){(h/J)}_{c}=1.000\phantom{\rule{0.16em}{0ex}}03(8)(h/J)c​=1.00003(8) in the thermodynamic limit, and we estimate the critical initial slip exponent $\ensuremath{\theta}=0.3734(2)$ and the static exponent $\ensuremath{\beta}/\ensuremath{\nu}=0.1251(2)$ by analyzing data on chains of length L=32−−256L=32--256L=32−−256 and 48--256, respectively. For the two-dimensional square-lattice system, the critical coupling ratio is given by 3.04451(7)3.044\phantom{\rule{0.16em}{0ex}}51(7)3.04451(7) in the thermodynamic limit, while the critical exponents are $\ensuremath{\theta}=0.209(4)$ and $\ensuremath{\beta}/\ensuremath{\nu}=0.518(1)$ estimated by data on systems of size L=24L=24L=24--64 and 32--64, respectively. Remarkably, the critical initial slip exponents obtained in both models are notably distinct from their classical counterparts due to the essential differences between classical and quantum dynamics. The short-time critical dynamics and the imaginary-time relaxation QMC approach can be readily adapted to various models.