Estimating Relative Free Energies from a Single Simulation of the Initial State

Author(s):  
Alan E. Mark ◽  
Heiko Schäfer ◽  
Haiyan Liu ◽  
Wilfred van Gunsteren
Keyword(s):  
1986 ◽  
Vol 64 (8) ◽  
pp. 1638-1642 ◽  
Author(s):  
Urmila Mandal ◽  
Kaushik Das ◽  
Kiron Kumar Kundu

Rate constants of acid-catalyzed hydrolysis of sucrose (S) to D-glucose and L-fructose have been determined at 25 °C by optical rotation measurements in aqueous mixtures of protophobic protic glycerol (GL), protophilic protic urea (UH), aprotic dioxane (D), and dipolar aprotic dimethyl sulphoxide (DMSO). Transfer free energies of the substrate sucrose, [Formula: see text] have also been determined in the solvents from solubility measurements. These values as well as those of H+, as obtained earlier by use of the widely used tetraphenylarsonium tetraphenylboron (TATB) reference electrolyte assumption, yielded transfer free energies of the transition state. The observed log (ks/kw) – composition profiles reveal that the rates increase monotonically in GL–water mixtures, that decrease more or less monotonically in UH– and DMSO–water mixtures, and decrease up to 10 mol% D in D–water mixtures, beyond which the values tend to increase. Examination of [Formula: see text]–composition profiles for the different species in each case indicates that the initial and transition state solvation get more or less compensated and the observed rates are dictated by the increased solvation of H+ in aqueous UH, DMSO, and D co-solvent systems. But in GL–water mixtures the decreased solvation of the transition state compared with the initial state is overcome by the decreased solvation of H+, thus resulting in the gradual enhancement of the rates of the reaction. The observed linearity of the correlative plots of −δ(ΔG≠) [= RT ln (ks/kw)] vs. [Formula: see text] with distinctly different slopes in the two cases also substantiates the relative importance of H+ solvation in dictating the rates of the reaction in these widely different aqueous co-solvents.


1986 ◽  
Vol 64 (2) ◽  
pp. 300-307 ◽  
Author(s):  
Urmila Mandal ◽  
Sumita Sen ◽  
Kaushik Das ◽  
Kiron Kumar Kundu

Rate constants (ks) of alkaline fading of crystal violet (CV+) have been determined at 25 °C by spectrophotometric measurements in aqueous mixtures of some protic, aprotic, and dipolar aprotic cosolvents. Transfer free energies of the substrate (CV+), [Formula: see text], were also determined in some of the solvent systems from solubility measurements of the chloride salt, and by subtracting [Formula: see text] obtained earlier by use of the tetraphenylarsonium tetraphenylboron (TATB) extrathermodynamic assumption. This helped determine transfer free energies of the transition state (X≠), [Formula: see text] values of lyate ion (S−) based on the TATB assumption are already known for all of these solvent systems. The observed log (ks/kw) – composition profiles reveal that the relative solvation of the reacting species rather than the dielectric constant of the solvents dictates the complex variation of the rates of the reaction in these solvent systems. Correlation of [Formula: see text] with [Formula: see text] indicates that the reaction is largely controlled by the relative solvation of S− in most of the cases. But analysis of [Formula: see text] – composition profiles for some of the solvent systems reveals that the non-compensation of the [Formula: see text] contributions of initial-state substrate and of the transition-state complex, which may be considered to be an outer-sphere complex [CV+](S−), is also in accord with what is expected from the relative solvating characteristics of the cosolvents as guided by their respective physico-chemical properties.


Author(s):  
F. I. Grace

An interest in NiTi alloys with near stoichiometric composition (55 NiTi) has intensified since they were found to exhibit a unique mechanical shape memory effect at the Naval Ordnance Laboratory some twelve years ago (thus refered to as NITINOL alloys). Since then, the microstructural mechanisms associated with the shape memory effect have been investigated and several interesting engineering applications have appeared.The shape memory effect implies that the alloy deformed from an initial shape will spontaneously return to that initial state upon heating. This behavior is reported to be related to a diffusionless shear transformation which takes place between similar but slightly different CsCl type structures.


2012 ◽  
Vol 8 (1) ◽  
pp. 37-48
Author(s):  
S. Chehaibi ◽  
K. Abrougui ◽  
F. Haouala

The effects of mechanical perforation densities by extracting soil cores through an aerator Vertidrain with a working width of 1.6 m and equipped with hollow tines spaced of 65 mm, were studied on a sandy soil of a grassy sward in the Golf Course El Kantaoui in Sousse (Tunisia). The mechanical aeration was performed at two densities: 250 and 350 holes/m2. The cone penetration resistance and soil water infiltration were measured. These parameters were performed at initial state before aeration (E0) and then on the 10th, 20th and 30th day after aeration. These results showed that perforation density of 350 holes/m2 had a positive effect on the soil by reducing its cone resistance to penetration compared to the initial state (Rp = 14.8 daN/cm2). At 5 cm depth the decrease in resistance to penetration was 34% and 43% on the 10th and 20th day after aeration, respectively. However, on the 30th day after aeration the soil resistance to penetration tended to grow and its value compared to the initial state decreased only by 21 and 26%, respectively, at 5 and 15 cm of depth only by 10% and 9% with 250 holes/m2 density. The soil water infiltration made a good improvement after aeration compared to the initial state. This parameter increased from 4.8 cm/h to 8.3, 10.9 and 13.1 cm/h with 250 holes/m2 density and to 10, 12.9 and 14.8 cm/h with 350 holes/m2 density on the 10th, 20th and 30th day following the aeration.


2014 ◽  
Vol 62 (1) ◽  
pp. 129-137
Author(s):  
A. Sawicki ◽  
J. Mierczyński

Abstract A basic set of experiments for the determination of mechanical properties of sands is described. This includes the determination of basic physical and mechanical properties, as conventionally applied in soil mechanics, as well as some additional experiments, which provide further information on mechanical properties of granular soils. These additional experiments allow for determination of steady state and instability lines, stress-strain relations for isotropic loading and pure shearing, and simple cyclic shearing tests. Unconventional oedometric experiments are also presented. Necessary laboratory equipment is described, which includes a triaxial apparatus equipped with local strain gauges, an oedometer capable of measuring lateral stresses and a simple cyclic shearing apparatus. The above experiments provide additional information on soil’s properties, which is useful in studying the following phenomena: pre-failure deformations of sand including cyclic loading compaction, pore-pressure generation and liquefaction, both static and caused by cyclic loadings, the effect of sand initial anisotropy and various instabilities. An important feature of the experiments described is that they make it possible to determine the initial state of sand, defined as either contractive or dilative. Experimental results for the “Gdynia” model sand are shown.


2020 ◽  
Author(s):  
Zhengqing Tong ◽  
Margaret S. Cheung ◽  
Barry D. Dunietz ◽  
Eitan Geva ◽  
Xiang Sun

The nonequilibrium Fermi’s golden rule (NE-FGR) describes the time-dependent rate coefficient for electronic transitions, when the nuclear degrees of freedom start out in a <i>nonequilibrium</i> state. In this letter, the linearized semiclassical (LSC) approximation of the NE-FGR is used to calculate the photoinduced charge transfer rates in the carotenoid-porphyrin-C<sub>60</sub> molecular triad dissolved in explicit tetrahydrofuran. The initial nonequilibrium state corresponds to impulsive photoexcitation from the equilibrated ground-state to the ππ* state, and the porphyrin-to-C<sub>60</sub> and the carotenoid-to-C<sub>60</sub> charge transfer rates are calculated. Our results show that accounting for the nonequilibrium nature of the initial state significantly enhances the transition rate of the porphyrin-to-C<sub>60</sub> CT process. We also derive the instantaneous Marcus theory (IMT) from LSC NE-FGR, which casts the CT rate coefficients in terms of a Marcus-like expression, with explicitly time-dependent reorganization energy and reaction free energy. IMT is found to reproduce the CT rates in the system under consideration remarkably well.


2020 ◽  
Author(s):  
Lucian Chan ◽  
Garrett Morris ◽  
Geoffrey Hutchison

The calculation of the entropy of flexible molecules can be challenging, since the number of possible conformers grows exponentially with molecule size and many low-energy conformers may be thermally accessible. Different methods have been proposed to approximate the contribution of conformational entropy to the molecular standard entropy, including performing thermochemistry calculations with all possible stable conformations, and developing empirical corrections from experimental data. We have performed conformer sampling on over 120,000 small molecules generating some 12 million conformers, to develop models to predict conformational entropy across a wide range of molecules. Using insight into the nature of conformational disorder, our cross-validated physically-motivated statistical model can outperform common machine learning and deep learning methods, with a mean absolute error ≈4.8 J/mol•K, or under 0.4 kcal/mol at 300 K. Beyond predicting molecular entropies and free energies, the model implies a high degree of correlation between torsions in most molecules, often as- sumed to be independent. While individual dihedral rotations may have low energetic barriers, the shape and chemical functionality of most molecules necessarily correlate their torsional degrees of freedom, and hence restrict the number of low-energy conformations immensely. Our simple models capture these correlations, and advance our understanding of small molecule conformational entropy.


2018 ◽  
Author(s):  
Maximiliano Riquelme ◽  
Alejandro Lara ◽  
David L. Mobley ◽  
Toon Vestraelen ◽  
Adelio R Matamala ◽  
...  

<div>Computer simulations of bio-molecular systems often use force fields, which are combinations of simple empirical atom-based functions to describe the molecular interactions. Even though polarizable force fields give a more detailed description of intermolecular interactions, nonpolarizable force fields, developed several decades ago, are often still preferred because of their reduced computation cost. Electrostatic interactions play a major role in bio-molecular systems and are therein described by atomic point charges.</div><div>In this work, we address the performance of different atomic charges to reproduce experimental hydration free energies in the FreeSolv database in combination with the GAFF force field. Atomic charges were calculated by two atoms-in-molecules approaches, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS). To account for polarization effects, the charges were derived from the solute's electron density computed with an implicit solvent model and the energy required to polarize the solute was added to the free energy cycle. The calculated hydration free energies were analyzed with an error model, revealing systematic errors associated with specific functional groups or chemical elements. The best agreement with the experimental data is observed for the MBIS atomic charge method, including the solvent polarization, with a root mean square error of 2.0 kcal mol<sup>-1</sup> for the 613 organic molecules studied. The largest deviation was observed for phosphor-containing molecules and the molecules with amide, ester and amine functional groups.</div>


Sign in / Sign up

Export Citation Format

Share Document