Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multipole moments”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

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

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

Neat, Simple, and Wrong: Debunking Electrostatic Fallacies Regarding Noncovalent Interactions

Multipole moments such as charge, dipole, and quadrupole are often invoked to rationalize intermolecular phenomena, but a low-order multipole expansion is rarely a valid description of electrostatics at the length scales that characterize nonbonded interactions. This is illustrated by examining several common misunderstandings rooted in erroneous electrostatic arguments. First, the notion that steric repulsion originates in Coulomb interactions is easily disproved by dissecting the interaction potential for Ar2. Second, the Hunter-Sanders model of π–π interactions, which is based on quadrupolar electrostatics, is shown to have no basis in accurate calculations. Third, curved “buckybowls” exhibit unusually large dipole moments, but these are ancillary to the forces that control their intermolecular interactions, as illustrated by two examples involving corannulene. Finally, the assumption that interactions between water and small anions are dictated by the dipole moment of H 2 O is shown to be false in the case of binary halide–water complexes. Furthermore, these examples present a compelling case that electrostatic explanations based on low-order multipole moments are very often counterfactual for nonbonded interactions at close range and should not be taken seriously in the absence of additional justification.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Chiral-Symmetric Higher-Order Topological Phases of Matter

Herein, we introduce novel higher-order topological phases of matter in chiral-symmetric systems (class AIII of the tenfold classification), most of which would be misidentified as trivial by current theories. These phases are protected by “multipole chiral numbers,” bulk integer topological invariants that in 2D and 3D are built from sublattice multipole moment operators, as defined herein. The integer value of a multipole chiral number indicates how many degenerate zero-energy states localize at each corner of a system. These higher-order topological phases of matter are generally boundary-obstructed and robust in the presence of chiral-symmetry-preserving disorder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Delaying the inevitable: tidal disruption in microstate geometries

Microstate geometries in string theory replace the black-hole horizon with a smooth geometric “cap” at the horizon scale. In geometries constructed using superstratum technology, this cap has the somewhat surprising property that induces very large tidal deformations on infalling observers that are far away from it. We find that this large-distance amplification of the tidal effects is also present in horizonless microstate geometries constructed as bubbling solutions, but can be tamed by suitably arranging the bubbles to reduce the strength of some of the gravitational multipole moments. However, despite this taming, these tidal effects still become large at a significant distance from the microstructure. This result suggests that an observer will not fall unharmed into the structure replacing the black hole horizon.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Intermolecular Interactions in Direct Air Capture Materials: Insights from Charge Density Analysis

Direct air capture (DAC) materials enable the removal of CO 2 from the atmosphere, but improving their efficiency requires a detailed understanding of the intermolecular interactions that govern CO 2 sorption and release. Here, we present an experimental electron density study of methylglyoxal-bis(iminoguanidine) (MGBIG), a promising DAC material, using high-resolution X-ray and neutron diffraction data combined with quantum crystallographic analysis. This approach bridges theoretical and experimental data by quantifying electron density distributions and revealing how hydrogen bonds stabilize CO 2 -derived carbonate phases and may influence the desorption behavior. We identify distinct hydrogen-bonding environments in two crystalline carbonate phases: P1, a transient kinetic product, and P3, a thermodynamically stable phase. Multipolar refinement and electrostatic potential and multipole moment calculations precisely map electron density distributions, revealing key hydrogen bonds involved in CO 2 capture. Topological analysis of electron density highlights a cooperative hydrogen-bonding network in the thermodynamically favored P3 phase, where enhanced electron density delocalization and water-mediated interactions contribute to a more stable lattice. Energetic analyses confirm that stronger hydrogen bonding networks enhance the stability of P3 with a binding energy of −607.0 kJ/mol and greater lattice stability (−847.3 kJ/mol) compared to P1 (−302.5 and −571.0 kJ/mol, respectively). Electrostatic potential maps further illustrate polarization patterns that may influence the stability of the binding of CO 2 and release conditions. These findings establish a direct experimental framework for linking electron density distributions to intermolecular interactions in DAC materials, providing a rational design strategy for optimizing sorbents with improved CO 2 capture efficiency and reduced energy demands.

Electron density↗
DOE OSTI2020

The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules

Abstract Maximum diversification of data is a central theme in building generalized and accurate machine learning (ML) models. In chemistry, ML has been used to develop models for predicting molecular properties, for example quantum mechanics (QM) calculated potential energy surfaces and atomic charge models. The ANI-1x and ANI-1ccx ML-based general-purpose potentials for organic molecules were developed through active learning; an automated data diversification process. Here, we describe the ANI-1x and ANI-1ccx data sets. To demonstrate data diversity, we visualize it with a dimensionality reduction scheme, and contrast against existing data sets. The ANI-1x data set contains multiple QM properties from 5 M density functional theory calculations, while the ANI-1ccx data set contains 500 k data points obtained with an accurate CCSD(T)/CBS extrapolation. Approximately 14 million CPU core-hours were expended to generate this data. Multiple QM calculated properties for the chemical elements C, H, N, and O are provided: energies, atomic forces, multipole moments, atomic charges, etc. We provide this data to the community to aid research and development of ML models for chemistry.

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

Quantum mechanical dataset of 836k neutral closed-shell molecules with up to 5 heavy atoms from C, N, O, F, Si, P, S, Cl, Br

Abstract We introduce the Vector-QM24 (VQM24) dataset comprehensively covering all possible neutral closed-shell small organic and inorganic molecules with up to five heavy (p-block) atoms: C, N, O, F, Si, P, S, Cl, Br. All valid stoichiometries, Lewis-rule-consistent graphs, and stable conformers (identified via GFN2-xTB) were enumerated combinatorially, yielding 577k conformational isomers spanning 258k constitutional isomers and 5,599 unique stoichiometries. DFT (ωB97X-D3/cc-pVDZ) optimizations were performed for all, and diffusion quantum Monte Carlo (DMC@PBE0(ccECP/cc-pVQZ)) energies are provided for 10,793 lowest-energy conformers with up to 4 heavy atoms. VQM24 includes structures, vibrational modes, rotational constants, thermodynamic properties (Gibbs free energies, enthalpies, ZPVEs, entropies, heat capacities), and electronic properties such as atomization, electron interaction, exchange-correlation, dispersion energies, multipole moments (dipole to hexadecapole), alchemical potentials, Mulliken charges, and wavefunctions. Machine learning models of atomization energies on this dataset reveal significantly higher complexity than QM9, with none achieving chemical accuracy. VQM24 offers a rigorous, high-fidelity benchmark for evaluating quantum machine learning models.

Science & Technology - Other Topics↗
DOE OSTI2021

Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water

Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. Finally, we propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.

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

Long-range model of vibrational autoionization in core-nonpenetrating Rydberg states of NO

In high orbital angular momentum (ℓ ≥ 3) Rydberg states, the centrifugal barrier hinders the close approach of the Rydberg electron to the ion-core. As a result, these core-nonpenetrating Rydberg states can be well described by a simplified model in which the Rydberg electron is only weakly perturbed by the long-range electric properties (i.e., multipole moments and polarizabilities) of the ion-core. We have used a long-range model to describe the vibrational autoionization dynamics of high-ℓ Rydberg states of nitric oxide (NO). In particular, our model explains the extensive angular momentum exchange between the ion-core and the Rydberg electron that had been previously observed in vibrational autoionization of f (ℓ = 3) Rydberg states. These results shed light on a long-standing mechanistic question around these previous observations and support a direct, vibrational mechanism of autoionization over an indirect, predissociation-mediated mechanism. In addition, our model correctly predicts newly measured total decay rates of g (ℓ = 4) Rydberg states because for ℓ ≥ 4, the non-radiative decay is dominated by autoionization rather than predissociation. We examine the predicted NO + ion rotational state distributions generated by vibrational autoionization of g states and discuss applications of our model to achieve quantum state selection in the production of molecular ions.

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

The completed SDSS-IV extended baryon oscillation spectroscopic survey: geometry and growth from the anisotropic void–galaxy correlation function in the luminous red galaxy sample

ABSTRACT We present an analysis of the anisotropic redshift-space void–galaxy correlation in configuration space using the Sloan Digital Sky Survey extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 luminous red galaxy (LRG) sample. This sample consists of LRGs between redshifts 0.6 and 1.0, combined with the high redshift z > 0.6 tail of the Baryon Oscillation Spectroscopic Survey Data Release 12 CMASS sample. We use a reconstruction method to undo redshift-space distortion (RSD) effects from the galaxy field before applying a watershed void-finding algorithm to remove bias from the void selection. We then perform a joint fit to the multipole moments of the correlation function for the growth rate fσ8 and the geometrical distance ratio DM/DH, finding $f\sigma _8(z_\rm {eff})=0.356\pm 0.079$ and $D_M/D_H(z_\rm {eff})=0.868\pm 0.017$ at the effective redshift $z_\rm {eff}=0.69$ of the sample. The posterior parameter degeneracies are orthogonal to those from galaxy clustering analyses applied to the same data, and the constraint achieved on DM/DH is significantly tighter. In combination with the consensus galaxy BAO and full-shape analyses of the same sample, we obtain fσ8 = 0.447 ± 0.039, DM/rd = 17.48 ± 0.23, and DH/rd = 20.10 ± 0.34. These values are in good agreement with the ΛCDM model predictions and represent reductions in the uncertainties of $13{{\ \rm per\ cent}}$, $23{{\ \rm per\ cent}}$, and $28{{\ \rm per\ cent}}$, respectively, compared to the combined results from galaxy clustering, or an overall reduction of 55 per cent in the allowed volume of parameter space.

Nadathur, Seshadri↗
DOE OSTI2022

Measurements of cosmic expansion and growth rate of structure from voids in the Sloan Digital Sky Survey between redshift 0.07 and 1.0

ABSTRACT We present measurements of the anisotropic cross-correlation of galaxies and cosmic voids in data from the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample, Baryon Oscillation Spectroscopic Survey (BOSS), and extended BOSS luminous red galaxy catalogues from SDSS Data Releases 7, 12, and 16, covering the redshift range 0.07 < z < 1.0. As in our previous work analysing voids in subsets of these data, we use a reconstruction method applied to the galaxy data before void finding in order to remove selection biases when constructing the void samples. We report results of a joint fit to the multipole moments of the measured cross-correlation for the growth rate of structure, fσ8(z), and the ratio DM(z)/DH(z) of the comoving angular diameter distance to the Hubble distance, in six redshift bins. For DM/DH, we are able to achieve a significantly higher precision than that obtained from analyses of the baryon acoustic oscillations and galaxy clustering in the same data sets. Our growth rate measurements are of lower precision but still comparable with galaxy clustering results. For both quantities, the results agree well with the expectations for a Λ cold dark matter model. Assuming a flat Universe, our results correspond to a measurement of the matter density parameter $\Omega _\mathrm{m}=0.337^{+0.026}_{-0.029}$. For more general models, the degeneracy directions obtained are consistent with and complementary to those from other cosmological probes. These results consolidate void–galaxy cross-correlation measurements as a pillar of modern observational cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Extraction of ground-state nuclear deformations from ultrarelativistic heavy-ion collisions: Nuclear structure physics context

The collective-flow-assisted nuclear shape-imaging method in ultrarelativistic heavy-ion collisions (UHICs) has recently been used to characterize nuclear collective states. In this paper, we assess the foundations of the shape-imaging technique employed in these studies. We argue that some current UHIC nuclear imaging techniques neglect fundamental aspects of spontaneous symmetry breaking and symmetry restoration in colliding ions and incorrectly infer one-body multipole moments from studies of nucleonic correlations. Therefore, the impact of this approach on nuclear structure research has been overstated. Conversely, efforts to incorporate existing knowledge on nuclear shapes into analysis pipelines can be beneficial for benchmarking tools and calibrating models used to extract information from ultrarelativistic heavy-ion experiments.

Nuclear data analysis & compilation↗
DOE OSTI2022

Bulk photospin effect: Calculation of electric spin susceptibility to second order in an electric field

Here, we compute the electric spin susceptibility of Bloch electrons with spin-orbit coupling to second order. We find that it is possible to generate a nonequilibrium spin polarization in the bulk of nonmagnetic inversion-symmetric materials using linearly polarized electric fields, but the process depends on interband coherence and produces heating. It may be possible to avoid heating with circular polarization in certain scenarios. The standard Edelstein effect and spin orientation effects are recovered in appropriate limits within the formalism. Finally, the electric spin susceptibility of metals has contributions proportional to spin multipole moments of the Fermi sea that dominate the low frequency spin response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Multimodal fission from self-consistent calculations

When multiple fission modes coexist in a given nucleus, distinct fragment yield distributions appear. Multimodal fission has been observed in a number of fissioning nuclei spanning the nuclear chart, and this phenomenon is expected to affect the nuclear abundances synthesized during the rapid neutron-capture process (𝑟-process). In this study, we generalize the previously proposed hybrid model for fission-fragment yield distributions to predict competing fission modes and estimate the resulting yield distributions. Here, our framework allows for a comprehensive large-scale calculation of fission-fragment yields suited for 𝑟-process nuclear network studies. Nuclear density functional theory is employed to obtain the potential energy and collective inertia tensor on a multidimensional collective space defined by mass multipole moments. Fission pathways and their relative probabilities are determined using the nudged elastic band method. Based on this information, mass and charge fission yields are predicted using the recently developed hybrid model. Fission properties of fermium isotopes are calculated in the axial quadrupole-octupole collective space for three energy density functionals (EDFs). Disagreement between the EDFs appears when multiple fission modes are present. Within our framework, the UNEDF⁢1 HFB EDF agrees best with experimental data. Calculations in the axial quadrupole-octupole-hexadecapole collective space improve the agreement with the experiment for SkM*. We also discuss the sensitivity of fission predictions on the choice of EDF for several superheavy nuclei. Fission-fragment yield predictions for nuclei with multiple fission modes are sensitive to the underlying EDF. For large-scale calculations in which a minimal number of collective coordinates is considered, UNEDF⁢1 HFB provides the best description of experimental data, though the sensitivity motivates robust quantification of the uncertainties of the theoretical model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Tidal Love numbers of analog black holes

Tidal Love numbers quantify the conservative static response of compact objects to external tidal fields, and are found to vanish exactly for asymptotically flat black holes in four-dimensional general relativity. Many aspects of the physics of black holes have an analog in the theory of supersonic acoustic flows, including the existence of an event horizon and associated phenomena, such as quasinormal modes and superradiance. In this paper, we investigate the tidal Love numbers of acoustic black holes in different number of dimensions. We find that they exhibit a number of similar properties as higher-dimensional general relativistic black holes, such as logarithmic running with radial distance, and vanishing tidal response for special multipole moments. As a result, we show that the latter is a consequence of ladder symmetries, analogous to those identified for black holes.

Classical black holes↗
DOE OSTI2024

Charge Conservation beyond Uniformity: Spatially Inhomogeneous Electromagnetic Response in Periodic Solids

Nonlinear electromagnetic response functions have reemerged as a crucial tool for studying quantum materials, due to recently appreciated connections between optical response functions, quantum geometry, and band topology. Most attention has been paid to responses to spatially uniform electric fields, relevant to low-energy optical experiments in conventional solid state materials. However, magnetic and magnetoelectric phenomena are naturally connected by responses to spatially varying electric fields due to Maxwell’s equations. Furthermore, in the emerging field of moiré materials, characteristic lattice scales are much longer, allowing spatial variation of optical electric fields to potentially have a measurable effect in experiments. In order to address these issues, we develop a formalism for computing linear and nonlinear responses to spatially inhomogeneous electromagnetic fields. Starting with the continuity equation, we derive an expression for the second-quantized current operator that is manifestly conserved and model independent. Crucially, our formalism makes no assumptions on the form of the microscopic Hamiltonian and so is applicable to model Hamiltonians derived from tight-binding or calculations. We then develop a diagrammatic Kubo formalism for computing the wave vector dependence of linear and nonlinear conductivities, using Ward identities to fix the value of the diamagnetic current order by order in the vector potential. We apply our formula to compute the magnitude of the Kerr effect at oblique incidence for a model of a moiré-Chern insulator and demonstrate the experimental relevance of spatially inhomogeneous fields in these systems. We further show how our formalism allows us to compute the (orbital) magnetic multipole moments and magnetic susceptibilities in insulators. Turning to nonlinear response, we use our formalism to compute the second-order transverse response to spatially varying transverse electric fields in our moiré-Chern insulator model, with an eye toward the next generation of experiments in these systems. Published by the American Physical Society 2024

Physics↗
DOE OSTI2024

Rapid Tuning of Synchrotron Surrogate Model at the Recycler Ring

The 8 GeV proton-storage Recycler Ring (RR) is essential for reaching megawatt beam intensity goals for the DUNE neutrino beam at Fermilab. Custom shims on each RR permanent magnet were designed to cancel manufacturing defects and bring magnetic fields to the design values. Remaining imperfections cause the observed tune variation vs energy to deviate from what is calculated using the design fields. Using the POUNDERS (“Practical Optimization Using No Derivatives for sums of Squares”) optimization method with Synergia in the loop, we demonstrate rapid convergence to a set of additive, higher-order multipole moments of these magnetic shims which reproduce that observed variation, and show that the convergence advantage grows with the parameter-space dimensionality.

43 PARTICLE ACCELERATORS↗