SEARCH · Engineering Papers

Results for “local minimum”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Multilevel Robustness for 2D Vector Field Feature Tracking, Selection and Comparison

Abstract Critical point tracking is a core topic in scientific visualization for understanding the dynamic behaviour of time‐varying vector field data. The topological notion of robustness has been introduced recently to quantify the structural stability of critical points, that is, the robustness of a critical point is the minimum amount of perturbation to the vector field necessary to cancel it. A theoretical basis has been established previously that relates critical point tracking with the notion of robustness, in particular, critical points could be tracked based on their closeness in stability, measured by robustness, instead of just distance proximity within the domain. However, in practice, the computation of classic robustness may produce artifacts when a critical point is close to the boundary of the domain; thus, we do not have a complete picture of the vector field behaviour within its local neighbourhood. To alleviate these issues, we introduce a multilevel robustness framework for the study of 2D time‐varying vector fields. We compute the robustness of critical points across varying neighbourhoods to capture the multiscale nature of the data and to mitigate the boundary effect suffered by the classic robustness computation. We demonstrate via experiments that such a new notion of robustness can be combined seamlessly with existing feature tracking algorithms to improve the visual interpretability of vector fields in terms of feature tracking, selection and comparison for large‐scale scientific simulations. We observe, for the first time, that the minimum multilevel robustness is highly correlated with physical quantities used by domain scientists in studying a real‐world tropical cyclone dataset. Such an observation helps to increase the physical interpretability of robustness.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Rye-soybean double-crop: planting method and N fertilization effects in the North Central US

Double-cropping winter rye cover crops (CC) with soybean in the North Central US could help with the global effort to sustainably intensify agriculture. Studies addressing the management of these systems are limited. Therefore, a field study was conducted from 2017 to 2019 in Central Iowa, US to evaluate winter rye CC biomass production, aboveground N accumulation, estimated economics, estimated within-field energy balance and estimated greenhouse gas (GHG) emissions under three N application rates (0, 60, 120 kg N ha –1 ) and three planting methods (pre- and post-harvest broadcast and post-harvest drilling). Averaged over N rates, all planting methods resulted in >5.0 Mg ha –1 year –1 rye aboveground biomass dry matter. Averaged over the 2-year study and compared with unfertilized treatments, applying 60 kg N ha –1 produced 1.1 Mg ha –1 more aboveground biomass (6.1 vs 5.0 Mg ha –1 ), accumulated 30 kg ha –1 more N in aboveground biomass (88 vs 58 kg N ha –1 ), and led to 20 GJ ha –1 more net energy. Biomass production was not significantly higher with 120 kg N ha –1 compared with the 60 kg N ha –1 rate. Even when accounting for an estimated 0.75 Mg ha –1 of above ground rye biomass left in the field after harvesting, more N was removed than applied at the 60 kg N ha –1 rate. The minimum rye prices over the 2-year study needed for double-cropping winter rye CC to be profitable (breakeven prices) averaged $117 and $104 Mg –1 for the 0 and 60 kg N ha –1 rates, which factors in estimated soybean yield reductions in 2019 compared with local averages but not off-site transportation. GHG emissions were estimated to increase approximately threefold between the unfertilized and 60 kg N ha –1 rates without considering bioenergy offsets. While environmental tradeoffs need further study, results suggest harvesting fertilized rye CC biomass before planting soybean is a promising practice for the North Central US to maximize total crop and net energy production.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Existence of a Nonzero Worst‐Case ACH for Short‐Term Exposure in Ventilated Indoor Spaces

A well‐ventilated room is essential to reduce the risk of airborne transmission. As such, the scientific community sets minimum limits on ventilation with the idea that increased ventilation reduces pathogen concentration and thus reduces the risk of transmission. In contrast, the upper limit on ventilation is usually determined by human comfort and the need to reduce energy consumption. While average pathogen concentration decreases with increased ventilation, local concentration depends on multiple factors and may not follow the same trend, especially within short exposure times over large separation distances. Here, we show through experiments and high‐fidelity simulations the existence of a worst‐case ventilation where local pathogen concentration increases near the receiving host. This occurs during the type of meetings that were recommended during the pandemic (and in some cases solely authorized) with reduced occupancy adhering to social distancing and short exposure times below 20 minutes. We maintain that for cases of high occupancy and long exposure time, increased ventilation remains necessary.

Construction & Building Technology↗

DOE OSTI2025

Dislocation-Grain Boundary Interaction Dataset for FCC Cu

Interactions between dislocations and grain boundaries play a major role in controlling the strength and ductility of structural materials. Experimentally, assessing and probing geometric and stress-based criteria at the local level for dislocation transmission through grain boundaries remains challenging. Therefore, there have been many efforts to systematically generate datasets of dislocation-grain boundary interactions (DGI) via computational models such as molecular dynamics simulations. So far, most DGI datasets have focused only on the subset of nominal minimum-energy grain boundary structures, which limits their applicability, especially to materials processed far from equilibrium. We present a comprehensive database of dislocation-grain boundary interactions for edge, screw, and 60° mixed dislocation with 330 <110> and 257 <112> symmetric tilt grain boundaries (total of 587) in FCC Cu consisting of 73 minimum-energy grain boundary structures and 514 metastable structures. The dataset contains the outcomes for 5234 unique interactions for various dislocation types, grain boundary structures, and applied shear stresses.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A numerical and experimental approach to compare the effect of sample thickness in small in-situ SEM and large ex-situ tensile testing in Alloy 709

ASTM standards for tensile tests define specific sample size requirements regardless of grain size of the material. However, sample size requirements for testing should be considered in conjunction with the number of grains within its reduced cross-sectional area. This is particularly important for in-situ Scanning Electron Microscope (SEM) tensile tests, as they have to be conducted on smaller samples. In this study, a comprehensive experimental and numerical evaluation of the effect of specimen thickness (and the resulting number of grains within the cross-section) were conducted using in-situ SEM tensile test (on submillimeter thick samples) and ex-situ tensile tests (on samples of 0.68–5.9 mm thickness) and the results are compared with FEM simulations outcome. All tests were conducted at room temperature and the results are correlated to the number of grains within the thickness of specimens. The tensile test results indicated that even though the 0.2% proof stress and the tensile strength do not vary for the current range of thicknesses (with number of grains in the cross-section varying from 13 to 118), a difference in necking mechanism exists. Post tensile strength, thinner samples undergo more shear failure and diagonal localized necking whereas thicker samples experience more diffuse necking indicated by decreased area of shear failure at the edges and increased area of dimpled tensile failure at the center of the sample. Further, FEM results complement the experimental findings by showing the formation of conjugated localized shear bands on the upper and lower surface in the 0.68 mm thick sample and shear bands combining to form singular shear bands in thicker samples. These results also confirm the validity of the in-situ SEM tensile tests conducted on thinner samples as long as the required minimum number of grains exist within the cross-section.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A high-pressure Raman study of FeTiO 3 ilmenite: Fermi resonance as a manifestation of Fe-Ti charge transfer

In this work, we investigated the 300 K high-pressure behavior of ilmenite using Raman spectroscopy to 54 GPa. Upon compression, we observe a Fermi resonance between the lowest frequency A g symmetry peaks (ν 4 and ν 5 ) between ~ 10 and ~ 30 GPa: bands that involve major components of Ti–O and Fe–O-related displacements, respectively. The peaks’ relative intensities switch at ~ 18 GPa and they also reach their minimum separation at ~ 20 GPa, indicating that their maximum resonance occurs between 18 and 20 GPa. The negative shift of the Ti–O-associated ν 4 vibration under compression is fully consistent with a shift in valence of Ti from 4 + to 3 + under compression. Anomalously small mode shifts of other, more localized vibrations are also consistent with a charge transfer from Fe to Ti under compression. At higher pressures, we have not found definitive evidence for a transition to the perovskite-structure at 300 K, which has been well characterized at high pressures and temperatures. At 40 GPa, we observe an apparent reversible disordering that persists up to our highest pressure. The 300 K mode shifts of the Raman active modes in FeTiO 3 under pressure are notably different from those of other ABO 3 compounds (where A = Mg, Mn and B = Ti, Si); in other ilmenite-structured compounds, the peaks shift at a faster rate and there has not been any observation of Fermi resonance. Thus, iron’s complex electronic structure, and its charge transfer with titanium, appears to play a primary role in the behavior of phonons in FeTiO 3 ilmenite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

On the energy landscape of symmetric quantum signal processing

Symmetric quantum signal processing provides a parameterized representation of a real polynomial, which can be translated into an efficient quantum circuit for performing a wide range of computational tasks on quantum computers. For a given polynomial f , the parameters (called phase factors) can be obtained by solving an optimization problem. However, the cost function is non-convex, and has a very complex energy landscape with numerous global and local minima. It is therefore surprising that the solution can be robustly obtained in practice, starting from a fixed initial guess Φ 0 that contains no information of the input polynomial. To investigate this phenomenon, we first explicitly characterize all the global minima of the cost function. We then prove that one particular global minimum (called the maximal solution) belongs to a neighborhood of Φ 0 , on which the cost function is strongly convex under the condition ‖ f ‖ ∞ = O ( d − 1 ) with d = d e g ( f ) . Our result provides a partial explanation of the aforementioned success of optimization algorithms.

Wang, Jiasu↗

DOE OSTI2023

Characterization of Orography-Influenced Riming and Secondary Ice Production and Their Effects on Precipitation Rates Using Radar Polarimetry and Doppler Spectra (CORSIPP-SAIL)

The Characterization of Orography-Influenced Riming and Secondary Ice Production and Their Effects on Precipitation Rates Using Radar Polarimetry and Doppler Spectra (CORSIPP) project was conducted to help improve the understanding of precipitation formation in orographically influenced terrain. Special focus is put on the two processes of riming and secondary ice production and their external drivers. Two instruments, a polarimetric W-band simultaneous transmission simultaneous reception (STSR) Doppler cloud radar manufactured by Radiometer Physics GmbH (RPG, instrument type RPG-FMCW-94-DP), from now on named LIMRAD94, and the video in situ snowfall sensor (VISSS), were deployed at the U. S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Surface Atmosphere Integrated Field Laboratory (SAIL) site in Gothic, Colorado between November 2022 and June 2023 during the second SAIL winter. Note that the exact dates of data availability differ between the instruments. Both instruments arrived at Gothic on November 2, 2022, 09:20 local time. VISSS, described by Maahn et al., is equipped with two camera systems with telecentric lenses. The two cameras are at a 90° angle to each other. This configuration allows for size-independent measurements by capturing images of hydrometeors from two sides at a high frame rate of 250 Hz. With a minimum detection size of 200 μm, VISSS provides valuable insights into particle size, number, shape, complexity, and fall velocity. The VISSS was deployed on the grassland next to the ARM facility with the amazing help of the ARM employees on site. The setup started on November 2, 2022, and was finished on November 5, 2022, without major problems. VISSS measurements were started on November 6, 2022. LIMRAD94 was installed on a scaffolding platform near Orehouse (Gothic) on November 9, 2022, with the great help of RMBL staff. LIMRAD94 was mounted on a cold temperature scanner prototype. After a short test of the setup on November 9, 2022, the digital control of the scanner elevation stopped working for (at that time) unknown reasons. All attempts to resolve the problem failed. This malfunction made it impossible to operate LIMRAD94 in scanning mode. The scanner was then manually moved to zenith pointing mode and between November 10 and November 15, 2022, vertical observations for polarimetric calibration were performed. On November 15, 2022, after the polarimetric calibration was applied, the scanner was manually moved to 40° elevation with azimuthal view towards the Ka-band ARM Zenith Radar (KAZR) and measurements were continued at constant elevation. Investigation of the scanner malfunction on February 6, 2023, by Benn Schmatz revealed a disengagement between the cogwheel of the elevation motor and the cogwheel moving the scanner in elevation. This mechanical issue was temporarily solved by re-engaging the cogwheels. This made the scanner operational again for about four weeks, until mechanical force disengaged the cogwheels again on March 15, 2023. This repeated scanner failure remained undetected for about three weeks until April 8, 2023; during this time the scanner was stuck at 72° elevation. However, the radar software continued to produce data files falsely indicating that the scanner was still operational. After the scanner failure was noticed, scanning was stopped again, and we returned to constant elevation measurements. On May 15, 2023, the cogwheels of the elevation motor were secured with additional screws sent by the manufacturer. At some point in May, the cogwheels of the azimuth motor were also disengaged by mechanical force, which still allowed for range height indicator (RHI) but no plan position indicator (PPI) scans in the last weeks of the campaign. The azimuth motor was repaired in Germany after the end of the campaign. Throughout the campaign, RMBL and ARM staff kept the radar and the VISSS free of snow.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Predictive Tools for Customizing Heat Treatment of Additively Manufactured Aerospace Components

Laser-bed powder fusion (LBPF) additive manufacturing is increasingly being used to produce components of complex geometries using the Ni-base superalloy Inconel 718. The composition and the microstructure of the alloy are currently well optimized for wrought components made using conventional manufacturing processes such as rolling, forging, extrusion, etc. The attractive mechanical properties of the alloy result from the underlying austenitic matrix with fine equiaxed grains, and a high density and uniform distribution of the precipitation hardening phase, γ". Heat treatment steps such as homogenization, solutioning and aging are well documented for the wrought alloy. However, when the same wrought alloy compositions are used for the additive manufacturing (AM) processes, the asprocessed microstructure is significantly different, because of the different thermal history associated with LBPF, including rapid solidification and multiple temperature excursions that lead to multiple re-melting and reheating in the solid state. Rapid solidification introduces potential non-equilibrium effects at the moving solid-liquid interfaces that impact the extent of solute segregation, as well as the morphology of the dendritic grains that form. In order to recover the target mechanical properties, AM components have to undergo post-process heat treatments. However, such heat treatments have to be custom designed for the AM process and the component geometry because of the expected vast differences in the microstructure at various locations of a component with complex geometry. The homogenization and precipitation steps should be optimized for the component so that target mechanical properties can be obtained throughout the part. The objective of this research is to utilize High Performance Computing in phase field simulations of microstructure evolution during post-processing of AM components. The physics-based modeling will be beneficial in reducing the experimental effort required for heat treatment process selection, optimization, and certification, thus leading to a significant reduction in energy consumption for AM and post-processing heat treatment. The optimization study will help identify heat treatments steps that are critical for development of a final desired microstructure with the minimum energy input. This combined with shortening of the production cycle (time-to-market) by reducing the number of failed parts (property targets), and reduction in the number of iterations for process optimization, will enable 30-40% savings in the energy costs. Phase field simulations of the degree of homogenization and the effect of local matrix composition on the nucleation and growth of competing precipitating phases were performed using the Microstructure Evolution Using Massively Parallel Phase Field Simulations code developed in-house at the Oak Ridge National Laboratory. The simulations were able to successfully capture the kinetics of nucleation and growth, and morphologies of various precipitating phases as a function of local matrix compositions and composition gradients characteristic of local microstructures arising from location-dependent variations in the thermal conditions. Future work will involve extending the simulations to a length scale consisting of multiple dendrites, so that the effect of homogenization on the coarsening of the dendrites can be simulated and used as an additional input to the optimization of the heat treatment process.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Influence of CO 2 -Regenerative Film Properties in Enhancing C 2+ Products Selectivity While Mitigating CO 2 Crossover

Zero-gap anion-exchange membrane electrode assembly (AEMEA) electrolyzers operating in alkaline media face challenges such as CO 2 crossover and salting-out. The bipolar membrane electrode assembly (BPMEA) electrolyzer, using DI water as the electrolyte, addresses both CO 2 crossover and salting-out issues. As cations, which are crucial in stabilizing the CO 2 R intermediates, are absent in the electrolyte, cations embedded in the AEM play an important role in dictating the CO 2 R selectivity and activity in BPMEA systems. So far, no systematic study has been conducted on the influence of cations embedded in AEMs on CO 2 selectivity and activity in BPMEA systems. Moreover, BPMEA systems impose an additional challenge: low stability due to delamination of the bipolar membrane caused by CO 2 regeneration at the membrane-membrane interface. To enhance the stability of the electrolyzer, a simple yet highly reproducible strategy for coating a porous CO 2 regenerative film on smooth Nafion 117 is demonstrated in this work, along with a systematic study of four different commercially available AEMs, composed of different cations and cation densities, for use in combination with Nafion 117 and copper catalyst at the cathode. We found that the PiperION membrane delivers selectivity and activity comparable to those of alkali-metal cations, owing to the enhanced local electric field resulting from the combined effects of a high positive charge on the N atom of the piperidinium cation and the high ion-exchange capacity. Further, we studied the influence of the thickness of the PiperION porous layer over Nafion 117 on CO 2 R selectivity and found that 70 μm is the minimum thickness to achieve maximum C 2+ products selectivity, reduced the CO 2 crossover to 5% from 25% at 4 SCCM and 150 mA cm –2 , and was stable for operation beyond 100 h. The selectivity of this system, compared with AEMEA, and stability outperformed both AEMEA and BPMEA. This study helps design more effective BPMs to inhibit CO 2 crossover while enabling stable and selective electrochemical CO 2 reduction to C 2+ hydrocarbons.

Cations↗

DOE OSTI2025

Low-Energy Isomers of the Magic Number H + (H 2 O) 21 Cluster

Electronic structure calculations are used to characterize low-energy isomers of H + (H 2 O) 21 . Eleven different classes of isomers, based on the (H 2 O) 20 pentagonal dodecahedron with the excess proton localized on the surface (as a hydrated hydronium ion) and the “extra” water molecule located in the interior of the cluster, are characterized. In 10 of these classes, the internal water molecule is engaged in six 5-membered rings, but in the remaining class, which is predicted to start at only 0.6 kcal/mol above the global minimum, the internal water is engaged in a 4- membered ring, an additional 6-membered ring, and four 5- membered rings. In addition, isomers with two 4-membered rings and two 6-membered rings on the cluster surface are predicted to start at only ∼1.3 kcal/mol above the lowest-energy dodecahedralbased structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Opportunities for wave energy in bulk power system operations

Wave energy resources have high, yet largely untapped potential as candidate generation technology. In this paper, we perform a data-driven analysis to characterize the impact of wave energy integration on bulk-scale power systems and market operations. Through data-driven sensitivity studies centered on an optimization-based production cost modeling formulation, our work characterizes the inflection point beyond which wave integration starts impacting power system operations, considering present day transmission infrastructure. Furthermore, our analysis also considers the joint effects of wave energy integration and system-wide transmission expansion. Finally, potential resilience scenarios such as wildfire-driven transmission contingencies and heat wave events are investigated, whereby the contributions of grid-integrated wave energy in alleviating the effects of the resilience events are analyzed. As our demonstration test bed, we consider a reduced-order network topology for the U.S. Western Interconnection with wave energy generation integrated at carefully selected sites across the coastal areas of Washington, Oregon, and northern California. Our results indicate that over a representative year of operations, wave energy integration systematically reduces locational marginal prices (LMPs) of energy and price volatility, especially during periods of high wave resource availability (winter months for the U.S. west coast). Average, maximum, and minimum of hourly LMPs over a typical year of operation was reduced by 2.95, 51.28, and 1.13 $\$$/MWh respectively (over a baseline scenario with no wave energy integration), when the selected network model had a total of 5000 MW wave power installed capacity during the representative year of study. The effects of wave energy integration can remain localized with existing transmission infrastructure (identified to be most pronounced in the Pacific Northwest region in the example we studied). However, with concurrent transmission expansion, the impacts of wave energy integration are likely to have a higher geographical spread. Our results also indicate that wave energy may be able to assist power system operations during resilience events such as major transmission contingencies and heat wave events, although such benefits might be dependent on factors such as proximity of affected area to wave resources, availability of adequate resource potential and adequate transmission capacity.

16 TIDAL AND WAVE POWER↗

DOE OSTI2025

Five Year Comparison of Mixing Height Determinations at the Savannah River Site

Air quality dispersion modeling is performed for the Savannah River Site (SRS) to demonstrate compliance with applicable regulations. The AMS/EPA Regulatory Model (AERMOD) modeling system is an EPA recommended model for air quality applications with a data preprocessor (AERMET) to incorporate meteorological data collected on site. AERMET parameterizes or calculates meteorological variables that are not directly measured onsite. One of the parameters estimated by AERMET is the atmospheric mixing height. While the mixing height is not currently a measurement input into AERMET, SRS has the capability to measure the local mixing height. The Savannah River National Laboratory (SRNL) operates a Vaisala CL31 Lidar Ceilometer which estimates mixing height from aerosol backscatter. This study compares the parameterized mixing height from AERMET to the ceilometer estimated mixing height for the current regulatory period at SRS incorporating data from 2015-2019. Results from this study showed the average daily minimum values (morning) from AERMET were an order of magnitude lower than the commonly used Holzworth (1972) method and the ceilometer estimated mixing heights. Additionally, on average, the ceilometer exhibited a daily maximum mixing height value that occurred 1-3 hours later than the AERMET estimated maximum. This difference is likely due to the nighttime atmospheric mixing height assumptions and calculations used by AERMET. The AERMET algorithm cuts off mixing height growth at sunset while the ceilometer data show ongoing evening convection typical of the southeastern United States. These results suggest that the AERMET parametrization scheme assumptions may not be representative of a forested landscape and evening convection which could account for more mixing overnight. The results obtained in this study are significant for air dispersion modeling applications for regulatory purposes and worker safety. Mixing height can impact model estimated pollutant concentrations. A greater mixing height will provide more volume for pollutant dispersion. This report documents efforts to quantify the dependence of mixing height inputs toward a conservative estimated pollutant concentration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Observational study of the South American low-level jet during the SALLJEX

The South American Low-Level Jet (SALLJ) east of the Andes is a narrow maximum northerly wind present just above the boundary layer. It is an important component of the tropical-extratropical exchange of heat and moisture in South America and favors deep moist convection in Southeastern South America (SESA). In the present study, the SALLJ characteristics during the 2002-2003 summer season are investigated using upper-air measurements collected during the SALLJ Experiment (SALLJEX). At most stations, the highest wind speed was observed below 3000 m between 06 and 12 UTC, mainly in association with a maximization of the northerly component of the wind while the zonal component reached its minimum. Here, an increase in the maximum wind speed from 10°S to Santa Cruz de la Sierra’s latitude and from there a weakening southward was found throughout the SALLJEX network. Additionally, SALLJEX observations allowed to locate SALLJ core purely from observations in a region just east of the Andes in agreement with previous studies. A new localized kinematic classification was developed and applied in this study to analyze the relationship between the SALLJ intensity and the thermodynamic properties of the lower layers of the atmosphere. Weak, moderate, and strong SALLJ categories were defined with a local threshold of the northerly wind speed maximum at each station. The strongest SALLJs were frequently observed at nighttime, while weaker SALLJs were likely to occur at any time of the day. A stabilization of the nocturnal boundary layer was observed between weak and moderate cases whereas a destabilization between moderate and strong cases were found near the time of the wind speed maximum (06 UTC) at Santiago del Estero and Resistencia stations. A deepening of the convective boundary layer and an increase of its mean temperature between weak and strong category during the previous afternoon to nocturnal SALLJs, reveal a possible contribution to the mechanisms that generate the maximum wind speed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Energetic Electron Precipitation Driven by Electromagnetic Ion Cyclotron Waves from ELFIN’s Low Altitude Perspective

We review comprehensive observations of electromagnetic ion cyclotron (EMIC) wave-driven energetic electron precipitation using data collected by the energetic electron detector on the Electron Losses and Fields InvestigatioN (ELFIN) mission, two polar-orbiting low-altitude spinning CubeSats, measuring 50-5000 keV electrons with good pitch-angle and energy resolution. EMIC wave-driven precipitation exhibits a distinct signature in energy-spectrograms of the precipitating-to-trapped flux ratio: peaks at >0.5 MeV which are abrupt (bursty) (lasting ~17 s, or ΔL ~ 0.56) with significant substructure (occasionally down to sub-second timescale). We attribute the bursty nature of the precipitation to the spatial extent and structuredness of the wave field at the equator. Multiple ELFIN passes over the same MLT sector allow us to study the spatial and temporal evolution of the EMIC wave - electron interaction region. Case studies employing conjugate ground-based or equatorial observations of the EMIC waves reveal that the energy of moderate and strong precipitation at ELFIN approximately agrees with theoretical expectations for cyclotron resonant interactions in a cold plasma. Using multiple years of ELFIN data uniformly distributed in local time, we assemble a statistical database of ~50 events of strong EMIC wave-driven precipitation. Most reside at L ~ 5 – 7 at dusk, while a smaller subset exists at L ~ 8 – 12 at post-midnight. The energies of the peak-precipitation ratio and of the half-peak precipitation ratio (our proxy for the minimum resonance energy) exhibit an L-shell dependence in good agreement with theoretical estimates based on prior statistical observations of EMIC wave power spectra. The precipitation ratio’s spectral shape for the most intense events has an exponential falloff away from the peak (i.e., on either side of ~ 1.45 MeV). It too agrees well with quasi-linear diffusion theory based on prior statistics of wave spectra. It should be noted though that this diffusive treatment likely includes effects from nonlinear resonant interactions (especially at high energies) and nonresonant effects from sharp wave packet edges (at low energies). Sub-MeV electron precipitation observed concurrently with strong EMIC wave-driven >1 MeV precipitation has a spectral shape that is consistent with efficient pitch-angle scattering down to ~ 200-300 keV by much less intense higher frequency EMIC waves at dusk (where such waves are most frequent). At ~100 keV, whistler-mode chorus may be implicated in concurrent precipitation. These results confirm the critical role of EMIC waves in driving relativistic electron losses. Nonlinear effects may abound and require further investigation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Models and Algorithms for Equilibrium Analysis of Mixed-Material Nucleic Acid Systems

Dynamic programming algorithms within the NUPACK software suite enable analysis of equilibrium base-pairing properties for complex and test tube ensembles containing arbitrary numbers of interacting nucleic acid strands. Currently, calculations are limited to single-material systems that are either all-RNA or all-DNA. Here, to enable analysis of mixed-material systems that are critical for modern applications in vitro, in situ, and in vivo, we develop physical models and dynamic programming algorithms that allow the material of the system to be specified at nucleotide resolution. Free energy parameter sets are constructed for both RNA/DNA and RNA/2'OMe-RNA mixed-material systems by combining available empirical mixed-material parameters with single-material parameter sets to enable treatment of the full complex and test tube ensembles. New dynamic programming recursions account for the material of each nucleotide throughout the recursive process. For a complex with N nucleotides, the mixed-material dynamic programming algorithms maintain the O(N 3 ) time complexity of the single-material algorithms, enabling efficient calculation of diverse physical quantities over complex and test tube ensembles (e.g., complex partition function, equilibrium complex concentrations, equilibrium base-pairing probabilities, minimum free energy secondary structure(s), and Boltzmann-sampled secondary structures) at a cost increase of roughly 2.0-3.5×. The results of existing single-material algorithms are exactly reproduced when applying the new mixed-material algorithms to single-material systems. Accuracy is significantly enhanced using mixed-material models and algorithms to predict RNA/DNA and RNA/2'OMe-RNA duplex melting temperatures from the experimental literature as well as RNA/DNA melt profiles from new experiments. In conclusion, mixed-material analyses can be performed online using the NUPACK web app (www.nupack.org) or locally using the NUPACK Python module.

2′OMe-RNA↗

DOE OSTI2022

Diagnosing the pedestal magnetic field and magnetohydrodynamics radial structure with pedestal–scrape of layer electron cyclotron emission radiation inversion in H-mode plasma (invited)

Forward modeling is used to interpret inversion patterns of the pedestal–Scrape of Layer (SOL) Electron Cyclotron Emission (ECE) in DIII-D H-mode experiments. The modeling not only significantly improves the ECE data interpretation quality but also leads to the potential measurements of (1) the magnetic field strength |B| at the separatrix, (2) the pedestal |B| evolution during an inter-Edge Localized Mode (ELM) period, and (3) the pedestal Magnetohydrodynamics (MHD) radial structure. The ECE shine-through effect leads to three types of pedestal–SOL radiation inversions that are discussed in this paper. The first type of inversion is the non-monotonic T e,rad profile with respect to the major radius. Using the ECE frequency at the minimum T e,rad , the inversion can be applied to measure the magnetic field |B| at the separatrix and calibrate the mapping of the ECE channels with respect to the separatrix. The second type of inversion refers to the opposite phase between the radiation fluctuations δT e,rad at the pedestal and SOL. This δT e,rad phase inversion is sensitive to density and temperature fluctuations at the pedestal foot and, thus, can be used to qualitatively measure the MHD radial structure. The third type of inversion appears when the pedestal and SOL T e,rad evolve in an opposite trend, which can be used to infer the pedestal |B| field change during an inter-ELM period. Finally, the bandwidth effect on measuring δT e,rad due to pedestal MHD is also investigated in the radiation modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Noncollinear magnetic configurations and substrate-mediated interactions in Mn trimers on the GaN (000$\bar{1}$) surface

Collinear and noncollinear calculations based on density functional theory are carried out to elucidate the magnetic ordering of Mn trimers on a GaN ($000\bar{1}$) substrate. These trimers had previously been observed in $3a \times 3a$ surface reconstructions through Mn deposition onto the N-polar face of wurtzite GaN($000\bar{1}$). In this work, we start off by studying the effect of spin orbit coupling for the case of monomers and dimers of Mn atoms on top of a GaN surface. Based on an effective spin Hamiltonian, we estimate the magnetic anisotropy energy (MAE) for those cases and found that it is four orders of magnitude weaker than the exchange magnetic coupling between Mn adatoms. In the Mn trimer case, the magnetic ground state has the Mn spins in-plane with the GaN surface in which the relative spin orientation within each trimer is noncollinear due to the competition between the two antiferromagnetic interactions that affect each Mn spin in the trimer, which leads to the found energy minimum with 120 degree angles between the spins. By exploring the nature and fundamental mechanisms for the magnetic interaction among the Mn trimers, we find that the surface states of the substrate play a key role, involving a Ruderman-Kittel-Kasuya-Yosida (RKKY)-type interaction. We report on an electron-mediated long-distance exchange coupling between localized magnetic moments on a GaN($000\bar{1}$) surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗