Path Optimization in Free Energy Calculations

Free energy calculations are a computational method for determining thermodynamic quantities, such as free energies of binding, via simulation.

Currently, due to computational and algorithmic limitations, free energy calculations are limited in scope.

In this work, we propose two methods for improving the efficiency of free energy calculations.

First, we expand the state space of alchemical intermediates, and show that this expansion enables us to calculate free energies along lower variance paths.

We use Q-learning, a reinforcement learning technique, to discover and optimize paths at low computational cost.

Second, we reduce the cost of sampling along a given path by using sequential Monte Carlo samplers.

We develop a new free energy estimator, pCrooks (pairwise Crooks), a variant on the Crooks fluctuation theorem (CFT), which enables decomposition of the variance of the free energy estimate for discrete paths, while retaining beneficial characteristics of CFT.

Combining these two advancements, we show that for some test models, optimal expanded-space paths have a nearly 80% reduction in variance relative to the standard path.

Additionally, our free energy estimator converges at a more consistent rate and on average 1.8 times faster when we enable path searching, even when the cost of path discovery and refinement is considered.

Understanding Co-translational Protein Targeting and Lithium Dendrite Formation through Free Energy Simulations and Coarse-Grained Models

We describe the application of alchemical free energy methods and coarse-grained models to study two key problems: (i) co-translational protein targeting and insertion to direct membrane proteins to the endoplasmic reticulum for proper localization and folding, (ii) lithium dendrite formation during recharging of lithium metal batteries. We show that conformational changes in the signal recognition particle, a central component of the protein targeting machinery, confer additional specificity during the the recognition of signal sequences. We then develop a three-dimensional coarse-grained model to study the long-timescale dynamics of membrane protein integration at the translocon and a framework for the calculation of binding free energies between the ribosome and translocon. Finally, we develop a coarse-grained model to capture the dynamics of lithium deposition and dissolution at the electrode interface with time-dependent voltages to show that pulse plating and reverse pulse plating methods can mitigate dendrite growth.

Dynamics of transformation and fragmentation of composite liquid nano-particles

Recent research on particle size distributions and particle concentrations near a busy road cannot be explained by the conventional mechanisms for particle evolution of combustion aerosols. Specifically they appear to be inadequate to explain the experimental observations of particle transformation and the evolution of the total number concentration. This resulted in the development of a new mechanism based on their thermal fragmentation, for the evolution of combustion aerosol nano-particles. A complex and comprehensive pattern of evolution of combustion aerosols, involving particle fragmentation, was then proposed and justified. In that model it was suggested that thermal fragmentation occurs in aggregates of primary particles each of which contains a solid graphite/carbon core surrounded by volatile molecules bonded to the core by strong covalent bonds. Due to the presence of strong covalent bonds between the core and the volatile (frill) molecules, such primary composite particles can be regarded as solid, despite the presence of significant (possibly, dominant) volatile component. Fragmentation occurs when weak van der Waals forces between such primary particles are overcome by their thermal (Brownian) motion. In this work, the accepted concept of thermal fragmentation is advanced to determine whether fragmentation is likely in liquid composite nano-particles. It has been demonstrated that at least at some stages of evolution, combustion aerosols contain a large number of composite liquid particles containing presumably several components such as water, oil, volatile compounds, and minerals. It is possible that such composite liquid particles may also experience thermal fragmentation and thus contribute to, for example, the evolution of the total number concentration as a function of distance from the source. Therefore, the aim of this project is to examine theoretically the possibility of thermal fragmentation of composite liquid nano-particles consisting of immiscible liquid v components. The specific focus is on ternary systems which include two immiscible liquid droplets surrounded by another medium (e.g., air). The analysis shows that three different structures are possible, the complete encapsulation of one liquid by the other, partial encapsulation of the two liquids in a composite particle, and the two droplets separated from each other. The probability of thermal fragmentation of two coagulated liquid droplets is discussed and examined for different volumes of the immiscible fluids in a composite liquid particle and their surface and interfacial tensions through the determination of the Gibbs free energy difference between the coagulated and fragmented states, and comparison of this energy difference with the typical thermal energy kT. The analysis reveals that fragmentation was found to be much more likely for a partially encapsulated particle than a completely encapsulated particle. In particular, it was found that thermal fragmentation was much more likely when the volume ratio of the two liquid droplets that constitute the composite particle are very different. Conversely, when the two liquid droplets are of similar volumes, the probability of thermal fragmentation is small. It is also demonstrated that the Gibbs free energy difference between the coagulated and fragmented states is not the only important factor determining the probability of thermal fragmentation of composite liquid particles. The second essential factor is the actual structure of the composite particle. It is shown that the probability of thermal fragmentation is also strongly dependent on the distance that each of the liquid droplets should travel to reach the fragmented state. In particular, if this distance is larger than the mean free path for the considered droplets in the air, the probability of thermal fragmentation should be negligible. In particular, it follows form here that fragmentation of the composite particle in the state with complete encapsulation is highly unlikely because of the larger distance that the two droplets must travel in order to separate. The analysis of composite liquid particles with the interfacial parameters that are expected in combustion aerosols demonstrates that thermal fragmentation of these vi particles may occur, and this mechanism may play a role in the evolution of combustion aerosols. Conditions for thermal fragmentation to play a significant role (for aerosol particles other than those from motor vehicle exhaust) are determined and examined theoretically. Conditions for spontaneous transformation between the states of composite particles with complete and partial encapsulation are also examined, demonstrating the possibility of such transformation in combustion aerosols. Indeed it was shown that for some typical components found in aerosols that transformation could take place on time scales less than 20 s. The analysis showed that factors that influenced surface and interfacial tension played an important role in this transformation process. It is suggested that such transformation may, for example, result in a delayed evaporation of composite particles with significant water component, leading to observable effects in evolution of combustion aerosols (including possible local humidity maximums near a source, such as a busy road). The obtained results will be important for further development and understanding of aerosol physics and technologies, including combustion aerosols and their evolution near a source.

Effect of Fe3O4 addition on removal of ammonium by zeolite NaA

Magnetic zeolite NaA with different Fe3O4 loadings was prepared by hydrothermal synthesis based on metakaolin and Fe3O4. The effect of added Fe3O4 on the removal of ammonium by zeolite NaA was investigated by varying the Fe3O4 loading, pH, adsorption temperature, initial concentration, adsorption time. Langmuir, Freundlich, and pseudo-second-order modeling were used to describe the nature and mechanism of ammonium ion exchange using both zeolite and magnetic zeolite. Thermodynamic parameters such as change in Gibbs free energy, enthalpy and entropy were calculated. The results show that all the selected factors affect the ammonium ion exchange by zeolite and magnetic zeolite, however, the added Fe3O4 apparently does not affect the ion exchange performance of zeolite to the ammonium ion. Freundlich model provides a better description of the adsorption process than Langmuir model. Moreover, kinetic analysis indicates the exchange of ammonium on the two materials follows a pseudosecond-order model. Thermodynamic analysis makes it clear that the adsorption process of ammonium is spontaneous and exothermic. Regardless of kinetic or thermodynamic analysis, all the results suggest that no considerable effect on the adsorption of the ammonium ion by zeolite is found after the addition of Fe3O4. According to the results, magnetic zeolite NaA can be used for the removal of ammonium due to the good adsorption performance and easy separation method from aqueous solution.

Adsorption thermodynamics and kinetics study for the self-assembly of 1,8,15,22-tetraaminophthalocyanatocobalt(II) on glassy carbon surface

Spontaneous adsorption of 1,8,15,22-tetraaminophthalocyanatocobalt(II) (4α-CoIITAPc) on glassy carbon (GC) electrode leads to the formation of a stable self-assembled monolayer (SAM). Since the SAM of 4α-CoIITAPc is redox active, its adsorption on GC electrode was followed by cyclic voltammetry. SAM of 4α-CoIITAPc on GC electrode shows two pairs of well-defined redox peaks corresponding to CoIII/CoII and CoIIIPc−1/CoIIIPc−2. The surface coverage (Γ) value, calculated by integrating the charge under CoII oxidation, was used to study the adsorption thermodynamics and kinetics of 4α-CoIITAPc on GC surface. Cyclic voltammetric studies show that the adsorption of 4α-CoIITAPc on GC electrode has reached the saturation coverage (Γs) within 3 h. The Γs value for the SAM of 4α-CoIITAPc on GC electrode was found to be 2.37 × 10−10 mol cm−2. Gibbs free energy (ΔGads) and adsorption rate constant (kad) for the adsorption of 4α-CoIITAPc on GC surface were found to be −16.76 kJ mol−1 and 7.1 M−1 s−1, respectively. The possible mechanism for the self-assembly of 4α-CoIITAPc on GC surface is through the addition of nucleophilic amines to the olefinic bond on the GC surface in addition to a meager contribution from π stacking. The contribution of π stacking was confirmed from the adsorption of unsubstituted phthalocyanatocobalt(II) (CoPc) on GC electrode. Raman spectra for the SAM of 4α-CoIITAPc on carbon surface shows strong stretching and breathing bands of Pc macrocycle, pyrrole ring and isoindole ring. Raman and CV studies suggest that 4α-CoIITAPc is adopting nearly a flat orientation or little bit tilted orientation.

An update on the thermodynamics of Ta2O5

Using a solid-state electrochemical cell incorporating yttria-doped thoria (YDT) as the electrolyte and a mixture of (Mn + MnO) as the reference electrode, standard Gibbs free energy of formation of beta-Ta2O5 has been determined as a function of temperature in the range (1000 to 1300) K. The solid-state electrochemical cell used can be represented as (-)Pt,Ta +Ta2O5//(Y2O3)ThO2//Mn + MnO, Pt(+) Combining the reversible e.m.f. of the cell with recent data on the free energy of formation of MnO, standard Gibbs free energy of formation of Ta2O5 from Ta metal and diatomic oxygen gas (O-2) in the temperature range (1000 to 1300) K is obtained: Delta fG degrees +/- 0.35/(kJ.mol(-1)) = -2004.376 + 0.40445(T/K). Because of the significant solid solubility of oxygen in tantalum, a small correction for the activity of Ta in the metal phase in equilibrium with Ta2O5 is applied. An analysis of the results obtained in this study and other free energy data reported in the literature by the "third law" method suggests the need for refining data for Ta2O5 reported in thermodynamic compilations. Used in the analysis is a revised value for standard entropy of Ta2O5 based on more recent low-temperature heat capacity measurements. An improved set of thermodynamic properties of ditantalum pentoxide (Ta2O5) are presented in the temperature range (298.15 to 2200) K. (C) 2008 Elsevier Ltd. All rights reserved.

Role of entropy in the stability of cobalt titanates

The standard molar Gibbs free energy of formation of Co2TiO4, CoTiO3,and CoTi2O5 as a function of temperature over an extended range (900 to 1675) K was measured using solid-state electrochemical cells incorporating yttria-stabilized zirconia as the electrolyte, with CoO as reference electrode and appropriate working electrodes. For the formation of the three compounds from their component oxides CoO with rock-salt and TiO2 with rutile structure, the Gibbs free energy changes are given by:Delta(f)G degrees((ox))(Co2TiO4) +/- 104/(J . mol(-1)) = -18865 - 4.108 (T/K)Delta(f)G degrees((ox))(CoTiO3) +/- 56/(J . mol(-1)) = -19627 + 2.542 (T/K) Delta(f)G degrees((ox))(CoTi2O5) +/- 52/(J . mol(-1)) = -6223 - 6.933 (T/K) Accurate values for enthalpy and entropy of formation were derived. The compounds Co2TiO4 with spinel structure and CoTi2O5 with pseudo-brookite structure were found to be entropy stabilized. The relatively high entropy of these compounds arises from the mixing of cations on specific crystallographic sites. The stoichiometry of CoTiO3 was confirmed by inert gas fusion analysis for oxygen. Because of partial oxidation of cobalt in air, the composition corresponding to the compound Co2TiO4 falls inside a two-phase field containing the spinet solid solution Co2TiO4-Co3O4 and CoTiO3. The spinel solid solution becomes progressively enriched in Co3O4 with decreasing temperature. (c) 2010 Elsevier Ltd. All rights reserved.

Thermodynamic properties of RhO2

A solid-state electrochemical cell, with yttria-stabilized zirconia as the electrolyte and pure O-2 gas at 0.1 MPa as the reference electrode, has been used to measure the oxygen chemical potential corresponding to the equilibrium between beta-Rh2O3 and RhO2 in the temperature range from 850 to 1050K. Using standard Gibbs energy of formation of beta-Rh2O3 available in the literature and the measured oxygen potential, the standard Gibbs free energy of formation of RhO2 is derived as a function of temperature: Delta G(f)degrees(RhO2)(+/- 71)/J mol(-1) = 238,418 + 179.89T Using an estimated value of Delta C-p degrees; for the formation reaction of RhO2 from its elements, the standard enthalpy of formation, standard entropy and isobaric heat capacity of RhO2 at 298.15 K are evaluated: Delta H-f degrees (298.15 K)(+/- 164)/kJ mol(-1) = -244.94, S degrees (298.15 K)(+/- 3.00)/J mol(-1) K-1 = 45.11 and C-p degrees(298.15 K)(+/- 2.6)1mol(-1) K-1 =64.28. (C) 2010 Elsevier B.V. All rights reserved.

Thermodynamic properties of strontium titanates Sr2TiO4, Sr3Ti2O7, Sr4Ti3O10, and SrTiO3

The chemical potentials of SrO in two-phase fields (TiO2 + SrTiO3) (SrTiO3 + Sr4Ti3O10) (Sr4Ti3O10 + Sr3Ti2O7) and (Sr3Ti2O7 + Sr2TiO4) of the pseudo-binary system (SrO + TiO2) have been measured in the temperature range (900 to 1250) K relative to pure SrO as the reference state using solid-state galvanic cells incorporating single crystal SrF2 as the electrolyte The cells were operated under pure oxygen at ambient pressure The standard Gibbs free energies of formation of strontium titanates SrTiO3 Sr4Ti3O10 Sr3Ti2O7 and Sr2TiO4 from their component binary oxides were derived from the reversible electromotive force (EMF) of the cells For the formation of the four compounds from their component oxides TiO2 with rutile structure and SrO the standard Gibbs free energy changes are given by Delta G((ox))(SrTiO3) +/- 89/(J mol(-1)) = -121878 + 3 881(T/K) Delta G((ox))(Sr4Ti3O10) +/- 284/(J mol(-1)) = -409197 + 14 749(T/K) Delta G((ox))(Sr3Ti2O7) +/- 190/(J mol(-1)) = -285827 + 10 022(T/K) Delta G((ox))(Sr2TiO4) +/- 110/(J mol(-1))= -159385 + 3 770(T/K) The reference state for solid TiO2 is the rutile form The results of this study are in good agreement with Gibbs free energy of formation data reported in the literature for SrTiO3 but differ significantly with data for Sr4Ti3O10 For Si3Ti2O7 and Si2TiO4 experimental measurements are not available in the literature for direct comparison with the results obtained in this study (C) 2010 Elsevier Ltd All rights reserved

Phase relationships in the system Ni-W-O and thermodynamic properties of NiWO,

The phase diagram of the Ni-W-O system at 1200 K was established by metallographic and X-ray identification of the phases present after equilibration at controlled oxygen potentials. The oxygen partial pressures over the samples were fixed by metered streams of CO+CO2 gas mixtures. There was only one ternary oxide, nickel tungstate (NiWO4), in the Ni-W-O system at a total pressure of 1 atm, and this compound decomposed to a mixture of Ni+WO2.72 on lowering the oxygen potential. The Gibbs' free energy of formation of NiWO4 was determined from the measurement of the e.m.f. of the solid oxide galvanic cell, Pt, Ni+NiWO4+WO2.72/CaO-ZrO2/Ni+NiO, Pt and thermodynamic properties of tungsten and nickel oxides available in the literature. For the reaction, NiO(s)+WO3(s)rarrNiWO4(s) DeltaG°=–10500–0.708 T (±250) cal mol–1.

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.

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.

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.

Nitride synthesis using ammonia and hydrazine - a thermodynamic panorama

Because of the wide variety of projected applications of ultrapure nitrides in advanced technologies, there is interest in developing new cost-effective methods of synthesis. Explored in this study is the use of ammonia and hydrazine for the synthesis of nitrides from oxides, sulfides and chlorides. Even when the standard Gibbs energy change for the nitridation reactions involved are moderately positive, the reaction can be made to proceed by lowering the partial pressure of the product gas below its equilibrium value. Use of a metastable form of precursor in the nanometric size range is an alternative method to facilitate nitridation. Ellingham-Richardson-Jeffes diagrams are used for a panoramic presentation of the driving force for each set of reactions as a function of temperature. Oxides are the least promising precursors for nitride synthesis; sulfides offer a larger synthetic window for many useful nitrides such as BN, AlN, InN, VN, TiN, ThN and Si3N4. The standard Gibbs free energy changes for reactions involving chlorides with either ammonia or hydrazine are much more negative. Hydrazine is a more powerful nitriding agent than ammonia. The metastability of hydrazine requires that it be introduced into a reactor through a water-cooled lance. The use of volatile halides with ammonia or hydrazine offers the potential for synthesis of pure and doped nanocrystalline nitrides. Nitride thin films can also be prepared by suitable adaptations of the chloride route. (C) 2002 Kluwer Academic Publishers.

Thermodynamic analysis of enzyme catalysed reactions: new insights into the Michaelis-Menten equation

A simple thermodynamic analysis of the well-known Michaelis-Menten equation (MME) of enzyme catalysis is proposed that employs the chemical potential mu to follow the Gibbs free energy changes attending the formation of the enzyme-substrate complex and its turnover to the product. The main conclusion from the above analysis is that low values of the Michaelis constant KM and high values of the turnover number k(cat) are advantageous: this supports a simple algebraic analysis of the MME, although at variance with current thinking. Available data apparently support the above findings. It is argued that transition state stabilisation - rather than substrate distortion or proximity - is the key to enzyme catalysis.