Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic systems”

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

Superconductivity in a topological lattice model with strong repulsion

The highly tunable nature of synthetic quantum materials-both in the solid-state and cold atom contexts-invites examining which microscopic ingredients aid in the realization of correlated phases of matter such as superconductors. Recent experimental advances in moiré materials suggest that unifying the features of the Fermi-Hubbard model and quantum Hall systems creates a fertile ground for the emergence of such phases. Here, we introduce the "double Hofstadter"model, a minimal 2D lattice model that incorporates exactly these features: Time-reversal symmetry, band topology, and strong repulsive interactions. By using infinite cylinder density matrix renormalization group methods (cylinder iDMRG), we investigate the ground state phase diagram of this model. We find that it hosts an interaction-induced quantum spin Hall insulator and demonstrate that weakly hole-doping this state gives rise to a superconductor at a finite circumference, with indications that this behavior persists on larger cylinders. Further, at the aforementioned circumference, the superconducting phase is surprisingly robust to perturbations including additional repulsive interactions in the pairing channel. By developing a technique to probe the superconducting gap function in iDMRG, we phenomenologically characterize the superconductor. Namely, we demonstrate that it is formed from the weak pairing of holes atop the quantum spin Hall insulator. Furthermore, we determine the pairing symmetry of the superconductor, finding it to be p-wave-reminiscent of the unconventional superconductivity reported in experiments on twisted bilayer graphene (TBG). Motivated by this, we elucidate structural similarities and differences between our model and those of TBG in its chiral limit. Finally, to provide a more direct experimental realization, we detail an implementation of our Hamiltonian in a system of cold fermionic alkaline-earth atoms in an optical lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Atomic ionization: sd energy imbalance and Perdew–Zunger self-interaction correction energy penalty in 3d atoms

To accurately describe the energetics of transition metal systems, density functional approximations (DFAs) must provide a balanced description of s- and d- electrons. One measure of this is the sd transfer error, which has previously been defined as E ( 3 d n − 1 4 s 1 ) − E ( 3 d n − 2 4 s 2 ) . Theoretical concerns have been raised about this definition due to its evaluation of excited-state energies using ground-state DFAs. A more serious concern appears to be strong correlation in the 4s 2 configuration. Here, we define a ground-state measure of the sd energy imbalance, based on the errors of s- and d-electron second ionization energies of the 3d atoms, that effectively circumvents the aforementioned problems. We find an improved performance as we move from the local spin density approximation (LSDA) to the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) to the regularized and restored Strongly Constrained and Appropriately Normed (r 2 SCAN) meta-GGA for first-row transition metal atoms. However, we find large (∼2 eV) ground-state sd energy imbalances when applying a Perdew–Zunger 1981 self-interaction correction. This is attributed to an “energy penalty” associated with the noded 3d orbitals. A local scaling of the self-interaction correction to LSDA results in a balance of s- and d-errors.

Science & Technology - Other Topics↗
DOE OSTI2019

Interatomic Potential Model Development: Finite-Temperature Dynamics Machine Learning

Developing an accurate interatomic potential model is a prerequisite for achieving reliable results from classical molecular dynamics (CMD) simulations; however, most of the potentials are biased as specific simulation purposes or conditions are considered in the parameterization. Here, for developing an unbiased potential, a finite-temperature dynamics machine learning (FTD-ML) approach is proposed, and its processes and feasibility are demonstrated using the Buckingham potential model and aluminum (Al) as an example. Compared with conventional machine learning approaches, FTD-ML exhibits three distinguished features: 1) FTD-ML intrinsically incorporates more extensive configurational and conditional space for enhancing the transferability of developed potentials; 2) FTD-ML employs various properties calculated directly from CMD, for ML model training and prediction validation against experimental data instead of first-principles data; 3) FTD-ML is much more computationally cost effective than first-principles simulations, especially when the system size increases over 10 3 atoms as employed in this research for ensuring reliable training data. The Al Buckingham potential developed by the FTD-ML approach exhibits good performance for general simulation purposes. Thus, the FTD-ML approach is expected to contribute to a fast development of interatomic potential model suitable for various simulation purposes and conditions, without limitation of model type, while maintaining experimental-level accuracy.

42 ENGINEERING↗
DOE OSTI2022

Simulation of multi-shell fullerenes using Machine-Learning Gaussian Approximation Potential

Multi-shell fullerenes ”buckyonions ” were simulated, starting from initially random configurations, using a density-functional-theory (DFT)-trained machine-learning carbon potential within the Gaussian Approximation Potential (GAP) Framework [Volker L. Deringer and Gábor Csányi, Phys. Rev. B 95, 094203 (2017)]. Fullerenes formed from seven different system sizes, ranging from 60 ~ 3774 atoms, were considered. The buckyonions are formed by clustering and layering starting from the outermost shell and proceeding inward. Inter-shell cohesion is partly due to interaction between delocalized π electrons protruding into the gallery. The energies of the models were validated ex post facto using density functional codes, VASP and SIESTA , revealing an energy difference within the range of 0.02 - 0.08 eV/atom after conjugate gradient energy convergence of the models was achieved with both methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

$\mathrm{Perturbo}$: A software package for ab initio electron–phonon interactions, charge transport and ultrafast dynamics

We report Perturbo is a software package for first-principles calculations of charge transport and ultrafast carrier dynamics in materials. The current version focuses on electron–phonon interactions and can compute phonon-limited transport properties such as the conductivity, carrier mobility and Seebeck coefficient. It can also simulate the ultrafast nonequilibrium electron dynamics in the presence of electron–phonon scattering. Perturbo uses results from density functional theory and density functional perturbation theory calculations as input, and employs Wannier interpolation to reduce the computational cost. It supports norm-conserving and ultrasoft pseudopotentials, spin–orbit coupling, and polar electron–phonon interactions for bulk and 2D materials. Hybrid MPI plus OpenMP parallelization is implemented to enable efficient calculations on large systems (up to at least 50 atoms) using high-performance computing. Taken together, Perturbo provides efficient and broadly applicable ab initio tools to investigate electron–phonon interactions and carrier dynamics quantitatively in metals, semiconductors, insulators, and 2D materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Two-step annealing induced structural rejuvenation: A cause for memory effect in metallic glasses

The structural origin of the memory effect is still elusive in glassy materials. Here, in this letter, the memory effect in a series of metallic glasses (MGs) is observed by two step annealing using differential scanning calorimetry (DSC) and molecular dynamics simulations. It demonstrates that the Kovacs peak can be reflected from not only the system potential energy but also the atomic dynamics, showing the memory effect could originate from the formation of more loose packing configurations around mobile atoms by structural rejuvenation in the second step of high temperature annealing. In contrast, the high-low temperature annealing protocol can only induce the monotonic relaxation. Therefore, our results pave a pathway to better understand the origin of memory effect in MGs and other glassy materials from their dynamics and local atomic structure.

Memory effect↗
DOE OSTI2021

How Accurate Are Approximate Density Functionals for Noncovalent Interaction of Very Large Molecular Systems?

Noncovalent intermolecular interactions are very important in many research areas. Therefore, it is vital to understand the extent to which approximate density functionals give a proper description of noncovalent interactions. Previous research has demonstrated that some approximate density functionals can predict usefully accurate interaction energies for many noncovalent systems; however, most of that work is limited to small and moderate-sized molecules. Very recently though, accurate benchmarks have become available for some very large molecules. Here, the present work applies 21 approximate density functionals to compute the binding energies of seven large molecular systems that have a number of atoms ranging from 200 to 910. The results are judged by comparison to the recently published CIM-DLPNO-CCSD(T) results, which are assumed to provide a reliable benchmark. The five most accurate methods among those tested are found to be PW6B95-D4, PW6B95-D3(BJ), revM11, M06-L, and MN15.

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

Excited State Properties of Point Defects in Semiconductors and Insulators Investigated with Time-Dependent Density Functional Theory

Here, we present a formulation of spin-conserving and spin-flip hybrid time-dependent density functional theory (TDDFT), including the calculation of analytical forces, which allows for efficient calculations of excited state properties of solid-state systems with hundreds to thousands of atoms. We discuss an implementation on both GPU- and CPU-based architectures along with several acceleration techniques. We then apply our formulation to the study of several point defects in semiconductors and insulators, specifically the negatively charged nitrogen-vacancy and neutral silicon-vacancy centers in diamond, the neutral divacancy center in 4H silicon carbide, and the neutral oxygen-vacancy center in magnesium oxide. Our results highlight the importance of taking into account structural relaxations in excited states in order to interpret and predict optical absorption and emission mechanisms in spin defects.

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

Modeling Interactions within and between Peptide Amphiphile Supramolecular Filaments

Many peptides are able to self-assemble into one-dimensional (1D) nanostructures, such as cylindrical fibers or ribbons of variable widths, but the relationship between the morphology of 1D objects and their molecular structure is not well understood. Here, we use coarse-grained molecular dynamics (CG-MD) simulations to study the nanostructures formed by self-assembly of different peptide amphiphiles (PAs). The results show that ribbons are hierarchical superstructures formed by laterally assembled cylindrical fibers. Simulations starting from bilayer structures demonstrate the formation of filaments, whereas other simulations starting from filaments indicate varying degrees of interaction among them depending on chemical structure. These interactions are verified by observations using atomic force microscopy of the various systems. The interfilament interactions are predicted to be strongest in supramolecular assemblies that display hydrophilic groups on their surfaces, while those with hydrophobic ones are predicted to interact more weakly as confirmed by viscosity measurements. The simulations also suggest that peptide amphiphiles with hydrophobic termini bend to reduce their interfacial energy with water, which may explain why these systems do not collapse into superstructures of bundled filaments. The simulations suggest that future experiments will need to address mechanistic questions about the self-assembly of these systems into hierarchical structures, namely, the preformation of interactive filaments vs equilibration of large assemblies into superstructures.

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

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm 2+ , Eu 2+ , Dy 2+ , Yb 2+ , Cf 2+ , and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6(DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics(MD) simulations have been performed on four divalent lanthanide halide salts(Sm 2+ , Eu 2+ , Dy 2+ , and Yb 2+ ) and one divalent actinide halide salt (Cf 2+ ) bound to three diastereoisomers of DCH18C6. Dy 2+ , Yb 2+ , Cf 2+ , DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm 2+ and Eu 2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln 2+ /An 2+ –O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln 2+ /An 2+ -N solvated in acetonitrile for systems containing iodide suggested that theN atom of acetonitrile was competitive with I – for cation interactions. Fluctuation between five and six CNs for Dy 2+ and Yb 2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm 2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Lastly, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

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

Strong Relativistic Effects in Lanthanide-Based Single-Molecule Magnets

Lanthanide-based single-molecule magnets (SMMs) are promising building blocks for quantum memory and spintronic devices. Designing lanthanide-based SMMs with long spin relaxation time requires a detailed understanding of their electronic structure, including the crucial role of the spin-orbit coupling (SOC). While traditional calculations of SOC using the perturbation theory applied to a solution of the non-relativistic Schrödinger equation are valid for light atoms, this approach is questionable for systems containing heavy elements such as lanthanides. We investigate the accuracy of the perturbation estimates of SOC by variationally solving the Dirac equation for the [DyO] + molecule, a prototype of a lanthanide-based SMM. We show that the energy splittings between the M J states involved in spin relaxation depend on the interplay between strong SOC and dynamic electron correlations. Here, we demonstrate that this interplay affects the resonances between the spin and vibrational transitions, and, therefore, the spin relaxation time.

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

Mechanistic Advantages of Organotin Molecular EUV Photoresists

Here, extreme ultraviolet (EUV)-induced radiation exposure chemistry in organotin-oxo systems, represented by the archetypal [(R-Sn) 12 O 14 (OH) 6 ](A) 2 cage, has been investigated with density functional theory. Upholding existing experimental evidence of Sn-C cleavage-dominant chemistry, computations have revealed that either electron attachment or ionization can single-handedly trigger tin-carbon bond cleavage, partially explaining the current EUV sensitivity advantage of metal oxide systems. We have revealed that tin atoms at different parts of the molecule react differently to ionization and electron attachment and have identified such selectivity as a result of local coordination chemistry instead of the macro geometry of the molecule. An ionization-deprotonation pathway has also been identified to explain the observed evolution of an anion conjugate acid upon exposure and anion mass dependence in resist sensitivity.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Origin of giant electric-field-induced strain in faulted alkali niobate films

A large electromechanical response in ferroelectrics is highly desirable for developing high-performance sensors and actuators. Enhanced electromechanical coupling in ferroelectrics is usually obtained at morphotropic phase boundaries requiring stoichiometric control of complex compositions. Recently it was shown that giant piezoelectricity can be obtained in films with nanopillar structures. Here, we elucidate its origin in terms of atomic structure and demonstrate a different system with a greatly enhanced response. This is in non-stoichiometric potassium sodium niobate epitaxial thin films with a high density of self-assembled planar faults. A giant piezoelectric coefficient of ~1900 picometer per volt is demonstrated at 1 kHz, which is almost double the highest ever reported effective piezoelectric response in any existing thin films. The large oxygen octahedral distortions and the coupling between the structural distortion and polarization orientation mediated by charge redistribution at the planar faults enable the giant electric-field-induced strain. Our findings demonstrate an important mechanism for realizing the unprecedentedly giant electromechanical coupling and can be extended to many other material functions by engineering lattice faults in non-stoichiometric compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A graph neural network-state predictive information bottleneck (GNN-SPIB) approach for learning molecular thermodynamics and kinetics

Molecular dynamics simulations offer detailed insights into atomic motions but face timescale limitations. Enhanced sampling methods have addressed these challenges but even with machine learning, they often rely on pre-selected expert-based features. Here, in this work, we present a Graph Neural Network-State Predictive Information Bottleneck (GNN-SPIB) framework, which combines graph neural networks and the state predictive information bottleneck to automatically learn low-dimensional representations directly from atomic coordinates. Tested on three benchmark systems, our approach predicts essential structural, thermodynamic and kinetic information for slow processes, demonstrating robustness across diverse systems. The method shows promise for complex systems, enabling effective enhanced sampling without requiring pre-defined reaction coordinates or input features.

Zou, Ziyue↗
DOE OSTI2020

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

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

Ultrafast state-selective tunneling in two-dimensional semiconductors with a phase- and amplitude-controlled THz-scanning tunneling microscope

THz-pulse driven scanning tunneling microscopy (THz-STM) enables access to the ultrafast quantum dynamics of low-dimensional material systems at simultaneous ultrafast temporal and atomic spatial resolution. State-selective tunneling requires precise amplitude and phase control of the THz pulses combined with quantitative near-field waveform characterization. Here, we employ our state-of-the-art THz-STM with multi-MHz repetition rates, efficient THz generation, and precisely tunable THz waveforms to investigate a single sulfur vacancy in monolayer MoS2. We demonstrate that 2D transition metal dichalcogenides (TMDs) are an ideal platform for near-field waveform sampling by THz cross-correlation. Furthermore, we determine the THz voltage via QEV scans, which measure the THz rectified charge Q as a function of THz field amplitude E and dc bias Vdc. Mapping the complex energy landscape of localized states with a resolution down to 0.01 electrons per pulse facilitates state-selective tunneling to the HOMO and LUMO orbitals of a charged sulfur vacancy.

Materials Science↗
DOE OSTI2025

Keldysh tuning of photoluminescence in a lead halide perovskite crystal

In 1964, Keldysh laid the groundwork for strong-field physics in atomic, molecular, and solid-state systems by delineating a ubiquitous transition from multiphoton absorption to quantum electron tunneling under intense AC driving forces. While both processes in semiconductors can generate carriers and result in photon emission through electron-hole recombination, the low quantum yields in most materials have hindered direct observation of the Keldysh crossover. Leveraging the large quantum yields of photoluminescence in lead halide perovskites, we show that we can not only induce bright light emission from extreme sub-bandgap light excitation but also distinguish between photon-induced and electric-field-induced processes. Our results are rationalized by the Landau-Dykhne formalism, providing insights into the non-equilibrium dynamics of strong-field light-matter interactions. These findings open new avenues for light upconversion and sub-bandgap photon detection, highlighting the potential of lead halide perovskites in advanced optoelectronic applications.

FOS: Physical sciences↗