Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MULTIPOLAR FIELD”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Magnetar giant flares in multipolar magnetic fields. II. Flux rope eruptions with current sheets

We propose a physical mechanism to explain giant flares and radio afterglows in terms of a magnetospheric model containing both a helically twisted flux rope and a current sheet (CS). With the appearance of a CS, we solve a mixed boundary value problem to get the magnetospheric field based on a domain decomposition method. We investigate properties of the equilibrium curve of the flux rope when the CS is present in background multipolar fields. In response to the variations at the magnetar surface, it quasi-statically evolves in stable equilibrium states. The loss of equilibrium occurs at a critical point and, beyond that point, it erupts catastrophically. New features show up when the CS is considered. In particular, we find two kinds of physical behaviors, i.e., catastrophic state transition and catastrophic escape. Magnetic energy would be released during state transitions. This released magnetic energy is sufficient to drive giant flares, and the flux rope would, therefore, go away from the magnetar quasi-statically, which is inconsistent with the radio afterglow. Fortunately, in the latter case, i.e., the catastrophic escape, the flux rope could escape the magnetar and go to infinity in a dynamical way. This is more consistent with radio afterglow observations of giant flares. We find that the minor radius of the flux rope has important implications for its eruption. Flux ropes with larger minor radii are more prone to erupt. We stress that the CS provides an ideal place for magnetic reconnection, which would further enhance the energy release during eruptions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Tailoring dielectric resonator geometries for directional scattering and Huygens’ metasurfaces

In this paper we describe a methodology for tailoring the design of metamaterial dielectric resonators, which represent a promising path toward low-loss metamaterials at optical frequencies. We first describe a procedure to decompose the far field scattered by subwavelength resonators in terms of multipolar field components, providing explicit expressions for the multipolar far fields. We apply this formulation to confirm that an isolated high-permittivity dielectric cube resonator possesses frequency separated electric and magnetic dipole resonances, as well as a magnetic quadrupole resonance in close proximity to the electric dipole resonance. We then introduce multiple dielectric gaps to the resonator geometry in a manner suggested by perturbation theory, and demonstrate the ability to overlap the electric and magnetic dipole resonances, thereby enabling directional scattering by satisfying the first Kerker condition. We further demonstrate the ability to push the quadrupole resonance away from the degenerate dipole resonances to achieve local behavior. These properties are confirmed through the multipolar expansion and show that the use of geometries suggested by perturbation theory is a viable route to achieve purely dipole resonances for metamaterial applications such as wave-front manipulation with Huygens’ metasurfaces. Our results are fully scalable across any frequency bands where high-permittivity dielectric materials are available, including microwave, THz, and infrared frequencies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The Multipolar Magnetic Field of the Millisecond Pulsar PSR J0030+0451

Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed at constraining the neutron star mass and radius, has inferred surface hot spots (the magnetic polar caps) that imply significantly nondipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole-plus-quadrupole components using a static vacuum field and force-free global magnetosphere models. Taking into account the compactness and observer angle values provided by Miller et al. and Riley et al., we compute geodesics from the observer plane to the polar caps to compute the resulting X-ray light curve. We explore, through Markov Chain Monte Carlo techniques, the detailed magnetic field configurations that can reproduce the observed X-ray light curve and have discovered degeneracies, i.e., diverse field configurations, which can provide sufficient descriptions of the NICER X-ray waveforms. Having obtained the force-free field structures, we then compute the corresponding synchronous γ-ray light curves following Kalapotharakos et al.; these we compare to those obtained by Fermi-LAT, to provide models consistent with both the X-ray and the γ-ray data, thereby restricting further the multipole field parameters. An essential aspect of this approach is the proper computation of the relative phase between the synchronous X- and γ-ray light curves. We conclude with a discussion of the broader implications of our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A nonlinear eigenvalue problem for self-similar spherical force-free magnetic fields

An axisymmetric force-free magnetic field B(r, θ) in spherical coordinates is defined by a function r sin θB{sub φ}=Q(A) relating its azimuthal component to its poloidal flux-function A. The power law r sin θB{sub φ}=aA|A|{sup 1/n}, n a positive constant, admits separable fields with A=(A{sub n}(θ))/(r{sup n}) , posing a nonlinear boundary-value problem for the constant parameter a as an eigenvalue and A{sub n}(θ) as its eigenfunction [B. C. Low and Y. Q Lou, Astrophys. J. 352, 343 (1990)]. A complete analysis is presented of the eigenvalue spectrum for a given n, providing a unified understanding of the eigenfunctions and the physical relationship between the field's degree of multi-polarity and rate of radial decay via the parameter n. These force-free fields, self-similar on spheres of constant r, have basic astrophysical applications. As explicit solutions they have, over the years, served as standard benchmarks for testing 3D numerical codes developed to compute general force-free fields in the solar corona. The study presented includes a set of illustrative multipolar field solutions to address the magnetohydrodynamics (MHD) issues underlying the observation that the solar corona has a statistical preference for negative and positive magnetic helicities in its northern and southern hemispheres, respectively; a hemispherical effect, unchanging as the Sun's global field reverses polarity in successive eleven-year cycles. Generalizing these force-free fields to the separable form B=(H(θ,φ))/(r{sup n+2}) promises field solutions of even richer topological varieties but allowing for φ-dependence greatly complicates the governing equations that have remained intractable. The axisymmetric results obtained are discussed in relation to this generalization and the Parker Magnetostatic Theorem. The axisymmetric solutions are mathematically related to a family of 3D time-dependent ideal MHD solutions for a polytropic fluid of index γ = 4/3 as discussed in the Appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Critical Height of the Torus Instability in Two-ribbon Solar Flares

We studied the background field for 60 two-ribbon flares of M-and-above classes during 2011–2015. These flares are categorized into two groups, i.e., eruptive and confined flares, based on whether a flare is associated with a coronal mass ejection or not. The background field of source active regions is approximated by a potential field extrapolated from the B {sub z} component of vector magnetograms provided by the Helioseismic and Magnetic Imager. We calculated the decay index n of the background field above the flaring polarity inversion line, and defined a critical height h {sub crit} corresponding to the theoretical threshold ( n {sub crit} = 1.5) of the torus instability. We found that h {sub crit} is approximately half of the distance between the centroids of opposite polarities in active regions and that the distribution of h {sub crit} is bimodal: it is significantly higher for confined flares than for eruptive ones. The decay index increases monotonously with increasing height for 86% (84%) of the eruptive (confined) flares but displays a saddle-like profile for the rest, 14% (16%), which are found exclusively in active regions of multipolar field configuration. Moreover, n at the saddle bottom is significantly smaller in confined flares than that in eruptive ones. These results highlight the critical role of background field in regulating the eruptive behavior of two-ribbon flares.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Rotating-Coil Measurement System for Small-Bore-Diameter Magnet Characterization

Rotating-coil measurement systems are widely used to measure the multipolar fields of particle accelerator magnets. This paper presents a rotating-coil measurement system that aims at providing a complete data set for the characterization of quadrupole magnets with small bore diameters (26 mm). The PCB magnetometer design represents a challenging goal for this type of transducer. It is characterized by an aspect ratio 30% higher than the state of the art, imposed by the reduced dimension of the external radius of the rotating shaft and the necessity of covering the entire magnet effective length (500 mm or higher). The system design required a novel design for the mechanical asset, also considering the innovation represented by the commercial carbon fiber tube, housing the PCB magnetometer. Moreover, the measurement system is based primarily on standard and commercially available components, with simplified control and post-processing software applications. The system and its components are cross-calibrated using a stretched-wire system and another rotating-coil system. The measurement precision is established in a measurement campaign performed on a quadrupole magnet characterized by an inner bore diameter of 45 mm.

42 ENGINEERING↗
DOE OSTI2015

Solvation of fluoro-acetonitrile in water by 2D-IR spectroscopy: A combined experimental-computational study

The solvent dynamics around fluorinated acetonitrile is characterized by 2-dimensional infrared spectroscopy and atomistic simulations. The lineshape of the linear infrared spectrum is better captured by semiempirical (density functional tight binding) mixed quantum mechanical/molecular mechanics simulations, whereas force field simulations with multipolar interactions yield lineshapes that are significantly too narrow. For the solvent dynamics, a relatively slow time scale of 2 ps is found from the experiments and supported by the mixed quantum mechanical/molecular mechanics simulations. With multipolar force fields fitted to the available thermodynamical data, the time scale is considerably faster—on the 0.5 ps time scale. The simulations provide evidence for a well established CF–HOH hydrogen bond (population of 25%) which is found from the radial distribution function g(r) from both, force field and quantum mechanics/molecular mechanics simulations.

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

On the nature of the extreme-ultraviolet late phase of solar flares

The extreme-ultraviolet (EUV) late phase of solar flares is a second peak of warm coronal emissions (e.g., Fe XVI) for many minutes to a few hours after the GOES soft X-ray peak. It was first observed by the EUV Variability Experiment on board the Solar Dynamics Observatory (SDO). The late-phase emission originates from a second set of longer loops (late-phase loops) that are higher than the main flaring loops. It is suggested to be caused by either additional heating or long-lasting cooling. In this paper, we study the role of long-lasting cooling and additional heating in producing the EUV late phase using the enthalpy based thermal evolution of loops model. We find that a long cooling process in late-phase loops can well explain the presence of the EUV late-phase emission, but we cannot exclude the possibility of additional heating in the decay phase. Moreover, we provide two preliminary methods based on the UV and EUV emissions from the Atmospheric Imaging Assembly on board SDO to determine whether or not additional heating plays a role in the late-phase emission. Using nonlinear force-free field modeling, we study the magnetic configuration of the EUV late phase. It is found that the late phase can be generated either in hot spine field lines associated with a magnetic null point or in large-scale magnetic loops of multipolar magnetic fields. In this paper, we also discuss why the EUV late phase is usually observed in warm coronal emissions and why the majority of flares do not exhibit an EUV late phase.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

AB INITIO PULSAR MAGNETOSPHERE: THREE-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF OBLIQUE PULSARS

We present “first-principles” relativistic particle-in-cell simulations of the oblique pulsar magnetosphere with pair formation. The magnetosphere starts to form with particles extracted from the surface of the neutron star. These particles are accelerated by surface electric fields and emit photons capable of producing electron–positron pairs. We inject secondary pairs at the locations of primary energetic particles whose energy exceeds the threshold for pair formation. We find solutions that are close to the ideal force-free magnetosphere with the Y-point and current sheet. Solutions with obliquities ≤40° do not show pair production in the open field line region because the local current density along the magnetic field is below the Goldreich–Julian value. The bulk outflow in these solutions is charge-separated, and pair formation happens in the current sheet and return current layer only. Solutions with higher inclinations show pair production in the open field line region, with high multiplicity of the bulk flow and the size of the pair-producing region increasing with inclination. We observe the spin-down of the star to be comparable to MHD model predictions. The magnetic dissipation in the current sheet ranges between 20% for the aligned rotator and 3% for the orthogonal rotator. Our results suggest that for low obliquity neutron stars with suppressed pair formation at the light cylinder, the presence of phenomena related to pair activity in the bulk of the polar region, e.g., radio emission, may crucially depend on the physics beyond our simplified model, such as the effects of curved spacetime or multipolar surface fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Electrical conductivity and magnetic dynamos in magma oceans of Super-Earths

Super-Earths are extremely common among the numerous exoplanets that have been discovered. The high pressures and temperatures in their interiors are likely to lead to long-lived magma oceans. If their electrical conductivity is sufficiently high, the mantles of Super-Earth would generate their own magnetic fields. With ab initio simulations, we show that upon melting, the behavior of typical mantle silicates changes from semi-conducting to semi-metallic. The electrical conductivity increases and the optical properties are substantially modified. Melting could thus be detected with high-precision reflectivity measurements during the short time scales of shock experiments. We estimate the electrical conductivity of mantle silicates to be of the order of 100 Ω -1 cm -1 , which implies that a magnetic dynamo process would develop in the magma oceans of Super-Earths if their convective velocities have typical values of 1 mm/s or higher. We predict exoplanets with rotation periods longer than 2 days to have multipolar magnetic fields.

58 GEOSCIENCES↗
DOE OSTI2022

Dynamic compression of water to conditions in ice giant interiors

Abstract Recent discoveries of water-rich Neptune-like exoplanets require a more detailed understanding of the phase diagram of H 2 O at pressure–temperature conditions relevant to their planetary interiors. The unusual non-dipolar magnetic fields of ice giant planets, produced by convecting liquid ionic water, are influenced by exotic high-pressure states of H 2 O—yet the structure of ice in this state is challenging to determine experimentally. Here we present X-ray diffraction evidence of a body-centered cubic (BCC) structured H 2 O ice at 200 GPa and ~ 5000 K, deemed ice XIX, using the X-ray Free Electron Laser of the Linac Coherent Light Source to probe the structure of the oxygen sub-lattice during dynamic compression. Although several cubic or orthorhombic structures have been predicted to be the stable structure at these conditions, we show this BCC ice phase is stable to multi-Mbar pressures and temperatures near the melt boundary. This suggests variable and increased electrical conductivity to greater depths in ice giant planets that may promote the generation of multipolar magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗
DOE OSTI2022

Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field

Protein–ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)–a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein–ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ~0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

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

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL–water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO 3 ] – anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV E , and liquid structures g(r)⁠. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid–water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid–water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ~0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. Furthermore, the ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Early Solar Wind and Dynamo Magnetic Field Topology Predictions for (16) Psyche and Other Asteroids

Abstract Asteroid (16) Psyche is a metal‐rich body that might record an ancient coherent magnetization if some relict crust or mantle is preserved. Herein, we use magnetohydrodynamic simulations to predict (16) Psyche's field topology for several distinct pathways: (i) an early solar wind‐induced magnetization imparted after a larger body was impacted, forming the present‐day asteroid, (ii) a core dynamo magnetization imparted in an asteroid that is either presently largely intact or was a rubble pile, and (iii) magnetization in the turbulent solar nebula disk. For pathway (i) we find the field to be predominantly dipolar and spin axis‐aligned. For pathway (ii) we find the field to be either dipolar and spin axis‐misaligned, or highly multipolar. We also find that (iii) a field produced earlier before the solar nebula cleared, would be highly multipolar. In cases (i) and (ii) we also place constraints on the field strength. Simple detection of a magnetic field without constraining its topology and temporal variability would be insufficient to confirm a remanent source, due to the influence of the present‐day solar wind, electromagnetic induction, and (16) Psyche's high obliquity. For sufficiently strong fields however, the field topology and orientation may reveal key observable consequences of the nature and history of (16) Psyche. Our framework is also broadly applicable to the study of magnetic fields from other asteroids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Experimental investigation of microwave interaction with magnetoplasma in miniature multipolar configuration using impedance measurements

A miniature microwave plasma source employing both radial and axial magnetic fields for plasma confinement has been developed for micro-propulsion applications. Plasma is initiated by launching microwaves via a short monopole antenna to circumvent geometrical cutoff limitations. The amplitude and phase of the forward and reflected microwave power is measured to obtain the complex reflection coefficient from which the equivalent impedance of the plasma source is determined. Effect of critical plasma density condition is reflected in the measurements and provides insight into the working of the miniature plasma source. A basic impedance calculation model is developed to help in understanding the experimental observations. From experiment and theory, it is seen that the equivalent impedance magnitude is controlled by the coaxial discharge boundary conditions, and the phase is influenced primarily by the plasma immersed antenna impedance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Beyond the electric-dipole approximation: A formulation and implementation of molecular response theory for the description of absorption of electromagnetic field radiation

We present a formulation of molecular response theory for the description of a quantum mechanical molecular system in the presence of a weak, monochromatic, linearly polarized electromagnetic field without introducing truncated multipolar expansions. The presentation focuses on a description of linear absorption by adopting the energy-loss approach in combination with the complex polarization propagator formulation of response theory. Going beyond the electric-dipole approximation is essential whenever studying electric-dipole-forbidden transitions, and in general, non-dipolar effects become increasingly important when addressing spectroscopies involving higher-energy photons. These two aspects are examined by our study of the near K-edge X-ray absorption fine structure of the alkaline earth metals (Mg, Ca, Sr, Ba, and Ra) as well as the trans-polyenes. In following the series of alkaline earth metals, the sizes of non-dipolar effects are probed with respect to increasing photon energies and a detailed assessment of results is made in terms of studying the pertinent transition electron densities and in particular their spatial extension in comparison with the photon wavelength. Along the series of trans-polyenes, the sizes of non-dipolar effects are probed for X-ray spectroscopies on organic molecules with respect to the spatial extension of the chromophore.

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

Multipolar Skyrmion Crystals in Non-Kramers Doublet Systems

Here, we study the Kondo lattice model of multipolar magnetic moments interacting with conduction electrons on a triangular lattice. Bond-dependent electron hoppings induce a compasslike anisotropy in the effective Ruderman-Kittel-Kasuya-Yosida interaction between multipolar moments. This unique anisotropy stabilizes multipolar skyrmion crystals at zero magnetic field. In a unit cell, the skyrmion fractionalizes into meron composites subject to the conservation of total topological charge. Diverse multipolar phases in the phase diagram give rise to novel spontaneous Hall response of conduction electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗