Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic phase”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High Capacity Step-Shaped Hydrogen Adsorption in Robust, Pore-Gating Zeolitic Imidazolate Frameworks (Final Technical Report)

The use of porous adsorbents for the densification of H 2 under relatively mild pressures and temperatures has been long pursued through various avenues of synthetic organic frameworks, namely metal–organic frameworks. However, the overwhelming majority of these materials are macroscopically rigid, exhibiting adsorption–desorption profiles that require appreciable energetic input for their full payload capacity to be stored and subsequently delivered. This reduced “deliverable” capacity is a core inefficiency in the use of porous adsorbents for H 2 densification. Our core mission is to determine a material-based strategy that fundamentally alters the adsorption¬–desorption profile in a manner that obviates these deleterious energetic requirements. Towards this, we have focused on metal–organic frameworks that exhibit “cooperative flexibility”, wherein the material undergoes reversible macroscopic phase changes in response to changes to external temperature and/or adsorbate pressure (i.e., concentration). Where these reversible switches occur between states of disparate enough accessible porosity, cooperatively flexible frameworks can exhibit “step-shaped” adsorption–desorption profiles that enable use of the delivery gaseous payload without the aforementioned additional energy requirements. However, before our work, such the observation of this phenomenon with very weakly adsorbing H 2 had not been observed with the necessary profile and P/T conditions desired for H 2 storage and delivery. During this project we have leveraged a deep understanding of the phase change in family of cooperatively flexible metal–organic frameworks to systematically alter their adsorption–desorption profiles to achieve step-shaped adsorption and desorption of H 2 in a precise P/T regime that enables a high deliverable capacity. Specifically, the thermodynamics of the adsorption/desorption induced phase changes in the baseline framework CdIF-13 (sod–Cd(benzimidazolate) 2 ) were tuned using the mixed-linker, or multivariate, approach to framework derivatization. Follow this, we have pursued a similar approach to tuning the adsorption–desorption profile of a second framework family, MIL-53(Al), which contains only Al(III) and theoretically exhibits a much higher deliverable capacity. In total, our work demonstrates that by first understanding the structural consequences of phase changes in these materials, their behavior can be systematically altered through ligand substitutions, tuning the pressure and temperature conditions at which H 2 is adsorbed and desorbed. Thus, enabling high deliverable capacities of H 2 with minimized energetic requirements.

08 HYDROGEN↗
DOE OSTI2020

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

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

Complexions at the iron-magnetite interface

Synthesizing distinct phases and controlling crystalline defects are key concepts in materials design. These approaches are often decoupled, with the former grounded in equilibrium thermodynamics and the latter in nonequilibrium kinetics. By unifying them through defect phase diagrams, we can apply phase equilibrium models to thermodynamically evaluate defects—including dislocations, grain boundaries, and phase boundaries—establishing a theoretical framework linking material imperfections to properties. Using scanning transmission electron microscopy (STEM) with differential phase contrast (DPC) imaging, we achieve the simultaneous imaging of heavy Fe and light O atoms, precisely mapping the atomic structure and chemical composition at the iron-magnetite (Fe/Fe 3 O 4 ) interface. We identify a well-ordered two-layer interface-stabilized phase state (referred to as complexion) at the Fe[001]/Fe 3 O 4 [001] interface. Using density-functional theory (DFT), we explain the observed complexion and map out various interface-stabilized phases as a function of the O chemical potential. The formation of complexions increases interface adhesion by 20% and alters charge transfer between adjacent materials, impacting transport properties. Our findings highlight the potential of tunable defect-stabilized phase states as a degree of freedom in materials design, enabling optimized corrosion protection, catalysis, and redox-driven phase transitions, with applications in materials sustainability, efficient energy conversion, and green steel production.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Characterizing Temporal Heterogeneity by Quantifying Nanoscale Fluctuations in Amorphous Fe‐Ge Magnetic Films

Abstract Equilibrium phase transitions are influenced by fluctuations and often discussed within the framework of the Gibbs free energy, wherein the exchange of energy between system and thermal bath is stationary and all regions of the sample exhibit the same phase. Presence of spatial heterogeneity in the magnetic structures such as pinning centers, domain walls, topological defects, etc. may cause temporal heterogeneity that modifies the nature of the magnetic phase transition. This study reports that interplay of nanoscale thermodynamics with spatio‐temporal heterogeneity gives rise to complex phase transition pathways in amorphous Fe x Ge 1‐x thin films with temperature and Fe‐concentration ( x ). Coherent resonant soft X‐ray scattering experiments that have simultaneous spatial, temporal, and spectral sensitivity show that the origin of helical to paramagnetic phase transition in amorphous Fe‐Ge thin films lies in the appearance of enhanced‐fluctuation spots deep inside the ordered state. The fluctuations are heterogeneous, starting over a small fraction of the domains that increases and becomes isotropic over the entire film as the temperature increases or the Fe‐concentration decreases. The fluctuating‐fraction, when normalized to magnetization for different Fe‐concentrations, follows a single power law behavior, suggesting that the nature of the transition can be described in terms of the underlying spatio‐temporal fluctuations.

Singh, Arnab↗
DOE OSTI2026

Computational Design of Interlayers for Thermally Stable Compositionally Graded Coatings on Nickel Alloys

To extend the service life of Ni-based superalloys, refractory metal coatings are often used. However, direct bonding between metals with dissimilar crystal structure promotes brittle intermetallic phase formation. This work presents a computational thermodynamic framework for high throughput design of functionally graded interlayers to suppress deleterious phases that may form at the interlayer. The Thermo-Calc software package was used to screen candidate metallic interlayer elements based on stability of solid-solution phases. Vanadium was identified as a promising interlayer due to its consistent suppression of intermetallic phases. Temperature-dependent phase diagram mapping between 600 and 1000 °C guided selection of a compositional pathway that significantly reduced intermetallic formation compared to directly joining the Ni-based and Nb refractory alloys. Time–temperature–transformation analysis was performed to assess whether equilibrium-predicted phases are kinetically accessible along regions of the graded path where non-solid-solution phases are not fully suppressed. The methodology was further applied to additional Ni-based alloy and coating systems, illustrating its transferability as an approach for rapid computational design of graded interlayers in dissimilar high-temperature materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Synthesis pathways to thin films of stable layered nitrides

Controlled synthesis of metastable materials away from equilibrium is of interest in materials chemistry. Thin-film deposition methods with rapid condensation of vapour precursors can readily synthesize metastable phases but often struggle to yield the thermodynamic ground state. Growing thermodynamically stable structures using kinetically limited synthesis methods is important for practical applications in electronics and energy conversion. Here, in this study, we reveal a synthesis pathway to thermodynamically stable, ordered layered ternary nitride materials, and discuss why disordered metastable intermediate phases tend to form. We show that starting from elemental vapour precursors leads to a 3D long-range-disordered MgMoN 2 thin-film metastable intermediate structure, with a layered short-range order that has a low-energy transformation barrier to the layered 2D-like stable structure. This synthesis approach is extended to ScTaN 2 , MgWN 2 and MgTa 2 N 3 , and may lead to the synthesis of other layered nitride thin films with unique semiconducting and quantum properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase Stability Through Machine Learning

Understanding the phase stability of a chemical system constitutes the foundation of materials science. Knowledge of the equilibrium state of a system under arbitrary thermodynamic conditions provides valuable information about the types of phases that are likely to be synthesized and how to get there. Accessing the phase diagram in a materials system provides one with the information necessary to design materials and microstructures with optimal properties. While the materials science community has long been focused on exploiting this knowledge to navigate the materials space, recent advances in machine learning (ML) and artificial intelligence (AI) have provided the community with novel ways of interrogating the materials thermodynamics space. Furthermore, this work presents some of the most recent advances in ML/AI applied to phase stability and thermodynamics of materials. Prof. John Morral always had a passion for understanding and teaching the fundamental characteristics of phase diagrams. This review is written to honor his memory.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Exploring the Origin of Anionic Redox Activity in Super Li-Rich Iron Oxide-Based High-Energy-Density Cathode Materials

Super alkali-rich materials (alkali:transition metal ratio of ≥2:1), such as Li 5 FeO 4 , exhibit the potential to realize anionic redox upon deep delithiation. Li 5 FeO 4 was shown to undergo reversible cycling between Li 4 FeO 3.5 and Li 3 FeO 3.5 with combined cation/anion redox, a remarkable capacity of 189 mAh/g, and no O 2 release. However, the impact of phase transformations on the reaction thermodynamics and the correlation between the structural changes and reaction reversibility remain unclear. In this study, we use first-principles calculations to examine the delithiation and (re)lithiation reactions of Li 5 FeO 4 . We show that the experimentally observed charge and discharge processes go through non-equilibrium pathways. Upon delithiation, the compound undergoes a phase transformation from Li 5 FeO 4 , with tetrahedrally coordinated (T d ) Fe ions, to a delithiated disordered rocksalt structure, with octahedral (O h ) Fe ions. Fe-ion migration has an asymmetric kinetic barrier that makes T d → O h migration facile, whereas the reverse process has a much larger barrier, explaining the difficulties in reaction reversibility. We further elucidate the transition metal and O redox sequences during the charge cycle and identify the complex electrochemistry associated with the dual participation of cationic redox (Fe 3+ /Fe 4+ ) and anionic redox (O 2– /O – /O o ). Armed with this knowledge, we conduct high-throughput screening of known alkali-rich transition metal oxides by evaluating their potential to enable reversible anionic redox, with multiple candidates proposed for further experimental trials. Furthermore, our work provides a useful guide for the further development of super alkali-rich anionic-redox-active electrodes for high-energy-density batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

The pervasive impact of critical fluctuations in liquid–liquid extraction organic phases

Liquid-liquid extraction is an essential chemical separation technique where polar solutes are extracted from an aqueous phase into a nonpolar organic solvent by amphiphilic extractant molecules. A fundamental limitation to the efficiency of this important technology is third phase formation, wherein the organic phase splits upon sufficient loading of polar solutes. The nanoscale drivers of phase splitting are challenging to understand in the complex hierarchically structured organic phases. In this study, we demonstrate that the organic phase structure and phase behavior are fundamentally connected in a way than can be understood with critical phenomena theory. For a series of binary mixtures of trialkyl phosphate extractants with linear alkane diluents, we combine small angle x-ray scattering and molecular dynamics simulations to demonstrate how the organic phase mesostructure over a wide range of compositions is dominated by critical concentration fluctuations associated with the critical point of the third phase formation phase transition. These findings reconcile many longstanding inconsistencies in the literature where small angle scattering features, also consistent with such critical fluctuations, were interpreted as reverse micellar-like particles. Altogether, this study shows how the organic phase mesostructure and phase behavior are intrinsically linked, deepening our understanding of both and providing a new framework for using molecular structure and thermodynamic variables to control mesostructure and phase behavior in liquid-liquid extraction.

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

Ab initio thermodynamics reveals the nanocomposite structure of ferrihydrite

Ferrihydrite is a poorly crystalline iron oxyhydroxide nanomineral that serves a critical role as the most bioavailable form of ferric iron for living systems. However, its atomic structure and composition remain unclear due in part to ambiguities in interpretation of X-ray scattering results. Prevailing models so far have not considered the prospect that at the level of individual nanoparticles multiple X-ray indistinguishable phases could coexist. Using ab initio thermodynamics we show that ferrihydrite is likely a nanocomposite of distinct structure types whose distribution depends on particle size, temperature, and hydration. Nanoparticles of two contrasting single-phase ferrihydrite models of Michel and Manceau are here shown to be thermodynamically equivalent across a wide range of temperature and pressure conditions despite differences in their structural water content. Higher temperature and water pressure favor the formation of the former, while lower temperature and water pressure favor the latter. For aqueous suspensions at ambient conditions, their coexistence is maximal for particle sizes up to 12 nm. The predictions inform and help resolve different observations in various experiments.

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

Electrochemical Control for Corrosion in Molten Chlorides During CSP Plant Operation

The Liquid Pathway of the Concentrating Solar Power Generation 3 (CSP Gen3) program proposed low-cost molten chloride salt for energy storage. However, online corrosion control was identified as a remining major risk of the Liquid Pathway approach. This project addressed that risk. Electrochemical solutions for corrosion mitigation during CSP plant operation were investigated and their feasibility and scalability were evaluated. The leading cause of corrosion in molten chloride salt systems was identified as corrosive impurities that form within the salt upon exposure to trace amounts of air and moisture. Leveraging electrochemistry, reduction/oxidation reactions can be employed to remove these corrosive impurities. In Phase 1 of this project, a bench-scale batch electrochemical reactor was designed, fabricated, and used to assess the kinetics and thermodynamics of electrochemical salt purification. In Phase 2, a laboratory-scale flow reactor was designed, fabricated, and used to assess the efficacy of the electrochemical method under flowing conditions. Results show that under proposed operating conditions for the Liquid Pathway Gen3 Pilot Plant, the electrochemical method is significantly more effective at removing impurities than alternative chemical and thermal methods, and that the electrochemical method produces less harmful byproducts. A key advance made in the course of this project was the development of a 2-electrode method for electrochemical purification that is more scalable than previously developed 3 electrode methods. This novel method is based on Magnesium (Mg) electrowinning. A provisional patent based on this invention has been submitted (USPTO Application No. 63/480,355). Additional key advances made during this project include assessment of the effect of dissimilar alloys on corrosion, kinetic and thermodynamic evaluation of thermolysis reactions of impurities within the molten salt, characterization of byproducts of purification reactions, and generation of IP focused on isolating value-added products using molten salt-based electrochemistry that could be deployed to valorize the process (USPTO Application No. 63/478,806). Ultimately, this project represented a step toward feasibility of Liquid Pathway Gen3 CSP. The method developed under this project could significantly reduce capital expenses and operating costs and increase plant profitability by enabling use of less expensive alloys, decreasing maintenance, and increasing plant longevity. Key focus areas for follow-on work have been identified as 1) evaluation of the efficacy of the electrochemical method under turbulent conditions in a larger flow system, such as the FASTR loop, 2) development of methods for removal of purification byproducts, 3) modeling pilot and industrial scale performance of electrochemical salt purification during plant operation and 4) further assessment of the effect of impurities on salt vapor phase.

14 SOLAR ENERGY↗
DOE OSTI2025

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗
DOE OSTI2023

A general non-Fourier Stefan problem formulation that accounts for memory effects

The Stefan problem is the classical model of a melting phase change. In heterogeneous systems, such phase changes can exhibit non-Fourier (anomalous) behaviors, where the advance of the melt interface does not follow the expected time scaling. These situations can be modeled by replacing the derivatives, in the governing partial differential equations, with fractional order derivatives. In particular, replacing the time derivatives leads to non-Fourier models that account for memory effects in the system. In this work, by using appropriate time convolution integrals, a general thermodynamic balance statement for melting phase problems, explicitly accounting for memory effects, is developed. From this balance, a general model formulation applicable to problems involving melting over a temperature range (i.e., a mushy region) is derived. A key component in this model is the representation of memory effects through the use of fractional derivative based constitutive models of the enthalpy and heat flux. Further, on shrinking the mushy region to a single isotherm, a general sharp interface melting model is obtained. Here, in contrast to the classic Stefan problem, the fractional derivatives induce a natural regularization, such that the constitutive models for enthalpy and heat flux are continuous at the melt interface; a result confirmed through numerical simulation. To further support the theoretical findings, a physical example of a non-Fourier Stefan problem is presented. Overall the development and results in this paper underscore the importance of explicitly relating the development of fractional calculus models to the appropriate thermodynamic balance statements.

42 ENGINEERING↗
DOE OSTI2020

Surface science of shape-selective metal nanocrystal synthesis from first-principles: Growth of Cu nanowires and nanocubes

The authors present results from first-principles density functional theory aimed at understanding the aqueous solution-phase growth of fivefold twinned copper nanowires and single-crystal nanocubes capped by hexadecylamine (HDA). The role of solution-phase chloride, present in the Cu salt or as an additive, is emphasized. Using ab initio thermodynamics, the authors delineate the range of solution-phase conditions, characterized by the chemical potentials of chloride and HDA, under which Cu nanowires and nanocubes can be grown. The authors discuss the likelihood of thermodynamic and/or kinetic nanostructures for various solution-phase concentrations. Their results are in good agreement with experiments and indicate that methods and insights developed for surface science in gas-phase or vacuum conditions can yield much insight into liquid-phase systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

New High-Entropy Perovskite Oxides with Increased Reducibility and Stability for Thermochemical Hydrogen Generation

This project aims to design, synthesize, and test a transformative class of High-Entropy Perovskite Oxides (HEPOs) as redox oxides to enable thermochemical hydrogen generation with improved stability, kinetics, and efficiency. These developed HEPOs are expected to demonstrate improved kinetics with oxygen surface exchange coefficient ( k > 7.5×10 -4 cm/s) in Budget Period (BP) 1, retain its structural stability in a broad range of oxygen non-stoichiometry (Δδ > 0.15) at a low operating reduction temperature of T red < 1400°C in BP 2, and deliver a H 2 yield of over 400 µmol per gram of oxide and high stability with less than 20% degradation after at least 50 cycles in BP 3. This project is feasible due to the unique thermodynamic properties (simultaneously increased reducibility and phase stability) and kinetic characters (stability against particle coarsening and potentially enhanced oxygen transport and surface reaction kinetics) of such HEPOs, and it is enabled by a unique active learning computational design approach. Computational studies have been conducted to investigate the oxygen vacancy formation in complex perovskite systems. * Accurate prediction of V O .. concentration with disordered A-site cations in Fe-based high-entropy perovskite oxides * Combined MC/DFT computation elucidates the mechanism of Co preference on the redox due to the strain introduced by local distortion. In this project, we explored a large number (~150) of perovskite compositions, which are listed in Tables 2 – 4). * All perovskite specimens have been synthesized through a high-throughput high-energy ball milling process, followed by the conventional sintering process. * XRD, SEM/EDS and TGA were performed to confirm the crystal structure, phase homogeneity and oxygen non-stoichiometry for compositionally complex perovskite oxides (CCPOs). * 110 compositions show single-phase from XRD * Unusual aliovalent doping effects in medium-entropy perovskite compositions. * V-shape relation between Δδ vs. x (La 1-x Sr x )(Mn 1/3 Fe 1/3 Ti 1/3 )O 3 * The highest reported hydrogen production for the CCPOs made in this project ( T re = 1350 ºC 30 min, T Ox = 1100 ºC 30 min) * B-site mixing (La 0.8 Sr 0.2 )(Mn 0.2 Fe 0.2 Co 0.4 Al 0.2 )O 3 : 89.97 ± 2.73 mmol H2 /mol oxide (395 ± 10 μmol/g oxide ) (i) No phase transformation during reactions when Co molar ratio is less than 61% (ii) Balance between intrinsic kinetics (oxygen surface exchange) and thermodynamics (oxygen non-stoichiometry) (iii) Preference of Co identified by in-situ XPS * A-site mixing (La 1/6 Pr 1/6 Nd 1/6 Gd 1/6 Ba 1/6 Sr1/6)MnO 3 : 98.48 mmol H2 /mol oxid e (~415 μmol/g oxide )

08 HYDROGEN↗
DOE OSTI2023

Thermodynamic modeling of the Nb-Ni system with uncertainty quantification using PyCalphad and ESPEI

Here, the Nb–Ni system is remodeled with uncertainty quantification (UQ) using software tools of PyCalphad and ESPEI (the Extensible, Self-optimizing Phase Equilibria Infrastructure) with the presently implemented capability of modeling site fraction based on Wyckoff positions. The five- and three-sublattice models are used to model the topologically close pack (TCP) μ-Nb 7 Ni 6 and δ-NbNi 3 phases according to their Wyckoff positions. The inputs for CALPHAD-based thermodynamic modeling include the thermochemical data as a function of temperature predicted by first-principles and phonon calculations based on density functional theory (DFT), ab initio molecular dynamics (AIMD) simulations, together with phase equilibrium and site fraction data in the literature. In addition to phase diagram and thermodynamic properties, the CALPHAD-based predictions of site fractions of Nb in μ-Nb 7 Ni 6 agree well with experimental data. Furthermore, the UQ estimation using the Markov Chain Monte Carlo (MCMC) method as implemented in ESPEI is applied to study the uncertainty of site fraction in μ-Nb 7 Ni 6 and enthalpy of mixing (ΔH mix ) in liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗