Search NASA⌕ Search

SEARCH · Search NASA

Results for “MULTIPOLES”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Non-perturbative Approach to Computing Seismic Normal Modes in Rotating Planets

In this work, a continuous Galerkin method based approach is presented to compute the seismic normal modes of rotating planets. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core using a polynomial filtering eigensolver. The relevant elastic-gravitational system of equations, including the Coriolis force, is subjected to a mixed finite-element method, while self-gravitation is accounted for with the fast multipole method. Our discretization utilizes fully unstructured tetrahedral meshes for both solid and fluid regions. The relevant eigenvalue problem is solved by a combination of several highly parallel and computationally efficient methods. We validate our three-dimensional results in the non-rotating case using analytical results for constant elastic balls, as well as numerical results for an isotropic Earth model from standard “radial” algorithms. We also validate the computations in the rotating case, but only in the slowly-rotating regime where perturbation theory applies, because no other independent algorithms are available in the general case. The algorithm and code are used to compute the point spectra of eigenfrequencies in several Earth and Mars models studying the effects of heterogeneity on a large range of scales.

58 GEOSCIENCES↗
DOE OSTI2019

Magnetic and Thermal Design of HTS Quadrupole Magnet for Newly Developed Superconducting Proton Cyclotron Beam Line

A magnetic and thermal design of high-temperature superconducting (HTS) quadrupole magnet for newly developed superconducting proton cyclotron was presented in this research. With superconducting technology, the design can reduce the magnet size and remove the heat loads more efficiently. Calculations are conducted with finite element method (FEM) to study the quadrupole magnet’s magnetic and thermal properties. For the magnet cold-mass and cryostat system, the heat load is mainly generated by conduction and radiation. Multilayer thermal insulation and G10 supports are used to restrict them. The magnetic field distribution in the HTS coils is calculated to consider the critical current degradation of the YBCO tapes. The effects on the field gradient quality in relation to the pole tip profile and end chamfer are analyzed. The sixth multipole component in the integral field can be improved with end chamfer. Finally, a 20 T/m HTS quadrupole magnet with integral field uniformity less than 0.05% out to 75% of the inscribed radius is proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Front-End Electron Transfer Dissociation Coupled to a 21 Tesla FT-ICR Mass Spectrometer for Intact Protein Sequence Analysis

High resolution mass spectrometry is a key technology for in-depth protein characterization. High-field Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) enables high-level interrogation of intact proteins in the most detail to date. However, an appropriate complement of fragmentation technologies must be paired with FTMS to provide comprehensive sequence coverage, as well as characterization of sequence variants, and post-translational modifications. Here we describe the integration of front-end electron transfer dissociation (FETD) with a custom-built 21 tesla FT-ICR mass spectrometer, which yields unprecedented sequence coverage for proteins ranging from 2.8 to 29 kDa, without the need for extensive spectral averaging (e.g., ~60% sequence coverage for apo-myoglobin with four averaged acquisitions). The system is equipped with a multipole storage device separate from the ETD reaction device, which allows accumulation of multiple ETD fragment ion fills. Consequently, an optimally large product ion population is accumulated prior to transfer to the ICR cell for mass analysis, which improves mass spectral signal-to-noise ratio, dynamic range, and scan rate. We find a linear relationship between protein molecular weight and minimum number of ETD reaction fills to achieve optimum sequence coverage, thereby enabling more efficient use of instrument data acquisition time. Finally, real-time scaling of the number of ETD reactions fills during method-based acquisition is shown, and the implications for LC-MS/MS top-down analysis are discussed. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Mobilities of polydisperse hard spheres near a no-slip wall

In the article, we have derived an analytical formulation for far-field hydrodynamic interactions among unequal size hard spheres near a no-slip wall and have implemented the formulation in a Stokesian dynamics model to simulate a suspension of polydisperse particles in a semi-bounded domain. The formulation is based on the multipole expansion of the boundary integral of Stokes flow, and the mobility tensors are deduced from Fáxen’s law together with Green’s function for Stokes flow near a no-slip wall. Lubrication approximation is incorporated to account for the close-distance interactions between any two particles and also between the particles and the wall. The implementation is validated against previous formulations for equal size particles and against a boundary-element code for unequal size particles. The code can be used to simulate the Stokesian or Brownian interaction of unequal particles in presence of a no-slip wall. In the current study, we applied the Stokesian dynamics model to investigate the trajectories of two unequal hard spheres and their redistribution during sedimentation parallel to a wall. We further used it to demonstrate the cases of many-particle sedimentation toward or parallel to a wall.

catalysis (heterogeneous)↗
DOE OSTI2018

The fast linear accelerator modeling engine for FRIB online model service

Commissioning of a large accelerator facility like FRIB needs support from an online beam dynamics model. Considering the new physics challenges of FRIB such as modeling of non-axisymmetric superconducting RF cavities and multi-charge state acceleration, there is no readily available online beam tuning code. The design code of FRIB super-conducting linac, IMPACT-Z, is not suitable for online tuning because of its code design and running speed. Therefore, the Fast Linear Accelerator Modeling Engine (FLAME), specifically designed to fulfill FRIB’s online modeling challenges, is proposed. The physics model of FLAME, especially its novel way of modeling non-axisymmetric superconducting RF cavities using a multipole expansion based thin-lens kick model, is discussed in detail, and the benchmark results against FRIB design code is presented, after which the software design strategy of FLAME and its execution speed is presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Reformulation of plasma spectral line broadening with computer simulations

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating radiator–plasma interaction obtained from computer simulations. These calculations usually described the radiator by a localized charge distribution of negligible size and the plasma by an electric potential plus inhomogeneities at the radiator nucleus. Using an alternative approach, it is shown that starting from the full Coulomb interaction and then introducing approximations leads to expressions that naturally take advantage of readily available data from atomic packages. Applying this alternative approach to screened Coulomb interactions in so-called trivial molecular dynamics to account for neglected plasma correlations leads to a modified multipole expansion compared to the standard method. Issues associated with the ad hoc introduction of the screened Coulomb potential are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The j and k dependencies of the line coupling and line mixing effects: Theoretical studies of the relaxation matrices of N 2 -broadened CH 3 D

Line coupling and line mixing effects in parallel and perpendicular bands of CH 3 D perturbed by N 2 have been studied. The work focuses on exhibiting the j and k dependencies of these two processes. The calculations were based on a previously reported anisotropic intermolecular potential including both the long-range multipole, induction, and dispersion forces and a short-range atom-atom model. Throughout this work, it is shown that components with L 1 = 3 of the atom-atom model are dominant contributions to the diffusion operator. As a consequence, in comparison with other molecular systems such as the CH 3 Cl-N 2 and CH 3 I-N 2 , theoretically predicted line coupling and line mixing effects exhibit completely new j and k dependencies. In general, the theoretically calculated halfwidths and intra-doublets’ off-diagonal elements of the relaxation matrix are in reasonable agreement with measurements.

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

Simulation of beam–beam interaction with crab cavities for $\mathrm{LHC}$ upgrade

The crab cavities are a critical component in the high luminosity LHC upgrade project to compensate the luminosity loss from large crossing angle collision. However, these crab cavities will not be perfect in the real accelerator. In this paper, we studied the effects of crab cavity imperfections on colliding beam luminosity lifetime degradation for the LHC upgrade through detailed numerical simulations. Our simulation results suggest that the white noise jitter in the crab cavity RF phase and voltage rms amplitudes should be kept below a few 10 -5 for a good luminosity lifetime, while with frequency-dependent jitter, the amplitudes should be kept below a few 10 -4 for a good lifetime. As a result, the RF multipole errors in the current crab cavity designs are small enough and would not cause extra luminosity degradation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Investigation of the isoscalar monopole response in the proton-rich nucleus 14 O

Deuteron inelastic scattering on 14 O was measured in inverse kinematics using an active-target time projection chamber and a magnetic spectrograph. The experimental technique enabled precise measurements of deuteron recoiling particles in coincidence with beam-like fragments detected in the spectrograph focal plane. The double differential cross section was reconstructed for scattering angles of 3–6 degrees and excitation energies up to 26 MeV. The monopole strength distribution was obtained from the data using a multipole decomposition analysis. The results were compared to quasiparticle random-phase approximation (QRPA) and generator coordinate method (GCM) calculations. The QRPA calculation accurately describes experimental data in the energy range of 13 to 26 MeV. GCM calculations assuming a 12 C (g . s .) + p + p cluster configuration were used to determine the 0 + strength in 14 O below 13 MeV. The monopole transition strength of these cluster states provides a good description of the experimental distribution in the 9–11 MeV region, while the 0$^{+}_{2}$ state accounts for only a small fraction of the experimental strength around 6 MeV.

Active target↗
DOE OSTI2023

Guiding centre motion for particles in a ponderomotive magnetostatic end plug

The Hamiltonian dynamics of a single particle in a rotating plasma column, interacting with an magnetic multipole is perturbatively solved for up to second order, using the method of Lie transformations. First, the exact Hamiltonian is expressed in terms of canonical action-angle variables, and then an approximate integrable Hamiltonian is introduced, using another set of actions and angles, which describe the centre of oscillation for the particle. The perturbation introduces an effective ponderomotive potential, which to leading order is positive. At the second order, the pseudopotential consists of a sum of terms of the Miller form, and can have either sign. Additionally, at second order, the ponderomotive interaction introduces a modification to the particle effective mass, when considering the motion along the column axis. It is found that particles can be axially confined by the ponderomotive potentials, but acquire radial excursions which scale as the confining potential. The radial excursions of the particle along its trajectory are investigated, and a condition for the minimal rotation frequency for which the particle remains radially confined is derived. Last, we comment on the changes to the aforementioned solution to the pseudopotentials and particle trajectory in the case of resonant motion, that is, a motion which has the same periodicity as the perturbation.

Rubin, T. (ORCID:0000000164850096)↗
DOE OSTI2019

Dispersion Forces: Neither Fluctuating Nor Dispersing

Dispersion forces are ubiquitous in chemistry, yet they are often misunderstood. This article provides background into why they are called dispersion forces and explains how to describe them in terms of time-independent quantum mechanics. Here, the article also describes the breakdown of the multipole series that is often used to describe dispersion forces, and it comments on the representation of dispersion forces in molecular mechanics. In general, the article tries to provide a proper framework for understanding and describing dispersion forces.

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

Polarizable Water Potential Derived from a Model Electron Density

A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.

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

pyEF: A Python Framework for QM and QM/MM Atom-Wise Electric Field Analysis

We introduce pyEF, a software package for computing molecular electric fields, electrostatic interaction energies, and electrostatic potentials from quantum mechanical (QM) atom-centered multipole expansions with atom-wise decomposable contributions. We demonstrate the computational efficiency and accuracy of this QM-derived electric field evaluation tool through several tests. To assess the influence of the underlying QM method and charge partitioning scheme on these electrostatic quantities, we analyze over 250 configurations of an acetone solute molecule in five solvents of variable polarity. We find that electric field calculations are highly sensitive to the choice of charge partitioning method. Even among real-space charge schemes, acetone Stark tuning rates differ by up to a factor of 2. Benchmarking computed solvent dipole moments against experimental bulk values, we conclude that the CM5, ADCH, and Hirshfeld-I charge schemes most reliably capture solvent electrostatics and therefore provide a more faithful foundation for computing electric fields. When constructed from these real-space charges, electric fields are nearly insensitive to basis set size and monotonically increase in magnitude with higher Fock exchange. We also demonstrate efficient convergence of QM electrostatics when more distant molecules are represented solely by MM point charges, reducing computational overhead. Leveraging these findings, we demonstrate the use of pyEF to deduce environmental effects on a transition metal complex from a Ga 4 L 6 12– nanocage and quantify the dominant role of organic linkers in orchestrating electrostatic preorganization.

electric fields↗
DOE OSTI2019

Probing Molecular Chirality by Orbital-Angular-Momentum-Carrying X-ray Pulses

A twisted X-ray beam with orbital angular momentum is employed in a theoretical study to probe molecular chirality. A nonlocal response description of the matter-field coupling is adopted to account for the field short wavelength and the structured spatial profile. In this study, we use the minimal-coupling Hamiltonian, which implicitly takes into account the multipole contributions to all orders. The combined interactions of the spin and orbital angular momentum of the X-ray beam give rise to circular-helical dichroism signals, which are stronger than ordinary circular dichroism signals, and may serve as a useful tool for the study of molecular chirality in the X-ray regime.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Second Linear Response Theory and the Analytic Calculation of Excited-State Properties

We present a method based on second linear response time-dependent density func- tional theory (TDDFT) to calculate permanent and transition multipoles of excited states, which are required to compute excited-state absorption/emission spectra, multi- photon optical processes, among others. In previous work, we examined computations based on second linear response theory in which linear response TDDFT was employed twice. In contrast, the present methodology only requires information from a single linear response calculation to compute the excited-state properties. These are eval- uated analytically through various algebraic operations involving electron repulsion integrals and excitation vectors. The present derivation focuses on full many-body wave functions instead of single orbitals, as in our previous approach. Here, we test the proposed method by applying it to several diatomic and triatomic molecules. This shows that the computed excited-state dipoles are consistent with respect to reference equation-of-motion coupled cluster calculations.

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

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

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

Understanding Chemical Bonding in Alloys and the Representation in Atomistic Simulations

Alloys are widely used in catalysts and structural materials. The nature of chemical bonding and the origin of alloy formation energies, defect energies, and interfacial properties have not been well understood to date but are critical to material performance. In this contribution, we explain the polar nature of chemical bonding and an implementation in classical and reactive atomistic simulations to understand such properties more quantitatively. Electronegativity differences between metal atoms lead to polar bonding, and exothermic alloy formation energies are related to charge transfer between the different elements. These differences can be quantified by atomic charges using pairwise charge increments, determined by matching the computed alloy formation energy to experimentally measured alloy formation energies using pair potentials for the pure metals. The polar character of alloys is comparable to organic molecules and partially ionic minerals, for example, AlNi and AlNi3 alloys assume significant atomic charges of ±0.40e and +0.60e/–0.20e, respectively. The subsequent analysis of defect sites and defect energies using force-field-based calculations shows excellent agreement with calculations using density functional theory and embedded atom models (EAM). The formation of vacancy and antisite defects is characterized by a redistribution of charge in the first shell of neighbor atoms in the classical models whereby electroneutrality is maintained and charge increments correlate with differences in electronegativity. The proposed atomic charges represent internal dipole and multipole moments, consistent with existing definitions for organic and inorganic compounds and with the extended Born model (Heinz, H.; Suter, U. W. J. Phys. Chem. B2004,108 (47), 18341–18352). The method can be applied to any alloy and has a reproducibility of ±10%. In contrast, quantum mechanical charge schemes remain associated with deviations exceeding ±100%. The atomic charges for alloys provide a simple initial measure for the internal electronic structure, surface adsorption of molecules, and reactivity in catalysis and corrosion. The models are compatible with the Interface force field (IFF), CHARMM, AMBER, OPLS-AA, PCFF, CVFF, and GROMOS for reliable atomistic simulations of alloys and their interfaces with minerals and electrolytes from the nanometer scale to the micrometer scale.

Chemistry↗
DOE OSTI2025

Completely Multipolar Model for Many-Body Water–Ion and Ion–Ion Interactions

This work constructs an advanced force field, the Completely Multipolar Model (CMM), to quantitatively reproduce each term of an energy decomposition analysis (EDA) for aqueous solvated alkali metal cations and halide anions and their ion pairings. We find that all individual EDA terms remain well-approximated in the CMM for ion-water and ion-ion interactions, except for polarization, which shows errors due to the partial covalency of ion interactions near their equilibrium. We quantify the onset of the dative bonding regime by examining the change in molecular polarizability and Mayer bond indices as a function of distance, showing that partial covalency manifests by breaking the symmetry of atomic polarizabilities while strongly damping them at short-range. This motivates an environment-dependent atomic polarizability parameter that depends on the strength of the local electric field experienced by the ions to account for strong damping, with anisotropy introduced by atomic multipoles. The resulting CMM model for ions provides accurate dimer surfaces and three-body polarization and charge transfer compared to EDA, and shows excellent performance on various ion benchmarks including vibrational frequencies and cluster geometries.

Heindel, Joseph P↗