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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
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 OSTI2014

Beyond first-order finite element schemes in micromagnetics

Magnetization dynamics in ferromagnetic materials is ruled by the Landau–Lifshitz–Gilbert equation (LLG). Reliable schemes must conserve the magnetization norm, which is a nonconvex constraint, and be energy-decreasing unless there is pumping. Some of the authors previously devised a convergent finite element scheme that, by choice of an appropriate test space – the tangent plane to the magnetization – reduces to a linear problem at each time step. The scheme was however first-order in time. We claim it is not an intrinsic limitation, and the same approach can lead to efficient micromagnetic simulation. We show how the scheme order can be increased, and the nonlocal (magnetostatic) interactions be tackled in logarithmic time, by the fast multipole method or the non-uniform fast Fourier transform. Our implementation is called feeLLGood. A test-case of the National Institute of Standards and Technology is presented, then another one relevant to spin-transfer effects (the spin-torque oscillator)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Accuracy and efficiency in computing electrostatic potential for an ion channel model in layered dielectric/electrolyte media

This paper will investigate the numerical accuracy and efficiency in computing the electrostatic potential for a finite-height cylinder, used in an explicit/implicit hybrid solvation model for ion channel and embedded in a layered dielectric/electrolyte medium representing a biological membrane and ionic solvents. A charge locating inside the cylinder cavity, where ion channel proteins and ions are given explicit atomistic representations, will be influenced by the polarization field of the surrounding implicit dielectric/electrolyte medium. Two numerical techniques, a specially designed boundary integral equation method and an image charge method, will be investigated and compared in terms of accuracy and efficiency for computing the electrostatic potential. The boundary integral equation method based on the three-dimensional layered Green's functions provides a highly accurate solution suitable for producing a benchmark reference solution, while the image charge method is found to give reasonable accuracy and highly efficient and viable to use the fast multipole method for interactions of a large number of charges in the atomistic region of the hybrid solvation model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A harmonic polynomial cell (HPC) method for 3D Laplace equation with application in marine hydrodynamics

We propose a new efficient and accurate numerical method based on harmonic polynomials to solve boundary value problems governed by 3D Laplace equation. The computational domain is discretized by overlapping cells. Within each cell, the velocity potential is represented by the linear superposition of a complete set of harmonic polynomials, which are the elementary solutions of Laplace equation. By its definition, the method is named as Harmonic Polynomial Cell (HPC) method. The characteristics of the accuracy and efficiency of the HPC method are demonstrated by studying analytical cases. Comparisons will be made with some other existing boundary element based methods, e.g. Quadratic Boundary Element Method (QBEM) and the Fast Multipole Accelerated QBEM (FMA-QBEM) and a fourth order Finite Difference Method (FDM). To demonstrate the applications of the method, it is applied to some studies relevant for marine hydrodynamics. Sloshing in 3D rectangular tanks, a fully-nonlinear numerical wave tank, fully-nonlinear wave focusing on a semi-circular shoal, and the nonlinear wave diffraction of a bottom-mounted cylinder in regular waves are studied. The comparisons with the experimental results and other numerical results are all in satisfactory agreement, indicating that the present HPC method is a promising method in solving potential-flow problems. The underlying procedure of the HPC method could also be useful in other fields than marine hydrodynamics involved with solving Laplace equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Solution of 3D inverse scattering problems by combined inverse equivalent current and finite element methods

An approach combining boundary integral and finite element methods is introduced for the solution of three-dimensional inverse electromagnetic medium scattering problems. Based on the equivalence principle, unknown equivalent electric and magnetic surface current densities on a closed surface are utilized to decompose the inverse medium problem into two parts: a linear radiation problem and a nonlinear cavity problem. The first problem is formulated by a boundary integral equation, the computational burden of which is reduced by employing the multilevel fast multipole method (MLFMM). Reconstructed Cauchy data on the surface allows the utilization of the Lorentz reciprocity and the Poynting's theorems. Exploiting these theorems, the noise level and an initial guess are estimated for the cavity problem. Moreover, it is possible to determine whether the material is lossy or not. In the second problem, the estimated surface currents form inhomogeneous boundary conditions of the cavity problem. The cavity problem is formulated by the finite element technique and solved iteratively by the Gauss–Newton method to reconstruct the properties of the object. Regularization for both the first and the second problems is achieved by a Krylov subspace method. The proposed method is tested against both synthetic and experimental data and promising reconstruction results are obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

An integral equation formulation for rigid bodies in Stokes flow in three dimensions

Here we present a new derivation of a boundary integral equation (BIE) for simulating the three-dimensional dynamics of arbitrarily-shaped rigid particles of genus zero immersed in a Stokes fluid, on which are prescribed forces and torques. Our method is based on a single-layer representation and leads to a simple second-kind integral equation. It avoids the use of auxiliary sources within each particle that play a role in some classical formulations. We use a spectrally accurate quadrature scheme to evaluate the corresponding layer potentials, so that only a small number of spatial discretization points per particle are required. The resulting discrete sums are computed in $\mathscr{O}(n)$ time, where n denotes the number of particles, using the fast multipole method (FMM). The particle positions and orientations are updated by a high-order time-stepping scheme. We illustrate the accuracy, conditioning and scaling of our solvers with several numerical examples.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2013

Nuclear Data Sheets for A = 75

The experimental nuclear spectroscopic data for known nuclides of mass number 75 (Fe, Co, Ni, Cu, Zn, Ga, Ge, As, Se, Br, Kr, Rb, Sr) have been evaluated and presented together with Adopted properties for levels and γ rays. New high-spin data are available for {sup 75}Ga, and {sup 75}Rb; and lifetime data for high-spin states in {sup 75}Br and {sup 75}Kr. For ΔJ =1, M1+E2 transitions in two rotational bands in {sup 75}Kr, several B(E2)(W.u.) values are anomalously high, deviating by 2–3 σ values from currently accepted RUL(E2) = 300. In the opinion of the evaluators, there is need to remeasure level lifetimes and multipole mixing ratios in {sup 75}Kr to resolve this serious discrepancy. New precise single-particle transfer cross section data are available for {sup 75}Ga, {sup 75}Ge, {sup 75}As and {sup 75}Se from several different reactions (2009Ka06,2008Sc03); these data give information for occupancy of valence neutron orbitals in the ground states of target nuclides: {sup 76}Ge, {sup 76}Se and {sup 78}Se. No significant new data since the 1999 NDS for A = 75 have been reported for {sup 75}As and {sup 75}Se. No data are yet available for excited states in {sup 75}Co, {sup 75}Ni and {sup 75}Sr. For {sup 75}Fe, only the isotopic identification is made with one observed event. The radioactive decay schemes of {sup 75}Co and {sup 75}Ni are unknown while those for {sup 75}Rb and {sup 75}Sr are incomplete. This work supersedes the data presented in the previous NDS evaluation of A = 75 published by 1999Fa05.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Nuclear Data Sheets for A = 248

Detailed level schemes, decay schemes, and the experimental data on which they are based are presented for all nuclei with mass number A=248. The experimental reaction and decay data are evaluated; inconsistencies and discrepancies are noted; and adopted values are given for various level properties (including spin and parity, halflife, and configuration), γ properties (including energy, intensity, and multipole character) and α properties (including energy, intensity, and hindrance factor). This evaluation replaces the A=248 evaluation published by A. Akovali in Nuclear Data Sheets 87, 249 (1999)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2017

Aberration corrected STEM by means of diffraction gratings

In the past 15 years, the advent of aberration correction technology in electron microscopy has enabled materials analysis on the atomic scale. This is made possible by precise arrangements of multipole electrodes and magnetic solenoids to compensate the aberrations inherent to any focusing element of an electron microscope. In this paper, we describe an alternative method to correct for the spherical aberration of the objective lens in scanning transmission electron microscopy (STEM) using a passive, nanofabricated diffractive optical element. This holographic device is installed in the probe forming aperture of a conventional electron microscope and can be designed to remove arbitrarily complex aberrations from the electron's wave front. In this work, we show a proof-of-principle experiment that demonstrates successful correction of the spherical aberration in STEM by means of such a grating corrector (GCOR). Our GCOR enables us to record aberration-corrected high-resolution high-angle annular dark field (HAADF-) STEM images, although yet without advancement in probe current and resolution. Finally, improvements in this technology could provide an economical solution for aberration-corrected high-resolution STEM in certain use scenarios.

47 OTHER INSTRUMENTATION↗
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↗