
2
sharp domain wall, as shown in Fig. 1. We find that sim-
ilar to the classical case [22,23], as the time evolves, the
sharp domain wall will blur and expand into an inter-
facial region. Let us focus on the behavior of the local
parameter at the position x, with xdenotes the distance
to the initial domain wall. When xis far away from the
interfacial region, the local order parameter decays as
M(τ, x)∝τ−β/νz . As time elapses, the interfacial region
spreads to the position x. Accordingly, the order param-
eter changes to decay as M(τ, x)∝τ−β1/νz . Here, β1is
a purely dynamic exponent, since it has no equilibrium
counterpart, similar to the classical case [22,23]. A full
scaling form is then developed to explain this behavior.
We take the 1D and 2D quantum Ising models as ex-
amples to verify this scaling theory. From the numerical
results, we find that the scaling function is an analytical
function and β1seems not an independent critical expo-
nent. Instead, it satisfies β1/νz =β/νz + 1, in analogy
to the classical case [22,23].
τ
Spin down Spin up
Interfacial
…
FIG. 1. Sketch of the imaginary-time relaxation dynamics
from a semi-ordered initial state. A sharp domain wall sepa-
rates the spin-up and -down domains in the initial state. As
time evolves, the domain wall extends to an interfacial region.
For the imaginary-time relaxation dynamics, the evo-
lution of the wave function |ψ(τ)iobeys the imaginary-
time Schr¨odinger equation −∂τ|ψ(τ)i=H|ψ(τ)iwith
the normalization condition hψ(τ)|ψ(τ)i= 1 [27,28].
The formal solution of the Schr¨odinger equation is given
by |ψ(τ)i= exp(−τH)|ψ(τ0)i/Z, in which τ0is the ini-
tial time of the evolution and Z=kexp(−τ H)|ψ(τ0)ik,
with k · k being the modulo operation. Since both the
imaginary-time dynamics and the classical dynamics de-
scribed by Model A, which refers to a set of models with
no conservation laws [29], are purely dissipative dynam-
ics, one expects that their dynamic scaling behaviors near
the critical point are similar [7]. For instance, from the
saturated ordered initial state, the order parameter M
follows the same power-law decay with the evolution time
with the exponent β/νz for both classical model A dy-
namics and quantum imaginary-time relaxation dynam-
ics.
Here we study the influence of a different type of initial
state to the quantum imaginary-time relaxation dynam-
ics. The initial state is set as a semi-ordered state with
two fully-ordered domains with opposite spin directions,
separated by a sharp domain wall, as shown in Fig. 1.
Since the translational symmetry is broken by the ini-
tial state, it is expected that the evolution behavior of
Mdepends on the distance to the initial domain wall x.
Near the domain wall, the spin at small xfeels a stronger
spin-flip intention from the other domain where the spins
orientated in the opposite direction than its homogeneous
ordered environment. Accordingly, one anticipates that
at the critical point, for small x, the local order parame-
ter should follows
M(x, τ)∝τ−β1/νz,(1)
with β1being an additional critical exponent which is
larger than β, since the other domain lures the spin at
xto flip. Note that similar to the classical case [22,23],
here β1is a purely nonequilibrium critical exponent, since
it has no static counterpart.
In contrast, for large x, the dynamic scaling behavior of
M(x, τ) is much richer as a result of the spread of the ef-
fects induced by the domain wall, as illustrated in Fig. 1.
In the short-time stage, the domain wall region is too far
away to control the spin at xand thus the local order pa-
rameter Mdecays according to M(x, τ)∝τ−β/νz , sim-
ilar to the case with homogeneous ordered initial state.
Therefore, this stage is referred to as the ‘homogeneous
region’. As time passes by, the domain wall extends into
an ‘interfacial region’ with growing size. When the lo-
cation xis absorbed into this interfacial region, Mwill
evolve following Eq. (1).
Scaling theory. For the relaxation critical dynamics
with homogeneous initial state, a characteristic scaling
behavior is the appearance of the critical initial slip char-
acterized by an independent critical exponent θ. A na-
ture question is whether or not β1is another independent
exponent. To answer this question, we develop a full scal-
ing form to describe the whole imaginary-time relaxation
process with the semi-ordered initial state. In analogy
with the classical situation [22,23], the scaling form of
the local order parameter Mat a quantum critical point
is given by
M(τ, x) = τ−β/νz f(xτ−1/z ),(2)
in which f(xτ−1/z) is the scaling function. By compar-
ing Eq. (2) with the imaginary-time relaxation scaling
theory with the homogeneous initial state [7], one finds
that the initial order parameter is absent in Eq. (2). The
reason is that the initial state keeps invariant under the
renormalization transformation, as illustrated in Fig. 2.
Similarly, the initial correlation is also not included since
both the initial correlation length and correlation time
are zero.
It is expected that the scaling behaviors discussed
above should be covered by the full scaling form Eq. (2).
This gives some constraints on the scaling function
f(xτ−1/z ): (i) f(xτ −1/z) should be an odd function of
xsince Mshould change sign on switching the spin