Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic simulations”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Study of vacuum ultraviolet emission in helium and helium/nitrogen mixtures

Helium is frequently used as a working medium for the generation of plasmas and is capable of energetic photon emissions. These energetic photon emissions are often attributed to the formation of helium excimer and subsequent photon emission. When the plasma device is exposed to another gas, such as nitrogen, this energetic photon emission can cause photoionization and further ionization wave penetration into the additional gas. Often ignored are the helium resonance emissions that are assumed to be radiation trapped and therefore not pertinent to photoionization. Here, experimental evidence for the presence of helium atomic emission in a pulsed discharge at ten’s of Torr is shown. Simulations of a discharge in similar conditions agree with the experimental measurements. In this context, the role of atomic and molecular helium light emission on photoionization of molecular nitrogen in an ionization wave is studied using a kinetic modeling approach that accounts for radiation dynamics in a developing low-temperature plasma. Three different mixtures of helium at a total pressure of 250 Torr are studied in simulation. Photoionization of the nitrogen molecule by vacuum ultraviolet helium emission is used as the only seed source ahead of the ionization front. It is found that even though radiation trapped, the atomic helium emission lines are the significant source of photoionization of nitrogen. The significant effect of radiation trapped photon emission on ionization wave dynamics demonstrates the need to consider these radiation dynamics in plasma reactors where self-absorbed radiation is ignored.

42 ENGINEERING↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

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

Hybrid Basis and Multi-Center Grid Method for Strong-Field Processes

We present a time-dependent framework that combines a hybrid basis, consisting of Gaussian-type orbitals (GTOs) and finite-element discrete-variable representation (FEDVR) functions, with a multicenter grid to simulate strong-field and attosecond dynamics in atoms and molecules. The method incorporates the construction of the orthonormal hybrid basis, the evaluation of electronic integrals, a unitary time-propagation scheme, and the extraction of optical and photoelectron observables. Its accuracy and robustness are benchmarked on one-electron systems such as atomic hydrogen and the dihydrogen cation (H$^+_2$) through comparisons with essentially-exact reference results for bound-state energies, high-harmonic generation spectra, photoionization cross sections, and photoelectron momentum distributions. This work establishes the groundwork for its integration with quantum-chemistry methods, which is already operational but will be detailed in future work, thereby enabling ab initio simulations of correlated polyatomic systems in intense ultrafast laser fields.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Outer midplane neutral density measurements and H-mode fueling studies in NSTX-U

Neutral density profiles are measured on the outboard midplane of the National Spherical Torus Experiment Upgrade (NSTX-U) using a two dimensional camera (ENDD, edge neutral density diagnostic) filtered for deuterium Balmer α (D α ) emission interpreted via simulations using the Monte Carlo neutral transport code DEGAS 2. Deuterium atomic densities n D are calculated by inverting the line-integrated D α brightness and using local measurements of electron density ne and temperature T e to determine atomic rate coefficients. The assumptions used in the derivation of n D from D α emissivity are validated using DEGAS 2 to estimate contributions to emissivity due to electron impact excitation and molecular processes. In this work, experimental measurements and DEGAS 2 simulations are compared over a database of L- and H-mode discharges, showing good agreement in D α emissivity profiles. Residual disagreement between experiment and simulations is suggestive of the possible role of intermittent transport and uncertainties in the molecular data. DEGAS 2 simulations are further used to complement the ENDD diagnostic, extracting neutral (atomic and molecular) densities at locations where the ENDD measurement are not made or where the assumptions used in the ENDD analysis are not valid. Deuterium atomic densities and ionization profiles on the outer midplane are compared for L-mode and H-mode discharges in NSTX-U. One-way coupling of DEGAS 2 to UEDGE multi-fluid simulations constrained by experimental data is used to study edge fueling and neutral penetration in NSTX-U discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Plasma-assisted atomic layer etching of single-crystal diamond

Applications of near-surface nitrogen-vacancy (NV) centers in diamond are often limited by surface defects created during processing. Understanding and controlling plasma-induced surface damage is important for preserving the optical and spin properties of diamond NV centers. We report molecular dynamics simulations of a novel form of plasma-aided atomic layer etching of diamond. In this proposed scheme, the initial surface modification step consists of Ar + ion bombardment. This creates an amorphous layer at the surface, the thickness of which is controlled by the ion energy. Amorphization is known to help smooth the surface, at least locally, and would also serve to sputter clean an initially contaminated surface. A second step impacts the amorphous carbon layer with O + between about 1 and 5 eV. Simulations show that this energy range will remove the amorphous carbon but will not etch the underlying diamond. Though this energy range is not trivial to achieve in a conventional low-temperature plasma, potential methods for creating these low-energy O + ions are discussed. In addition to a-C etching, the O + impacts are predicted to remove any isolated (100) diamond terraces by selectively attacking the edges of the terraces. The proposed atomic layer etching (ALE) approach reverses the conventional ALE sequence in which the surface modification step is usually chemical modification (oxidation in this case), followed by a removal step using Ar + impacts. This plasma ALE procedure is predicted to create a diamond surface that is atomically flat and defect free.

Draney, J. S. [Princeton University, NJ (United St↗
DOE OSTI2020

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2020

Identifying Dense NiSe 2 /CoSe 2 Heterointerfaces Coupled with Surface High-Valence Bimetallic Sites for Synergistically Enhanced Oxygen Electrocatalysis

Constructing heterostructures with abundant interfaces is essential for integrating the multiple functionalities in single entities. Herein, the synthesis of NiSe 2 /CoSe 2 heterostructures with different interfacial densities via an innovative strategy of successive ion injection is reported. The resulting hybrid electrocatalyst with dense heterointerfaces exhibits superior electrocatalytic properties in an alkaline electrolyte, superior to other benchmarks and precious metal catalysts. Advanced synchrotron techniques, post structural characterizations, and density functional theory (DFT) simulations reveal that the introduction of atomic-level interfaces can lower the oxidation overpotential of bimetallic Ni and Co active sites (whereas Ni$^{2+}$ can be more easily activated than Co$^{2+}$) and induce the electronic interaction between the core selenides and surface in situ generated oxides/hydroxides, which play a critical role in synergistically reducing energetic barriers and accelerating reaction kinetics for catalyzing the oxygen evolution. Hence, the heterointerface structure facilitates the catalytic performance enhancement via increasing the intrinsic reactivity of metallic atoms and enhancing the synergistic effect between the inner selenides and surface oxidation species. This work not only complements the understanding on the origins of the activity of electrocatalysts based on metal selenides, but also sheds light on further surface and interfacial engineering of advanced hybrid materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular level insight into non-bilayer structure formation in thylakoid membranes: a molecular dynamics study

In oxygenic photosynthetic organisms, the light reactions are performed by protein complexes embedded in the lipid bilayer of thylakoid membranes (TMs). The organization of the bulk lipid molecules into bilayer structures provide optimal conditions for the build-up of the proton motive force (pmf) and its utilization for ATP synthesis. However, the lipid composition of TMs is dominated by the non-bilayer lipid species monogalactosyl diacylglycerol (MGDG), and functional plant TMs, besides the bilayer, contain large amounts of non-bilayer lipid phases. Bulk lipids have been shown to be associated with lumenal, stromal-side and marginal-region proteins and proposed to play roles in the self-assembly and photoprotection of the photosynthetic machinery. Furthermore, it has recently been pointed out that the generation and utilization of pmf for ATP synthesis according to the ‘protet’ or protonic charge transfer model Kell (Biochim Biophys Acta Bioenerg 1865(4):149504, 2024), requires high MGDG content Garab (Physiol Plant 177(2):e70230, 2025). In this study, to gain better insight into the structural and functional roles of MGDG, we employed all atom and coarse-grained molecular dynamics simulations to explore how temperature, hydration levels and varying MGDG concentrations affect the structural and dynamic properties of bilayer membranes constituted of plant thylakoid lipids. Our findings reveal that MGDG promotes increased membrane fluidity and dynamic fluctuations in membrane thickness. MGDG-rich stacked bilayers spontaneously formed inverted hexagonal phases; these transitions were enhanced at low hydration levels and at elevated but physiologically relevant temperatures. It can thus be inferred that MGDG plays important roles in heat and drought stress mechanisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Helium causing disappearance of a/2<111> dislocation loops in binary Fe-Cr ferritic alloys

In this work, single and dual-beam self-ion irradiations were performed at 500°C on ultra-high purity Fe14%Cr alloy to ~0.33 displacements-per-atom (dpa), and 0 or 3030 atomic-parts-per-million (appm) helium/dpa, respectively. Using transmission electron microscopy, we reveal that helium can drastically modify the dislocation loop Burgers vector in Fe-Cr alloys. Helium co-implantation caused complete disappearance of a/2<111> type dislocation loops, and the microstructure consisted of only a<100> loops. Conversely, a/2<111> type loops were predominant without He co-implantation. The total loop density remained largely unaffected. The results strikingly contrast literature asserting that helium stabilizes a/2<111> type loops in bcc Fe alloys, based on low temperature irradiations. Collectively analyzing the results with literature suggest that the small positive interaction between helium and self-interstitial atoms (SIA) in Fe predicted by atomistic simulations maybe insufficient to holistically explain the dislocation loop microstructure development in presence of helium. Helium-SIA positive binding inadvertently implies elevated a/2<111> loop fraction and higher loop densities that the present results contradict. Helium induced high cavity density causing a preferential loss of highly glissile <111> clusters, leaving the matrix saturated with <100> type clusters is proposed as a potential mechanism. Further, the in-situ irradiations combined with Burgers vector analysis strengthened the evidence of Cr-induced dislocation loop mobility reduction that appears to stabilize the a/2<111> type loops and causes higher loop densities in Fe-Cr alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sputter deposition of high electrical resistivity Au-Ta alloy coatings on rotating substrates

Sputter deposition of gold-tantalum alloy coatings is a key process for manufacturing hohlraums for magnetically-assisted inertial confinement fusion implosions. In this report we describe direct current magnetron sputter deposition of ~ 10 - μ~ 10 - μm-thick films of Au-80 at.% Ta onto rotating sphero-cylindrical hohlraum and planar Si witness substrates. Emphasis is given to how film microstructure and properties are affected by main deposition parameters, including argon working gas pressure, substrate bias, and the source composition (a single alloyed target compared to co-sputtering from two elemental targets). Experimental findings are correlated with distributions of landing energies and incident angles of depositing species calculated by Monte Carlo simulations of ballistic collisions and gas phase atomic transport. Deposition conditions characterized by low energetics of depositing species favor the formation of a β-Ta-like phase. Implications of these results to hohlraum fabrication are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Programming Amphiphilic Peptoid Oligomers for Hierarchical Assembly and Inorganic Crystallization

Natural organisms make a wide variety of exquisitely complex, nano-, micro-, and macroscale structured materials in an energy-efficient and highly reproducible manner. During these processes, the information-carrying biomolecules (e.g., proteins, peptides, and carbohydrates) enable (1) hierarchical organization to assemble scaffold materials and execute high-level functions and (2) exquisite control over inorganic materials synthesis, generating biominerals whose properties are optimized for their functions. Inspired by nature, significant efforts have been devoted to developing functional materials that can rival those natural molecules by mimicking in vivo functions using engineered proteins, peptides, DNAs, sequence-defined synthetic molecules (e.g., peptoids), and other biomimetic polymers. Among them, peptoids, a new type of synthetic mimetics of peptides and proteins, have received particular attention because they combine the merits of both synthetic polymers (e.g., high chemical stability and efficient synthesis) and biomolecules (e.g., sequence programmability and biocompatibility). The lack of both chirality and hydrogen bonds in their backbone results in a highly designable peptoid-based system with reduced structural complexity and side chain-chemistry-dominated properties. Here in this Account, we present our recent efforts in this field by programming amphiphilic peptoid sequences for (1) the controlled self-assembly into different hierarchically structured nanomaterials with favorable properties and (2) manipulating inorganic (nano)crystal nucleation, growth, and assembly into superstructures. First, we designed a series of amphiphilic peptoids with controlled side chain chemistries that self-assembled into 1D highly stiff and dynamic nanotubes, 2D membrane-mimetic nanosheets, hexagonally patterned nanoribbons, and 3D nanoflowers. These crystalline nanostructures exhibited sequence-dependent properties and showed promise for different applications. The corresponding peptoid self-assembly pathways and mechanisms were also investigated by leveraging in situ atomic force microscopy studies and molecular dynamics simulations, which showed precise sequence dependency. Second, inspired by peptide- and protein-controlled formation of hierarchical inorganic nanostructures in nature, we developed peptoid-based biomimetic approaches for controlled synthesis of inorganic materials (e.g., noble metals and calcite), in which we took advantage of the substantial side chain chemistry of peptoids and investigated the relationship between the peptoid sequences and the morphology and growth kinetics of inorganic materials. For example, to overcome the challenges (e.g., complexity of protein- and peptide-folding, poor thermal and chemical stabilities) facing the area of protein- and peptide-controlled synthesis of inorganic materials, we recently reported the design of sequence-defined peptoids for controlled synthesis of highly branched plasmonic gold particles. Moreover, we developed a rule of thumb for designing peptoids that predictively enabled the morphological evolution from spherical to coral-shaped gold nanoparticles (NPs). With this Account, we hope to stimulate the research interest of chemists and materials scientists and promote the predictive synthesis of functional and robust materials through the design of sequence-defined synthetic molecules.

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

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

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

Modular Approach for the Synthesis of Bottlebrush Diblock Copolymers from Poly(Glycidyl Methacrylate)-block-Poly(Vinyldimethylazlactone) Backbones

A strategy that enables the facile synthesis of bottlebrush block copolymers with flexible backbones was developed. A demonstration of the strategy’s utility was carried out by grafting end-functionalized polymethylmethacrylate (PMMA) and polystyrene (PS) to the dually reactive block copolymer, poly(glycidyl methacrylate)-block-poly(vinyldimethylazlactone) (PGMA-b-PVDMA). Five different bottlebrush diblock copolymers were investigated by size-exclusion chromatography (SEC), 1H NMR, Fourier transform infrared (FT-IR), differential scanning calorimetry (DSC), X-ray scattering methods, atomic force microscopy (AFM), rheology and computational simulations using molecular dynamics (MD), and self-consistent field theory (SCFT). A relationship between the glass transition temperature and the fraction of chain ends was demonstrated by DSC and highlights the potential of this synthetic method to tailor polymer properties. All five samples were found to be in a disordered phase exhibiting multiscale structures revealed by two broad peaks in small-angle X-ray scattering (SAXS) that can be attributed to graft-to-graft and backbone-to-backbone density correlations using MD simulations. The SCFT-based simulations justify the observation of a disordered phase due to its stabilization by the grafts. Additionally, this modular approach can be easily extended to other grafts, including responsive, conducting, and charged polymers with the prerequisite end groups. The versatility and ease of assembling these functional bottlebrushes constitute a powerful “toolbox” method for the rapid and scalable synthesis of novel bottlebrush block copolymers with desired properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hydrogen bonded trimesic acid networks on Cu(111) reveal how basic chemical properties are imprinted in HR-AFM images

High resolution non-contact atomic force microscopy measurements characterize assemblies of trimesic acid molecules on Cu(111) and the link group interactions, providing the first fingerprints utilizing CO-based probes for this widely studied paradigm for hydrogen bond driven molecular self assembly. The enhanced submolecular resolution offered by this technique uniquely reveals key aspects of the competing interactions. Accurate comparison between full-density-based modeled images and experiment allows to identify key structural elements in the assembly in terms of the electron-withdrawing character of the carboxylic groups, interactions of those groups with Cu atoms in the surface, and the valence electron density in the intermolecular region of the hydrogen bonds. This study of trimesic acid assemblies on Cu(111) combining high resolution atomic force microscopy measurements with theory and simulation forges clear connections between fundamental chemical properties of molecules and key features imprinted in force images with submolecular resolution.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Structure, bonding and electronic characteristics of amorphous Se

Ovonic threshold switching (OTS) selector can effectively improve the storage density and suppress the leakage current of advanced phase-change memory. As a prototypical OTS material, amorphous GeSe is widely investigated. But the attentions paid to amorphous Se (i.e., the functional constituent in amorphous GeSe) is very limited up to now. Here, in this work, we have explored the structure, bonding and electronic characteristics of amorphous Se using ab initio molecular dynamics simulations. The results reveal that the Se atoms in amorphous Se tend to form the 2-coordinated configurations, and they connect with each other to form the long chains. The fraction of vibrational density of state located in high frequency range is relatively large, and the formation energy of Se-Se bond is as large as 4.44 eV, hinting that the Se-Se bonds in chains possess a high stability. In addition, the mid-gap state related to the OTS behavior is also found in the amorphous Se despite of the small proportion. Our findings enrich the knowledge of the amorphous Se, which contributes to the applications of Se-based OTS selectors.

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

Assessing the roles of Cu- and Ag-deficient layers in chalcopyrite-based solar cells through first principles calculations

Chalcopyrites are a demonstrated material platform for realizing efficient thin-film photovoltaics, with the most well known Cu(In,Ga)Se 2 (CIGS)-based solar cells exceeding 23%. Several factors, including flexibility in tuning the absorber bandgap, enhanced surface treatments, and the electrically benign nature of common defects are responsible for the existing high performance and future promise in chalcopyrite-based photovoltaic devices. The introduction of Cu-poor phases (also known as ordered-vacancy compounds or OVCs) between the absorber and buffer layers in CIGS solar cells is known to enhance device performance; however, the overall properties and role of OVCs remain poorly understood. Using first principles calculations based on the density functional theory with screened hybrid functionals, here we explore the electronic structure and stability of OVCs and their band offsets with defect-free chalcopyrite layers in Cu- and Ag-based compounds (ABX 2 where A=Cu, Ag; B=In, Ga, Al; and X=S, Se). Using AB 3 X 5 and AB 5 X 8 stoichiometries as model OVC systems, we report on the variation of the bandgap with the A/B ratio and discuss the trends in other Cu- and Ag-based chalcopyrites beyond CuInSe 2 . We find that the valence and conduction bands are lower in energy in OVCs with respect to the parent ABX 2 chalcopyrite owing to a reduced p–d interaction between X and A atoms. We additionally perform device-level simulations to assess the implications of the results, finding that the valence band offsets of OVCs are favorable, while the conduction band offsets of chalcopyrites beyond CuInSe 2 -based absorbers may be detrimental in conventional solar cell device designs.

14 SOLAR ENERGY↗
DOE OSTI2021

Cartesian message passing neural networks for directional properties: Fast and transferable atomic multipoles

The message passing neural network (MPNN) framework is a promising tool for modeling atomic properties but is, until recently, incompatible with directional properties, such as Cartesian tensors. We propose a modified Cartesian MPNN (CMPNN) suitable for predicting atom-centered multipoles, an essential component of ab initio force fields. The efficacy of this model is demonstrated on a newly developed dataset consisting of 46 623 chemical structures and corresponding high-quality atomic multipoles, which was deposited into the publicly available Molecular Sciences Software Institute QCArchive server. We show that the CMPNN accurately predicts atom-centered charges, dipoles, and quadrupoles and that errors in the predicted atomic multipoles have a negligible effect on multipole–multipole electrostatic energies. The CMPNN is accurate enough to model conformational dependencies of a molecule’s electronic structure. This opens up the possibility of recomputing atomic multipoles on the fly throughout a simulation in which they might exhibit strong conformational dependence.

Glick, Zachary L. (ORCID:0000000309002849)↗
DOE OSTI2025

Artificial intelligence for advanced functional materials: exploring current and future directions

This perspective addresses the topic of harnessing the tools of artificial intelligence (AI) for boosting innovation in functional materials design and engineering as well as discovering new materials for targeted applications in energy storage, biomedicine, composites, nanoelectronics or quantum technologies. It gives a current view of experts in the field, insisting on challenges and opportunities provided by the development of large materials databases, novel schemes for implementing AI into materials production and characterization as well as progress in the quest of simulating physical and chemical properties of realistic atomic models reaching the trillion atoms scale and with near ab initio accuracy.

36 MATERIALS SCIENCE↗