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

An adaptable multiple power source for mass spectrometry and other scientific instruments

Power supplies are commonly used in the operation of many types of scientific equipment, including mass spectrometers and ancillary instrumentation. A generic modern mass spectrometer comprises an ionization source, such as electrospray ionization (ESI), ion transfer devices such as ion funnels and multipole ion guides, and ion signal detection apparatus. Very often such platforms include, or are interfaced with ancillary elements in order to manipulate samples before or after ionization. In order to operate such scientific instruments, numerous direct current (DC) channels and radio frequency (RF) signals are required, along with other controls such as temperature regulation. In particular, DC voltages in the range of ±400 V, along with MHz range RF signals with peak-to-peak amplitudes in the hundreds of volts range are commonly used to transfer ionized samples under vacuum. Additionally, an ESI source requires a high voltage (HV) DC source capable of producing several thousand volts and heaters capable of generating temperatures up to 300°C. All of these signals must be properly synchronized and managed in order to carry out ion trapping, accumulation and detection.

47 OTHER INSTRUMENTATION↗
DOE OSTI2015

Hyper-Rayleigh scattering in centrosymmetric systems

Hyper-Rayleigh scattering (HRS) is an incoherent mechanism for optical second harmonic generation. The frequency-doubled light that emerges from this mechanism is not emitted in a laser-like manner, in the forward direction; it is scattered in all directions. The underlying theory for this effect involves terms that are quadratic in the incident field and involves an even-order optical susceptibility (for a molecule, its associated hyperpolarizability). In consequence, HRS is often regarded as formally forbidden in centrosymmetric media. However, for the fundamental three-photon interaction, theory based on the standard electric dipole approximation, representable as E1{sup 3}, does not account for all experimental observations. The relevant results emerge upon extending the theory to include E1{sup 2}M1 and E1{sup 2}E2 contributions, incorporating one magnetic dipolar or electric quadrupolar interaction, respectively, to a consistent level of multipolar expansion. Both additional interactions require the deployment of higher orders in the multipole expansion, with the E1{sup 2}E2 interaction analogous in rank and parity to a four-wave susceptibility. To elicit the correct form of response from fluid or disordered media invites a tensor representation which does not oversimplify the molecular components, yet which can produce results to facilitate the interpretation of experimental observations. The detailed derivation in this work leads to results which are summarized for the following: perpendicular detection of polarization components both parallel and perpendicular to the pump radiation, leading to distinct polarization ratio results, as well as a reversal ratio for forward scattered circular polarizations. The results provide a route to handling data with direct physical interpretation, to enable the more sophisticated design of molecules with sought nonlinear optical properties.

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

Mean force on a finite-sized rigid particle, droplet, or bubble in a viscous compressible medium

In this paper, a force formulation to compute the axial acoustic mean second-order force on finite-sized compressible and rigid particles is presented. The flow inside and outside the spherical inclusion is considered viscous and compressible. Other than for volumetric pulsations of the bubble/droplet, the sphericity of the inclusion is maintained (taken to be unity). A far-field derivation approach has been used to compute the force due to standing and traveling waves; and the force is expressed as a multipole expansion (infinite series). In case of a bubble and a rigid particle, there exist three length scales that govern the mean second-order force: mean radius of the spherical inclusion (R 0 ), wavelength of the incoming acoustics (λ), and the momentum diffusion thickness of the ambient fluid (δ o ). While R 0 and λ are arbitrary, we assume the viscous length scale is negligibly small compared to the acoustic wavelength. In case of a droplet, however, the following additional parameters (inside to outside fluid ratios) also play a role: density ratio ($\sim\atop{p}$), viscosity ratio ($\sim\atop{μ}$, and speed of sound ratio ($\sim\atop{c}$). The force expression yields the correct behavior in several limiting cases considered: (i) inviscid bubble and droplets with R 0 /λ « 1, (ii) inviscid bubbles with finite R 0 /λ, and (iii) finite size rigid immovable particles. In general, while the monopole alone is sufficient to capture the force for small bubbles, higher-order terms are found to be important when R 0 /λ ≥ 0.02. In addition to reporting similar behavior for droplets, we study the effect arising from $\sim\atop{p}$, $\sim\atop{μ}$, $\sim\atop{c}$, and δ o /R 0 on the mean second-order force.

42 ENGINEERING↗
DOE OSTI2015

Non-linear non-local molecular electrodynamics with nano-optical fields

The interaction of optical fields sculpted on the nano-scale with matter may not be described by the dipole approximation since the fields may vary appreciably across the molecular length scale. Rather than incrementally adding higher multipoles, it is advantageous and more physically transparent to describe the optical process using non-local response functions that intrinsically include all multi-poles. In this work, we present a semi-classical approach for calculating non-local response functions based on the minimal coupling Hamiltonian. The first, second, and third order response functions are expressed in terms of correlation functions of the charge and the current densities. Finally, this approach is based on the gauge invariant current rather than the polarization, and on the vector potential rather than the electric and magnetic fields.

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

Low energy, high power hydrogen neutral beam for plasma heating

A high power, relatively low energy neutral beam injector was developed to upgrade of the neutral beam system of the gas dynamic trap device and C2-U experiment. The ion source of the injector produces a proton beam with the particle energy of 15 keV, current of up to 175 A, and pulse duration of a few milliseconds. The plasma emitter of the ion source is produced by superimposing highly ionized plasma jets from an array of four arc-discharge plasma generators. A multipole magnetic field produced with permanent magnets at the periphery of the plasma box is used to increase the efficiency and improve the uniformity of the plasma emitter. Multi-slit grids with 48% transparency are fabricated from bronze plates, which are spherically shaped to provide geometrical beam focusing. The focal length of the Ion Optical System (IOS) is 3.5 m and the initial beam diameter is 34 cm. The IOS geometry and grid potentials were optimized numerically to ensure accurate beam formation. The measured angular divergences of the beam are ±0.01 rad parallel to the slits and ±0.03 rad in the transverse direction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

The consequences of improperly describing oscillator strengths beyond the electric dipole approximation

The interaction between a quantum mechanical system and plane wave light is usually modeled within the electric dipole approximation. This assumes that the intensity of the incident field is constant over the length of the system and transition probabilities are described in terms of the electric dipole transition moment. For short wavelength spectroscopies, such as X-ray absorption, the electric dipole approximation often breaks down. Higher order multipoles are then included to describe transition probabilities. The square of the magnetic dipole and electric quadrupole are often included, but this results in an origin-dependent expression for the oscillator strength. The oscillator strength can be made origin-independent if all terms through the same order in the wave vector are retained. We will show the consequences and potential pitfalls of using either of these two expressions. It is shown that the origin-dependent expression may violate the Thomas-Reiche-Kuhn sum rule and the origin-independent expression can result in negative transition probabilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Puzzles in quarkonium hadronic transitions with two pion emission

The anomalously large rates of some hadronic transitions from quarkonium are studied using QCD multipole expansion (QCDME) in the framework of a constituent quark model which has been successful in describing hadronic phenomenology. The hybrid intermediate states needed in the QCDME method are calculated in a natural extension of our constituent quark model based on the Quark Confining String (QCS) scheme. Some of the anomalies are explained due to the presence of an hybrid state with a mass near the mass of the decaying resonance whereas other are justified by the presence of molecular components in the wave function. Some unexpected results are pointed out.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Multilevel summation with B-spline interpolation for pairwise interactions in molecular dynamics simulations

The multilevel summation method for calculating electrostatic interactions in molecular dynamics simulations constructs an approximation to a pairwise interaction kernel and its gradient, which can be evaluated at a cost that scales linearly with the number of atoms. The method smoothly splits the kernel into a sum of partial kernels of increasing range and decreasing variability with the longer-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate in the context of dynamics and minimization, because it can produce continuous gradients. This article explores the use of B-splines to increase the accuracy of the multilevel summation method (for nonperiodic boundaries) without incurring additional computation other than a preprocessing step (whose cost also scales linearly). To obtain accurate results efficiently involves technical difficulties, which are overcome by a novel preprocessing algorithm. Numerical experiments demonstrate that the resulting method offers substantial improvements in accuracy and that its performance is competitive with an implementation of the fast multipole method in general and markedly better for Hamiltonian formulations of molecular dynamics. The improvement is great enough to establish multilevel summation as a serious contender for calculating pairwise interactions in molecular dynamics simulations. In particular, the method appears to be uniquely capable for molecular dynamics in two situations, nonperiodic boundary conditions and massively parallel computation, where the fast Fourier transform employed in the particle–mesh Ewald method falls short.

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

Electronic structure of an anticancer drug DC81 and its interaction with DNA base pairs

The drug, 8-Hydroxy-7-methoxy-pyrrolo-[2,1-c][1,4] benzodiazepine-5-one, commonly christened as DC81 belongs to the pyrrolo-[2,1-c][1,4]benzodiazepine (PBDs) family. It is a member of the group of naturally occurring antitumour antibiotics produced by various Streptomyces species. The antitumour activity of DC81 is attributed to its sequence specific interaction with G-C rich DNA region in particular, for Pu-G-Pu motifs. In the present paper, physico-chemical properties DC81 have been carried out using an ab-initio method, HF/6-31G(d,p) with GAMESS program. MEP, HOMO and LUMO surfaces have been scanned. Ionization potential, electron affinity, electronegativity, global hardness and softness of the drug have been calculated. Further, drug-DNA interactions have been examined using modified second order perturbation theory along with multicentred-multipole expansion technique. Results have been discussed in the light of other theoretical and experimental observations. Efforts have been made to elucidate the binding patterns and thereby biological properties of the drug.

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

VR-SCOSMO: A smooth conductor-like screening model with charge-dependent radii for modeling chemical reactions

To better represent the solvation effects observed along reaction pathways, and of ionic species in general, a charge-dependent variable-radii smooth conductor-like screening model (VR-SCOSMO) is developed. This model is implemented and parameterized with a third order density-functional tight binding quantum model, DFTB3/3OB-OPhyd, a quantum method which was developed for organic and biological compounds, utilizing a specific parameterization for phosphate hydrolysis reactions. Unlike most other applications with the DFTB3/3OB model, an auxiliary set of atomic multipoles is constructed from the underlying DFTB3 density matrix which is used to interact the solute with the solvent response surface. The resulting method is variational, produces smooth energies, and has analytic gradients. As a baseline, a conventional SCOSMO model with fixed radii is also parameterized. The SCOSMO and VR-SCOSMO models shown have comparable accuracy in reproducing neutral-molecule absolute solvation free energies; however, the VR-SCOSMO model is shown to reduce the mean unsigned errors (MUEs) of ionic compounds by half (about 2-3 kcal/mol). The VR-SCOSMO model presents similar accuracy as a charge-dependent Poisson-Boltzmann model introduced by Hou et al. [J. Chem. Theory Comput. 6, 2303 (2010)]. VR-SCOSMO is then used to examine the hydrolysis of trimethylphosphate and seven other phosphoryl transesterification reactions with different leaving groups. Two-dimensional energy landscapes are constructed for these reactions and calculated barriers are compared to those obtained from ab initio polarizable continuum calculations and experiment. Results of the VR-SCOSMO model are in good agreement in both cases, capturing the rate-limiting reaction barrier and the nature of the transition state.

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

Transfer of energy between a pair of molecules near a plasmonic core-shell nanoparticle: Tunability and sensing

Our model is applied to the calculation of interaction energy between a pair of dipolar molecules (point dipoles) in the vicinity of a nanoshell monomer with core-shell structure, based on the dipole quasi-electrostatic theory of classical electrodynamics and using the Drude and Maxwell-Garnett model. In other words, this work discusses the intermolecular energy transfer from a donor molecule to an acceptor molecule near a spherical nanoparticle that is important for practical applications like sensing. It is shown that the proximity of plasmonic nanoparticles can have a strong effect on the energy transfer between molecules. In addition to the influence of the size, composition, embedding medium, and the filling fraction of doped particles on the interaction energy, the contribution of the dipolar, quadrupolar, octupolar, hexadecapolar, triakontadipolar, and higher order multipole interactions is presented and analyzed. Briefly, we will show that it is possible to achieve enhanced energy transfer by manipulation of different parameters as mentioned above.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Understanding the many-body expansion for large systems. III. Critical role of four-body terms, counterpoise corrections, and cutoffs

Papers I and II in this series [R. M. Richard et al., J. Chem. Phys. 141, 014108 (2014); K. U. Lao et al., ibid. 144, 164105 (2016)] have attempted to shed light on precision and accuracy issues affecting the many-body expansion (MBE), which only manifest in larger systems and thus have received scant attention in the literature. Many-body counterpoise (CP) corrections are shown to accelerate convergence of the MBE, which otherwise suffers from a mismatch between how basis-set superposition error affects subsystem versus supersystem calculations. In water clusters ranging in size up to (H2O)37, four-body terms prove necessary to achieve accurate results for both total interaction energies and relative isomer energies, but the sheer number of tetramers makes the use of cutoff schemes essential. To predict relative energies of (H2O)20 isomers, two approximations based on a lower level of theory are introduced and an ONIOM-type procedure is found to be very well converged with respect to the appropriate MBE benchmark, namely, a CP-corrected supersystem calculation at the same level of theory. Results using an energy-based cutoff scheme suggest that if reasonable approximations to the subsystem energies are available (based on classical multipoles, say), then the number of requisite subsystem calculations can be reduced even more dramatically than when distance-based thresholds are employed. The end result is several accurate four-body methods that do not require charge embedding, and which are stable in large basis sets such as aug-cc-pVTZ that have sometimes proven problematic for fragment-based quantum chemistry methods. Even with aggressive thresholding, however, the four-body approach at the self-consistent field level still requires roughly ten times more processors to outmatch the performance of the corresponding supersystem calculation, in test cases involving 1500–1800 basis functions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

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

Monte Carlo transport analysis to assess intensity dependent response of a carbon-doped GaN photoconductor

In this report evaluation of the photoresponse in wurtzite GaN photoconductive switches is presented based on kinetic Monte Carlo simulations. The focus is on electron transport physics and assessment of high frequency operation. The roles of GaN band structure, Pauli exclusion, and treatment of internal fields based on the fast multipole method are all comprehensively included. The implementation was validated through comparisons of velocity-field characteristics for GaN with computational results in the literature. Photocurrent widths of less than ~7 ps for the 1 μm device can be expected, which translates into a 100 GHz upper bound. Photocurrent pulse compression below the laser full width at half maxima at high applied fields are predicted based on the interplay of space-charge effects and the negative differential velocity characteristics of GaN.

36 MATERIALS SCIENCE↗
DOE OSTI2021

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

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

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

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

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

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

Flowing plasma rearrangement in the presence of static perturbing fields

Charged particles interacting with electromagnetic waves have a portion of their energy tied up in wave-driven oscillations. When these waves are localized to the exhaust of linear magnetic confinement systems, this ponderomotive effect can be utilized to enhance particle confinement. The same effect can be derived for particles moving via an E×B drift into a region of a static perturbation to the electromagnetic fields which has a large wave vector component in the direction of the motion. In this work, we use a simplified slab model to self-consistently solve for the electromagnetic fields within the fluid flowing plasma of a static flute-like (k∥=0) perturbation and evaluate the resulting ponderomotive potential. We find that two types of perturbations can exist within the flowing plasma, which are an O wave and an X wave in the frame moving with the fluid. In the case of tenuous plasma, these perturbations are magnetostatic or electrostatic multipole-analog perpendicular to the guiding magnetic field in the lab frame, respectfully. For denser plasmas, the O wave-like perturbation is screened at the electron skin depth scale, and the X wave-like perturbation is a combination of a similar perpendicular electric perturbation and parallel magnetic perturbation. The ponderomotive potential generated in the X wave-like case is gyrofrequency-dependent and can be used as either potential barriers or potential wells, depending on the direction of the flow velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗