Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamic entropy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Metastable sound speed in gas-liquid mixtures

A new method of calculating speed of sound for two-phase flow is presented. The new equation assumes no phase change during the propagation of an acoustic disturbance and assumes that only the total entropy of the mixture remains constant during the process. The new equation predicts single-phase values for the speed of sound in the limit of all gas or all liquid and agrees with available two-phase, air-water sound speed data. Other expressions used in the two-phase flow literature for calculating two-phase, metastable sound speed are reviewed and discussed. Comparisons are made between the new expression and several of the previous expressions -- most notably a triply isentropic equation as used, a triply isentropic equation as used, among others, by Karplus and by Wallis. Appropriate differences are pointed out and a thermodynamic criterion is derived which must be satisfied in order for the triply isentropic expression to be thermodynamically consistent. This criterion is not satisfied for the cases examined, which included two-phase nitrogen, air-water, two-phase parahydrogen, and steam-water. Consequently, the new equation derived is found to be superior to the other equations reviewed.

Bursik, J. W.↗
NASA NTRS1975

GASP: A computer code for calculating the thermodynamic and transport properties for ten fluids: Parahydrogen, helium, neon, methane, nitrogen, carbon monoxide, oxygen, fluorine, argon, and carbon dioxide

A FORTRAN IV subprogram called GASP is discussed which calculates the thermodynamic and transport properties for 10 pure fluids: parahydrogen, helium, neon, methane, nitrogen, carbon monoxide, oxygen, fluorine, argon, and carbon dioxide. The pressure range is generally from 0.1 to 400 atmospheres (to 100 atm for helium and to 1000 atm for hydrogen). The temperature ranges are from the triple point to 300 K for neon; to 500 K for carbon monoxide, oxygen, and fluorine; to 600 K for methane and nitrogen; to 1000 K for argon and carbon dioxide; to 2000 K for hydrogen; and from 6 to 500 K for helium. GASP accepts any two of pressure, temperature and density as input conditions along with pressure, and either entropy or enthalpy. The properties available in any combination as output include temperature, density, pressure, entropy, enthalpy, specific heats, sonic velocity, viscosity, thermal conductivity, and surface tension. The subprogram design is modular so that the user can choose only those subroutines necessary to the calculations.

Hendricks, R. C.↗
DOE OSTI2025

Hybrid Doping Strategy with High‐Entropy Cu/Fe Surface Modification and Zr Bulk Incorporation for Ni‐Rich Cathodes

A hybrid doping strategy combining Zr 4+ bulk doping with high-entropy Cu 2+ /Fe 3+ surface doping is developed to enhance the structural and interfacial stability of Ni-rich layered oxide cathodes. Cu and Fe are selectively introduced at the particle surface via a surface-selective ion-exchange process, forming a ≈15 nm Fe-rich layer while preserving the layered framework. Compared to the pristine cathode, the hybrid sample exhibits significantly improved electrochemical performance in both half-cell and full-cell configurations. In half-cells, the hybrid retains 88.5% and 90.2% after 100 cycles at 1C under 4.6 and 4.5 V, respectively. During high-voltage full-cell cycling, the hybrid cathode maintains over 80% capacity retention, whereas the pristine counterpart retains less than 10% under identical conditions over the same cycling period. XPS, EELS, and DEMS analyses confirm improved oxygen retention, suppressed gas evolution, and stable surface chemistry, while DFT calculations indicate enhanced Me–O bonding in the selected Fe 0.75 Cu 0.25 (Mn 1/16 Co 2/16 Ni 13/16 )O 2 surface composition, which is identified through DFT-calculated mixing energy reaching a minimum at this ratio, indicating the most thermodynamically favorable configuration. In conclusion, these results demonstrate the effectiveness of this hybrid doping strategy in mitigating coupled degradation pathways in Ni-rich cathodes.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Characterization of Interactions between CTX-M-15 and Clavulanic Acid, Desfuroylceftiofur, Ceftiofur, Ampicillin, and Nitrocefin

Cefotaximase-Munich (CTX-M) extended-spectrum beta-lactamases (ESBLs) are commonly associated with Gram-negative, hospital-acquired infections worldwide. Several beta-lactamase inhibitors, such as clavulanate, are used to inhibit the activity of these enzymes. To understand the mechanism of CTX-M-15 activity, we have determined the crystal structures of CTX-M-15 in complex with two specific classes of beta-lactam compounds, desfuroylceftiofur (DFC) and ampicillin, and an inhibitor, clavulanic acid. The crystal structures revealed that Ser70 and five other residues (Lys73, Tyr105, Glu166, Ser130, and Ser237) participate in catalysis and binding of those compounds. Based on analysis of steady-state kinetics, thermodynamic data, and molecular docking to both wild-type and S70A mutant structures, we determined that CTX-M-15 has a similar affinity for all beta-lactam compounds (ceftiofur, nitrocefin, DFC, and ampicillin), but with lower affinity for clavulanic acid. A catalytic mechanism for tested β-lactams and two-step inhibition mechanism of clavulanic acid were proposed. CTX-M-15 showed a higher activity toward DFC and nitrocefin, but significantly lower activity toward ampicillin and ceftiofur. The interaction between CTX-M-15 and both ampicillin and ceftiofur displayed a higher entropic but lower enthalpic effect, compared with DFC and nitrocefin. DFC, a metabolite of ceftiofur, displayed lower entropy and higher enthalpy than ceftiofur. This finding suggests that compounds containing amine moiety (e.g., ampicillin) and the furfural moiety (e.g., ceftiofur) could hinder the hydrolytic activity of CTX-M-15.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Towards absolutely stable ergodicity breaking in two and three dimensions

We propose physically reasonable systems capable of avoiding ergodicity at infinite time in the thermodynamic limit, even with generic perturbations and when coupled to a heat bath. In two dimensions, the rainbow loop soup has (stretched) exponentially numerous absolutely stable nonergodic states with diverging energy but vanishing energy density. In three dimensions the rainbow membrane soup has (stretched) exponentially numerous nonergodic states with diverging energy barriers, leading to infinite-time robust ergodicity breaking that even survives coupling to a nonzero temperature heat bath. We describe our results in the language of exact emergent symmetries and demonstrate how the systems avoid common instabilities. Furthermore, our construction naturally connects to quantum dimer models, topologically ordered systems, the group word construction, and Hamiltonians whose low-energy eigenstates exhibit anomalous entanglement entropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The u 3 si 2 -H system

U 3 Si 2 is of interest to the nuclear industry as a candidate fuel material due to its high uranium density and high thermal conductivity. However, it has been observed to react with hydrogen, resulting in material decrepitation. As a result, it is important to understand the thermodynamics of the U 3 Si 2 -H system. In this study, the thermodynamics of the hydrogen absorption reaction of USi were determined experimentally using Sievert’s gas absorption and related to crystallographic evolution with hydrogen content using X-ray diffraction. Experimentally-determined thermodynamic parameters were compared with results from density functional theory modeling. Results from this study were also compared with those determined in previous work. Sievert’s gas absorption results were used to develop the pressure-composition-temperature (PCT) curves of the U 3 Si 2 -H system. It was found that the hydride phase exhibited a maximum stoichiometry between U 3 Si 2 H 1.8 and U 3 Si 2 H 2 . The two-phase region for hydride formation from U 3 Si 2 exhibited a miscibility gap with a critical temperature between 623 and 673 K, as calculated from the PCT curves. Analysis of the PCT curves also showed that both the enthalpy and entropy of the hydrogen absorption reaction increased with hydrogen content but were lower than the values for uranium trihydride formation from uranium metal. The enthalpy of reaction for hydrogen absorption was calculated to range between -86.9 and -94.8 kJ mol -1 , while the entropy of reaction was calculated to range between 101.9 and 138.8 J mol -1 K -1 . Furthermore, DFT modeling of the thermoydnamic stability of the USi hydride phases yielded a decomposition temperature of U 3 Si 2 H 2 of approximately 600 K, which was consistent with the experimental results. Similarly, the DFT-calculated enthalpy and entropy of reaction to form USiH were determined to be -106.5kJ mol and 121.8J mol -1 K -1 , respectively, which were both in close agreement with the experimentally-determined values.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Denoising diffusion probabilistic models for generative alloy design

Inverse material design is an extremely challenging optimization task made difficult by, in part, the highly nonlinear relationship linking performance with composition. Quantitative approaches have improved significantly owing to advances in high throughput experimentation and computational thermodynamics. However, existing physics-based tools are mostly forward models; input a chemistry and obtain a prediction. More recently the materials community has leveraged advances in the machine learning community to establish novel inverse design frameworks. Very recently denoising diffusion probabilistic models have been shown to be extremely powerful generators producing synthetic data of various modalities e.g. images, text, audio, tables, etc.. In this work a novel framework for alloy design and optimization is proposed leveraging these class of models. Five key generative tasks are demonstrated (1) unconditional generation (2) composition conditioned generation (3) property conditioned generation (4) multi-feedstock conditioned generation and (5) generative optimization. These methods were tested on three case studies: high entropy alloy design, superalloy binder jet additive manufacturing, and in-situ dual-feedstock wire-arc additive manufacturing. Results indicate that the established models are extremely flexible, expressive, and robust. The architecture’s flexibility and training procedure empower the model to learn complex intra-compositional and composition-property relationships. Furthermore, the probabilistic nature of these models makes them well suited for addressing solution non-uniqueness and tackling uncertainty quantification tasks. While the fidelity and quantity of the underlying training data is paramount, we envision that future alloy design frameworks will make extensive use of these kinds of machine learning models as “search” tools bolstering the utility of experimental and computational approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model

We report molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. The SLIC/CDC multiscale model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). The SLIC/CDC model predicts Gibbs energies of solvation for a database of 500 solutes in water with average accuracy better than 1 kcal/mol both for experimental measurements and for explicit-solvent molecular dynamics simulations. The separate SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent accessible surface area (SASA) nonpolar energy mode. The SLIC/SASA model predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems and better than 1.6 kcal/mol average accuracy in ionic liquids. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Accurate prediction of oxygen vacancy concentration with disordered A-site cations in high-entropy perovskite oxides

Abstract Entropic stabilized ABO 3 perovskite oxides promise many applications, including the two-step solar thermochemical hydrogen (STCH) production. Using binary and quaternary A-site mixed {A}FeO 3 as a model system, we reveal that as more cation types, especially above four, are mixed on the A-site, the cell lattice becomes more cubic-like but the local Fe–O octahedrons are more distorted. By comparing four different Density Functional Theory-informed statistical models with experiments, we show that the oxygen vacancy formation energies ( $${E}_{V}^{f}$$ E V f ) distribution and the vacancy interactions must be considered to predict the oxygen non-stoichiometry ( δ ) accurately. For STCH applications, the $${E}_{V}^{f}$$ E V f distribution, including both the average and the spread, can be optimized jointly to improve Δ δ (difference of δ between the two-step conditions) in some hydrogen production levels. This model can be used to predict the range of water splitting that can be thermodynamically improved by mixing cations in {A}FeO 3 perovskites.

08 HYDROGEN↗
DOE OSTI2023

Probing supersymmetric black holes with surface defects

It has long been conjectured that the large N deconfinement phase transition of $\mathcal{N}$ = 4 SU(N) super-Yang-Mills corresponds via AdS/CFT to the Hawking-Page transition in which black holes dominate the thermal ensemble, and quantitative evidence of this has come through the recent matching of the superconformal index of $\frac{1}{16}$-BPS states to the supersymmetric black hole entropy. Here we introduce the half-BPS Gukov-Witten surface defect as a probe of the superconformal index, which also serves as an order parameter for the deconfinement transition. This can be studied directly in field theory as a modification of the usual unitary matrix model or in the dual description as a D3-brane probe in the background of a (complex) supersymmetric black hole. Using a saddle point approximation, we determine our defect index in the large N limit as a simple function of the chemical potentials and show independently that it is reproduced by the renormalized action of the brane in the black hole background. Along the way, we also comment on the Cardy limit and the thermodynamics of the D3-brane in the generalized ensemble. The defect index sharply distinguishes between the confining and the deconfining phases of the gauge theory and thus is a supersymmetric non-perturbative order parameter for these large N phase transitions which deserves further investigation. Finally, our work provides an example where the properties of a black hole coupled to an external system can be analyzed precisely.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Differentiating Hund from Mott physics in a three-band Hubbard-Hund model: Temperature dependence of spectral, transport, and thermodynamic properties

We study the interplay between Mott physics, driven by Coulomb repulsion U, and Hund physics, driven by Hund's coupling J, for a minimal model for Hund metals, the orbital-symmetric three-band Hubbard-Hund model (3HHM) for a lattice filling of 1/3. Hund-correlated metals are characterized by spin-orbital separation (SOS), a Hund's-rule-induced two-stage Kondo-type screening process, in which spin screening occurs at much lower energy scales than orbital screening. By contrast, in Mott-correlated metals, lying close to the phase boundary of a metal-insulator transition, the SOS window becomes negligibly small and the Hubbard bands are well separated. Using dynamical mean-field theory and the numerical renormalization group as real-frequency impurity solver, we identify numerous fingerprints distinguishing Hundness from Mottness in the temperature dependence of various physical quantities. Furthermore, These include ARPES-type spectra, the local self-energy, static local orbital and spin susceptibilities, resistivity, thermopower, and lattice and impurity entropies. Our detailed description of the behavior of these quantities within the context of a simple model Hamiltonian will be helpful for distinguishing Hundness from Mottness in experimental and theoretical studies of real materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗
NASA NTRS1979

Comparison between thermochemical and phase stability data for the quartz-coesite-stishovite transformations

Phase stability and elasticity data have been used to calculate the Gibbs free energy, enthalpy, and entropy changes at 298 K and 1 bar associated with the quartz-coesite and coesite-stishovite transformations in the system SiO2. For the quartz-coesite transformation, these changes disagree by a factor of two or three with those obtained by calorimetric techniques. The phase boundary for this transformation appears to be well determined by experiment; the discrepancy, therefore, suggests that the calorimetric data for coesite are in error. Although the calorimetric and phase stability data for the coesite-stishovite transformation yield the same transition pressure at 298 K, the phase-boundary slopes disagree by a factor of two. At present, it is not possible to determine which of the data are in error. Thus serious inconsistencies exist in the thermodynamic data for the polymorphic transformations of silica.

Weaver, J. S.↗
DOE OSTI2023

On the black hole/string transition

We discuss aspects of the possible transition between small black holes and highly excited fundamental strings. We focus on the connection between black holes and the self gravitating string solution of Horowitz and Polchinski. This solution is interesting because it has non-zero entropy at the classical level and it is natural to suspect that it might be continuously connected to the black hole. Surprisingly, we find a different behavior for heterotic and type II cases. For the type II case we find an obstruction to the idea that the two are connected as classical solutions of string theory, while no such obstruction exists for the heterotic case. We further provide a linear sigma model analysis that suggests a continuous connection for the heterotic case. We also describe a solution generating transformation that produces a charged version of the self gravitating string. This provides a fuzzball-like construction of near extremal configurations carrying fundamental string momentum and winding charges. We provide formulas which are exact in α' relating the thermodynamic properties of the charged and the uncharged solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Molecular Origins of the Barriers to Proton Transport in Acidic Aqueous Solutions

The self-consistent iterative multistate empirical valence bond (SCI-MS-EVB) method is used to analyze the structure, thermodynamics, and dynamics of hydrochloric acid solutions. The reorientation timescales of irreversible proton transport are elucidated by simulating 0.43 M, 0.85 M, 1.68 M, and 3.26 M HCl solutions at 270 K, 285 K, 300 K, 315 K, and 330 K. The results indicate increased counter-ion pairing with increasing concentration, which manifests itself via a reduced hydronium oxygen–chloride (O*-Cl) structuring in the radial distribution functions. Increasing ionic concentration also reduces the diffusion of the hydrated excess protons, principally by reducing the contribution of the Grotthuss proton hopping (shuttling) mechanism to the overall diffusion process. In agreement with prior experimental findings, a decrease in the activation energy of reorientation timescales was also observed, which is explicitly explained using activated rate theory and an energy-entropy decomposition of the state-averaged radial distribution functions. Furthermore, these results provide atomistic verification of suggestions from recent two-dimensional infrared spectroscopy experiments that chloride anions (as opposed to hydrated excess protons) create entropic barriers to proton transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1982

Thermodynamic and transport combustion properties of hydrocarbons with air. Part 1: Properties in SI units

Thermodynamic and transport combustion properties were calculated for a wide range of conditions for the reaction of hydrocarbons with air. Three hydrogen-carbon atom ratios (H/C = 1.7, 2.0, 2.1) were selected to represent the range of aircraft fuels. For each of these H/C ratios, combustion properties were calculated for the following conditions: Equivalence ratio: 0, 0.25, 0.5, 0.75, 1.0, 1.25 Water - dry air mass ratio: 0, 0.03 Pressure, kPa: 1.01325, 10.1325, 101.325, 1013.25, 5066.25 (or in atm: 0.01, 0.1, 1, 10, 50) Temperature, K: every 10 degrees from 200 to 900 K; every 50 degrees from 900 to 3000 K Temperature, R: every 20 degrees from 360 to 1600 R; very 100 degrees from 1600 to 5400 R. The properties presented are composition, density, molecular weight, enthalphy, entropy, specific heat at constant pressure, volume derivatives, isentropic exponent, velocity of sound, viscosity, thermal conductivity, and Prandtl number. Property tables are based on composites that were calculated by assuming both: (1) chemical equilibrium (for both homogeneous and heterogeneous phases) and (2) constant compositions for all temperatures. Properties in SI units are presented in this report for the Kelvin temperature schedules.

Gordon, S.↗
DOE OSTI2026

Nanodomain Formation and Temperature-Dependent Diffusion in Deep Eutectic Solvents Revealed by Single-Molecule Tracking

Deep eutectic solvents (DESs) are typically regarded as homogeneous liquids; however, recent work shows that many exhibit nanoscale structural heterogeneity. Most studies attribute these nanoscale features to short-range chemical interactions. It is still unclear whether a long-range physical mechanism also plays a role. Here, in this study, we examined the nanoscale structure in two hydrophobic DESs, 1:3 tetrabutylammonium bromide: l-menthol (DES-butyl) and 1:3 tetraoctylammonium bromide: l-menthol (DES-octyl). The notation 1:3 represents the molar ratio of the hydrogen bond acceptors to hydrogen bond donors used in the synthesis of the DESs. Single-molecule tracking (SMT) coupled with maximum entropy method (MEM) analysis was used to measure the number of diffusion populations of a dilute concentration of an added fluorescent probe. The presence of more than one population of diffusion coefficients indicates the existence of multiple local environments for the fluorescent probe (i.e., nanoscale structures in the DES). DES-butyl showed a relatively narrow diffusion coefficient distribution centered at 0.55 μm 2 /s, whereas DES-octyl displayed two distinct diffusing populations at 20 °C, with diffusion coefficients of 0.12 μm 2 /s and 0.53 μm 2 /s for the slow and fast populations, respectively. As DES-octyl was heated, the slow-diffusing population steadily diminished and disappeared above ∼30 °C, indicating that the nanodomains present at lower temperatures collapse as the liquid becomes more thermodynamically mixed. This temperature-dependent homogenization is consistent with a physical mechanism of nanostructure formation, for example, liquid–liquid phase separation (LLPS), wherein the structure is not driven solely by specific chemical interactions. The SMT-MEM results suggest that a long-range physical mechanism is the most plausible origin of the measured nanoscale structure in DES-octyl.

Opare-Addo, Jemima [Ames Laboratory (AMES), Ames, ↗
DOE OSTI2021

Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics

Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the “spectral gap optimization of order parameters (SGOOP)” approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. Here, we believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗