Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular dynamics (MD)”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Co-scheduling Ensembles of In Situ Workflows

Molecular dynamics (MD) simulations are widely used to study large-scale molecular systems. HPC systems are ideal platforms to run these studies, however, reaching the necessary simulation timescale to detect rare processes is challenging, even with modern supercomputers. To overcome the timescale limitation, the simulation of a long MD trajectory is replaced by multiple short-range simulations that are executed simultaneously in an ensemble of simulations. Analyses are usually co-scheduled with these simulations to efficiently process large volumes of data generated by the simulations at runtime, thanks to in situ techniques. Executing a workflow ensemble of simulations and their in situ analyses requires efficient co- scheduling strategies and sophisticated management of computational resources so that they are not slowing down each other. In this paper, we propose an efficient method to co-schedule simulations and in situ analyses such that the makespan of the workflow ensemble is minimized. We present a novel approach to allocate resources for a workflow ensemble under resource constraints by using a theoretical framework modeling the workflow ensemble’s execution. We evaluate the proposed approach using an accurate simulator based on the WRENCH simulation framework on various workflow ensemble configurations. Results demonstrate the significance of co-scheduling simulations and in situ analyses that couple data together to benefit from data locality, in which inefficient scheduling decisions can lead to slowdown in makespan up to a factor of 30.

Do, Tu Mai Anh↗
DOE OSTI2026

Single-Molecule Electron Transport in Peptoids

Peptoids are structural analogs of peptides in which side chains are appended to the backbone nitrogen rather than the α-carbon. The sequence-defined modularity of peptoids enables precise control over structure−function relationships, enabling applications in energy storage and biomedical materials. Despite recent progress, the role of sequence and conformation on electron transport in peptoid molecules is not fully understood. Here, we synthesize a library of peptoid oligomers and characterize their molecular electronic properties using the scanning tunneling microscope-break junction (STM-BJ) technique. Our results show well-defined electron transport behavior for peptoid sequences containing aromatic side groups lacking hydrogen bonds (H-bonds) and without chemical substitutions at the N−C α position. This behavior fundamentally differs from electron transport in peptides, where H-bond interactions give rise to higher conductance states. All-atom molecular dynamics (MD) simulations are used to understand the conformational heterogeneity of peptoids, and molecular conformations obtained from MD simulations are used in quantum mechanical calculations based on the nonequilibrium Green’s function−density functional theory (NEGF-DFT) formalism. In all cases, computational results are in reasonable qualitative agreement with experiments. Our work demonstrates that the conductance behavior of peptoids depends on monomer identity, including side-chain aromaticity and substitution at the N−C α position. Overall, this work provides new insights into the structure−function relationships governing electron transport in peptoid-based materials and establishes design rules for peptoid-based molecular junctions.

Charge transport↗
DOE OSTI2023

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose–Einstein statistics into an “effective temperature” in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

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

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Development of a deep potential model for F and CF 2 etching of Si and SiO 2

An understanding of plasma-surface interactions at increasingly smaller scales is invaluable for the development of novel technologies and processing techniques. Molecular dynamics (MD) simulations can provide insights into atomic-scale interactions, though they are restricted by the availability of interatomic potentials. Machine learning methods, such as Deep Potential Molecular Dynamics (DeepMD), provide a systematic framework for the development of accurate and flexible ab initio-based models. In this work, we develop DeepMD models for the ion-enhanced etching of Si and SiO 2 by F and CF 2 radicals. We employ an active learning process to expand the data set on which the model is trained and demonstrate its effect on the model accuracy. The DeepMD results are compared to data from classical MD simulations and experiments. Physical sputtering yields of SiO 2 by Ar + ions show good agreement with previous simulation results using conventional interatomic potentials, though the predicted depth profiles are different. Etching yields are calculated as a function of ion energy and neutral to ion flux ratio for the Ar + ion-enhanced etching of SiO 2 and Si by F atoms, as well as for etching of SiO 2 by CF 2 radicals, showing reasonable agreement with experimental data. Finally, an ion-enhanced surface kinetic model is fitted to the DeepMD etch yields, and the fitted parameters are compared to quantities computed directly from DeepMD simulations. This study illustrates how molecular dynamics simulations using machine learning potentials can provide an accurate model of etching processes relevant to device manufacturing.

Kounis-Melas, Andreas [Princeton Univ., NJ (United↗
DOE OSTI2021

Enabling particle applications for exascale computing platforms

The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four “sub-motifs”: short-range particle–particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed particle hydrodynamics (SPH) methods), long-range particle–particle interactions (e.g., electrostatic MD and gravitational N-body), particle-in-cell (PIC) methods, and linear-scaling electronic structure and quantum molecular dynamics (QMD) algorithms. Our crosscutting co-designed technologies fall into two categories: proxy applications (or “apps”) and libraries. Proxy apps are vehicles used to evaluate the viability of incorporating various types of algorithms, data structures, and architecture-specific optimizations and the associated trade-offs; examples include ExaMiniMD, CabanaMD, CabanaPIC, and ExaSP2. Libraries are modular instantiations that multiple applications can utilize or be built upon; CoPA has developed the Cabana particle library, PROGRESS/BML libraries for QMD, and the SWFFT and fftMPI parallel FFT libraries. Success is measured by identifiable “lessons learned” that are translated either directly into parent production application codes or into libraries, with demonstrated performance and/or productivity improvement. The libraries and their use in CoPA’s ECP application partner codes are also addressed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhanced Ion Transport and Molecular Packing Stability in Asymmetric 2D Nanostructured π‐Conjugated Thieno[3,2‐b]Thiophene‐Based Liquid Crystal

Organic semiconductors based on liquid crystal (LC) molecules have attracted increasing interest. Here, in this work, two linear LCs based on 2,5‐bis(thien‐2‐yl)thieno[3,2‐b]thiophene (BTTT) mesogen are designed and synthesized, including BTTT/dEO3 with two symmetrically attached tri(ethylene oxide) groups and BTTT/mEO6 with one asymmetrically attached hexa(ethylene oxide) group. These two molecules have comparable functional‐group compositions but different molecular geometries, leading to their moderately different material performances. Both LCs show smectic mesophases with relatively low transition temperatures as confirmed by differential scanning calorimetry and polarized optical microscopy. A combination of experimental grazing incidence wide‐angle X‐ray scattering and molecular dynamics (MD) simulations reveals a herringbone packing motif of BTTT segments in both LCs while a smaller molecular tilt angle in BTTT/mEO6. Ionic conductivities are measured by doping LCs with different amounts of ionic dopants, lithium bis(trifluoromethanesulfonyl)imide (LiTFSI). BTTT/mEO6 shows better smectic phase stability to higher LiTFSI doping ratios. Both LCs exhibit similar ionic conductivities in the smectic phases, but BTTT/mEO6 outperforms BTTT/dEO3 by a factor of three in the amorphous phase at higher temperatures. MD simulations, performed to examine the ion solvation environment, reveal that BTTT/mEO6 is more efficient in coordinating Li‐ions and screening their interactions with TFSI‐ions which further promote ionic transport.

grazing incidence wide-angle X-ray scattering↗
DOE OSTI2024

Characterize and control the multi-range structure of solution-phase systems with resonant x-ray scattering

This work aims at demonstrating the feasibility and impact of Anomalous X-ray Scattering (AXS) in characterizing the multi-range structures of solution-state systems. We will show preliminary investigation of long-range correlations in concentrated aqueous electrolytes with a combination of AXS and molecular dynamics (MD) simulations. We will also start exploring the capability of AXS to capture intra- and inter-molecular structure of dilute molecular systems, and specifically its sensitivity to the chemical environment surrounding active metal sites. We anticipate that the development of this method will help us to understand ion solvation and transport, which affect the performances of, for instance, electrocatalytic cells, as well as to identify/control the intrinsic structural factors that lead to selectivity, efficiency, and stability control during catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quasi-Classical Trajectory Calculation of Rate Constants Using an Ab Initio Trained Machine Learning Model (aML-MD) with Multifidelity Data

Machine learning (ML) provides a great opportunity for the construction of models with improved accuracy in classical molecular dynamics (MD). However, the accuracy of a ML trained model is limited by the quality and quantity of the training data. Generating large sets of accurate ab initio training data can require significant computational resources. Furthermore, inconsistent or incompatible data with different accuracies obtained using different methods may lead to biased or unreliable ML models that do not accurately represent the underlying physics. Recently, transfer learning showed its potential for avoiding these problems as well as for improving the accuracy, efficiency, and generalization of ML models using multifidelity data. In this work, ab initio trained ML-based MD (aML-MD) models are developed through transfer learning using DFT and multireference data from multiple sources with varying accuracy within the Deep Potential MD framework. Further, the accuracy of the force field is demonstrated by calculating rate constants for the H + HO 2 → H 2 + 3 O 2 reaction using quasi-classical trajectories. We show that the aML-MD model with transfer learning can accurately predict the rate constants while reducing the computational cost by more than five times compared to the use of more expensive quantum chemistry training data sets. Hence, the aML-MD model with transfer learning shows great potential in using multifidelity data to reduce the computational cost involved in generating the training set for these potentials.

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

How Rigid Are Anthranilamide Molecular Electrets?

As important as molecular electrets are for electronic materials and devices, conformational fluctuations strongly impact their macrodipoles and intrinsic properties. Herein, we employ molecular dynamics (MD) simulations with the polarizable charge equilibrium (PQEq) method to investigate the persistence length (L P ) of molecular electrets composed of anthranilamide (Aa) residues. The PQEq-MD dissipates the accepted static notions about Aa macromolecules, and L P represents the shortest Aa rigid segments. The classical model with a single L P value does not describe these oligomers. Introducing multiple L P values for the same macromolecule follows the observed trends and discerns the enhanced rigidity in their middle sections from the reduced stiffness at their terminal regions. Furthermore, L P distinctly depends on solvent polarity. The Aa oligomers maintain extended conformations in nonpolar solvents with L P exceeding 4 nm, while in polar media, increased conformational fluctuations reduce L P to about 2 nm. These characteristics set key guidelines about the utility of Aa conjugates for charge-transfer systems within organic electronics and energy engineering.

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

Solvent-induced membrane stress in biofuel production: molecular insights from small-angle scattering and all-atom molecular dynamics simulations

The disruptive effect of organic solvents on microbial membranes represents a significant challenge to the economical production of green fuels and value-added chemicals from lignocellulosic feedstocks. One route to overcoming this challenge is to engineer microbes with membranes capable of resisting organic solvent stresses. In this regard, it is useful to understand the mechanisms by which organic solvents disrupt typical biomembranes. In this study, molecular dynamics (MD) simulation, complemented by small-angle X-ray and neutron scattering (SANS/SAXS), provide a molecular-scale view of the disruption of a microbial model membrane by 1-butanol and tetrahydrofuran (THF), two common water–organic cosolvent mixtures of importance in biofuel production. Solvent interactions at the interface between the head-group and fatty acid tail regions lead to more dramatic membrane changes than interactions solely at the head-groups or tails. Although both organic solvents are found to partition into the membrane, the depth of solvent penetration into the membrane is quite different. Specifically, 1-butanol localizes near the interface between the lipid heads and tails at low concentrations, but partitions into both the head and tail regions at high concentrations. In contrast, THF, overall, partitions less than 1-butanol and prefers the lipid tail regions. Importantly, the presence of 1-butanol near the head/tail interface introduces drastic membrane changes not seen with THF. The organic solvent interactions with the lipids lead to membrane thinning and fluidization, but more so for 1-butanol than for THF. These results suggest that an aim for the future engineering of robust membranes could be to design lipid head groups that reduce the accumulation of organic solvents at the head–tail interface and that rational designs need also be cognizant of the different solvent-specific mechanisms responsible for membrane disruption.

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

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Diffusion characteristics of water molecules in a lamellar structure formed by triblock copolymers

The distribution and diffusion of water molecules are playing important roles in determining self-assembly and transport properties of polymeric systems. Small-angle neutron scattering (SANS) experiments and molecular dynamics (MD) simulation have been applied to understand the distribution of water molecules and their dynamics in the lamellar membrane formed by Pluronic L62 block copolymers. Penetration of water molecules into the polyethylene oxide (PEO) layers of the membranes has been estimated using scattering length density (SLD) profiles obtained from SANS measurements, which agree well with the molecular distribution observed from MD simulations. The water diffusion coefficient at different regions of the lamellar membrane was further investigated using MD simulation. The diffusion characteristic shows a transition from normal to anomalous diffusion as the position of the water molecule changes from the bulk to PEO and to the polypropylene oxide (PPO) layer. Here we find that water molecules within the PEO or PPO layers follow subdiffusive dynamics, which can be interpreted by the model of fractional Brownian motion.

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

The Single Residue K12 Governs the Exceptional Voltage Sensitivity of Mitochondrial Voltage-Dependent Anion Channel Gating

The voltage-dependent anion channel (VDAC) is a β-barrel channel of the mitochondrial outer membrane (MOM) that passively transports ions, metabolites, polypeptides, and single-stranded DNA. VDAC responds to a transmembrane potential by “gating,” i.e. transitioning to one of a variety of low-conducting states of unknown structure. The gated state results in nearly complete suppression of multivalent mitochondrial metabolite (such as ATP and ADP) transport, while enhancing calcium transport. Voltage gating is a universal property of β-barrel channels, but VDAC gating is anomalously sensitive to transmembrane potential. Here, we show that a single residue in the pore interior, K12, is responsible for most of VDAC’s voltage sensitivity. Using the analysis of over 40 μs of atomistic molecular dynamics (MD) simulations, we explore correlations between motions of charged residues inside the VDAC pore and geometric deformations of the β-barrel. Residue K12 is bistable; its motions between two widely separated positions along the pore axis enhance the fluctuations of the β-barrel and augment the likelihood of gating. Single channel electrophysiology of various K12 mutants reveals a dramatic reduction of the voltage-induced gating transitions. The crystal structure of the K12E mutant at a resolution of 2.6 Å indicates a similar architecture of the K12E mutant to the wild type; however, 60 μs of atomistic MD simulations using the K12E mutant show restricted motion of residue 12, due to enhanced connectivity with neighboring residues, and diminished amplitude of barrel motions. We conclude that β-barrel fluctuations, governed particularly by residue K12, drive VDAC gating transitions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Deep potential molecular dynamics simulations of low-temperature plasma-surface interactions

Machine learning approaches to potential generation for molecular dynamics (MD) simulations of low-temperature plasma-surface interactions could greatly extend the range of chemical systems that can be modeled. Empirical potentials are difficult to generalize to complex combinations of multiple elements with interactions that might include covalent, ionic, and metallic bonds. This work demonstrates that a specific machine learning approach, Deep Potential Molecular Dynamics (DeepMD), can generate potentials that provide a good model of plasma etching in the Si-Cl-Ar system. Comparisons are made between MD results using DeepMD models and empirical potentials, as well as experimental measurements. Pure Si properties predicted by the DeepMD model are in reasonable agreement with experimental results. Simulations of Si bombardment by Ar + ions demonstrate the ability of the DeepMD method to predict sputtering yields as well as the depth of the amorphous-crystalline interface. Etch yields as a function of flux ratio and ion energy for simultaneous Cl 2 and Ar + impacts are in good agreement with previous simulation results and experiment. Predictions of etch yields and etch products during plasma-assisted atomic layer etching of Si-Cl 2 -Ar are shown to be in good agreement with MD predictions using empirical potentials and with experiment. Finally, good agreement was also seen with measurements for the spontaneous etching of Si by Cl atoms at 300 K. Further, the demonstration that DeepMD can reproduce results from MD simulations using empirical potentials is a necessary condition to future efforts to extend the method to a much wider range of systems for which empirical potentials may be difficult or impossible to obtain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗