Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduce order modeling”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Improving noisy free-energy measurements by adding more noise

Estimating free-energy differences using nonequilibrium work relations, such as the Jarzynski equality, is hindered by poor convergence when work fluctuations are large. For systems governed by overdamped Langevin dynamics, we propose the counterintuitive approach of adding noise in order to increase the precision of such calculations. Here, by introducing additional stochastic fluctuations to the system and rescaling its potential energy accordingly, we leave the thermodynamics of the system unchanged while increasing its relaxation rate. For a given time-dependent protocol this modification reduces the dissipated reduced work, leading to more accurate free-energy estimates. The method is designed to be used in experiment, and we illustrate its operation using computer simulations applied to two model systems. However, the regime of applicability of this strategy is likely limited, because it requires control of the system's potential energy in a way that is feasible in only a few experimental settings.

Whitelam, Stephen [Lawrence Berkeley National Labo↗
DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

INTEGRATED WORKFLOW MANAGEMENT FOR PARTICLE ACCELERATOR SIMULATION

Supercomputing systems are used for a wide range of computationally demanding tasks in many fields of science and engineering. They play a key role in numerical simulation, in which mathematical models are computed in order to simulate the behavior of physical systems. Scientists and engineers that use supercomputers for numerical simulation often have their productivity limited by the need to manually organize and manage extremely large amounts of data that are often produced and consumed by the software programs run on these systems. Recognizing these limitations, Kitware Inc. (Clifton Park, NY) and SLAC National Accelerator Laboratory (Menlo Park, CA) are developing an advanced software platform that can reduce the cognitive overhead required by knowledge workers when using supercomputers for numerical simulation. Phase I of the project is complete and includes the development of new capabilities for organizing simulation project files, improvements to the user interface and overall usability, and deployment of a “middle tier” server to sit between user desktop machines and supercomputers to offload much of the data management workload. The project also developed prototype software for executing sequences of numerical simulations, and a prototype for migrating supercomputing software to cloud-based computing systems to provide a potential alternative to supercomputers with different logistical and price-to-performance tradeoffs.

Tourtellott, John↗
DOE OSTI2020

TRUST Testing Procedure Documentation

The TRUST testbed was designed and fabricated to conduct vibration testing and validate the behavior of hyper-elastic foams in a pre-loaded condition. The main goals of this testbed design were to ensure symmetry, proper instrumentation alignment and mounting, simplified geometry to make the data analysis validation process more straightforward, and to apply and monitor the applied pre-load. Based off of the test parameters, a shaker system was purchased. The baseplate was designed to mount directly to the bolt holes provided on the top of the shaker table. It was determined that the specimens will be cylindrical and vary in thickness from 2-10mm. The baseplate riser, center-mass, and testbed cap were designed to be cylindrical to make alignment about the center axis more straightforward. It was determined to fabricate the parts previously mentioned using Aluminum 6061-T6 due to its stiffness in order to reduce the likelihood of the stiff material impacting the data being collected for the foam specimens. The design concept used to apply the pre-load to the testbed was a through bolt/nut configuration. A through hole was included in each piece of the testbed, specimens, and instrumentation large enough to allow for 1/8” clearance between the through bolt and each of these pieces. This was imperative to avoid causing friction and complicating the validation model. A load cell was then purchased and used to measure the pre-load being applied to the testbed. It was also required to collect data for 3 axes at 3 different locations 120° apart on the center-mass. This configuration was required to be placed on the cylindrical face and toward the top and bottom of the center-mass. Slots were milled out in the middle of each sensor block with tight tolerances for mounting each of the accelerometers at 120° apart. This same approach was taken to mount 2 accelerometers on the baseplate at 180° apart. Holes were milled out of the corners of each accelerometer mounting slot to relieve the corner and allow for a flat edge at the back of each slot. Sensor blocks were mounted into the center-mass using 6-32 bolts at the required locations. This initial testbed design was fabricated at the TA53 machine shop.

42 ENGINEERING↗
DOE OSTI2023

Robust highly durable solid oxide fuel cell cathodes – Improved materials compatibility & self-regulating surface chemistry

Solid oxide fuel cells (SOFCs) are electrochemical conversion devices that directly transform hydrogen or hydrocarbon fuels to electricity, with energy efficiencies as high as 90%, coupled with reduced emissions. Several factors, however, remain to be addressed when considering scale-up of SOFC technology, including the need to overcome decreased performance due to sluggish rates of the oxygen reduction reaction (ORR) at the cathode under reduced temperatures and susceptibility to degradation in performance from surface poisoning e.g. from chromia, while limiting the use of critical raw materials (lanthanides and transition metals) present in high performing mixed ionic electronic conducting electrodes like (La,Sr)CoO 3 (LSC). In this project we explored the key descriptors for determining ORR activity in SOFC electrodes and tried to recover performance degradation by applying them to SOFC electrodes. In order to do this, we first selected a model mixed ionic electronic conducting (MIEC) oxide, Pr-doped CeO 2 (Pr 0.1 Ce 0.9 O 2-δ , PCO), which is a chemically stable fluorite and free of inherent poison sources (e.g. Sr segregation in LSC) that potentially react with external impurities such as Cr-species vaporized from the interconnect. The three approaches originally planned in this project are as follows: 1) evaluation of scavenger exsolution characteristics, 2) study of scavengers gettering efficacy following Cr and Si poisoning and 3) integration of new compositions into porous electrodes. Among them, exceptional progress has been made in 2) and 3), especially understanding the role of surface infiltrants in impacting electrode performance and degradation of PCO materials. We found that the Smith acidity scale for binary oxides serves as a powerful descriptor for tuning and predicting the oxygen exchange kinetics on MIEC PCO surfaces. As a result, with infiltration with binary oxides, ranging from strongly basic (Li 2 O) to strongly acidic (SiO 2 ) onto the surface of porous PCO, it was possible to systematically vary the chemical surface exchange coefficient (k chem ) by 6 orders of magnitude! L i2 O increased k chem by nearly 1,000 times over that of pristine PCO, while SiO 2 decreased k chem by nearly the same factor. Strikingly, although the pre-exponential of k chem scales linearly with the acidity of the infiltrated binary oxide, there is nearly no change in the activation energy. With this insight, we attributed the origin of these dramatic changes in k chem values to the systematic increase and decrease in the surface electron density induced by infiltrated binary oxides. More interestingly, although both Cr 2 O 3 and SiO 2 were determined to be acidic by Smith, suggesting that this feature could likely be the primary reason that these compounds serve to poison the ORR on SOFC cathodes, the effect of poisoning could be subsequently tuned by adding multiple infiltrants and controlling their relative surface acidities. We also systematically examined the effect of serial infiltration of basic and acidic oxides. It turned out that serial infiltration of Li not only recovers approximately 20-fold degraded k chem of PCO by acidic Cr 2 O 3 but its k chem is enhanced even beyond that of the non-infiltrated PCO by more than three orders of magnitude. This was further verified with a screen-printing PCO symmetric cell in terms of the electrode performance (area-specific resistance, ASR) related to approach 3). These observations point to acidity as a key descriptor not only in tuning and predicting the ORR activity of SOFC cathodes that largely determines the overall performance of SOFC, but in mitigating and reactivating poisoned electrode performance. This work provides novel guidelines for making the electrode performance much more active and robust in SOFCs, which can further be applied to all applications requiring oxygen exchange reaction, such as electrolyzers, permeation membranes and gas sensors.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

The Installation of Direct Water-Cooling Systems to Reduce Cooling Energy Requirements for High-Performance Computing Centers

A large cluster of High-Performance Computing (HPC) equipment at the Lawrence Livermore National Laboratory in California was retrofitted with an Asetek cooling system. The Asetek system is a hybrid scheme with water-cooled cold plates on high-heat-producing components in the information technology (IT) equipment, and with the remainder of the heat being removed by conventional air-cooling systems. In order to determine energy savings of the Asetek system, data were gathered and analyzed two ways: using top-down statistical models, and bottom-up engineering models. The cluster, “Cabernet”, rejected its heat into a facilities cooling water loop which in turn rejected the heat into the same chilled water system serving the computer-room air handlers (CRAHs) that provided the air-based cooling for the room. Because the “before” and “after” cases both reject their heat into the chilled water system, the only savings is due to reduction in CRAH fan power. The top-down analysis showed a 4% overall energy savings for the data center (power usage effectiveness (PUE) —the ratio of total data center energy to IT energy— dropped from 1.60 to 1.53, lower is better); the bottom-up analysis showed a 3% overall energy savings (PUE from 1.70 to 1.66) and an 11% savings for the Cabernet system by itself (partial PUE of 1.51). Greater savings, on the order of 15-20%, would be possible if the chilled water system was not used for rejecting the heat from the Asetek system. About 37% of the heat from the Cab system was rejected to the cooling water, lower than at other installations.

Earni, Shankar↗
DOE OSTI2020

Arsenic stabilization performance of a novel starch-modified Fe-Mn binary oxide colloid

Arsenic (As) is an environmentally hazardous contaminant and can be a serious threat to human health. In China, remediation needs for a large number of As-contaminated sites renders a strong demand for efficient remedial reagents and cost effective approaches. In this study, a novel starch-modified Fe-Mn binary oxide (SFM), an amorphous colloidal material, has been synthesized as a remedial reagent and its As stabilization performance has been evaluated. A set of laboratory batch experiments were carried out by directly adding SFM of different dosages into three contaminated soils. Results demonstrated that SFM could transform As in soil from non-specifically and specifically sorbed fractions to the more stable form bounded to amorphous iron hydrous oxides, thus reducing the As concentration in TCLP leachates by 42.0-93.2%. Results from adsorption tests and microscopic analysis indicated that the interactions between SFM and As are mainly controlled by adsorption, oxidation, and precipitation processes. SFM has abundant surface hydroxyl groups, with excellent adsorption properties for both As(V) and As(III), with the maximum adsorption capacities of 160.63 and 284.64 mg/g respectively at pH 7.0. The adsorption process closely fitted pseudo second-order kinetics and Freundlich isotherm model. SFM could increase soil Eh and oxidize As(III) to As(V), which facilitated the As stabilization in soilColloidal iron-based material directly used for stabilization in As contaminated soils is reported here for the first time. Starch modification improves both the reactivity and mobility of SFM in soil. Our findings propose an efficient and convenient reagent for As remediation in soil.

Yan, Xiulan↗
DOE OSTI2025

Renewable Activated Carbon Sorbent for the Desulfurization of Liquid Fuels

In recent years, significant efforts have been devoted to the reduction of sulfur levels in transportation fuels in order to reach regulatory limits and reduce the emission of harmful SO x into the atmosphere. In particular, adsorptive desulfurization (ADS) has the potential to produce zero-sulfur fuels without the need for a high energy intensity process, high H 2 pressure, and long process times as in the case of other technologies such as hydrodesulfurization. In this work, we have demonstrated the effectiveness of renewable activated carbons derived from food waste (FWAC) for the ADS of model jet and diesel fuels. The optimal FWAC materials were those fabricated to maximize the micropore volume, providing available sites for the adsorption of dibenzothiophene (DBT) and dimethyldibenzothiophene (DMDBT). The FWAC had a greater sulfur adsorption capacity than commercial AC, Y zeolite sorbents, and AC derived from other biomass sources including miscanthus, coconut shell, and walnut shell. Elemental analysis suggests that the inorganic impurities inherent in food waste, notably K, Ca, P, and Na, may contribute to its improved sulfur adsorption compared to other AC materials. Microscopy and X-ray diffraction studies further demonstrated the presence of inorganic species on FWAC that may provide active sites for the chemisorption of sulfur molecules.

Adsorption↗
DOE OSTI2024

Using Molecular Structure to Tune Intrachain and Interchain Charge Transport in Indacenodithiophene-Based Copolymers

Here, in this work, we compare two structurally near-amorphous rigid-rod polymers-poly(indacenodithiophene-co-benzothiadiazole), p(IDT-BT), and poly(indacenodithiophene-co-benzopyrollodione), p(IDT-BPD)-with orders of magnitude different mobilities to understand the effect charge carrier intrachain delocalization has on electronic transport. Quantum chemical calculations show that p(IDT-BPD) has a barrier to torsion that is significantly lower than that of p(IDT-BT) and is thus more likely to have reduced conjugation lengths. We utilize absorption and photoluminescence spectroscopy to characterize energetic disorder and show that p(IDT-BPD) has higher energetic disorder. Charge modulation spectroscopy (CMS) and model calculations are used to show that charge carriers are substantially delocalized in p(IDT-BT) and occupy near-uniform energetic environments. We find that mobility activated hopping barriers are similar in these two materials. Electronic structure calculations show that both intrachain and interchain couplings of monomer units are poor enough in p(IDT-BPD) that charge carriers collapse to single IDT units and transport via a through-space tunneling mechanism. This work highlights the remarkable charge transport properties of p(IDT-BT) by showing that high mobilities are achievable on device-relevant length scales with only 1D carrier delocalization.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A semi-automated algorithm for designing stellarator divertor and limiter plates and application to HSX

We present a semi-automated algorithm for designing three-dimensional divertor or limiter plates targeting low heat loads. The algorithm designs the plates in two stages: firstly, the parallel heat flux distribution is caught on vertically-inclined plates at one or several toroidal locations. Secondly, the power per unit area is reduced by stretching, tilting and bending the plates toroidally. Heat transport is modelled using the EMC3-Lite code, which uses an anisotropic diffusion model. We apply this scheme to HSX, a medium-sized stellarator located at the University of Wisconsin–Madison. Starting from the current machine with an extended vessel wall, we construct plates which are able to effectively catch and spread the heat for three different magnetic configurations. The scheme has a computational cost in the order of tens of CPU-minutes, making it a powerful tool for semi-automated plasma-facing component design in three-dimensional environments.

anisotropic diffusion↗
DOE OSTI2020

Machine-learning prediction for quasiparton distribution function matrix elements

There have been rapid developments in the direct calculation in lattice QCD (LQCD) of the Bjorken-x dependence of hadron structure through large-momentum effective theory (LaMET). LaMET overcomes the previous limitation of LQCD to moments (that is, integrals over Bjorken x) of hadron structure, allowing LQCD to directly provide the kinematic regions where the experimental values are least known. LaMET requires large-momentum hadron states to minimize its systematics and allow us to reach small- x reliably. This means that very fine lattice spacing to minimize lattice artifacts at order $(P_za)^n$ will become crucial for next-generation LaMET-like structure calculations. Furthermore, such calculations require operators with long Wilson-link displacements, especially in finer lattice units, increasing the communication costs relative to that of the propagator inversion. In this work, we explore whether machine-learning algorithms can make predictions of correlators to reduce the computational cost of these LQCD calculations. We consider two algorithms, gradient-boosting decision tree and linear models, applied to LaMET data, the matrix elements needed to determine the kaon and $η_s$ unpolarized parton distribution functions (PDFs), meson distribution amplitude (DA), and the nucleon gluon PDF. We find that both algorithms can reliably predict the target observables with different prediction accuracy and systematic errors. The predictions from smaller displacement z to larger ones work better than those for momentum p due to the higher correlation among the data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Mitigating the noise of DESI mocks using analytic control variates

In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological N-body simulations provide a reference through which we can test our models, but their output suffers from sample variance on large scales. Fortunately, this is the regime where accurate analytic approximations exist. To reduce the variance, which is key to making optimal use of these simulations, we can leverage the accuracy and precision of such analytic descriptions using Control Variates (CV). The power of control variates stems from utilizing inexpensive but highly correlated surrogates of the statistics one wishes to measure. The stronger the correlation between the surrogate and the statistic of interest, the larger the variance reduction delivered by the method. We apply two control variate formulations to mock catalogs generated in anticipation of upcoming data from the Dark Energy Spectroscopic Instrument (DESI) to test the robustness of its analysis pipeline. Our CV-reduced measurements offer a factor of 5-10 improvement in the measurement error compared with the raw measurements. We explore the relevant properties of the galaxy samples that dictate this reduction and comment on the improvements we find on some of the derived quantities relevant to Baryon Acoustic Oscillation (BAO) analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A bi-level data-driven framework for fault-detection and diagnosis of HVAC systems

Long-term operation of heating, ventilation, and air conditioning (HVAC) systems will eventually lead to a range of HVAC system failures, resulting in excessive energy consumption and maintenance costs. Here, to avoid HVAC malfunctioning, fault detection diagnostic (FDD) is utilized as a common practice. Machine learning methods have lately received considerable interest for FDD analysis of HVAC systems due to their high detection accuracy. Meanwhile, HVAC malfunctions are regarded as rare occurrences, hence normal operating data samples are much more accessible than data samples in faulty and malfunctioning conditions. The dominating frequency of normal operation in HVAC datasets has also led to heavily biased classification algorithms within the literature. Moreover, the focus of previous literature has been on increasing the accuracy of the models which leads to a high number of false positives (misleading alarms) in the system. In order to enhance the performance of diagnostic procedures and fill the mentioned gaps, this study proposes a novel data-driven framework. A bi-level machine learning framework is developed for diagnosing faults in air handling units (AHUs) and rooftop units (RTUs) based on principal component analysis (PCA), time series anomaly detection, and random forest (RF). It is shown that PCA can reduce the dataset dimension with one principal component accounting for 95% of data variance. Also, the random forest could classify the faults with 89% precision for single-zone AHU, 85% precision for RTU, and 79% for multi-zone AHU. By proposing this framework, three persistent challenges are addressed: (I) minimizing false positives; (II) accounting for data imbalance; and (III) normal condition monitoring of equipment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Toward shell model interactions with credible uncertainties

Background: The nuclear shell model is a powerful framework for predicting nuclear structure observables, but relies on interaction matrix elements fit to experimental data as its inputs. Extending the shell model's applicability, particularly toward dripline nuclei, requires efficient fitting methods and credible uncertainty quantification. Traditional approaches face computational challenges and may underestimate uncertainties. Purpose: We develop and test a framework combining eigenvector continuation and Markov chain Monte Carlo to efficiently fit shell model interaction matrix elements and quantify their uncertainties. Methods: Eigenvector continuation is used to emulate shell model calculations, reducing computational costs. The emulator enables Markov chain Monte Carlo sampling to optimize interaction matrix elements and rigorously assess parametric uncertainties. Here, the framework is benchmarked using the USDB interaction in the 𝑠⁢𝑑 shell. Results: The emulator reproduces the USDB interaction with negligible error, validating its use in shell model fitting applications. However, we find that to obtain credible predictive intervals, the model defect of the shell model itself, rather than experimental or emulator error, must be taken into account in order to obtain credible uncertainties. Conclusions: The proposed framework provides an efficient and rigorous approach for fitting shell model interactions and quantifying uncertainties. Further, the normality assumption used in the past appears sufficient to describe the distribution of interaction matrix elements. However, it is crucial to account for model correlations to avoid underestimating uncertainties.

Nuclear forces↗
DOE OSTI2020

Evaluation of Leading Modes of Climate Variability in the CMIP Archives

The adequate simulation of internal climate variability is key for our understanding of climate as it underpins efforts to attribute historical events, predict on seasonal and decadal time scales, and isolate the effects of climate change. Here the skill of models in reproducing observed modes of climate variability is assessed, both across and within the CMIP3, CMIP5, and CMIP6 archives, in order to document model capabilities, progress across ensembles, and persisting biases. A focus is given to the well-observed tropical and extratropical modes that exhibit small intrinsic variability relative to model structural uncertainty. These include El Niño–Southern Oscillation (ENSO), the Pacific decadal oscillation (PDO), the North Atlantic Oscillation (NAO), and the northern and southern annular modes (NAM and SAM). Significant improvements are identified in models’ representation of many modes. Canonical biases, which involve both amplitudes and patterns, are generally reduced across model generations. For example, biases in ENSO-related equatorial Pacific sea surface temperature, which extend too far westward, and associated atmospheric teleconnections, which are too weak, are reduced. Stronger tropical expression of the PDO in successive CMIP generations has characterized their improvement, with some CMIP6 models generating patterns that lie within the range of observed estimates. For the NAO, NAM, and SAM, pattern correlations with observations are generally higher than for other modes and slight improvements are identified across successive model generations. Finally, for ENSO and PDO spectra and extratropical modes, changes are small compared to internal variability, precluding definitive statements regarding improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

On Ohm’s law in reduced plasma fluid models

Drift-reduced MHD models are widely used to study magnetised plasma phenomena, in particular for magnetically confined fusion applications, as well as in solar and astrophysical research. This letter discusses the choice of Ohm's law in these models, the resulting dispersion relations for the dynamics parallel to the magnetic field, and the implications for numerical simulations. We find that if electron pressure is included in Ohm's law, then both electromagnetic and finite electron mass effects must also be included in order to obtain physical dispersion relations. A simple modification to the plasma vorticity is also found which improves handling of low density regions, of particular relevance to the simulation of the boundary region of magnetised plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of the phase transition in random unitary circuits with measurements

Here, we present a theory of the entanglement transition tuned by measurement strength in qudit chains evolved by random unitary circuits and subject to either weak or random projective measurements. The transition can be understood as a nonanalytic change in the amount of information extracted by the measurements about the initial state of the system, quantified by the Fisher information. To compute the von Neumann entanglement entropy $\textit{S}$ and the Fisher information $\mathcal{F}$, we apply a replica method based on a sequence of quantities $\tilde{S}^{(n)}$ and $\mathcal{F}^{(n)}$ that depend on the $\textit{n}$th moments of density matrices and reduce to $\textit{S}$ and $\mathcal{F}$ in the limit $\textit{n}$ → 1. These quantities with $\textit{n}$ ≥ 2 are mapped to free energies of a classical spin model with $\textit{n}$! internal states in two dimensions with specific boundary conditions. In particular, $\tilde{S}^{(n)}$ is the excess free energy of a domain wall terminating on the top boundary, and $\mathcal{F}^{(n)}$ is related to the magnetization on the bottom boundary. Phase transitions occur as the spin models undergo ordering transitions in the bulk. Taking the limit of large local Hilbert space dimension $\textit{q}$ followed by the replica limit $\textit{n}$ → 1, we obtain the critical measurement probability $p_c$ = 1/2 and identify the transition as a bond percolation in the 2D square lattice in this limit. Finally, we show there is no phase transition if the measurements are allowed in an arbitrary nonlocal basis, thereby highlighting the relation between the phase transition and information scrambling. We establish an explicit connection between the entanglement phase transition and the purification dynamics of a mixed state evolution and discuss implications of our results to experimental observations of the transition and simulability of quantum dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Unraveling the magnetic ground state in the alkali-metal lanthanide oxide Na 2 ⁢Pr⁡O 3

Here, a comprehensive set of muon spin spectroscopy and neutron scattering measurements supported by ab initio and model Hamiltonian simulations have been used to investigate the magnetic ground state of Na 2 ⁢PrO 3 . μSR reveals a Néel antiferromagnetic order below T N ~ 4.9K, with a small static magnetic moment m static ≤ 0.22 μ B /Pr collinearly aligned along the c axis. Inelastic neutron measurements reveal the full spectrum of crystal field excitations and confirm that the Pr 4+ ground-state wave function deviates significantly from the Γ 7 limit that is relevant to the Kitaev model. Single- and two-magnon excitations are observed in the ordered state below T N =4.6K and are well described by nonlinear spin wave theory from the Néel state using a magnetic Hamiltonian with Heisenberg exchange J=1 meV and symmetric anisotropic exchange Γ/J=0.1, corresponding to an XY model. Intense two magnon excitations are accounted for by g-factor anisotropy g z /g ± = 1.29. A fluctuating moment δm 2 = 0.57 ⁢(22) ⁢μ$^2_B$/Pr extracted from the energy and momentum integrated inelastic neutron signal is reduced from expectations for a local J = 1/2 moment with average g factor g avg ≈ 1.1. Together, the results demonstrate that the small moment in Na 2 ⁢PrO 3 arises from crystal field and covalency effects and the material does not exhibit significant quantum fluctuations.

36 MATERIALS SCIENCE↗