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
NASA NTRS1993

Point charge representation of multicenter multipole moments in calculation of electrostatic properties

Distributed Point Charge Models (PCM) for CO, (H2O)2, and HS-SH molecules have been computed from analytical expressions using multi-center multipole moments. The point charges (set of charges including both atomic and non-atomic positions) exactly reproduce both molecular and segmental multipole moments, thus constituting an accurate representation of the local anisotropy of electrostatic properties. In contrast to other known point charge models, PCM can be used to calculate not only intermolecular, but also intramolecular interactions. Comparison of these results with more accurate calculations demonstrated that PCM can correctly represent both weak and strong (intramolecular) interactions, thus indicating the merit of extending PCM to obtain improved potentials for molecular mechanics and molecular dynamics computational methods.

Non-NASA Center↗
NASA NTRS1992

Cumulative atomic multipole moments complement any atomic charge model to obtain more accurate electrostatic properties

The quality of several atomic charge models based on different definitions has been analyzed using cumulative atomic multipole moments (CAMM). This formalism can generate higher atomic moments starting from any atomic charges, while preserving the corresponding molecular moments. The atomic charge contribution to the higher molecular moments, as well as to the electrostatic potentials, has been examined for CO and HCN molecules at several different levels of theory. The results clearly show that the electrostatic potential obtained from CAMM expansion is convergent up to R-5 term for all atomic charge models used. This illustrates that higher atomic moments can be used to supplement any atomic charge model to obtain more accurate description of electrostatic properties.

NASA Discipline Number 52-20↗
NASA NTRS1986

Exact formulas for multipole moments using Slater-type molecular orbitals

A triple infinite sum of formulas expressed as an expansion in Legendre polynomials is generated by use of computer algebra to represent the potential from the midpoint of two Slater-type orbitals; the charge density that determines the potential is given as the product of the two orbitals. An example using 1s orbitals shows that only a few terms are needed to obtain four-figure accuracy. Exact formulas are obtained for multipole moments by means of a careful study of expanded formulas, allowing an 'extrapolation to infinity'. This Loewdin alpha-function approach augmented by using a C matrix to characterize Slater-type orbitals can be readily generalized to all cases.

Jones, H. W.↗
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↗
NASA NTRS1994

Multipole moments of stellar oscillation modes

The oscillating mass 2(exp l)-pole moment, M(sub nl), of a star in a given (normalized) oscillation mode determines the energy that can be absorbed by the mode in a tidal interaction and the power radiated by the mode in gravitational waves, both of which are proportional to (absolute value of M(sub nl))(exp 2). The coefficients in the expansion of the vector fields del(r(exp l)Y(sub lm)(theta, phi)) in terms of the displacement fields of modes of given l and m are proportional to M(sub nl). This expansion leads to a sum rule sum over n(absolute value of M(sub nl))(exp 2) = constant. For stars of weak to moderate central condensation (such as neutron stars), the f-mode is well approximated by the vector field being expanded, and therefore it takes the lion's share of the sum. Thus the multipole moments of all other modes must be small. In there numerical evaluation, it is necessary to know the shape of the eigenfunctions quite precisly, since a small f-mode contamination can significantly increase the obtained values. This contamination occurs in some `hybrid' numerical computations of neutron star oscillations with relativistic equilibrium stars and Newtonian dynamics (e.g., McDermott et al. 1988). In this case, it is due to a slight inconsistency in the models and leads to a large overestimate of the power radiated in gravitational waves by modes other than the f-mode, although their oscillation periods are nearly unaffected.

Reisenegger, Andreas↗
NASA NTRS2007

Anatomy of the Binary Black Hole Recoil: A Multipolar Analysis

We present a multipolar analysis of the recoil velocity computed in recent numerical simulations of binary black hole coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including 1 = 4 are sufficient to accurately reproduce the final recoil velocity (= 98%) and that only a few dominant modes contribute significantly to it (2 95%). We describe how the relative amplitude, and more importantly, the relative phase, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ring-down phases. We also find that the numerical results can be reproduced, to a high level of accuracy, by an effective Newtonian formula for the multipole moments obtained by replacing in the Newtonian formula the radial separation with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes. Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of '.anti-kick" for each of the simulations. Lastly, we apply this multipolar analysis to understand the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes.

Schnittman, Jeremy↗
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↗
NASA NTRS1991

Multipolar magnetic fields in neutron stars

The small periods and period derivatives of millisecond pulsars permit interesting upper bounds to be placed on the strength of multipole components in their magnetic fields if the loss of rotational energy may be approximated by vacuum electromagnetic radiation formulas. The surface field strength of millisecond pulsars may be as high as about 10 to the 12th G if the dominant multipole component is l not less than 3. Significant multipole moments may help explain the unusual pulse morphology of millisecond pulsars. If the magnetic fields in gamma-ray burst sources (old neutron stars?) have the same evolutionary history as in millisecond pulsars, the about 10 to the 12th G local fields observed in several gamma-ray bursts must be due to multipole moments with l not less than 3.

Krolik, Julian H.↗
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↗
NASA NTRS1993

Multipole correction of atomic monopole models of molecular charge distribution. I. Peptides

The defects in atomic monopole models of molecular charge distribution have been analyzed for several model-blocked peptides and compared with accurate quantum chemical values. The results indicate that the angular characteristics of the molecular electrostatic potential around functional groups capable of forming hydrogen bonds can be considerably distorted within various models relying upon isotropic atomic charges only. It is shown that these defects can be corrected by augmenting the atomic point charge models by cumulative atomic multipole moments (CAMMs). Alternatively, sets of off-center atomic point charges could be automatically derived from respective multipoles, providing approximately equivalent corrections. For the first time, correlated atomic multipoles have been calculated for N-acetyl, N'-methylamide-blocked derivatives of glycine, alanine, cysteine, threonine, leucine, lysine, and serine using the MP2 method. The role of the correlation effects in the peptide molecular charge distribution are discussed.

NASA Discipline Exobiology↗
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↗