Search NASA⌕ Search

SEARCH · Search NASA

Results for “Higher order calculations”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The solubility of La hydroxide and stability of La 3+ and La hydroxyl complexes at acidic to mildly acidic pH from 25 to 250 °C

The mobility of rare earth elements (REE) in natural hydrothermal systems can be assessed using geochemical modeling, which requires reliable thermodynamic data of relevant aqueous species. In this study, we evaluate the controls of pH and temperature on La speciation and the role of hydroxyl complexes in REE transport at hydrothermal conditions. Batch-type hydrothermal solubility experiments were conducted using synthetic La hydroxide powders equilibrated in perchloric acid-based aqueous solutions at temperatures between 150 and 250 °C and starting pH of 2 to 5. The La hydroxide solubility is retrograde with temperature and displays a strong pH dependence with a decrease in La concentrations from acidic to mildly acidic pH spanning between 3 and 5 orders of magnitude (e.g. log La molality of –2.5 to –7.2 at 250 °C). Thermodynamic optimizations using GEMSFITS allow to retrieve the standard partial molal Gibbs energies for the La 3+ aqua ion and the formation constants for the La hydroxyl species (i.e., LaOH 2+ , La(OH) 2 + , La(OH) 3 0 ) between 25 and 250 °C. A comparison between the experimentally derived thermodynamic properties with the calculated values from the Helgeson-Kirkham-Flowers equation of state parameters indicates an increased divergence with temperature. Discrepancies in standard partial molal Gibbs energies range between ~ 1 – 12 kJ/mol and result in a predicted La hydroxide solubility differing by up to 3 orders of magnitude at 250 °C. Speciation calculations indicate a higher stability of La 3+ and LaOH 2+ over the other La hydroxyl species in the studied pH range of 3.4 to 6. Here, the optimized thermodynamic properties for La aqueous species have important implications for modeling the solubility of REE minerals such as monazite and the mobility of REE in hydrothermal systems.

58 GEOSCIENCES↗
DOE OSTI2025

Progress in modeling hydrogen assisted ammonia oxidation with new experiments and a further reconciliation of the NH 3 + OH rate constant

Hydrogen-assisted oxidation of ammonia in a premixed, laminar flow tubular reactor under reducing conditions was investigated experimentally and through chemical kinetic modeling. Due to its impact on the competition for OH among ammonia and hydrogen, the rate constant for NH 3 + OH (R1) was determined through state-of-the-art theoretical kinetics calculations, employing composite energies that include the effects of higher order electronic excitations on the electronic energies along the variational reaction path and treating the limitations in the kinetics posed by the passage through a hydrogen-bonded complex. The resulting rate constant was in close agreement with the recent experimental value from Zaczek et al. (2025), settling a long-term dispute about the high-temperature value of k 1 and confirming within 20% the value previously used in modeling. The chemical kinetic model, with no other changes, captured well measured concentrations of NH 3 , H 2 , NO, and N 2 O from flow reactor oxidation of NH 3 /H 2 at slightly reducing conditions over a range of temperature (900-1350 K) and NH 3 /H 2 ratios (0.5-2.0). Comparison of the present results with reported data from a non-premixed setup indicates that for laminar flow tubular reactors, the reactor configuration may have implications for the observed H 2 consumption due to the possibility of preferential oxidation during mixing.

Ab initio theory↗
DOE OSTI2026

The Complexation Properties of Self-Defensive Microgel-Modified Antimicrobial Surfaces

The complexation of cationic antimicrobials with polyanionic microgels on a biomaterial surface can render that surface self-defensive against bacteria by killing those bacteria which physically contact the antimicrobial-loaded microgels. This killing has been attributed to the contact-driven transfer of antimicrobial from a microgel to a challenging bacterium, though much remains unknown about this process. Here, in this study, we use a combination of experiments and computational modeling to identify key aspects of the complexation phenomena which influence the self-defensive properties. We synthesize poly(acrylic acid) (PAA) microgels (∼2–5 μm diameter) via membrane emulsification and electrostatically deposit them onto polycaprolactone (PCL) coupons or onto glass to form a discontinuous submonolayer. Subsequent microgel loading with colistin or with Sub5 antimicrobial peptide (AMP) causes microgel deswelling. Under physiological conditions Sub5 remains stably sequestered whereas colistin is quickly released. Coarse-grained molecular dynamics (CGMD) simulations confirm stronger Sub5/PAA complexation. CGMD calculations also indicate that Sub5 forms dimers and higher-order structures, a prediction confirmed experimentally by Small-Angle X-ray Scattering (SAXS). Supramolecular structure entropically enhances the complexation strength because of enhanced counterion release per complexation event, and this finding can help identify other antimicrobials well suited for such a nonelutive yet self-defensive strategy. CGMD simulations also show that Sub5 has a higher complexation strength with the Staphylococcus aureus membrane than it does with PAA, confirming that there is a thermodynamic driving force for antimicrobial transfer. Such self-defensive surfaces significantly reduce S. aureus colonization (over 90% reduction relative to unmodified controls) in an in vitro hematogenous contamination model and remain cyto-compatible as evidenced by mesenchymal stem cell spreading and proliferation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantum Hamiltonian for the surface charge density on a ring torus and radiative decay of plasmons

Photon scattering and the ensuing excitation of surface plasmons in single rings, single composite rings, and many-ring systems was recently shown to provide dispersion and field distributions [K. V. Garapati et al., J. Phys. Commun. 2, 015031 (2018)] of potential for specific applications in particle and molecular trapping [R. Alaee et al., Appl. Phys. Lett. 109, 141102 (2016); M. Salhi et al., Phys. Rev. A 92, 033416 (2015)], in addition to metamaterials and sensing. Furthermore, following photon or electron interactions with metallic nanorings, both radiative and nonradiative decay channels are important in the consideration of the nanoparticle as a photon or phonon source, and in related applications. In this work, we quantize the electronic normal modes of the ring torus and obtain the radiative decay of plasmons. Due to a geometry-related complexity, we employ a perturbation approach to obtain analytical expressions for the radiative decay channel for a vacuum bounded single solid nanoring. In quantizing the fields, the frequency spectrum of the charge density normal modes of the nanoring is obtained and shown to agree with the exact quasistatic plasmon dispersion relations. Higher-order corrections beyond the presented zero- and first-order calculations may be obtained following the presented results. The results are of potential interest in quantum sensing such as entangling photons and plasmons, or plasmons and trapped molecules.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nonperturbative phonon scatterings and the two-channel thermal transport in Tl 3 VSe 4

We review the role of nonperturbative phonon scattering in strongly anharmonic materials having ultralow lattice thermal conductivity with unusual temperature dependence. We take Tl 3 VSe 4 as an example and investigate its lattice dynamics using perturbation theory (PT) up to the fourth order and molecular dynamics (MD) with a machine-learning potential. We find distinct differences of phonon linewidth between PT and MD in the whole Brillouin zone. The comparison between the theoretical phonon linewidths and experiments suggests that PT severely underestimates the phonon scatterings, even when the fourth-order anharmonicity is included. Moreover, we extend our calculations to higher temperatures and evaluate the two-channel thermal conductivity based on the unified theory developed by Simoncelli et al. [Nat. Phys. 15, 809 (2019)]. We find a crucial coherence contribution to the total thermal conductivity at high temperatures. Our results pave the path for future studies of phonon properties and lattice thermal conductivities of strongly anharmonic crystals beyond the conventional PT realm.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles in $\mathrm {PbPb}$ collisions at $\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 2.76$ and 5.02$\,\text {TeV}$

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective anisotropic flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in $\mathrm {PbPb}$ collisions at nucleon-nucleon center-of-mass energies $\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 2.76$ and 5.02$\,\text {TeV}$ with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Higher-order topology induced by structural buckling

Abstract Higher-order topological insulator (HOTI) states, such as two-dimension (2D) HOTI featured with topologically protected corner modes at the intersection of two gapped crystalline boundaries, have attracted much recent interest. However, the physical mechanism underlying the formation of HOTI states is not fully understood, which has hindered our fundamental understanding and discovery of HOTI materials. Here we propose a mechanistic approach to induce higher-order topological phases via structural buckling of 2D topological crystalline insulators (TCIs). While in-plane mirror symmetry is broken by structural buckling, which destroys the TCI state, the combination of mirror and rotation symmetry is preserved in the buckled system, which gives rise to the HOTI state. We demonstrate that this approach is generally applicable to various 2D lattices with different symmetries and buckling patterns, opening a horizon of possible materials to realize 2D HOTIs. The HOTIs so generated are also shown to be robust against buckling height fluctuation and in-plane displacement. A concrete example is given for the buckled $\beta $-Sb monolayer from first-principles calculations. Our finding not only enriches our fundamental understanding of higher-order topology, but also opens a new route to discovering HOTI materials.

42 ENGINEERING↗
DOE OSTI2023

Interplay between magnetism and band topology in the kagome magnets R Mn 6 Sn 6

Kagome-lattice magnets R Mn 6 Sn 6 recently emerged as a new platform to exploit the interplay between magnetism and topological electronic states. Some of the most exciting features of this family are the dramatic dependence of the easy magnetization direction on the rare-earth specie, despite other magnetic and electronic properties being essentially unchanged, and the kagome geometry of the Mn planes that in principle can generate flat bands and Dirac points; gapping of the Dirac points by spin-orbit coupling has been suggested recently to be responsible for the observed anomalous Hall response in the member TbMn 6 Sn 6 . In this paper, we address both issues with ab initio calculations. We have discovered the significant role played by higher-order crystal-field parameters and rare-earth magnetic anisotropy constants in these systems. Further, we demonstrate that the microscopic origin of rare-earth magnetic anisotropy can also be quantified and understood at various levels: ab initio, phenomenological, and analytical. In particular, using a simple and physically transparent analytical model based on perturbation theory, we are able to explain, with full quantitative agreement, the evolution of rare-earth magnetic anisotropy across the series. We analyze in detail the topological properties of Mn-dominated bands and demonstrate how they emerge from the multiorbital planar kagome model. We further show that, despite this fact, most of the topological features at the Brillouin zone corner K are strongly 3D and therefore cannot explain the observed quasi-2D anomalous Hall effect, while the most pronounced quasi-2D dispersion are too far removed from the Fermi level. By employing self-consistent calculations with ab initio many-body approaches, we demonstrate that the exchange-correlation effects beyond the density functional theory for itinerant Mn- d electrons do not significantly alter the obtained electronic and magnetic structure. Therefore, we conclude that, contrary to previous claims, the most pronounced 2D kagome-derived topological band features bear little relevance to transport in RMn 6 Sn 6 , albeit they may possibly be brought to focus by electron or hole doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Lattice disorder effect on magnetic ordering of iron arsenides

This study investigates magnetic ordering temperature in nano- and mesoscale structural features in an iron arsenide. Although magnetic ground states in quantum materials can be theoretically predicted from known crystal structures and chemical compositions, the ordering temperature is harder to pinpoint due to potential local lattice variations that calculations may not account for. In this work we find surprisingly that a locally disordered material can exhibit a significantly larger Néel temperature (T N ) than an ordered material of precisely the same chemical stoichiometry. Here, a EuFe 2 As 2 crystal, which is a ‘122’ parent of iron arsenide superconductors, is found through synthesis to have ordering below T N = 195 K (for the locally disordered crystal) or T N = 175 K (for the ordered crystal). In the higher T N crystals, there are shorter planar Fe-Fe bonds [2.7692(2) Å vs. 2.7745(3) Å], a randomized in-plane defect structure, and diffuse scattering along the [00 L] crystallographic direction that manifests as a rather broad specific heat peak. For the lower T N crystals, the a-lattice parameter is larger and the in-plane microscopic structure shows defect ordering along the antiphase boundaries, giving a larger T N and a higher superconducting temperature (T c ) upon the application of pressure. First-principles calculations find a strong interaction between c-axis strain and interlayer magnetic coupling, but little impact of planar strain on the magnetic order. Neutron single-crystal diffraction shows that the low-temperature magnetic phase transition due to localized Eu moments is not lattice or disorder sensitive, unlike the higher-temperature Fe sublattice ordering. This study demonstrates a higher magnetic ordering point arising from local disorder in 122.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stochastic Inflation at NNLO

Stochastic Inflation is an important framework for understanding the physics of de Sitter space and the phenomenology of inflation. In the leading approximation, this approach results in a Fokker-Planck equation that calculates the probability distribution for a light scalar field as a function of time. Despite its successes, the quantum field theoretic origins and the range of validity for this equation have remained elusive, and establishing a formalism to systematically incorporate higher order effects has been an area of active study. In this paper, we calculate the next-to-next-to-leading order (NNLO) corrections to Stochastic Inflation using Soft de Sitter Effective Theory (SdSET). In this effective description, Stochastic Inflation manifests as the renormalization group evolution of composite operators. The leading impact of non-Gaussian quantum fluctuations appears at NNLO, which is presented here for the first time; we derive the coefficient of this term from a two-loop anomalous dimension calculation within SdSET. We solve the resulting equation to determine the NNLO equilibrium distribution and the low-lying relaxation eigenvalues. In the process, we must match the UV theory onto SdSET at one-loop order, which provides a non-trivial confirmation that the separation into Wilson-coefficient corrections and contributions to initial conditions persists beyond tree level. Furthermore, these results illustrate how the naive factorization of time and momentum integrals in SdSET no longer holds in the presence of logarithmic divergences. It is these effects that ultimately give rise to the renormalization group flow that yields Stochastic Inflation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Evaluating Shallow Convection Parameterization Assumptions With a q t – w Quadrant Analysis

Abstract Uncertainties associated with the response of shallow clouds to global warming remain a great challenge for climate projection. Due to the small size of these clouds, parameterizations are required to represent them in both current‐ and next‐generation climate models. We present a quantitative evaluation of several important assumptions used in both mass‐flux and Assumed‐PDF Higher‐Order Closure parameterizations of shallow convection. We use large‐eddy simulations (LESs) of four different shallow convection regimes as benchmarks and apply a q t (total moisture)– w (vertical velocity) quadrant analysis technique to identify the “coherent structures” of moist and dry up/down‐drafts. The statistics of these coherent structures and the environment are then used to evaluate commonly used assumptions in mass‐flux parameterizations related to (a) the convective downdrafts and (b) the subplume variability of convective updrafts. For quantitative evaluation of the PDF closure used in AP‐HOCs, we perform offline calculations with the PDF closure used in the Cloud Layers Unified By Bi‐normals and the Simplified Higher‐Order Closure schemes by supplying the closure with statistical quantities directly calculated from the LESs. Then, the impact of parameterization assumptions embedded in the PDF closure in the same two categories is analyzed and understood through comparisons of the output from the closure with the statistics of the moist and dry up/down‐drafts from the LESs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems

Molten salt electrolysis can recover the actinides from spent nuclear fuels, and it involves a eutectic LiCl-KCl in molten form as an electrolyte. During reprocessing, the concentration of fission products such as La, Nd, Pr, etc., increases and affects the recovery efficiency of the electrolyte. In this work, thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems was carried out using the CALPHAD (Calculation of Phase Diagrams) approach for the first time. The thermodynamic functions for the pure salts were taken from the SGTE (Scientific Group Thermodata Europe) Substances (SSUB) database. The experimental thermochemical and phase equilibria data available in the literature were used as input for the assessment of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems. The model parameters for the KCl-LiCl system were adjusted to include the new Gibbs energy descriptions for the pure salts. In addition, the sublattice model for the liquid phase in the LiCl-PrCl 3 system was modified and reassessed to ensure the model compatibility for higher-order extrapolation. There is a good agreement between the experimental and calculated thermochemical and phase diagram data for all the optimized constituent binaries and the ternary system. Furthermore, this work will be beneficial for determining the solubility limit of PrCl 3 in molten LiCl-KCl electrolytes and their thermodynamic properties for improving the efficiency of the pyrochemical process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Time-interval distributions in a digital gamma multi-channel analyzer at extreme input rates

Here, this study explores the relationship between experimentally observed inter-event time-interval distributions (TIDs) and measured dead times for pulses generated by a HPGe gamma detector and processed by a digital spectrometer. The system utilizes a fast and slow channel, pile-up rejector, trapezoid filtering, and flash analog-to-digital converter. The experimentally derived TIDs were compared with theory for validation. The results demonstrate that the theoretical model reliably describes measured TIDs up to 40% dead time. However, significant distortion effects become increasingly pronounced at higher input rates. It appears that the deviation between the measured and calculated time-interval distributions take the shape of higher-order convolutions of the TID, which represent the time difference between counts for more than two successive events. In this work, theoretical functions for the TID are expanded to reproduce measurements up to 90%. This refinement in the interpretation and treatment of the measured TIDs provides improved accuracy and precision in the prediction of the true event rate and measured dead time in the counter. Even though it reflects occasional count loss that is specific to the set-up, it is expected that a similar adjustment may be applicable to other detection systems as well.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A precise measurement of the Z -boson double-differential transverse momentum and rapidity distributions in the full phase space of the decay leptons with the ATLAS experiment at $\sqrt{s}=8$ TeV

This paper presents for the first time a precise measurement of the production properties of the Z boson in the full phase space of the decay leptons. This is in contrast to the many previous precise unfolded measurements per formed in the fiducial phase space of the decay leptons. The measurement is obtained from proton–proton collision data collected by the ATLAS experiment in 2012 at $\sqrt{s}$ = 8 TeV at the LHC and corresponding to an integrated luminosity of 20.2 fb -1 . The results, based on a total of 15.3 million Z boson decays to electron and muon pairs, extend and improve a previous measurement of the full set of angular coefficients describing Z-boson decay. The double-differential cross section distributions in Z-boson transverse momentum pT and rapidity y are measured in the pole region, defined as 80 < m ℓℓ < 100 GeV, over the range |y| < 3.6. The total uncertainty of the normalised cross-section measurements in the peak region of the pT distribution is dominated by statistical uncertainties over the full range and increases as a function of rapidity from 0.5–1.0% for |y| < 2.0 to 2 - 7% at higher rapidities. The results for the rapidity-dependent transverse momentum distributions are compared to state of-the-art QCD predictions, which combine in the best cases approximate N 4 LL resummation with N 3 LO fixed-order per turbative calculations. The differential rapidity distributions integrated over pT are even more precise, with accuracies from 0.2–0.3% for |y| < 2.0 to 0.4–0.9% at higher rapidities, and are compared to fixed-order QCD predictions using the most recent parton distribution functions. The agreement between data and predictions is quite good in most cases.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Patterned dielectric back contact design for GaAs thermophotovoltaic devices

-Patterned-dielectric back contact structures in optoelectronic devices are designed to boost the reflectance of light from the device back surface while retaining a low-resistance pathway for electrical conductance. Their reduced light absorption at near- and sub-bandgap photon energies leads to improved luminescence in light-emitting diodes, greater photon recycling, voltage, and efficiency in photovoltaic cells, and greater recuperation of unabsorbed sub-bandgap light in thermophotovoltaic (TPV) systems. However, diffraction from the patterned features can deflect incident light in propagation directions that lead to light trapping and parasitic absorption in the cell. Here, in this article, we use rigorous coupled-wave analysis (RCWA) to study three-dimensional diffractive scattering of electromagnetic waves by periodic metal point-contact gratings on 1.42-eV GaAs TPV cells, to analyze their effect on unwanted sub-bandgap absorption in order to achieve higher TPV system efficiency. Solutions of Maxwell's equations calculated using RCWA are compared to measured sub-bandgap reflectance in experimental GaAs TPV devices with varying metal point-contact diameters and spacing. Modeling and experiments indicate decreased total reflectance due to these diffractive effects for a small point contact diameter of 1 μm, and this effect is much stronger at higher contact coverage fractions.

14 SOLAR ENERGY↗
DOE OSTI2020

Exciton energy transfer reveals spectral signatures of excited states in clusters

Electronic excitation and concomitant energy transfer leading to Penning ionization in argon–acetylene clusters generated in a supersonic expansion are investigated with synchrotron-based photoionization mass spectrometry and electronic structure calculations. Spectral features in the photoionization efficiency of the mixed argon–acetylene clusters reveal a blue shift from the 2 P 1/2 and 2 P 3/2 excited states of atomic argon. Analysis of this feature suggests that excited states of argon clusters transfer energy to acetylene, resulting in its ionization and successive evaporation of argon. Theoretically calculated Ar n (n = 2–6) cluster spectra are in excellent agreement with experimental observations, and provide insight into the structure and ionization dynamics of the clusters. A comparison between argon–acetylene and argon–water clusters reveals that argon solvates water better, allowing for higher-order excitons and Rydberg states to be populated. Furthermore, these results are explained by theoretical calculations of respective binding energies and structures.

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

Field redefinitions can be nonlocal

We revisit the lore establishing the allowed space of field redefinitions and show that there are essentially no restrictions. Our conclusions hold to all orders in perturbation theory and for any dispersion relation. Field redefinitions can be nonlocal, symmetry breaking, or in certain cases have explicit dependence on spacetime. We address field redefinitions that can be resummed into the propagator, which demonstrates how to perform perturbative calculations away from the minimum in field space. Field redefinitions are used to derive higher-order Schwinger-Dyson equations, which imply multiparticle soft theorems. Non-standard field redefinitions are showcased using both relativistic and nonrelativistic examples.

effective field theories↗
DOE OSTI2024

Measurements of Lund subjet multiplicities in 13 TeV proton-proton collisions with the ATLAS detector

This Letter presents a differential cross-section measurement of Lund subjet multiplicities, suitable for testing current and future parton shower Monte Carlo algorithms. This measurement is made in dijet events in 140 fb -1 of $\sqrt{s}$ =13 TeV proton–proton collision data collected with the ATLAS detector at CERN's Large Hadron Collider. The data are unfolded to account for acceptance and detector-related effects, and are then compared with several Monte Carlo models and to recent resummed analytical calculations. The experimental precision achieved in the measurement allows tests of higher-order effects in QCD predictions. Most predictions fail to accurately describe the measured data, particularly at large values of jet transverse momentum accessible at the Large Hadron Collider, indicating the measurement's utility as an input to future parton shower developments and other studies probing fundamental properties of QCD and the production of hadronic final states up to the TeV-scale.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗