Search NASA⌕ Search

SEARCH · Search NASA

Results for “Moments”

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 55 records · Page 3

Quantum closures for neutrino moment transport

A computationally efficient method for calculating the transport of neutrino flavor in simulations is to use angular moments of the neutrino one-body reduced density matrix, i.e., “quantum moments.” As with any moment-based radiation transport method, a closure is needed if the infinite tower of moment evolution equations is truncated. We derive a general parametrization of a quantum closure and the limits the parameters must satisfy in order for the closure to be physical. We then derive from multiangle calculations the evolution of the closure parameters in two test cases which we then progressively insert into a moment evolution code and show how the parameters affect the moment results until the full multiangle results are reproduced. This parametrization paves the way to setting prescriptions for genuine quantum closures adapted to neutrino transport in a range of situations.

79 ASTRONOMY AND ASTROPHYSICS↗

Impurity-induced moment freezing in NaFe x⁢ Ru 1–x O 2

We report the impact of magnetic impurity substitution on the quantum disordered magnetic ground state of NaRuO 2 . Local S = 5/2 moments are introduced into the frustrated triangular lattice of J eff = 1/2 Ru moments via Fe substitution in NaFe x⁢ Ru 1–x O 2 , and the evolution of the magnetic ground state is reported. Local spin freezing associated with conventional spin-glass behavior is observed upon Fe substitution, marking an impurity-induced freezing of the primarily dynamic magnetic ground state in NaRuO 2 . Furthermore, local Fe moments induce a Curie-Weiss magnetic behavior in the uniform magnetic susceptibility, and the local moment magnitude is best described by dynamic Ru moments polarized about impurity sites. Furthermore, our results establish an impurity-doping phenomenology consistent with inherently dynamic moments in NaRuO 2 that are pinned by local magnetic impurities, similar to “swiss cheese” models of impurity-substituted copper oxides.

36 MATERIALS SCIENCE↗

Bayesian Inference for the Seismic Moment Tensor Using Regional Waveforms and Teleseismic- P Polarities with a Data-Derived Distribution of Velocity Models and Source Locations

The largest source of uncertainty in any source inversion is the velocity model used in the transfer function that relates observed ground motion to the seismic moment tensor. However, standard inverse procedure often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. Here, we incorporate this uncertainty into an estimation of the seismic moment tensor using a data-derived distribution of velocity models based on complementary geophysical data sets, including thickness constraints, velocity profiles, gravity data, surface-wave group velocities, and regional body-wave travel times. The data-derived distribution of velocity models is then used as a prior distribution of Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional and teleseismic-P waveforms. The use of multiple data sets is important for gaining resolution to different components of the moment tensor. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and interpreted in terms of the most probable source type.

58 GEOSCIENCES↗

Calculation of the First Moment of Energy Using D-T Reactivity Formalisms Under the Maxwell-Boltzmann Distribution (Pt. 1)

Nuclear fusion science is an example of a scientific field with a rich history of expert involvement and scientific publications, which together, form an expert-knowledge base. One example of this history is the utilization of published reaction rates from a variety of authors. Investigators for Deuterium-Tritium (D-T) ion fusion can choose from using frequently cited methods: the Bosch and Hale reactivity, thermonuclear reaction rates from Caughlan and Fowler, and the reactivity evaluation from Miley, Towner & Ivich which forms the basis of the Naval Research Lab (NRL) formulary. There are other choices available. Each of the reactivity formulations considered here, are based upon the Maxwell- Boltzmann velocity distribution for D-T fusion ion reactants. Numerical methods for computer codes simulating hot, energetic plasmas, include tabulations of the reactivity, and the first moment of energy. The purpose of this study is to answer the question: what is the first moment of energy, and how has it been formulated? This report is part of a series of reports by the authors on D-T fusion-reaction formalism. The present focus is on defining the mathematical relationship of the first moment of D-T fusion ion kinetic energy, $\langle$E$\rangle$, with the fusion crosssection, fusion reactivity and its derivative with respect to ion-temperature. Three variants of the first moment $\langle$E$\rangle$ are analytically developed and explored: 1) constant cross-section, 2) a normalized first moment, and 3) a particular function of the first moment. Much of the mathematical development is relegated to the Appendices, for a concise presentation of topics and results. A research component of this discussion is our comparison of some of the reactivity and its derivative results between evaluations (1) versus (2) and (3), which the authors have not found in any other publication to date. This investigation is useful for validation, verification (V&V) and uncertainty quantification (UQ) as it compares the work of several authors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Majorana versus Dirac Constraints on the Neutrino Dipole Moments

Massive neutrinos are guaranteed to have nonzero electromagnetic moments and, since there are at least three neutrino species, these dipole moments define a matrix. Here, we estimate the current upper bounds on all independent neutrino electromagnetic moments, concentrating on Earth-bound experiments and measurements with solar neutrinos, including the very recent results reported by XENONnT. We make no simplifying assumptions and compare the hypotheses that neutrinos are Majorana fermions or Dirac fermions. In particular, we fully explore constraints in the Dirac-neutrino parameter space. Majorana and Dirac neutrinos are different; for example, the upper bounds on the magnitudes of the elements of the dipole moment matrix are weaker for Dirac neutrinos, relative to Majorana neutrinos. The potential physics reach of next-generation experiments also depends on the nature of the neutrino. We find that a next-generation experiment two orders of magnitude more sensitive to the neutrino electromagnetic moments via $\nu_{\mu}$ elastic scattering may discover that the neutrino electromagnetic moments are nonzero if the neutrinos are Dirac fermions. Instead, if the neutrinos are Majorana fermions, such a discovery is ruled out by existing solar neutrino data, unless there are more than three light neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗

Dilute magnetic impurity-induced effective phonon magnetic moment in Fe-doped monolayer MoS 2

Realization of large effective phonon magnetic moment in monolayer MoS 2 has established an important route for exploring intriguing magnetic phenomena in a nonmagnetic material. The sizable coupling between the orbital transition and the circularly polarized phonon results in the large effective phonon magnetic moment. In this work, using magneto-Raman spectroscopy, we investigate substitutional doping of magnetic atoms as a tuning knob of the electronic and phononic properties of MoS 2 . We show that Fe-doping polarizes the spin of the conduction bands and introduces a localized Fe band underneath the conduction band. As a result, an additional orbital transition between the Mo 4d and Fe 3d states emerges, producing an orbital-phonon hybridized mode at 283 cm −1 . Our magnetic field dependent measurements demonstrate that this new mode carries 2.8 $\mu_{\mathrm{B}}$ effective phonon magnetic moment, which is comparable to that of the undoped MoS 2 . Moreover, even though a long-range magnetic order is absent in Fe-doped MoS 2 , the local magnetic moment of Fe modifies the nature of the spin fluctuation, producing monotonically increasing quasielastic scattering spectral weight as temperature decreases. Our results highlight two-dimensional dilute magnetic semiconductors synthesized by substitutional doping as a promising material platform to manipulate the phonon magnetic moment through orbital-phonon coupling.

36 MATERIALS SCIENCE↗

Tuning the melting point of selected ionic liquids through adjustment of the cation’s dipole moment.

In prior research with thermally robust salts [Cassity et al., PCCP, 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. Here, we also test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation–anion interactions, which are recovered upon melting. If, however, the dipole- dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimization of foreground moment deprojection for semi-blind CMB polarization reconstruction

Abstract Upcoming Cosmic Microwave Background (CMB) experiments, aimed at measuring primordial CMB polarization B-modes, require exquisite control of instrumental systematics and Galactic foreground contamination. Blind minimum-variance techniques, like the Needlet Internal Linear Combination (NILC), have proven effective in reconstructing the CMB polarization signal and mitigating foregrounds and systematics across diverse sky models without suffering from foreground mismodelling errors. Still, residual foreground contamination from NILC may bias the recovered CMB polarization at large angular scales when confronted with the most complex foreground scenarios.By adding constraints to NILC to deproject statistical moments of the Galactic emission, the Constrained Moment ILC (cMILC) method has been demonstrated to further enhance foreground subtraction, albeit with an associated increase in overall noise variance. Faced with this trade-off between foreground bias reduction and overall variance minimization, there is still no recipe on which moments to deproject and which are better suited for blind variance minimization. To address this, we introduce the optimized cMILC (ocMILC) pipeline, which performs full automated optimization of the required number and set of foreground moments to deproject, pivot parameter values, and deprojection coefficients across the sky and angular scales, depending on the actual sky complexity, available frequency coverage, and experiment sensitivity. The optimal number of moments for deprojection, before paying significant noise penalty, is determined through a data diagnosis inspired by the Generalized NILC (GNILC) method.Validated on B-mode simulations of thePICOspace mission concept with four challenging foreground models, ocMILC exhibits lower Galactic foreground contamination compared to NILC and cMILC at all angular scales, with limited noise penalty. This multi-layer optimization enables the ocMILC pipeline to achieve unbiased posteriors of the tensor-to-scalar ratio, regardless of foreground complexity.

Astronomy & Astrophysics↗

Dark Energy Survey Year 3 results: cosmology with moments of weak lensing mass maps – validation on simulations

ABSTRACT We present a simulated cosmology analysis using the second and third moments of the weak lensing mass (convergence) maps. The second moment, or variances, of the convergence as a function of smoothing scale contains information similar to standard shear two-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The analysis is geared towards the third year (Y3) data from the Dark Energy Survey (DES), but the methodology can be applied to other weak lensing data sets. We present the formalism for obtaining the convergence maps from the measured shear and for obtaining the second and third moments of these maps given partial sky coverage. We estimate the covariance matrix from a large suite of numerical simulations. We test our pipeline through a simulated likelihood analyses varying 5 cosmological parameters and 10 nuisance parameters and identify the scales where systematic or modelling uncertainties are not expected to affect the cosmological analysis. Our simulated likelihood analysis shows that the combination of second and third moments provides a 1.5 per cent constraint on S8 ≡ σ8(Ωm/0.3)0.5 for DES Year 3 data. This is 20 per cent better than an analysis using a simulated DES Y3 shear two-point statistics, owing to the non-Gaussian information captured by the inclusion of higher order statistics. This paper validates our methodology for constraining cosmology with DES Year 3 data, which will be presented in a subsequent paper.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of point spread function higher moments error on weak gravitational lensing

ABSTRACT Weak gravitational lensing is one of the most powerful tools for cosmology, while subject to challenges in quantifying subtle systematic biases. The point spread function (PSF) can cause biases in weak lensing shear inference when the PSF model does not match the true PSF that is convolved with the galaxy light profile. Although the effect of PSF size and shape errors – i.e. errors in second moments – is well studied, weak lensing systematics associated with errors in higher moments of the PSF model require further investigation. The goal of our study is to estimate their potential impact for LSST weak lensing analysis. We go beyond second moments of the PSF by using image simulations to relate multiplicative bias in shear to errors in the higher moments of the PSF model. We find that the current level of errors in higher moments of the PSF model in data from the Hyper Suprime-Cam survey can induce a ∼0.05 per cent shear bias, making this effect unimportant for ongoing surveys but relevant at the precision of upcoming surveys such as LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

First-principles electronic structure in second-moment calculation of mode frequencies: Failure of quasiharmonic approximation in silicon

Using first-principles electronic structure calculations (local density approximation), we present second-moment calculations of the temperature dependence of the vibrational normal mode frequencies in Si up to 1500 K . The method is based on simple ensemble averages of displacements and forces and is easy to implement. These efficient calculations incorporate all interactions nonperturbatively, implicitly including all high order force constants, in contrast to typical perturbative approaches. Also in this work, we propose and apply to our classical calculations an estimated quantum correction that could be used for classical approaches other than the moments method. We compare our results to available experiments, including recent inelastic neutron scattering data and conclude that the second-moment gives the frequency shift to within a factor of 2 in contrast to quasiharmonic approaches that fail by an order of magnitude to predict the anharmonic shift. Noting the relative simplicity and speed of the calculations, we conclude that the second-moment approximation reasonably accounts for anharmonicity in a system that has been very closely studied experimentally and may reliably be applied in an exploratory way to more complex systems. Finally, we suggest that further improvement on the results could be obtained by incorporating the fourth moment, which can be implemented with little additional cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spontaneous Hall effect enhanced by local Ir moments in epitaxial Pr 2 Ir 2 O 7 thin films

Rare-earth pyrochlore iridates ( RE 2 Ir 2 O 7 ) consist of two interpenetrating cation sublattices, the RE with highly frustrated magnetic moments, and the iridium with extended conduction orbitals significantly mixed by spin-orbit interactions. The coexistence and coupling of these two sublattices create a landscape for discovery and manipulation of quantum phenomena such as the topological Hall effect, massless conduction bands, and quantum criticality. Thin films allow extended control of the material system via symmetry-lowering effects such as strain. While bulk Pr 2 Ir 2 O 7 shows a spontaneous hysteretic Hall effect below 1.5 K, we observe the effect at elevated temperatures up to 15 K in epitaxial thin films on (111) yttria-stabilized zirconia (YSZ) substrates synthesized via solid-phase epitaxy. Similar to the bulk, the lack of observable long-range magnetic order in the thin films points to a topological origin. We use synchrotron-based element-specific x-ray diffraction and x-ray magnetic circular dichroism to compare powders and thin films to attribute the spontaneous Hall effect in the films to localization of the Ir moments. Further, we link the thin-film Ir local moments to lattice distortions absent in the bulklike powders. We conclude that the elevated-temperature spontaneous Hall effect is caused by the topological effect originating either from the Ir or Pr sublattice, with interaction strength enhanced by the Ir local moments. This spontaneous Hall effect with weak net moment highlights the effect of vanishingly small lattice distortions as a means to discover topological phenomena in metallic frustrated magnetic materials.

36 MATERIALS SCIENCE↗

New method for detecting fast neutrino flavor conversions in core-collapse supernova models with two-moment neutrino transport

Fast-pairwise neutrino oscillations potentially affect many aspects of core-collapse supernovae (CCSNe): the explosion mechanism, neutrino signals, and nucleosynthesis in the ejecta. This particular mode of collective neutrino oscillations has a deep connection to the angular structure of neutrinos in momentum space; for instance, the appearance of electron-neutrino lepton number (ELN) angular crossings in momentum space is a good indicator of the occurrence of flavor conversions. However, many multidimensional CCSN simulations are carried out with approximate neutrino transport (such as two-moment methods), which limits access to the angular distributions of neutrinos, i.e., inhibits ELNcrossing searches. In this paper, we develop a new method of searching for ELN crossing in these CCSN simulations. The required data is the zeroth and first angular moments of neutrinos and the matter profile, all of which are available in CCSN models with the two-moment method. One of the novelties of our new method is the use of a ray-tracing neutrino transport to determine ELNs in the direction of the stellar center. It is designed to compensate for shortcomings of the crossing searches with only two angular moments. We assess the capability of the method by carrying out a detailed comparison with results of full Boltzmann neutrino transport in 1D and 2D CCSN models. In this work, we find that the ray-tracing neutrino transport improves the accuracy of crossing searches; indeed, the appearance/disappearance of the crossings is accurately detected even in the region of forward-peaked angular distributions. The new method is computationally cheap and has the benefit of efficient parallelization; hence, it will be useful for ELN-crossing searches in any CCSN models that employ two-moment neutrino transport.

79 ASTRONOMY AND ASTROPHYSICS↗

Gate-Tunable Phonon Magnetic Moment in Bilayer Graphene

Here we develop a first-principles quantum scheme to calculate the phonon magnetic moment in solids. As a showcase example, we apply our method to study gated bilayer graphene, a material with strong covalent bonds. According to the classical theory based on the Born effective charge, the phonon magnetic moment in this system should vanish, yet our quantum mechanical calculations find significant phonon magnetic moments. Furthermore, the magnetic moment is highly tunable by changing the gate voltage. Our results firmly establish the necessity of the quantum mechanical treatment, and identify small-gap covalent materials as a promising platform for studying tunable phonon magnetic moment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗

Propagation Pattern for Moment Representation of the Lattice Boltzmann Method

A propagation pattern for the moment representation of the regularized lattice Boltzmann method (LBM) in three dimensions is presented. Using effectively lossless compression, the simulation state is stored as a set of moments of the lattice Boltzmann distribution function, instead of the distribution function itself. An efficient cache-aware propagation pattern for this moment representation has the effect of substantially reducing both the storage and memory bandwidth required for LBM simulations. This article extends recent work with the moment representation by expanding the performance analysis on central processing unit (CPU) architectures, considering how boundary conditions are implemented, and demonstrating the effectiveness of the moment representation on a graphics processing unit (GPU) architecture.

42 ENGINEERING↗