Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Real-space quadrature: A convenient, efficient representation for multipole expansions

Multipoles are central to the theory and modeling of polarizable and nonpolarizable molecular electrostatics. This has made a representation in terms of point charges a highly sought after goal, since rotation of multipoles is a bottleneck in molecular dynamics implementations. All known point charge representations are orders of magnitude less efficient than spherical harmonics due to either using too many fixed charge locations or due to nonlinear fitting of fewer charge locations. We present the first complete solution to this problem—completely replacing spherical harmonic basis functions by a dramatically simpler set of weights associated to fixed, discrete points on a sphere. This representation is shown to be space optimal. It reduces the spherical harmonic decomposition of Poisson’s operator to pairwise summations over the point set. As a corollary, we also shows exact quadrature-based formulas for contraction over trace-free supersymmetric 3D tensors. Moreover, multiplication of spherical harmonic basis functions translates to a direct product in this representation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

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

Spin induced multipole moments for the gravitational wave flux from binary inspirals to third Post-Newtonian order

Using effective field theory techniques we calculate the source multipole moments needed to obtain the spin contributions to the power radiated in gravitational waves from inspiralling compact binaries to third Post-Newtonian order (3PN). The multipoles depend linearly and quadratically on the spins and include both spin(1)spin(2) and spin(1)spin(1) components. The results in this paper provide the last missing ingredient required to determine the phase evolution to 3PN including all spin effects which we will report in a separate paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Spontaneous breaking of multipole symmetries

Multipole symmetries are of interest both as a window on fracton physics and as a crucial ingredient in realizing new universality classes for quantum dynamics. Here we address the question of whether and when multipole symmetries can be spontaneously broken, both in thermal equilibrium and at zero temperature. In this work, we derive generalized Mermin-Wagner arguments for the total or partial breaking of multipolar symmetry groups and generalized Imry-Ma arguments for the robustness of such multipolar symmetry breaking to disorder. We present both general results and explicit examples. Our results should be directly applicable to quantum dynamics with multipolar symmetries and also provide a useful stepping stone to understanding the robustness of fracton phases to thermal fluctuations, quantum fluctuations, and disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Mathematical and Numerical Aspects of the Adaptive Fast Multipole Poisson-Boltzmann Solver

This paper summarizes the mathematical and numerical theories and computational elements of the adaptive fast multipole Poisson-Boltzmann (AFMPB) solver. We introduce and discuss the following components in order: the Poisson-Boltzmann model, boundary integral equation reformulation, surface mesh generation, the nodepatch discretization approach, Krylov iterative methods, the new version of fast multipole methods (FMMs), and a dynamic prioritization technique for scheduling parallel operations. For each component, we also remark on feasible approaches for further improvements in efficiency, accuracy and applicability of the AFMPB solver to large-scale long-time molecular dynamics simulations. Lastly, the potential of the solver is demonstrated with preliminary numerical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Rogue Mode Shileding in NSLS-II Multipole Vacuum Chambers

Modes with transverse electric field (TE-modes) in the NSLS-II multipole vacuum chamber can be generated at frequencies above 450MHz due to its geometric dimensions. Since the NSLS-II BPM system monitors signals within 10 MHz band at RF frequency of 500 MHz, frequencies of higher-order modes (HOM) can be generated within the transmission band of the band pass filter. In order to avoid systematic errors in the NSLS-II BPM system, we introduced frequency shift of HOMs by using RF metal shielding located in the antechamber slot. We demonstrated numerical modeling and experimental studies of the spurious TE modes in the NSLS-II vacuum chambers with antechamber slot. Calculated frequencies of TE-modes in considered chambers with and without RF shielding were verified experimentally. Flexible BeCu RF shielding inside each chamber at proper location shifts frequencies of H{sub 10p}-modes above {approx}900MHz, except chambers S6 odd and even. These chambers need special attention because of synchrotron radiation from downstream magnets. S6 odd multipole vacuum chamber needs to be measured and the RF shielding length has to be optimized. RF shielding looks adequate for baseline design. Fifty percent of open space provides adequate pumping speed.

25mm434mhz↗
DOE OSTI2015

Effects of Crab Cavities' Multipole Content in an Electron-Ion Collider

The impact on the beam dynamics of the Medium Energy Electron-Ion Colider (MEIC) due to the multipole content of the 750 MHz crab cavity was studied using thin multipole elements for 6D phase space particle tracking in ELEGANT. Target values of the sextupole component for the cavity’s field expansion were used to perform preliminary studies on the proton beam stability when compared to the case of pure dipole content of the rf kicks. Finally, important effects on the beam sizes due to non-linear components of the crab cavities’ fields were identified, and some criteria for their future study were proposed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

A multipole-expanded effective field theory for vortex ring-sound interactions

The low-energy dynamics of a zero temperature superfluid or of the compressional modes of an ordinary fluid can be described by a simple effective theory for a scalar field — the superfluid ‘phase’. However, when vortex lines are present, to describe all interactions in a local fashion one has to switch to a magnetic-type dual two-form description, which comes with six degrees of freedom (in place of one) and an associated gauge redundancy, and is thus considerably more complicated. Here we show that, in the case of vortex rings and for bulk modes that are much longer than the typical ring size, one can perform a systematic multipole expansion of the effective action and recast it into the simpler scalar field language. In a sense, in the presence of vortex rings the non-single valuedness of the scalar can be hidden inside the rings, and thus out of the reach of the multipole expansion. As an application of our techniques, we compute by standard effective field theory methods the sound emitted by an oscillating vortex ring.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Spinor representations for fields with any spin: Lorentz tensor basis for operators and covariant multipole decomposition

This paper discusses a framework to parametrize and decompose operator matrix elements for particles with higher spin (j > 1/2) using chiral representations of the Lorentz group, i.e. the (j, 0) and (0, j) representations and their parity-invariant direct sum. Unlike traditional approaches that require imposing constraints to eliminate spurious degrees of freedom, these chiral representations contain exactly the 2j + 1 components needed to describe a spin-j particle. The central objects in the construction are the t-tensors, which are generalizations of the Pauli four-vector σ μ for higher spin. For the generalized spinors of these representations, we demonstrate how the algebra of the t-tensors allows to formulate a generalization of the Dirac matrix basis for any spin. For on-shell bilinears, we show that a set consisting exclusively of covariant multipoles of order 0 ≤ m ≤ 2j forms a complete basis. We provide explicit expressions for all bilinears of the generalized Dirac matrix basis, which are valid for any spin value. As a byproduct of our derivations we present an efficient algorithm to compute the t-tensor matrix elements. The formalism presented here paves the way to use a more unified approach to analyze the non-perturbative QCD structure of hadrons and nuclei across different spin values, with clear physical interpretation of the resulting distributions as covariant multipoles.

Angular momentum of light↗
DOE OSTI2016

Static electric and magnetic multipole susceptibilities for Dirac one-electron atoms in the ground state

We present tabulated data for several families of static electric and magnetic multipole susceptibilities for hydrogenic atoms with nuclear charge numbers from the range 1⩽Z⩽137. Atomic nuclei are assumed to be point-like and spinless. The susceptibilities considered include the multipole electric polarizabilities α{sub EL→EL} and magnetizabilities (magnetic susceptibilities) χ{sub ML→ML} with 1⩽L⩽4 (i.e., the dipole, quadrupole, octupole and hexadecapole ones), the electric-to-magnetic cross-susceptibilities α{sub EL→M(L−1)} with 2⩽L⩽5 and α{sub EL→M(L+1)} with 1⩽L⩽4, the magnetic-to-electric cross-susceptibilities χ{sub ML→E(L−1)} with 2⩽L⩽5 and χ{sub ML→E(L+1)} with 1⩽L⩽4 (it holds that χ{sub ML→E(L∓1)}=α{sub E(L∓1)→ML}), and the electric-to-toroidal-magnetic cross-susceptibilities α{sub EL→TL} with 1⩽L⩽4. Numerical values are computed from general exact analytical formulas, derived by us elsewhere within the framework of the Dirac relativistic quantum mechanics, and involving generalized hypergeometric functions {sub 3}F{sub 2} of the unit argument.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

An improved fast multipole method for electrostatic potential calculations in a class of coarse-grained molecular simulations

This paper presents a novel algorithm to approximate the long-range electrostatic potential field in the Cartesian coordinates applicable to 3D coarse-grained simulations of biopolymers. In such models, coarse-grained clusters are formed via treating groups of atoms as rigid and/or flexible bodies connected together via kinematic joints. Therefore, multibody dynamic techniques are used to form and solve the equations of motion of such coarse-grained systems. In this article, the approximations for the potential fields due to the interaction between a highly negatively/positively charged pseudo-atom and charged particles, as well as the interaction between clusters of charged particles, are presented. These approximations are expressed in terms of physical and geometrical properties of the bodies such as the entire charge, the location of the center of charge, and the pseudo-inertia tensor about the center of charge of the clusters. Further, a novel substructuring scheme is introduced to implement the presented far-field potential evaluations in a binary tree framework as opposed to the existing quadtree and octree strategies of implementing fast multipole method. Using the presented Lagrangian grids, the electrostatic potential is recursively calculated via sweeping two passes: assembly and disassembly. In the assembly pass, adjacent charged bodies are combined together to form new clusters. Then, the potential field of each cluster due to its interaction with faraway resulting clusters is recursively calculated in the disassembly pass. The method is highly compatible with multibody dynamic schemes to model coarse-grained biopolymers. Since the proposed method takes advantage of constant physical and geometrical properties of rigid clusters, improvement in the overall computational cost is observed comparing to the tradition application of fast multipole method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Doping Bottleneck in Hematite: Multipole Clustering by Small Polarons

We report highly effective doping in transition metal oxides is critical to fundamentally overcome low carrier conductivity due to small polaron formation and reach their ideal efficiency for energy conversion applications. However, the optimal doping concentration in polaronic oxides such as hematite has been extremely low, for example, less than a percent, which hinders the benefits of doping for practical applications. In this work, we investigate the underlying mechanism of low optimal doping concentration with group IV (Ti, Zr, and Hf) and XIV (Si, Ge, Sn, and Pb) dopants from first-principles calculations. We find that novel dopant-polaron clustering occurs even at very low dopant concentrations and resembles electric multipoles. These multipoles can be very stable at room temperature and are difficult to fully ionize compared to separate dopants, and thus they are detrimental to carrier concentration improvement. This allows us to uncover mysteries of the doping bottleneck in hematite and provide guidance for optimizing doping and carrier conductivity in polaronic oxides toward highly efficient energy conversion applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

H-Atom Assignment and Sb–O Bonding of [Mes 3 SbOH][O 3 SPh] Confirmed by Neutron Diffraction, Multipole Modeling, and Hirshfeld Atom Refinement

Neutron wavelength-resolved Laue diffraction experiments permit accurate refinement of the H-atom positions and anisotropic displacement parameters of [Mes 3 SbOH][O 3 SPh]. A multipole-based charge density refinement and a topological analysis of the refined electron density were also performed. Hirshfeld atom refinement (HAR) recovers the neutron-determined H-atom parameters, and the quantum-mechanical electron density used in HAR recovers the electron density topology from the refined multipole model. We report these results confirm that [Mes 3 SbOH][O 3 SPh] does indeed feature a hydroxystibonium cation with a nominal Sb–O single bond and not a stibine oxide with an Sb=O/Sb + –O – bond.

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

High-order elastic multipoles as colloidal atoms

Achieving and exceeding diversity of colloidal analogs of chemical elements and molecules as building blocks of matter has been the central goal and challenge of colloidal science ever since Einstein introduced the colloidal atom paradigm. Recent advances in colloids assembly have been achieved by exploiting the machinery of DNA hybridization but robust physical means of defining colloidal elements remain limited. Here we introduce physical design principles allowing us to define high-order elastic multipoles emerging when colloids with controlled shapes and surface alignment are introduced into a nematic host fluid. Combination of experiments and numerical modeling of equilibrium field configurations using a spherical harmonic expansion allow us to probe elastic multipole moments, bringing analogies with electromagnetism and a structure of atomic orbitals. We illustrate that, at least in view of the symmetry of the “director wiggle wave functions,” diversity of elastic colloidal atoms can far exceed that of known chemical elements.

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

Higher-order electric multipole contributions to retarded non-additive three-body dispersion interaction energies between atoms: Equilateral triangle and collinear configurations

The theory of molecular quantum electrodynamics (QED) is used to calculate higher electric multipole contributions to the dispersion energy shift between three atoms or molecules arranged in a straight line or in an equilateral triangle configuration. As in two-body potentials, three-body dispersion interactions are viewed in the QED formalism to arise from exchange of virtual photons between coupled pairs of particles. By employing an interaction Hamiltonian that is quadratic in the electric displacement field means that third-order perturbation theory can be used to yield the energy shift for a particular combination of electric multipole polarizable species, with only six time-ordered diagrams needing to be summed over. Specific potentials evaluated include dipole-dipole-quadrupole (DDQ), dipole-quadrupole-quadrupole (DQQ), and dipole-dipole-octupole (DDO) terms. For the geometries of interest, near-zone limiting forms are found to exhibit an R{sup −11} dependence on separation distance for the DDQ interaction, and an R{sup −13} behaviour for DQQ and DDO shifts, agreeing with an earlier semi-classical computation. Retardation weakens the potential in each case by R{sup −1} in the far-zone. It is found that by decomposing the octupole moment into its irreducible components of weights-1 and -3 that the former contribution to the DDO potential may be taken to be a higher-order correction to the leading triple dipole energy shift.

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

Novel multipole Wien filter as three-dimensional spin manipulator

Spin polarized electron beam is often used in material characterizations which relates to magnetism as well as in the high energy particle physics. The manipulation of the spin polarization toward the arbitrary direction is indispensable in such studies. In the present work, a novel multipole Wien filter is proposed as the three-dimensional spin manipulator, and a prototype 8-pole Wien filter is developed. It is applied to spin polarized low energy electron microscopy, and the variation of the magnetic contrast with managing the spin polarization is evaluated. It is confirmed that the novel multipole Wien filter can manipulate the spin polarization three-dimensionally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Exchange splitting of the interaction energy and the multipole expansion of the wave function

The exchange splitting J of the interaction energy of the hydrogen atom with a proton is calculated using the conventional surface-integral formula J{sub surf}[Φ], the volume-integral formula of the symmetry-adapted perturbation theory J{sub SAPT}[Φ], and a variational volume-integral formula J{sub var}[Φ]. The calculations are based on the multipole expansion of the wave function Φ, which is divergent for any internuclear distance R. Nevertheless, the resulting approximations to the leading coefficient j{sub 0} in the large-R asymptotic series J(R) = 2e{sup −R−1}R(j{sub 0} + j{sub 1}R{sup −1} + j{sub 2}R{sup −2} + ⋯) converge with the rate corresponding to the convergence radii equal to 4, 2, and 1 when the J{sub var}[Φ], J{sub surf}[Φ], and J{sub SAPT}[Φ] formulas are used, respectively. Additionally, we observe that also the higher j{sub k} coefficients are predicted correctly when the multipole expansion is used in the J{sub var}[Φ] and J{sub surf}[Φ] formulas. The symmetry adapted perturbation theory formula J{sub SAPT}[Φ] predicts correctly only the first two coefficients, j{sub 0} and j{sub 1}, gives a wrong value of j{sub 2}, and diverges for higher j{sub n}. Since the variational volume-integral formula can be easily generalized to many-electron systems and evaluated with standard basis-set techniques of quantum chemistry, it provides an alternative for the determination of the exchange splitting and the exchange contribution of the interaction potential in general.

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

Linearly scaling and almost Hamiltonian dielectric continuum molecular dynamics simulations through fast multipole expansions

Hamiltonian Dielectric Solvent (HADES) is a recent method [S. Bauer et al., J. Chem. Phys. 140, 104103 (2014)] which enables atomistic Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric solvent continua. Such simulations become rapidly impractical for large proteins, because the computational effort of HADES scales quadratically with the number N of atoms. If one tries to achieve linear scaling by applying a fast multipole method (FMM) to the computation of the HADES electrostatics, the Hamiltonian character (conservation of total energy, linear, and angular momenta) may get lost. Here, we show that the Hamiltonian character of HADES can be almost completely preserved, if the structure-adapted fast multipole method (SAMM) as recently redesigned by Lorenzen et al. [J. Chem. Theory Comput. 10, 3244-3259 (2014)] is suitably extended and is chosen as the FMM module. By this extension, the HADES/SAMM forces become exact gradients of the HADES/SAMM energy. Their translational and rotational invariance then guarantees (within the limits of numerical accuracy) the exact conservation of the linear and angular momenta. Also, the total energy is essentially conserved—up to residual algorithmic noise, which is caused by the periodically repeated SAMM interaction list updates. These updates entail very small temporal discontinuities of the force description, because the employed SAMM approximations represent deliberately balanced compromises between accuracy and efficiency. The energy-gradient corrected version of SAMM can also be applied, of course, to MD simulations of all-atom solvent-solute systems enclosed by periodic boundary conditions. However, as we demonstrate in passing, this choice does not offer any serious advantages.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗