Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

X-Ray Imagery as the Record of All Data of Interest in Hypervelocity Impact Fragment Studies

Laboratory study of hypervelocity spacecraft fragmentation has traditionally involved the collection and analysis of fragments that were caught in deceleration material surrounding the impact. This process has typically involved the disintegration of the catchment material either through chemical dissolution, or through physical excavation to recover the fragments. Due to the scale of the three impact tests—the Satellite Orbital Debris Characterization Impact Test (SOCIT), the DebriSat satellite impact test, and the DebrisLV launch vehicle impact test—the latter two using more than 12 cubic meters of polyurethane foam to capture the fragments, hese projects have used x-ray imagery to precisely locate and thus, to more efficiently extract fragments in the soft-catch material. Three years into the DebriSat fragment extraction process, a side study was initiated to explore what additional information could be discerned from the x-rays, with significant results. This study was instrumental to a rapid replacement and retooling as the project was forced to replace the x-ray system around which the extraction process had been based. The revised process continues to map the debris for extraction. The project has, in parallel, systematically addressed the limits/tolerances of what x-rays can reveal about size, shape, density, mass, velocity, energy, and deformation/damage of the fragment during the deceleration in the catchment material while replicating the original extraction mapping function. All of these features have been optimized or have sufficient understanding to characterize the basic factors that will define a complete data set extracted solely from x-ray imagery. It is an ideal time to develop such a process, with extracted fragments providing “ground truth” against image-only data, and abundant available imagery of the same fragments under both the prior and replacement x-ray technologies, which have several fundamentally different characteristics. This paper addresses the types and quality of hypervelocity fragmentation data that can be and has been extracted from x-rays. It further addresses the question of whether and under what circumstances future hypervelocity experiments can use x-ray methods to largely—or to completely—avoid the extraction process in recording all appropriate results. Lastly, this paper addresses lessons learned and how future efforts can be further optimized.

John B. Bacon↗
DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computing x-ray absorption spectra from linear-response particles atop optimized holes

State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn–Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. Herein, this is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ~0.6 eV rms error for the K-edges of C–F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.

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

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗
DOE OSTI2026

Correction to “Toward Chemical Accuracy for Chemi- and Physisorption with an Efficient Density Functional”

In the original manuscript, the value of parameter c was not explicitly provided. We wish to clarify that c was refitted to the hydrogen atom exchange energy, using the updated value of κ. With the updated κ = 0.2501, the optimized value of the parameter is c = 0.391596228. These values supersede any implicit assumptions about c in the original manuscript, and readers are encouraged to use these parameters for reproducing the results presented therein.

Adsorption↗
DOE OSTI2020

Time-Resolved XUV Opacity Measurements of Warm Dense Aluminum

The free-free opacity in plasmas is fundamental to our understanding of energy transport in stellar interiors and for inertial confinement fusion research. However, theoretical predictions in the challenging dense plasma regime are conflicting and there is a dearth of accurate experimental data to allow for direct model validation. Here we present time-resolved transmission measurements in solid-density Al heated by an XUV free-electron laser. We use a novel functional optimization approach to extract the temperaturedependent absorption coefficient directly from an oversampled pool of single-shot measurements, and find a pronounced enhancement of the opacity as the plasma is heated to temperatures of order of the Fermi energy. Plasma heating and opacity enhancement are observed on ultrafast timescales, within the duration of the femtosecond XUV pulse. We attribute further rises in the opacity on ps timescales to melt and the formation of warm dense matter

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Inelastic Neutron Scattering from Thin Film Biaxially Oriented Polyethylene Terephthalate

Recent interest in emerging processes for polymer manufacturing and bio-based chemistries for direct chemical recycling/upcycling has motivated new research focused on a deeper understanding of atomic-scale polymer properties and how they influence macroscopic phenomena. Uncovering the fundamental properties of polymers that give rise to macroscopic behavior could enable new pathways for improved recyclability or utilization of alternative “greener” polymer analogues. In this study, the neutron vibrational spectrum was measured for a film of biaxially oriented polyethylene terephthalate (BoPET) using inelastic neutron scattering (INS), to investigate the relationship between the structure and dynamics of a widely used polymer. Compared to conventional spectroscopic techniques, the use of INS is advantageous for polymeric materials due to the absence of selection rules (i.e., all transitions are allowed), broad-band energy range, and considerable sensitivity to hydrogen modes. In order to distinguish the vibrational modes caused by trans and gauche rotational isomerism, the normal modes of vibration were calculated from a density functional theory-optimized structure of crystalline PET (cPET), representative of the all-trans state, and compared with INS from “highly crystalline” PET powder. Although in- and out-of-plane wagging of hydrogens on the ring structure exhibit significant contribution to both BoPET and cPET spectra, the wagging, rocking, and twisting modes of hydrogen on the ethylene glycol group are, in most cases, conformation-specific. These results were further rationalized by investigating the role of hyperconjugation in stabilizing both conformations using the natural bond order method. Here, through comparison of experimental and calculated INS results, this work provides the fundamental basis for discovering the role of structure and dynamics in shaping the macroscopic properties of PET and polymer analogues.

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

Data-driven approach to parameterize SCAN + U for an accurate description of 3 d transition metal oxide thermochemistry

Semilocal density-functional theory (DFT) methods exhibit significant errors for the phase diagrams of transition-metal oxides that are caused by an incorrect description of molecular oxygen and the large self-interaction error in materials with strongly localized electronic orbitals. Empirical and semiempirical corrections based on the DFT + U method can reduce these errors, but the parameterization and validation of the correction terms remains an on-going challenge. Here we develop a systematic methodology to determine the parameters and to statistically assess the results by considering interlinked thermochemical data across a set of transition metal compounds. We consider three interconnected levels of correction terms: (1) a constant oxygen binding correction, (2) Hubbard-U correction, and (3) DFT/DFT + U compatibility correction. The parameterization is expressed as a unified optimization problem. We demonstrate this approach for 3d transition metal oxides, considering a target set of binary and ternary oxides. With a total of 37 measured formation enthalpies taken from the literature, the dataset is augmented by the reaction energies of 1710 unique reactions that were derived from the formation energies by systematic enumeration. To ensure a balanced dataset across the available data, the reactions were grouped by their similarity using clustering and suitably weighted. The parameterization is validated using leave-one-out cross validation (CV), a standard technique for the validation of statistical models. We apply the methodology to the strongly constrained and appropriately normed (SCAN) density functional. Based on the CV score, the error of binary (ternary) oxide formation energies is reduced by 40% (75%) to 0.10 (0.03) eV/atom. A simplified correction scheme that does not involve SCAN/SCAN + U compatibility terms still achieves an error reduction of 30% (25%). The method and tools demonstrated here can be applied to other classes of materials or to parameterize the corrections to optimize DFT + U performance for other target physical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Towards atom-level understanding of metal oxide catalysts for the oxygen evolution reaction with machine learning

Abstract Green hydrogen production is crucial for a sustainable future, but current catalysts for the oxygen evolution reaction (OER) suffer from slow kinetics, despite many efforts to produce optimal designs, particularly through the calculation of descriptors for activity. In this study, we develop a dataset of density functional theory calculations of bulk and surface perovskite oxides, and adsorption energies of OER intermediates, which includes compositions up to quaternary and facets up to (555). We demonstrate that per-site properties of perovskite oxides such as Bader charge or band center can be tuned through element substitution and faceting, and develop a machine learning model that accurately predicts these properties directly from the local chemical environment. We leverage these per-site properties to identify promising perovskites with high theoretical OER activity. The identified design principles and promising materials provide a roadmap for closing the gap between current artificial catalysts and biological enzymes such as photosystem II.

Chemistry↗
DOE OSTI2020

Quantum Chemistry-Informed Active Learning to Accelerate the Design and Discovery of Sustainable Energy Storage Materials

Here we employed Density Functional Theory (DFT) to compute oxidation potentials of 1,400 homobenzylic ether molecules to search for the ideal sustainable redoxmer design. The generated data were used to construct an active learning model based on Bayesian optimization (BO) that targets candidates with desired oxidation potentials utilizing only a minimal number of DFT calculations. The active learning model demonstrated not only significant efficiency improvement over the random selection approach but also robust capability in identifying desired candidates in an untested set of 112,000 homobenzylic ether molecules. Our findings highlight the efficacy of quantum chemistry-informed active learning to accelerate the discovery of materials with desired properties from a vast chemical space.

25 ENERGY STORAGE↗
NASA NTRS2019

High-Performance Computing Optimization for Aladyn – Adaptive Neural Network Molecular Dynamics Mini-Application

This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.

Yamakov, Vesselin I.↗
DOE OSTI2021

Dark Energy Survey Year 3 results: Optimizing the lens sample in a combined galaxy clustering and galaxy-galaxy lensing analysis

We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 metacalibration sample for the sources. We explore easily reproducible selections based on magnitude cuts in i-band as a function of (photometric) redshift, zphot, and benchmark the potential gains against those using the well-established redMaGiC [E. Rozo et al., Mon. Not. R. Astron. Soc. 461, 1431 (2016)MNRAA40035-871110.1093/mnras/stw1281] sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MagLim sample, satisfies i<4zphot+18 and has ∼30% wider redshift distributions but ∼3.5 times more galaxies than redMaGiC. Assuming a wCDM model (i.e. with a free parameter for the dark energy equation of state) and equivalent scale cuts to mitigate nonlinear effects, this leads to 40% increase in the figure of merit for the pair combinations of Ωm, w, and σ8, and gains of 16% in σ8, 10% in Ωm, and 12% in w. Similarly, in ΛCDM, we find an improvement of 19% and 27% on σ8 and Ωm, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, i<22.2, has ∼7 times more galaxies and ∼20% wider redshift distributions compared to MagLim, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Rational optimization of substituted α–MnO 2 cathode for aqueous zinc–ion battery

Density functional theory (DFT) is utilized to explore the effects of increasing concentrations of vanadium (V) and chromium (Cr) substitution on the discharge of α-MnO 2 cathode in hydrated zinc-ion batteries (ZIBs). During H + -intercalation and Zn 2+ -intercalation, Cr-substitution proves to be more effective than V-substitution in promoting discharge behaviors. Transitioning from Mn 0.875 Cr 0.125 O 2 , Mn 0.75 Cr 0.25 O 2 to Mn 0.625 Cr 0.375 O 2 is found to consistently enhance the discharge voltage, along with improved tunnel structure retention and volume expansion suppression. In comparison, the promoting effect of increasing V-substitution is relatively small at initial discharge stages and leads to degradation at later stages, primarily due to an increased concentration of unstable Mn 2+ ion. The superior effect of Cr-substitution is attributed to the unique atomic and electronic structures of substituted Cr 4+ and reduced Cr 3+ ions during discharge. These ions serve as active electron acceptors to limit the formation of Mn 3+ and Mn 2+ ions, and as anchors to stabilize the α-MnO 2 framework and intercalated H + /Zn 2+ ions, respectively. Finally, our study highlights fine-tuning through substitution to enhance the performance of α-MnO 2 -based cathode materials in ZIBs.

25 ENERGY STORAGE↗
DOE OSTI2021

Net Electronic Charge as an Effective Electronic Descriptor for Oxygen Release and Transport Properties of SrFeO 3 -based Oxygen Sorbents

Perovskite oxides, as oxygen sorbents, exhibit excellent potential in thermochemical redox applications such as chemical looping air separation (CLAS), resulting from their excellent redox properties and high tunability. The structural and compositional flexibility of perovskites, while making these mixed oxides highly versatile, also poses challenges for their design and optimization. Moreover, the elevated operating temperature and dynamic oxygen partial pressure swings during redox processes add additional complications for oxygen sorbent design. From a fundamental standpoint, finding simple yet effective descriptors for the effect of cation substitutions on oxygen sorbents’ redox and oxygen transport properties are of significant importance for rational optimization of perovskite-based oxygen sorbents. In the present study, a series of A/B-site-doped SrFeO 3-δ perovskites are investigated using density functional theory (DFT) + U simulations. From these simulations, an effective electronic structural descriptor, the net electronic charge (Δe) of oxygen anions, is extracted based on its excellent correlations with vacancy formation energy, experimentally derived oxygen capacity, and oxygen vacancy migration barrier. Overall, these findings provide an effective tool to correlate the oxygen sorbents’ redox performances with a simple yet physically meaningful descriptor for rational design and optimization of perovskites in the context of CLAS.

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

Exploring the Effects of Node Topology, Connectivity, and Metal Identity on the Binding of Nerve Agents and Their Hydrolysis Products in Metal–Organic Frameworks

Recent studies have shown that metal–organic frameworks (MOFs) built from hexanuclear M(IV) oxide cluster nodes are effective catalysts for nerve agent hydrolysis, where the properties of the active sites on the nodes can strongly influence the reaction energetics. The connectivity and metal identity of these M6 nodes can be easily tuned, offering extensive opportunities for computational screening to predict promising new materials. Thus, we used density functional theory (DFT) to examine the effects of node topology, connectivity, and metal identity on the binding energies of multiple nerve agents and their corresponding hydrolysis products. By computing an optimization metric based on the relative binding strengths of key hydrolysis reaction species (water, agent, and bidentate-bound products), we predicted optimal M6 nodes for hydrolyzing specific nerve agent and simulant molecules, where our results are in qualitative agreement with observed experimental trends. This analysis highlighted the notion that no single metal or node topology is optimal for all possible organophosphates, suggesting that MOFs should be selected based on the agent of interest. Using the large amount of data generated from our DFT calculations, we then derived quantitative structure–activity relationship (QSAR) models to help explain the complex trends observed in the binding energies. Through linear regression, we identified the most important descriptors for describing the binding of nerve agents and their hydrolysis products to M6 nodes. These results suggested that both molecular and node properties, including both structural and chemical features, collectively contribute to the binding energetics. By performing a thorough statistical analysis, we showed that our QSAR models are capable of making quantitatively accurate binding energy predictions for nerve agents and their hydrolysis products in a wide variety of M(IV)-MOFs. In conclusion, the insights gained herein can be used to guide future experiments for the synthesis of MOFs with enhanced catalytic activity for organophosphate hydrolysis.

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

2D in-Plane Ordered MXene Nanosheets Derived from (Mo 2/3 Er 1/3 ) 2 AlC Rare-Earth i-MAX for Energy Storage Applications

MXenes have become one of the most versatile families of two-dimensional (2D) materials due to their high conductivity, hydrophilicity, and remarkable electrochemical performance. This has stimulated intense efforts to design and synthesize MXenes, including structurally unique in-plane ordered 2D MXenes called i-MXenes. Here, we have synthesized the quaternary rare earth (RE)-based i-MAX phase (Mo 2/3 Er 1/3 ) 2 AlC using an arc melting method, and the corresponding 2D i-MXene was then obtained through a LiF/HCl soft etching process. Literature studies have shown that Al and the RE element are etched out during the etching process, leading to the formation of pure vacancy-ordered Mo1.33C 2D i-MXene. However, our investigation reveals that upon exposure to a fluorine solution, the i-MAX phase forms RE fluoride impurities, which are challenging to remove through HCl−DI water washing and persist in the final product, resulting in impure Mo 1.33 C@Er i-MXene. These results were confirmed by various characterizations such as X-ray diffraction, Raman spectroscopy, X-ray photoelectron spectroscopy, and scanning transmission electron microscopy. Although the Mo 1.33 C@Er electrode showed a 24-fold increase in specific capacitance compared to its parent i-MAX phase, it still exhibited a high charge-transfer resistance arising from the insulating nature of RE fluoride byproducts, which adversely influence the overall capacitance behavior of the synthesized 2D Mo 1.33 C@Er i-MXenes. This study contributes to identifying pathways for the preparation of pure 2D i-MXenes from RE-based i-MAX phases and developing improved synthesis methods. With additional process optimization, the 2D i-MXene holds a strong potential for electrochemical energy storage applications. Additionally, the electronic structures of Mo 1.33 C were theoretically studied using first-principles density functional theory calculations, which revealed that pristine Mo 1.33 C is metallic, and this metallic nature is preserved even with −O, −F, and mixed functionalization.

chemical structure↗
DOE OSTI2024

Factors influencing quantum evaporation of helium from polar semiconductors from first principles

While there is much indirect evidence for the existence of dark matter (DM), to date it has evaded detection. Current efforts focus on DM masses over ~GeV—to push the sensitivity of DM searches to lower masses, new DM targets and detection schemes are needed. In this work, we focus on the latter—a novel detection scheme recently proposed to detect 10–100 meV phonons in polar target materials. Previous work showed that well-motivated models of DM can interact with polar semiconductors to produce an athermal population of phonons. This new sensing scheme proposes that these phonons then facilitate quantum evaporation of 3 He from a van der Waals film deposited on the target material. However, a fundamental understanding of the underlying process is still unclear, with several uncertainties related to the precise rate of evaporation and how it can be controlled. In this work, we use ab initio density functional theory calculations to compare the adsorption energies of helium atoms on a polar target material, sodium iodide, to understand the underlying evaporation physics. We explore the role of surface termination, monolayer coverage, and elemental species on the rate of He evaporation from the target material. Using this, we discuss the optimal target features for He-evaporation experiments and their range of tunability through chemical and physical modifications such as applied field and surface termination.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comparing the Expense and Accuracy of Methods to Simulate Atomic Vibrations in Rubrene

Atomic vibrations can inform about materials properties from hole transport in organic semiconductors to correlated disorder in metal–organic frameworks. Currently, there are several methods for predicting these vibrations using simulations, but the accuracy–efficiency tradeoffs have not been examined in depth. In this work, rubrene is used as a model system to predict atomic vibrational properties using six different simulation methods: density functional theory, density functional tight binding, density functional tight binding with a Chebyshev polynomial-based correction, a trained machine learning model, a pretrained machine learning model called ANI-1, and a classical forcefield model. The accuracy of each method is evaluated by comparison to the experimental inelastic neutron scattering spectrum. All methods discussed here show some accuracy across a wide energy region, though the Chebyshev-corrected tight-binding method showed the optimal combination of high accuracy with low expense. We then offer broad simulation guidelines to yield efficient, accurate results for inelastic neutron scattering spectrum prediction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗