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.

Interaction of sheep liver apo-serine hydroxymethyltransferase with pyridoxal-5'-phosphate: A physicochemical, kinetic, and thermodynamic study

Sheep liver serine hydroxymethyltransferase (EC 2.1.2.1) is a homotetramer of M(r) 213,000 requiring pyridoxal-5'-phosphate (PLP) as cofactor, Removal of PLP from the holoenzyme converted the enzyme to the apo form which, in addition to being inactive, was devoid of the characteristic absorption spectrum. Upon the addition of PLP to the apoenzyme, complete activity was restored and the visible absorption spectrum with a maximum at 425 nm was regained. The interaction of PLP with the apoenzyme revealed two phases of reaction with pseudo-first-order rate constants of 20 +/- 5 s(-1) and 12.2 +/- 2.0 x 10(-3) s(-1), respectively. However, addition of PLP to the apoenzyme did not cause gross conformational changes as evidenced by circular dichroic and fluorescence spectroscopy. Although conformationally apoenzyme and holoenzyme were indistinguishable, they had distinct apparent melting temperatures of 51 +/- 2 and 58 +/- 2 degrees C, respectively, and the reconstituted holoenzyme was thermally as stable as the native holoenzyme. These results suggested that there was no apparent difference in the secondary structure of holoenzyme, apoenzyme, and reconstituted holoenzyme, However, sedimentation analysis of the apoenzyme revealed the presence of two peaks of S-20,S-w values of 8.7 +/- 0.5 and 5.7 +/- 0.3 S, respectively. A similar pattern was observed when the apoenzyme was chromatographed on a calibrated Sephadex G-150 column. The first peak corresponded to the tetrameric form (M(r) 200,000 +/- 15,000) while the second peak had a M(r) of 130,000 +/- 10,000. Reconstitution experiments revealed that only the tetrameric form of the apoenzyme could be converted into an active holoenzyme while the dimeric form could not be reconstituted into an active enzyme. These results demonstrate that PLP plays an important role in maintaining the structural integrity of the enzyme by preventing the dissociation of the enzyme into subunits, in addition to its function in catalysis. (C) 1996 Academic Press, Inc.

KFVS (Kinetic Flux Vector Splitting) on moving grids for unsteady aerodynamics

Accurate numerical solutions to the problems in fluid-structure (aeroelasticity) interaction are becoming increasingly important in recent years. The methods based on FCD (Fixed Computational Domain) and ALE (Alternate Lagrangian Eulerian) to solve such problems suffer from numerical instability and loss of accuracy. They are not general and can not be extended to the flowsolvers on unstructured meshes. Also, global upwind schemes can not be used in ALE formulation thus leads to the development of flow solvers on moving grids. The KFVS method has been shown to be easily amenable on moving grids required in unsteady aerodynamics. The ability of KFMG (Kinetic Flux vector splitting on Moving Grid) Euler solver in capturing shocks, expansion waves with small and very large pressure ratios and contact discontinuities has been demonstrated.

Potential-distance relationships of clay-water systems considering the Stern theory

This paper deals with the use of Stem theory as applied to a clay-water electrolyte system, which is more realistic to understand the force system at micro level man the Gouy-Chapman theory. The influence of the Stern layer on potential-distance relationship has been presented quantitatively for certain specified clay-water systems and the results are compared with the Gouy-Chapman model. A detailed parametric study concerning the number of adsorption spots on the clay platelet, the thickness of the Stern layer, specific adsorption potential and the value of dielectric constant of the pore fluid in the Stern layer, was carried out. This study investigates that the potential obtained at any distance using the Stern theory is higher than that obtained by the Gouy-Chapman theory. The hydrated size of the ion is found to have a significant influence on the potential-distance relationship for a given clay, pore fluid characteristics and valence of the exchangeable ion.

Realtime kinetic analysis of antigen-antibody interaction using solid phase binding: Transformation of hCG-monoclonal antibody complex

Kinetic constants of MAb-hCG interactions have been determined using solid phase binding of I-125[hCG] to immobilized MAb. While association has been shown to follow the expected pattern, dissociation consists of at least two reversible steps, one with a rate constant of 0.0025 min(-1), and a second with a rate constant of 0.00023 min(-1). Validity of affinity constant measurements in the light of the complex reaction kinetics is discussed, A comparison between the method of surface plasmon resonance technology (BIAcore) and solid phase binding (SPB) for determination of kinetic parameters shows that SPB provides not only a cost-effective approach for determination of realtime kinetic parameters of macromolecular ligand-ligate interaction but also a method with several advantages over the BIAcore system in investigating the mechanism of antigen-antibody interaction.

The structural theory and physical organic chemistry

Woolley's revolutionary proposal that quantum mechanics does not sanction the concept of ''molecular structure'' - which is but only a ''metaphor'' - has fundamental implications for physical organic chemistry. On the one hand, the Uncertainty Principle limits the precision with which transition state structures may be defined; on the other, extension of the structure concept to the transition state may be unviable. Attempts to define transition states have indeed caused controversy. Consequences for molecular recognition, and a mechanistic classification, are also discussed.

Instability, intermittency, and multiscaling in discrete growth models of kinetic roughening

We show by numerical simulations that discretized versions of commonly studied continuum nonlinear growth equations (such as the Kardar-Parisi-Zhangequation and the Lai-Das Sarma-Villain equation) and related atomistic models of epitaxial growth have a generic instability in which isolated pillars (or grooves) on an otherwise flat interface grow in time when their height (or depth) exceeds a critical value. Depending on the details of the model, the instability found in the discretized version may or may not be present in the truly continuum growth equation, indicating that the behavior of discretized nonlinear growth equations may be very different from that of their continuum counterparts. This instability can be controlled either by the introduction of higher-order nonlinear terms with appropriate coefficients or by restricting the growth of pillars (or grooves) by other means. A number of such ''controlled instability'' models are studied by simulation. For appropriate choice of the parameters used for controlling the instability, these models exhibit intermittent behavior, characterized by multiexponent scaling of height fluctuations, over the time interval during which the instability is active. The behavior found in this regime is very similar to the ''turbulent'' behavior observed in recent simulations of several one- and two-dimensional atomistic models of epitaxial growth.

Nonequilibrium-phase transition and 'specific-heat' singularity in the kinetic Ising model: A Monte Carlo study

The nonequilibrium-phase transition has been studied by Monte Carlo simulation in a ferromagnetically interacting (nearest-neighbour) kinetic Ising model in presence of a sinusoidally oscillating magnetic field. The ('specific-heat') temperature derivative of energies (averaged over a full cycle of the oscillating field) diverge near the dynamic transition point.

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).

Modified Kirchhoffs theory of plates including transverse shear deformations

A primary flexure problem defined by Kirchhoff theory of plates in bending is considered. Significance of auxiliary function introduced earlier in the in-plane displacements in resolving Poisson-Kirchhoffs boundary conditions paradox is reexamined with reference to reported sixth order shear deformation theories, in particular, Reissner's theory and Hencky's theory. Sixth order modified Kirchhoff's theory is extended here to include shear deformations in the analysis. (C) 2011 Elsevier Ltd. All rights reserved.

Solvation dynamics in monohydroxy alcohols: Agreement between theory and different experiments

Recently three different experimental studies on ultrafast solvation dynamics in monohydroxy straight-chain alcohols (C-1-C-4) have been carried out, with an aim to quantify the time constant (and the amplitude) of the ultrafast component. The results reported are, however, rather different from different experiments. In order to understand the reason for these differences, we have carried out a detailed theoretical study to investigate the time dependent progress of solvation of both an ionic and a dipolar solute probe in these alcohols. For methanol, the agreement between the theoretical predictions and the experimental results [Bingemann and Ernsting J. Chem. Phys. 1995, 102, 2691 and Horng et al. J: Phys, Chern, 1995, 99, 17311] is excellent. For ethanol, propanol, and butanol, we find no ultrafast component of the time constant of 70 fs or so. For these three liquids, the theoretical results are in almost complete agreement with the experimental results of Horng et al. For ethanol and propanol, the theoretical prediction for ionic solvation is not significantly different from that of dipolar solvation. Thus, the theory suggests that the experiments of Bingemann and Ernsting and those of Horng et al. studied essentially the polar solvation dynamics. The theoretical studies also suggest that the experimental investigations of Joo et al. which report a much faster and larger ultrafast component in the same series of solvents (J. Chem. Phys. 1996, 104, 6089) might have been more sensitive to the nonpolar part of solvation dynamics than the polar part. In addition, a discussion on the validity of the present theoretical approach is presented. In this theory the ultrafast component arises from almost frictionless inertial motion of the individual solvent molecules in the force field of its neighbors.

Some recent advances in the theory of homogeneous isotropic turbulence

We review some advances in the theory of homogeneous, isotropic turbulence. Our emphasis is on the new insights that have been gained from recent numerical studies of the three-dimensional Navier Stokes equation and simpler shell models for turbulence. In particular, we examine the status of multiscaling corrections to Kolmogorov scaling, extended self similarity, generalized extended self similarity, and non-Gaussian probability distributions for velocity differences and related quantities. We recount our recent proposal of a wave-vector-space version of generalized extended self similarity and show how it allows us to explore an intriguing and apparently universal crossover from inertial- to dissipation-range asymptotics.

Strong shock approximation in the new theory of shock dynamics

This paper investigates the propagation of a strong shock into an inhomogeneous medium using the new theory of shock dynamics. The equations are simple to solve and involve no trial-and-error method commonly used in this case. The results compare favourably with earlier results obtained in the case of self-similar flows, which arise as a special case of this theory.