First-Principles Atomistic Thermodynamics and Configurational Entropy
Ontology highlight
ABSTRACT:
SUBMITTER: Sutton C
PROVIDER: S-EPMC7793851 | biostudies-literature | 2020 Jan
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Understanding the structural tendencies of nanoconfined water is of great interest for nanoscience and biology, where nano/micro-sized objects may be separated by very few layers of water. Here we investigate the properties of ice confined to a quasi-2D monolayer by a featureless, chemically neutral potential, in order to characterize its intrinsic behaviour. We use density-functional theory simulations with a non-local van der Waals density functional. An ab initio random structure search reveals all the energetically competitive monolayer configurations to belong to only two of the previously-identified families, characterized by a square or honeycomb hydrogen-bonding network, respectively. We discuss the modified ice rules needed for each network, and propose a simple point dipole 2D lattice model that successfully explains the energetics of the square configurations. All identified stable phases for both networks are found to be non-polar (but with a topologically non-trivial texture for the square) and, hence, non-ferroelectric, in contrast to previous predictions from a five-site empirical force-field model. Our results are in good agreement with very recently reported experimental observations.
| S-EPMC4700474 | biostudies-literature
Project description:Recently, the design of lightweight high entropy alloys (HEAs) with a mass density lower than 5 g/cm3 has attracted much research interest in structural materials. We applied a first principles-based high-throughput method to design lightweight HEAs in single solid-solution phase. Three lightweight quinary HEA families were studied: AlBeMgTiLi, AlBeMgTiSi and AlBeMgTiCu. By comprehensively exploring their entire compositional spaces, we identified the most promising compositions according to the following design criteria: the highest stability, lowest mass density, largest elastic modulus and specific stiffness, along with highest Pugh's ratio. We found that HEAs with the topmost compositions exhibit a negative formation energy, a low density and high specific Young's modulus, but a low Pugh's ratio. Importantly, we show that the most stable composition, Al0.31Be0.15Mg0.14Ti0.05Si0.35 is energetically more stable than its metallic compounds and it significantly outperforms the current lightweight engineering alloys such as the 7075 Al alloy. These results suggest that the designed lightweight HEAs can be energetically more stable, lighter, and stiffer but slightly less ductile compared to existing Al alloys. Similar conclusions can be also drawn for the AlBeMgTiLi and AlBeMgTiCu. Our design methodology and findings serve as a valuable tool and guidance for the experimental development of lightweight HEAs.
| S-EPMC10728166 | biostudies-literature
Project description:The restriction of a small molecule's motion on binding to a protein causes a loss of configurational entropy, and thus a penalty in binding affinity. Some energy models used in computer-aided ligand design neglect this entropic penalty, whereas others account for it based on an expected drop in the number of accessible rotamers upon binding. However, the validity of the physical assumptions underlying the various approaches is largely unexamined. The present study addresses this issue by using Mining Minima calculations to analyze the association of amprenavir with HIV protease. The computed loss in ligand configurational entropy is large, contributing approximately 25 kcal/mol (4.184 kJ/kcal) to DeltaG degrees. Most of this loss results from narrower energy wells in the bound state, rather than a drop in the number of accessible rotamers. Coupling among rotation/translation and internal degrees of freedom complicates the decomposition of the entropy change into additive terms. The results highlight the potential to gain affinity by designing conformationally restricted ligands and have implications for the formulation of energy models for ligand scoring.
| S-EPMC1780070 | biostudies-literature
Project description:Entropy and the second law of thermodynamics are fundamental concepts that underlie all natural processes and patterns. Recent research has shown how the entropy of a landscape mosaic can be calculated using the Boltzmann equation, with the entropy of a lattice mosaic equal to the logarithm of the number of ways a lattice with a given dimensionality and number of classes can be arranged to produce the same total amount of edge between cells of different classes. However, that work seemed to also suggest that the feasibility of applying this method to real landscapes was limited due to intractably large numbers of possible arrangements of raster cells in large landscapes. Here I extend that work by showing that: (1) the proportion of arrangements rather than the number with a given amount of edge length provides a means to calculate unbiased relative configurational entropy, obviating the need to compute all possible configurations of a landscape lattice; (2) the edge lengths of randomized landscape mosaics are normally distributed, following the central limit theorem; and (3) given this normal distribution it is possible to fit parametric probability density functions to estimate the expected proportion of randomized configurations that have any given edge length, enabling the calculation of configurational entropy on any landscape regardless of size or number of classes. I evaluate the boundary limits (4) for this normal approximation for small landscapes with a small proportion of a minority class and show it holds under all realistic landscape conditions. I further (5) demonstrate that this relationship holds for a sample of real landscapes that vary in size, patch richness, and evenness of area in each cover type, and (6) I show that the mean and standard deviation of the normally distributed edge lengths can be predicted nearly perfectly as a function of the size, patch richness and diversity of a landscape. Finally, (7) I show that the configurational entropy of a landscape is highly related to the dimensionality of the landscape, the number of cover classes, the evenness of landscape composition across classes, and landscape heterogeneity. These advances provide a means for researchers to directly estimate the frequency distribution of all possible macrostates of any observed landscape, and then directly calculate the relative configurational entropy of the observed macrostate, and to understand the ecological meaning of different amounts of configurational entropy. These advances enable scientists to take configurational entropy from a concept to an applied tool to measure and compare the disorder of real landscapes with an objective and unbiased measure based on entropy and the second law.
| S-EPMC7512816 | biostudies-literature
Project description:Although geophysical observations of mantle regions that suggest the presence of partial melt have often been interpreted in light of the properties of basaltic liquids erupted at the surface, the seismic and rheological consequences of partial melting in the upper mantle depend instead on the properties of interstitial basaltic melt at elevated pressure. In particular, basaltic melts and glasses display anomalous mechanical softening upon compression up to several GPa, suggesting that the relevant properties of melt are strongly pressure-dependent. A full understanding of such a softening requires study, under compression, of the atomic structure of primitive small-degree basaltic melts at their formation depth, which has proven to be difficult. Here we report multiNMR spectra for a simplified basaltic glass quenched at pressures up to 5 GPa (corresponding to depths down to ∼150 km). These data allow quantification of short-range structural parameters such as the populations of coordination numbers of Al and Si cations and the cation pairs bonded to oxygen atoms. In the model basaltic glass, the fraction of [5,6]Al is ∼40% at 5 GPa and decreases to ∼3% at 1 atm. The estimated fraction of nonbridging oxygens at 5 GPa is ∼84% of that at ambient pressure. Together with data on variable glass compositions at 1 atm, these results allow us to quantify how such structural changes increase the configurational entropy of melts with increasing density. We explore how configurational entropy can be used to explain the anomalous mechanical softening of basaltic melts and glasses.
| S-EPMC7486762 | biostudies-literature
Project description:Configurational entropy is an important parameter in amorphous systems. It is involved in the thermodynamic considerations, plays an important role in the molecular mobility calculations through its appearance in the Adam-Gibbs equation and provides information on the solubility increase of an amorphous form compared to its crystalline counterpart. This paper presents a calorimetric method which enables the scientist to quickly determine the values for the configurational entropy at any temperature and obtain the maximum of information from these measurements.
| S-EPMC3986718 | biostudies-literature
Project description:We present a procedure to determine temperature-dependent thermodynamic properties of crystalline materials from density functional theory molecular dynamics (DFT-MD). Finite temperature properties (structural, thermal, and mechanical properties) of the phases (ground state monoclinic B33, martensitic B19', and austenitic B2) of the shape memory alloy NiTi are investigated. Fluctuation formulas and numerical derivatives are used to evaluate mechanical and thermal properties. A modified version of thermodynamic upsampling is used to enable simulations with lower DFT convergence thresholds. In addition, a thermodynamic integration expression is developed to compute free energies from isobaric DFT-MD simulations that accounts for volume changes. Structural parameters, elastic constants, volume expansion, and specific heats as a function of temperature are evaluated. Phase transitions between B2 and B19' and between B19' and B33 are characterized according to their thermal energy, entropy, and free energy differences as well as their latent heats. Anharmonic effects are shown to play a large role in both stabilizing the austenite B2 phase as well as suppressing the martensitic phase transition. The quasiharmonic approximation to the free energy results in large errors in estimating the martensitic transition temperature by neglecting these large anharmonic components.
| S-EPMC6476555 | biostudies-literature
Project description:We use first-principles density functional theory total energy and linear response phonon calculations to compute the Helmholtz and Gibbs free energy as a function of temperature, pressure, and cell volume in the flexible metal-organic framework material MIL-53(Cr) within the quasiharmonic approximation. GGA and metaGGA calculations were performed, each including empirical van der Waals (vdW) forces under the D2, D3, or D3(BJ) parameterizations. At all temperatures up to 500 K and pressures from -30 MPa to 30 MPa, two minima in the free energy versus volume are found, corresponding to the narrow pore (np) and large pore (lp) structures. Critical positive and negative pressures are identified, beyond which there is only one free energy minimum. While all results overestimated the stability of the np phase relative to the lp phase, the best overall agreement with experiment is found for the metaGGA PBEsol+RTPSS+U+J approach with D3 or D3(BJ) vdW forces. For these parameterizations, the calculated free energy barrier for the np-lp transition is only 3 to 6 kJ per mole of Cr4(OH)4(C8H4O4)4.
| S-EPMC5497523 | biostudies-literature
Project description:The thermodynamics of disordering in Cu3Au have been investigated by measuring the heat capacity of samples with different degrees of long- and short-range order between T = 5 and 720 K using relaxation and differential scanning calorimetry. The heat capacities of L12-ordered and fcc-disordered samples show similar behaviour at low temperatures (<300 K). They deviate positively from the linear combination of the end-member heat capacities between ∼30 and 160 K. However, small differences between the two samples exist, as the disordered sample has a larger heat capacity producing a vibrational entropy of disordering of ∼0.05 R. At temperatures higher than 300 K, the heat capacity of the ordered sample shows a prominent lambda-type anomaly at 675 K due to the diffusive L12-fcc phase transition. When starting these measurements with disordered samples, ordering effects are observed between 400 and 620 K, and the disordering reaction is observed at 660 K. Evaluation of the data gives an enthalpy and entropy of disordering at 683 K of 2.0 kJ mol-1 and 0.39 R, respectively. However, these values increase with increasing temperature, thereby reducing the short-range order. Because the vibrational and configurational disordering effects become active at different temperature regimes, i.e., the vibrational effects at low temperatures (T ≪ 300 K) and the sum of both effects at higher temperatures (T > 300 K), they have been successfully separated.
| S-EPMC4394139 | biostudies-literature
Project description:Advances in simulation techniques and computing hardware have created a substantial overlap between the timescales accessible to atomic-level simulations and those on which the fastest-folding proteins fold. Here we demonstrate, using simulations of four variants of the human villin headpiece, how simulations of spontaneous folding and unfolding can provide direct access to thermodynamic and kinetic quantities such as folding rates, free energies, folding enthalpies, heat capacities, ?-values, and temperature-jump relaxation profiles. The quantitative comparison of simulation results with various forms of experimental data probing different aspects of the folding process can facilitate robust assessment of the accuracy of the calculations while providing a detailed structural interpretation for the experimental observations. In the example studied here, the analysis of folding rates, ?-values, and folding pathways provides support for the notion that a norleucine double mutant of villin folds five times faster than the wild-type sequence, but following a slightly different pathway. This work showcases how computer simulation has now developed into a mature tool for the quantitative computational study of protein folding and dynamics that can provide a valuable complement to experimental techniques.
| S-EPMC3497772 | biostudies-literature