Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytical mathematics”

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

Enhanced Twist-Averaging Technique for Magnetic Metals: Applications Using Quantum Monte Carlo

In this study, we propose an improved twist-averaging (TA) scheme for quantum Monte Carlo methods that use converged Kohn–Sham or Hartree–Fock orbitals as the reference. This TA technique is tailored to sample the Brillouin zone of magnetic metals, although it naturally extends to nonmagnetic (NM) conducting systems. The proposed scheme aims to reproduce the reference magnetization and achieves charge neutrality by construction, thus avoiding the large energy fluctuations and the postprocessing needed to correct the energies. It shows the most robust convergence of total energy and magnetism to the thermodynamic limit (TDL) when compared to four other TA schemes. Diffusion Monte Carlo applications are shown on NM Al and ferromagnetic α-Fe. The cohesive energy of Al in the TDL shows an excellent agreement with the experimental result. Furthermore, the magnetic moments in α-Fe exhibit rapid convergence with an increasing number of twists.

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

Efficient Treatment of Large Active Spaces through Multi-GPU Parallel Implementation of Direct Configuration Interaction

In this study, we have extended our graphical processing unit (GPU)-accelerated direct configuration interaction program to multiple devices, reducing iteration times for configuration spaces of 165 million determinants to only 3 s using NVIDIA P100 GPUs. Similar improvements in the one- and two-particle reduced density matrix formation allow for fast analytical energy gradients and electronic properties. Our parallel algorithm enables the calculation of arbitrarily large configuration spaces (limited only by available system memory), with iteration times of 13 min for an active space of 18 electrons in 18 orbitals (2.4 billion determinants) using six consumer grade NVIDIA 1080Ti GPUs. These advances enable routine molecular dynamics simulations, geometry optimizations, and absorption spectrum calculations for molecules with large configuration spaces, a task that has heretofore required massive computational effort. In this work, we demonstrate the utility of our program by generating the absorption spectrum for diphenyl acetylene at the floating occupation molecular orbital complete active space configuration interaction level of theory. Lastly, several active spaces were investigated to assess the dependence of spectral features on orbital space dimension.

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

Photoelectron Spectroscopic and ab Initio Computational Studies of the Anion, HThO –

The synergetic combination of anion photoelectron spectroscopy and high-level relativistic coupled-cluster calculations was employed to study the anion, HThO-. The atomic connectivity of this anion was found to be HThO- and not ThOH-. Vibrational and electronic energy spacings in the HThO- photoelectron spectrum were measured and calculated, with good agreement between them being found. Computations yielded electronic energies and equilibrium structures as well as enabling orbital analyses. The adiabatic electron affinity of HThO was determined to be 1.297 ± 0.035 eV.

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

Toward the Minimal Floating Operation Count Cholesky Decomposition of Electron Repulsion Integrals

As quantum chemistry calculations deal with molecular systems of increasing size, the memory requirement to store electron-repulsion integrals (ERIs) greatly outpaces the physical memory available in computing hardware. Here, the Cholesky decomposition of ERIs provides a convenient yet accurate technique to reduce the storage requirement of integrals. Recent developments of a two-step algorithm have drastically reduced the memory operation (MOP) count, leaving the floating operation (FLOP) count as the last frontier of cost reduction in the Cholesky ERI algorithm. In this report, we introduce a dynamic integral tracking, reusing, and compression/elimination protocol embedded in the two-step Cholesky ERI method. Benchmark studies suggest that this technique becomes particularly advantageous when the basis set consists of many computationally expensive high-angular-momentum basis functions. With this dynamic ERI improvement, the Cholesky ERI approach proves to be a highly efficient algorithm with minimal FLOP and MOP count.

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

Projector-Based Quantum Embedding for Molecular Systems: An Investigation of Three Partitioning Approaches

Projector-based embedding is a relatively recent addition to the collection of methods that seek to utilize chemical locality to provide improved computational efficiency. This work considers the interactions between the different proposed procedures for this method and their effects on the accuracy of the results. The interplay between the embedded background, projector type, partitioning scheme, and level of atomic orbital (AO) truncation are investigated on a selection of reactions from the literature. The Huzinaga projection approach proves to be more reliable than the level-shift projection when paired with other procedural options. Active subsystem partitioning from the subsystem projected AO decomposition (SPADE) procedure proves slightly better than the combination of Pipek-Mezey localization and Mulliken population screening (PMM). Along with these two options, a new partitioning criteria is proposed based on subsystem von Neumann entropy and the related subsystem orbital occupancy. This new method overlaps with the previous PMM method, but the screening process is computationally simpler. Finally, AO truncation proves to be a robust option for the tested systems when paired with the Huzinaga projection, with satisfactory results being acquired at even the most severe truncation level.

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

Comparison of Machine Learning Approaches for Prediction of the Equivalent Alkane Carbon Number for Microemulsions Based on Molecular Properties

The chemical properties of oils are vital in the design of microemulsion systems. The hydrophilic–lipophilic difference equation used to predict microemulsions’ phase behavior expresses the oils’ physiochemical properties as the equivalent alkane carbon number (EACN). The experimental determination of EACN requires knowledge of the temperature dependence of the microemulsion system and the effects of different surfactant concentrations. Thus, the experimental determination is time-intensive and tedious, requiring days to months for proper separations. Furthermore, the experiments require high purity of chemicals because microemulsions are sensitive to impurities. Our work focuses on the quick and reliable predictions of the EACN with machine learning (ML) models. Due to the immaturity of ML chemical predictions, we compare three graph neural networks (GNNs) and a gradient-boosted tree algorithm, known as XGBoost. The GNNs use the molecular structures represented as simplified molecular-input line-entry system (SMILES) codes for the initial input, which allows us to assess whether geometry optimization is necessary for reliable results. The XGBoost model also begins with the SMILES representations of the molecules but uses molecular descriptors instead of geometry optimizations. As a result, the best model tested (crystal graph convolutional neural network with Merck molecular force field-94) has an error of 1.15 EACN units of the true EACN for unknown data with the errors skewed toward zero and an R² score of 0.9

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

Accurate Prediction of pKb in Amines: Validation of the CAM-B3LYP/6-311+G(d,p)/SMD Model

Amines play several key roles in chemistry and biology and are involved in numerous industrial processes, often with significant economic impacts. Recently, amines are also garnering interest as catalysts for polymer synthesis and for CO 2 fixation, incentivizing the need to rapidly design and screen new amino compounds. Hence, developing reliable methods to predict their physicochemical properties, e.g., the base dissociation constant (pKb), is pivotal. Here, a density functional theory (DFT)-based approach was employed to compute the pKb of substituted amines, exploring the impact of several key parameters, including (i) the number of explicit water molecules at the reaction center, (ii) the van der Waals (vdW) surface, and (iii) solvent polarizability. In previous work, it was determined that including two explicit water molecules at the reaction center resulted in highly accurate pKb estimates for primary amines. Here, we find that including a third water molecule at the reaction center is essential for accurate pKb for secondary and tertiary amines. The revised methodology was then applied to a wider selection of amines, obtaining a minimum average error (MAE) < 0.4. In conclusion, this result represents an extension of our “easy-to-use method,” a simple and direct DFT approach exploiting CAM-B3LYP/SMD/6-311G+(d,p) to compute pKb without post facto modifications.

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

Infrared Spectroscopy of Li + Solvation in EmimBF 4 and in Propylene Carbonate: Ab Initio Molecular Dynamics and Experiment

Infrared (IR) spectra of solutions of the lithium salt LiBF 4 in the ionic liquid 1-ethyl-3-methylimidazolium tetrafluoroborate (EmimBF 4 ) and in the organic solvent propylene carbonate (PC) are studied via infrared spectroscopy and ab initio molecular dynamics (AIMD) simulations. The measurements show that the major effects of LiBF 4 in both solutions, compared to their neat counterparts, are the appearance of a new broad band in the 300-450 cm -1 frequency region and a broadening of the IR structure in the 900-1200 cm -1 region with the development of a new peak at 980 cm -1 . Computational analysis indicates that hindered translational motions of Li + in its solvation cage are mainly responsible for the former, while the latter is due to Li + -induced structural changes and accompanying vibrational frequency shifts of constituent ions and molecules of the solutions. In addition, molecular motions in these and lower frequency regions are generally correlated, disclosing the collective nature of the vibrational dynamics, which involve multiple ions/molecules. Herein, detailed analysis of these features via AIMD simulations of the spectrum and its components arising from auto- and cross-correlations of motions of constituent molecular species, combined with generalized normal modes of the solutions and normal modes of small Li + -containing clusters, is presented. Furthermore, other minor spectral changes caused by the lithium salt as well as the interaction-induced effect on IR spectra are also discussed.

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

Experimental Investigation of Thermodynamic Stabilization in Boron Imidazolate Frameworks (BIFs) Synthesized by Mechanochemistry

This study experimentally explores the energetics for the formation of boron–imidazolate frameworks (BIFs), which are synthesized by mechanochemistry. The topologically similar frameworks employ the same tetratopic linker based on tetrakis(imidazolyl)boric acid but differ in the monovalent cation metal nodes. This permits assessment of the stabilizing effect of metal nodes in frameworks with sodalite (SOD) and diamondoid (dia) topologies. The enthalpy of formation from endmembers (metal oxide and linker), which define thermodynamic stability of the structures, has been determined by use of acid solution calorimetry. The results show that heavier metal atoms in the node promote greater energetic stabilization of denser structures. Overall, in BIFs the relation between cation descriptors (ionic radius and electronegativity) and thermodynamic stability depends on framework topology. Thermodynamic stability increases with the metallic character of the cation employed as the metal node, independent of the framework topology. In conclusion, the results suggest unifying aspects for thermodynamic stabilization across MOF systems.

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

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

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

Reducing Qubit Requirements while Maintaining Numerical Precision for the Variational Quantum Eigensolver: A Basis-Set-Free Approach

We present a basis-set-free approach to the variational quantum eigensolver using an adaptive representation of the spatial part of molecular wave functions.Our approach directly determines system-specific representations of qubit Hamiltonians while fully omitting globally defined basis sets. In this work, we use directly determined pair-natural orbitals on the level of second-order perturbation theory. This results in compact qubit Hamiltonians with high numerical accuracy. We demonstrate initial applications with compact Hamiltonians on up to 22 qubits where conventional representation would for the same systems require 40–100 or more qubits. Here, we further demonstrate reductions in the quantum circuits through the structure ofthe pair-natural orbitals.

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

Emulating Quantum Entangled Biphoton Spectroscopy Using Classical Light Pulses

We show that for a class of quantum light spectroscopy (QLS) experiments using n = 0, 1, 2, ··· classical light pulses and an entangled photon pair (a biphoton state) where one photon acts as a reference without interacting with the matter sample, identical signals can be obtained by replacing the biphotons with classical-like coherent states of light, where these are defined explicitly in terms of the parameters of the biphoton states. An input-output formulation of quantum nonlinear spectroscopy is used to prove this equivalence. We demonstrate the equivalence numerically by comparing a classical pump–quantum probe experiment with the corresponding classical pump–classical probe experiment. This analysis shows that understanding the equivalence between entangled biphoton probes and carefully designed classical-like coherent state probes leads to quantum-inspired classical experiments that yield equivalent signals and provides insights for the future design of QLS experiments that could provide a true quantum advantage.

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

Scaling Theory of Neutral Sequence-Specific Polyampholytes

The primary sequence of charged monomers in a polyampholyte has a profound effect on its conformational properties. In this work, theory and simulations are used to predict how the sequence influences the single-chain behavior of globally neutral polyampholytes (PAs) under salt-free conditions. We consider PAs with Markov statistics of charges, where charge “blockiness” is defined by the correlation λ along the chain. These PAs cover a wide spectrum of primary sequences ranging from alternating (λ = –1) to ideally random (λ = 0) to diblock (λ → 1). In a theta solvent, sufficiently long PAs of any primary sequence form globules, but their internal structure, as well as the position and sharpness of the coil–globule transition, crucially depends on λ. Further, using scaling arguments, we demonstrate that the excess charge of the globule correlation volume (blob) increases with increasing λ, thereby inducing stronger Coulomb attractions between neighboring blobs of opposite charge and leading to higher density PA globules. Scaling predictions are supported by the results of the random phase approximation (RPA) and coarse-grained molecular dynamics simulations. A generalization to good and poor solvents is presented that enables the construction of the resulting scaling diagram of PA conformations in terms of solvent quality and charge blockiness, thereby providing a comprehensive molecular understanding of the role of the PA primary sequence in the chain structure.

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

Amidinate- and Dithiolene-Based Silicon Complexes

Reactions of the amidinato-silylene chloride PhC( t BuN) 2 SiCl (1) with imidazole-based dithione dimer 2, lithium dithiolene radical 3, and dithiolate dimer 4 result in the synthesis of a series of silicon complexes 5-7, respectively, containing both amidinato and dithiolene ligands. 7 is the first structurally characterized silicon(II) dithiolene complex. The structural and bonding characteristics of 5-7 have been probed by both experimental and theoretical methods.

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

Stochasticity in the Switching of Nanodisks for Probabilistic Computing

Stochasticity in magnetic nanodevices is an essential characteristic for harnessing these devices to computing based on population coding or the building blocks of probabilistic computing, p-bits. A magnetic tunneling junction (MTJ) consisting of a patterned magnetic element is considered a promising computing unit in the concept of artificial neurons and p-bits. A comprehensive understanding of the stochasticity in the switching of patterned magnetic elements is crucial for realizing MTJ-based probabilistic computing technology. In the present work, the stochastic behavior in the switching process of a perpendicularly magnetized Co/Pt disk within an array was directly observed utilizing full-field soft X-ray microscopy. Within 50 repeated hysteretic cycles, the stochastic magnetization switching of individual Co/Pt disks within disk arrays is identified. We found that the stochasticity in the magnetization switching of disks considerably depends on the disk size. The stochasticity initially decreases as the disk radius gets bigger from 125 to 375 nm (region I), then increases with further enlarging the disk size to 625 nm (region II). The variance of thermal fluctuation relevant to the disk size and the multilevel switching within a disk are severely involved in the observed size-dependent stochasticity. Furthermore, this work provides the way for controlling the stochasticity in the switching of nanopatterned elements, which is a key aspect of MTJ-based probabilistic computing.

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

Machine Learning-Guided Identification of PET Hydrolases from Natural Diversity

The enzymatic depolymerization of poly(ethylene terephthalate) (PET) is emerging as a leading chemical recycling technology for waste polyester. As part of this endeavor, new candidate enzymes identified from natural diversity can serve as useful starting points for enzyme evolution and engineering. In this study, we improved upon HMM searches by applying an iterative machine learning strategy to identify 400 putative PET-degrading enzymes (PET hydrolases) from naturally occurring homologs. Using high-throughput (HTP) experimental techniques, we successfully expressed and purified >200 enzyme candidates and assayed them for PET hydrolysis activity as a function of pH, temperature, and substrate crystallinity. From this library, we discovered 91 previously unknown PET hydrolases, 35 of which retain activity at pH 4.5 on crystalline material, which are conditions relevant to developing more efficient commercial processes. Notably, four enzymes showed equal to or higher activity than LCC-ICCG, a benchmark PET hydrolase, at this challenging condition in our screening assay, and 11 of which have pH optima <7. Using these data, we identified regions of PETases statistically correlated to activity at lower pH. We additionally investigated the effect of condition-specific activity data on trained machine learning predictors and found a precision (putative hit rate) improvement of up to 30% compared to a Hidden Markov Model alone. Our findings show that by pointing enzyme discovery toward conditions of interest with multiple rounds of experimental and machine learning, we can discover large sets of active enzymes and explore factors associated with activity at those conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗