Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microscopic simulation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantitative Nonlinear Optical Polarimetry with High Spatial Resolution

Nonlinear optical microscopy such as in the optical second-harmonic generation (SHG) modality has become a popular tool today for probing materials in the physical and biological sciences. While imaging and spectroscopy are widely used in the microscopy mode, nonlinear polarimetry, which can shed light on materials’ symmetry and microstructure, is relatively underdeveloped. This is partly because quantitative analytical modeling of the optical SHG response for anisotropic crystals and films largely assumes low-numerical aperture (NA) focusing of light, where the plane-wave approximation is sufficient. Tight focusing provides unique benefits in revealing out-of-plane polarization responses, which cannot be detected by near-plane-wave illumination at normal incidence. Here, we outline a method for quantitatively analyzing SHG polarimetry measurements obtained under high-NA focusing within a microscope geometry. Experiments and simulations of a variety of standard samples, from single crystals to thin films, are in good agreement, including measured and simulated spatial SHG maps of ferroelectric domains. A solution to the inverse problem is demonstrated, where the spatial distribution of an SHG tensor with unknown tensor coefficient magnitudes is determined by experimentally measured polarimetry. The ability to extract the out-of-plane component of the nonlinear polarization in normal incidence is demonstrated, which can be valuable for high-resolution polarimetry of 2D materials, thin films, heterostructures, and uniaxial crystals with a strong out-of-plane response.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deuteration removes quantum dipolar defects from KDP crystals

Abstract Dielectric properties of the hydrogen-bonded ferroelectric crystal KH 2 PO 4 (KDP) differ significantly from those of KD 2 PO 4 (DKDP). It is well established that deuteration affects the interplay of hydrogen-bond switches and heavy ion displacements that underlie the emergence of macroscopic polarization, but a detailed microscopic model is missing. We show that all-atom path integral molecular dynamics simulations can predict the isotope effects, revealing the microscopic mechanism that differentiates KDP and DKDP. Proton tunneling generates phosphate configurations that do not contribute to the polarization. At low temperatures, these quantum dipolar defects are substantial in KDP but negligible in DKDP. These intrinsic defects explain why KDP has lower spontaneous polarization and transition entropy than DKDP. The prominent role of quantum fluctuations in KDP is related to the unusual strength of the hydrogen bonds and should be equally important in other crystals of the KDP family, which exhibit similar isotope effects.

Yang, Bingjia (ORCID:0000000340749553)↗
DOE OSTI2024

Insights into the structure and dynamics of K + ions at the muscovite–water interface from machine learning potential simulations

The surfaces of many minerals are covered by naturally occurring cations that become partially hydrated and can be replaced by hydronium or other cations when the surface is exposed to water or an aqueous solution. These ion exchange processes are relevant to various chemical and transport phenomena, yet elucidating their microscopic details is challenging for both experiments and simulations. Here, in this work, we make a first step in this direction by investigating the behavior of the native K + ions at the interface between neat water and the muscovite mica (001) surface with ab-initio -based machine learning molecular dynamics and enhanced sampling simulations. Our results show that the desorption of the surface K + ions in pure ion-free water has a significant free energy barrier irrespective of their local surface arrangement. In contrast, facile K + diffusion between mica’s ditrigonal cavities characterized by different Al/Si orderings is observed. This behavior suggests that the K + ions may favor a dynamic disordered surface arrangement rather than complete desorption when exposed to deionized water.

Ab-initio molecular dynamics↗
DOE OSTI2026

Controlling Deformation in Al/Ti: How Interface Roughness and Orientation Drive Bimetal Mechanics

The microstructural characteristics and morphology of interfaces in metals can be crucial in governing the initiation of plasticity and early deformation mechanisms under extreme conditions. During high-strain-rate deformation, these interfaces significantly affect dislocation nucleation, twinning, and other mechanisms that directly impact material strength and failure. Despite their importance, a substantial knowledge gap remains between the observed macroscopic material behavior and the underlying role of bimetal interfaces in plasticity initiation. Here, to address this gap, large-scale molecular dynamics simulations are performed on Al/Ti bimetal structures to examine the effect of interface characteristics on the onset of plasticity under uniaxial compression. Specifically, this study investigates how interface roughness (flat vs waveform interfaces) modifies the initiation of plastic events in the microstructure. Atomistic simulations indicate that interface roughness (a microscopic behavior) reduces the stress required for dislocation nucleation, thereby reducing the peak stress relative to a flat interface. For the square interface, the step height strongly influences plasticity initiation by setting the separation of locally flat regions. Varying the interface rotation relative to the loading direction (macroscopic behavior) reveals that the peak stress for both flat and waveform interfaces initially decreases and then increases with rotation, accompanied by a shift in dominant mechanisms—from Al twinning (flat) to interface sliding and Ti-dominated phase transformation and twinning at lower angles. The change in the stress-strain slope during initial compression reflects how rotation alters the resolved shear stress and activates different slip systems. Overall, the study provides valuable insights into the role of bimetal interfaces in controlling plasticity initiation and early deformation pathways for Al/Ti under extreme loading conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Using machine learning to jointly harness the strength of microscopic, fundamental-science driven and macroscopic, application-driven experiments

The PARADIGM project aims at accelerating progress in science by quantitatively answering the following question: What is the optimal combination of fundamental-science and application driven experiments to maximally reduce pertinent data uncertainties? Hence, we are bridging between microscopic experiments and data, and macroscopic simulations and experiments. Answering this question entails solving a high-dimensional and complex optimization problem which we solve with machine learning techniques.

LANSCE↗
DOE OSTI2025

Simulating spin dynamics of supersolid states in a quantum Ising magnet

Motivated by a recent experimental study on the quantum Ising magnet K 2 Co(SeO 3 ) 2 that presented spectroscopic evidence of zero-field supersolidity (Chen et al., arXiv:2402.15869), we simulate the excitation spectrum of the corresponding microscopic XXZ model for the compound, using the recently developed excitation ansatz for infinite projected entangled-pair states. Here, we map out the ground state phase diagram and compute the dynamical spin structure factors across a range of magnetic field strengths, focusing especially on the two supersolid phases found near zero and saturation fields. Our simulated excitation spectra for the zero-field supersolid “Y” phase are in excellent agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Electrolyte-Dependent, “Microscopically Irreversible” H-Atom Transfer Kinetics of Ce-Based Metal–Organic Framework, Ce-MOF-808

Redox reactions at the interface of metal oxides and protic electrolytes almost always involve protons and electrons in equal amounts. Given the stoichiometry, these proton-coupled electron transfer (PCET) reactions are thermochemically equivalent to net H-atom transfer (HAT) reactions. The correlation between the chemical nature of solid catalysts and HAT kinetics has been employed for decades as the design principle for energy-relevant reactions (e.g., reactions of 2H + /H 2 ). More recently, chemists have experimentally determined that a change in liquid electrolytes that alters the microenvironment at the redox-active sites has an equally profound impact on electrocatalysis involving PCET/HAT. Yet, precise correlations between the chemical nature of electrolytes and the PCET kinetics are, to date, rare in the literature. Herein, we report our findings using the Ce-based metal−organic framework, Ce-MOF-808, as a model system. Each Ce 6 (μ 3 −O) 4 (μ 3 − OH) 4 (OH) 6 (H 2 O) 6 node of this MOF undergoes a 1H + /1e − redox reaction. Using chronoamperometry and the Cottrell analysis, we have determined that the PCET hopping kinetics within the pores of Ce-MOF-808 can change by orders of magnitude by altering the buffer species and the proton activity of the electrolyte. Furthermore, in all buffers, reductive reactions were ∼3−10 times faster in kinetics than the reverse oxidative reaction with the same electrochemical driving force, suggesting that the system, at first glance, violates the principle of microscopic reversibility. Isothermal titration calorimetry (ITC) and computational simulations corroborated that the buffer-node binding thermodynamics are quite distinct, depending on the chemical nature of the buffer and the oxidation state of the node. Together, these results suggest that the substrate and the product during the oxidative vs reductive reaction of Ce-MOF-808 are chemically different species, which explains the apparent ‘microscopic irreversibility.’ Thus, the rational modulation of electrolytes can dramatically enhance PCET kinetics, even though the solid electrodes remain identical. Implications of these findings are contrasted with the electrochemical/electrocatalytic behavior of other redox-active MOFs, heterogeneous catalysts, and enzymatic systems at the solid−liquid interface.

Ce-based MOF↗
DOE OSTI2025

Semiconductor‐Electrocatalyst Contacts: Theory, Experiment, and Applications to Solar Water Photoelectrolysis

The goal of the work supported by this program was to use new types of photoelectrochemical experiments, supported by simulation and theory, to develop an improved microscopic picture of electron transfer in combined semiconductor|catalyst|solution (sem|cat|sol) systems. An emphasis was placed on developing techniques which can directly characterize electronic interfacial phenomena without first assuming a junction model. We then aimed to illustrate the design principles that enhance interfacial processes and enable high-efficiency devices. Our work thus provides deep insight into the processes that, in general, use light to drive chemical reactions.

14 SOLAR ENERGY↗
DOE OSTI2025

Machine-learning-informed scattering correlation analysis of sheared colloids

We have carried out theoretical analysis, Monte Carlo simulations and machine-learning analysis to quantify microscopic rearrangements of dilute dispersions of spherical colloidal particles from coherent scattering intensity. Both monodisperse and polydisperse dispersions of colloids were created and underwent a rearrangement consisting of an affine simple shear and non-affine rearrangement using the Monte Carlo method. We calculated the coherent scattering intensity of the dispersions and the correlation function of intensity before and after the rearrangement and generated a large data set of angular correlation functions for varying system parameters, including number density, polydispersity, shear strain and non-affine rearrangement. Singular value decomposition of the data set shows the feasibility of machine-learning inversion from the correlation function for the polydispersity, shear strain and non-affine rearrangement using only three parameters. A Gaussian process regressor is then trained on the data set and can retrieve the affine shear strain, non-affine rearrangement and polydispersity with relative errors of 3%, 1% and 6%, respectively. Altogether, our model provides a framework for quantitative studies of both steady and non-steady microscopic dynamics of colloidal dispersions using coherent scattering methods.

Gaussian process regression↗
DOE OSTI2026

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

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

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2025

Terahertz-field activation of polar skyrons

Unraveling collective modes arising from coupled degrees of freedom is crucial for understanding complex interactions in solids and developing new functionalities. Unique collective behaviors emerge when two degrees of freedom, ordered on distinct length scales, interact. Polar skyrmions, three-dimensional electric polarization textures in ferroelectric superlattices, disrupt the lattice continuity at the nanometer scale with nontrivial topology, leading to previously unexplored collective modes. Here, using terahertz-field excitation and femtosecond x-ray diffraction, we discover subterahertz collective modes, dubbed “skyrons”, which appear as swirling patterns of atomic displacements functioning as atomic-scale gearsets. The key to activating skyrons is the use of the THz field that couples primarily to skyrmion domain walls. Momentum-resolved time-domain measurements of diffuse scattering reveal an avoided crossing in the dispersion relation of skyrons. Atomistic simulations and dynamical phase-field modeling provide microscopic insights into the three-dimensional crystallographic and polarization dynamics. The amplitude and dispersion of skyrons are demonstrated to be controlled by sample temperature and electric-field bias. The discovery of skyrons and their coupling with terahertz fields opens avenues for ultrafast control of topological polar structures.

ferroelectrics↗
DOE OSTI2026

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗
DOE OSTI2024

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

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

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

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

Efficiently Measuring 𝑑-Wave Pairing and Beyond in Quantum Gas Microscopes

Understanding the mechanism of high-temperature superconductivity is among the most important problems in physics, one for which quantum simulation can provide new insights. However, it remains challenging to characterize superconductivity in existing cold-atom quantum simulation platforms. Here, in this work, we introduce a protocol for measuring a broad class of observables in fermionic quantum gas microscopes, including long-range superconducting pairing correlations (after a repulsive-to-attractive mapping). The protocol only requires global controls followed by site-resolved particle number measurements—capabilities that have been already demonstrated in multiple experiments—and is designed by analyzing the Hilbert space of fermions on two sites. The protocol is sample-efficient and we further optimize our pulses for robustness to experimental imperfections such as lattice inhomogeneity. Our Letter introduces a general tool for manipulating quantum states on optical lattices, enhancing their ability to tackle problems such as high-temperature superconductivity.

Hubbard model↗
DOE OSTI2026

Microwave Imaging of Edge Conductivity in Graphene at Charge Neutrality and Quantum Hall States

Here, we report local conductivity imaging of edge states in monolayer graphene by millikelvin microwave impedance microscopy. At the charge-neutrality point, as the magnetic field increases, the local conductivity at the edge drops to zero more slowly than in the bulk. This behavior is consistent with the calculated spatial profile of the charge gap in the canted antiferromagnetic phase. For comparison, we also perform microwave imaging of integer quantum Hall states away from neutrality, which host dissipationless chiral edge channels. The evolution of the edge signal as a function of the bulk gap is fundamentally different between the Landau level filling factor 𝜈 = 0 and |𝜈| ≥ 1 integer quantum Hall states, which can be qualitatively explained by numerical simulations and theoretical analysis. Our results provide a comprehensive microscopic picture of the edge and bulk states as the Fermi level moves across the unique Landau-level spectrum of graphene.

Hall effect↗
DOE OSTI2025

Unraveling Exciton Trap Dynamics and Nonradiative Loss Pathways in Quantum Dots via Atomistic Simulations

Surface defects in colloidal quantum dots are a major source of nonradiative losses, yet the microscopic mechanisms underlying exciton trapping and recombination remain elusive. Here, we develop a model Hamiltonian based on atomistic electronic calculations to investigate exciton dynamics in CdSe/CdS core/shell QDs containing a single hole trap introduced by an unpassivated sulfur atom. By systematically varying the defect depth and reorganization energy, we uncover how defect-induced excitonic states mediate energy relaxation pathways. Our simulations reveal that a single localized defect can induce a rich spectrum of excitonic states, leading to multiple dynamical regimes, from slow, energetically off-resonant trapping to fast, cascaded relaxation through in-gap defect states. Crucially, we quantify how defect-induced polaron shifts and exciton-phonon couplings govern the balance between efficient radiative emission and rapid nonradiative decay. These insights clarify the microscopic origin of defect-assisted loss channels and suggest pathways for tailoring QD optoelectronic properties via surface and defect engineering.

Defects↗