Crossover dynamics at large metastability in gas-liquid nucleation

We have developed an alternate description of dynamics of nucleation in terms of an extended set of order parameters. The order parameters consist of an ordered set of kth largest clusters, ordered such that k = 1 is the largest cluster in the system, k = 2 is the second largest cluster, and so on. We have derived an analytic expression for the free energy for the kth largest cluster, which is in excellent agreement with the simulated results. At large supersaturation, the free energy barrier for the growth of the kth largest cluster disappears and the nucleation becomes barrierless. The major success of this extended theoretical formalism is that it can clearly explain the observed change in mechanism at large metastability P. Bhimalapuram et al., Phys. Rev. Lett. 98, 206104 (2007)] and the associated dynamical crossover. The classical nucleation theory cannot explain this crossover. The crossover from activated to barrierless nucleation is found to occur at a supersaturation where multiple clusters cross the critical size. We attribute the crossover as the onset of the kinetic spinodal. We have derived an expression for the rate of nucleation in the barrierless regime by modeling growth as diffusion on the free energy surface of the largest cluster. The model reproduces the slower increase in the rate of growth as a function of supersaturation, as observed in experiments.

Thermodynamic and structural consequences of changing a sulfur atom to a methylene group in the M13Nle mutation in ribonuclease-S

Two fragments of pancreatic ribonuclease A, a truncated version of S-peptide (residues 1-15) and S-protein (residues 21-124), combine to give a catalytically active complex. We have substituted the wild-type residue at position 13, methionine (Met), with norleucine (Nle), where the only covalent change is the replacement of the sulfur atom with a methylene group. The thermodynamic parameters associated with the binding of this variant to S-protein, determined by titration calorimetry in the temperature range 10-40 degrees C, are reported and compared to values previously reported [Varadarajan, R., Connelly, P. R., Sturtevant, J. M., & Richards, F. M. (1992) Biochemistry 31, 1421-1426] for other position 13 analogs. The differences in the free energy and enthalpy of binding between the Met and Nle peptides are 0.6 and 7.9 kcal/mol at 25 degrees C, respectively. These differences are slightly larger than, but comparable to, the differences in the values for the Met/Ile and Met/Leu pairs. The structure of the mutant complex was determined to 1.85 Angstrom resolution and refined to an R-factor of 17.4% The structures of mutant and wild-type complexes are practically identical although the Nle side chain has a significantly higher average B-factor than the corresponding Met side chain. In contrast, the B-factors of the atoms of the cage of residues surrounding position 13 are all somewhat lower in the Nle variant than in the Met wild-type. Thus, the large differences in the binding enthalpy appear to reside entirely in the difference in chemical properties or dynamic behavior of the -S- and -CH2- groups and not in differences in the geometry of the side chains or the internal cavity surface. In addition, a novel method of obtaining protein stability data by means of isothermal titration calorimetry is introduced.

Thickness-dependent fcc-hcp phase transformation in polycrystalline titanium thin films

Polycrystalline Ti thin films are shown to gradually transform from face-centered cubic (fcc) to hexagonal close-packed structure (hcp) with increasing film thickness. Diffraction stress analysis revealed that the fcc phase is formed in a highly compressive hcp matrix (>= 2 GPa), the magnitude of which decreases with increasing film thickness. A correlation between stress and crystallographic texture vis-a-vis the fcc-hcp phase transformation has been established. The total free energy change of the system upon phase transformation calculated using the experimental results shows that the fcc-hcp transformation is theoretically possible in the investigated film thickness regime (144-720 nm) and the hcp structure is stable for films thicker than 720 nm, whereas the fcc structure can be stabilized in Ti films much thinner than 144 nm. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

A thermodynamic criterion for selection of gas compositions for diamond deposition

Isoactivity lines for carbon with respect to diamond as the standard state have been calculated in the ternary system C-H-O at 1223 K to identify the diamond deposition domain. The gas composition is calculated by suppressing the formation of all condensed forms of carbon using the SOLGASMIX free-energy minimization program. Thirty six gas species were included in the calculation. From the gas composition, isoactivity lines are computed using recent data on the Gibbs energy of diamond. Except for activities less than 0.1, the isoactivity lines are almost linear on the C-H-O ternary diagram. Gas compositions which generate activity of diamond ranging from 1 to 100 at 1223 K fall inside a narrow wedge originating from the point representing CO. This wedge is very similar to the revised lens-shaped diamond growth domain identified by Bachman et al., using inputs from experiment. The small difference between the calculated and observed domains may be attributed to variation in the supersaturation required for diamond deposition with gas composition. The diamond solubility in the gas phase along the isoactivity line for a(di)=100 and P=6.7 kPa exhibits a minimum at 1280 K, which is close to the optimum temperature found experimentally. At higher supersaturations, non-diamond forms of carbon, including amorphous varieties, are expected. The results suggest that thermodynamic calculations can be useful for locating diamond growth domains in more complex CVD systems containing halogens, for which very little experimental data is available.

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.

Thermodynamics of water entry in hydrophobic channels of carbon nanotubes

Experiments and computer simulations demonstrate that water spontaneously fills the hydrophobic cavity of a carbon nanotube. To gain a quantitative thermodynamic understanding of this phenomenon, we use the recently developed two phase thermodynamics method to compute translational and rotational entropies of confined water molecules inside single-walled carbon nanotubes and show that the increase in energy of a water molecule inside the nanotube is compensated by the gain in its rotational entropy. The confined water is in equilibrium with the bulk water and the Helmholtz free energy per water molecule of confined water is the same as that in the bulk within the accuracy of the simulation results. A comparison of translational and rotational spectra of water molecules confined in carbon nanotubes with that of bulk water shows significant shifts in the positions of the spectral peaks that are directly related to the tube radius. (C) 2011 American Institute of Physics.

Quantum and non-Markovian effects in the electron transfer reaction dynamics in the Marcus inverted region

Electron transfer reactions in large molecules may often be coupled to both the polar solvent modes and the intramolecular vibrational modes of the molecule. This can give rise to a complex dynamics which may in some systems, like betaine, be controlled more by vibrational rather than by solvent effects. Additionally, a significant contribution from an ultrafast relaxation component in the solvation dynamics may enhance the complexity. To explain the wide range of behavior that has been observed experimentally, Barbara et al. recently proposed that a model of an electron transfer reaction should minimally consist of a low-frequency classical solvent mode (X), a low-frequency vibrational mode (Q), and a high-frequency quantum mode (q) (J. Phys. Chem. 1991, 96, 3728). In the present work, a theoretical study of this model is described. This study generalizes earlier work by including the biphasic solvent response and the dynamics of the low-frequency vibrational mode in the presence of a delocalized, extended reaction zone. A novel Green's function technique has been developed which allowed us to study the non-Markovian dynamics on a multidimensional surface. The contributions from the high-frequency vibrational mode and the ultrafast component in the non-Markovian solvent dynamics are found to be primarily responsible for the dramatic increase in charge transfer rate over the prediction of the classical theories that neglect both these factors. These, along with a large coupling between the reactant and the product states, may combine to render the electron transfer rate both very large and constant over a wide range of solvent relaxation rates. A study on the free energy gap dependence of the electron transfer rate reveals that the rates are sensitive to changes in the quantum frequency particularly when the free energy gap is very large.

Representation of thermodynamic properties of ternary systems and its application to the system silver-gold-copper at 1350 K

The article presents a generalized analytical expression for description of the integral excess Gibbs free energy of mixing of a ternary system. Twelve constants of the equation are assessed by the least mean squares regressional analysis of the experimental integral excess data of the constituent binaries; three ternary parameters are evaluated by a regressional analysis based on the partial experimental data of a component of the ternary system. The assessed values of the ternary parameters describe the nature of the ternary interaction in the system. Activities and isoactivities of the components in the Ag-Au-Cu system at 1350 K are calculated and found to be in good agreement with the experimental data. This analytical treatment is particularly useful to ternary systems where the thermodynamic data are available from different sources.

Role of Water in the Enzymatic Catalysis: Study of ATP plus AMP -> 2ADP Conversion by Adenylate Kinase

The catalytic conversion ATP + AMP -> 2ADP by the enzyme adenylate kinase (ADK) involves the binding of one ATP. molecule to the LID domain and one AMP molecule to the NMP domain. The latter is followed by a. phosphate transfer and then the release of two ADP molecules. We have computed a novel two-dimensional configurational free energy surface (2DCFES), with one reaction coordinate each for the LID and the NMP domain motions, while considering explicit water interactions. Our computed 2DCFES clearly reveals the existence of a stable half-open half-closed (HOHC) intermediate stale of the enzyme. Cycling of the enzyme through the HOHC state reduces the conformational free energy barrier for. the reaction by about 20 kJ/mol. We find that the stability of the HOHC state (missed in all earlier studies with implicit solvent model) is largely because of the increase of specific interactions of the polar amino acid side chains with water, particularly with the arginine and the histidine residues. Free energy surface of the LID domain is rather rugged, which can conveniently slow down LID's conformational motion, thus facilitating a new substrate capture after the product release in the catalytic cycle.

Fluctuation-dissipation relation for nonlinear Langevin equations

It is shown that the fluctuation-dissipation theorem is satisfied by the solutions of a general set of nonlinear Langevin equations with a quadratic free-energy functional (constant susceptibility) and field-dependent kinetic coefficients, provided the kinetic coefficients satisfy the Onsager reciprocal relations for the irreversible terms and the antisymmetry relations for the reversible terms. The analysis employs a perturbation expansion of the nonlinear terms, and a functional integral calculation of the correlation and response functions, and it is shown that the fluctuation-dissipation relation is satisfied at each order in the expansion.

Activities in the spinel solid solution FeXMg1−XAl2O4

Abstract: Activities in the spinel solid solution FexMg1-xAl2O4 saturated with alpha-Al2O3 have been measured for the compositional range 0 < X < 1 between 1100 and 1350 K using a bielectrolyte solid-state galvanic cell, which may be represented as Pt, Fe + FexMg1-xAl2O4 + alpha-Al2O3//(Y2O3)ThO2/ (CaO)ZrO2//Fe + FeAl2O4 + alpha-Al2O3, Pt Activities of ferrous and magnesium aluminates exhibit small negative deviations from Raoult's law. The excess free energy of mixing of the solid solution is a symmetric function of composition and is independent of temperature: Delta G(E) = -1990 X(1 - X J/mol. Theoretical analysis of cation distribution in spinel solid solution also suggests mild negative deviations from ideality. The lattice parameter varies linearly with composition in samples quenched from 1300 K. Phase relations in the FeO-MgO-Al2O3 system at 1300 K are deduced from the results of this study and auxiliary thermodynamic data from the literature. The calculation demonstrates the influence of intracrystalline ion exchange equilibrium between nonequivalent crystallographic sites in the spinel structure on intercrystalline ion exchange equilibrium between the monoxide and spinel solid solutions (tie-lines). The composition dependence of oxygen partial pressure at 1300 K is evaluated for three-phase equilibria involving the solid solutions Fe + FexMg1-xAl2O4 + alpha-Al2O3 and Fe + FeyMg1-yO + FexMg1-xAl2O4. Dependence of X, denoting the composition of the spinel solid solution, on parameter Y, characterizing the composition of the monoxide solid solution with rock salt structure, in phase fields involving the two solid solutions is elucidated. The tie-lines are slightly skewed toward the MgAl2O4 corner.

Relation between integral and partial properties of ternary solutions along different composition trajectories

The relations between partial and integral properties of ternary solutions along composition trajectories suggested by Kohler, Colinet and Jacob, and along an arbitrary path are derived. The chemical potentials of the components are related to the slope of integral free energy by expressions involving the binary compositions generated by the intersections of the composition trajectory with the sides of the ternary triangle. Only along the Kohler composition trajectory it is possible to derive the integral free energy from the variation of the chemical potential of a single component with composition or vice versa. Along all other paths the differential of the integral free energy is related to two chemical potentials. The Gibbs-Duhem integration proposed by Darken for the ternary system uses the Kohler isogram. The relative merits of different limits for integration are discussed.

Representation of thermodynamic properties of quaternary systems using interaction parameters

Integral excess free energy of a quaternary system has been expressed in terms of the MacLaurin infinite series. The series is subjected to appropriate boundary conditions and each of the derivatives correlated to the corresponding interaction coefficients. The derivation of the partial functions involves extensive summation of various infinite series pertaining to the first order and quaternary parameters to remove any truncational error. The thermodynamic consistency of the derived partials has been established based on the Gibbs-Duhem relations. The equations are used to interpret the thermodynamic properties of the Fe-Cr-Ni-N system.

Statistical thermodynamics of lattice models in zeolites: Implications of local versus global mean field interactions

The statistical thermodynamics of adsorption in caged zeolites is developed by treating the zeolite as an ensemble of M identical cages or subsystems. Within each cage adsorption is assumed to occur onto a lattice of n identical sites. Expressions for the average occupancy per cage are obtained by minimizing the Helmholtz free energy in the canonical ensemble subject to the constraints of constant M and constant number of adsorbates N. Adsorbate-adsorbate interactions in the Brag-Williams or mean field approximation are treated in two ways. The local mean field approximation (LMFA) is based on the local cage occupancy and the global mean field approximation (GMFA) is based on the average coverage of the ensemble. The GMFA is shown to be equivalent in formulation to treating the zeolite as a collection of interacting single site subsystems. In contrast, the treatment in the LMFA retains the description of the zeolite as an ensemble of identical cages, whose thermodynamic properties are conveniently derived in the grand canonical ensemble. For a z coordinated lattice within the zeolite cage, with epsilon(aa) as the adsorbate-adsorbate interaction parameter, the comparisons for different values of epsilon(aa)(*)=epsilon(aa)z/2kT, and number of sites per cage, n, illustrate that for -1 0. We compare the isotherms predicted with the LMFA with previous GMFA predictions [K. G. Ayappa, C. R. Kamala, and T. A. Abinandanan, J. Chem. Phys. 110, 8714 (1999)] (which incorporates both the site volume reduction and a coverage-dependent epsilon(aa)) for xenon and methane in zeolite NaA. In all cases the predicted isotherms are very similar, with the exception of a small steplike feature present in the LMFA for xenon at higher coverages. (C) 1999 American Institute of Physics. [S0021-9606(99)70333-8].

Thermodynamic modeling of quaternary Al-Ga-In-As system

The present research describes the modeling of the thermodynamic properties of the liquid Al-Ga-In-As alloys at 1073 and 1173 K, and investigates the solid-liquid equilibria in the systems. The isothermal molar excess free energy function for the liquid alloys is represented in terms of 37 parameters pertaining to six of the constituent binaries, four ternaries and the quaternary interactions in the system. The corresponding solid alloys which consist of AlAs, GaAs and InAs are assumed to be quasi-regular ternary solutions. The solidus and liquidus compositions are calculated at 1073 and 1173 K using the derived values of the partial components for the solid and liquid alloys at equilibrium. They are in good agreement with those of the experimentally determined values available in the literature. (C) 1999 Elsevier Science S.A. All rights reserved.