Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Internal state variables”

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 469 records · Page 26

Nonlinearity of the post-spinel transition and its expression in slabs and plumes worldwide

Phase transitions in the mantle control its internal dynamics and structure. The post-spinel transition marks the upper–lower mantle boundary, where ringwoodite dissociates into bridgmanite plus ferropericlase, and its Clapeyron slope regulates mantle flow across it. This interaction has previously been assumed to have no lateral spatial variations, based on the assumption of a linear post-spinel boundary in pressure and temperature. Here we present laser-heated diamond anvil cell experiments with synchrotron X-ray diffraction to better constrain this boundary, especially at higher temperatures. Combining our data with results from the literature, and using a global analysis based on machine learning, we find a pronounced nonlinearity in the post-spinel boundary, with its slope ranging from –4 MPa/K at 2100 K, to –2 MPa/K at 1950 K, and to 0 MPa/K at 1600 K. Changes in temperature over time and space can therefore cause the post-spinel transition to have variable effects on mantle convection and the movement of subducting slabs and upwelling plumes.

58 GEOSCIENCES↗

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasizes the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M qi) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in turns of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, Derck↗

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasises the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters, but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M (dot) q(sub i)) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in terms of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, D.↗

The Global Ocean Observing System

A Global Ocean Observing System (GOOS) should be established now with international coordination (1) to address issues of global change, (2) to implement operational ENSO forecasts, (3) to provide the data required to apply global ocean circulation models, and (4) to extract the greatest value from the one billion dollar investment over the next ten years in ocean remote sensing by the world's space agencies. The objectives of GOOS will focus on climatic and oceanic predictions, on assessing coastal pollution, and in determining the sustainability of living marine resources and ecosystems. GOOS will be a complete system including satellite observations, in situ observations, numerical modeling of ocean processes, and data exchange and management. A series of practical and economic benefits will be derived from the information generated by GOOS. In addition to the marine science community, these benefits will be realized by the energy industries of the world, and by the world's fisheries. The basic oceanic variables that are required to meet the oceanic and predictability objectives of GOOS include wind velocity over the ocean, sea surface temperature and salinity, oceanic profiles of temperature and salinity, surface current, sea level, the extent and thickness of sea ice, the partial pressure of CO2 in surface waters, and the chlorophyll concentration of surface waters. Ocean circulation models and coupled ocean-atmosphere models can be used to evaluate observing system design, to assimilate diverse data sets from in situ and remotely sensed observations, and ultimately to predict future states of the system. The volume of ocean data will increase enormously over the next decade as new satellite systems are launched and as complementary in situ measuring systems are deployed. These data must be transmitted, quality controlled, exchanged, analyzed, and archived with the best state-of-the-art computational methods.

Kester, Dana↗

Modeled Results of Four Residential Energy Efficiency Measure Packages for Deriving Advanced Building Construction Research Targets

The Advanced Building Construction (ABC) Initiative from the U.S. Department of Energy Building Technologies Office is working to accelerate industrialized construction innovations for decarbonizing buildings. To inform performance and cost targets for research under the ABC Initiative, this analysis used the ResStock™ tool to evaluate the energy savings, utility bill impacts, and carbon emissions impacts of four simulated upgrade packages with specific target performance levels on a large sample of residential dwelling units (approximately 550,000) representative of the U.S. housing stock.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Differences in cluster and internal wake effects from mesoscale and large-eddy simulations off the US East Coast

Mesoscale simulations are increasingly used to estimate wake effects within and between large wind farms, despite limited validation for large-scale wake effects. This study evaluates the capabilities and limitations of mesoscale simulations in capturing wake-induced impacts on wind turbine power production through a direct comparison with large-domain large-eddy simulations (LESs) for three planned offshore wind farms under realistic atmospheric conditions and a range of atmospheric stabilities. We assess mesoscale performance in replicating wake characteristics behind single and multiple turbine clusters and quantify the resulting variability in mean turbine power. Results show that mesoscale Weather Research and Forecasting simulations with the Fitch wind farm parameterization capture key features of the velocity deficit downstream of both single and multiple wind farms, with mean root-mean-square errors near 5 % and good agreement with stability-driven wake behavior. However, in these simulations, the mesoscale Fitch parameterization underestimates power losses from internal wake effects, particularly when turbines align with the prevailing wind direction or under stable stratification. In these conditions, individual wakes persist and dominate downstream power deficits. The coarse resolution of the mesoscale simulations limits their ability to resolve individual wind turbine wakes that drive power fluctuations within wind farms. Nonetheless, mesoscale simulations can yield accurate estimates of combined wake losses from internal and cluster effects across some wind direction sectors, where errors in wake representation may cancel each other out. These findings underscore the strengths of mesoscale simulations for capturing broader wake patterns while highlighting their limitations for modeling turbine-level power losses. Future work should explore hybrid modeling approaches to capture both long-range cluster wake propagation and localized internal wake dynamics.

17 WIND ENERGY↗

A high resolution, gridded product for vapor pressure deficit using Daymet

Vapor pressure deficit (VPD) is a critical variable in assessing drought conditions and evaluating plant water stress. Gridded products of global and regional VPD are not freely available from satellite remote sensing, model reanalysis, or ground observation datasets. We present two versions of the first gridded VPD product for the Continental US and parts of Northern Mexico and Southern Canada (CONUS+) at a 1 km spatial resolution and daily time step. We derived VPD from Daymet maximum daily temperature and average daily vapor pressure and scale the estimates based on (1) climate determined by the Köppen-Geiger classifications and (2) land cover determined by the International Geosphere-Biosphere Programme. Ground-based VPD data from 253 AmeriFlux sites representing different climate and land cover classifications were used to improve the Daymet-derived VPD estimates for every pixel in the CONUS+ grid to produce the final datasets. We evaluated the Daymet-derived VPD against independent observations and reanalysis data. The CONUS+ VPD datasets will aid in investigating disturbances including drought and wildfire, and informing land management strategies.

54 ENVIRONMENTAL SCIENCES↗

One‐at‐a‐Time Parameter Perturbation Ensemble of the Community Land Model, Version 5.1

Comprehensive land models are subject to significant parametric uncertainty, which can be hard to quantify due to the large number of parameters and high model computational costs. We constructed a large parameter perturbation ensemble (PPE) for the Community Land Model version 5.1 with biogeochemistry configuration (CLM5.1-BGC). We performed more than 2,000 simulations perturbing 211 parameters across six forcing scenarios. This provides an expansive data set, which can be used to identify the most influential parameters on a wide range of output variables globally, by biome, or by plant functional type. We found that parameter effects can exceed scenario effects and that a small number of parameters explains a large fraction of variance across our ensemble. The most important parameters can differ regionally and also based on the forcing scenario. The software infrastructure developed for this experiment has greatly reduced the human and computer time needed for CLM PPEs, which can facilitate routine investigation of parameter sensitivity and uncertainty, as well as automated calibration.

Kennedy, Daniel [NSF National Center for Atmospher↗

Standardising the “Gregory method” for calculating equilibrium climate sensitivity

The equilibrium climate sensitivity (ECS) – the equilibrium global mean temperature response to a doubling of atmospheric CO 2 – is a high-profile metric for quantifying the Earth system's response to human-induced climate change. A widely applied approach to estimating the ECS is the “Gregory method” (Gregory et al., 2004), which uses an ordinary least squares (OLS) regression between the net radiative flux, N, and surface air temperature anomalies, ΔT, from a 150 year experiment in which atmospheric CO 2 concentrations are quadrupled. The ECS is determined by extrapolating the linear fit to N=0, i.e. the ΔT-intercept, indicating the point at which the system is back in equilibrium. This method has been used to compare ECS estimates across the CMIP5 and CMIP6 ensembles and will likely be a key diagnostic for CMIP7. Despite its widespread application, there is little consistency or transparency between studies in how the climate model data is processed prior to the regression, leading to potential discrepancies in ECS estimates. We identify 32 alternative data processing pathways, varying by differences in global mean weighting, net radiative flux variable, anomaly calculation method, and linear regression fit. Using 44 CMIP6 models, we systematically assess the impact of these choices on ECS estimates and calculate uncertainty ranges using two bootstrap approaches. While the inter-model ECS range is insensitive to the data processing pathway, individual outlier models exhibit notable differences. Approximating a model's native grid cell area (if irregular) with cosine of the latitude can decrease the ECS by 11 %, the choice of N-variable can change the ECS by 6 %, and some anomaly calculation methods can introduce spurious temporal correlations in the processed data. Beyond data processing choices, we also evaluate an alternative linear regression method – total least squares (TLS) – which has a more statistically robust basis than OLS. However, for consistency with previous literature, and given TLS may reduce the ECS compared to OLS (by up to 24 %), thereby making a known bias in the Gregory method worse, we do not feel there is sufficient clarity to recommend a transition to TLS in all cases. To improve reproducibility and comparability in future studies, we recommend a standardised Gregory method: weighting the global mean by cell area, using the top of the atmosphere (as opposed to the top of model) N-variable, and calculating anomalies by first applying a rolling average to the preindustrial control timeseries then subtracting from the raw CO 2 quadrupling experiment. This approach accounts for model drift while reducing noise in the data to best meet the pre-conditions of the linear regression. While CMIP6 results of the multi-model mean ECS appear insensitive to these processing choices, similar assumptions may not hold for CMIP7, underscoring the need for standardised data preparation in future climate sensitivity assessments.

Geosciences↗

Providing Satellite-Assisted Particulate Matter (SAPM) Estimates for the Models, In situ, and Remote sensing of Aerosols (MIRA) Working Group

The Models, In situ, and Remote sensing of Aerosols (MIRA) Working Group is an international forum that promotes collaboration amongst researchers of these three atmospheric aerosol communities. MIRA currently consists of a collection of five interdisciplinary and independently funded Topic groups with focused goals and generally involve requests for additional scientific datasets. Participants of the MIRA Satellite-Assisted Particulate Matter (SAPM) Topic group study particulate matter with diameters smaller than 2.5 microns (PM2.5), because of its substantial role in air pollution and the resulting negative impacts on human health. SAPM seeks to provide intercomparisons of different methods and techniques for obtaining surface PM2.5 proxies using spaceborne passive and active remote sensors, aerosol models, and in situ observations. The ultimate motivation of SAPM is to provide greater PM2.5 coverage than in situ ground stations, which are limited in some areas of the contiguous United States and large regions throughout the globe. In this poster presentation, we provide an overview of the PM2.5 estimation techniques applied by current SAPM team members, as well as the strengths and limitations of each technique. These approaches include spaceborne lidar (Cloud-Aerosol Lidar with Orthogonal Polarization) only, combined spaceborne lidar (Cloud Aerosol Transport System) and a model (Goddard Earth Observing System), and using aerosol optical depth (AOD) from a spaceborne passive remote sensor (Moderate Resolution Imaging Spectroradiometer) to explore city trends in AOD and relate those to city trends in PM2.5. Also in this poster presentation, we highlight results of our SAPM study of India, a country characterized by high levels of PM2.5 concentrations. We show the temporal and spatial variability of PM2.5 over India derived using various approaches and validation with ground-based in situ PM2.5 observations. The SAPM Topic group actively seeks additional collaborators, including those working with in situ aerosol measurements, and is interested in aerosol datasets to improve and/or validate the PM2.5 estimates.

Travis D Toth↗

Indicators of Global Climate Change 2023: annual update of key indicators of the state of the climate system and human influence

Intergovernmental Panel on Climate Change (IPCC) assessments are the trusted source of scientific evidence for climate negotiations taking place under the United Nations Framework Convention on Climate Change (UNFCCC). Evidence-based decision-making needs to be informed by up-to-date and timely information on key indicators of the state of the climate system and of the human influence on the global climate system. However, successive IPCC reports are published at intervals of 5–10 years, creating potential for an information gap between report cycles. We follow methods as close as possible to those used in the IPCC Sixth Assessment Report (AR6) Working Group One (WGI) report. We compile monitoring datasets to produce estimates for key climate indicators related to forcing of the climate system: emissions of greenhouse gases and short-lived climate forcers, greenhouse gas concentrations, radiative forcing, the Earth's energy imbalance, surface temperature changes, warming attributed to human activities, the remaining carbon budget, and estimates of global temperature extremes. The purpose of this effort, grounded in an open-data, open-science approach, is to make annually updated reliable global climate indicators available in the public domain. As they are traceable to IPCC report methods, they can be trusted by all parties involved in UNFCCC negotiations and help convey wider understanding of the latest knowledge of the climate system and its direction of travel. The indicators show that, for the 2014–2023 decade average, observed warming was 1.19 [1.06 to 1.30] °C, of which 1.19 [1.0 to 1.4] °C was human-induced. For the single-year average, human-induced warming reached 1.31 [1.1 to 1.7] °C in 2023 relative to 1850–1900. The best estimate is below the 2023-observed warming record of 1.43 [1.32 to 1.53] °C, indicating a substantial contribution of internal variability in the 2023 record. Human-induced warming has been increasing at a rate that is unprecedented in the instrumental record, reaching 0.26 [0.2–0.4] °C per decade over 2014–2023. This high rate of warming is caused by a combination of net greenhouse gas emissions being at a persistent high of 53±5.4 Gt CO 2 e yr -1 over the last decade, as well as reductions in the strength of aerosol cooling. Despite this, there is evidence that the rate of increase in CO 2 emissions over the last decade has slowed compared to the 2000s, and depending on societal choices, a continued series of these annual updates over the critical 2020s decade could track a change of direction for some of the indicators presented here.

54 ENVIRONMENTAL SCIENCES↗

Spatial Study 2021: Sample-Based Surface Water Chemistry and Organic Matter Characterization across Watersheds in the Yakima River Basin, Washington, USA (v3)

This dataset supports a broader study examining the drivers of spatial variability in sediment respiration rates in the Yakima River Basin. The dataset provides geochemistry and organic matter characterization data generated from samples collected during the same two-week period at 47 sites within multiple rivers throughout the Yakima River Basin in Washington, USA. Related sensor data are published at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1892054. This data package was originally published September 2022. It was updated May 2023 (modified files) and November 2024 (new and modified files). See the change history section in the readme for more details. This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) dissolved inorganic carbon (DIC), dissolved organic carbon (DOC; reported as non-purgeable organic carbon; NPOC), total nitrogen (TN), total suspended solids (TSS), ions, and benzene polycarboxylic acid (BPCA) concentration and stable isotope data; (5) averaged values from water chemistry data; (6) surface water sampling protocol; (7) sensor protocol (8) readme; (9) methods codes; (10) international generic-sample number (IGSN) mapping file; and (11) folder of high resolution characterization of organic matter via 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). This folder contains two subfolders, one containing the .xml data files and the other containing instructions for using Formultitude (https://github.com/PNNL-Comp-Mass-Spec/Formultitude) and an R script to process the data based on the user's specific needs. All files are .csv, .pdf, .R, .ref, or .xml. We thank the United States Forest Service, Washington Department of Natural Resources, Washington Department of Fish and Wildlife, Washington State Parks, Confederated Tribes and Bands of the Yakama Nation, and Cowiche Canyon Conservancy for access to field locations where these samples were collected. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗

Radiative absorption enhancements by black carbon controlled by particle-to-particle heterogeneity in composition

Black carbon (BC) absorbs solar radiation, leading to a strong but uncertain warming effect on climate. A key challenge in modeling and quantifying BC’s radiative effect on climate is predicting enhancements in light absorption that result from internal mixing between BC and other aerosol components. Modeling and laboratory studies show that BC, when mixed with other aerosol components, absorbs more strongly than pure, uncoated BC; however, some ambient observations suggest more variable and weaker absorption enhancement. We show that the lower-than-expected enhancements in ambient measurements result from a combination of two factors. First, the often used spherical, concentric core-shell approximation generally overestimates the absorption by BC. Second, and more importantly, inadequate consideration of heterogeneity in particle-to-particle composition engenders substantial overestimation in absorption by the total particle population, with greater heterogeneity associated with larger model–measurement differences. We show that accounting for these two effects—variability in per-particle composition and deviations from the core-shell approximation—reconciles absorption enhancement predictions with laboratory and field observations and resolves the apparent discrepancy. Furthermore, our consistent model framework provides a path forward for improving predictions of BC’s radiative effect on climate.

54 ENVIRONMENTAL SCIENCES↗

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗

Application of Simulated Annealing and Related Algorithms to TWTA Design

Simulated Annealing (SA) is a stochastic optimization algorithm used to search for global minima in complex design surfaces where exhaustive searches are not computationally feasible. The algorithm is derived by simulating the annealing process, whereby a solid is heated to a liquid state and then cooled slowly to reach thermodynamic equilibrium at each temperature. The idea is that atoms in the solid continually bond and re-bond at various quantum energy levels, and with sufficient cooling time they will rearrange at the minimum energy state to form a perfect crystal. The distribution of energy levels is given by the Boltzmann distribution: as temperature drops, the probability of the presence of high-energy bonds decreases. In searching for an optimal design, local minima and discontinuities are often present in a design surface. SA presents a distinct advantage over other optimization algorithms in its ability to escape from these local minima. Just as high-energy atomic configurations are visited in the actual annealing process in order to eventually reach the minimum energy state, in SA highly non-optimal configurations are visited in order to find otherwise inaccessible global minima. The SA algorithm produces a Markov chain of points in the design space at each temperature, with a monotonically decreasing temperature. A random point is started upon, and the objective function is evaluated at that point. A stochastic perturbation is then made to the parameters of the point to arrive at a proposed new point in the design space, at which the objection function is evaluated as well. If the change in objective function values (Delta)E is negative, the proposed new point is accepted. If (Delta)E is positive, the proposed new point is accepted according to the Metropolis criterion: rho((Delta)f) = exp((-Delta)E/T), where T is the temperature for the current Markov chain. The process then repeats for the remainder of the Markov chain, after which the temperature is decremented and the process repeats. Eventually (and hopefully), a near-globally optimal solution is attained as T approaches zero. Several exciting variants of SA have recently emerged, including Discrete-State Simulated Annealing (DSSA) and Simulated Tempering (ST). The DSSA algorithm takes the thermodynamic analogy one step further by categorizing objective function evaluations into discrete states. In doing so, many of the case-specific problems associated with fine-tuning the SA algorithm can be avoided; for example, theoretical approximations for the initial and final temperature can be derived independently of the case. In this manner, DSSA provides a scheme that is more robust with respect to widely differing design surfaces. ST differs from SA in that the temperature T becomes an additional random variable in the optimization. The system is also kept in equilibrium as the temperature changes, as opposed to the system being driven out of equilibrium as temperature changes in SA. ST is designed to overcome obstacles in design surfaces where numerous local minima are separated by high barriers. These algorithms are incorporated into the optimal design of the traveling-wave tube amplifier (TWTA). The area under scrutiny is the collector, in which it would be ideal to use negative potential to decelerate the spent electron beam to zero kinetic energy just as it reaches the collector surface. In reality this is not plausible due to a number of physical limitations, including repulsion and differing levels of kinetic energy among individual electrons. Instead, the collector is designed with multiple stages depressed below ground potential. The design of this multiple-stage collector is the optimization problem of interest. One remaining problem in SA and DSSA is the difficulty in determining when equilibrium has been reached so that the current Markov chain can be terminated. It has been suggested in recent literature that simulating the thermodynamic properties opecific heat, entropy, and internal energy from the Boltzmann distribution can provide good indicators of having reached equilibrium at a certain temperature. These properties are tested for their efficacy and implemented in SA and DSSA code with respect to TWTA collector optimization.

Radke, Eric M.↗

Fast and Invertible Simplicial Approximation of Magnetic‐Following Interpolation for Visualizing Fusion Plasma Simulation Data

We introduce a fast and invertible approximation for fusion plasma simulation data represented as 2D planar meshes with connectivities approximating magnetic field lines along the toroidal dimension in deformed 3D toroidal spaces. Scientific variables (e.g., density and temperature) in these fusion data are interpolated following a complex magnetic-field-line-following scheme in the toroidal space represented by a cylindrical coordinate system. This deformation in the 3D space poses challenges for root-finding and interpolation. To this end, we propose a novel paradigm for visualizing and analyzing such data based on a newly developed algorithm for constructing a 3D simplicial mesh within the deformed 3D space. Our algorithm generates a tetrahedral mesh that connects the 2D meshes using tetrahedra while adhering to the constraints on node connectivities imposed by the magnetic field-line scheme. Specifically, we first divide the space into smaller partitions to reduce complexity based on the input geometries and constraints on connectivities. Then, we independently search for a feasible tetrahedralization of each partition, considering nonconvexity. We demonstrate our method with two X-Point Gyrokinetic Code (XGC) simulation datasets on the International Thermonuclear Experimental Reactor (ITER) and Wendelstein 7-X (W7-X), and use an ocean simulation dataset to substantiate broader applicability of our method. An open source implementation of our algorithm is available at https://github.com/rcrcarissa/DeformedSpaceTet.

Ren, Congrong [The Ohio State Univ., Columbus, OH ↗

Towards provision of regularly updated climate data from the Coupled Model Intercomparison Project

The Coupled Model Intercomparison Project (CMIP) is a flagship of the World Climate Research Programme (WCRP). CMIP has become a recognised ‘brand’ in climate circles evolving over the last thirty years from a targeted research activity by a small number of climate modelling centres intercomparing their Earth System Model (ESM) simulations to a broad international coordinated research effort (Durack et al, 2025). CMIP is organized as a research activity leveraging funded and in-kind contributions from experts within modelling centres and the broader scientific community supported more recently by a fully-funded International Project Office. Within CMIP, Model Intercomparison Projects (MIPs) are community-designed to understand past, present and future climate. CMIP data provides a valuable resource for climate research and is routinely used to assess model representation of climate processes and test scientific hypotheses in the context of model uncertainty and (forced and internal) variability as evident from its prolific use in scientific publications1 . The impact relies on enabling infrastructure (most prominently via the Earth System Grid Federation (ESGF)), which allows sharing of simulation output, provision of the boundary conditions used in each simulation, and definition of the data standards that are essential to facilitating wide use of the data. The impact is supplemented by the wide-ranging scrutiny to which model simulations are subjected. Beyond its use in research, CMIP data is a key resource for communities producing derived climate information from downscaling and impact studies, such as the Coordinated Regional Downscaling Experiment (CORDEX; Gutowski et al., 2016) and the Intersectoral Impacts MIP (ISIMIP; Frieler et al., 2024). Government, academic and commercial entities also increasingly rely on CMIP and its downstream data for climate risk assessments and climate services (for example, Copernicus Climate Change Service and World Bank portal). This means that, although CMIP is a research activity, it increasingly serves a secondary and very relevant role as a provider of climate data – a long-recognised dichotomy (Stevens, 2024). Research and applications have distinct needs, with the former requiring flexibility and generality and the latter consistency. Here we explain how the design of the research activity has been adapted to reduce the burdens imposed by applications and how the research infrastructure might evolve to further enable scientific inquiry. We propose one possible approach to consistently providing model information and projections for applications in the future.

Environmental sciences↗

Hybrid Quantum–Classical Graph Transformers for Efficient Sentiment Analysis

Quantum Machine Learning (QML) offers a promising paradigm that leverages quantum computing principles to develop efficient and expressive models for learning from complex and structured data. Recent advances in natural language processing (NLP) and artificial intelligence (AI) have demonstrated capabilities in understanding, generating, and reasoning over linguistic and multimodal information. In this work, we present the Quantum Graph Transformer (QGT), a hybrid quantum–classical architecture that extends graph transformer capabilities through quantum self-attention. The QGT models variable-length sentences as token graphs, where both the embedding encoding and the self-attention mechanisms are implemented using parameterized quantum circuits (PQCs), enabling efficient contextual learning with significantly fewer trainable parameters. We train QGT using both fully connected and 𝑘 -nearest-neighbor graph structures and evaluate it on five benchmark sentiment-classification datasets. Experimental results show that QGT consistently achieves higher or comparable accuracy to existing quantum NLP models and outperforms a Classical Graph Transformer (CGT) baseline with identical architecture, achieving 29.4 × fewer parameters while requiring 3–5 × fewer samples to reach comparable performance. These findings highlight the potential of graph-based quantum models as scalable and data-efficient architectures for natural language understanding.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗