Measurement of Gibbs energies of formation of CoF2 and MnF2 using a new composite dispersed solid electrolyte

Gibbs energies of formation of CoF2 and MnF2 have been measured in the temperature range from 700 to 1100 K using Al2O3-dispersed CaF2 solid electrolyte and Ni+NiF2 as the reference electrode. The dispersed solid electrolyte has higher conductivity than pure CaF2 thus permitting accurate measurements at lower temperatures. However, to prevent reaction between Al2O3 in the solid electrolyte and NiF2 (or CoF2) at the electrode, the dispersed solid electrolyte was coated with pure CaF2, thus creating a composite structure. The free energies of formation of CoF2 and MnF2 are (± 1700) J mol−1; {fx37-1} The third law analysis gives the enthalpy of formation of solid CoF2 as ΔH° (298·15 K) = −672·69 (± 0·1) kJ mol−1, which compares with a value of −671·5 (± 4) kJ mol−1 given in Janaf tables. For solid MnF2, ΔH°(298·15 K) = − 854·97 (± 0·13) kJ mol−1, which is significantly different from a value of −803·3 kJ mol−1 given in the compilation by Barinet al.

Phase formation and stability of alloy phases in free nanoparticles: some insights

This paper explores phase formation and phase stability in free nanoparticles of binary alloys. A procedure for estimating the size and composition dependent free energies incorporating the contributions from the interfaces has been presented. Both single phase solid solution and two phase morphology containing interphase interfaces have been considered. A free energy scenario has been evaluated for two binary alloy systems Ag-Ni and Ag-Cu to predict the microstructure of the alloy nanoparticles at different size ranges as a function of composition. Both Ag-Cu and Ag-Ni systems exhibit wide bulk immiscibility. Ag-Ni nanoparticles were synthesized using the wet chemical synthesis technique whereas Ag-Cu nanoparticles were synthesized using laser ablation of a Ag-Cu target immersed in distilled water. Microstructural and compositional characterization of Ag-Ni and Ag-Cu nanoparticles on a single nanoparticle level was conducted using transmission electron microscopy. Nanoparticle microstructures observed from the microscopic investigation have been correlated with thermodynamic calculation results. It is shown that the observed two phase microstructure consisting of Ag-Ni solid solution in partial decomposed state coexisting with pure Ag phases in the case of Ag-Ni nanoparticles can be only be rationalized by invoking the tendency for phase separation of an initial solid solution with increase in nanoparticle size. Smaller sized Ag-Ni nanoparticles prefer a single phase solid solution microstructure. Due to an increase in particle size during the synthesis process the initial solid solution decomposes into an ultrafine scale phase separated microstructure. We have shown that it is necessary to invoke critical point phenomenon and wetting transition in systems showing a critical point that leads to phase separated Ag-Ni nanoparticles providing a catalytic substrate for the nucleation of equilibrium Ag over it. In the case of the Ag-Cu system, we report the experimental observation of a core shell structure at small sizes. This can be rationalized in terms of a metastable solid solution. It is argued that the nucleation barrier can prevent the formation of biphasic morphology with an internal interface. In such a situation, demixing of the solid solution can bring the system to a lower energy configuration. This has lead to the observed core-shell morphology in the Ag-Cu system during room temperature synthesis.

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.

Evaluation of the surface free energy of plant surfaces: toward standardizing the procedure

Plant surfaces have been found to have a major chemical and physical heterogeneity and play a key protecting role against multiple stress factors. During the last decade, there is a raising interest in examining plant surface properties for the development of biomimetic materials. Contact angle measurement of different liquids is a common tool for characterizing synthetic materials, which is just beginning to be applied to plant surfaces. However, some studies performed with polymers and other materials showed that for the same surface, different surface free energy values may be obtained depending on the number and nature of the test liquids analyzed, materials' properties, and surface free energy calculation methods employed. For 3 rough and 3 rather smooth plant materials, we calculated their surface free energy using 2 or 3 test liquids and 3 different calculation methods. Regardless of the degree of surface roughness, the methods based on 2 test liquids often led to the under- or over-estimation of surface free energies as compared to the results derived from the 3-Liquids method. Given the major chemical and structural diversity of plant surfaces, it is concluded that 3 different liquids must be considered for characterizing materials of unknown physico-chemical properties, which may significantly differ in terms of polar and dispersive interactions. Since there are just few surface free energy data of plant surfaces with the aim of standardizing the calculation procedure and interpretation of the results among for instance, different species, organs, or phenological states, we suggest the use of 3 liquids and the mean surface tension values provided in this study.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Free energy reconstruction from nonequilibrium single-molecule pulling experiments

Laser tweezers and atomic force microscopes are increasingly used to probe the interactions and mechanical properties of individual molecules. Unfortunately, using such time-dependent perturbations to force rare molecular events also drives the system away from equilibrium. Nevertheless, we show how equilibrium free energy profiles can be extracted rigorously from repeated nonequilibrium force measurements on the basis of an extension of Jarzynski's remarkable identity between free energies and the irreversible work.

Direct measurement of salt-bridge solvation energies using a peptide model system: implications for protein stability.

The solvation energies of salt bridges formed between the terminal carboxyl of the host pentapeptide AcWL- X-LL and the side chains of Arg or Lys in the guest (X) position have been measured. The energies were derived from octanol-to-buffer transfer free energies determined between pH 1 and pH 9. 13C NMR measurements show that the salt bridges form in the octanol phase, but not in the buffer phase, when the side chains and the terminal carboxyl group are charged. The free energy of salt-bridge formation in octanol is approximately -4 kcal/mol (1 cal = 4.184 J), which is equal to or slightly larger than the sum of the solvation energies of noninteracting pairs of charged side chains. This is about one-half the free energy that would result from replacing a charge pair in octanol with a pair of hydrophobic residues of moderate size. Therefore, salt bridging in octanol can change the favorable aqueous solvation energy of a pair of oppositely charged residues to neutral or slightly unfavorable but cannot provide the same free energy decrease as hydrophobic residues. This is consistent with recent computational and experimental studies of protein stability.

Relative basicities of free base porphyrins; understanding the role of macrocyclic distortion

Porphyrins have been the center of numerous investigations in different areas of chemistry, geochemistry, and the life sciences. In nature the conformation of the porphyrin macrocycle varies, depending on the function of its apoenzyme. It is believed that the conformation of the porphyrin ring is necessary for the enzyme to achieve its function and modify its reactivity. It is important to understand how the conformation of the porphyrin ring will influence its properties. ^ In synthetic porphyrins particular conformations and ring deformations can be achieved by peripheral substitution, metallation, core substitution, and core protonation among other alterations of the macrocycle. The macrocyclic distortions will affect the ring current, the ability of pyrroles to intramolecularly hydrogen bond and the relative basicity of each of the porphyrins. To understand these effects different theoretical models are used. The ground state structure of each of 19 free base porphyrins is determined using molecular mechanics (MM+) and semiempirical methods (PM3). The energetics of deformation of the macrocyclic core is calculated by carrying out single point energy calculations for the conformation achieved by each synthetic compound. Enthalpies of solution and enthalpies of protonation of 10 porphyrins with varying degrees of macrocyclic deformation and varying electron withdrawing groups in the periphery are determined using solution calorimetry. Using Hess's Law, the relative basicity of each of the different free base porphyrins is calculated. NMR results are described, including the determination of free energies of activation of ring tautomerization and hydrogen bonding for several compounds. It was found that in the absence of electronic effects, the greater macrocyclic deformation, the greater the basicity of the porphyrins. This basicity is attenuated by the presence of electron withdrawing groups and ability to of the macrocycle to intramolecularly hydrogen bond. ^

Relative basicities of free base porphyrins : understanding the role of macrocyclic distortion

Porphyrins have been the center of numerous investigations in different areas of chemistry, geochemistry, and the life sciences. In nature the conformation of the porphyrin macrocycle varies, depending on the function of its apoenzyme. It is believed that the conformation of the porphyrin ring is necessary for the enzyme to achieve its function and modify its reactivity. It is important to understand how the conformation of the porphyrin ring will influence its properties. In synthetic porphyrins particular conformations and ring deformations can be achieved by peripheral substitution, metallation, core substitution, and core protonation among other alterations of the macrocycle. The macrocyclic distortions will affect the ring current, the ability of pyrroles to intramolecularly hydrogen bond and the relative basicity of each of the porphyrins. To understand these effects different theoretical models are used. The ground state structure of each of 19 free base porphyrins is determined using molecular mechanics (MM+) and semiempirical methods (PM3). The energetics of deformation of the macrocyclic core is calculated by carrying out single point energy calculations for the conformation achieved by each synthetic compound. Enthalpies of solution and enthalpies of protonation of 10 porphyrins with varying degrees of macrocyclic deformation and varying electron withdrawing groups in the periphery are determined using solution calorimetry. Using Hess's Law, the relative basicity of each of the different free base porphyrins is calculated. NMR results are described, including the determination of free energies of activation of ring tautomerization and hydrogen bonding for several compounds. It was found that in the absence of electronic effects, the greater macrocyclic deformation, the greater the basicity of the porphyrins. This basicity is attenuated by the presence of electron withdrawing groups and ability to of the macrocycle to intramolecularly hydrogen bond.

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.

Experimental evidence for competitive N–O and O–C bond homolysis in gasphase alkoxyamines

The extensive use of alkoxyamines in controlled radical polymerisation and polymer stabilisation is based on rapid cycling between the alkoxyamine (R1R2NO–R3) and a stable nitroxyl radical (R1R2NO•) via homolysis of the labile O–C bond. Competing homolysis of the alkoxyamine N–O bond has been predicted to occur for some substituents leading to production of aminyl and alkoxyl radicals. This intrinsic competition between the O–C and N–O bond homolysis processes has to this point been difficult to probe experimentally. Herein we examine the effect of local molecular structure on the competition between N–O and O–C bond cleavage in the gas phase by variable energy tandem mass spectrometry in a triple quadrupole mass spectrometer. A suite of cyclic alkoxyamines with remote carboxylic acid moieties (HOOC–R1R2NO–R3) were synthesised and subjected to negative ion electrospray ionisation to yield [M – H]− anions where the charge is remote from the alkoxyamine moiety. Collision-induced dissociation of these anions yield product ions resulting, almost exclusively, from homolysis of O–C and/or N–O bonds. The relative efficacy of N–O and O–C bond homolysis was examined for alkoxyamines incorporating different R3 substituents by varying the potential difference applied to the collision cell, and comparing dissociation thresholds of each product ion channel. For most R3 substituents, product ions from homolysis of the O–C bond are observed and product ions resulting from cleavage of the N–O bond are minor or absent. A limited number of examples were encountered however, where N–O homolysis is a competitive dissociation pathway because the O–C bond is stabilised by adjacent heteroatom(s) (e.g., R3 = CH2F). The dissociation threshold energies were compared for different alkoxyamine substituents (R3) and the relative ordering of these experimentally determined energies is shown to correlate with the bond dissociation free energies, calculated by ab initio methods. Understanding the structure-dependent relationship between these rival processes will assist in the design and selection of alkoxyamine motifs that selectively promote the desirable O–C homolysis pathway.

Spinel-Corundum Phase Equilibria in the Systems Mn-Cr-Al-O and Co-Cr-Al-O at 1373 K

The tie lines delineating ion-exchange equilibria between MCr2O4-MAl2O4 spinel solid solution, where M is either Mn or Co, and Cr2O3-Al2O3 solid solution with the corundum structure were determined at 1373 K by electron microprobe and E0AX point count analysis of the oxide phases equilibrated with metallic Co and Au-5% Mn. The component activities in the spinel solid solutions are derived from the tie lines and the thernodynamic data for Cr2O3-Al2O3 soiid solutions available hi the literature. The Gibbs free energies of mixing calculated from the experimental data are discussed in relation to the values derived from the cation distribution a.odel based on the site preference energies and assuming random mixing on both tetrahedral and octahedral sites. Positive deviations from ideality observed in this study suggest a miscibility gap for both series of spinel solid solutions at low temperatures in the absence of oxidation.