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.

Multifractal analysis, a method to investigate the morphology of materials

Using the multifractal formalism, we discuss the results obtained to characterized the morphology of polymer alloys and granular discontinuous metallic thin films. In the first case we have found a correlation between the multifractality and the mechanical properties of the alloys. In the second case, we have found that it is possible to measure the differences between the morphology of thin films induced by a growth process on a subtrate and that of percolation clusters of the classical theory of percolation.

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.

Mean first passage time approach to the problem of optimal barrier subdivision for Kramer's escape rate

We consider the effect of subdividing the potential barrier along the reaction coordinate on Kramer's escape rate for a model potential, Using the known supersymmetric potential approach, we show the existence of an optimal number of subdivisions that maximizes the rate, We cast the problem as a mean first passage time problem of a biased random walker and obtain equivalent results, We briefly summarize the results of our investigation on the increase in the escape rate by placing a blow-torch in the unstable part of one of the potential wells. (C) 1999 Elsevier Science B.V. All rights reserved.

Stochastic dynamics modeling of the protein sequence length distribution in genomes: implications for microbial evolution

In this paper, we report an analysis of the protein sequence length distribution for 13 bacteria, four archaea and one eukaryote whose genomes have been completely sequenced, The frequency distribution of protein sequence length for all the 18 organisms are remarkably similar, independent of genome size and can be described in terms of a lognormal probability distribution function. A simple stochastic model based on multiplicative processes has been proposed to explain the sequence length distribution. The stochastic model supports the random-origin hypothesis of protein sequences in genomes. Distributions of large proteins deviate from the overall lognormal behavior. Their cumulative distribution follows a power-law analogous to Pareto's law used to describe the income distribution of the wealthy. The protein sequence length distribution in genomes of organisms has important implications for microbial evolution and applications. (C) 1999 Elsevier Science B.V. All rights reserved.

The 2D Coulomb gas on a 1D lattice

The statistical mechanics of a two-dimensional Coulomb gas confined to one dimension is studied, wherein hard core particles move on a ring. Exact self-duality is shown for a version of the sine-Gordon model arising in this context, thereby locating the transition temperature exactly. We present asymptotically exact results for the correlations in the model and characterize the low- and high-temperature phases. Numerical simulations provide support to these renormalization group calculations. Connections with other interesting problems, such as the quantum Brownian motion of a panicle in a periodic potential and impurity problems, are pointed out.

Glassy behavior in neural network models of associative memory

Neural network models of associative memory exhibit a large number of spurious attractors of the network dynamics which are not correlated with any memory state. These spurious attractors, analogous to "glassy" local minima of the energy or free energy of a system of particles, degrade the performance of the network by trapping trajectories starting from states that are not close to one of the memory states. Different methods for reducing the adverse effects of spurious attractors are examined with emphasis on the role of synaptic asymmetry. (C) 2002 Elsevier Science B.V. All rights reserved.

The extended Enskog operator for simple fluids with continuous potentials: single particle and collective properties

We generalized the Enskog theory originally developed for the hard-sphere fluid to fluids with continuous potentials, such as the Lennard–Jones. We derived the expression for the k and ω dependent transport coefficient matrix which enables us to calculate the transport coefficients for arbitrary length and time scales. Our results reduce to the conventional Chapman–Enskog expression in the low density limit and to the conventional k dependent Enskog theory in the hard-sphere limit. As examples, the self-diffusion of a single atom, the vibrational energy relaxation, and the activated barrier crossing dynamics problem are discussed.

Understanding glassy dynamics from the free-energy landscape

Many interesting features of the dynamics of simple liquids near the glass transition may be understood in terms of properties of the free-energy landscape obtained from numerical studies of a model free-energy functional. Main results obtained from this approach are summarized and a list of references to relevant publications is provided. (C) 2002 Elsevier Science B.V. All rights reserved.

Utilizing Model Nanostructured Porous Inorganic Compounds Such as Silica to Beneficially Confine “Bio”-Relevant Molecules and Demonstrate their Potential in Therapeutics

Abstract | The importance of well-defined inorganic porous nanostructured materials in the context of biotechnological applications such as drug delivery and biomolecular sensing is reviewed here in detail. Under optimized conditions, the confinement of “bio”-relevant molecules such as pharmaceutical drugs, enzymes or proteins inside such inorganic nanostructures may be remarkably beneficial leading to enhanced molecular stability, activity and performance. From the point of view of basic research, molecular confinement inside nanostructures poses several formidable and intriguing problems of statistical mechanics at the mesoscopic scale. The theoretical comprehension of such non-trivial issues will not only aid in the interpretation of observed phenomena but also help in designing better inorganic nanostructured materials for biotechnological applications.

Dynamics of end-to-end loop formation: A flexible chain in the presence of hydrodynamic interaction

Based on the Wilemski-Fixman approach G. Wilemski, M. Fixman, J. Chem. Phys. 60 (1974) 866], we show that, for a flexible chain in theta solvent, hydrodynamic interaction treated with a pre-averaging approximation makes ring closing faster if the chain is not very short. We also show that the ring closing time for a long chain with hydrodynamic interaction in theta solvent scales with the chain length (N) as N-1.5, in agreement with the previous renormalization group calculation based prediction by Freidman and O'Shaughnessy B. Friedman, B. O'Shaughnessy, Phys. Rev. A 40 (1989) 5950]. (C) 2012 Elsevier B.V. All rights reserved.

Reducing convergence times of self-propelled swarms via modified nearest neighbor rules

Vicsek et al. proposed a biologically inspired model of self-propelled particles, which is now commonly referred to as the Vicsek model. Recently, attention has been directed at modifying the Vicsek model so as to improve convergence properties. In this paper, we propose two modification of the Vicsek model which leads to significant improvements in convergence times. The modifications involve an additional term in the heading update rule which depends only on the current or the past states of the particle's neighbors. The variation in convergence properties as the parameters of these modified versions are changed are closely investigated. It is found that in both cases, there exists an optimal value of the parameter which reduces convergence times significantly and the system undergoes a phase transition as the value of the parameter is increased beyond this optimal value. (C) 2012 Elsevier B.V. All rights reserved.

The force distribution function of an oscillator model of polymer stretching at constant velocity

Recent simulations of the stretching of tethered biopolymers at a constant speed v (Ponmurugan and Vemparala, 2011 Phys. Rev. E 84 060101(R)) have suggested that for any time t, the distribution of the fluctuating forces f responsible for chain deformation is governed by a relation of the form P(+ f)/ P(- f) = expgamma f], gamma being a coefficient that is solely a function of v and the temperature T. This result, which is reminiscent of the fluctuation theorems applicable to stochastic trajectories involving thermodynamic variables, is derived in this paper from an analytical calculation based on a generalization of Mazonka and Jarzynski's classic model of dragged particle dynamics Mazonka and Jarzynski, 1999 arXiv:cond-\textbackslashmat/9912121v1]. However, the analytical calculations suggest that the result holds only if t >> 1 and the force fluctuations are driven by white rather than colored noise; they further suggest that the coefficient gamma in the purported theorem varies not as v(0.15)T-(0.7), as indicated by the simulations, but as vT(-1).

Yielding and large deviations in micellar gels: a model

We present a simple model that can be used to account for the rheological behaviour observed in recent experiments on micellar gels. The model combines attachment detachment kinetics with stretching due to shear, and shows well-defined jammed and flowing states. The large-deviation function (LDF) for the coarse-grained velocity becomes increasingly non-quadratic as the applied force F is increased, in a range near the yield threshold. The power fluctuations are found to obey a steady-state fluctuation relation (FR) at small F. However, the FR is violated when F is near the transition from the flowing to the jammed state although the LDF still exists; the antisymmetric part of the LDF is found to be nonlinear in its argument. Our approach suggests that large fluctuations and motion in a direction opposite to an imposed force are likely to occur in a wider class of systems near yielding.

Chapman-Enskog analysis for finite-volume formulation of lattice Boltzmann equation

The classical Chapman-Enskog expansion is performed for the recently proposed finite-volume formulation of lattice Boltzmann equation (LBE) method D.V. Patil, K.N. Lakshmisha, Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh, J. Comput. Phys. 228 (2009) 5262-5279]. First, a modified partial differential equation is derived from a numerical approximation of the discrete Boltzmann equation. Then, the multi-scale, small parameter expansion is followed to recover the continuity and the Navier-Stokes (NS) equations with additional error terms. The expression for apparent value of the kinematic viscosity is derived for finite-volume formulation under certain assumptions. The attenuation of a shear wave, Taylor-Green vortex flow and driven channel flow are studied to analyze the apparent viscosity relation.