Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic stability behaviors”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamics of Three-Dimensional Shock-Wave/Boundary-Layer Interactions

Advances in measuring and understanding separated, nominally two-dimensional (2D) shock-wave/turbulent-boundary-layer interactions (STBLI) have triggered recent campaigns focused on three-dimensional (3D) STBLI, which display far greater configuration diversity. Nonetheless, unifying properties emerge for semi-infinite interactions, taking the form of conical asymptotic behavior where shock-generator specifics become insignificant. The contrast between 2D and 3D separation is substantial; the skewed vortical structure of 3D STBLI reflects the essentially 2D influence of the boundary layer on the 3D character of the swept shock. As with 2D STBLI, conical interactions engender prominent spectral content below that of the turbulent boundary layer. However, the uniform separation length scale, which is crucial to normalizing the lowest-frequency dynamics in 2D STBLI, is absent. Comparatively, the spectra of 3D STBLI are more representative of the mid-frequency, convective, shear-layer dynamics in 2D, while phenomena associated with 2D separation-shock breathing are muted. Asymptotic behavior breaks down in many regions important to 3D-STBLI dynamics, occurring in a configuration-dependent manner. Aspects of inceptive regions near shock generators and symmetry planes are reviewed. Focused efforts toward 3D modal and nonmodalanalyses, moving-shock/boundary-layer interactions, fluid/structure interactions, and flow control are suggested as directions for future work.

42 ENGINEERING↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗
DOE OSTI2022

Depletion-Driven Assembly of Polymer-Coated Nanocrystals

Depletion-driven assembly has been widely studied for micron-sized colloids, but questions remain at the nanoscale where the governing physics are impacted by stabilizing surface ligands or wrapping polymers, whose dimensions are similar to those of the colloidal core and the depletant. Here, we probe how wrapping colloidal tin-doped indium oxide nanocrystals with polymers affects their depletion-induced interactions and assembly in solutions of polyethylene glycol. Copolymers of polyacrylic acid grafted with polyethylene oxide provide nanocrystal wrappings with different effective polymer graft densities and molecular weights. (Ultra) small angle X-ray scattering, coarse-grained molecular dynamics simulation, and molecular thermodynamic theory were combined to analyze how depletant size and polymer wrapping characteristics affect depletion interactions, structure, and phase behavior. The results show how depletant molecular weight, as well as surface density and molecular weight of polymer grafts, set thresholds for assembly. These signatures are unique to depletion-driven assembly of nanoscale colloids and mirror phase behaviors of grafted nanoparticle–polymer composites. Optical and rheological responses of depletion-driven assemblies of nanocrystals with different polymer shell architectures were probed to learn how their structural differences impact properties. Finally, we discuss how these wrinkles for depletion-driven assembly at the nanoscale may provide fresh opportunities for designing responsive depletion interactions and dynamically reconfigurable materials.

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

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗
DOE OSTI2021

Recycled polyethylene and crumb rubber composites modified asphalt with improved aging resistance and thermal stability

Crumb rubbers (CR) and recycled Polyethylene (PE) are commonly used as modifiers to improve the asphalt properties. However, the addition of these modifiers individually may result in phase separation of PE and inadequate degradation of CR, which may not simultaneously satisfy the rheological and aging resistance of the binder blends. In this study, a type of CR/PE composites manufactured through melt extrusion procedures was adopted for the improvement of the modified binder blends. The modified binder blends after short-term (RTFO) and long-term (PAV) aging were collected for dynamic shear rheometer and bending beam rheometer tests. FTIR tests were also used to quantify the aging resistance, while differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) were applied to evaluate the thermal behaviors of the binder blends. Optical microscopy tests were used to explore the modification mechanism and visualize the microstructures of original and PAV aged binder blends. Test results showed that CR/PE composites could significantly improve the aging resistance of the binder blends. The dispersed PE particles still interacted with each other and meanwhile worked with the diffused rubbery phase, which could form a similar TPEs network. The oxidation and high-temperature aging of such binder blends would strengthen the network that resisted the aging severity of the binder blends. Furthermore, it was also revealed that CR/PE modified asphalt could sustain more heat before a rapid degradation and the well-dispersed PE particles could interact with the diffused CR particles to form a network to resist aging.

02 PETROLEUM↗
DOE OSTI2024

Constructing Crown Ether-Based Supramolecular Adhesives with Ambient-Temperature Applicable and Durable Adhesion

Supramolecular adhesives are attracting considerable attentions owing to their dynamic and reversible bonding, while the complex curing conditions and poor performance under extreme conditions severely restrict their applicability. Herein, a series of tetra-crown ether-terminated polyethylene glycol (TCE-PEG) adhesives are presented in which plenty of supramolecular interactions, including Van der Waals, hydrogen bonding, π–π staking, metal coordination, and electrostatic ones, provide strong adhesion on multiple types of substrates. Highly relevant for practical applications, ultra-strong adhesivity of TCE-PEG can be achieved without requiring additional ultraviolet radiation, heat, or chemical treatment at ambient temperature. It also exhibits durable adhesive behavior, excellent underwater stability, and low-temperature tolerance. Especially, by incorporating ionic polymer fragments of polyethyleneimine hydrochloride, additional hydrogen bonding between secondary amines and crown ether rings can be induced to enhance the lap-shear strength to 7.21 MPa due to the formation of supramolecular cross-linked network. Finally, this work pioneers a rather unique supramolecular approach for easily formed, ultra-strong crown-ether-based adhesives bearing great potential for applications in critical environments of limited heat source and unallowed solvent usage.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhanced Frequency Support Scheme of Generic Inverter-Based Resource Models for Renewable-Dominated Power Grids

The frequency response of SG-dominated power grids is predictable ahead of an occurrence of a frequency event because the frequency response of SGs is consistent, and it can be inferred from the swing equation [1]. However, increasing the portion of IBRs in an SG-dominated power grid might make the characteristics of the conventional power grids no longer valid because this changing resource mix affects grid dynamics and controls [2]. Thus, maintaining these characteristics greatly benefits the control and operation of the power grids with high penetration of IBRs. To maintain these characteristics in IBR-dominated power grids, IBRs should have frequency response capability similar to that of an SG. The WECC modeling validation subcommittee has developed generic IBR models for large system planning [3]-[5]. These models can represent various vendors' dynamic behavior for WTG, PV, and ESS [5]. The current generic IBR models approved by WECC can provide frequency response only from droop control loops in REPC models [6], [7]. The contribution of the loops is proportional to the frequency deviation from the nominal frequency. Thus, it presents an insufficient contribution to arrest frequency variation compared to the frequency response of SGs because it allows a high ROCOF in the early stage of frequency events. This shortfall will become greater as the PL of IBRs increases in power grids. Controller enhancement for the generic IBR models is required to secure the frequency stability under high PL of IBRs as in the SG-dominated power grids. This paper proposes a control extension for the generic IBR models to enhance the frequency support capabilities and discusses the classification of frequency support for the different types of IBR considering their operating constraints. An inertial control scheme is implemented in the REPC and REEC models of the generic IBR models to achieve these objectives. The inertial control scheme includes the following stages: Control area data acquisition, inertia time constant estimation, IBR-related constraint check, IBR contribution determination, and inertial response provision. In the scheme, a REPC acquires control area data from a system operator and estimates a total inertia time constant for the control area the applicable IBR power plant belongs. Then, the estimated inertial time constant is transferred to each IBR controller—REEC—within the power plant. Each REEC checks the availability of applicable IBR for inertial response participation. If the IBR is available, the REEC amplifies the estimated inertial time constant to utilize it for inertial response provision. In this way, the proposed inertial response scheme extends the functionality of the generic IBR models to provide SG-like frequency response within their constraints. Various scenarios considering different IBR types, IBR penetration levels, and frequency control schemes were simulated and compared in an IEEE 39-bus system using PSCAD simulator to verify the effectiveness of the proposed scheme.

Kim, Jinho↗
DOE OSTI2025

Probing the Self-Assembly dynamics of cellulose nanocrystals by X-ray photon correlation spectroscopy

Charge-stabilized colloidal cellulose nanocrystals (CNCs) can self-assemble into higher-ordered chiral nematic structures by varying the volume fraction. The assembly process exhibits distinct dynamics during the isotropic to liquid crystal phase transition, which can be elucidated using X-ray photon correlation spectroscopy (XPCS). Anionic CNCs were dispersed in propylene glycol (PG) and water spanning a range of volume fractions, encompassing several phase transitions. Coupled with traditional characterization techniques, XPCS was conducted to monitor the dynamic evolution of the different phases. Additionally, simulated XPCS results were obtained using colloidal rods and compared with the experimental data, offering additional insights into the dynamic behavior of the system. The results indicate that the particle dynamics of CNCs undergo a stepped decay in three stages during the self-assembly process in PG, coinciding with the observed phases. The phase transitions are associated with a total drop of Brownian diffusion rates by four orders of magnitude, a decrease of more than a thousand times slower than expected in an ideal system of repulsive Brownian rods. Given the similarity in the phase behaviors in CNCs dispersed in PG and in water, we hypothesize that these dynamic behaviors can be extrapolated to other polar solvent environments. Importantly, these findings represent the direct measurement of CNC dynamics using XPCS, underscoring the feasibility of directly assessing the dynamic behavior of other rod-like colloidal suspensions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Unexpected electrochemical behavior of an anolyte redoxmer in flow battery electrolytes: solvating cations help to fight against the thermodynamic–kinetic dilemma

Redoxmers are redox-active molecules that can store energy in electrolytes for redox flow batteries (RFBs), and their electrochemical properties are significantly affected by the choice of supporting electrolytes. Herein, we use 2,1,3-benzothiadiazole (BzNSN) as a model system to scrutinize the supporting electrolyte impact. By systemically varying the components of supporting salts, BzNSN not only shows substantial redox potential shifts but also exhibits varying electrochemical stabilities. Specifically, changing the size of cations can effectively alter the coordination between the supporting salt and BzNSN species. From Li + , Na + , K + , to NEt 4 + , the redox potential of BzNSN shifts negatively, from -1.63 V to -1.82 V vs. Ag/Ag + . Molecular dynamics and density functional theory simulations revealed that smaller cations, like Li + , are closer to the charged BzNSN when coordinated, implying stronger coordination, while larger cations, like K + and NEt 4 + , are farther away. Interestingly, the large cation electrolytes also lead to much improved electrochemical stability, evidenced by the extraordinarily enhanced kinetic lifetime from electron paramagnetic resonance measurement. This study demonstrates the first example of tuning an anolyte redoxmer toward a concurrent improvement of lowered redox potentials AND enhanced calendar lives via solvation means, which is usually constrained by the thermodynamic–kinetic relation.

2,1,3-benzothiadiazole↗
DOE OSTI2023

Granger Causality for prediction in Dynamic Mode Decomposition: Application to power systems

Here, the dynamic mode decomposition (DMD) technique extracts the dominant modes characterizing the innate dynamical behavior of the system within the measurement data. For appropriate identification of dominant modes from the measurement data, the DMD algorithm necessitates ensuring the quality of the input measurement data sequences. On that account, for validating the usability of the dataset for the DMD algorithm, the paper proposed two conditions: Persistence of excitation (PE) and the Granger Causality Test (GCT). The virtual data sequences are designed with the hankel matrix representation such that the dimensions of the subspace spanning the essential system modes are increased with the addition of new state variables. The PE condition provides the lower bound for the trajectory length, and the GCT provides the order of the model. Satisfying the PE condition enables estimating an approximate linear model, but the predictability with the identified model is only assured with the temporal causation among data searched with GCT. The proposed methodology is validated with the application for coherency identification (CI) in a multi-machine power system (MMPS), an essential phenomenon in transient stability analysis. The significance of PE condition and GCT is demonstrated through various case studies implemented on 22 bus six generator system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Unexpected Kinetic Solvent Effects Enhance Activity and Selectivity in Biphasic Systems

Biphasic dehydration of fructose to 5-hydroxymethylfurfural (HMF) has shown unprecedented increases in productivity, but a mechanistic understanding is lacking. Herein, we couple fast experimental reaction kinetics, multiscale modeling (phase behavior, classical molecular dynamics(MD), and quantum mechanics/molecular mechanics MD), in situ sampling, and IR and 13 C-NMR spectroscopy to elucidate the complex effects of nonpolar extracting organic solvents on the kinetics of fructose dehydration. We show that these organic solvents can reach significant mutual solubility with water at reaction temperatures, enabling the partition of the sugar and catalyst into the extracting phase. In the organic-rich environment, the dehydration of fructose proceeds faster and more selectively than in water due to increased relative abundance of the reactive furanose isomer, enhanced water–catalyst–substrate interactions driven by nanophase separation, and higher product stability stemming from preferential solvation. Furthermore, we demonstrate that these solvent effects impact other critical biphasic reactions in biomass upgrading and provide qualitative principles for solvent selection.

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

Effect of iodides on thermal behavior and phase partitioning in LiCl-KCl

Liquid-fueled molten salt reactors (MSRs) are designed to operate with fissile materials and, ultimately, fission products dissolved in the primary molten salt coolant. Understanding the speciation and transport of iodine—a high-yield fission product—is essential because this element’s accidental release poses significant environmental concerns due to its capacity to be readily absorbed by the human thyroid gland. Here, we report the impact of iodide species (LiI and KI) on phase transitions, phase distribution, and phase stability in LiCl–KCl-eutectic salt mixtures. The study employed a combination of computational and experimental techniques, including thermodynamic FactSage calculations, differential scanning calorimetry, and high-temperature X-ray diffraction. The results indicate that the presence of iodide (10–25 wt%) significantly alters the melting behavior of the LiCl–KCleutectic system. Adding 10 wt% LiI has a more-pronounced effect than 10 wt% KI, as LiI converts to KI, leading to formation of LiCl, thereby, altering the LiCl-KCl ratio which significantly affects the melting temperature of the mixture. Furthermore, the evolution of crystalline structure, solid-fraction composition, and the dynamics of mixed-halide solid–liquid partitioning as a function of temperature indicate the potential for selective iodide separation from chloride-salt mixtures via solid–liquid separation techniques. Overall, the presented findings provide valuable insights that are beneficial for the design and operation of MSRs, as well as for the safe handling and effective processing of used nuclear fuel using advanced pyrochemical techniques.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Phase instability-coupled fracture behavior in garnet LLZO solid electrolytes: a machine learning-enabled atomistic study

Fracture in the garnet-type solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) poses a critical threat to both the performance and safety of solid-state batteries. To unravel the coupled chemomechanical processes that govern fracture in LLZO under external loading, we carry out large-scale molecular-dynamics simulations with a validated machine-learning force field. Our results demonstrate that triaxial stresses at crack flanks trigger a localized cubic-to-tetragonal phase transformation, which is accompanied by Li-ion rearrangement. The emergent tetragonal domains feature lattice contraction normal to the fracture plane, imposing coherent misfit strains that provide an additional driving force for further crack propagation. Crucially, introducing Li deficiencies stabilizes the cubic phase, postponing the phase transition and thereby delaying fracture initiation. These findings highlight the role of intrinsic phase instability in dictating LLZO's fracture resistance and its critical connection to local Li concentration. This chemomechanical coupling points toward targeted strategies to enhance the mechanical robustness of garnet electrolytes, including tuning Li content, ensuring dopant homogeneity, and refining processing protocols.

Materials science↗
DOE OSTI2021

Behavior of solitary waves of coupled nonlinear Schrödinger equations subjected to complex external periodic potentials with odd-$\mathcal{PT}$ symmetry

In this work, we discuss the response of both moving and trapped solitary wave solutions of a two-component nonlinear Schrödinger system in 1 + 1 dimensions to an odd-$\mathcal{PT}$ external periodic complex potential. The dynamical behavior of perturbed solitary waves is explored by conducting numerical simulations of the nonlinear system and using a collective coordinate variational approximation. We present case examples corresponding to choices of parameter values and initial conditions involved therein. The results of the collective coordinate approximation are compared against numerical simulations where we observe qualitatively good agreement between the two. Unlike the case for a single-component solitary wave in a complex periodic $\mathcal{PT}$-symmetric potential, the collective coordinate equations do not have a small oscillation regime, and initially the height of the two components changes in opposite directions often causing instability. We find that the dynamic stability criteria we have used in the one-component case are a good indicator for the onset of dynamic instabilities in the present setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

In situ probing of structure and deagglomeration of SnO 2 colloids via small-angle X-ray scattering

Transforming dry nanopowders into stable colloidal dispersions remains challenging due to the cohesive forces between the nanoparticles (NPs) which promote agglomeration. Effective dispersion and deagglomeration of these agglomerates is a critical process in the formulation and preparation of nanoparticle-based functional materials via the colloidal route. Understanding the deagglomeration dynamics provides information to improve microstructural quality in various applications by enabling engineering of agglomerate size, structure and morphology. However, the deagglomeration process dynamics with respect to the evolution of the fractal agglomerate structures, particularly for very small NPs, is still poorly understood and requires further investigation. This study employs in situ small-angle X-ray scattering (SAXS) to investigate the sonication-induced deagglomeration of SnO 2 NPs. Electrostatically stabilized SnO 2 colloids with varying primary particle size (6–21 nm) are investigated in a specifically designed in situ cell using synchrotron-based SAXS to study the influence of sonication time and intensity on the nanoscaled agglomerates. Complete structural analysis via SAXS reveals a direct correlation between changes in agglomerate size and structure, size-dependent deagglomeration behavior and a dependence on the overall energy introduced during sonication into the dispersion regardless of the actual power as well as ultrasonic process parameters in case of SnO 2 NPs. The results suggest that control of the dispersion process during ultrasonic deagglomeration results in tailoring agglomerates with respect to size and structure.

Deagglomeration↗
DOE OSTI2022

Structural and dynamic properties of solvated hydroxide and hydronium ions in water from ab initio modeling

Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water from ab initio molecular dynamics (AIMD) has been a challenging task. The difficulty mainly comes from a lack of accurate and efficient exchange–correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange–correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice for describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. In conclusion, our AIMD simulations indicate that hydroxide diffuses more slowly than hydronium in water, which is consistent with the experimental results.

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

Interface-Driven Chiral Magnetism in Ultrathin Metallic Ferromagnets: Towards Skyrmion Spintronics

Magnetic skyrmions are localized chiral magnetic textures in the form of nanoscale vortices or bubbles that are topologically protected from being ‘unwound’. Their topological nature gives rise to rich behaviors including ordered lattice formation, emergent electrodynamics and robust current-driven displacement at remarkably low current densities. However, magnetic skyrmions have so far been restricted to just a few materials and observed only at low temperatures, limiting the experimental accessibility and technological application of these unique topological objects. This project aimed to realize magnetic skyrmions at room temperature in a new class of engineered chiral ferromagnets that exploit broken mirror symmetry at interfaces to generate helical magnetic order. Ultrathin magnetic heterostructures were engineered and patterned into laterally-confined nanostructures in which skyrmions were stabilized, manipulated, and detected for the first time under ambient conditions. Static and dynamic properties were examined using advanced x-ray imaging to provide a fundamental understanding of topological magnetic phases. The interactions between magnetic skyrmions and electron charge and spin currents were studied to identify fundamental new physics and enable device applications. These experimental studies, supported by detailed modeling, will help launch a new subfield of skyrmion-based spintronics in which individual skyrmions can be used to encode, store, and transport information in high-performance, low power memory and logic devices.

36 MATERIALS SCIENCE↗