Shear alignment of a disordered lamellar mesophase

The shear alignment of an initially disordered lamellar phase is examined using lattice Boltzmann simulations of a mesoscopic model based on a free-energy functional for the concentration modulation. For a small shear cell of width 8 lambda, the qualitative features of the alignment process are strongly dependent on the Schmidt number Sc = nu/D (ratio of kinematic viscosity and mass diffusion coefficient). Here, lambda is the wavelength of the concentration modulation. At low Schmidt number, it is found that there is a significant initial increase in the viscosity, coinciding with the alignment of layers along the extensional axis, followed by a decrease at long times due to the alignment along the flow direction. At high Schmidt number, alignment takes place due to the breakage and reformation of layers because diffusion is slow compared to shear deformation; this results in faster alignment. The system size has a strong effect on the alignment process; perfect alignment takes place for a small systems of width 8 lambda and 16 lambda, while a larger system of width 32 lambda does not align completely even at long times. In the larger system, there appears to be a dynamical steady state in which the layers are not perfectly aligned-where there is a balance between the annealing of defects due to shear and the creation due to an instability of the aligned lamellar phase under shear. We observe two types of defect creation mechanisms: the buckling instability under dilation, which was reported earlier, as well as a second mechanism due to layer compression.

Metastable extension of the fluorite phase field in Y2O3---ZrO2 and its effect on grain growth

Pure Y2O3 and Y2O3---ZrO2 solid solutions have been prepared by melt atomization and by pyrolysis of nitrate solutions. Extended solubility is readily achieved in both techniques for the entire composition range investigated: melts with 0–30% ZrO2 and precursors with 0–50% ZrO2. However, solidification of under cooled droplets yields almost exclusively single phase powders with the structure of cubic yttria (D53). In contrast, the pyrolysis route leads to a sequence of metastable microstructures beginning with a nanocrystalline disordered fluorite-based (C1) solid solution. Further heating leads to the evolution of much larger (micron size) flake crystals with a {001} texture, concurrent with partial ordering of the oxygen ions to the sites occupied in the D53 structure. The driving force for ordering and the rate of grain growth decrease with increasing ZrO2 addition. Abrupt heating to high temperatures or electron irradiation can induce ordering without substantial grain growth. There is no significant reduction in porosity during the recrystallization, which with the other observations suggests that grain growth is driven by the free energy available for the ordering transformation from fluorite to the yttria structure. This route offers opportunities for single crystal thin film development at relatively low processing temperatures.

High-temperature stability of Ca2ZrSi4O12

The standard Gibbs free energy of formation of orthorhombic Ca2ZrSi4O12 from component oxides ZrO2 (monoclinic), CaO (rock salt), and SiO2 (quartz) has been determined in the temperature range 973 to 1273 K using a solid-state cell incorporating single-crystal CaF2 as the electrolyte: Delta G(f) degrees = -219930 + 11.77T (+/- 1500) J.mol(-1) This is the only quantitative information now available on the stability of Ca2ZrSi4O12.

Metastable extension of the fluorite phase field in Y2O3---ZrO2 and its effect on grain growth

Pure Y2O3 and Y2O3-ZrO2 solid solutions have been prepared by melt atomization and by pyrolysis of nitrate solutions. Extended solubility is readily achieved in both techniques for the entire composition range investigated: melts with 0-30% ZrO2 and precursors with 0-50% ZrO2. However, solidification of under cooled droplets yields almost exclusively single phase powders with the structure of cubic yttria (D5(3)). In contrast, the pyrolysis route leads to a sequence of metastable microstructures beginning with a nanocrystalline disordered fluorite-based (C1) solid solution. Further heating leads to the evolution of much larger (micron size) flake crystals with a {001} texture, concurrent with partial ordering of the oxygen ions to the sites occupied in the D5(3) structure. The driving force for ordering and the rate of grain growth decrease with increasing ZrO2 addition. Abrupt heating to high temperatures or electron irradiation can induce ordering without substantial grain growth. There is no significant reduction in porosity during the recrystallization, which with the other observations suggests that grain growth is driven by the free energy available for the ordering transformation from fluorite to the yttria structure. This route offers opportunities for single crystal thin film development at relatively low processing temperatures.

Phase relations in the system Cu─Ho─O and stability of Cu2Ho2O5

The phase relations in the system Cu-Ho-O have been determined at 1300 K using X-ray diffraction, optical microscopy, and electron microprobe analysis of samples equilibrated in evacuated quartz ampules and in pure oxygen. Only one ternary compound, Cu2Ho2O5, was found to be stable. The Gibbs free energy of formation of this compound has been measured using the solid-state cell Pt,Cu2O + Cu2Ho2O5 + Ho2O3/(Y2O3)ZrO2/CuO + Cu2O,Pt in the temperature range of 973 to 1350 K. For the formation of Cu2Ho2O5 from its binary component oxides, 2CuO(s) + Ho2O3(S) --> Cu2Ho2O5(s) DELTAG-degrees = 11190 - 13.8T(+/- 120) J-mol-1 Since the formation is endothermic, CU2Ho2O5 becomes thermodynamically unstable with respect to CuO and Ho2O3 below 810 K. When the oxygen partial pressure over Cu2Ho2O5 is lowered, it decomposes according to the reaction 2Cu2Ho2O5(s) --> 2Ho2O3(s) + 2Cu2O(S) + O2(g) for which the equilibrium oxygen potential is given by DELTAmu(O2) = - 238510 + 160.2T(+/- 450) J.mol-1 The decomposition temperature at an oxygen partial pressure of 1.52 x 10(4) Pa was measured using a combined DTA-TGA apparatus. Based on these results, an oxygen potential diagram for the system Cu-Ho-O at 1300 K is presented.

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.

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.