Sparse nonlinear models of chaotic electroconvection

Convection is a fundamental fluid transport phenomenon, where the large-scale motion of a fluid is driven, for example, by a thermal gradient or an electric potential. Modelling convection has given rise to the development of chaos theory and the reduced-order modelling of multiphysics systems; however, these models have been limited to relatively simple thermal convection phenomena. In this work, we develop a reduced-order model for chaotic electroconvection at high electric Rayleigh number. The chaos in this system is related to the standard Lorenz model obtained from Rayleigh-Benard convection, although our system is driven by a more complex three-way coupling between the fluid, the charge density, and the electric field. Coherent structures are extracted from temporally and spatially resolved charge density fields via proper orthogonal decomposition (POD). A nonlinear model is then developed for the chaotic time evolution of these coherent structures using the sparse identification of nonlinear dynamics (SINDy) algorithm, constrained to preserve the symmetries observed in the original system. The resulting model exhibits the dominant chaotic dynamics of the original high-dimensional system, capturing the essential nonlinear interactions with a simple reduced-order model.

Kinetics of motile solitons in nematic liquid crystals

The generation of spatially localized, soliton-like hydrodynamic disturbances in microscale fluidic systems is an intriguing challenge. Herein, we introduce nonequilibrium solitons in nematic liquid crystals stimulated by an electric field. These dynamic solitons are robust as long as the electric field is maintained. Interestingly, their kinetic behaviours depend on the field condition-Tuning of the amplitude and frequency of the applied electric field alters the solitons to self-assemble into lattice ordering like physical particles or to command them to various dynamic states. Our key property to the realisation is the electrohydrodynamic instability due to the coupling between the fluid elasticity and the background convection. This paper describes a new mechanism for realising dynamic solitons in fluid systems on the basis of the electrohydrodynamic phenomena.

Slow diffusive structure in Nikolaevskii turbulence

Weak turbulence has been investigated in nonlinear-nonequilibrium physics to understand universal characteristics near the transition point of ordered and disordered states. Here the one-dimensional Nikolaevskii turbulence, which is a mathematical model of weak turbulence, is studied theoretically. We calculate the velocity field of the Nikolaevskii turbulence assuming a convective structure and carry out tagged-particle simulations in the flow to clarify the Nikolaevskii turbulence from the Lagrangian description. The tagged particle diffuses in the disturbed flow and the diffusion is superdiffusive in an intermediate timescale between ballistic and normal-diffusive scale. The diffusion of the slow structure is characterized by the power law for the control parameter near the transition point of the Nikolaevskii turbulence, suggesting that the diffusive characteristics of the slow structure remain scale invariant. We propose a simplified model, named two-scale Brownian motion, which reveals a hierarchy in the Nikolaevskii turbulence.

Three-dimensional electroconvective vortices in cross flow

This study focuses on the three-dimensional (3D) electrohydrodynamic flow instability between two parallel electrodes driven by unipolar charge injection with and without cross flow. Lattice Boltzmann method with a two-relaxation time model is used to compute flow patterns. In the absence of cross flow, the base-state solution is hydrostatic, and the electric field is one-dimensional. With strong charge injection and high electrical Rayleigh number, the system exhibits electroconvective vortices. Disturbed by perturbation patterns, such as rolling pattern, square pattern, and hexagon pattern, the flow develops corresponding to the most unstable mode. The growth rate and pattern transitions are studied using dynamic mode decomposition of the transient numerical solutions. The interactions between cross flow and electroconvective vortices lead to suppression and disappearance of structures with velocity components in the direction of cross flow, while the other components are not affected. Surprisingly, the transition from a 3D to a 2D flow pattern enhances the convective charge transport, marked by an increase in the electric Nusselt number. Hysteresis in the 3D to 2D transition is characterized by the nondimensional parameter Y, a ratio of the electrical force term to the viscous term in the momentum equation.

Correlations between the leading Lyapunov vector and pattern defects for chaotic Rayleigh-Bénard convection

We probe the effectiveness of using topological defects to characterize the leading Lyapunov vector for a high-dimensional chaotic convective flow field. This is accomplished using large-scale parallel numerical simulations of Rayleigh-Bénard convection for experimentally accessible conditions. We quantify the statistical correlations between the spatiotemporal dynamics of the leading Lyapunov vector and different measures of the flow field pattern's topology and dynamics. We use a range of pattern diagnostics to describe the flow field structures which includes many of the traditional diagnostics used to describe convection as well as some diagnostics tailored to capture the dynamics of the patterns. We use the ideas of precision and recall to build a statistical description of each pattern diagnostic's ability to describe the spatial variation of the leading Lyapunov vector. The precision of a diagnostic indicates the probability that it will locate a region where the Lyapunov vector is larger than a threshold value. The recall of a diagnostic indicates its ability to locate all of the possible spatial regions where the Lyapunov vector is above threshold. By varying the threshold used for the Lyapunov vector magnitude, we generate precision-recall curves which we use to quantify the complex relationship between the pattern diagnostics and the spatiotemporally varying magnitude of the leading Lyapunov vector. We find that pattern diagnostics which include information regarding the flow history outperform pattern diagnostics that do not. In particular, an emerging target defect has the highest precision of all of the pattern diagnostics we have explored.

Time-Dependent Diffusion Coefficients for Chaotic Advection due to Fluctuations of Convective Rolls

The properties of chaotic advection arising from defect turbulence, that is, weak turbulence in the electroconvection of nematic liquid crystals, were experimentally investigated. Defect turbulence is a phenomenon in which fluctuations of convective rolls arise and are globally disturbed while maintaining convective rolls locally. The time-dependent diffusion coefficient, as measured from the motion of a tagged particle driven by the turbulence, was used to clarify the dependence of the type of diffusion on coarse-graining time. The results showed that, as coarse-graining time increases, the type of diffusion changes from superdiffusion → subdiffusion → normal diffusion. The change in diffusive properties over the observed timescale reflects the coexistence of local order and global disorder in the defect turbulence.

Electrically driven three-dimensional solitary waves as director bullets in nematic liquid crystals

Electric field-induced collective reorientation of nematic molecules is of importance for fundamental science and practical applications. This reorientation is either homogeneous over the area of electrodes, as in displays, or periodically modulated, as in electroconvection. The question is whether spatially localized three-dimensional solitary waves of molecular reorientation could be created. Here we demonstrate that the electric field can produce particle-like propagating solitary waves representing self-trapped "bullets" of oscillating molecular director. These director bullets lack fore-aft symmetry and move with very high speed perpendicularly to the electric field and to the initial alignment direction. The bullets are true solitons that preserve spatially confined shapes and survive collisions. The solitons are topologically equivalent to the uniform state and have no static analogs, thus exhibiting a particle-wave duality. Their shape, speed, and interactions depend strongly on the material parameters, which opens the door for a broad range of future studies.

Electroconvection in one-dimensional liquid crystal cells

We investigate the alternating current (ac) -driven electroconvection (EC) in one-dimensional cells (1DCs) under the in-plane switching mode. In 1DCs, defect-free EC can be realized. In the presence and absence of external multiplicative noise, the features of traveling waves (TWs), such as their Hopf frequency f_{H} and velocity, are examined in comparison with those of conventional two-dimensional cells (2DCs) accompanying defects of EC rolls. In particular, we show that the defects significantly contribute to the features of the TWs. Additionally, owing to the defect-free EC in the 1DCs, the effects of the ac and noise fields on the TW are clarified. The ac field linearly increases f_{H}, independent of the ac frequency f. The noise increases f_{H} monotonically, but f_{H} does not vary below a characteristic noise intensity V_{N}^{*}. In addition, soliton-like waves and unfamiliar oscillation of EC vortices in 1DCs are observed, in contrast to the localized EC (called worms) and the oscillation of EC rolls in 2DCs.

Confined Electroconvective and Flexoelectric Instabilities Deep in the Freedericksz State of Nematic CB7CB

We report wormlike flexoelectric structures evolving deep in the Freedericksz state of a nematic layer of the liquid crystal cyanobiphenyl-(CH2)₇-cyanobiphenyl. They form in the predominantly splay-bend thin boundary layers and are built up of solitary flexoelectric domains of the Bobylev-Pikin type. Their formation is possibly triggered by the gradient flexoelectric surface instability that remains optically discernible up to unusually high frequencies. The threshold voltage at which the worms form scales as square root of the frequency; in their extended state, worms often appear as labyrinthine structures on a section of loops that separate regions of opposite director deviation. Such asymmetric loops are also derived through pincement-like dissociation of ring-shaped walls. Formation of isolated domains of bulk electroconvection precedes the onset of surface instabilities. In essence, far above the Freedericksz threshold, the twisted nematic layer behaves as a combination of two orthogonally oriented planar half-layers destabilized by localized flexoelectric distortion.

Traveling waves and worms in ac-driven electroconvection under external multiplicative noise

In the presence of external multiplicative noise, ac-driven electroconvection (EC) in a nematic liquid crystal is investigated. Noise-induced traveling waves (TWs) including localized ones (worms) are found with a typical, stationary wave. Three kinds of TWs are classified by their dynamic features (e.g., noise-intensity-dependent Hopf frequency and space-time map). Moreover, ac frequency-dependent threshold voltages of EC are examined in high noise intensities causing abnormal charge redistribution of the EC cell, and the roles of ac and noise fields with respect to TWs are elucidated in successive pattern evolutions. The mechanism of TWs is discussed in terms of a locally perturbed dynamic conductivity due to the noise field additionally applied to the EC; such a conductivity can be related to a weak-electrolyte model for a Hopf bifurcation to a TW.

Lagrangian chaos and particle diffusion in electroconvection of planar nematic liquid crystals

Two types of spatiotemporal chaos in the electroconvection of nematic liquid crystals, such as defect turbulence and spatiotemporal intermittency, have been statistically investigated according to the Lagrangian picture. Here fluctuations are traced using the motion of a single particle driven by chaotic convection. In the defect turbulence (fluctuating normal rolls), a particle is mainly trapped in a roll but sometimes jumps to a neighboring roll. Its activation energy is then obtained from the jumping (hopping) rate. This research clarifies that diffusion in the defect turbulence regime in electroconvection can be regarded as a kind of hopping process. The spatiotemporal intermittency appears as a coexistent state of ordered grid domains and turbulent domains. The motion of a single particle shows weak and strong diffusion, respectively, in the ordered and turbulent domains. The diffusion characteristics intermittently change from one to another with certain durations as the domains change. This research has found that the distribution function of the duration that a particle remains in an ordered area has a power-law decay for which the index is different from that obtained by the Eulerian measurement.

Localized vegetation patterns, fairy circles, and localized patches in arid landscapes

We investigate the formation of localized structures with varying widths in one- and two-dimensional systems. The mechanism of stabilization is attributed to strongly nonlocal coupling mediated by a Lorentzian type of kernel. We show that, in addition to stable dips found recently [see, e.g. Fernandez-Oto et al., Phys. Rev. Lett. 110, 174101 (2013)], there are stable localized peaks which appear as a result of strongly nonlocal coupling. We applied this mechanism to arid ecosystems by considering a prototype model of a Nagumo type. In one dimension, we study the front connecting the stable uniformly vegetated state to the bare one under the effect of strongly nonlocal coupling. We show that strongly nonlocal coupling stabilizes both-dip and peak-localized structures. We show analytically and numerically that the width of a localized structure, which we interpret as a fairy circle, increases strongly with the aridity parameter. This prediction is in agreement with published observations. In addition, we predict that the width of localized patch decreases with the degree of aridity. Numerical results are in close agreement with analytical predictions.

Localized states in the conserved Swift-Hohenberg equation with cubic nonlinearity

The conserved Swift-Hohenberg equation with cubic nonlinearity provides the simplest microscopic description of the thermodynamic transition from a fluid state to a crystalline state. The resulting phase field crystal model describes a variety of spatially localized structures, in addition to different spatially extended periodic structures. The location of these structures in the temperature versus mean order parameter plane is determined using a combination of numerical continuation in one dimension and direct numerical simulation in two and three dimensions. Localized states are found in the region of thermodynamic coexistence between the homogeneous and structured phases, and may lie outside of the binodal for these states. The results are related to the phenomenon of slanted snaking but take the form of standard homoclinic snaking when the mean order parameter is plotted as a function of the chemical potential, and are expected to carry over to related models with a conserved order parameter.

Depinning, front motion, and phase slips

Pinning and depinning of fronts bounding spatially localized structures in the forced complex Ginzburg-Landau equation describing the 1:1 resonance is studied in one spatial dimension, focusing on regimes in which the structure grows via roll insertion instead of roll nucleation at either edge. The motion of the fronts is nonlocal but can be analyzed quantitatively near the depinning transition.

Steady state-Hopf mode interactions at the onset of electroconvection in the nematic liquid crystal Phase V

We report on a new mode interaction found in electroconvection experiments on the nematic liquid crystal mixture Phase V in planar geometry. The mode interaction (codimension two) point occurs at a critical value of the frequency of the driving AC voltage. For frequencies below this value the primary pattern-forming instability at the onset voltage is an oblique stationary instability involving oblique rolls, and above this value it is an oscillatory instability giving rise to normal traveling rolls (oriented perpendicular to and traveling in the director direction). The transition has been confirmed by measuring the roll angle and the dominant frequency of the time series, as both quantities exhibit a discontinuous jump across zero when the AC frequency is varied near threshold. The globally coupled system of Ginzburg–Landau equations that qualitatively describe this mode interaction is constructed, and the resulting normal form, in which slow spatial variations of the mode amplitudes are ignored, is analyzed. This analysis shows that the Ginzburg–Landau system provides the adequate theoretical description for the experimentally observed phenomenon. The experimentally observed patterns at and higher above the onset allow us to narrow down the range of the parameters in the normal form.

Diagnosis of spatiotemporal chaos in wave envelopes of a nematic electroconvection pattern

In this paper we report and analyze complex spatiotemporal dynamics recorded in electroconvection in the nematic liquid crystal I52, driven by an ac voltage slightly above the onset value. The instability mechanism creating the pattern is an oscillatory (Hopf) instability, giving rise to two pairs of counterpropagating rolls traveling in oblique directions relative to the unperturbed director axis. If a system of nonlinear partial differential equations shows the same set of unstable modes, the pattern above the onset is represented in a weakly nonlinear analysis as a superposition of the traveling rolls in terms of wave envelopes varying slowly in space and time. Motivated by this procedure, we extract slowly varying envelopes from the space-time data of the pattern, using a four-wave demodulation based on Fourier analysis. In order to characterize the spatiotemporal dynamics, we apply a variety of diagnostic methods to the envelopes, including the calculation of mean intensities and correlation lengths, global and local Karhunen-Loève decompositions in Fourier space and physical space, the location of holes, the identification of coherent vertical structures, and estimates of Lyapunov exponents. The results of this analysis provide strong evidence that our pattern exhibits extensive spatiotemporal chaos. One of its main characteristics is the presence of coherent structures of low and high intensities extended in the vertical (parallel to the director) direction.

Self-organization of topological defects due to applied constraints

While topological defects are essential to our understanding of self-organizing periodic systems, little is known about how these systems respond when their defects are subjected to geometrical constraints. In an experiment on spatially modulated thermal convection patterns, we observe that applied geometrical constraints bind topological defects into robust self-localized structures that evolve through aggregation, annihilation, and self-replication. We demonstrate that this unexpected cooperative response to the modulation is a natural consequence of three generic elements: phase locking, symmetry breaking, and spatial resonance. Our work provides new insights into the interplay between order, chaos, and control in self-organizing systems.

Curvature fields, topology, and the dynamics of spatiotemporal chaos

The curvature field is measured from tracer-particle trajectories in a two-dimensional fluid flow that exhibits spatiotemporal chaos and is used to extract the hyperbolic and elliptic points of the flow. These special points are pinned to the forcing when the driving is weak, but wander over the domain and interact in pairs at stronger driving, changing the local topology of the flow. Their behavior reveals a two-stage transition to spatiotemporal chaos: a gradual loss of spatial and temporal order followed by an abrupt onset of topological changes.

Formation of traveling waves in nematics due to material parameter ramps

We propose a simple one-dimensional linear model to describe the formation of localized traveling waves due to material parameter ramps in nematic liquid crystals. We assume that due to some external perturbation the material parameters (conductivity, dielectric constants, viscosity, elasticity) of the liquid crystal become slowly varying functions of position. Temperature gradient in localized regions induced by a laser beam could be one such perturbation. We obtain a 4x4 system of electrohydrodynamic equations [partial derivative equations (PDE) with respect to time and position]. At first we assume that all parameters change from their nonperturbed value by the same ramp function. Then, to reveal the parameters which slowly change are predominant in the formation of localized traveling waves, we obtain four more systems of equations. Each time we assume that just one of the material parameters changes while the others are held constant. Accordingly, by reducing 4x4 systems of equations into one equation and using a WKB-like approach for nonuniform media, we obtain the relevant dispersion relations. We show that ramps of elasticity and dielectric parameters play the dominant role in the formation of localized traveling waves in nonuniform nematic media.

Quantitative and qualitative characterization of zigzag spatiotemporal chaos in a system of amplitude equations for nematic electroconvection

It has been suggested by experimentalists that a weakly nonlinear analysis of the recently introduced equations of motion for the nematic electroconvection by M. Treiber and L. Kramer [Phys. Rev. E 58, 1973 (1998)] has the potential to reproduce the dynamics of the zigzag-type extended spatiotemporal chaos and localized solutions observed near onset in experiments [M. Dennin, D. S. Cannell, and G. Ahlers, Phys. Rev. E 57, 638 (1998); J. T. Gleeson (private communication)]. In this paper, we study a complex spatiotemporal pattern, identified as spatiotemporal chaos, that bifurcates at the onset from a spatially uniform solution of a system of globally coupled complex Ginzburg-Landau equations governing the weakly nonlinear evolution of four traveling wave envelopes. The Ginzburg-Landau system can be derived directly from the weak electrolyte model for electroconvection in nematic liquid crystals when the primary instability is a Hopf bifurcation to oblique traveling rolls. The chaotic nature of the pattern and the resemblance to the observed experimental spatiotemporal chaos in the electroconvection of nematic liquid crystals are confirmed through a combination of techniques including the Karhunen-Loeve decomposition, time-series analysis of the amplitudes of the dominant modes, statistical descriptions, and normal form theory, showing good agreement between theory and experiments.

Onset of electroconvection of homeotropically aligned nematic liquid crystals

We present experimental measurements near the onset of electroconvection (EC) of homeotropically aligned nematic liquid crystals Phase 5A and MBBA. A voltage of amplitude square root 2V0 and frequency f was applied. With increasing V0, EC occurred after the bend Freedericksz transition. We found supercritical bifurcations to EC that were either stationary bifurcations or Hopf bifurcations to traveling convection rolls, depending on the sample conductances. Results for the onset voltages Vc, the critical wave numbers kc, the obliqueness angles thetac, and the traveling-wave (Hopf) frequencies at onset omegac over a range of sample conductances and driving frequencies are presented and compared, to the extent possible, with theoretical predictions. For the most part good agreement was found. However, the experiment revealed some unusual results for the orientations of the convection rolls relative to the direction selected by the Freedericksz domain.

Spatiotemporal chaos in electroconvection of a homeotropically aligned nematic liquid crystal

We present patterns of electroconvection (EC) for the homeotropically aligned nematic liquid crystal MBBA. A voltage V = square root of 2V0 sin(2pift) was applied. With increasing V0, the bend Freedericksz transition at VF was followed by the onset of EC at Vc > VF. We found four distinct pattern types. First, a primary supercritical Hopf bifurcation to traveling waves (TW's) of convection rolls occurred. The structure factor S(k) of this state reflected the azimuthal anisotropy of the underlying Freedericksz state. For f < fL approximately 75 Hz there was a superposition of two oblique-roll modes (pattern I). These patterns were chaotic in space and time. For larger f the patterns consisted of chaotic TW normal rolls (pattern II). Here the chaos was attributable to the motion of dislocations and domain walls between left- and right-traveling waves. A secondary bifurcation yielded pattern III; it had no dominant TW frequency but had broadband chaotic dynamics dominated by the motion of dislocations. This pattern type had been referred to by others as a "chevron pattern;" its structure factor still revealed azimuthal anisotropy. Finally, at somewhat larger identical with epsilon = V2/Vc2 -1 a highly disordered pattern IV with defect dynamics was found. This state had been studied before by Kai and co-workers and was referred to by them as "phase turbulence." It had a structure factor that was (within our resolution) invariant under rotation. For patterns I, II, and III, S(k) contained crescent-shaped peaks. The peak shape was qualitatively different from the case of planar EC where the structure factor has an elliptical cross section. We present measurements of the widths 1/xik and 1/xitheta in the radial (k) and the azimuthal (theta) directions. For small epsilon (patterns I and II) we found that xik was consistent with the usual Ginzburg-Landau scaling xik approximately epsilon(-nuk) with nuk approximately 1/2. However, for xitheta we found xitheta approximately epsilon(-nutheta) with nutheta approximately 3/4. Presumably this anomalous scaling of xitheta is associated with the Goldstone mode of homeotropic EC. We also show data for the height S0 of the structure factor that are consistent with S0 approximately epsilonbeta with beta approximately -0.5, implying that S0 diverges at onset. This differs from the case of domain chaos in rotating Rayleigh-Bénard convection where experiment is consistent with beta = 1/2 and thus with a vanishing S0. The difference between the shape of the structure-factor cross section and between the exponents, for the present case, for planar EC, and for domain chaos suggests that there are different universality classes for spatiotemporal chaos.

Use of dynamic schlieren interferometry to study fluctuations during free diffusion

We used a form of schlieren interferometry to measure the mean-squared amplitude and temporal autocorrelation function of concentration fluctuations driven by the presence of a gradient during the free diffusion of a urea solution into water. By taking and processing sequences of images separated in time by less than the shortest correlation time of interest, we were able to simultaneously measure dynamics at a number of different wave vectors. The technique is conceptually similar to the shadowgraph method, which has been used to make similar measurements, but the schlieren method has the advantage that the transfer function is wave-vector independent rather than oscillatory.

Charge transport scaling in turbulent electroconvection

We describe a local-power-law scaling theory for the mean dimensionless electric current Nu in turbulent electroconvection. The experimental system consists of a weakly conducting, submicron-thick liquid-crystal film supported in the annulus between concentric circular electrodes. It is driven into electroconvection by an applied voltage between its inner and outer edges. At sufficiently large voltage differences, the flow is unsteady and electric charge is turbulently transported between the electrodes. Our theoretical development, which closely parallels the Grossmann-Lohse model for turbulent thermal convection, predicts the local-power law Nu approximately F(gamma)R(gamma)P(delta). R and P are dimensionless numbers that are similar to the Rayleigh and Prandtl numbers of thermal convection, respectively. The dimensionless function F(gamma), which is specified by the model, describes the dependence of Nu on the aspect ratio gamma. We find that measurements of Nu are consistent with the theoretical model.

Influence of excitation wave forms and frequencies on the fundamental time symmetry of the system dynamics, studied in nematic electroconvection

The dynamics of periodically driven nematic electroconvection, a classical dissipative pattern forming system, is studied experimentally and theoretically. We demonstrate that for certain excitation wave forms, the system's dynamic response can be periodic with the excitation or subharmonic, depending on the periodicity of the excitation as control parameter, while for some classes of wave forms, a subharmonic response seems to be principally excluded. In particular, we describe influences of frequency and time symmetry of triangular excitation wave forms. Two intrinsically different routes for the transition to subharmonic dynamics are observed. The time characteristics of the system variables are determined by numerical solution of appropriate model equations and a Floquet analysis. Experimental data are compared to calculations of the model system of two coupled linear differential equations. Results of experiment and model are in excellent agreement.

Statistics of defect trajectories in spatio-temporal chaos in inclined layer convection and the complex Ginzburg-Landau equation

For spatio-temporal chaos observed in numerical simulations of the complex Ginzburg-Landau equation (CGL) and in experiments on inclined-layer convection (ILC) we report numerical and experimental data on the statistics of defects and of defect loops. These loops consist of defect trajectories in space-time that are connected to each other through the pairwise annihilation or creation of the associated defects. While most such loops are small and contain only a few defects, the loop distribution functions decay only slowly with the quantities associated with the loop size, consistent with power-law behavior. For the CGL, two of the three power-law exponents are found to agree, within our computational precision, with those from previous investigations of a simple lattice model. In certain parameter regimes of the CGL and ILC, our results for the single-defect statistics show significant deviations from the previously reported findings that the defect dynamics are consistent with those of random walkers that are created with fixed probability and annihilated through random collisions.

Dynamics of laser-induced electroconvection pulses

We first report that, for planar nematic 4-methoxy-benzilidene-4-butylaniline (MBBA), the electroconvection threshold voltage has a nonmonotonic temperature dependence, with a well-defined minimum, and a slope of about -0.12 V/degrees C near room temperature at 70 Hz. Motivated by this observation, we have designed an experiment in which a weak continuous-wave absorbed laser beam with a diameter comparable to the pattern wavelength generates a locally supercritical region, or pulse, in dye-doped MBBA. Working 10-20 % below the laser-free threshold voltage, we observe a steady-state pulse shaped as an ellipse with the semimajor axis oriented parallel to the nematic director, with a typical size of several wavelengths. The pulse is robust, persisting even when spatially extended rolls develop in the surrounding region, and displays rolls that counterpropagate along the director at frequencies of tenths of Hz, with the rolls on the left (right) side of the ellipse moving to the right (left). Systematic measurements of the sample-voltage dependence of the pulse amplitude, spatial extent, and frequency show a saturation or decrease when the control parameter (evaluated at the center of the pulse) approaches approximately 0.3. We propose that the model for these pulses should be based on the theory of control-parameter ramps, supplemented with new terms to account for the advection of heat away from the pulse when the surrounding state becomes linearly unstable. The advection creates a negative feedback between the pulse size and the efficiency of heat transport, which we argue is responsible for the attenuation of the pulse at larger control-parameter values.

Dependence of domain wall dynamics on background wave number

We report on the growth of domains of standing waves in electroconvection in a nematic liquid crystal, focusing on the evolution of domain walls. An ac voltage is applied to the system, forming an initial state that consists of traveling striped patterns with two different orientations, zig and zag rolls. The standing waves are generated by suddenly applying a periodic modulation of the amplitude of the applied voltage that is approximately resonant with the traveling frequency of the pattern. By varying the modulation frequency, we are able to vary the steady-state, average wave number. We characterize the evolution of the domain walls as a function of the average background wave number by measuring the total area and length of domain walls present in the system as a function of time. We find that as the background wave number is varied away from the "natural" wave number for the pattern, the evolution of the domain walls occurs at a faster rate.

Penta-hepta defect chaos in a model for rotating hexagonal convection

In a model for rotating non-Boussinesq convection with mean flow, we identify a regime of spatiotemporal chaos that is based on a hexagonal planform and is sustained by the induced nucleation of dislocations by penta-hepta defects. The probability distribution function for the number of defects deviates substantially from the usually observed Poisson-type distribution. It implies strong correlations between the defects in the form of density-dependent creation and annihilation rates of defects. We extract these rates from the distribution function and also directly from the defect dynamics.

Patterns of electroconvection in the nematic liquid crystal N4

Electroconvection using the liquid crystal N4 is studied as a function of two control parameters: the applied frequency and the applied voltage. As a function of voltage, there is a rich series of bifurcations that takes the system from stationary rolls to chaos. As a function of the frequency, the initial pattern changes from stationary oblique rolls at low frequencies to stationary normal rolls at higher frequencies. There is also a change in the secondary bifurcations. In particular, we observe that the bimodal-varicose instability is replaced by the skewed-varicose instability as the applied frequency is increased. Comparisons with theoretical predictions are made.

Modulation of localized states in electroconvection

We report on the effects of temporal modulation of the driving force on a particular class of localized states, known as worms, that have been observed in electroconvection in nematic liquid crystals. The worms consist of the superposition of traveling waves and have been observed to have unique, small widths, but to vary in length. The transition from the pure conduction state to worms occurs via a backward bifurcation. A possible explanation of the formation of the worms has been given in terms of coupled amplitude equations. Because the worms consist of the superposition of traveling waves, temporal modulation of the control parameter is a useful probe of the dynamics of the system. We observe that temporal modulation increases the average length of the worms and stabilizes worms below the transition point in the absence of modulation.

Temporal and spatial properties of fluctuations below a supercritical primary bifurcation to traveling oblique-roll electroconvection

We present measurements of thermally-induced oblique-roll traveling-wave (TW) fluctuations below the supercritical primary bifurcation to electroconvection (EC) in the nematic liquid crystal 4-ethyl-2-fluoro-4'-[2-(trans-4-pentylcyclohexyl)ethyl]-biphenyl (I52). First we analyze time sequences of one-dimensional shadowgraph images taken parallel to the director to obtain the TW frequency omega and the fluctuation lifetime tau. Within our resolution we find that omega is independent of epsilon [triple bond] V/V(c)-1 (V is the applied voltage amplitude and V(c) its value at the onset of convection). Contrary to linear theory, the relaxation rate 1/tau remains finite at the bifurcation. Next we present the analysis of temporally uncorrelated two-dimensional shadowgraph images of the fluctuations for several values of the electrical conductivity sigma. We fitted an anisotropic two-dimensional Lorentzian function, corresponding to oblique-roll EC, to the time-averaged structure factors S(k) derived from the images. This yielded information about the components of the mean wave vector k(0) and about the correlation length xi as a function of sigma and epsilon. The angle of obliqueness theta of the roll patterns was independent of sigma but decreased anomalously as epsilon approached zero. The modulus k(0) of k(0) depended on sigma. It also showed an anomalous reduction close to onset. The anomalous epsilon dependence of k(0) and theta disagrees with linear theory, which predicts a smooth, essentially linear dependence on epsilon, and presumably is caused by nonlinear interactions between the fluctuations.

Fluctuation and dissipation in liquid-crystal electroconvection

In this experiment a steady-state current is maintained through a liquid-crystal thin film. When the applied voltage is increased through a threshold, a phase transition is observed to a convective state characterized by the chaotic motion of rolls. Above the threshold, an increase in power consumption is observed that is manifested by an increase in the mean conductivity. A sharp increase in the ratio of the power fluctuations to the mean power dissipated is observed above the transition. This ratio is compared to the predictions of the fluctuation theorem of Gallavotti and Cohen using an effective temperature associated with the rolls' chaotic motion.

Nonequilibrium defect-unbinding transition: defect trajectories and loop statistics

In a Ginzburg-Landau model for parametrically driven waves, a transition between a state of ordered and one of disordered spatiotemporal defect chaos is found. To get insight into the breakdown of the order, the defect trajectories are tracked in detail. Since the defects are created and annihilated in pairs, the trajectories form loops in space-time. The probability distribution functions for the size of the loops and the number of defects involved in them undergo a transition from exponential decay in the ordered regime to a power-law decay in the disordered regime. These power laws are also found in a simple lattice model of randomly created defect pairs that diffuse and annihilate upon collision.

Complex Ginzburg-Landau equation in the presence of walls and corners

We investigate the influence of walls and corners (with Dirichlet and Neumann boundary conditions) in the evolution of two-dimensional autooscillating fields described by the complex Ginzburg-Landau equation. Analytical solutions are found, and arguments provided, to show that Dirichlet walls introduce strong selection mechanisms for the wave pattern. Corners between walls provide additional synchronization mechanisms and associated selection criteria. The numerical results fit well with the theoretical predictions in the parameter range studied.

Domain coarsening in electroconvection

We report on experimental measurements of the growth of regular domains evolving from an irregular pattern in electroconvection. The late-time growth of the domains is consistent with the size of the domains scaling as t(n). We use two isotropic measurements of the domain size: the structure factor and the domain wall length. Measurements using the structure factor are consistent with t(1/5) growth. Measurements using the domain wall length are consistent with t(1/4) growth. One source of this discrepancy is the fact that the distribution of local wave numbers is approximately independent of the domain size. In addition, we measure the anisotropy of the growing domains.

Temporal modulation of the control parameter in electroconvection in the nematic liquid crystal I52

I report on the effects of a periodic modulation of the control parameter on electroconvection in the nematic liquid crystal I52. Without modulation, the primary bifurcation from the uniform state is a direct transition to a state of spatiotemporal chaos. This state is the result of the interaction of four degenerate traveling modes: right and left zig and zag rolls. Periodic modulations of the driving voltage at approximately twice the traveling frequency are used. For a large enough modulation amplitude, standing waves that consist of only zig or zag rolls are stabilized. The standing waves exhibit regular behavior in space and time. Therefore, modulation of the control parameter represents a method of eliminating spatiotemporal chaos. As the modulation frequency is varied away from twice the traveling frequency, standing waves that are a superposition of zig and zag rolls, i.e., standing rectangles, are observed. These results are compared with existing predictions based on coupled complex Ginzburg-Landau equations.

Direct observation of a twist mode in electroconvection

I report on the direct observation of a uniform twist mode of the director field in electroconvection in the nematic liquid crystal I52. Recent theoretical work suggests that such a uniform twist mode of the director field is responsible for a number of secondary bifurcations in both electroconvection and thermal convection in nematics. I show here evidence that the proposed mechanisms are consistent with being the source of the previously reported stationary oblique roll pattern identified as the SO2 state of electroconvection in the liquid crystal I52. The same mechanisms also contribute to a tertiary Hopf bifurcation that I observe in electroconvection in the liquid crystal I52. There are quantitative differences between the experiment and calculations that only include the twist mode. These differences suggest that a complete description must include effects described by the weak-electrolyte model of electroconvection.

Deviations from linear theory for fluctuations below the supercritical primary bifurcation to electroconvection

Over two decades ago it was predicted that nonlinear interactions between thermally driven fluctuations in dissipative nonlinear nonequilibrium systems lead to deviations from mean-field theory. Here we report experimental observations of such deviations as a supercritical primary bifurcation is approached. We measured the mean-square director-angle fluctuations <straight theta(2)> below the bifurcation to electroconvection of two different nematic liquid crystals. For epsilon(mf) identical withV2/V(2)(c,mf)-1 less, similar-0.1 ( V is the applied voltage) we find <straight theta(2)> approximately |epsilon(mf)|(-gamma) with gamma given by linear theory (LT). Closer to the bifurcation there are deviations from LT with a smaller gamma and with V(2)(c)>V(2)(c,mf).

Temporally harmonic oscillons in newtonian fluids

Stationary, highly localized (oscillon) structures are observed in a Newtonian fluid when nonlinear surface waves are parametrically excited with two frequencies. Oscillons have a characteristic structure, that of periodically self-focusing jets. In contrast to previously observed oscillons in highly non-Newtonian media, these states are temporally harmonic with the forcing. For wave amplitudes greater than a critical value, they nucleate from an initial pattern via a hysteretic bifurcation, and can therefore be localized on a background of patterns with a variety of different spatial symmetries.

Defect chaos of oscillating hexagons in rotating convection

Using coupled Ginzburg-Landau equations, the dynamics of hexagonal patterns with broken chiral symmetry are investigated, as they appear in rotating non-Boussinesq or surface-tension-driven convection. We find that close to the secondary Hopf bifurcation to oscillating hexagons the dynamics are well described by a single complex Ginzburg-Landau equation (CGLE) coupled to the phases of the hexagonal pattern. At the band center these equations reduce to the usual CGLE and the system exhibits defect chaos. Away from the band center a transition to a frozen vortex state is found.

Frequency-induced structure transition of nematic electroconvection in twist cells

We investigate electroconvection of nematic liquid crystals in planar 90 degrees twist cells in the dielectric regime. These cells provide competing boundary conditions for the nematic director at both electrodes. When the ratio of cell gap and roll width is not too small (in thin cells and at low frequencies), the convection rolls form along the director in the cell middle, which is diagonal to the anchoring directions. This is in accordance with known behavior in the conduction regime and with the assumption of a bulk instability (rolls traversing the cell). When the ratio of roll periodicity to cell gap is very small in the high-frequency range and the thick cells, a regular convection pattern sets in with wave vectors directed along the two alignment directions. We suggest the interpretation that two independent systems of convection rolls form localized near the electrode plates. It is very likely that a similar behavior occurs in nontwisted cells where it is difficult to identify experimentally.

Bifurcation to worms in electroconvection

The primary bifurcation to electroconvection of the liquid crystal 4-ethyl-2-fluoro-4'-[2-(trans-4-pentylcyclohexyl)ethyl]-biphenyl (I52) with planar alignment leads to localized structures of convection rolls known as "worms" when the conductivity of the fluid is relatively small. Worms coexist with the conduction state. They have a unique small width in the direction perpendicular to the director and a varying, usually much greater, length parallel to the director. Previous experiments had not determined whether the bifurcation to worms is supercritical or subcritical. We estimated the voltage V(c) corresponding to the stability limit of the conduction state by measuring the mean-square amplitude of the thermally induced fluctuations below onset and extrapolating to V(c). We found that worms appear already well below V(c). Thus the bifurcation is subcritical. Measurements of the lifetime of the conduction state below V(c) gave information about the voltage V(s) corresponding to the saddle node below which no worms form. We measured V(c) and epsilon(s)=V(s)2/V(c)2-1 as a function of the conductance sigma for a cell of thickness 24 microm and found for our sample that epsilon(s) approaches zero from negative values near sigma approximately 1.2 x 10(-8) Omega(-1) m(-1) as sigma increases. For larger sigma we found the bifurcation to be supercritical. We have been unable to determine so far whether the experimentally observed transition with decreasing sigma from a supercritical to a subcritical bifurcation occurs via a tricritical bifurcation, or whether the worm saddle node is disconnected from the primary supercritical bifurcation line as suggested by theory.