SEARCH · Engineering Papers

Results for “Theory of nucleation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Modeling Framework to Predict Melting Dynamics at Microstructural Defects in TNT-HMX High Explosive Composites

Many high explosive (HE) formulations are composite materials whose microstructure is understood to impact functional characteristics. Interfaces are known to mediate the formation of hot spots that control their safety and initiation. Here, to study such processes at molecular scales, we developed all-atom force fields (FFs) for Octol, a prototypical HE formulation comprised of TNT (2,4,6-trinitrotoluene) and HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). We extended a FF for TNT and recasted it in a form that can be readily combined with a well-established FF for HMX. The resulting FF was extensively validated against experimental results and density functional theory calculations. We applied the new combined TNT-HMX FF to predict and rank surface and interface energies, which indicate that there is an energetic driver for coarsening of microstructural grains in TNT-HMX composites. Finally, we assess the impact of several microstructural environments on the dynamic melting of TNT crystal under ultrafast thermal loading. We find that both free surfaces and planar material interfaces are effective nucleation points for TNT melting. However, MD simulations show that TNT crystal is prone to superheating by at least 50 K on subnanosecond time scales and that the degree of superheating is inversely correlated with surface and interface energy. The modeling framework presented here will enable future studies on hot spot formation processes in accident scenarios that are governed by strong coupling between microstructural interfaces, material mechanics, momentum and energy transport, phase transitions, and chemistry.

36 MATERIALS SCIENCE↗

DOE OSTI2026

The Fate of Nucleated Black Holes in de Sitter Quantum Gravity

The Euclidean Nariai geometry has long been proposed as the instanton describing the nucleation of maximal-mass black holes in de Sitter space. We place this interpretation on firmer footing by showing that, once an observer is included, the gravitational path integral produces the imaginary phase required for a transition rate. As a warmup, we revisit the Hawking-Moss instanton and, as a byproduct, find that scalar fields can enhance black-hole nucleation, suggesting a quantum-gravity bound on scalar potentials with de Sitter solutions. We then study the subsequent semiclassical evolution of the nucleated black hole. We show that the previously claimed "anti-evaporation" channel is unphysical, arising from a quantum state with singular horizons. In a smooth state, the black hole instead undergoes standard thermal Hawking evaporation. We verify explicit agreement with the no-boundary state and argue that this evaporation is not subject to large quantum-gravity corrections. The nucleated black hole thus evaporates completely back to the maximally-entropic empty de Sitter vacuum, making the full process a Boltzmann fluctuation.

FOS: Physical sciences↗

DOE OSTI2025

Attribution of heterogeneous stress distributions in low-grain polycrystals under conditions leading to damage

In high-purity polycrystalline metallic materials, voids tend to favor grain boundaries as nucleation sites due to the elevated stress states produced by granular interactions and the weakened grain boundary from the relative atomic disorder. To quantify the key factors of this elevated stress state, simple compression of a small multi-grain cylinder of body-centered cubic tantalum was simulated using a single crystal plasticity model that incorporates non-Schmid effects. Four increasingly complex synthetic microstructures were created to tractably incorporate grain boundary interactions, and a statistically significant number of combinations were performed by varying the initial crystallographic orientations of the microstructure. Most of these simulations produce the maximum von Mises stress on a grain boundary and less frequently at the multi-grain junctions. To build a statistical model for the maximum von Mises stress at the grain boundary, physically based features that could contribute to the elevated stress state were selected. Then, a learning algorithm based on information theory was used to identify which of these features contributed the most information to the data set. The identified features include a grain’s propensity to accommodate both elastic and plastic deformations and their directional components. The misalignment of the direction of each grain’s mechanical response was found to be strongly correlated to the magnitude of the stress near the grain boundary. For all of the synthetic microstructures, the statistical models produce a residual distribution that is nearly Gaussian with a variance of, at most, 10% of the prior distribution. The successful performance of the statistical model implies the correct identification of the physical features that cause severe stress localization in polycrystalline materials. The statistical models constructed here can be used to formulate a physically motivated void nucleation model which is sensitive to a microstructure’s propensity to produce elevated stress states. As a result, these statistical models also enable the design of material microstructures, in which the crystallographic orientation is chosen to resist void nucleation.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Exploring poly-crystallization in semiconductors through assumption-less growth simulations: CdTe/CdS case study

Crystal growth is a complex process with far-reaching implications for high-performance materials across various fields. Recent advancements in structural analysis methods such as polyhedral template matching, which allows semiconductor-specific analysis, coupled with simulation technology, have enabled the comprehensive study of crystallization dynamics in semiconductors. However, the exploration of polycrystalline semiconductors created with minimal external intervention of the crystallization processes is relatively uncharted in comparison with metals. In this study, we employ molecular dynamics to simulate the growth of polycrystalline CdTe/CdS with the assumptions of classical mechanics, a Stillinger–Weber potential, an amorphous substrate, and common vapor growth conditions to allow the polycrystalline structures to evolve naturally. Post-simulation, we identify and analyze impactful structures and events, comparing them to theory and experiment to gain insight into various modes of crystallization dynamics. Two research questions guided the study: (1) How realistic are assumption-less simulated polycrystalline semiconductor structures? (2) To what extent can the approach provide insight into crystallization? The simulations, performed with minimal external control, yield polycrystalline structures mirroring experimental findings. The analysis reveals key crystallization insights, such as the role of amorphous atoms in the transition from nucleation to grain growth and the transformative impact of single events, such as dislocations, on crystallization dynamics. The method paves the way for reproducing and analyzing realistic polycrystalline semiconductor structures with minimal simulation assumptions across various growth modes.

Abdullah, Sharmin↗

DOE OSTI2022

Low temperature CVD growth of WSe 2 enabled by moisture-assisted defects in the precursor powder

Two-dimensional transition metal dichalcogenides (TMDs) have been proposed for a wide variety of applications, such as neuromorphic computing, flexible field effect transistors, photonics, and solar cells, among others. However, for most of these applications to be feasible, it is necessary to integrate these materials with the current existing silicon technology. Although chemical vapor deposition is a promising method for the growth of high-quality and large-area TMD crystals, the high temperatures necessary for the growth make this technique incompatible with the processes used in the semiconductor industry. Herein, we demonstrate the possibility of low-temperature growth of TMDs, using tungsten selenide (WSe 2 ) as a model, by simply using moisture-assisted defective tungsten oxide (WO 3 ) precursor powders during the growth of these materials. Density functional theory calculations reveal the mechanism by which moisture promotes the defect formation on the precursor crystal structure and how it dictates the reduction of the temperature of the growth. The results were compared with the standard growth at high temperatures and with a precursor mixture with alkali salts to show the high quality of the WSe 2 grown at temperatures as low as 550 °C. To conclude, the work improves the understanding of nucleation and growth mechanisms of WSe 2 at low temperatures and provides a useful strategy for the growth of TMDs at temperatures required for the back-end-of-line compatibility with current silicon technology.

2D materials↗

NASA NTRS2004

Notes on Interface Growth Kinetics 50 Years after Burton, Cabrera and Frank

This overview is devoted to some unresolved basic problems in crystal growth kinetics. The density wave approach to propagation of a spatially diffuse interface between a growing crystal and its simple (e.g., metallic) melt is discussed is Section 2. This approach allows for the calculation of kinetic coefficients and is an alternative to the localized interface concept in which each atom belongs to either a solid or a liquid. Sections 3 and 4 deals mainly with layer growth from solution. Mutual retardation of the growth steps via their bulk and d a c e diffusion fields is the major subject. The influence of solution flow on step bunching (Section 4) suggests the essential influence of bulk diffusion on the surface morphology. The flow within the solution boundary layer enhances step-step interaction, influences the step bunching process and the resulting step pattern morphology on the growing surface. Recent experiments on the rates at which strongly polygonized steps on protein and small molecule crystals propagate during growth from solution are analyzed in Section 5 . We have shown that the step segments may be "singular" and that "one-dimensional nucleation" may be the rate limiting stage for the segments that are shorter or comparable in length to the thermodynamically equilibrium interlink distance. In this case, the reciprocal dependence of the segment propagation rate on the segment length that follow from the Gibbs-Thompson law, should be replaced by an abrupt switch from zero to a finite constant velocity. Until recently, the Kossel crystal remained the only model used in crystal growth theory. In such Kossel Gibbs-Thomson law, should be replaced by an abrupt switch &om zero to a finite constant velocity. crystals, all kinks at the steps are identical and the kink rate is a linear function of the supersaturation. In the non-Kossel crystals, there may be several kink configurations characterized by different geometries and energies. These configurations should appear in a specific sequence when each new lattice unit cell is filled. As a result of such a cooperative interaction within the unit cell, a non-linear dependence of the kink rate on the vapor pressure or solution concentration in excess over the equilibrium value should be expected.

Chernov, A. A.↗

DOE OSTI2025

Observation of nanoparticle coalescence during core-shell metallic nanowire growth in colloids via nanoscale imaging

The surface morphology and shape of crystalline nanowires significantly influence their functional properties, including phonon transport, electrocatalytic performance, to name but a few. However, the kinetic pathways driving these morphological changes remain underexplored due to challenges in real-space and real-time imaging at single-particle and atomic resolutions. This study investigates the dynamics of shell (Au, Pd, Pt, Fe, Cu, Ni) deposition on AuAg alloy seed nanowires during core-shell formation. By using chiral/non-chiral seed nanowires, advanced imaging techniques, including liquid-phase transmission electron microscopy (LPTEM), cryogenic TEM, and three-dimensional electron tomography, a three-step deposition process is revealed: heterogeneous nucleation, nanoparticle attachment, and coalescence. It is found that colloidal Ostwald ripening, metal reactivity, and deposition amount modulate nanoparticle size and surface roughness, shaping final morphologies. Noble metal nanoparticles (Au, Ag, Pd, Pt) coalesce with seed nanowire along the 〈111〉 direction, distinct from that of other metals. These findings are consistent across different metals, including Ru, Cu, Fe, and Ni, highlighting the hypothesis of these processes in nanowire formation. These findings enhance traditional crystallographic theories and provide a framework for designing nanowire morphology. Additionally, our imaging techniques may be applied to investigate phenomena like electrodeposition, dendrite growth in batteries, and membrane deformation.

Yang, Dahai↗

NASA NTRS1993

Theory of Metastable State Relaxation in a Gravitational Field for Non-Critical Binary Systems with Non-Conserved Order Parameter

A new mathematical ansatz is developed for solution of the time-dependent Ginzburg-Landau nonlinear partial differential equation describing metastable state relaxation in binary (solute+solvent) non-critical solutions with non-conserved scalar order parameter in presence of a gravitational field. It has been demonstrated analytically that in such systems metastability initiates heterogeneous solute redistribution which results in the formation of a non-equilibrium singly-periodic spatial solute structure in the new solute-rich phase. The critical radius of nucleation and the induction time in these systems are gravity-dependent. It has also been proved that metastable state relaxation in vertical columns of supersaturated non-critical binary solutions leads to formation of the solute concentration gradient. Analytical expression for this concentration gradient is found and analysed. It is concluded that gravity can initiate phase separation (nucleation or spinodal decomposition).

Izmailov, Alexander F.↗

DOE OSTI2021

Emergent unitarity in de Sitter from matrix integrals

We study Jackiw-Teitelboim gravity with positive cosmological constant as a model for de Sitter quantum gravity. We focus on the quantum mechanics of the model at past and future infinity. There is a Hilbert space of asymptotic states and an infinite-time evolution operator between the far past and far future. This evolution is not unitary, although we find that it acts unitarily on a subspace up to non-perturbative corrections. These corrections come from processes which involve changes in the spatial topology, including the nucleation of baby universes. There is significant evidence that this 1+1 dimensional model is dual to a 0+0 dimensional matrix integral in the double-scaled limit. So the bulk quantum mechanics, including the Hilbert space and approximately unitary evolution, emerge from a classical integral. We find that this emergence is a robust consequence of the level repulsion of eigenvalues along with the double scaling limit, and so is rather universal in random matrix theory.

2D gravity↗

NASA NTRS2000

Mechanisms for the Crystallization of ZBLAN

A number of research teams have observed that glass forming melts that are solidified in low-g exhibit enhanced glass formation. This project will examine one of these glasses, the heavy metal fluoride glass ZBLAN. A four year ground based research program has been approved to examine the crystallization of ZBLAN glasses with the purpose of testing a theory for the crystallization of ZBLAN glass. The theory could explain the general observations of enhanced glass formation of other glasses melted and solidified in low-g. Fluid flow in 1-g results from buoyancy forces and surface tension driven convection. This fluid flow can introduce shear in undercooled liquids in 1-g. In low-g it is known that fluid flows are greatly reduced so that the shear rate in fluids in low-g are extremely low. It is believed that fluids may have some weak structure in the absence of flow. Even very small shear rates could cause this structure to collapse in response to the shear. A general result would be shear thinning of the fluid. The hypothesis of this research is that: Shear thinning in undercooled liquids increases the rate of nucleation and crystallization of glass forming melts. Shear of the melt can be reduced in low-g enhancing undercooling and glass formation. Samples will be melted and quenched in 1-g under quiescent conditions at a number of controlled cooling rates to determine times and temperatures of crystallization and heated at controlled heating rates to determine kinetic crystallization parameters. Experiments will also be performed on the materials while under controlled vibration conditions and compared with the quiescent experiments in order to evaluate the effect of shear in the liquid on crystallization kinetics. After the experimental parameters are well known, experiments will be repeated under low-g (and 2-g) conditions on the KC-135 aircraft during low-g parabolic maneuvers. The results will determine the effects of shear on crystallization. Our experimental setups will be designed with low-g experiments in mind and will be tested as breadboard low-g experiments. It is very likely that the thermal analysis instrumentation can be adapted to be run in the microgravity glovebox facilities. Critical space experiments may result to test the theory at longer low-g time experiments in space.

Ethridge, Edwin C.↗

DOE OSTI2023

Direct Imaging of Antiferromagnet-Ferromagnet Phase Transition in van der Waals Antiferromagnet CrSBr

The advent of van der Waals (vdW) ferromagnetic (FM) and antiferromagnetic (AFM) materials offers unprecedented opportunities for spintronics and magneto-optic devices. Combining magnetic Kerr microscopy and density functional theory calculations, the AFM-FM transition is investigated and a surprising abnormal magneto-optic anisotropy in vdW CrSBr associated with different magnetic phases (FM, AFM, or paramagnetic state) is discovered. This unique magneto-optic property leads to different anisotropic optical reflectivity from different magnetic states, permitting direct imaging of the AFM Néel vector orientation and the dynamic process of the AFM-FM transition within a magnetic field. Using Kerr microscopy, not only the domain nucleation and propagation process is imaged but also the intermediate spin-flop state in the AFM-FM transition is identified. In conclusion, the unique magneto-optic property and clear identification of the dynamics process of the AFM-FM phase transition in CrSBr demonstrate the promise of vdW magnetic materials for future spintronic technology.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Aggregation of Size-Selected Oxide Clusters Deposited onto Au(111)

Kinetic Monte Carlo (kMC) simulations along with density functional theory (DFT) calculations were used to investigate the aggregation of size-selected Nb 3 O y (y = 5, 6, 7) clusters deposited onto the Au(111) surface. Recent STM experiments showed that the cluster binding sites and sizes of the cluster assemblies on the Nb 3 O y /Au(111) surfaces strongly depend on the stoichiometry of the clusters, i.e., the oxygen-to-niobium ratio. To better understand the origins of these differences, kMC simulations of the nucleation and growth of cluster assemblies were performed using energy barriers for diffusion and intercluster interactions estimated from DFT calculations of cluster binding and dimerization energies, respectively. Comparisons of the kMC simulations with STM images of the as-deposited Nb 3 O y /Au(111) surfaces at RT and after high temperature annealing were used to further optimize the energetics and gauge the importance of nearest neighbor interactions. The kMC simulations demonstrate that the assembly of Nb 3 O y clusters on Au(111) are largely controlled by the magnitude of the barriers for diffusion and interparticle-bond formation, while changes at higher temperatures are sensitive to the binding energies between nearest neighbors. Simulations for the Nb 3 O 5 and Nb 3 O 6 clusters, which exhibit smaller cluster assembly sizes in STM, required larger diffusion barriers as well as different barriers for interparticle binding, which reflected differences in DFT calculated dimerization energies. Furthermore, the results demonstrate the effectiveness of combined DFT and kMC calculations for understanding how the stoichiometry affects the aggregation of small oxide clusters on a metal surface.

36 MATERIALS SCIENCE↗

NASA NTRS1977

The interaction of bubbles with solidification interfaces

The behavior of bubbles at a dendritic solidification interface was studied during the coasting phase of a sounding rocket flight. Sequential photographs of the gradient freeze experiment showed nucleation, growth and coalescence of bubbles at the moving interface during both the low-gravity and one-gravity tests. In the one-gravity test the bubbles were observed to detach from the interface and float to the top of the melt. However, in the low-gravity tests no bubble detachment from the interface or steady state bubble motion occurred and large voids were grown into the crystal. These observations are discussed in terms of the current theory of thermal migration of bubbles and in terms of their implications on the space processing of metals.

Papazian, J. M.↗

DOE OSTI2026

Domain-wall-mediated polarization switching in ferroelectric AlScN: Strain relief and field-dependent dynamics

While scandium-doped aluminum nitride (AlScN) exhibits robust ferroelectricity and excellent thermal stability, its utility is limited by an exceptionally high coercive field (𝐸 𝑐 ) for polarization switching. Unraveling the atomistic switching dynamics is therefore critical for tailoring 𝐸 𝑐 . Here, in this study, we combine density functional theory and machine-learning molecular dynamics to elucidate the polarization switching mechanisms in AlScN over various Sc concentrations and applied electric fields. We find that excessive lattice strain strictly prohibits collective polarization switching, but the preexisting domain walls relieve strain and lead to a distinct switching dynamics—dictating a field-dependent switching mechanism. At low electric fields, switching occurs via gradual domain-wall propagation consistent with the Kolmogorov-Avrami-Ishibashi model. In contrast, high fields stimulate additional nucleation, driving a rapid, homogeneous reversal process described by the simultaneous nonlinear nucleation and growth model. These findings highlight the critical role of domain-wall dynamics and suggest domain engineering as a viable strategy to tailor coercive fields in AlScN and related ferroelectrics.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Closing the loop between in-situ stress complexity and EGS fracture complexity

The goal of the project is to employ a combination of high-fidelity simulations and true-triaxial block fracturing tests at high temperature to explore the intricate relationship between in-situ stress and hydraulic fracture patterns and better characterize the in-situ stress at Utah FORGE. Laboratory experiments are used to validate numerical models and better understand the impact of various parameters (e.g., temperature, well orientation, stress, etc) on the hydraulic fracture nucleation process in EGS. The project will have a significant impact as it will (1) improve the characterization of the in-situ stress field at FORGE; (2) demonstrate how high-fidelity modeling tools can be employed to analyze data for EGS; (3) provide a unique set of high-temperature hydraulic fracturing experimental results for EGS to help identifying the most important components for in-situ stress estimation and to validate the current in-situ stress estimation theory and establish a strong foundation for proposing a new one; (4) provide a validated set of numerical tools in an open source simulation framework (GEOS) that can be used to model EGS.

15 GEOTHERMAL ENERGY↗

NASA NTRS1980

A theory for the origin of a self-replicating chemical system. I - Natural selection of the autogen from short, random oligomers

A general theory is presented for the origin of a self-replicating chemical system, termed an autogen, which is capable of both crude replication and translation (protein synthesis). The theory requires the availability of free energy and monomers to the system, a significant background low-yield synthesis of kinetically stable oligopeptides and oligonucleotides, the localization of the oligomers, crude oligonucleotide selectivity of amino acids during oligopeptide synthesis, crude oligonucleotide replication, and two short peptide families which catalyze replication and translation, to produce a localized group of at least one copy each of two protogenes and two protoenzymes. The model posits a process of random oligomerization, followed by the random nucleation of functional components and the rapid autocatalytic growth of the functioning autogen to macroscopic amounts, to account for the origin of the first self-replicating system. Such a process contains steps of such high probability and short time periods that it is suggested that the emergence of an autogen in a laboratory experiment of reasonable time scale may be possible.

White, D. H.↗

DOE OSTI2023

Effect of the Precursor Chemistry on the Crystallization of Triple Cation Mixed Halide Perovskites

Triple cation, mixed halide perovskite compositions have been reported to be more thermally stable, exhibit fewer phase impurities, and show higher power conversion efficiency and better reproducibility than single cation perovskites. In this work, we explain the formation of Cs0.05FA0.81MA0.14Pb(I0.85Br0.15)3 via a multimodal in situ study combining structural information from synchrotron grazing-incidence wide-angle X-ray scattering (GIWAXS) and optical properties from photoluminescence (PL) spectroscopy with density functional theory calculations (DFT). The focus here is on the effects of the solvent and antisolvent during crystallization. The predominantly used solvents N,N-dimethylformamide (DMF), dimethyl sulfoxide (DMSO), and the antisolvent chlorobenzene (CB) as well as the solvent-antisolvent-precursor interactions are investigated. Given the high elemental complexity and mutual interdependencies between solvent, antisolvent, and perovskite precursors, we found significant differences in the crystallization pathways. DMF-pure precursors show the formation of the DMF-containing intermediate phase and the nucleation of compositionally distinct perovskite phases, while when DMSO is added, only crystalline α- and δ-phases were found. In addition, the presence of DMSO helps the formation of α-perovskite. Coordination energy and bond order (BO) calculations support our experimental findings. Dripping of CB induces nucleation at room temperature, slows the α-phase formation rate, and appears to reduce the nucleation radius. These findings provide novel insights into solvent, antisolvent, and perovskite precursor interactions and their formation pathways. The complexity of interactions between solvents and reagents highlights the importance of understanding these effects to further improve the reproducibility and optimize processing conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Cu–Ni Oxidation Mechanism Unveiled: A Machine Learning-Accelerated First-Principles and in Situ TEM Study

Here, the development of accurate methods for determining how alloy surfaces spontaneously restructure under reactive and corrosive environments is a key, long-standing, grand challenge in materials science. Using machine learning-accelerated density functional theory and rare-event methods, in conjunction with in situ environmental transmission electron microscopy (ETEM), we examine the interplay between surface reconstructions and preferential segregation tendencies of CuNi(100) surfaces under oxidation conditions. Our modeling approach predicts that oxygen-induced Ni segregation in CuNi alloys favors Cu(100)-O c(2 × 2) reconstruction and destabilizes the Cu(100)-O (2√2 × √2)R45° missing row reconstruction (MRR). In situ ETEM experiments validate these predictions and show Ni segregation followed by NiO nucleation and growth in regions without MRR, with secondary nucleation and growth of Cu 2 O in MRR regions. Our approach based on combining disparate computational components and in situ ETEM provides a holistic description of the oxidation mechanism in CuNi, which applies to other alloy systems.

36 MATERIALS SCIENCE↗