SEARCH · Search NASA

Results for “surrogate models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

930 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗

DOE OSTI2026

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

On-Orbit Calibration Assessment of NOAA-21 VIIRS Thermal Emissive Bands and Implications for JPSS-4 VIIRS

The Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard the NOAA-21 (N21) satellite has been successfully operating for over three years collecting valuable scientific measurements. A large suite of weather forecasting models and scientific research applications are supported with the VIIRS measurements from N21 combined with those acquired from the VIIRS instrument onboard the S-NPP and NOAA-20 spacecraft. Among the 22 VIIRS spectral bands, there are seven thermal emissive bands (TEB) covering the 3.7 to 12.2 µm spectral range at two different spatial resolutions. The VIIRS TEB detectors are calibrated onboard using a blackbody at a controlled temperature and a space view for background signal measurement. In this paper, we discuss the on-orbit performance of the N21 VIIRS TEB using various parameters such as the detector gain, noise, and offset. This on-orbit assessment provides critical insights into instrument behavior and calibration performance that will inform pre-launch testing and post-launch validation strategies for the upcoming JPSS-4 VIIRS mission scheduled for launch in 2027. The methodologies and lessons learned from characterizing the N21 TEB performance will enable more efficient and comprehensive radiometric assessment of the next VIIRS instrument. Specific focus is placed on identifying performance trends and anomaly signatures observed during the N21 commissioning phase that can enhance the JPSS4 instrument checkout procedures and accelerate the transition to routine on-orbit operations. The established performance baselines and uncertainty estimates will serve as an acceptance criterion for the JPSS-4 commissioning activities, ensuring mission readiness and data continuity for operational users.

Amit Angal↗

DOE OSTI2026

Topographic effects on reflected acoustic waves from the OSIRIS-REx reentry observed from stratospheric balloons

During long-distance sound propagation in planetary atmospheres, acoustic waves may reflect off the air/surface interface one or more times. For low sound frequencies and flat interfaces, the incident and reflected wave tend to be nearly identical. However, this may not be the case when the downgoing acoustic wave encounters topography. Here, we describe a set of direct and reflected acoustic signals recorded on free-flying balloons during the hypersonic entry of the OSIRIS-REx sample return capsule (SRC). In two of the three cases presented here, an impulsive reflected arrival similar in form to the direct sonic boom of the SRC is observed, followed by a diffuse coda. In contrast, one of the floating stations lacked an impulsive reflection entirely, with only the coda present. We use full waveform modeling to show how reflection in the presence of complex topography can explain coda seen in all three examples as well as the lack of impulsive arrival on the third. Our results indicate that the complex signals often observed in long range acoustic propagation could be due, in part, to interactions with topography during transmission.

Lees, Jonathan M. [University of North Carolina, C↗

DOE OSTI2026

Characterization of Flashback and Flame-Holding in a Jet-in-Crossflow Mixing Configuration With Methane–Hydrogen Fuel Blends

An approach that combines experimental and numerical analyses has been implemented to characterize the fundamentals of flashback events and flame-holding phenomena during high-hydrogen combustion in a jet-in-crossflow (JICF) configuration. Such flame dynamics are visualized experimentally using nanosecond (ns)-based hydroxyl planar laser-induced fluorescence (OH-PLIF) and chemiluminescence diagnostics techniques. The JICF burner has an optically accessible premixing tube allowing the optical diagnostics. The testing was conducted for varied premixer velocities (V) and equivalence ratio (ϕ) for 90–100% (H 2 , by mole) H 2 /CH 4 reactant mixtures at atmospheric temperature and pressure conditions. Two distinct flashback events were identified – conventional rich flashback and lean flashback, recorded while increasing ϕ and decreasing ϕ, respectively. The cause of lean flashback was attributed to lower momentum flux ratio which bends the jet sharply, closer toward the injection plane. The mean OH-PLIF images characterized the flame-holding behavior, where the flame was found to be stabilized on the leeward side only or on both windward and leeward sides as a lifted flame near the fuel port. A large eddy simulation (LES) with detailed chemistry combustion modeling approach was implemented along with an OH* submechanism, and it showed qualitative agreement with the integrated line-of-sight chemiluminescence results as well as the planar OH-PLIF measurements.

chemiluminescence diagnostic↗

DOE OSTI2026

Eddy Covariance Theory: A Review

Eddy covariance (EC), the gold standard for measuring ecosystem scale gas and heat exchanges, has transformed our understanding of the breathing of the biosphere, and thus global change biology. Despite numerous methodological improvements and insights gained from the technique, the community faces persistent challenges that have been present since the first EC measurements. Here, we review the theoretical developments underpinning EC. We present theoretical developments in four important areas that have relevance to EC measurements of the net ecosystem exchanges (NEE) of gases and heat from a single tower: (i) measuring the total vertical flux density, (ii) flux attenuation, (iii) coordinate rotations, and (iv) energy balance closure. Persistent problems with EC measurements, such as the inability to close the energy budget, led us to identify two priorities for revisiting the theory underlying: (i) sensible heat flux calculations, and (ii) constraining the mean vertical wind velocity. We present a framework for improved calculation of sensible heat flux derived from first principles of fluid mechanics and thermodynamics that considers coupled heat and mass transfer so that conservation of both is obeyed. These refinements are motivated by the need for unbiased measurements of energy and mass transfer between the land surface and atmosphere for ecosystem research and to validate satellite observations and land surface models.

ecosystem fluxes↗

NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

pragmatic method↗

DOE OSTI2024

Osmotic and phoretic competition explains chemotaxic assembly and sorting

Microscale objects responding to chemical gradients by migrating toward or away from a preferred species is a simple yet constitutive mechanism by which transport occurs in biological organisms. Synthetic chemotaxis provides key physical descriptions of simplified systems that can be used in biological models, or in the creation of advanced responsive material systems. In this article, we provide a quantitative framework for understanding synthetic chemotaxis of microparticles which involves a competition between phoresis and osmosis. We present separate quantitative measurements of phoresis and osmosis acting on individual taxing particles, finding that phoresis follows the long-predicted v ∼ 1 / r 2 scaling while the osmotic contribution depends on the geometry and details of the system, and must be solved on a case-by-case basis. Through this, we are able to develop a more accurate picture of particle transport at the single particle level. Equipped with this approach, we go on to describe how high concentrations of particles in a symmetric chemical gradient grow close-packed hives that reach a steady-state size tunable through light intensity or particle size. Last, we demonstrate that mixed particles experiencing the same chemical gradient will selectively migrate toward or away depending on the nature of the particle surface, thereby locally sorting out a particular species. We anticipate these results will be important in describing both biological and synthetic chemotaxis in phoretic systems and should bring a wealth of studies that take advantage of competing osmotic flows to illicit unexpected dynamic active behavior.

Science & Technology - Other Topics↗

NASA NTRS2026

Proactive Wildfire Management: A Remote Sensing and Multimodal CNN-MLP Architecture for Ignition Risk Forecasting

As the frequency and intensity of wildfires increase, with fire seasons now starting earlier and ending later than they have over the past decades, current monitoring systems, such as lookout towers and satellites, are hindered by cloud cover, low-resolution imagery, and static data gaps that fail to track vegetation moisture levels fast enough to catch rapid pre-ignition changes. This report proposes a Machine Learning-enabled Wildfire Ignition Prediction framework that combines satellite monitoring with dynamic and high-resolution remote sensing from Unmanned Aerial Vehicle (UAV) swarms. The method would use multispectral and thermal data from the Landsat program to create a baseline for vegetation health, calculating a two-band Enhanced Vegetation Index (EVI2) and the moisture content of the vegetation. These inputs will later be fused with microscale UAV weather data, including thermal hotspots found through thick canopies, hyperspectral chemical signatures of pre-visual combustion, and local weather streams. The multispectral satellite, multispectral Light Detection and Ranging (LiDAR), and thermal data would then be processed through a Convolutional Neural Network (CNN), alongside a Multilayer Perceptron (MLP) for the micro-weather telemetry. The outputs of these networks would be fused into a single feature representation and passed through a final prediction network to generate real-time ignition risk scores and hotspot alerts. Model performance would be assessed using standard classification metrics, including a Receiver Operating Characteristic - Area Under the Curve (ROC AUC) and F1 score. This system would allow first responders to identify high-risk zones and intervene before ignition occurs, improving emergency response time compared to current approaches.

machine learning↗

NASA NTRS2014

High Efficiency Centrifugal Compressor for Rotorcraft Applications

The report "High Efficiency Centrifugal Compressor for Rotorcraft Applications" documents the work conducted at UTRC under the NRA Contract NNC08CB03C, with cost share 2/3 NASA, and 1/3 UTRC, that has been extended to 4.5 years. The purpose of this effort was to identify key technical barriers to advancing the state-of-the-art of small centrifugal compressor stages; to delineate the measurements required to provide insight into the flow physics of the technical barriers; to design, fabricate, install, and test a state-of-the-art research compressor that is representative of the rear stage of an axial-centrifugal aero-engine; and to acquire detailed aerodynamic performance and research quality data to clarify flow physics and to establish detailed data sets for future application. The design activity centered on meeting the goal set outlined in the NASA solicitation-the design target was to increase efficiency at higher work factor, while also reducing the maximum diameter of the stage. To fit within the existing Small Engine Components Test Facility at NASA Glenn Research Center (GRC) and to facilitate component re-use, certain key design parameters were fixed by UTRC, including impeller tip diameter, impeller rotational speed, and impeller inlet hub and shroud radii. This report describes the design effort of the High Efficiency Centrifugal Compressor stage (HECC) and delineation of measurements, fabrication of the compressor, and the initial tests that were performed. A new High-Efficiency Centrifugal Compressor stage with a very challenging reduction in radius ratio was successfully designed, fabricated and installed at GRC. The testing was successful, with no mechanical problems and the running clearances were achieved without impeller rubs. Overall, measured pressure ratio of 4.68, work factor of 0.81, and at design exit corrected flow rate of 3 lbm/s met the target requirements. Polytropic efficiency of 85.5 percent and stall margin of 7.5 percent were measured at design flow rate and speed. The measured efficiency and stall margin were lower than pre-test CFD predictions by 2.4 percentage points (pt) and 4.5 pt, respectively. Initial impressions from the experimental data indicated that the loss in the efficiency and stall margin can be attributed to a design shortfall in the impeller. However, detailed investigation of experimental data and post-test CFD simulations of higher fidelity than pre-test CFD, and in particular the unsteady CFD simulations and the assessment with a wider range of turbulence models, have indicated that the loss in efficiency is most likely due to the impact of unfavorable unsteady impeller/diffuser interactions induced by diffuser vanes, an impeller/diffuser corrected flow-rate mismatch (and associated incidence levels), and, potentially, flow separation in the radial-to-axial bend. An experimental program with a vaneless diffuser is recommended to evaluate this observation. A subsequent redesign of the diffuser (and the radial-to-axial bend) is also recommended. The diffuser needs to be redesigned to eliminate the mismatching of the impeller and the diffuser, targeting a slightly higher flow capacity. Furthermore, diffuser vanes need to be adjusted to align the incidence angles, to optimize the splitter vane location (both radially and circumferentially), and to minimize the unsteady interactions with the impeller. The radial-to-axial bend needs to be redesigned to eliminate, or at least minimize, the flow separation at the inner wall, and its impact on the flow in the diffuser upstream. Lessons were also learned in terms of CFD methodology and the importance of unsteady CFD simulations for centrifugal compressors was highlighted. Inconsistencies in the implementation of a widely used two-equation turbulence model were identified and corrections are recommended. It was also observed that unsteady simulations for centrifugal compressors require significantly longer integration times than what is current practice in industry.

Gorazd Medic↗

NASA NTRS1998

Behavior of 2-Cell CPV Ni-H2 Battery During Pulse Discharge

A study was carried out to determine the transient voltage behavior of the 2-cell CPV nickel-hydrogen battery with the objective of using the results as a basis for mathematical modeling. The 2-cell CPV battery which is manufactured by Eagle Picher, Inc. for the GOES program yields 18.5 Ah at C/2 rate of discharge at 10 C with a mid-discharge voltage of 2.514 V. The capacity increased with decrease of temperature and a maximum capacity of 22 Ah was obtained at -5 C. The pulse tests consisted of obtaining the voltage profile in the first 20 milliseconds of the one minute pulse discharge at 37 A and pulse discharge was repeated as a function of state-of-charge. The pulse test at 10 C and 20 C provided voltage profiles with the expected decrease in voltage as the pulse was applied. The end of pulse voltage decreased with the state-of-charge. The battery voltage was above 2V at the end of the one-minute pulse at 8 % state-of-charge at 10 C. The voltage profile during the 37 A pulse discharge consists of an initial drop in voltage which was independent of the state of charge. The invariability in the value for the initial drop in voltage with state of charge is a very important observation. The results show that towards the end of discharge the dominant resistance is not ohmic in nature. It could be hypothesized from the nature of the voltage transients that the dominant mechanism towards the end of discharge is proton diffusion. The study also shows that the dominant resistance in the voltage plateau during discharge is activation polarization.

Vaidyanathan, Hari↗

NASA NTRS2024

Artemis I Orion ESM Propulsion System Engine Performance

NASA's Orion spacecraft transports humans and cargo into cislunar space for the Artemis program. The European Service Module (ESM), supplied by ESA and its European industry partners, provides Orion with power and in-space propulsion. The Orion-ESM propulsion system is a bipropellant hypergolic propulsion system using monomethyl hydrazine (MMH) and nitrogen tetroxide (MON-3). Primary translational propulsion is provided by the Orbital Maneuvering System Engine(OMS-E), with backup translational propulsion provided by eight Auxiliary thrusters (AUX). Attitude control and small translational maneuvers are provided by twenty four Reaction Control System (RCS) engines. The 2022 Artemis I mission was the first integrated flight test of the Orion-ESM spacecraft and its propulsion system. The OMS-E used on Artemis I was a refurbished Space Shuttle OMS-E that previously flew on nineteen missions ranging from STS-41G in 1984 to STS-112 in 2002. The Auxiliary engines are modified Aerojet Rocketdyne R4D-11 engines produced specifically for the Orion program. The RCS engines are Ariane Group engines originally used for the Automated Transfer Vehicle (ATV) program. This paper will discuss the unique operational requirements for each engine on Orion and the development and qualification effort sat both the engine and system-level that were completed to enable a successful Artemis I mission. Next the paper will evaluate the in-flight performance of the engines during the Artemis I mission showing nominal performance as expected. Additionally, comparisons to models will be presented showing very good correlation. Finally, the paper will address the plan for the engines on future Orion missions and the evolution of the system operation.

Liquid Propulsion Systems↗

NASA NTRS2024

Artemis I Orion ESM Propulsion System Engine Performance

Reuse↗

NASA NTRS2012

General Purpose Data-Driven System Monitoring for Space Operations

Modern space propulsion and exploration system designs are becoming increasingly sophisticated and complex. Determining the health state of these systems using traditional methods is becoming more difficult as the number of sensors and component interactions grows. Data-driven monitoring techniques have been developed to address these issues by analyzing system operations data to automatically characterize normal system behavior. The Inductive Monitoring System (IMS) is a data-driven system health monitoring software tool that has been successfully applied to several aerospace applications. IMS uses a data mining technique called clustering to analyze archived system data and characterize normal interactions between parameters. This characterization, or model, of nominal operation is stored in a knowledge base that can be used for real-time system monitoring or for analysis of archived events. Ongoing and developing IMS space operations applications include International Space Station flight control, spacecraft vehicle system health management, launch vehicle ground operations, and fleet supportability. As a common thread of discussion this paper will employ the evolution of the IMS data-driven technique as related to several Integrated Systems Health Management (ISHM) elements. Thematically, the projects listed will be used as case studies. The maturation of IMS via projects where it has been deployed or is currently being integrated to aid in fault detection will be described. The paper will also explain how IMS can be used to complement a suite of other ISHM tools, providing initial fault detection support for diagnosis and recovery

Space Propulsion↗

NASA NTRSDate not supplied

Mass Economy Evaluation for Integrated ECLSS and Propulsion Architecture

As missions in Low Earth Orbit (LEO) lengthen and extend to deep space, minimizing resupply needs becomes vital for sustaining crewed operations. Traditional life support systems depend on consumables resupplied from Earth, a method that is increasingly impractical for missions beyond LEO, such as lunar outposts or Mars transit. Long-duration missions require more efficient, autonomous systems that can recycle essential resources, particularly water and oxygen, to minimize the frequency and mass of resupply missions. The Environmental Control and Life Support System (ECLSS) is essential to such missions, with the International Space Station (ISS) serving as a testbed for advanced water recovery and partial oxygen recycling via physico-chemical methods. Yet, ECLSS and propulsion subsystems generally operate independently, despite overlapping requirements and potential areas for synergy. For instance, ECLSS byproducts, water, CO₂, and hydrogen, could be repurposed for propulsion, potentially reducing dedicated propellant mass and increasing overall system efficiency. One promising approach is to develop shared-resource architectures that integrate ECLSS with propulsion systems. This study examines the potential of such integration through the Sabatier CO₂ reduction process, focusing on water management as a key factor in system mass trade-offs. The Sabatier reaction produces water and methane from metabolic CO₂ and electrolytic hydrogen, partially closing the life support loop and providing methane, which could serve as a propellant. This integration could minimize waste, reduce resupply requirements, and enhance mission mass efficiency. A dynamic modeling framework will be used to simulate resource flows over long missions, capturing interactions between life support and propulsion. By comparing integrated versus separate system configurations, the study aims to quantify mass benefits and penalties, informing future habitat designs and trade studies for missions prioritizing autonomy and mass efficiency.

ECLSS↗

NASA NTRS2026

Micromechanical Design of Carbon Nanotube Ribbon Reinforced Polymer Composite Materials

Lightweight materials are an important component of the design of aerospace structures. Carbon nanotube materials have been considered for this purpose due to the strength and stiffness of individual nanotubes, and the commercial availability of bulk formats such as fibers. These fibers can have a ribbon cross section which results in a different design space for their composites relative to traditional reinforcements which have a round cross section. This work applies brick-and-mortar micromechanical models and classical lamination theory with an inverse approach to investigate the design space of these composites. Using this approach, the influences of fiber geometry and axial and transverse mechanical properties are mapped. Finally, a sensitivity study is performed and the relative impacts of ±10% variations in the constituent material and geometric properties are ranked. Lamina axial moduli were found to range from a maximum of 3x to 1x minimum relative to a target quasi-isotropic laminate modulus depending on the anisotropy and shear modulus of the lamina. The fiber targets depended strongly on the fiber volume fraction in the lamina and the fiber axial modulus target was found to range from 4.8x to 3.2x the quasi-isotropic laminate target. The sensitivity analysis found that the largest driver of performance was the volume fraction, followed by the fiber axial modulus. While bio-based brick-and-mortar composites, such as nacre, can benefit from reinforcement aspect ratios above 10, for carbon nanotube ribbon (or carbon fiber)/polymer composites the sensitivity study indicated that the optimal cross-sectional aspect ratio was relatively smaller, potentially less than three.

Micromechanics↗

NASA NTRS2026

Towards an Improved Understanding of the Antarctic Coastal Zone and Its Contribution to Future Global Sea Level

Understanding the coastal zone of the Antarctic Ice Sheet (AIS), where it interacts with the Southern Ocean and warmer air masses, is crucial for predicting Antarctica's influence on the global climate and sea level. This region has multiple tipping mechanisms that could trigger large, rapid, and potentially irreversible changes in the AIS, the Southern Ocean and their global connections in the coming centuries. The AIS remains the largest source of uncertainty in future sea-level projections. Bed topography beneath the ice shelves and the coastal ice sheet is not yet well documented, and is a major source of this uncertainty. This review assesses current knowledge of the coastal zone and highlights methods to investigate it, including aerogeophysical surveys, ground- and ship-based measurements, satellite observations, and computer modeling. An ensemble analysis of published bed topography data sets identifies significant data gaps and their regional distribution, framed in the context of current ice-sheet behavior and potential instability. We propose scientific priorities and guidelines for future aerogeophysical surveys, advocating for a comprehensive, coordinated international effort to build a next-generation data set of Antarctic bed properties. Such an initiative would significantly advance understanding of the role of coastal processes in ice-sheet dynamics, reducing uncertainties in sea-level rise projections and improving predictions of future ocean and climate changes.

Kenichi Matsuoka↗

DOE OSTI2025

Photon–photon chemical thermodynamics of frequency conversion processes in highly multimode systems

Abstract Frequency generation in highly multimode nonlinear optical systems is inherently a complex process, giving rise to an exceedingly convoluted landscape of evolution dynamics. While predicting and controlling the global conversion efficiencies in such nonlinear environments has long been considered impossible, here, we formally address this challenge even in scenarios involving a very large number of spatial modes. By utilizing fundamental notions from optical statistical mechanics, we develop a universal theoretical framework that effectively treats all frequency components as chemical reactants/products, capable of undergoing optical thermodynamic reactions facilitated by a variety of multi-wave mixing effects. These photon–photon reactions are governed by conservation laws that directly determine the optical temperatures and chemical potentials of the ensued chemical equilibria for each frequency species. In this context, we develop a comprehensive stoichiometric model and formally derive an expression that relates the chemical potentials to the optical stoichiometric coefficients, in a manner akin to atomic/molecular chemical reactions. This advancement unlocks new predictive capabilities that can facilitate the optimization of frequency generation in highly multimode photonic arrangements, surpassing the limitations of conventional schemes that rely exclusively on nonlinear optical dynamics. Notably, we identify a universal regime of Rayleigh–Jeans thermalization where an optical reaction at near-zero optical temperatures can promote the complete and entropically irreversible conversion of light to the fundamental mode at a target frequency. Our theoretical results are corroborated by numerical simulations in settings where second-harmonic generation, sum-frequency generation and four-wave mixing processes can manifest.

Optics↗