Gas-liquid nucleation at large metastability: unusual features and a new formalism

Nucleation at large metastability is still largely an unsolved problem, even though it is a problem of tremendous current interest, with wide-ranging practical value, from atmospheric research to materials science. It is now well accepted that the classical nucleation theory (CNT) fails to provide a qualitative picture and gives incorrect quantitative values for such quantities as activation-free energy barrier and supersaturation dependence of nucleation rate, especially at large metastability. In this paper, we present an alternative formalism to treat nucleation at large supersaturation by introducing an extended set of order parameters in terms of the kth largest liquid-like clusters, where k = 1 is the largest cluster in the system, k = 2 is the second largest cluster and so on. At low supersaturation, the size of the largest liquid-like cluster acts as a suitable order parameter. At large supersaturation, the free energy barrier for the largest liquid-like cluster disappears. We identify this supersaturation as the one at the onset of kinetic spinodal. The kinetic spinodal is system-size-dependent. Beyond kinetic spinodal many clusters grow simultaneously and competitively and hence the nucleation and growth become collective. In order to describe collective growth, we need to consider the full set of order parameters. We derive an analytic expression for the free energy of formation of the kth largest cluster. The expression predicts that, at large metastability (beyond kinetic spinodal), the barrier of growth for several largest liquid-like clusters disappears, and all these clusters grow simultaneously. The approach to the critical size occurs by barrierless diffusion in the cluster size space. The expression for the rate of barrier crossing predicts weaker supersaturation dependence than what is predicted by CNT at large metastability. Such a crossover behavior has indeed been observed in recent experiments (but eluded an explanation till now). In order to understand the large numerical discrepancy between simulation predictions and experimental results, we carried out a study of the dependence on the range of intermolecular interactions of both the surface tension of an equilibrium planar gas-liquid interface and the free energy barrier of nucleation. Both are found to depend significantly on the range of interaction for the Lennard-Jones potential, both in two and three dimensions. The value of surface tension and also the free energy difference between the gas and the liquid phase increase significantly and converge only when the range of interaction is extended beyond 6-7 molecular diameters. We find, with the full range of interaction potential, that the surface tension shows only a weak dependence on supersaturation, so the reason for the breakdown of CNT (with simulated values of surface tension and free energy gap) cannot be attributed to the supersaturation dependence of surface tension. This remains an unsettled issue at present because of the use of the value of surface tension obtained at coexistence.

Influence of poling methods on the orientational dynamics of 2-methyl-4-nitro-aniline in poly(methyl methacrylate) studied by second harmonic generation

The magnitude and stability of the induced dipolar orientation of 2-methyl-4-nitroaniline (MNA)/poly(methyl methacrylate) (PMMA) guest/host system is investigated. The chromophores are aligned using both the corona discharge and contact electrode poling techniques. The magnitude of order parameter (also an indicator for the second order nonlinear susceptibility) is measured by recording absorbances of the poled (by the two different techniques) and unpoled PMMA films at different concentrations of MNA. Under the same conditions the corona poling technique creates a higher alignment of molecules along the field direction. The time dependence of the second harmonic intensity of the MNA/PMMA film prepared by the two techniques can be described by a Kohlrausch-Williams-Watts stretched exponential. The temperature dependence of the decay time constant is found to generally follow a modified Williams-Landel-Ferry (WLF) or Vogel-Tamann-Fulcher (VTF) equation. The glass transition temperature seems to be the single most important parameter for determining the relaxation time tau(T).

NMR in liquid crystals spinning at and near magic angle

NMR spectra of liquid crystalline phases and the molecules dissolved therein, spinning at and near the magic angle provide information on the director dynamics and the order parameter. The studies on the dynamics of the liquid crystal director for sample spinning near magic angle in mesophases with positive and negative diamagnetic susceptibility anisotropies (Delta chi) and their mixtures with near-zero macroscopic diamagnetic susceptibility anisotropies have been reported. In systems with weakly positive Delta chi, the director has been observed to switch from an orientation parallel to the spinning axis at low rotational speeds to one perpendicular to the spinning axis at high rotational speeds, when the angle theta, the axis of rotation makes with the magnetic field is smaller than the magic angle theta(m). For systems with a small negative Delta chi, similar director behaviour has been observed for theta greater than theta(m). At magic angle, the spectra under slow spinning speeds exhibit a centre band and side bands at integral values of the spinning speeds. The intensities of the spinning side bands have been shown to contain information on the sign and the magnitude of the order parameter(s). The results are discussed with illustrative examples. Results on the orientation of the chemical shielding tensor obtained from a combination of the NMR studies in the solid and the liquid crystalline states, have been described.

Nonequilibrium phase transition in the kinetic Ising model: Critical slowing down and the specific-heat singularity

The nonequilibrium dynamic phase transition, in the kinetic Ising model in the presence of an oscillating magnetic field has been studied both by Monte Carlo simulation and by solving numerically the mean-field dynamic equation of motion for the average magnetization. In both cases, the Debye ''relaxation'' behavior of the dynamic order parameter has been observed and the ''relaxation time'' is found to diverge near the dynamic transition point. The Debye relaxation of the dynamic order parameter and the power law divergence of the relaxation time have been obtained from a very approximate solution of the mean-field dynamic equation. The temperature variation of appropriately defined ''specific heat'' is studied by the Monte Carlo simulation near the transition point. The specific heat has been observed to diverge near the dynamic transition point.

Nonequilibrium phase transition in the kinetic Ising model: Divergences of fluctuations and responses near the transition point

The nonequilibrium dynamic phase transition in the kinetic Ising model in the presence of an oscillating magnetic field is studied by Monte Carlo simulation. The fluctuation of the dynamic older parameter is studied as a function of temperature near the dynamic transition point. The temperature variation of appropriately defined ''susceptibility'' is also studied near the dynamic transition point. Similarly, the fluctuation of energy and appropriately defined ''specific heat'' is studied as a function of temperature near the dynamic transition point. In both cases, the fluctuations (of dynamic order parameter and energy) and the corresponding responses diverge (in power law fashion) near the dynamic transition point with similar critical behavior (with identical exponent values).

Phenomenological aspects of phase fluctuation in high Tc superconductors: Cooper pair droplets

We point out how fluctuation of the phase of the superconducting order parameter can play a key role in our understanding of high Te superconductors. A simple universal criterion is given which illustrates why all oxide superconductors in contrast to classical superconductors ought to behave as a lattice of cooper pairs. T-c is to be thought of as the temperature of phase coherence or the temperature above which the lattice of Cooperpair 'melts' into a phase of Cooper-pair droplets that starts forming at T approximate to T-* . This is the pseudo-gap region. Quantum fluctuation of the phase predicts a superconductor to insulator phase transition for all underdoped materials.

Nematogens with more flexible chains than aromatic rings in the core

Mesogens containing four rings in the main core can accommodate one terminal and two nearby lateral chains on each outside aromatic ring. These compounds containing six chains present an enantiotropic nematic range which is influenced by the rigidity of the links. The conformational behaviour of the first methyleneoxy group within the chains was investigated by one and two dimensional C-13 NMR. The sign of the jump in chemical shifts when entering the nematic phase indicates the folding of each lateral branch. Dipolar oscillations during cross-polarization contact provide the values of the bond order parameter. The two First lateral fragments do not behave in the same way, demonstrating the influence of the fragment along which the chain is back: folded.

Phase separation in a simple model with dynamical asymmetry

We perform computer simulations of a Cahn-Hilliard model of phase separation that has dynamical asymmetry between the two coexisting phases. The dynamical asymmetry is incorporated by considering a mobility function that is order parameter dependent. Simulations of this model reveal morphological features similar to those observed in viscoelastic phase separation. In the early stages, the minority phase domains form a percolating structure that shrinks with time, eventually leading to the formation of disconnected regions that are characterized by the presence of random interfaces as well as isolated droplets. The domains grow as L(t)similar to t(1/3) in the very late stages. Although dynamical scaling is violated in the area shrinking regime, it is restored at late times. However, the form of the scaling function is found to depend on the extent of dynamical asymmetry. [S1063-651X(99)12101-9].

Dynamics of strain bifurcations in a magnetostrictive ribbon

We develop a coupled nonlinear oscillator model involving magnetization and strain to explain several experimentally observed dynamical features exhibited by forced magnetostrictive ribbon. Here we show that the model recovers the observed period-doubling route to chaos as function of the dc field for a fixed ac field and quasiperiodic route to chaos as a function of the ac field, keeping the dc field constant. The model also predicts induced and suppressed chaos under the influence of an additional small-amplitude near-resonant ac field. Our analysis suggests rich dynamics in coupled order-parameter systems such as magnetomartensitic and magnetoelectric materials.

Transport across a junction of topological insulators and a superconductor

We study transport across a line junction lying between two orthogonal topological insulator surfaces and a superconductor which can have either s-wave (spin-singlet) or p-wave (spin-triplet) pairing symmetry. The junction can have three time-reversal invariant barriers on three sides. We compute the charge and the spin conductance across such a junction and study their behaviors as a function of the bias voltage applied across the junction and the three parameters used to characterize the barrier. We find that the presence of topological insulators and a superconductor leads to both Dirac- and Schrodinger-like features in charge and spin conductances. We discuss the effect of bound states on the superconducting side of the barrier on the conductance; in particular, we show that for triplet p-wave superconductors, such a junction may be used to determine the spin state of its Cooper pairs. Our study reveals that there is a nonzero spin conductance for some particular spin states of the triplet Cooper pairs; this is an effect of the topological insulators which break the spin rotation symmetry. Finally, we find an unusual satellite peak (in addition to the usual zero bias peak) in the spin conductance for p-wave symmetry of the superconductor order parameter.

Solid-liquid transition in polydisperse Lennard-Jones systems

We study melting of a face-centered crystalline solid consisting of polydisperse Lennard-Jones spheres with Gaussian polydispersity in size. The phase diagram reproduces the existence of a nearly temperature invariant terminal polydispersity (delta(t) similar or equal to 0.11), with no signature of reentrant melting. The absence of reentrant melting can be attributed to the influence of the attractive part of the potential upon melting. We find that at terminal polydispersity the fractional density change approaches zero, which seems to arise from vanishingly small compressibility of the disordered phase. At constant temperature and volume fraction the system undergoes a sharp transition from crystalline solid to the disordered amorphous or fluid state with increasing polydispersity. This has been quantified by second- and third-order rotational invariant bond orientational order, as well as by the average inherent structure energy. The translational order parameter also indicates a similar sharp structural change at delta similar or equal to 0.09 in case of T* = 1.0, phi = 0.58. The free energy calculation further supports the sharp nature of the transition. The third-order rotationally invariant bond order shows that with increasing polydispersity, the local cluster favors a more icosahedral arrangement and the system loses its local crystalline symmetry. Interestingly, the value of structure factor S(k) of the amorphous phase at delta similar or equal to 0.10 (just beyond the solid-liquid transition density at T* = 1) becomes 2.75, which is below the value of 2.85 required for freezing given by the empirical Hansen-Verlet rule of crystallization, well known in the theory of freezing.

Nucleation of a Stable Solid from Melt in the Presence of Multiple Metastable Intermediate Phases: Wetting, Ostwald's Step Rule, and Vanishing Polymorphs

In many systems, nucleation of a stable solid may occur in the presence of other (often more than one) metastable phases. These may be polymorphic solids or even liquid phases. Sometimes, the metastable phase might have a lower free energy minimum than the liquid but higher than the stable-solid-phase minimum and have characteristics in between the parent liquid and the globally stable solid phase. In such cases, nucleation of the solid phase from the melt may be facilitated by the metastable phase because the latter can ``wet'' the interface between the parent and the daughter phases, even though there may be no signature of the existence of metastable phase in the thermodynamic properties of the parent liquid and the stable solid phase. Straightforward application of classical nucleation theory (CNT) is flawed here as it overestimates the nucleation barrier because surface tension is overestimated (by neglecting the metastable phases of intermediate order) while the thermodynamic free energy gap between daughter and parent phases remains unchanged. In this work, we discuss a density functional theory (DFT)-based statistical mechanical approach to explore and quantify such facilitation. We construct a simple order-parameter-dependent free energy surface that we then use in DFT to calculate (i) the order parameter profile, (ii) the overall nucleation free energy barrier, and (iii) the surface tension between the parent liquid and the metastable solid and also parent liquid and stable solid phases. The theory indeed finds that the nucleation free energy barrier can decrease significantly in the presence of wetting. This approach can provide a microscopic explanation of the Ostwald step rule and the well-known phenomenon of ``disappearing polymorphs'' that depends on temperature and other thermodynamic conditions. Theory reveals a diverse scenario for phase transformation kinetics, some of which may be explored via modem nanoscopic synthetic methods.

Hard-core bosons in a zig-zag optical superlattice

We study a system of hard-core bosons at half-filling in a one-dimensional optical superlattice. The bosons are allowed to hop to nearest-and next-nearest-neighbor sites. We obtain the ground-state phase diagram as a function of microscopic parameters using the finite-size density-matrix renormalization-group method. Depending on the sign of the next-nearest-neighbor hopping and the strength of the superlattice potential the system exhibits three different phases, namely the bond-order (BO) solid, the superlattice induced Mott insulator (SLMI), and the superfluid (SF) phase. When the signs of both hopping amplitudes are the same (the unfrustratedase), the system undergoes a transition from the SF to the SLMI at a nonzero value of the superlattice potential. On the other hand, when the two amplitudes differ in sign (the frustrated case), the SF is unstable to switching on a superlattice potential and also exists only up to a finite value of the next-nearest-neighbor hopping. This part of the phase diagram is dominated by the BO phase which breaks translation symmetry spontaneously even in the absence of the superlattice potential and can thus be characterized by a bond-order parameter. The transition from BO to SLMI appears to be first order.

Strong-coupling expansion for ultracold bosons in an optical lattice at finite temperatures in the presence of superfluidity

We develop a strong-coupling (t << U) expansion technique for calculating the density profile for bosonic atoms trapped in an optical lattice with an overall harmonic trap at finite temperature and finite on-site interaction in the presence of superfluid regions. Our results match well with quantum Monte Carlo simulations at finite temperature. We also show that the superfluid order parameter never vanishes in the trap due to the proximity effect. Our calculations for the scaled density in the vacuum-to-superfluid transition agree well with the experimental data for appropriate temperatures. We present calculations for the entropy per particle as a function of temperature which can be used to calibrate the temperature in experiments. We also discuss issues connected with the demonstration of universal quantum critical scaling in the experiments.

Fluctuating micro-heterogeneity in water-tert-butyl alcohol mixtures and lambda-type divergence of the mean cluster size with phase transition-like multiple anomalies

Water-tert-butyl alcohol (TBA) binary mixture exhibits a large number of thermodynamic and dynamic anomalies. These anomalies are observed at surprisingly low TBA mole fraction, with x(TBA) approximate to 0.03-0.07. We demonstrate here that the origin of the anomalies lies in the local structural changes that occur due to self-aggregation of TBA molecules. We observe a percolation transition of the TBA molecules at x(TBA) approximate to 0.05. We note that ``islands'' of TBA clusters form even below this mole fraction, while a large spanning cluster emerges above that mole fraction. At this percolation threshold, we observe a lambda-type divergence in the fluctuation of the size of the largest TBA cluster, reminiscent of a critical point. Alongside, the structure of water is also perturbed, albeit weakly, by the aggregation of TBA molecules. There is a monotonic decrease in the tetrahedral order parameter of water, while the dipole moment correlation shows a weak nonlinearity. Interestingly, water molecules themselves exhibit a reverse percolation transition at higher TBA concentration, x(TBA) approximate to 0.45, where large spanning water clusters now break-up into small clusters. This is accompanied by significant divergence of the fluctuations in the size of largest water cluster. This second transition gives rise to another set of anomalies around. Both the percolation transitions can be regarded as manifestations of Janus effect at small molecular level. (C) 2014 AIP Publishing LLC.