Exact and Truncated Dynamics in Nonequilibrium Field Theory

Nonperturbative dynamics of quantum fields out of equilibrium is often described by the time evolution of a hierarchy of correlation functions, using approximation methods such as Hartree, large N, and nPI-effective action techniques. These truncation schemes can be implemented equally well in a classical statistical system, where results can be tested by comparison with the complete nonlinear evolution obtained by numerical methods. For a 1+1 dimensional scalar field we find that the early-time behaviour is reproduced qualitatively by the Hartree dynamics. The inclusion of direct scattering improves this to the quantitative level. We show that the emergence of nonthermal temperature profiles at intermediate times can be understood in terms of the fixed points of the evolution equations in the Hartree approximation. The form of the profile depends explicitly on the initial ensemble. While the truncated evolution equations do not seem to be able to get away from the fixed point, the full nonlinear evolution shows thermalization with a (surprisingly) slow relaxation.Comment: 30 pages with 12 eps figures, minor changes; to appear in Phys.Rev.

The approach to thermalization in the classical phi^4 theory in 1+1 dimensions: energy cascades and universal scaling

We study the dynamics of thermalization and the approach to equilibrium in the classical phi^4 theory in 1+1 spacetime dimensions. At thermal equilibrium we exploit the equivalence between the classical canonical averages and transfer matrix quantum traces of the anharmonic oscillator to obtain exact results for the temperature dependence of several observables, which provide a set of criteria for thermalization. We find that the Hartree approximation is remarkably accurate in equilibrium. The non-equilibrium dynamics is studied by numerically solving the equations of motion in light-cone coordinates for a broad range of initial conditions and energy densities.The time evolution is described by several stages with a cascade of energy towards the ultraviolet. After a transient stage, the spatio-temporal gradient terms become larger than the nonlinear term and a stage of universal cascade emerges.This cascade starts at a time scale t_0 independent of the initial conditions (except for very low energy density). Here the power spectra feature universal scaling behavior and the front of the cascade k(t) grows as a power law k(t) sim t^alpha with alpha lesssim 0.25. The wake behind the cascade is described as a state of Local Thermodynamic Equilibrium (LTE) with all correlations being determined by the equilibrium functional form with an effective time dependent temperatureTeff(t) which slowly decreases as sim t^{-alpha}.Two well separated time scales emerge while Teff(t) varies slowly, the wavectors in the wake with k < k(t) attain LTE on much shorter time scales.This universal scaling stage ends when the front of the cascade reaches the cutoff at a time t_1 sim a^{-1/alpha}. Virialization starts to set much earlier than LTE. We find that strict thermalization is achieved only for an infinite time scale.Comment: relevance for quantum field theory discussed providing validity criteria. To appear in Phys. Rev.

Nonlinear dynamics of phase separation in thin films

We present a long-wavelength approximation to the Navier-Stokes Cahn-Hilliard equations to describe phase separation in thin films. The equations we derive underscore the coupled behaviour of free-surface variations and phase separation. We introduce a repulsive substrate-film interaction potential and analyse the resulting fourth-order equations by constructing a Lyapunov functional, which, combined with the regularizing repulsive potential, gives rise to a positive lower bound for the free-surface height. The value of this lower bound depends on the parameters of the problem, a result which we compare with numerical simulations. While the theoretical lower bound is an obstacle to the rupture of a film that initially is everywhere of finite height, it is not sufficiently sharp to represent accurately the parametric dependence of the observed dips or `valleys' in free-surface height. We observe these valleys across zones where the concentration of the binary mixture changes sharply, indicating the formation of bubbles. Finally, we carry out numerical simulations without the repulsive interaction, and find that the film ruptures in finite time, while the gradient of the Cahn--Hilliard concentration develops a singularity.Comment: 26 pages, 20 figures, PDFLaTeX with RevTeX4 macros. A thorough analysis of the equations is presented in arXiv:0805.103

Waiting and Residence Times of Brownian Interface Fluctuations

We report on the residence times of capillary waves above a given height $h$ and on the typical waiting time in between such fluctuations. The measurements were made on phase separated colloid-polymer systems by laser scanning confocal microscopy. Due to the Brownian character of the process, the stochastics vary with the chosen measurement interval $\Delta t$. In experiments, the discrete scanning times are a practical cutoff and we are able to measure the waiting time as a function of this cutoff. The measurement interval dependence of the observed waiting and residence times turns out to be solely determined by the time dependent height-height correlation function $g(t)$. We find excellent agreement with the theory presented here along with the experiments.Comment: 5 figure

Bubbles and Filaments: Stirring a Cahn-Hilliard Fluid

The advective Cahn-Hilliard equation describes the competing processes of stirring and separation in a two-phase fluid. Intuition suggests that bubbles will form on a certain scale, and previous studies of Cahn-Hilliard dynamics seem to suggest the presence of one dominant length scale. However, the Cahn-Hilliard phase-separation mechanism contains a hyperdiffusion term and we show that, by stirring the mixture at a sufficiently large amplitude, we excite the diffusion and overwhelm the segregation to create a homogeneous liquid. At intermediate amplitudes we see regions of bubbles coexisting with regions of hyperdiffusive filaments. Thus, the problem possesses two dominant length scales, associated with the bubbles and filaments. For simplicity, we use use a chaotic flow that mimics turbulent stirring at large Prandtl number. We compare our results with the case of variable mobility, in which growth of bubble size is dominated by interfacial rather than bulk effects, and find qualitatively similar results.Comment: 20 pages, 27 figures. RevTeX

Dissipation in equations of motion of scalar fields

The methods of non-equilibrium quantum field theory are used to investigate the possibility of representing dissipation in the equation of motion for the expectation value of a scalar field by a friction term, such as is commonly included in phenomenological inflaton equations of motion. A sequence of approximations is exhibited which reduces the non-equilibrium theory to a set of local evolution equations. However, the adiabatic solution to these evolution equations which is needed to obtain a local equation of motion for the expectation value is not well defined; nor, therefore, is the friction coefficient. Thus, a non-equilibrium treatment is essential, even for a system that remains close to thermal equilibrium, and the formalism developed here provides one means of achieving this numerically.Comment: 17 pages, 5 figure

Particle Production and Effective Thermalization in Inhomogeneous Mean Field Theory

As a toy model for dynamics in nonequilibrium quantum field theory we consider the abelian Higgs model in 1+1 dimensions with fermions. In the approximate dynamical equations, inhomogeneous classical (mean) Bose fields are coupled to quantized fermion fields, which are treated with a mode function expansion. The effective equations of motion imply e.g. Coulomb scattering, due to the inhomogeneous gauge field. The equations are solved numerically. We define time dependent fermion particle numbers with the help of the single-time Wigner function and study particle production starting from inhomogeneous initial conditions. The particle numbers are compared with the Fermi-Dirac distribution parametrized by a time dependent temperature and chemical potential. We find that the fermions approximately thermalize locally in time.Comment: 16 pages + 6 eps figures, some clarifications and two references added, typos corrected; to appear in Phys.Rev.

Symmetry Breaking and False Vacuum Decay after Hybrid Inflation

We discuss the onset of symmetry breaking from the false vacuum in generic scenarios in which the mass squared of the symmetry breaking (Higgs) field depends linearly with time, as it occurs, via the evolution of the inflaton, in models of hybrid inflation. We show that the Higgs fluctuations evolve from quantum to classical during the initial stages. This justifies the subsequent use of real-time lattice simulations to describe the fully non-perturbative and non-linear process of symmetry breaking. The early distribution of the Higgs field is that of a smooth classical gaussian random field, and consists of lumps whose shape and distribution is well understood analytically. The lumps grow with time and develop into ``bubbles'' which eventually collide among themselves, thus populating the high momentum modes, in their way towards thermalization at the true vacuum. With the help of some approximations we are able to provide a quasi-analytic understanding of this process.Comment: 33 pages, 16 figures, LaTeX, uses revtex. Version to be published in Phys. Rev. with minor change

Kinetic pathways of the Nematic-Isotropic phase transition as studied by confocal microscopy on rod-like viruses

We investigate the kinetics of phase separation for a mixture of rodlike viruses (fd) and polymer (dextran), which effectively constitutes a system of attractive rods. This dispersion is quenched from a flow-induced fully nematic state into the region where the nematic and the isotropic phase coexist. We show experimental evidence that the kinetic pathway depends on the overall concentration. When the quench is made at high concentrations, the system is meta-stable and we observe typical nucleation-and-growth. For quenches at low concentration the system is unstable and the system undergoes a spinodal decomposition. At intermediate concentrations we see the transition between both demixing processes, where we locate the spinodal point.Comment: 11 pages, 6 figures, accepted in J. Phys.: Condens. Matter as symposium paper for the 6th Liquid Matter Conference in Utrech

The Development of Equilibrium After Preheating

We present a fully nonlinear study of the development of equilibrium after preheating. Preheating is the exponentially rapid transfer of energy from the nearly homogeneous inflaton field to fluctuations of other fields and/or the inflaton itself. This rapid transfer leaves these fields in a highly nonthermal state with energy concentrated in infrared modes. We have performed lattice simulations of the evolution of interacting scalar fields during and after preheating for a variety of inflationary models. We have formulated a set of generic rules that govern the thermalization process in all of these models. Notably, we see that once one of the fields is amplified through parametric resonance or other mechanisms it rapidly excites other coupled fields to exponentially large occupation numbers. These fields quickly acquire nearly thermal spectra in the infrared, which gradually propagates into higher momenta. Prior to the formation of total equilibrium, the excited fields group into subsets with almost identical characteristics (e.g. group effective temperature). The way fields form into these groups and the properties of the groups depend on the couplings between them. We also studied the onset of chaos after preheating by calculating the Lyapunov exponent of the scalar fields.Comment: 15 pages, 23 figure