SEARCH · Engineering Papers

Results for “second order phase transitions”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Theory for anomalous terahertz emission in striped cuprate superconductors

Recent experiments in the doped cuprates La 2−𝑥 ⁢Ba 𝑥⁢ CuO 4 have revealed the emission of anomalous terahertz radiation after impulsive optical excitation. Here, we theoretically investigate the nonlinear electrodynamics of such striped superconductors and explore the origin of the observed radiation. We argue that photoexcitation is converted into a photocurrent by a second-order optical nonlinearity, which is activated by the breaking of inversion symmetry in certain stripe configurations. We point out the importance of including umklapp photocurrents modulated at the stripe periodicity itself, which impulsively drive surface Josephson plasmons and lead to a resonant structure of outgoing radiation, consistent with the experiments. We speculate on the utility of the proposed mechanism in the context of generating tunable terahertz radiation.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Intermediate-Term Thermal Aging Effect Evaluation for Grade 92 and 316L at The LWR Relevant Temperature

Life extension of the existing nuclear reactors imposes accumulated damages, such as higher fluences and longer period of corrosion, to structural materials, which would result in significant challenges to the traditional reactor materials such as type 304 and 316 stainless steels. Advanced alloys with superior radiation resistance will increase safety margins, design flexibility, and economics for not only the life extension of the existing fleet but also new builds with advanced reactor designs. The Electric Power Research Institute (EPRI) teamed up with Department of Energy (DOE) on the Advanced Radiation Resistant Materials (ARRM) program, aiming to develop and test degradation resistant alloys from current commercial alloy specifications by 2021 to a new advanced alloy with superior degradation resistance in light water reactor (LWR)-relevant environments by 2024. Based on a comprehensive microstructure and property screening performed in Phase-1 of the ARRM program, a total of five alloys, together with 316L and X-750 as references, were down-selected for Phase-2 neutron irradiation studies. Because thermal aging could exert a synergistic effect on neutron irradiation due to the low neutron damage rate on the order of 10 –7 displacements per atom per second (dpa/s), Grade 92 (one of the five down-selected alloys) and two heats of 316L were selected in this task to study the effect of aging at 350°C for 12.5–12.7 kh on microstructure and mechanical properties. The aging time is approximately corresponding to 5 dpa neutron irradiation, which can be served as a reference for the 5-dpa-irradiated samples of the alloys to help understand the independent neutron damage influence on microstructure and mechanical properties. Optical microscopy and scanning electron microscopy were used for microstructural characterization. Hardness, tensile, Charpy impact toughness, and fracture toughness in the ductile regime of the aged samples were examined. The aging of Grade 92 led to the formation of many Laves phase in sizes of ~100–200 nm and resulted in some reduction in hardness and yield/ultimate tensile strength with some increases in uniform and total plastic elongations, which may have helped the enhancement of Charpy impact toughness, e.g., ~4 J increase in upper-shelf energy and 20.5°C reduction in ductile-brittle transition temperature compared with the unaged condition. The fracture toughness of the aged Grade 92 showed decent toughness of ~303 MPa√m (K Jq ) with ~87 tearing modulus at 22°C, which reduced to ~241 MPa√m (~20% reduction) with ~75 tearing modulus at 300°C. The two heats of 316L are differentiated by their amounts of δ-ferrite, i.e., ~1 vol% in heat T1103 and ~4 vol% in heat N5B8. The two heats of 316L were subjected to 15% cold work (CW) as this condition is often used in nuclear reactors to help trapping radiation-induced defects and thus postpone the steady state swelling stage. The aging did not result in any noticeable microstructure changes, except for possible segregations at grain boundaries that need to be further investigated. The aging resulted in slight reduction in hardness of 316L-T1103 and slight increase in hardness of 316L-N5B8, which is consistent with their tensile testing results. The 15%CW might have introduced some inhomogeneity, resulting in large standard deviations in hardness, which were reduced after the aging. The aged 316L followed the same trend as the unaged condition for the strength and elongation results from the tensile testing. Minimum elongations appeared within ~300–600°C, above which the elongations seem to increase. Unlike the negligible or minor changes in microstructure, hardness and tensile results of the two heats of 316L, their fracture toughness showed noticeable difference. The aged 316L-T1103 showed good toughness of ~350 MPa√m with ~80 tearing modulus at 22°C, which reduced to ~271 MPa√m (~23% reduction) with ~83 tearing modulus at 300°C. In contrast, the aged 316L-N5B8 showed lower toughness of ~244 MPa√m with ~53 tearing modulus at 22°C, which reduced to ~180 MPa√m (~26% reduction) with ~31 tearing modulus at 300°C. The preliminary results indicate that the presence of high volume of δ-ferrite would noticeable impair fracture toughness although it may not influence hardness and tensile properties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Constraining the hadronic spectrum and repulsive interactions in a hadron resonance gas via fluctuations of conserved charges

We simultaneously incorporate two common extensions of the hadron resonance gas model, namely the addition of extra, unconfirmed resonances to the particle list and the excluded volume repulsive interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow to constrain them separately. In particular, ratios of second-order susceptibilities like $\chi_{11}^{BQ}/\chi_2^B$ and $\chi_{11}^{BS}/\chi_2^B$ are sensitive only to the baryon spectrum, while fourth-to-second order ratios like $\chi_4^B/\chi_2^B$, $\chi_{31}^{BS}/\chi_{11}^{BS}$, or $\chi_{31}^{BQ}/\chi_{11}^{BQ}$ are mainly determined by repulsive interactions. Analysis of the available lattice results suggests the presence of both the extra states in the baryon-strangeness sector and the repulsive baryonic interaction, with indications that hyperons have a smaller repulsive core than non-strange baryons. Furthermore, the modified hadron resonance gas model presented here significantly improves the description of lattice QCD susceptibilities at chemical freeze-out and can be used for the analysis of event-by-event fluctuations in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Updated design of the CMB polarization experiment satellite LiteBIRD

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Unconventional fieldlike spin torques in CrPt 3

The topological semimetal CrPt 3 has potential for generating unconventional spin torques due to its ferrimagnetic ordering, topological band structure, and high anomalous Hall effect. CrPt 3 exhibits ferrimagnetic behavior only in its chemically ordered phase and is paramagnetic in its chemically disordered phase. By controlling the growth and annealing temperatures, we prepare epitaxial films of both chemically ordered and chemically disordered phases of CrPt 3 , allowing us to investigate the effect of magnetic ordering on unconventional-torque generation. We use angle-dependent spin-torque-ferromagnetic-resonance and second-harmonic Hall measurements to probe the spin torques generated from epitaxial CrPt 3 in CrPt 3 /Cu/Ni 81 ⁢Fe 19 heterostructures. With current applied along specific directions with respect to the crystal order, we reveal unconventional spin torques in both ordered and disordered films. When current flows parallel to the [$1\bar{1}1$] and [ $\bar{1}11$] directions, we observe an unconventional fieldlike torque that is opposite in sign for the two directions. Finally, our calculations reveal that this unconventional torque originates from an indirect nonlocal spin-orbit torque due to spin scattering at the CrPt 3 /Cu interface, in addition to symmetry breaking at this interface.

ferromagnetic resonance↗

DOE OSTI2022

Fluid-Kinetic Coupling: Advanced Discretizations for Simulations on Emerging Heterogeneous Architectures (LDRD FY20-0643)

Plasma physics simulations are vital for a host of Sandia mission concerns, for fundamental science, and for clean energy in the form of fusion power. Sandia's most mature plasma physics simulation capabilities come in the form of particle-in-cell (PIC) models and magnetohydrodynamics (MHD) models. MHD models for a plasma work well in denser plasma regimes when there is enough material that the plasma approximates a fluid. PIC models, on the other hand, work well in lower-density regimes, in which there is not too much to simulate; error in PIC scales as the square root of the number of particles, making high-accuracy simulations expensive. Real-world applications, however, almost always involve a transition region between the high-density regimes where MHD is appropriate, and the low-density regimes for PIC. In such a transition region, a direct discretization of Vlasov is appropriate. Such discretizations come with their own computational costs, however; the phase-space mesh for Vlasov can involve up to six dimensions (seven if time is included), and to apply appropriate homogeneous boundary conditions in velocity space requires meshing a substantial padding region to ensure that the distribution remains sufficiently close to zero at the velocity boundaries. Moreover, for collisional plasmas, the right-hand side of the Vlasov equation is a collision operator, which is non-local in velocity space, and which may dominate the cost of the Vlasov solver. The present LDRD project endeavors to develop modern, foundational tools for the development of continuum-kinetic Vlasov solvers, using the discontinuous Petrov-Galerkin (DPG) methodology, for discretization of Vlasov, and machine-learning (ML) models to enable efficient evaluation of collision operators. DPG affords several key advantages. First, it has a built-in, robust error indicator, allowing us to adapt the mesh in a very natural way, enabling a coarse velocity-space mesh near the homogeneous boundaries, and a fine mesh where the solution has fine features. Second, it is an inherently high-order, high-intensity method, requiring extra local computations to determine so-called optimal test functions, which makes it particularly suited to modern hardware in which floating-point throughput is increasing at a faster rate than memory bandwidth. Finally, DPG is a residual-minimizing method, which enables high-accuracy computation: in typical cases, the method delivers something very close to the $L^2$ projection of the exact solution. Meanwhile, the ML-based collision model we adopt affords a cost structure that scales as the square root of a standard direct evaluation. Moreover, we design our model to conserve mass, momentum, and energy by construction, and our approach to training is highly flexible, in that it can incorporate not only synthetic data from direct-simulation Monte Carlo (DSMC) codes, but also experimental data. We have developed two DPG formulations for Vlasov-Poisson: a time-marching, backward-Euler discretization and a space-time discretization. We have conducted a number of numerical experiments to verify the approach in a 1D1V setting. In this report, we detail these formulations and experiments. We also summarize some new theoretical results developed as part of this project (published as papers previously): some new analysis of DPG for the convection-reaction problem (of which the Vlasov equation is an instance), a new exponential integrator for DPG, and some numerical exploration of various DPG-based time-marching approaches to the heat equation. As part of this work, we have contributed extensively to the Camellia open-source library; we also describe the new capabilities and their usage. We have also developed a well-documented methodology for single-species collision operators, which we applied to argon and demonstrated with numerical experiments. We summarize those results here, as well as describing at a high level a design extending the methodology to multi-species operators. We have released a new open-source library, MLC, under a BSD license; we include a summary of its capabilities as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Hydrolysis Reaction Pathways of Thorium Oxide Nanoclusters

Density functional theory benchmarked by correlated molecular orbital theory is used to develop a fundamental and predictive understanding of the interaction of thorium oxide nanoclusters with gas phase water to provide insight into nuclear-waste storage, production of thorium nuclear reactor fuels, and reprocessing of spent fuel. The structures of Th n O 2n (n = 3 – 6) clusters and their interactions with water have been studied at the B3LYP, MP2, and CCSD(T) levels. Hydrolysis is initiated by the formation of Lewis acid-base adducts, with relative H 2 O binding energies (physisorption) ranging from −15 kcal/mol to −22 kcal/mol. The initial H 2 O physisorption energy is ca. −21 kcal/mol regardless of the cluster size and is consistent with the experimentally obtained initial adsorption energy on a thorium dioxide surface. The physisorption enthalpies for additional water molecules can be affected by the presence of terminal groups OH groups generated by proton transfer to a Th-O near the site of adsorption. The hydrolysis products (chemisorption) form either bridging or terminal hydroxides. More exothermic hydrolysis steps were predicted for the formation of terminal hydroxides as compared to the formation of bridging hydroxides. Here, the calculated transition state barriers for transfer of protons from bound water complexes to form the chemisorption products are very low. Overall, water readily reacts with thorium oxide clusters preferring hydroxide products over hydrated complexes. First and second order fits were predicted for the combined physisorption and chemisorption energies for the hydrolysis of thorium oxide clusters. Finally, ionization energies and electron affinities were calculated as were HOMO-LUMO gaps to provide additional insights into the properties of the thorium oxide and hydroxide clusters.

Adsorption↗

DOE OSTI2024

Magnetotransport of Sm 2 Ir 2 O 7 across the pressure-induced quantum-critical phase boundary

Rare-earth pyrochlore iridates host two interlocking magnetic sublattices of corner-sharing tetrahedra and can harbour a unique combination of frustrated moments, exotic excitations and highly correlated electrons. They are also the first systems predicted to display both topological Weyl semimetal and axion insulator phases. We have measured the transport and magnetotransport properties of single-crystal Sm 2 Ir 2 O 7 up to and beyond the pressure-induced quantum critical point for all-in-all-out (AIAO) Ir order at p c = 63 kbar previously identified by resonant X-ray scattering and close to which Weyl semimetallic behavior has been previously predicted. Our findings overturn the accepted expectation that the suppression of AIAO order should lead to metallic conduction persisting down to zero temperature. Instead, the resistivity-minimum temperature, which tracks the decrease in the AIAO ordering temperature for pressures up to 30 kbar, begins to increase under further application of pressure, pointing to the presence of a second as-yet unidentified mechanism leading to non-metallic behavior. The magnetotransport does track the suppression of Ir magnetism, however, with a strong hysteresis observed only within the AIAO phase boundary, similar to that found for Ho 2 Ir 2 O 7 and attributed to plastic deformation of Ir domains. Around pc we find the emergence of a new type of electronic phase, characterized by a negative magnetoresistance with small hysteresis at the lowest temperatures, and hysteresis-free positive magnetoresistance above approximately 5K. The temperature dependence of our low-temperature transport data are found to be best described by a model consistent with a Weyl semimetal across the entire pressure range.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Terahertz-Driven Local Dipolar Correlation in a Quantum Paraelectric

Light-induced ferroelectricity in quantum paraelectrics is a new avenue of achieving dynamic stabilization of hidden orders in quantum materials. In this Letter, we explore the possibility of driving a transient ferroelectric phase in the quantum paraelectric KTaO 3 via intense terahertz excitation of the soft mode. We observe a long-lived relaxation in the terahertz-driven second harmonic generation (SHG) signal that lasts up to 20 ps at 10 K, which may be attributed to light-induced ferroelectricity. Through analyzing the terahertz-induced coherent soft-mode oscillation and finding its hardening with fluence well described by a single-well potential, we demonstrate that intense terahertz pulses up to 500 kV=cm cannot drive a global ferroelectric phase in KTaO 3 . Instead, we find the unusual long-lived relaxation of the SHG signal comes from a terahertz-driven moderate dipolar correlation between the defect-induced local polar structures. We discuss the impact of our findings on current investigations of the terahertz-induced ferroelectric phase in quantum paraelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Low-field magnetic anomalies in single crystals of the A-type square-lattice antiferromagnet EuGa 4

The body-centered-tetragonal antiferromagnet EuGa 4 was recently identified as a Weyl nodal-line semimetal that exhibits the topological Hall effect below its reported antiferromagnetic (AFM) ordering temperature T N = 15 – 16.5 K which we find to be T N = 16.4 ( 2 ) K. The Eu +2 ions are located at the corners and body centers of the unit cells. EuGa 4 exhibits A-type antiferromagnetic order below T N , where the Eu 2+ spin-7/2 moments are ferromagnetically aligned in the ab plane with the Eu moments in adjacent Eu planes along the c axis aligned antiferromagnetically. Low-field magnetization versus field M(H ab ) data at T = 2 K with the field aligned in the ab plane are reported that exhibit anomalous positive curvature up to a critical field H c1 at which a second-order transition occurs with H c1 ≈ 0.85 kOe for H∥ [1,1,0] and ≈ 4.8 kOe for H∥[1,0,0]. For larger fields, the linear behavior M ab = χ(T N ) H a b is followed until the critical field H c is reached at which all moments become aligned with the applied field. A theory is formulated for T = 0 K that fits the observed M(H ab ) behavior at T = 2 K well, where domains of A -type AFM order with fourfold rotational symmetry occur in the AFM state in zero field. The moments in the four domains reorient to become almost perpendicular to H ab at H c1 , followed by increasing canting of all moments toward the field with increasing field up to H c which is reported to be 71 kOe. A first-order transition in M(H ab ) at H ab = H c1 is predicted by the theory for T = 0 K when H ab is at a small angle from the [1,0,0] direction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

DOE OSTI2020

Time-Resolved Spectroscopy of Insulator-Metal Transitions: Exploring Low-Energy Dynamics in Strongly Correlated Systems. Final Report

Transition metal oxides (TMOs) have attracted much attention in recent years. The strong coupling between charge, spin, and orbital degrees of freedom creates novel functional properties. The research at the interplay between ferromagnetism, ferroelectricity and conductivity provides a promising way to the development of future spintronic devices. The interface characterization of multifunctional oxide heterostructures is crucial to understand deeply the magnetic response, spin injection and spin transport. In year one, we characterized optically the magnetic property with the method of magnetization-induced second-harmonic generation (MSHG), which is a sensitive technique to selectively probe the interface where the spatial inversion symmetry is broken. The interface of n-type BaTiO 3 (BTO)/La 0.67 Sr 0.33 MnO 3 (LSMO) forms a Schottky barrier and displays a good diode effect. By applying an external electric field across the sample, we observe the magnetoelectric (ME) coupling effect to modulate the interface magnetization. The voltage dependent magnetic contrast sharply vanishes at positive voltage (applied to LSMO contact) and reveals a ferromagnetic (FM)-to- antiferromagnetic (AFM) transition occurring at the interface, while the magneto-optic Kerr effect (MOKE) indicates that the bulk property is not affected. The strain mediated mechanism and ferroelectric (FE) polarization induced charge medicated mechanism are both discussed, but they do not play an important role in the observed magnetization transition. A new mechanism is proposed - the injected minority spins through strong Hund's interaction with the local magnetic moments at positive voltage weaken the double-exchange interaction of nearby Mn e g electrons, thereby reducing the ferromagnetic (FM) ordering. The dominant superexchange interaction of localized t 2g electrons finally leads to the AFM configuration. By reducing the electron carrier concentration of BTO, a positive voltage shift of the transition is observed. The oxygen richer sample with lower carrier concentration requires a higher voltage to increase the hole doping level of LSMO to induce the interfacial AFM phase. By inserting a thin (1nm) La 0.5 Ca 0.5 MnO 3 layer between BTO and LSMO, the interface magnetic response from LSMO is enhanced. In year two, a non-multiferroic heterostructure SrTiO 3 (STO)/ La 0.5 Ca 0.5 MnO 3 (LCMO)/La 0.67 Sr 0.33 MnO 3 (LSMO) was studied. In this system, the interface FM-AFM transition is observed at negative voltage. The STO layer may act as a hole-donating layer, which induces the hole transfer to LSMO making it overdoped, while the LCMO interlayer at the 0.5 hole doping level displays complicated CE type AFM phase. Thus the AFM phase is dominant at the interface and electron accumulation is required to induce the FM phase by applying a negative voltage. For oxygen rich sample, the magnetic MSHG loop reappears at high positive voltage (around +2.5 V). At high doping level (x > 0.5), LSMO exhibits A-type AFM order, consisting of AFM superexchange coupling between adjacent layers and intralayer FM double-exchange coupling. The less tensile strain of STO and LSMO due to a smaller lattice mismatch as compared to BTO would result in a weaker interlayer coupling . Thus, interfacial A-type AFM order can occur in STO/LCMO/LSMO heterostructure at high hole doping concentration. The study of ME coupling of multifunctional oxide heterostructures is definitely important to the development and performance of future spintronic devices. The results presented here regarding the interfacial magnetic transition controlled by an external electric field offer some new perspective to study the spin transport of magnetic tunnel junctions.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Ab Initio Study of Vibronic and Magnetic 5f-to-5f and Dipole-Allowed 5f-to-6d and Charge-Transfer Transitions in [UX 6 ] n – (X = Cl, Br; n = 1, 2)

The absorption spectra of the octahedral [UX 6 ] n– (X = Cl, Br; n = 1, 2) complexes in the near-infrared (NIR) and UV–vis spectral regions were studied theoretically, using a relativistic restricted active space second-order perturbation theory (RASPT2) wavefunction framework, with the spin–orbit (SO) coupling treated by state interaction, in conjunction with Kohn–Sham density functional theory calculations for determining the vibrational normal modes. The electric-dipole-allowed ligand-to-metal charge-transfer (LMCT) and 5f-to-6d transitions, and the electric-dipole-forbidden 5f-to-5f ligand field (LF) transitions, are thereby obtained within the same theoretical framework. For the 5f-to-5f LF transitions, the observed absorption intensity is mostly due to vibronic coupling with low-energy electric-dipole-allowed transitions, but in some cases, the magnetic dipole intensity of the purely electronic transition has comparable intensity to the vibronic transitions. Here, experimental LF spectra of 5f 2 open-shell systems have been reported decades back, but ab initio calculations of their vibronic intensity have not yet been reported in the literature. Although the LF spectra for the 5f 2 systems can be assigned in detail, based on the calculations, the spectra are very complex and the underlying electronic states are strongly multiconfigurational. Therefore, the usefulness of the LF spectra beyond serving as a “fingerprint” of the LF and the metal oxidation state appears to be limited.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Tailoring Cementitious Materials Towards Value-Added Use of Large CO 2 Volumes

Hydraulic cements with alternative chemistries were developed for large-volume and value-added use of carbon dioxide. The hydraulic cements were processed using the energy-efficient mechanochemical technique at room temperature and atmospheric pressure. Carbon dioxide was captured directly from combustion emissions during processing. For this purpose, mechanochemical processing of hydraulic cements was accomplished under a flow of combustion emissions prior to the release of emissions to the atmosphere. The process removed a significant fraction of carbon dioxide from combustion emissions. The hydraulic cements captured carbon dioxide at about 10% of their weight. Two hydraulic cement chemistries were developed, and their mechanochemical processing was successfully scaled-up. The characteristic feature of one chemistry was its relatively low (near-neutral) pH where the integration of carbon dioxide yielded clear value. The second cement chemistry was based on alkali activation of industrial wastes. This chemistry could make value-added use of carbon dioxide, but it had to be refined to control the pH drop caused by CO 2 integration. These two cements render binding effects upon hydration by forming a combination of stable carbonates and aluminosilicates or phosphates. Efforts to develop cement chemistries based solely on carbonates were not successful. The mechanochemical process was found to integrate carbon dioxide into the alternative cement chemistries in the form of disordered and metastable carbonates. During hydration reactions, the disordered/metastable carbonates are either transformed into stable carbonate phases with desired binding efforts, or carbonates get integrated into the primary inorganic binders (hydrates). When compared with mechanochemical processing in pure carbon dioxide, mechanochemical processing in combustion emissions produced hydraulic cements with improved engineering properties. The mechanochemical process was scaled-up, and its variables as well as the raw materials formulations were optimized for implementation at pilot scale. Scale-up was found to enhance the carbon capture potential and the engineering qualities of the resulting hydraulic cements. This was because scale-up raises the intensity of mechanical energy input to raw materials. Certain mechanochemical phenomena cannot be induced, irrespective of the cumulative mechanical energy input, unless the intensity of impact is raised above a minimum level that cannot be achieved in laboratory-scale implementation of the process. Due to this effect, the duration of the mechanochemical process as well as its energy demand could be reduced significantly (by an order of magnitude) upon transition from laboratory to pilot scale, with the hydraulic cements produced at pilot scale offering engineering properties that were superior to those realized in laboratory-scale mechanochemical processing. The hydraulic cements produced at pilot scale via mechanochemical processing under a flow of combustion emissions were thoroughly characterized. They were found to meet standard requirements for ‘General Use’ hydraulic cements. They were compatible with the industrial-scale concrete production and construction practices that are used with the currently prevalent hydraulic cements. The new hydraulic cements with integrated carbon dioxide were found to offer distinct advantages over the currently prevalent Portland cement in terms of net carbon footprint and energy content. The combined raw materials and energy costs of the new hydraulic cements are competitive, and major cost savings can be realized because of the simplified production process that significantly lowers the capital investment in cement production plants. Mechanochemical processing of the new hydraulic cements under a flow of combustion emissions can be implemented using some existing components of cement manufacturing plants; this facilitates adoption of the technology by the cement industry.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Depinning, melting, and sliding of generalized Wigner crystals in Moiré systems

We numerically examine the depinning, sliding, and melting of commensurate and incommensurate Wigner crystals on two-dimensional hexagonal periodic substrates near fillings of 1 / 3 , 1 / 2 , and 2 / 3 to model the dynamics of generalized Wigner crystals in moiré heterostructures. At low temperatures where thermal fluctuations are irrelevant, for commensurate fillings we find a strongly pinned state that depins into a sliding crystal, while at incommensurate fillings the depinning threshold is strongly reduced. For fillings below 1 / 2 , the depinning occurs in a two-step process. Above the first depinning threshold, there is an extended range of drives where the conduction occurs via the sliding of antikinks along the charge stripes, followed by a second threshold where all the charges begin to slide. At finite temperatures, the external driving reduces the effective melting temperature at commensurate fillings and enhances creep at incommensurate fillings. We show that depinning into different sliding states, such as moving fluids or moving crystals, produces nonlinear features in the transport curves. We also show that transport is asymmetric on either side of a commensurate filling due to the different dynamics for interstitials versus holes in the commensurate structure. A variety of sliding states should be accessible for generalized Wigner crystals at finite temperatures even for low drives. If experiments can realize sliding dynamics in moiré systems, it would open another class of commensurate and incommensurate phases for study. Additionally, the sliding dynamics and antikink flow could provide unique functionalities for charge transport in moiré heterostructures. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Noncollinear polar magnet Fe 2 (SeO 3 ) 3 (H 2 O) 3 with inequivalent Fe 3+ sites

The emergence of novel magnetic states becomes more likely when the inversion symmetry of the crystal field, relative to the center between two spins, is broken. We propose that placing magnetic spins in inequivalent sites in a polar lattice can promote a realization of nontrivial magnetic states and associated magnetic properties. To test our hypothesis, we study Fe 2 (SeO 3 )(H 2 O) 3 as a model system that displays two distinct Fe(1) and Fe(2) magnetic sites in a polar structure (R3c space group). At low fields μ 0 H≤ 0.06 T, the material undergoes an antiferromagnetic ordering with T N1 = 77 K and a second transition at T N2 ≈ 4 K. At μ 0 H≥ 0.06 T and 74 K ≤T≤ 76 K, a positive entropy change of ∼0.12 mJ mol −1 K −1 can be associated with a metamagnetic transition to possibly nontrivial spin states. At zero field, Fe(1) is nearly fully ordered at T≈ 25 K, while Fe(2) features magnetic frustration down to T = 4 K. The magnetic ground state, a result corroborated by single-crystal neutron diffraction and 57 Fe Mössbauer spectroscopy, is a noncollinear antiparallel arrangement of ferrimagnetic Fe(1)–Fe(2) dimers along the c-axis. The results demonstrate that placing distinct magnetic sites in a polar crystal lattice can enable a new pathway to modifying spin, orbital, and lattice degrees of freedom for unconventional magnetism.

Crystal lattices↗

DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗