Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global optimization”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An adaptive stochastic sequential quadratic programming with differentiable exact augmented lagrangians

In this study, we consider solving nonlinear optimization problems with a stochastic objective and deterministic equality constraints. We assume for the objective that its evaluation, gradient, and Hessian are inaccessible, while one can compute their stochastic estimates by, for example, subsampling. We propose a stochastic algorithm based on sequential quadratic programming (SQP) that uses a differentiable exact augmented Lagrangian as the merit function. To motivate our algorithm design, we first revisit and simplify an old SQP method Lucidi developed for solving deterministic problems, which serves as the skeleton of our stochastic algorithm. Based on the simplified deterministic algorithm, we then propose a non-adaptive SQP for dealing with stochastic objective, where the gradient and Hessian are replaced by stochastic estimates but the stepsizes are deterministic and prespecified. Finally, we incorporate a recent stochastic line search procedure Paquette and Scheinberg into the non-adaptive stochastic SQP to adaptively select the random stepsizes, which leads to an adaptive stochastic SQP. The global "almost sure" convergence for both non-adaptive and adaptive SQP methods is established. Numerical experiments on nonlinear problems in CUTEst test set demonstrate the superiority of the adaptive algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

High-pressure ignition delay time measurements of a four-component gasoline surrogate and its high-level blends with ethanol and methyl acetate

Ignition delay times of a gasoline surrogate (iso-octane/toluene/n-heptane/1-hexene at 55/25/15/5% by liquid volume, developed in McCormick et al. 2017) with high-level bioblendstock-gasoline surrogate blends (50% and 85% biofuel by liquid volume of each ethanol and methyl acetate) were collected behind reflected shock waves. Post-reflected-shock temperatures ranged from 968 to 1361 K at pressures of about 4, 10, 25, and 50 atm. Data were collected for real fuel–air mixtures at fuel lean and stoichiometric conditions (φ = 0.5, 1.0) with focus on the more-practical, fuel-lean conditions. Ignition delay times were measured from OH* chemiluminescence around 307 nm using an endwall diagnostic. The RON and MON of the surrogate were 90.3 and 84.7, respectively. Ethanol was chosen due to its wide use in flex-fuel vehicles and availability, while also providing an increased octane rating. Methyl acetate was chosen for its especially high octane rating to investigate the effect of an extreme case. To validate the gasoline surrogate, the data are compared to real gasoline (RD387) and several gasoline surrogate experiments from the literature. Using the wide range of pressures studied, the gasoline surrogate’s pressure dependence was quantified to account for test-to-test variations using regression analysis. Similarly, a global correlation for gasoline and its surrogates was developed using all available data from the literature. Two modern chemical kinetics models targeting gasoline and its surrogates are compared to the ignition delay time measurements. These new high-pressure, high-bioblendstock concentration tests provide required chemical kinetic data for optimizing fuel and engine design.

42 ENGINEERING↗
DOE OSTI2022

Biogeochemical Equilibrium Responses to Maximal Productivity in High Nutrient Low Chlorophyll Regions

Whether ocean iron fertilization in high nutrient low chlorophyll (HNLC) regions is a practical geoengineering strategy to combat climate warming has been debated because the fate of absorbed carbon by the ocean remains unclear. Here, we use an optimized oceanic biogeochemical inverse model to explore the equilibrium responses to increased productivity in major HNLC regions. We find that an overall downward shift of nutrients and carbon from the surface and intermediate to the deep waters with maximal productivity in the HNLC Southern Ocean (SO), North Pacific (NP) and eastern equatorial Pacific (EP). The efficiency of soft tissue pump is increased by 37.8%, 2.7% and 5.4% in the SO, EP and NP simulation while atmospheric CO 2 concentration is drawn down by 45.1, 3.4 and 6.7 ppm, respectively. Meanwhile, global mean oxygen is decreased by 21.1%, 1.3% and 2.4% in the SO, EP and NP simulation, accompanied by an expansion of oxygen minimum zones by 400%, 26% and 63%, respectively. Our study indicates that iron fertilization has the potential to be an effective geoengineering approach to mitigate climate warming but also has a strong impact on the development of ocean hypoxia.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Modeling Ring Current Proton Fluxes Using Artificial Neural Network and Van Allen Probe Measurements

Abstract Terrestrial ring current dynamics are a critical part of the near‐space environment, in that they directly drive geomagnetic field variations that control particle drifts, and define geomagnetic storms. The present study aims to specify a global and time‐varying distribution of ring current proton using geomagnetic indices and solar wind parameters with their history as input. We train an artificial neural network (ANN) model to reproduce proton fluxes measured by the Radiation Belt Storm Probes Ion Composition Experiment instrument onboard Van Allen Probes. By choosing optimal feature parameters and their history length, the model results show a high correlation and a small error between model specifications and satellite measurements. The modeled results well capture energy‐dependent proton dynamics in association with geomagnetic storms, including inward radial diffusion, acceleration and decay. Our ANN model produces proton fluxes with their corresponding 3D spatiotemporal variations, capturing the latitudinal distribution and local time asymmetry that are consistent with observations and that can further inform theory.

Li, Jinxing↗
DOE OSTI2024

Multi-omics analysis reveals the dynamic interplay between Vero host chromatin structure and function during vaccinia virus infection

The genome folds into complex configurations and structures thought to profoundly impact its function. The intricacies of this dynamic structure-function relationship are not well understood particularly in the context of viral infection. To unravel this interplay, here we provide a comprehensive investigation of simultaneous host chromatin structural (via Hi-C and ATAC-seq) and functional changes (via RNA-seq) in response to vaccinia virus infection. Over time, infection significantly impacts global and local chromatin structure by increasing long-range intra-chromosomal interactions and B compartmentalization and by decreasing chromatin accessibility and inter-chromosomal interactions. Local accessibility changes are independent of broad-scale chromatin compartment exchange (~12% of the genome), underscoring potential independent mechanisms for global and local chromatin reorganization. While infection structurally condenses the host genome, there is nearly equal bidirectional differential gene expression. Despite global weakening of intra-TAD interactions, functional changes including downregulated immunity genes are associated with alterations in local accessibility and loop domain restructuring. Therefore, chromatin accessibility and local structure profiling provide impactful predictions for host responses and may improve development of efficacious anti-viral counter measures including the optimization of vaccine design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

xesn: Echo state networks powered by Xarray and Dask

Xesn is a Python package that allows scientists to easily design Echo State Networks (ESNs) for forecasting problems. ESNs are a Recurrent Neural Network architecture introduced by Jaeger (2001) that are part of a class of techniques termed Reservoir Computing. One defining characteristic of these techniques is that all internal weights are determined by a handful of global, scalar parameters, thereby avoiding problems during backpropagation and reducing training time significantly. Because this architecture is conceptually simple, many scientists implement ESNs from scratch, leading to questions about computational performance. Xesn offers a straightforward, standard implementation of ESNs that operates efficiently on CPU and GPU hardware. The package leverages optimization tools to automate the parameter selection process, so that scientists can reduce the time finding a good architecture and focus on using ESNs for their domain application. Importantly, the package flexibly handles forecasting tasks for out-of-core, multi-dimensional datasets, eliminating the need to write parallel programming code. Xesn was initially developed to handle the problem of forecasting weather dynamics, and so it integrates naturally with Python packages that have become familiar to weather and climate scientists such as Xarray (Hoyer & Hamman, 2017). However, the software is ultimately general enough to be utilized in other domains where ESNs have been useful, such as in signal processing (Jaeger & Haas, 2004).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

National Smart Manufacturing Strategic Plan: To Facilitate More Rapid Development, Deployment and Adoption of Smart Manufacturing Technologies

Smart manufacturing technologies provide real-time data and insight to improve the productivity, efficiency, and competitiveness of U.S. manufacturing, creating the potential for new jobs in the manufacturing sector. These technologies can support U.S. manufacturers’ ability to increase throughput and energy efficiency, and decrease waste, defects, and costs. A wide range of manufacturing and industrial subsectors, particularly energy-intensive and energy-dependent industries, have the potential to benefit from smart manufacturing technologies. There are technical improvements still needed to reduce the costs and barriers (trained workforce and upskilling, software-hardware integration, cost, and technical barriers to deployment of advanced sensors, computing, and communication technologies for existing manufacturing assets) to the adoption of smart manufacturing technologies, especially to small and medium-sized manufacturers, and subsequently, increase the overall adoption rate of smart manufacturing technologies. This report outlines the Department of Energy’s (DOE) strategic plan to accelerate the development and implementation of smart manufacturing technologies in the United States and the actions DOE has taken to use these technologies in smart manufacturing for the United States. DOE’s plan to facilitate more rapid development, deployment and widespread adoption of smart manufacturing technologies derive from the 2018 Strategy for American Leadership in Advanced Manufacturing. The strategy outlines a vision for America’s leadership in advanced manufacturing including the development of intelligent manufacturing systems to optimize manufacturing facilities and support the manufacturing of transformative materials. DOE’s Clean Energy Smart Manufacturing Innovation Institute (CESMII), a Manufacturing USA Institute, focuses on accelerating the development and adoption of advanced sensors, controls, platforms, and models needed for smart manufacturing. The objective is to enhance U.S. manufacturing productivity and global competitiveness through the research and development of these technologies. Smart manufacturing has the potential to improve the overall performance, energy productivity, and efficiency of manufacturing, fostering the economic competitiveness of the manufacturing sector.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Environment-sensitivity functions for gross primary productivity in light use efficiency models

The sensitivity of photosynthesis to environmental changes is essential for understanding carbon cycle responses to global climate change and for the development of modeling approaches that explains its spatial and temporal variability. We collected a large variety of published sensitivity functions of gross primary productivity (GPP) to different forcing variables to assess the response of GPP to environmental factors. These include the responses of GPP to temperature; vapor pressure deficit, some of which include the response to atmospheric CO 2 concentrations; soil water availability (W); light intensity; and cloudiness. These functions were combined in a full factorial light use efficiency (LUE) model structure, leading to a collection of 5600 distinct LUE models. Each model was optimized against daily GPP and evapotranspiration fluxes from 196 FLUXNET sites and ranked across sites based on a bootstrap approach. The GPP sensitivity to each environmental factor, including CO 2 fertilization, was shown to be significant, and that none of the previously published model structures performed as well as the best model selected. From daily and weekly to monthly scales, the best model's median Nash-Sutcliffe model efficiency across sites was 0.73, 0.79 and 0.82, respectively, but poorer at annual scales (0.23), emphasizing the common limitation of current models in describing the interannual variability of GPP. Although the best global model did not match the local best model at each site, the selection was robust across ecosystem types. The contribution of light saturation and cloudiness to GPP was observed across all biomes (from 23% to 43%). Temperature and W dominates GPP and LUE but responses of GPP to temperature and W are lagged in cold and arid ecosystems, respectively. The findings of this study provide a foundation towards more robust LUE-based estimates of global GPP and may provide a benchmark for other empirical GPP products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dual-Fuel Ammonia Equivalence Ratio Sweep Data

Ammonia (NH3) has garnered significant interest as an alternative fuel for meeting international emissions reduction mandates in sectors with high weight and distance requirements, such as shipping. Technical barriers and unanswered questions remain on the combustion strategies that can maximize ammonia utilization and minimize emissions. Prior research studies at the US Department of Energy’s Oak Ridge National Laboratory have shown strong performance with NH3 under dual-fuel mode using conventional diesel combustion (CDC) manifold air pressure settings. Diesel airflow was initially used to simplify retrofitting (no turbocharger modification), which resulted in air-fuel equivalence ratios (λ) greater than 1.5. To characterize potential improvements in dual-fuel NH3 combustion performance at richer in-cylinder conditions, a global λ sweep compared the use of early (E-pilot) and late (L-pilot) single diesel injections. The experiments were conducted at 1200 RPM and 12.8 ± 0.2 bar (75 % load), and λ was varied by decreasing the commanded air flow to the engine at greater than 90 % ammonia energy substitution level. A diesel injection timing sweep was conducted for both the injection strategies at fixed λ, and the timing with the lowest engine-out N2O emissions was identified. The results indicated an optimal balance between CO2,eq and thermal efficiency benefits both E-pilot and L-pilot injection strategy cases compared with CDC at a λ of 1.4. The indicated nitrogen-based emissions exhibited a strong correlation to the ratio of CA5–50 and ignition delay for L-pilot, but no apparent trend emerged for the E-pilot injection strategy at the tested boundary conditions. This dataset includes the raw experimental data and documentation of the experiment conditions and methods.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Engineering acid-base active sites on an N-doped biochar catalyst for microwave-assisted biodiesel production from waste cooking oil

Acid-base bifunctional catalysts play a crucial role in converting high-free fatty acid (FFA) feedstocks into biodiesel by enabling the simultaneous esterification and transesterification of triglycerides in a single step, thereby minimizing saponification and streamlining the process. In this study, we report the design and development of a novel N-doped acid-base bifunctional biochar catalyst synthesized from renewable biomass precursors- glucose and chitosan. The catalyst’s acidic sites were introduced via surface functionalization with o-sulfobenzoic acid cyclic anhydride, grafting -SO₃H groups onto the carbon framework, while basic sites were generated from amine functionalities (-NH) derived from chitosan, yielding a robust and truly bifunctional surface. The catalyst was characterized and used for the microwave-assisted transesterification of waste cooking oil, and process optimization using response surface methodology (RSM-CCD) yielding a maximum biodiesel yield of 98.8 ± 0.4% under optimal conditions (methanol-to-oil molar ratio 20:1, catalyst loading 5.7 wt%, temperature 99 °C, and reaction time 46.6 min). The catalyst demonstrated good stability, maintaining a yield of 92.2 ± 0.4% after eight reuse cycles. A Life Cycle Assessment (LCA) conducted for 1 kg of biodiesel production revealed a global warming potential of 0.85 kg CO 2 eq and human toxicity potential of 1.59 kg 1,4-DB eq, with the largest environmental contributions arising from catalyst preparation and reuse. In conclusion, this work highlights a green and circular approach to biodiesel production, combining renewable catalyst design, waste feedstock utilization, and quantitative sustainability assessment to guide future developments in sustainable catalysis and biofuel engineering.

Bifunctional catalyst↗
DOE OSTI2020

Multi-year incubation experiments boost confidence in model projections of long-term soil carbon dynamics

Abstract Global soil organic carbon (SOC) stocks may decline with a warmer climate. However, model projections of changes in SOC due to climate warming depend on microbially-driven processes that are usually parameterized based on laboratory incubations. To assess how lab-scale incubation datasets inform model projections over decades, we optimized five microbially-relevant parameters in the Microbial-ENzyme Decomposition (MEND) model using 16 short-term glucose (6-day), 16 short-term cellulose (30-day) and 16 long-term cellulose (729-day) incubation datasets with soils from forests and grasslands across contrasting soil types. Our analysis identified consistently higher parameter estimates given the short-term versus long-term datasets. Implementing the short-term and long-term parameters, respectively, resulted in SOC loss (–8.2 ± 5.1% or –3.9 ± 2.8%), and minor SOC gain (1.8 ± 1.0%) in response to 5 °C warming, while only the latter is consistent with a meta-analysis of 149 field warming observations (1.6 ± 4.0%). Comparing multiple subsets of cellulose incubations (i.e., 6, 30, 90, 180, 360, 480 and 729-day) revealed comparable projections to the observed long-term SOC changes under warming only on 480- and 729-day. Integrating multi-year datasets of soil incubations (e.g., > 1.5 years) with microbial models can thus achieve more reasonable parameterization of key microbial processes and subsequently boost the accuracy and confidence of long-term SOC projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Towards Scaling Law Analysis For Spatiotemporal Weather Data

Compute-optimal scaling laws are relatively well studied for NLP and CV, where objectives are typically single-step and targets are comparatively homogeneous. Weather forecasting is harder to characterize in the same framework: autoregressive rollouts compound errors over long horizons, outputs couple many physical channels with disparate scales and predictability, and globally pooled test metrics can disagree sharply with per-channel, late-lead behavior implied by short-horizon training. We extend neural scaling analysis for autoregressive weather forecasting from single-step training loss to long rollouts and per-channel metrics. We quantify (1) how prediction error is distributed across channels and how its growth rate evolves with forecast horizon, (2) if power law scaling holds for test error, relative to rollout length when error is pooled globally, and (3) how that fit varies jointly with horizon and channel for parameter, data, and compute-based scaling axes. We find strong cross-channel and cross-horizon heterogeneity: pooled scaling can look favorable while many channels degrade at late leads. We discuss implications for weighted objectives, horizon-aware curricula, and resource allocation across outputs.

Kiefer Jr, Alexander [ORNL] (ORCID:000000025398874↗
DOE OSTI2024

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗
DOE OSTI2025

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗
DOE OSTI2021

Humate's Role in Terrestrial Greenhouse Gas and Water Reduction

Carbon dioxide (CO 2 ) is considered the sole culprit for global warming; however, nitrous oxide (N 2 O), a greenhouse gas (GHG) with approximately 300 times more global warming potential than CO 2 , accounts for 6% of the GHG emissions in the United States. Seventy five percent of N 2 O emissions come from synthetic nitrogen (N) fertilizer usage in the agriculture sector primarily due to excess fertilization. Numerous studies have shown that changes in soil management practices, specifically optimizing N fertilizer use and amending soil with organic and humate materials can reverse soil damage and improve a farmer's or land reclamation company's balance sheet. Soil restoration is internationally recognized as one of the lowest cost GHG abatement opportunities available. Profitability improves in two ways: (1) lower operating costs resulting from lower input costs (water and fertilizer); and (2) increased revenue by participation in emerging GHG offsets markets, and water quality trading markets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The 8RH 2 Process for Producing Clean Hydrogen with Autothermal Reforming and Carbon Capture

The primary objective of this project is to demonstrate the commercial viability and environmental benefits of an 8 Rivers Hydrogen (8RH 2 ) process, a proprietary technology developed to optimally produce high-purity blue H 2 with up to 99% carbon capture. This 8RH 2 process is particularly notable for its integration of a low-energy cryogenic CO 2 separation system that not only maximizes hydrogen recovery but also minimizes the carbon footprint of the operation. This aligns with global environmental goals and supports the United States' push toward energy sustainability and independence.

03 NATURAL GAS↗
DOE OSTI2025

Design and synthesis of biobased superhydrophobic biochar catalyst derived from Citrus sinensis for biodiesel production using inedible oil feedstocks

In an one-pot in situ trans/esterification of low-grade feedstocks, the simultaneous presence of triglycerides (TAGs) and free fatty acids (FFAs) presents a dual catalytic challenge. While TAGs undergo transesterification to form biodiesel and glycerol, FFAs present in the feedstocks are esterified simultaneously, producing water as a by-product. This water can deactivate acid catalysts by interacting with their active sites, reducing catalytic efficiency and reusability. To overcome this, we developed a hydrophobic sulfonic acid-functionalized biochar catalyst derived from Citrus sinensis (orange peel) via sulfonation with H 2 SO 4 and subsequent silylation using hexamethyldisilazane (HMDS). The optimized catalyst exhibited a high surface area (355.09 m 2 g −1 ), sulfur content (5.33 wt%), and strong hydrophobicity (water contact angle: 154°). Compared to non-hydrophobic analogs, it showed enhanced activity and reusability (up to 10 cycles). Using response surface methodology with central composite design, a biodiesel yield of 99.1 ± 0.4% was achieved under optimal conditions. Life cycle assessment was performed to evaluate the environmental impacts of biodiesel production utilising the synthesized catalyst, considering 1000 kg of biodiesel produced as 1 functional unit. In conclusion, the recorded results showed the cumulative abiotic depletion of fossil resources over the entire biodiesel production process as 87 243.423 MJ and global warming potential as 4103.494 kg CO 2 equivalent.

Biochar↗
DOE OSTI2022

Extreme Low-Temperature Stress Affects Nutritional Quality of Amino Acids in Rice

Global climate change has increased the frequency of extreme climate events, and their effects on the nutritional quality, especially on amino acids in rice, have not been quantified. The data from a 3-year low temperature stress (LTS) experiment including two rice varieties (Huaidao 5 and Nanjing 46), seven minimum/maximum temperature levels (one optimal 21/27°C and six LTS levels from 17/23 to 6/12°C), and three LTS durations (3, 6, and 9 days) after flowering, revealed significant interactive effects of LTS at different stages, durations, and temperature levels on the content and accumulation of amino acids. LTS increased rice total amino acid content, while decreasing its accumulation, with higher sensitivities to LTS at the flowering stage than at the grain filling stage. In most treatments, the lysine (the first limiting amino acid) and phenylalanine content were increased under LTS at early and peak flowering stages but decreased at the grain filling stage in both varieties, and only leucine content was increased at all three stages after flowering, while the content of other essential amino acids differed among the two varieties. With an increase of 1°C·d per day in the accumulated cold degree days, the relative content of the essential amino acids was increased by 0.01–0.41%, depending on the rice variety and growth stage. Our results suggest that LTS can improve nutritional quality of amino acids of rice grains in terms of amino acids content, especially at flowering stage. These results provide critical insights for assessing the potential impact of extreme climates on the nutrient quality of rice under future climate change.

Kang, Min↗