Results for “Hysteresis model” — Page 6

DOE OSTI2023

Effective Landau-type model of a Hf x Zr 1 - x O 2 -graphene nanostructure

Here, to describe the charge-polarization coupling in the nanostructure formed by a thin Hf x Zr 1-x O 2 film with a single-layer graphene as a top electrode, we develop the “effective” Landau-Ginzburg-Devonshire model. This approach is based on the parametrization of the Landau expansion coefficients for the polar (FE) and antipolar (AFE) orderings in thin Hf 1-x Zr x O 2 films from a limited number of polarization-field curves and hysteresis loops. The Landau expansion coefficients are nonlinearly dependent on the film thickness h and Zr/[Hf+Zr] ratio x, in contrast to h-independent and linearly x-dependent expansion coefficients of a classical Landau energy. We explain the dependence of the Landau expansion coefficients by the strong nonmonotonic dependence of the Hf 1-x Zr x O 2 film polar properties on the film thickness, grain size and surface energy. The proposed Landau free energy with five “effective” expansion coefficients, which are interpolation functions of x and h, describes the continuous transformation of polarization dependences on applied electric field and hysteresis loop shapes induced by the changes of x and h in the range 0 < x < 1 and 5 nm < h < 35 nm. Using the effective free energy, we demonstrated that the polarization of Hf 1-x Zr x O 2 film influences strongly on the graphene conductivity, and the full correlation between the distribution of polarization and charge carriers in graphene is revealed. In accordance with our modeling, the polarization of the (5 – 25) nm thick Hf 1-x Zr x O 2 films, which are in the ferroelectric-like or antiferroelectric-like states for the chemical compositions 0.35 ≤ x ≤ 0.95, determine the concentration of carriers in graphene and can control its field dependence. The result can be promising for creation of next generation Si-compatible nonvolatile memories and graphene-ferroelectric FETs, because the working voltages applied to the Hf 1-x Zr x O 2 film (which acts as a gate) can be relatively low (less than 2 V). These low voltages are sufficient to induce the pronounced hysteresis of ferroelectric polarization in the Hf 1-x Zr x O 2 gate, which, due to the strong electric coupling, induces the hysteresis of the graphene charge.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2021

Propagation of priors for more accurate and efficient spectroscopic functional fits and their application to ferroelectric hysteresis

Multi-dimensional spectral-imaging is a mainstay of the scanning probe and electron microscopies, micro-Raman, and various forms of chemical imaging. In many cases, individual spectra can be fit to a specific functional form, with the model parameter maps, providing direct insight into material properties. Since spectra are often acquired across a spatial grid of points, spatially adjacent spectra are likely to be similar to one another; yet, this fact is almost never used when considering parameter estimation for functional fits. On datasets tried here, we show that by utilizing proximal information, whether it be in the spatial or spectral domains, it is possible to improve the reliability and increase the speed of such functional fits by ~2-3x, as compared to random priors. We explore and compare three distinct new methods: (1) spatially averaging neighborhood spectra, and propagating priors based on functional fits to the averaged case, (2) hierarchical clustering-based methods where spectra are grouped hierarchically based on response, with the priors propagated progressively down the hierarchy, and (3) regular clustering without hierarchical methods with priors propagated from fits to cluster means. Our results highlight that utilizing spatial and spectral neighborhood information is often critical for accurate parameter estimation in noisy environments, which we show for ferroelectric hysteresis loops acquired on a prototypical PbTiO3 thin film with piezoresponse spectroscopy. This method is general and applicable to any spatially measured spectra where functional forms are available. Examples include exploring the superconducting gap with tunneling spectroscopy, using the Dynes formula, or current-voltage curve fits in conductive atomic force microscopy mapping. Here we explore the problem for ferroelectric hysteresis, which, given its large parameter space, constitutes a more difficult task than, for example, fitting current-voltage curves with a Schottky emission formula.

42 ENGINEERING↗

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 OSTI2022

Analysis of Postulated Accident Scenarios for the Transformational Challenge Reactor

The Transformational Challenge Reactor (TCR) is a helium-cooled, yttrium-hydride-moderated reactor that was designed for the U.S. Department of Energy Office of Nuclear Energy. A key objective of the TCR was to employ advanced manufacturing techniques in a nuclear system and demonstrate their potential for revolutionizing the nuclear reactor design process. One purpose of the present work is to demonstrate the safety of the TCR under postulated accidents. Based on RELAP5-3D and COMSOL analyses, the TCR remained below all current safety limits and far below the expected failure limits for the core materials. Another purpose of this work is to provide useful insights and recommendations regarding the application of RELAP5-3D to gas-cooled or other advanced reactors. A novel approach was implemented for simultaneously modeling conduction and radiation in RELAP5-3D, which was found to provide reasonable predictions of radial core, vessel, and ex-vessel heat transfer during postulated events. A multicode approach was also applied, in which high-fidelity COMSOL calculations were used to tune the radial heat transfer parameters in RELAP5-3D. The tuned RELAP5-3D model demonstrated comparable peak temperature predictions as COMSOL, despite a coarse treatment of the core in RELAP5-3D consisting of only two lumped heat structures. This high-fidelity tuning approach enabled enhanced accuracy as well as minimal complexity within the RELAP5-3D model, even for complex fuel geometric designs as in the TCR. Finally, investigations were made into the potential for flow reversal during a pressurized loss-of-forced-flow event in the TCR. The TCR is designed with downward helium flow through the core during normal operation. The RELAP5-3D model predicted that this downward flow would persist, without flow reversal, up to several days after the circulator trip. This was attributed to natural circulation hysteresis effects as have been noted in similar thermofluidic systems. Although flow stagnation and eventual reversal did not lead to unsafe TCR conditions, interesting spatial effects were observed which may have safety relevance for other reactor system designs and coolant types that are designed for downward core flow during normal operation, warranting closer investigation of the flow reversal phenomenon.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Antiprotons and Elementary Particles over a Solar Cycle: Results from the Alpha Magnetic Spectrometer

We present results over an 11-year Solar cycle of cosmic antiprotons based on 1.1 × 10 6 events in the rigidity range from 1.00 to 41.9 GV. The $\bar{𝑝}$ fluxes exhibit distinct properties. The magnitude of the $\bar{𝑝}$ flux temporal variation is significantly smaller than those of 𝑝, 𝑒 − , and 𝑒 + . A hysteresis between the $\bar{𝑝}$ fluxes and the 𝑝 fluxes is observed, whereas the $\bar{𝑝}$ and 𝑒 − fluxes show a linear correlation. With a model-independent analysis, we found a universal relation between the shape of the rigidity spectrum and the magnitude of flux temporal variation over an 11-year Solar cycle for both positively and negatively charged particles. The simultaneous results on $\bar{𝑝}$ and 𝑝, 𝑒 − , and 𝑒 + provide unique information for understanding particle transport in the Solar System as a function of mass, charge, and spectral shape.

cosmic ray composition & spectra↗

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 $p_c$ 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 5 K. 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 OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Domain wall enabled steep slope switching in MoS2 transistors towards hysteresis-free operation

Abstract The device concept of ferroelectric-based negative capacitance (NC) transistors offers a promising route for achieving energy-efficient logic applications that can outperform the conventional semiconductor technology, while viable operation mechanisms remain a central topic of debate. In this work, we report steep slope switching in MoS 2 transistors back-gated by single-layer polycrystalline PbZr 0.35 Ti 0.65 O 3 . The devices exhibit current switching ratios up to 8 × 10 6 within an ultra-low gate voltage window of $$V_{{{\mathrm{g}}}} = \pm \! 0.5$$ V g = ± 0.5 V and subthreshold swing (SS) as low as 9.7 mV decade −1 at room temperature, transcending the 60 mV decade −1 Boltzmann limit without involving additional dielectric layers. Theoretical modeling reveals the dominant role of the metastable polar states within domain walls in enabling the NC mode, which is corroborated by the relation between SS and domain wall density. Our findings shed light on a hysteresis-free mechanism for NC operation, providing a simple yet effective material strategy for developing low-power 2D nanoelectronics.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Machine learning methods for probabilistic locked-mode predictors in tokamak plasmas

A rotating tokamak plasma can interact resonantly with the external helical magnetic perturbations, also known as error fields. This can lead to locking and then to disruptions. We leverage machine learning (ML) methods to predict the locking events. We use a coupled third-order nonlinear ordinary differential equation model to represent the interaction of the magnetic perturbation and the plasma rotation with the error field. This model is sufficient to describe qualitatively the locking and unlocking bifurcations. Here, we explore using ML algorithms with the simulation data and experimental data, focusing on the methods that can be used with sparse datasets. These methods lead to the possibility of the avoidance of locking in real-time operations. We describe the operational space in terms of two control parameters: the magnitude of the error field and the rotation frequency associated with the momentum source that maintains the plasma rotation. The outcomes are quan- tified by order parameters that completely characterize the state, whether locked or unlocked. We use unsupervised ML methods to classify locked/unlocked states and note the usefulness of a certain normalization of the order parameters. Three supervised ML classifiers are used in suite to estimate the probability of locking in the region of control parameter space with hysteresis, i.e., the set of control parameters for which both locked and unlocked states can exist. The results show that a neural network gives the best estimate of the locking probability. An analogy of the present locking model with the van der Waals equation of state is also provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Permeability Hysteresis From Microchannels Opening During Dissolution/Reprecipitation Cycle

Abstract Permeability is a critical parameter for geological resources characterization. Its evolution with respect to porosity is particularly interesting and many research initiatives focus on deriving such relationships, to understand some hydraulic impacts of microstructure alteration. Permeability evolution from chemical reactions, for instance, can become complex as flow channels may open during rock dissolution. In this contribution, we show that this phenomenon can lead to irregular porosity‐permeability curves and permeability hysteresis after reprecipitation. Current approaches describing permeability as a simple function of porosity can therefore not capture this behavior, and we advocate instead the use of dynamic modeling for such scenarios. We demonstrate our approach by modeling a dissolution/precipitation cycle for a unit cell pore channel and quantify the process at larger scale on three different rock samples, whose microstructures are reconstructed from segmented micro‐Computerised Tomography scans.

Lesueur, Martin↗

DOE OSTI2020

Room Temperature Quasi-static Characterization and Constitutive Model Parametrization of Flexible Polyurethane Foams of Different Densities Loaded in Different Orientations

This report describes the efforts to characterize and model General Plastics TF6070 and EF4000 flexible polyurethane foams under room temperature, large deformation quasi-static cyclic mechanical loading conditions. Densities from three to fifteen pound per cubic foot (PCF) are examined, which correspond to relative densities of approximately 4 to 20%. These foams are open cell structured and flexible at room temperature with a glass transition transition less than -30°C, and they fully recover their original shape when unloaded (at room temperature). Uniaxial compression tests were conducted with accompanying lateral image series for Digital Image Correlation (DIC) analysis with the goal of extracting transverse strain responses. Due to difficulties with DIC analysis at large strains, lateral strains were instead extracted for each test via edge tracking. The experimental results exhibit a nonlinear elastic response and anisotropic material behavior (particularly for the lower densities). Some hysteresis is observed that is different between the first and subsequent cycles of deformation indicating both a small degree of permanent damage (reduced stiffness during reloading) and viscoelasticity. These inelastic mechanisms are not considered in the modeling and calibration in this report. This work considers only the rate independent, room temperature foam behavior. Individual foam densities were calibrated for loading in two directions, parallel and perpendicular to the foam bubble rise direction, since the mechanical behavior is different in these two directions. The Flex Foam constitutive model was used for all parameterizations despite the fact that the model is isotropic. A review of the constitutive model is given as well as necessary data reduction procedures to parameterize it for each foam density and orientation are discussed. Finally, two different parameterizations are developed that take the undeformed foam density as an input that span all densities considered. These two parameterized models represent foams loaded in the rise and transverse directions respectively. We summarize the assumptions and limitations of the parameterizations provided in this report to guide analysis choices with them. All parameterizations presented herein have the following traits, room temperature, rate independent, damage-free, and non-dissipative . Isotropy (even if they are representing anisotropic data). Supplied Sierra Solid Mechanics Flex Foam Model Inputs are in units: pounds, inches, Celsius, and seconds

36 MATERIALS SCIENCE↗

DOE OSTI2024

Parametric reduced-order modeling for component-oriented treatment and localized nonlinear feature inclusion

Abstract We propose coupling a physics-based reduction framework with a suited response decomposition technique to derive a component-oriented reduction (COR) approach, which is suitable for assembly systems featuring localized nonlinearities. Dependencies on influencing parameters are injected into the reduced-order model (ROM), thus ensuring robustness and validity over a domain of parametric inputs, while capturing nonlinear effects. The implemented approach employs individual component modes to capture localized features while additionally relying on reduced modes of a global nature to approximate the system’s dynamics accurately. The global modes are derived from a linear monolithic system, defined as a result of a coordinate separation scheme, which permits the proposed COR-ROM to naturally couple the response between linear and nonlinear subdomains. The derived low-order representation utilizes a proper orthogonal decomposition projection and is additionally reinforced with the inclusion of a hyper-reduction technique to capture the underlying high-fidelity model response while providing accelerated computations. The resulting approach is exemplified in the synthetic case studies of a four-story shear frame with multiple nonlinear regions driven by hysteresis and a large-scale kingpin connection featuring plasticity.

Vlachas, Konstantinos (ORCID:000000029124219X)↗

DOE OSTI2025

Proximity ferroelectricity in wurtzite heterostructures

Proximity ferroelectricity is an interface-associated phenomenon in electric-field-driven polarization reversal in a non-ferroelectric polar material induced by one or more adjacent ferroelectric materials. Here, in this paper, we report proximity ferroelectricity in wurtzite ferroelectric heterostructures. In the present case, the non-ferroelectric layers are AlN and ZnO, whereas the ferroelectric layers are Al 1−x B x N, Al 1−x Sc x N and Zn 1−x Mg x O. The layered structures include nitride–nitride, oxide–oxide and nitride–oxide stacks that feature two-layer (asymmetric) and three-layer (symmetric) configurations. Ferroelectric switching in both layers is validated by multimodal characterization methods, including polarization hysteresis, anisotropic chemical etching, second harmonic generation, piezo response force microscopy, electromechanical testing and atomic resolution polarization orientation imaging in real space by scanning transmission electron microscopy. We present a physical switching model in which antipolar nuclei originate in the ferroelectric layer and propagate towards the internal non-ferroelectric interface. The domain wall leading edge produces elastic and electric fields that extend beyond the interface at close proximity, reducing the switching barrier in the non-ferroelectric layer, and allowing complete domain propagation without breakdown. Density functional theory calculations of polymorph energies, reversal barriers and domain wall energies support this model. Proximity ferroelectricity enables polarization reversal in wurtzites without the chemical or structural disorder that accompanies elemental substitution, opening new questions and opportunities regarding interface-based ferroelectricity.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Quantum Barkhausen noise induced by domain wall cotunneling

Most macroscopic magnetic phenomena (including magnetic hysteresis) are typically understood classically. Here, we examine the dynamics of a uniaxial rare-earth ferromagnet deep within the quantum regime, so that domain wall motion, and the associated hysteresis, is initiated by quantum nucleation, which then grows into large-scale domain wall motion, which is observable as an unusual form of Barkhausen noise. We observe noncritical behavior in the resulting avalanche dynamics that only can be explained by going beyond traditional renormalization group methods or classical domain wall models. We find that this “quantum Barkhausen noise” exhibits two distinct mechanisms for domain wall movement, each of which is quantum-mechanical, but with very different dependences on an external magnetic field applied transverse to the spin (Ising) axis. These observations can be understood in terms of the correlated motion of pairs of domain walls, nucleated by cotunneling of plaquettes (sections of domain wall), with plaquette pairs correlated by dipolar interactions; this correlation is suppressed by the transverse field. Similar macroscopic correlations may be expected to appear in the hysteresis of other systems with long-range interactions.

Science & Technology - Other Topics↗

DOE OSTI2021

Hysteresis curves reveal the microscopic origin of cooperative CO 2 adsorption in diamine-appended metal–organic frameworks

Diamine-appended metal-organic frameworks (MOFs) of the form Mg 2 (dobpdc)(diamine) 2 adsorb CO 2 in a cooperative fashion, exhibiting an abrupt change in CO 2 occupancy with pressure or temperature. This change is accompanied by hysteresis. While hysteresis is suggestive of a first-order phase transition, we show that hysteretic temperature-occupancy curves associated with this material are qualitatively unlike the curves seen in the presence of a phase transition; they are instead consistent with CO 2 chain polymerization, within one-dimensional channels in the MOF, in the absence of a phase transition. Here, our simulations of a microscopic model reproduce this dynamics, providing a physical understanding of cooperative adsorption in this industrially important class of materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Vertical Hydrologic Exchange Flows Control Methane Emissions from Riverbed Sediments

CH 4 emissions from inland waters are highly uncertain in the current global CH 4 budget, especially for streams, rivers, and other lotic systems. Previous studies have attributed the strong spatiotemporal heterogeneity of riverine CH 4 to environmental factors such as sediment type, water level, temperature, or particulate organic carbon abundance through correlation analysis. However, a mechanistic understanding of the basis for such heterogeneity is lacking. Here, in this study, we combine sediment CH 4 data from the Hanford reach of the Columbia River with a biogeochemical-transport model to show that vertical hydrologic exchange flows (VHEFs), driven by the difference between river stage and groundwater level, determine CH 4 flux at the sediment–water interface. CH 4 fluxes show a nonlinear relationship with the magnitude of VHEFs, where high VHEFs introduce O 2 into riverbed sediments, which inhibit CH 4 production and induce CH 4 oxidation, and low VHEFs cause transient reduction in CH 4 flux (relative to production) due to reduced advective CH 4 transport. In addition, VHEFs lead to the hysteresis of temperature rise and CH 4 emissions because high river discharge caused by snowmelt in spring leads to strong downwelling flow that offsets increasing CH 4 production with temperature rise. Our findings reveal how the interplay between in-stream hydrologic flux besides fluvial-wetland connectivity and microbial metabolic pathways that compete with methanogenic pathways can produce complex patterns in CH 4 production and emission in riverbed alluvial sediments.

54 ENVIRONMENTAL SCIENCES↗