SEARCH · Search NASA
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.
A prominent goal in quantum chemistry is to solve the molecular electronic structure problem for ground state energy with high accuracy.
Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.
Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.
First, high-resolution sub-Doppler infrared spectroscopic results for cyclopentyl radical (C 5 H 9 ) are reported on the α-CH stretch fundamental with suppression of spectral congestion achieved by adiabatic cooling to T rot ≈ 19(4) K in a slit jet expansion. Surprisingly, cyclopentyl radical exhibits a rotationally assignable infrared spectrum, despite 3N – 6 = 36 vibrational modes and an upper vibrational state density (ρ ≈ 40–90 #/cm –1 ) in the critical regime (ρ ≈ 100 #/cm –1 ) necessary for onset of intramolecular vibrational relaxation (IVR) dynamics. Such high-resolution data for cyclopentyl radical permit detailed fits to a rigid-rotor asymmetric top Hamiltonian, initial structural information for ground and vibrationally excited states, and opportunities for detailed comparison with theoretical predictions. Specifically, high level ab initio calculations at the coupled-cluster singles, doubles, and perturbative triples (CCSD(T))/ANO0, 1 level are used to calculate an out-of-plane bending potential, which reveals a C 2 symmetry double minimum 1D energy surface over a C 2v transition state. Here, the inversion barrier [V barrier ≈ 3.7(1) kcal/mol] is much larger than the effective moment of inertia for out-of-plane bending, resulting in localization of the cyclopentyl wavefunction near its C 2 symmetry equilibrium geometry and tunneling splittings for the ground state too small (<1 MHz) to be resolved under sub-Doppler slit jet conditions. The persistence of fully resolved high-resolution infrared spectroscopy for such large cyclic polyatomic radicals at high vibrational state densities suggests a “deceleration” of IVR for a cycloalkane ring topology, much as low frequency torsion/methyl rotation degrees of freedom have demonstrated a corresponding “acceleration” of IVR processes in linear hydrocarbons.
EuCd 2 As 2 is a remarkably complex magnetic semimetal that may behave as a topological insulator or host two pairs of Weyl points, depending on the growth conditions and the final magnetic state. Both antiferromagnetic (AFM) and ferromagnetic (FM) forms have been grown, and we show here, using 151 Eu Mössbauer spectroscopy, that the differences between the AFM and FM forms extend well beyond their ground state magnetic structures. Whereas the AFM form undergoes a conventional AFM → paramagnetic transition on warming, the FM form passes through a complex incommensurate modulated state before becoming paramagnetic.
A new iterative solver is proposed to efficiently calculate the ground state electronic structure in density functional theory calculations. This algorithm is particularly useful for simulating physical systems considered difficult to converge by standard solvers, in particular metallic systems. Here, the effectiveness of the proposed algorithm is demonstrated on various applications.
We use the entanglement negativity, a bipartite measure of entanglement in mixed quantum states, to study how multipartite entanglement constrains the real-space structure of the ground state wavefunctions of (2 + 1)-dimensional topological phases. We focus on the (Abelian) Laughlin and (non-Abelian) Moore-Read states at filling fraction ν = 1/m. We show that a combination of entanglement negativities, calculated with respect to specific cylinder and torus geometries, determines a necessary condition for when a topological state can be disentangled, i.e., factorized into a tensor product of states defined on cylinder subregions. This condition, which requires the ground state to lie in a definite topological sector, is sufficient for the Laughlin state. On the other hand, we find that a general Moore-Read ground state cannot be disentangled even when the disentangling condition holds.
Heavy-fermion family exhibits fascinating and often puzzling properties due to the presence of open-shell f ions and the complexity of the associated charge, orbital, and spin degrees of freedom. SmB 6 is a prototypical heavy-fermion compound that is electrically insulating but yet it displays quantum oscillations, which are a telltale signature of the metallic state. Adding to the enigma is the possibility that SmB 6 is a topological Kondo insulator. Here, by treating the spin degree of freedom on an equal footing with other degrees of freedom using the parameter-free strongly-constrained and appropriately-normed (SCAN) density functional, we explore the ground-state electronic structure of SmB6. A number of competing magnetic phases lying very closely in energy are found, indicating the key role of spin fluctuations in the material. The computed band structure, crystal-field splittings in the f-electron complex, the heavy effective electron mass at the Fermi energy, and the large specific heat are all in good agreement with the corresponding experimental results. In particular, our predicted FS explains the experimentally observed bulk quantum oscillations as well as the low electrical conductivity of SmB 6 . The topological Kondo state of SmB6 is shown to be robust regardless of its magnetic configuration. The excellent performance of SCAN in heavy-fermion systems is explained in terms of its ability to treat self-interaction errors and symmetry breaking within the framework of the density functional theory. Our study provides a new approach for modeling heavy-fermion materials
The multipartite entanglement structure for the ground states of two-dimensional (2D) topological phases is an interesting albeit not well-understood question. Utilizing the bulk-boundary correspondence, the calculation of tripartite entanglement in 2D topological phases can be reduced to that of the vertex state, defined by the boundary conditions at the interfaces between spatial regions. In this paper, we use the conformal interface technique to calculate entanglement measures in the vertex state, which include area-law terms, corner contributions, and topological pieces, and a possible additional order-one contribution. Furthermore, this explains our previous observation of the Markov gap ℎ = $\frac{c}{3}$ ln2 in the three-vertex state, and generalizes this result to the p-vertex state, general rational conformal field theories, and more choices of subsystems. Finally, we support our prediction by numerical evidence, finding precise agreement.
The electronic structure and magnetic properties of a newly predicted stable ternary compound La 18 Co 28 Pb 3 are studied using electronic structure analysis. The ground state of this compound is ferromagnetic, with three positions of nonequivalent magnetic Co atoms. A strong dependence of magnetic properties on volume shows that this system is situated near the point of magnetic instability. The coexistence of high- and low-spin ferromagnetic states as a function of volume near equilibrium was discovered. A corresponding spin tunneling splitting was estimated. The stability of the theoretically predicted magnetic ground state was tested by varying the Hubbard parameter. The thermal spin fluctuations were added to estimate the paramagnetic moment and a Curie temperature. Finally, the necessity of experimental verification of the obtained results is emphasized.
The CLAS Physics Database (CLAS DB) has been developed in collaboration between SINP at MSU and Hall B at Jefferson Lab. CLAS DB contains most experimental results on inclusive, semi-inclusive, and exclusive reactions on proton and nuclei obtained with the CLAS detector. The substantial part of these experimental results represents the only data sets available in the world. The results and analysis tools provided by the CLAS DB are of particular importance for studies of the structure of the ground and excited hadron states, offering insight into strong nonperturbative QCD dynamics, which underlie the generation of hadrons from quarks and gluons.
Abstract High-frequency seismic data on Mars are dominated by wind-generated lander vibrations, which are radiated partially to the subsurface. Autocorrelation functions (ACFs) of seismic data on Mars filtered between 1 and 5 Hz show clear phases at ∼1.3, ∼2.6, and ∼3.9 s. Daily temporal changes of their arrival times (dt/t) correlate well with the daily changes of ground temperature, with ∼5% daily variation and ∼50 min apparent phase delay. The following two mechanisms could explain the observations: (1) the interference of two predominant spectral peaks at ∼3.3 and ∼4.1 Hz, assumed to be both lander resonance modes, generate the apparent arrivals in the ACFs; (2) the interference of the lander vibration and its reflection from an interface ∼200 m below the lander generate the 3.3 Hz spectral peak and ∼1.3 s arrival in the ACFs. The driving mechanism of the resolved dt/t that most likely explains the ∼50 min delay is thermoelastic strain at a near-surface layer, affecting the lander–ground coupling and subsurface structures. The two outlined mechanisms suggest, respectively, up to ∼10% changes in ground stiffness at 1–5 Hz and ∼15% velocity changes in the top ∼20 m layer. These are upper bound values considering also other possible contributions. The presented methodology and results contribute to analysis of ACFs with limited data and the understanding of subsurface materials on Mars.
In this work, a total of four solar panel layouts, applied to Agrivoltaics, are investigated. This study presents two hypotheses: 1. Prefabrication - Lower total cost can be achieved by constructing a pre-assembled, pre-wired solar array in a factory, with automation, then rapidly deploying it in the field. 2. Tension Structures - Using tension structures in a way to suspend an overhead solar array also reduces costs when compared to conventional support methods. Designs are presented for frames to hold the panels. Configurations are advanced to connect the frames into arrays. Supports and foundations to hold the arrays in the desired layout are calculated, selected and presented. Cost models are prepared for the four layouts, which include definition of the business enterprise, the manufacturing operation, and the required facility for production at scale. Cost models include, detail material and task/labor take-offs for both manufacturing and installation. Comparisons with conventional system adapted to Agrivoltaics reveal that the two hypotheses do not hold true. A prefabricated array requires additional material that is not needed in the base system and the cost of this additional material is not outweighed by the cost savings of reduced field installation. The cost of wire rope for a tension structure does result in an economy of material, however, the cost of end connections and tensioning devices adds significantly to the cost of the overall tension structure. Foundation loads with tension structures, especially the large wind uplift seen with a fixed tilt solar array, impose significant constraints and high costs. The following conclusions are recommended for further consideration: 1. Best Case Agrivoltaic Configuration - The best value/lowest cost configuration for agrivoltaics is the installation of solar panels mounted on a single axis tracker in either a 2-in-portrait or 2-in landscape arrangement. This is useful for grazing lands and for staple crops which require high light levels. These applications represent the vast majority of the potential agrivoltaic market. Important elements in this configuration are the specific geometric layout coordinated with the various field operations and a control system linked between the farm equipment and single axis tracker. The control system is to monitor location and orientation, then actively tilt the solar panels or brake the field equipment to provide clearance and avoid collisions. 2. Long Span Applications - The use of tension structures for an overhead fixed tilt array is beneficial where long spans are a necessity, especially where the base support in the project location has consolidated rock near the surface or some other solid structure. One application is spanning wide irrigation canals. Another may be as an installation at the top deck of parking structures. 3. Direct Ground Mount Configuration – Having skids at the base of a pre-fabricated array that extend when the array is deployed, could be a successful design. This could be a useful strategy for dense, flat, ground-mount configurations like that employed by a couple of successful commercial operations (5G Maverick and Erthos). Arrays would be orientated north and south, deployed very low on prepared grade. Arrays would be aligned closely next to each other and anchored with small ground screws or simply with ballast. This same product arrangement could be very useful for rapid deployment and set-up of solar arrays for temporary deployments, disaster response and the military operations
The COVID-19 pandemic underscored the promise of monoclonal antibody-based prophylactic and therapeutic drugs, especially where protective candidates can be rapidly identified or developed for emerging biothreats and escape variants. Current cutting-edge technologies for this purpose still rely on pathogen-exposed convalescent volunteers and a large screening effort to find a proverbial needle in a haystack. Computational design of protective antibodies based on pre-existing templates skips those requirements and allows for greater control over the breadth and target epitope, while also co-optimizing for potency and developability or other biophysical characteristics. We approached this problem by building and expanding an in vitro experimental rapid antibody production and characterization pipeline to support development of an autonomous, closed loop, active learning software system based on structural simulation and ground truth experimental data to design and evaluate antibody antigen interactions. Starting from early in the pandemic, we targeted SARS-CoV-2, the causative agent of COVID-19, by re-purposing neutralizing antibodies against SARS-CoV-1 that had been identified in the wake of that outbreak in the early 2000’s. We successfully re-targeted three different anti-SARS-CoV-1 antibodies to neutralize SARS-CoV-2 in vitro, where the antibodies were generated externally and tested through conventional binding and neutralization assays internally or with collaborators. As antibodies were identified from the blood of humans infected with SARS-CoV-2, we shifted to collaborate with academic partners to develop improved versions of their human-derived antibodies. This work reached its most important stage in rapid response to the emergence of the Omicron variant of concern (VOC) in late 2021. In a matter of weeks, enabled by on demand innovation to our screening pipeline, we computationally designed and experimentally characterized derivative antibodies of COV2-2130, one of two antibodies from Vanderbilt that form the basis of the AstraZeneca Evusheld prophylactic drug product. This drug product suffers a serious loss of efficacy against Omicron BA.1 and BA.1.1, the first Omicron strains. Due to tight integration of computational design and experimental evaluation, we were able to identify a pair of designs with potent neutralization of the main targets Omicron BA.1 and BA.1.1; but also the earlier Delta variant, and subsequent Omicron strains including BA.2, BA.4, BA.5, and BA.2.75, demonstrating that our multi-target design process can, by its nature, produce robust antibody designs that strictly improve over the parental antibody. These results, recognized by a 2022 Director’s Science and Technology award, have enabled the follow-on GUIDE program, to commence in FY23. While earlier design campaigns were substantially outsourced, we have engineered better and faster processes internally to better compliment, calibrate, and speed computational designs. As part of the follow-on GUIDE program, we will stand up a rapid and high-throughput antibody production and characterization facility staffed with the expertise and capabilities to foster our current collaboration across PLS and ENG as well as other partnerships toward computational design of biologics.
The photocycloaddition (PCA) of chalcones represents an important reaction pathway for accessing substituted cyclobutanes, which is a molecular framework with utility in synthetic chemistry, materials science, and medicine. In the past, our group has demonstrated the utility of the large cavity of γ-CD as a container for encapsulating two photo reactants for directing the PCA of several classes of aryl alkenes with high stereo- and regioselectivity: the cavitand-mediated photodimerization (CMP) approach. The CMP of chalcones reported in this work further demonstrates the effectiveness of this approach as high yields of dimers were observed in the photoreactions, while they were non-reactive in the solid state and yielded only the isomerization product in homogeneous media. The γ-CD CMP of chalcones yielded predominantly dimerized products in very good to high yields (>70%), composed of a mixture of three dimers in different proportions with syn HH as the major product. Computational analysis of the ground state complex structures revealed a strong correlation between the stability of the complex and predominance of the stereoisomer in the mixture. Further insights were deduced from temperature-dependence studies, which showed a shift in dimer selectivity tending towards a single stereoisomer.
Trichomes show 47 morphological phenotypes, while literature reports only two root hair phenotypes in all plants. However, could hair-like structures exist below-ground in a similar wide range of morphologies like trichomes? Genetic mutants and root hair stress phenotypes point to the possibility of uncharacterized morphological variation existing belowground. For example, such root hairs in Arabidopsis (Arabidopsis thaliana) can be wavy, curled, or branched. We found hints in the literature about hair-like structures that emerge before root hairs belowground. As such, these early emerging hair structures can be potential exceptions to the contrasting morphological variation between trichomes and root hairs. Here, in this work, we show a previously unreported ‘hooked’ hair structure growing below-ground in common bean. The unique ‘hooking’ shape distinguishes the ‘hooked hair’ morphologically from root hairs. Currently, we cannot fully characterize the phenotype of our observation due to the lack of automated methods for phenotyping root hairs. This phenotyping bottleneck also handicaps the discovery of more morphology types that might exist below-ground as manual screening across species is slower than computer-assisted high-throughput screening.
Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.
QCD gives rise to a rich spectrum of excited baryon states. Understanding their internal structure is important for many areas of nuclear physics, such as nuclear forces, dense matter, and neutrino-nucleus interactions. Generalized parton distributions (GPDs) are an established tool for characterizing the QCD structure of the ground-state nucleon. They are used to create 3D tomographic images of the quark/gluon structure and quantify the mechanical properties such as the distribution of mass, angular momentum, and forces in the system. Transition GPDs extend these concepts to N → N* transitions and can be used to characterize the 3D structure and mechanical properties of baryon resonances. They can be probed in high-momentum-transfer exclusive electroproduction processes with resonance transitions e + N → e' + M + N*, such as deeply-virtual Compton scattering (M = γ) or meson production (M = π, K, etc.), and in related photon/hadron-induced processes. This White Paper describes a research program aiming to explore baryon resonance structure with transition GPDs. This includes the properties and interpretation of the transition GPDs, theoretical methods for structures and processes, first experimental results from JLab 12 GeV, future measurements with existing and planned facilities (JLab detector and energy upgrades, COMPASS/AMBER, EIC, EicC, J-PARC, LHC ultraperipheral collisions), and the theoretical and experimental developments needed to realize this program.