Scaling behavior of two-time correlations in a twisted nematic liquid crystal

Phase-ordering kinetics in the Allen-Cahn (Model A) class: Universal aspects elucidated by electrically induced transition in liquid crystals

The two-dimensional (2D) Ising model is the statistical physics textbook example for phase transitions and their kinetics. Quenched through the Curie point with Glauber rates, the late-time description of the ferromagnetic domain coarsening finds its place at the scalar sector of the Allen-Cahn (or Model A) class, which encompasses phase-ordering kinetics endowed with a nonconserved order parameter. Resisting exact results sought for theoreticians since Lifshitz's first account in 1962, the central quantities of 2D Model A-most scaling exponents and correlation functions-remain known up to approximate theories whose disparate outcomes urge experimental assessment. Here we perform such assessment based on a comprehensive study of the coarsening of 2D twisted nematic liquid crystals whose kinetics is induced by a superfast electrical switching from a spatiotemporally chaotic (disordered) state to a two-phase concurrent, equilibrium one. Tracking the dynamics via optical microscopy, we first show the sharp evidence of well-established Model A aspects, such as the dynamic exponent z=2 and the dynamic scaling hypothesis, to then move forward. We confirm the Bray-Humayun theory for Porod's regime describing intradomain length scales of the two-point spatial correlators and show that their nontrivial decay beyond the Porod's scale can be captured in a free-from-parameter fashion by Gaussian theories, namely the Ohta-Jasnow-Kawasaki (OJK) and Mazenko theories. Regarding time-related statistics, we corroborate the aging hypothesis in Model A systems, which includes the collapse of two-time correlators into a master curve whose format is, actually, best accounted for by a solution of the local scaling invariance theory: the same solution that fits the 2D nonconserved Ising model correlator along with the Fisher-Huse conjecture. We also suggest the true value for the local persistence exponent in Model A class, in disfavor of the exact outcome for the diffusion and OJK equations. Finally, we observe a fractal morphology for persistence clusters and extract their universal dimension. Given its accuracy and possibilities, this experimental setup may work as a prototype to address further universality issues in the realm of nonequilibrium systems.

Universal dynamics of coarsening during polymer-polymer thin-film spinodal dewetting kinetics

The dewetting dynamics of a supported bilayer polymer thin film on a solid substrate is investigated using grazing incidence x-ray photon correlation spectroscopy. We find that the top layer dewets via the spinodal mechanism. The kinetics of the dewetting is studied by monitoring the time evolution of the surface diffuse x-ray scattering intensity. We study the time evolution of fluctuations about the average surface structure by measuring the two-time x-ray intensity fluctuation correlation functions. Using these two-time correlation functions we quantify the crossover from early-time diffusive dynamics to hydrodynamics. The early diffusive regime satisfies dynamic universality. The two-time correlation functions also quantify the onset of hydrodynamic effects. The hydrodynamic regime is observed during the spinodal dewetting process as these interactions are not screened.

Direct test of the Gaussian auxiliary field ansatz in nonconserved order parameter phase ordering dynamics

The assumption that the local order parameter is related to an underlying spatially smooth auxiliary field, u(r[over ⃗],t), is a common feature in theoretical approaches to non-conserved order parameter phase separation dynamics. In particular, the ansatz that u(r[over ⃗],t) is a Gaussian random field leads to predictions for the decay of the autocorrelation function which are consistent with observations, but distinct from predictions using alternative theoretical approaches. In this paper, the auxiliary field is obtained directly from simulations of the time-dependent Ginzburg-Landau equation in two and three dimensions. The results show that u(r[over ⃗],t) is equivalent to the distance to the nearest interface. In two dimensions, the probability distribution, P(u), is well approximated as Gaussian except for small values of u/L(t), where L(t) is the characteristic length-scale of the patterns. The behavior of P(u) in three dimensions is more complicated; the non-Gaussian region for small u/L(t) is much larger than that in two dimensions but the tails of P(u) begin to approach a Gaussian form at intermediate times. However, at later times, the tails of the probability distribution appear to decay faster than a Gaussian distribution.

Dislocations and metastable chevrons in the electroconvective inplane normal roll state of a bent core nematic liquid crystal

We report experimental results on the formation, dynamics, and annihilation of edge dislocations of opposite topological charge in the electroconvective inplane vortex state of a bent core nematic liquid crystal. The approach of paired, oppositely charged defects toward each other is a two-step process. Near constant velocity at large separation and accelerated motion close to annihilation are found, as in the case of nematic rolls belonging to standard electroconvection. Periodic arrays of dislocations of alternating polarity form upon a sudden, strong elevation of the control parameter. Chevron structures that appear between undulatory defect chains are metastable, and their decay with time is accompanied by an exponential reduction in the dislocation density. The initial periodicity of defect chains also drops exponentially with increasing field.

Monte Carlo study of the phase transition in the critical behavior of the Ising model with shear

The critical behavior of the Ising model with nonconserved dynamics and an external driving field mimicking a shear profile is analyzed by studying its dynamical evolution in the short-time regime. Starting from high-temperature disordered configurations (fully disordered configurations, FDC), the critical temperature Tc is determined when the order parameter, defined as the absolute value of the transversal spin profile, exhibits a power-law behavior with an exponent that is a combination of some of the critical exponents of the transition. For each value of the shear field magnitude, labeled as gamma, Tc has been estimated and two stages have been found: (1) a growing stage at low values of gamma, where Tc approximately gammapsi and psi=0.52(3), and (2) a saturation regime at large gamma. The same values of Tc(gamma) were found studying the dynamical evolution from the ground-state configuration with all spins pointing in the same direction. By combining the exponents of the corresponding power laws obtained from each initial configuration, the set of critical exponents was calculated. These values, at large external field magnitude, define a critical behavior different from that of the Ising model and of other driven lattice gases.

Influence of thermal fluctuations on the geometry of interfaces of the quenched Ising model

We study the role of the quench temperature Tf in the phase-ordering kinetics of the Ising model with single spin flip in d=2,3 . Equilibrium interfaces are flat at Tf=0 , whereas at Tf>0 they are curved and rough (above the roughening temperature in d=3 ). We show, by means of scaling arguments and numerical simulations, that this geometrical difference is important for the phase-ordering kinetics as well. In particular, while the growth exponent z=2 of the size of domains L(t) approximately t 1/z is unaffected by Tf, other exponents related to the interface geometry take different values at Tf=0 or Tf>0 . For Tf>0 a crossover phenomenon is observed from an early stage where interfaces are still flat and the system behaves as at Tf=0 , to the asymptotic regime with curved interfaces characteristic of Tf>0 . Furthermore, it is shown that the roughening length, although subdominant with respect to L(t) , produces appreciable correction to scaling up to very long times in d=2 .

Dynamics of the Bose-Einstein condensation: analogy with the collapse dynamics of a classical self-gravitating Brownian gas

We consider the dynamics of a gas of free bosons within a semiclassical Fokker-Planck equation for which we give a physical justification. In this context, we find a striking similarity between the Bose-Einstein condensation in the canonical ensemble, and the gravitational collapse of a gas of classical self-gravitating Brownian particles. The paper is mainly devoted to the complete study of the Bose-Einstein "collapse" within this model. We find that at the Bose-Einstein condensation temperature Tc, the chemical potential mu(t) vanishes exponentially with a universal rate that we compute exactly. Below Tc, we show analytically that square root mu(t) vanishes linearly in a finite time t coll. After t coll, the mass of the condensate grows linearly with time and saturates exponentially to its equilibrium value for large time. We also give analytical results for the density scaling functions, for the corrections to scaling, and for the exponential relaxation time. Finally, we find that the equilibration time (above Tc) and the collapse time T coll(below Tc) both behave like -T -3 c ln|T-Tc|, near Tc.

Effect of shear on persistence in coarsening systems

We analytically study the effect of a uniform shear flow on the persistence properties of coarsening systems. The study is carried out within the anisotropic Ohta-Jasnow-Kawasaki (OJK) approximation for a system with nonconserved scalar order parameter. We find that the persistence exponent theta has a nontrivial value: theta = 0.5034 in space dimension d = 3, and theta = 0.2406 for d = 2, the latter being exactly twice the value found for the unsheared system in d = 1. We also find that the autocorrelation exponent lambda is affected by shear in d = 3 but not in d = 2.

Effect of the Elastic Constant Anisotropy on Disclination Interaction in the Nematic Polymers

In this work, disclination interaction behavior in relation to Frank elastic constant anisotropy in nematics has been studied. A large number of (+(1/2), -(1/2)) disclination pairs are revealed by spontaneous band texture in a semiflexible copolyester. The pairs show no preferential relative orientation, with the intervening fields showing intermediate patterns. A two-dimensional tensor lattice model considering unequal elastic constants is applied to simulate the interaction behavior and patterns of disclination pairs in the presence of elastic anisotropy. Scaling laws for disclination density rho(t) as a function of time step t with different elastic anisotropy are obtained as t(-nu). The value of the exponent nu decreases as elastic anisotropy is increased. Obviously, elastic anisotropy slows the texture coarsening. The simulations also show that angular forces arise in the presence of elastic anisotropy and change the patterns of pairs during the texture coarsening. When disclination density is considerably decreased, some +(1/2) disclinations start to rotate to the energetically favored patterns depending on the sign of the elastic anisotropy. As a result of the disclination rotation, the distribution of patterns of pairs continues to change during the annihilation. However, disclination pairs are influenced not only by elastic anisotropy but also by disclination interaction during the whole annihilation. Therefore, in a real system, the dependence of pairs on elastic anisotropy is not as strong as the theoretical prediction for an isolated pair, and the full pattern range of disclination pairs can be observed.

Monte Carlo results for the Ising model with shear

We study the kinetics of domain growth in the Ising model with nonconserved dynamics under the action of a stochastic driving field that mimics the action of a shear flow. At late times, we found multistriped configurations with constant transversal size and linear growth in the direction of the flow. In cases with weak shear, a regime characterized by the decreasing of the transversal size is found that could correspond to previous theoretical investigations. This behavior is confirmed by the analysis of the structure factor patterns.

Autocorrelation exponent of conserved spin systems in the scaling regime following a critical quench

We study the autocorrelation function of a conserved spin system following a quench at the critical temperature. Defining the correlation length L(t) approximately t(1/z), we find that for times t' and t satisfying L(t')<<L(t)<<L(t')(phi) well inside the scaling regime, the spin autocorrelation function behaves like s(t)s(t') approximately L(t')(-(d-2+eta))[L(t')/L(t)](lambda(')(c)). For the O(n) model in the n-->infinity limit, we show that lambda(')(c)=d+2 and phi=z/2. We give a heuristic argument suggesting that this result is, in fact, valid for any dimension d and spin vector dimension n. We present numerical simulations for the conserved Ising model in d=1 and d=2, which are fully consistent with the present theory.

Experimental study of spinodal decomposition in a 1D conserved order parameter system

We study the zigzag instability coarsening of splay-bend walls formed in a nematic liquid crystal under external fields. The vertexes of zigzag can be considered as kinks in a one-dimensional order parameter system and the geometrical constraints associated with the necessary equal length sum of zig and zag segments impose a conserved quantity in this Cahn-Hilliard-type problem. In the late stage of coarsening, the characteristic length of the system L(t) shows a logarithmic increase in time and the dynamical scaling law holds. We then try to extract the nontrivial asymptotic scaling exponent lambda of the two-time correlation function, defined by lim(<phi(0,t)phi(0,t('))> approximately [L(t)/L(t('))](-lambda). The scaling exponents with respective time references, t(')=32 and 64 s, after quench are found to be lambda approximately 2 which is larger than the value with respective time reference t(')=0, predicted by numerical simulation.

Coarsening dynamics of phase-separating systems

When a system such as a binary liquid is cooled rapidly from a homogeneous phase into a two-phase region, domains of the two equilibrium phases form and grow ('coarsen') with time. In the absence of an external drive, such as gravity or an imposed shear flow, a dynamical-scaling regime emerges in which the domain morphology is statistically self-similar at different times, up to an overall length-scale (coarsening scale) that grows with time. In the first part of the paper, the scaling phenomenology will be reviewed and the time-dependence of the coarsening scale will be discussed in the context of a number of different physical systems and scaling regimes. In the second part, the influence of an external drive, in particular a shear flow, will be addressed and recent developments reviewed. Interesting open questions include the late-time behaviour under shear and whether the coarsening continues indefinitely or is ultimately arrested by the shear flow.

Experimental study of coarsening dynamics of the zigzag wall in a nematic liquid crystal with negative dielectric anisotropy

When a homeotropically aligned nematic liquid crystal cell is placed above two permanent magnets forming a magnetic quadrupole, a straight splay-bend wall, or a so-called Ising wall, is formed. With a material of positive dielectric anisotropy, it has been shown that the application of an electric field perpendicular to the plates leads to a zigzag instability of the wall, exclusively related to the elastic anisotropy of the liquid crystal. In this case, the coarsening process of the zigzag is very slow, which in turn leads to experimental difficulties concerning its quantitative investigation. If a material of negative dielectric anisotropy is used under an electric field with low voltage and low frequency, two convective rolls appear along the Ising wall due to the charge focusing effect, which is also responsible, at a higher voltage in the homogenous tilted regions, for the appearance of Williams domains electrohydrodynamic instability. If the voltage is higher than a threshold value, the straight Ising wall spontaneously breaks into a zigzag shape and a fast coarsening of the zigzag proceeds, associated with the annihilation of two neighboring vertices. In the present paper, the coarsening dynamics of this system, which can be considered as a one-dimensional Ising situation, are investigated experimentally. At late times, the average width of the zigzag increases logarithmically with time. This finding is consistent with the theory and also with the numerical simulation of a one-dimensional Cahn-Hilliard situation having a conserved order parameter. The scaling analysis of size distribution of the Ising domain, the shape of the power spectrum, and of the correlation function of the Ising order parameter, as well as the number density correlation functions of kinks also confirms that the dynamical scaling law predicted for one-dimensional conservative systems holds for the coarsening process. As supposed from symmetry arguments, it is confirmed that this experiment constitutes a one-dimensional analog of spinodal decomposition.

Grain boundary pinning and glassy dynamics in stripe phases

We study numerically and analytically the coarsening of stripe phases in two spatial dimensions, and show that transient configurations do not achieve long ranged orientational order but rather evolve into glassy configurations with very slow dynamics. In the absence of thermal fluctuations, defects such as grain boundaries become pinned in an effective periodic potential that is induced by the underlying periodicity of the stripe pattern itself. Pinning arises without quenched disorder from the nonadiabatic coupling between the slowly varying envelope of the order parameter around a defect, and its fast variation over the stripe wavelength. The characteristic size of ordered domains asymptotes to a finite value R(g) approximately lambda(0)epsilon(-1/2)exp(absolute value of a/square root of epsilon), where epsilon<<1 is the dimensionless distance away from threshold, lambda(0) the stripe wavelength, and a a constant of order unity. Random fluctuations allow defect motion to resume until a new characteristic scale is reached, function of the intensity of the fluctuations. We finally discuss the relationship between defect pinning and the coarsening laws obtained in the intermediate time regime.

Numerical confirmation of late-time t(1/2) growth in three-dimensional phase ordering

Results for the late-time regime of phase ordering in three dimensions are reported, based on numerical integration of the time-dependent Ginzburg-Landau equation with nonconserved order parameter at zero temperature. For very large systems (700(3)) at late times, t> or =150, the characteristic length grows as a power law, R(t) approximately t(n), with the measured n in agreement with the theoretically expected result n=1/2 to within statistical errors. In this time regime R(t) is found to be in excellent agreement with the analytical result of Ohta, Jasnow, and Kawasaki [Phys. Rev. Lett. 49, 1223 (1982)]. At early times, good agreement is found between the simulations and the linearized theory with corrections due to the lattice anisotropy.

Ohta-jasnow-kawasaki approximation for nonconserved coarsening under shear

We analytically study coarsening dynamics in a system with nonconserved scalar order parameter, when a uniform time-independent shear flow is present. We use an anisotropic version of the Ohta-Jasnow-Kawasaki approximation to calculate the growth exponents in two and three dimensions: for d=3 the exponents we find are the same as expected on the basis of simple scaling arguments, that is, 3/2 in the flow direction and 1/2 in all the other directions, while for d=2 we find an unusual behavior, in that the domains experience an unlimited narrowing for very large times and a nontrivial dynamical scaling appears. In addition, we consider the case where an oscillatory shear is applied to a two-dimensional system, finding in this case a standard t(1/2) growth, modulated by periodic oscillations. We support our two-dimensional results by means of numerical simulations and we propose to test our predictions by experiments on twisted nematic liquid crystals.

Corrections to scaling in the phase-ordering dynamics of a vector order parameter

Corrections to scaling, associated with deviations of the order parameter from the scaling morphology in the initial state, are studied for systems with O(n) symmetry at zero temperature in phase-ordering kinetics. Including corrections to scaling, the equal time pair correlation function has the form C(r,t)=f0(r/L)+L(-omega)f1(r/L)+., where L is the coarsening length scale. The correction-to-scaling exponent omega and the correction-to-scaling function f1(x) are calculated for both nonconserved and conserved order parameter systems using the approximate Gaussian closure theory of Mazenko. In general omega is a nontrivial exponent which depends on both the dimensionality d of the system and the number of components n of the order parameter. Corrections to scaling are also calculated for the nonconserved one-dimensional XY model, where an exact solution is possible.