SEARCH · Search NASA
Results for “activity coefficient model”
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.
Substantial oxygen loss and chemical expansion in lithium-rich layered oxides at moderate delithiation
Delithiation of layered oxide electrodes triggers irreversible oxygen loss, one of the primary degradation modes in lithium-ion batteries. However, the delithiation-dependent mechanisms of oxygen loss remain poorly understood. Here we investigate the oxygen non-stoichiometry in Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ electrodes as a function of Li content by using cycling protocols with long open-circuit voltage steps at varying states of charge. Surprisingly, we observe substantial oxygen loss even at moderate delithiation, corresponding to 2.5, 4.0 and 7.6 ml O 2 per gram of Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ after resting at upper capacity cut-offs of 135, 200 and 265 mAh g −1 for 100 h. Our observations suggest an intrinsic oxygen instability consistent with predictions of high oxygen activity at intermediate potentials versus Li/Li + . In addition, we observe a large chemical expansion coefficient with respect to oxygen non-stoichiometry, which is about three times greater than those of classical oxygen-deficient materials such as fluorite and perovskite oxides. Furthermore, our work challenges the conventional wisdom that deep delithiation is a necessary condition for oxygen loss in layered oxide electrodes and highlights the importance of calendar ageing for investigating oxygen stability.
Characterizing the Variation and Covariation of Cloud Microphysical Properties and Implications for Simulation of Subgrid-scale Warm-Rain Processes in Earth System Models (Final DOE-ASR Report)
Warm marine boundary layer (MBL) clouds constitute an important component in the global climate system, and precipitation plays a central role in controlling the water budget, radiative effects, and lifetime of these MBL clouds. Unfortunately, because of the relatively coarse effective grid resolution of the current generation of Earth system models (ESMs), the variety of cloud microphysical processes occurring inside an ESM grid cell are often oversimplified or unconstrained by observations. For example, the warm rain processes (e.g., autoconversion and accretion) are usually parameterized as nonlinear functions of grid-mean cloud properties. Because of the nonlinear nature of these functions, neglecting variability within the ESM grid volume can lead to substantial biases in precipitation production, cloud cover, and surface radiative fluxes. In state-of-the art ESMs, the influence of subgrid-scale variability is represented as an enhancement factor (EF) coefficient to the autoconversion, and accretion rates calculated from the model variables. However, EF is typically taken to be a constant or even used as a knob to tune model cloud properties to match observations, an ad hoc approach that may yield a desired cloud outcome yet introduce compensating errors. In this project, we used the combination of in situ cloud microphysics measurements from the ACE-ENA field campaign and large-eddy simulations (LES) to characterize and understand subgrid-scale variations and co-variations of cloud microphysical properties and use the results to evaluate and improve the representation of subgrid warm-rain processes in ESMs, in particular the EF used to tune the autoconversion and accretion processes. In this final report, we summarize our research activities and main findings in Section 2, provide a list of publications (Section 3) and presentations (Section 4) resulted from our research, and briefly discuss the student activities supported by this project.
Paper #50 Preliminary Sizing of Active Magnetic Bearings for sCO2 Waste-Heat Recovery Application
The team of Southwest Research Institute (SwRI) and GE Vernova (GE) is executing a project developing conceptual designs of high-temperature active magnetic bearings (AMBs) for sCO2 machinery applications. The goals of the project are to develop conceptual designs for radial and thrust AMBs capable of operating in sCO2 environments up to 1000°F. A waste-heat recovery (WHR) application with hermetic machinery is used as a reference case to apply AMBs. This application was previously investigated to compare hermetic machinery enabled by process-lubricated bearings and high-speed motors/generators with traditional machinery configurations using oil-lubricated bearings, gearboxes, and grid-frequency electric machines. In this paper, AMBs are sized for the WHR application machinery. Load capacity estimates are made using first-principle-derived formulas that consider operation at elevated temperatures. Estimates of linearized AMB system coefficients are also made to aid future controller development. Rotor models are developed that account for AMB rotor geometry and housing space requirements. Rotordynamics studies are performed to identify reasonable target values for closed-loop AMB stiffness and closed-loop AMB damping. These target values will be used in future control studies.
Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers
Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.
Electro-physical properties of surface-barrier diodes on low-resistance n-Cd 0.96 Mn 0.04 Te 0.96 Se 0.04 crystals
The electro-physical properties of Cd 0.96 Mn 0.04 Te 0.96 Se 0.04 n-type conductivity single crystals with relatively low resistivity were studied. The determined resistivity of semiconductor single crystals Cd 0.96 Mn 0.04 Te 0.96 Se 0.04 n-type conductivity is equal to ρ≈2 Ω-cm at 293K. The activation energy of dark conductivity, measured from the temperature dependence of the resistivity of the In/Cd 0.96 Mn 0.04 Te 0.96 Se 0.04 /In structure, is equal to ΔE≈0.11eV. From optical measurements of the absorption coefficient from the energy of photons in the region of large α, the band-gap width of the single crystal was found, which is equal to E g =1.53eV. Diode structures of Ni/Cd 0.96 Mn 0.04 Te 0.96 Se 0.04 /In obtained by vacuum sputtering of Ni on the surface of a single crystal n-Cd 0.96 Mn 0.04 Te 0.96 Se 0.04 were manufactured and studied. Here, the rectification coefficient at a voltage of U=1V is equal to 10 8 . Within the framework of the Sah-Noys-Shockley carrier generation-recombination model, a quantitative match of the experimental measurements of the volt-ampere characteristics with the calculation was achieved.
First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport
Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.
Evaluation of a quasi-steady-state approximation of the cloud droplet growth equation (QDGE) scheme for aerosol activation in global models using multiple aircraft data over both continental and marine environments
This research introduces a numerically efficient aerosol activation scheme and evaluates it by using stratus and stratocumulus cloud data sampled during multiple aircraft campaigns in Canada, Chile, Brazil, and China. The scheme employs a quasi-steady-state approximation of the cloud droplet growth equation (QDGE) to efficiently simulate aerosol activation, the vertical profile of supersaturation, and the activated cloud droplet number concentration (CDNC) near the cloud base. The calculated maximum supersaturation values using the QDGE scheme were compared with multiple parcel model simulations under various aerosol and environmental conditions. The differences are all below 0.18 %, indicating good performance and accuracy of the QDGE scheme. We evaluated the QDGE scheme by specifying observed environmental thermodynamic variables and aerosol information from 31 cloud cases as input and comparing the simulated CDNC with cloud observations. The average of mean relative error ($\overline{MRE}$) of the simulated CDNC for cloud cases in each campaign ranges from 17.30 % in Brazil to 25.90 % in China, indicating that the QDGE scheme successfully reproduces observed variations in CDNC over a wide range of different meteorological conditions and aerosol regimes. Additionally, we carried out an error analysis by calculating the maximum information coefficient (MIC) between the MRE and input variables for the individual campaigns and all cloud cases. MIC values were then sorted by aerosol properties, pollution level, environmental humidity, and dynamic condition according to their relative importance to MRE. Based on the error analysis, we found that the magnitude of MRE is more relevant to the specification of input aerosol pollution level in marine regions and aerosol hygroscopicity in continental regions than to other variables in the simulation.
New Developments and Capabilities Within WEC-Sim
WEC-Sim is an open-source software for simulating wave energy converters and has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory and Sandia National Laboratories within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with the latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1, which spans a 2-year period from June 2020 to September 2022. New features described here include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.
Metabolic Costs of Walking with Weighted Vests
ABSTRACT Introduction The US Army Load Carriage Decision Aid (LCDA) metabolic model is used by militaries across the globe and is intended to predict physiological responses, specifically metabolic costs, in a wide range of dismounted warfighter operations. However, the LCDA has yet to be adapted for vest-borne load carriage, which is commonplace in tactical populations, and differs in energetic costs to backpacking and other forms of load carriage. Purpose The purpose of this study is to develop and validate a metabolic model term that accurately estimates the effect of weighted vest loads on standing and walking metabolic rate for military mission-planning and general applications. Methods Twenty healthy, physically active military-age adults (4 women, 16 men; age, 26 ± 8 yr old; height, 1.74 ± 0.09 m; body mass, 81 ± 16 kg) walked for 6 to 21 min with four levels of weighted vest loading (0 to 66% body mass) at up to 11 treadmill speeds (0.45 to 1.97 m·s −1 ). Using indirect calorimetry measurements, we derived a new model term for estimating metabolic rate when carrying vest-borne loads. Model estimates were evaluated internally byk-fold cross-validation and externally against 12 reference datasets (264 total participants). We tested if the 90% confidence interval of the mean paired difference was within equivalence limits equal to 10% of the measured walking metabolic rate. Estimation accuracy, precision, and level of agreement were also evaluated by the bias, standard deviation of paired differences, and concordance correlation coefficient (CCC), respectively. Results Metabolic rate estimates using the new weighted vest term were statistically equivalent (P< 0.01) to measured values in the current study (bias, −0.01 ± 0.54 W·kg −1 ; CCC, 0.973) as well as from the 12 reference datasets (bias, −0.16 ± 0.59 W·kg −1 ; CCC, 0.963). Conclusions The updated LCDA metabolic model calculates accurate predictions of metabolic rate when carrying heavy backpack and vest-borne loads. Tactical populations and recreational athletes that train with weighted vests can confidently use the simplified LCDA metabolic calculator provided as Supplemental Digital Content to estimate metabolic rates for work/rest guidance, training periodization, and nutritional interventions.
Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media
We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.
Determining Activity Coefficients of SOA from Isothermal Evaporation in a Laboratory Chamber
Nonideal molecular interactions in aerosol particles influence the partitioning of semivolatile organic compounds (SVOCs). However, few direct measurements exist that determine activity coefficients (γ), which quantify nonideality, for individual organic compounds in different atmospheric mixtures. By measuring the isothermal evaporation of SOA formed from multifunctional SVOCs in an environmental chamber, this study determines the mole-fraction-based γ of the bulk SOA in a variety of pre-existing chemically complex particles, approximated as if the bulk SOA behaved like a single SVOC. These multicomponent mixtures contain molecular structures commonly found in the atmosphere. When treated as a single liquid phase, γ of the bulk SOA tends to increase from ~1 to ~5 as the organic seeds and the SOA have more differing polarities. A high computed γ value of 74 for a wet ammonium sulfate-SOA system indicates phase separation. The γ for some individual species was also quantified based on gas-phase measurements and AIOMFAC model predictions. However, the bulk SOA γ cannot be explained by the simplified speciated SVOC–seed interactions, suggesting incomplete compositional understanding of the bulk SOA. In this work, these results demonstrate a method to quantify nonideal behavior and show that it can occur in multicomponent mixtures and therefore influence OA formation, evolution, and lifetimes.
Global Simulations of Phase State and Equilibration Time Scales of Secondary Organic Aerosols with GEOS-Chem
The phase state of secondary organic aerosols (SOA) can range from liquid through amorphous semisolid to glassy solid, which is important to consider as it influences various multiphase processes including SOA formation and partitioning, multiphase chemistry, and cloud activation. In this study, we simulate the glass transition temperature and viscosity of SOA over the globe using the global chemical transport model, GEOS-Chem. The simulated spatial distributions show that SOA at the surface exist as liquid over equatorial regions and oceans, semisolid in the midlatitude continental regions, and glassy solid over lands with low relative humidity. The predicted SOA viscosities are mostly consistent with the available measurements. In the free troposphere, SOA particles are mostly predicted to be semisolid at 850 hPa and glassy solid at 500 hPa, except over tropical regions including Amazonia, where SOA are predicted to be low viscous. Phase state also exhibits seasonal variation with a higher frequency of semisolid and solid particles in winter compared to warmer seasons. We calculate equilibration time scales of SOA partitioning (τ eq ) and effective mass accommodation coefficient (α eff ), indicating that τ eq is shorter than the chemical time step of GEOS-Chem of 20 min and αeff is close to unity for most locations at the surface level, supporting the application of equilibrium SOA partitioning. However, τ eq is prolonged and α eff is lowered over drylands and most regions in the upper troposphere, suggesting that kinetically limited growth would need to be considered for these regions in future large-scale model studies.
New Developments and Capabilities Within WEC-Sim: Preprint
WEC-Sim is an open-source software for simulating wave energy converters, which has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (Sandia) within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with our latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1 which spans over a two year period from June 2020 to September 2022. New features to be described will include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.
ANN-based ground motion model for Turkey using stochastic simulation of earthquakes
SUMMARY Turkey is characterized by a high level of seismic activity attributed to its complex tectonic structure. The country has a dense network to record earthquake ground motions; however, to study previous earthquakes and to account for potential future ones, ground motion simulations are required. Ground motion simulation techniques offer an alternative means of generating region-specific time-series data for locations with limited seismic networks or regions with seismic data gaps, facilitating the study of potential catastrophic earthquakes. In this research, a local ground motion model (GMM) for Turkey is developed using region-specific simulated records, thus constructing a homogeneous data set. The simulations employ the stochastic finite-fault approach and utilize validated input-model parameters in distinct regions, namely Afyon, Erzincan, Duzce, Istanbul and Van. To overcome the limitations of linear regression-based models, artificial neural network is used to establish the form of equations and coefficients. The predictive input parameters encompass fault mechanism (FM), focal depth (FD), moment magnitude (Mw), Joyner and Boore distance (RJB) and average shear wave velocity in the top 30 m (Vs30). The data set comprises 7359 records with Mw ranging between 5.0 and 7.5 and RJB ranging from 0 to 272 km. The results are presented in terms of spectral ordinates within the period range of 0.03–2.0 s, as well as peak ground acceleration and peak ground velocity. The quantification of the GMM uncertainty is achieved through the analysis of residuals, enabling insights into inter- and intra-event uncertainties. The simulation results and the effectiveness of the model are verified by comparing the predicted values of ground motion parameters with the observed values recorded during previous events in the region. The results demonstrate the efficacy of the proposed model in simulating physical phenomena.
A systematic analysis of Cu-based membrane-electrode assemblies for CO 2 reduction through multiphysics simulation
A multiphysics model is presented to study potential losses in Cu-MEAs and how various physical phenomena impact the product distribution of CO 2 reduction.
3D reconstruction of an inertial-confinement fusion implosion with neural networks using multiple heterogeneous data sources
3D asymmetries are major degradation mechanisms in inertial-confinement fusion implosions at the National Ignition Facility (NIF). These asymmetries can be diagnosed and reconstructed with the neutron imaging system (NIS) on three lines of sight around the NIF target chamber. Conventional tomographic reconstructions are used to reconstruct the 3D morphology of the implosion using NIS [Volegov et al., J. Appl. Phys. 127, 083301 (2020)], but the problem is ill-posed with only three imaging lines of sight. Asymmetries can also be diagnosed with the real-time neutron activation diagnostics (RTNAD) and the neutron time-of-flight (nToF) suite. Since the NIS, RTNAD, and nToF each sample a different part of the implosion using different physical principles, we propose that it is possible to overcome the limitations of too few imaging lines of sight by performing 3D reconstructions that combine information from all three heterogeneous data sources. This work presents a new machine learning-based reconstruction technique to do just this. By using a simple physics model and group of neural networks to map 3D morphologies to data, this technique can easily account for data of multiple different types. A simple proof-of-principle is presented, demonstrating that this technique can accurately reconstruct a hot-spot shape using synthetic primary neutron images and a hot-spot velocity vector. In particular, the hot-spot’s asymmetry, quantified as spherical harmonic coefficients, is reconstructed to within ±4% of the radius in 90% of test cases. In the future, this technique will be applied to actual NIS, RTNAD, and nToF data to better understand 3D asymmetries at the NIF.
Chemistry of Sugar Formation in the Gas Phase: Following the Activated Aldehyde
Sugars are produced by living organisms, and are required building blocks for life as we know it, which raises the foundational question of how sugars formed in a prebiotic environment. The abiotic formose reaction produces sugars from formaldehyde, but our understanding of its initiation step remains murky, with chemists invoking the concept of an “activated aldehyde” to seed this reaction. Singlet hydroxycarbenes, high-energy isomers of aldehydes, were recently reported to facilitate sugar formation under cold, nonaqueous conditions relevant to interstellar environments. Here, we generate singlet methylhydroxycarbene ( 1 CH 3 –C̈–OH) from the photodissociation of pyruvic acid and experimentally measure its gas-phase reaction with d 4 -acetaldehyde using multiplexed photoionization mass spectrometry. The C 4 H 4 D 4 O 2 isomer d 4 -acetoin is the sole product, which we kinetically link to the reactant CH 3 –C̈–OH, and attribute to a carbonyl-ene formation mechanism. We see no evidence of 3-hydroxybutanal, the C–H insertion product expected in carbene chemistry. Using automated exploration we calculate stationary points on the potential energy surface and report master equation rate coefficients from T = 20–600 K, providing quantitative kinetics of this fast reaction for use in chemical models. The prereactive complex in this reaction is stabilized by both hydrogen bonding and electrophilic carbene C═O interactions. These effects create a short-range dynamical bottleneck for the reaction besides the long- and midrange barrierless bottlenecks. Combined with recent reports of 1 HC̈OH production from methanol photodissociation and pyruvic acid production in cold irradiated ices, this work provides evidence that singlet hydroxycarbene + aldehyde chemistry is a feasible path to prebiotic sugar formation.