Time-dependent Order Parameter Research Articles

Random anisotropy magnet at finite temperature

We present finite-temperature Monte Carlo studies of a 2D random-anisotropy (RA) magnet on lattices containing one million spins. The correlated spin-glass state predicted by analytical theories is reproduced in simulations, as are the field-cooled and zero-field-cooled magnetization curves observed in experiments. The orientations of lattice spins begin to freeze when the temperature is lowered. The freezing transition is due to the energy barriers generated by the RA rather than due to random interactions in conventional spin-glasses. We describe freezing by introducing the time-dependent spin-glass order parameter q and the spin-melting time τ M defined via q = τ M/t above freezing, where t is the time of the experiment represented by the number of Monte Carlo steps.

Real-time simulation of light-driven spin chains on quantum computers

In this work, we study the real-time evolution of periodically driven (Floquet) systems on a quantum computer using IBM quantum devices. We consider a driven Landau-Zener model and compute the transition probability between the Floquet steady states as a function of time. We find that for this simple one-qubit model, Floquet states can develop in real time, as indicated by the transition probability between Floquet states. Next, we model light-driven spin chains and compute the time-dependent antiferromagnetic order parameter. We consider models arising from light coupling to the underlying electrons as well as those arising from light coupling to phonons. For the two-spin chains, the quantum devices yield time evolutions that match the effective Floquet Hamiltonian evolution for both models once readout error mitigation is included. For three-spin chains, zero-noise extrapolation yields a time dependence that follows the effective Floquet time evolution. Therefore, the current IBM quantum devices can provide information on the dynamics of small Floquet systems arising from light drives once error mitigation procedures are implemented.

Cumulants of the chiral order parameter in a nonequilibrium Bjorken expansion with a QCD critical point

To understand experimentally obtained net-proton number cumulants in the search for the QCD critical point, we study a dynamical model based on an effective quark-meson Lagrangian with chiral symmetry. We investigate the evolution of the expanding medium created in a heavy-ion collision using a spatially homogeneous fluid and a time-dependent order parameter, the sigma field evolved by a Langevin equation. We extract cumulants of the sigma field along a parametrized freeze-out curve and match the obtained freeze-out points to corresponding beam energies. These cumulants can be related to cumulants of the net-proton number through the sigma-proton coupling to provide a qualitative comparison to experimental data from STAR’s beam energy scan program. We demonstrate that the presence of the spinodal or mixed phase region around the first-order chiral phase transition allows for a wide interval of cumulants at the lowest beam energies.

Dynamical quantum phase transition in a bosonic system with long-range interactions

In this paper, we investigate the dynamical quantum phase transitions appearing in the Loschmidt echo and the time-dependent order parameter of a quantum system of harmonically coupled degenerate bosons as a function of the power-law decay $\sigma$ of long-range interactions. Following a sudden quench, the nonequilibrium dynamics of this system are governed by a set of nonlinear coupled Ermakov equations. To solve them, we develop an analytical approximation valid at late times. Based on this approximation, we show that the emergence of a dynamical quantum phase transition hinges on the generation of a finite mass gap following the quench, starting from a massless initial state. In general, we can define two distinct dynamical phases characterized by the finiteness of the post-quench mass gap. The Loschmidt echo exhibits periodical nonanalytic cusps whenever the initial state has a vanishing mass gap and the final state has a finite mass gap. These cusps are shown to coincide with the maxima of the time-dependent long-range correlations.

Fate of dynamical phases of a BCS superconductor beyond the dissipationless regime

The BCS model of an isolated superconductor initially prepared in a nonequilibrium state, predicts the existence of interesting dynamical phenomena in the time-dependent order parameter as decaying oscillations, persistent oscillations and overdamped dynamics. To make contact with real systems remains an open challenge as one needs to introduce dissipation due to the environment in a self-consistent computation. Here, we reach this goal with the use of the Keldysh formalism to treat the effect of a thermal bath. We show that, contrary to the dissipationless case, all dynamical phases reach the equilibrium order parameter in a characteristic time that depends on the coupling with the bath. Remarkably, as time evolves, the overdamped phase shows a fast crossover where the superconducting order parameter recovers to reach a state with a well-developed long range order that tends towards equilibrium with the damped Higgs mode oscillations. Our results provide a benchmark for the description of the dynamics of real out-of-equilibrium superconductors relevant for quantum technological applications.

Consequences of integrability breaking in quench dynamics of pairing Hamiltonians

We study the collisionless dynamics of two classes of nonintegrable pairing models. One is a BCS model with separable energy-dependent interactions, the other - a 2D topological superconductor with spin-orbit coupling and a band-splitting external field. The long-time quantum quench dynamics at integrable points of these models are well understood. Namely, the squared magnitude of the time-dependent order parameter $\Delta(t)$ can either vanish (Phase I), reach a nonzero constant (Phase II), or periodically oscillate as an elliptic function (Phase III). We demonstrate that nonintegrable models too exhibit some or all of these nonequilibrium phases. Remarkably, elliptic periodic oscillations persist, even though both their amplitude and functional form change drastically with integrability breaking. Striking new phenomena accompany loss of integrability. First, an extremely long time scale emerges in the relaxation to Phase III, such that short-time numerical simulations risk erroneously classifying the asymptotic state. This time scale diverges near integrable points. Second, an entirely new Phase IV of quasiperiodic oscillations of $|\Delta|$ emerges in the quantum quench phase diagrams of nonintegrable pairing models. As integrability techniques do not apply for the models we study, we develop the concept of asymptotic self-consistency and a linear stability analysis of the asymptotic phases. With the help of these new tools, we determine the phase boundaries, characterize the asymptotic state, and clarify the physical meaning of the quantum quench phase diagrams of BCS superconductors. We also propose an explanation of these diagrams in terms of bifurcation theory.

Emergence and destruction of macroscopic wave functions

The concept of the macroscopic wave function is a key for understanding macroscopic quantum phenomena. The existence of this object reflects a certain order, as it is present in a Bose-Einstein condensate when a single-particle orbital is occupied by a macroscopic number of bosons. We extend these ideas to situations in which a condensate is acted on by an explicitly time-dependent force. While one might assume that such a force would necessarily degrade any pre-existing order, we demonstrate that macroscopic wave functions can persist even under strong forcing. Our definition of the time-dependent order parameter is based on a comparison of the evolution of N-particle states on the one hand, and of states with particles on the other. Our simulations predict the possibility of an almost instantaneous dynamical destruction of a macroscopic wave function under currently accessible experimental conditions.

Localization transition, Lifschitz tails, and rare-region effects in network models.

Effects of heterogeneity in the suspected-infected-susceptible model on networks are investigated using quenched mean-field theory. The emergence of localization is described by the distributions of the inverse participation ratio and compared with the rare-region effects appearing in simulations and in the Lifschitz tails. The latter, in the linear approximation, is related to the spectral density of the Laplacian matrix and to the time dependent order parameter. I show that these approximations indicate correctly Griffiths phases both on regular one-dimensional lattices and on small-world networks exhibiting purely topological disorder. I discuss the localization transition that occurs on scale-free networks at γ=3 degree exponent.

On the broken time translation symmetry in macroscopic systems: Precessing states and off-diagonal long-range order

The broken symmetry state with off-diagonal long-range order (ODLRO), which is characterized by the vacuum expectation value of the operator of creation of the conserved quantum number Q, has the time-dependent order parameter. However, the breaking of the time translation symmetry is observable only if the charge Q is not strictly conserved and may decay. This dichotomy is resolved in systems with quasi-ODLRO. These systems have two well separated relaxation times: the relaxation time τ Q of the charge Q and the energy relaxation time τ E . If τ Q ≫ τ E , the perturbed system relaxes first to the state with the ODLRO, which persists for a long time and finally relaxes to the full equilibrium static state. In the limit τQ → ∞, but not in the strict limit case when the charge Q is conserved, the intermediate ODLRO state can be considered as the ground state of the system at fixed Q with the observable spontaneously broken time translation symmetry. Examples of systems with quasi-ODLRO are provided by superfluid phase of liquid 4He, Bose-Einstein condensation of magnons (phase coherent spin precession) and precessing vortices.

Approximate solution to the stochastic Kuramoto model

We study Kuramoto phase oscillators with temporal fluctuations in the frequencies. The infinite-dimensional system can be reduced in a Gaussian approximation to two first-order differential equations. This yields a solution for the time-dependent order parameter, which characterizes the synchronization between the oscillators. The known critical coupling strength is exactly recovered by the Gaussian theory. Extensive numerical experiments further show that the analytical results are very accurate below and sufficiently above the critical value. We obtain the asymptotic order parameter in closed form, which suggests a tighter upper bound for the corresponding scaling. As a last point, we elaborate the Gaussian approximation in complex networks with distributed degrees.

Reduced dynamics for delayed systems with harmonic or stochastic forcing

The analysis of nonlinear delay-differential equations (DDEs) subjected to external forcing is difficult due to the infinite dimensionality of the space in which they evolve. To simplify the analysis of such systems, the present work develops a non-homogeneous center manifold (CM) reduction scheme, which allows the derivation of a time-dependent order parameter equation in finite dimension. This differential equation captures the major dynamical features of the delayed system. The forcing is assumed to be small compared to the amplitude of the autonomous system, in order to cause only small variations of the fixed points and of the autonomous CM. The time-dependent CM is shown to satisfy a non-homogeneous partial differential equation. We first briefly review CM theory for DDEs. Then we show, for the general scalar case, how an ansatz that separates the CM into one for the autonomous problem plus an additional time-dependent order-two correction leads to satisfying results. The paper then details the application to a transcritical bifurcation subjected to single or multiple periodic forcings. The validity limits of the reduction scheme are also highlighted. Finally, we characterize the specific case of additive stochastic driving of the transcritical bifurcation, where additive white noise shifts the mode of the probability density function of the state variable to larger amplitudes.

Unifying treatment of nonequilibrium and unstable dynamics of cold bosonic atom system with time-dependent order parameter in Thermo Field Dynamics

The coupled equations which describe the temporal evolution of the Bose–Einstein condensed system are derived in the framework of nonequilibrium Thermo Field Dynamics. The key element is that they are not the naive assemblages of assumed equations, but are the self-consistent ones derived by appropriate renormalization conditions. While the order parameter is time-dependent, an explicit quasiparticle picture is constructed by a time-dependent expansion. Our formulation is valid even for the system with a unstable condensate, and describes the condensate decay caused by the Landau instability as well as by the dynamical one.

Simulation of non-equilibrium phase transition in a random magnetic system

The time-dependent phase transition of a random magnetic system is investigated using a dynamic Monte Carlo simulation. The present random system is a simple dilution of magnetic atoms with Ising spins in a square lattice. The time-dependent order parameters are calculated as a function of temperature for several concentrations of magnetic atoms. The results indicate that as the averaging time becomes shorter, the phase transition becomes more gradual. This tendency is most pronounced around the percolation concentration. These results are compared to results of a recent NMR spin echo experiment on MnxCd1-x(HCOO)2·2(NH2)2CO reported by Kubo et al. (2007). Although the present simulation includes essentially no adjustable parameter, the results are in satisfactory agreement.

Dynamics of texture for spinor condensates

We study a dynamics of texture for a two-component spinor bose condensate. This is carried out by adopting a time dependent Landau-Ginzburg Lagrangian for a spinor order parameter. By using a polar form of the spinor order parameter, we obtain a field equation for the texture. In particular we consider a one dimensional model in which we can obtain analytic forms for the textures in terms of elliptic functions of several kinds. We find that these solutions are characterized by a modulus parameter, and changes in this parameter cause structural changes of texture.

Generalization of Faraday’s Law to Include Nonconservative Spin Forces

The usual Faraday's Law E=-dPhi/dt determines an electromotive force E which accounts only for forces resulting from the charge of electrons. In ferromagnetic materials, in general, there exist nonconservative spin forces which also contribute to E. These might be included in Faraday's Law if the magnetic flux Phi is replaced by [Planck's constant/(-e)]gamma, where gamma is a Berry phase suitably averaged over the electron spin direction. These contributions to E represent the requirements of energy conservation in itinerant ferromagnets with time dependent order parameters.

Backbone Motions of Free and Pheromone-Bound Major Urinary Protein I Studied by Molecular Dynamics Simulation

Molecular motions of free and pheromone-bound mouse major urinary protein I, previously investigated by NMR relaxation, were simulated in 30 ns molecular dynamics (MD) runs. The backbone flexibility was described in terms of order parameters and correlation times, commonly used in the NMR relaxation analysis. Special attention was paid to the effect of conformational changes on the nanosecond time scale. Time-dependent order parameters were determined in order to separate motions occurring on different time scales. As an alternative approach, slow conformational changes were identified from the backbone torsion angle variances, and "conformationally filtered" order parameters were calculated for well-defined conformation states. A comparison of the data obtained for the free and pheromone-bound protein showed that some residues are more rigid in the bound form, but a larger portion of the protein becomes more flexible upon the pheromone binding. This finding is in general agreement with the NMR results. The higher flexibility observed on the fast (fs-ps) time scale was typically observed for the residues exhibiting higher conformational freedom on the ns time scale. An inspection of the hydrogen bond network provided a structural explanation for the flexibility differences between the free and pheromone-bound proteins in the simulations.

Non-linear field dynamics of reacting nuclear, atomic and molecular species

We develop a field theoretical formalism for describing the non-linear dynamics of interacting nuclear, atomic and molecular species. Specifically, four simple types of reaction are considered all of which involve up to two reactants and two products. The equations of motion for space- and time-dependent quantum order parameter fields, describing the distribution of reactant and product species, are derived from Heisenberg's equation using appropriate second-quantised Hamiltonians. For simplicity, both internal degrees of freedom and heat bath modes are being neglected at this stage. We solve a number of cases exactly both in the space-independent and space-dependent régimes. The results obtained possess a variety of interesting properties which include solitonic behaviour and propagating waves of various types indicating the existence of a wealth of non-linear phenomena.

Molecular dynamics study of structural changes in berlinite

Fractional coordinates and anisotropic temperature factors of atoms in berlinite, AlPO4 with the quartz topology, were successfully simulated in a molecular dynamics simulation (MDS) at high temperatures. Time-dependent and time-averaged atomic order parameters were analyzed in detail with the aid of spectral densities calculated from trajectory data. These parameters show characteristic behavior of the order-disorder regime for a structure change, where atoms spend most of the time oscillating around the 1-sites (or 2-sites) in the low temperature α-phase, but oscillate over both sites in the higher temperature α-phase and the β-phase. In the spectral density functions calculated for atom order parameters, a nearly zero-frequency excitation, which is accompanied by the emergence of large-scale 1 and 2 clusters, appears at the Γ point of the Brillouin zone below the transition point T o, and increases in intensity up to T o. A low-lying optic branch along Γ-M, which is strongly temperature dependent in the small q-region, is another characteristic of the spectral density functions for the β phase. The spectrum at Γ continuously reduces its frequency from 0.6 THz at temperatures far above T o to nearly 0 THz at temperatures approaching T o from above. The dynamical behavior in β berlinite rapidly but continuously changes from that in less oscillatory clusters in the vicinity of T o to that in the typical β phase at temperatures departing from T o.

Time Delay in the Kuramoto Model of Coupled Oscillators

We generalize the Kuramoto model of coupled oscillators to allow time-delayed interactions. New phenomena include bistability between synchronized and incoherent states, and unsteady solutions with time-dependent order parameters. We derive exact formulas for the stability boundaries of the incoherent and synchronized states, as a function of the delay, in the special case where the oscillators are identical. The experimental implications of the model are discussed for populations of chirping crickets, where the finite speed of sound causes communication delays, and for physical systems such as coupled phase-locked loops or lasers.

Slowing down in the three-dimensional three-state Potts glass with nearest neighbor exchange : A Monte Carlo study

,Static and dynamic properties of the Potts model on the simple cubic lattice with nearest neighbor ±Ĵ-interaction are obtained from Monte Carlo simulations in a temperature range where full thermal equilibrium still can be achieved (T/Ĵ ≥ 0.6). For a lattice size L = 16, in this range finite size effects are still negligible, but the data for the spin glass susceptibility agree with previous extrapolations based on finite size scaling of very small lattices. While the static properties are compatible with a zero temperature transition, they certainly do not prove it. Unlike the Ising spin glass, the decay of the time-dependent order parameter is compatible with a simple Kohlrausch function, q(t) ∝ exp[‒(t/τ) y(T) ], while a power law prefactor cannot be distinguished. The Kohlrausch exponent y(T) decreases from y ≈ 0.64 at T/Ĵ = 1.1 to y ≈ 0.37 at T/Ĵ = 0.6 however. The relaxation time τ is compatible with the exponential divergence postulated by McMillan for spin glasses at their lower critical dimension, lnτ ∝ T‒ σ but the exponent σ that can be extracted (σ ≈ 2.5) still differs significantly from the theoretical value, σ = 3. Thus the present results support the conclusion that the Potts spin glass in d = 3 dimensions differs qualitatively from the Ising spin glass.