SEARCH · Search NASA

Results for “Kinetic Modelling”

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

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS2024

Fundamental Experimental Tests and Modeling of LOX/CH4 Engines at High Pressures

In recent years, the space travel industry has grown exponentially, resulting in the need for a low-cost, efficient rocket engine fuel. As such, there has been renewed interest in utilizing liquid natural gas (LNG) as it is less likely to soot than kerosene-based fuels and has widespread availability and low cost of use as compared to the traditional RP-X and even liquid methane (LCH4) fuels due to the reduced need for refinement. However, current literature does not fully cover varying blends of natural gas and the effects that impurities have on natural gas/methane fuels, with most studies being confined to low-pressure applications. Therefore, new experimental ignition delay time measurements at rocket engine relevant pressures are essential to ensure current chemical kinetic models capture the ignition behavior at these elevated conditions and to determine the optimal purity levels for a reliable, cost-efficient, aerospace-grade LCH4/LNG fuel. The current work explores various blends of natural gas/methane fuel with various impurities, including higher hydrocarbons as well as nitrogen-carriers. A shock tube study was carried out to study the ignition delay times and carbon monoxide time histories of these natural gas blends utilizing chemiluminescence and laser absorption spectroscopy at 20, 50 and 100 atm over temperatures ranging from 1400-1700K and at an equivalence ratio of 1. The experimental data was then compared to two different chemical kinetic mechanisms: the industry-standard GRI 3.0 and the in-house UCF 2022. The data is used to improve chemical kinetic mechanisms and modeling of LNG fuels for rocket engine applications.

Methane↗

DOE OSTI2025

Active Site Dynamics in Molybdenum-Based Silica-Supported Olefin Metathesis Catalysts: Site Renewal and Decay Beyond the Chauvin Cycle

Heterogeneous olefin metathesis catalysts exhibit low active site densities and unpredictable kinetics due to dynamic active site formation and decay processes. Here, in this study, we establish a quantitative framework that captures active site generation, renewal, and decay in olefin metathesis over silica-supported molybdenum oxide catalysts, enabling a mechanistic explanation of catalytic behavior and strategies to achieve high, stable activity. Steady-state active site titrations reveal that 2,3-dimethyl-butene isomers (4MEs) cofeeding increases active site density by up to 4.3-fold, directly correlating with enhanced metathesis rates. Spectroscopic studies demonstrate that 4MEs facilitate Mo(VI) reduction to Mo(IV) and interact strongly with surface Si–OH groups, generating labile protons that promote active site formation via a 1,2-proton shift mechanism. Kinetic modeling indicates that ethylene acts as a decay promoter, shifting kinetic control away from the Chauvin cycle and suppressing metathesis activity. Comparative studies on catalysts with varying Mo loading reveal that promotion is most effective for dispersed molybdate species, with a decline at higher Mo loadings. These findings provide a unified mechanistic framework for heterogeneous olefin metathesis, offering new strategies to enhance active site accessibility, mitigate deactivation, and optimize catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

New physics in multi-lepton tau decays

Dark particles with lepton-flavor-violating couplings to the tau lepton can induce rare neutrinoless $τ$ decays with large final state multiplicities. We study models where transitions of the type $τ^\pm\to \ell^\pm\,ϕ$, with $ϕ$ a light new particle, initiate a chain of decays in the dark sector that terminate with decays into electrons, muons, or pions. These decay cascades appear as rare five or even seven-body $τ$ decays with multiple reconstructable resonances. We survey several representative models: kinetically mixed dark photon, gauged $L_i-L_j$ models, and other more exotic charge assignments such as chiral $U(1)'$ extensions of the Standard Model. The main new ingredient is the possibility of flavor violation at very high scales. In these models, a number of channels that have not yet been searched for experimentally, such as $τ\to 5μ$, $τ\to 3μ\,2e$, $τ\to μ\,4e$, and hadronic channels like $τ\to μ\,4π$, typically dominate over the previously-considered signatures such as $τ\to 3μ$. While some of the models, such as the gauged $L_i-L_j$ ones, also contain more challenging channels with missing energy due to decays to neutrinos, they can still be searched for via fully visible channels.

Ema, Yohei [Florida U.]↗

DOE OSTI2023

Re-evaluating cloud chamber constraints on depositional ice growth in cirrus clouds – Part 1: Model description and sensitivity tests

Abstract. Ice growth from vapor deposition is an important process for the evolution of cirrus clouds, but the physics of depositional ice growth at the low temperatures (<235 K) characteristic of the upper troposphere and lower stratosphere is not well understood. Surface attachment kinetics, generally parameterized as a deposition coefficient αD, control ice crystal habit and also may limit growth rates in certain cases, but significant discrepancies between experimental measurements have not been satisfactorily explained. Experiments on single ice crystals have previously indicated the deposition coefficient is a function of temperature and supersaturation, consistent with growth mechanisms controlled by the crystal's surface characteristics. Here we use observations from cloud chamber experiments in the Aerosol Interactions and Dynamics in the Atmosphere (AIDA) aerosol and cloud chamber to evaluate surface kinetic models in realistic cirrus conditions. These experiments have rapidly changing temperature, pressure, and ice supersaturation such that depositional ice growth may evolve from diffusion limited to surface kinetics limited over the course of a single experiment. In Part 1, we describe the adaptation of a Lagrangian parcel model with the Diffusion Surface Kinetics Ice Crystal Evolution (DiSKICE) model (Zhang and Harrington, 2014) to the AIDA chamber experiments. We compare the observed ice water content and saturation ratios to that derived under varying assumptions for ice surface growth mechanisms for experiments simulating ice clouds between 180 and 235 K and pressures between 150 and 300 hPa. We found that both heterogeneous and homogeneous nucleation experiments at higher temperatures (>205 K) could generally be modeled consistently with either a constant deposition coefficient or the DiSKICE model assuming growth on isometric crystals via abundant surface dislocations. Lower-temperature experiments showed more significant deviations from any depositional growth model, with different ice growth rates for heterogeneous and homogeneous nucleation experiments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Efficient kinetic thermal inverse modeling for organic material decomposition

The prevalent use of organic materials in manufacturing is a fire safety concern, and motivates the need for predictive thermal decomposition models. A critical component of predictive modeling is numerical inference of kinetic parameters from bench scale data. Currently, an active area of computational pyrolysis research focuses on identifying efficient, robust methods for optimization. This paper demonstrates that kinetic parameter calibration problems can successfully be solved using classical gradient-based optimization. We explore calibration examples that exhibit characteristics of concern: high nonlinearity, high dimensionality, complicated schemes, overlapping reactions, noisy data, and poor initial guesses. The examples demonstrate that a simple, non-invasive change to the problem formulation can simultaneously avoid local minima, avoid computation of derivative matrices, achieve a computational efficiency speedup of 10x, and make optimization robust to perturbations of parameter components. Techniques from the mathematical optimization and inverse problem communities are employed. By re-examining gradient-based algorithms, we highlight opportunities to develop kinetic parameter calibration methods that should outperform current methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Surface chemistry models for GaAs epitaxial growth and hydride cracking using reacting flow simulations

Hydride vapor phase epitaxy (HVPE) is a promising technology that can aid in the cost reduction of III-V materials and devices manufacturing, particularly high-efficiency solar cells for space and terrestrial applications. However, recent demonstrations of ultrafast growth rates (~500 µm/h) via uncracked hydrides are not well described by present models for the growth. Therefore, it is necessary to understand the kinetics of the growth process and its coupling with transport phenomena, so as to enable fast and uniform epitaxial growth. In this work, we derive a kinetic model using experimental data and integrate it into a computational fluid dynamics simulation of an HVPE growth reactor. We also modify an existing hydride cracking model that we validate against numerical simulations and experimental data. Here, we show that the developed growth model and the improved cracking model are able to reproduce experimental growth measurements of GaAs in an existing HVPE system.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

DOE OSTI2022

Formation of PAHs, phenol, benzofuran, and dibenzofuran in a flow reactor from the oxidation of ethylene, toluene, and n-decane

The necessity to reveal the formation mechanism of not only polycyclic aromatic hydrocarbons (PAHs) but also oxygenated PAHs (OPAHs) during combustion is increasing. Although many studies on PAHs have been conducted, fundamental studies investigating OPAH formation are still limited. Phenol, benzofuran, and dibenzofuran were selected as OPAHs in this study. We experimentally investigated fuel-rich oxidation in a flow reactor at atmospheric pressure, mean gas temperatures from 1050 to 1350 K, residence times from 0.2 to 1.5 s, and equivalence ratios from 3.0 to 12.0. Ethylene, toluene, and n-decane were used as fuels. Three kinds of OPAHs as well as 23 kinds of PAHs including monocyclic structures were quantitatively measured through a direct sampling using gas chromatography mass spectrometry. The results showed that concentrations of PAHs and OPAHs were strongly affected by the fuel type, with increasing concentrations in the following order: ethylene < n-decane < toluene. A chemical kinetic model for OPAH formation was developed based on a recent PAH growth model. The model showed that the predicted concentrations of OPAHs and PAHs were in reasonable agreement with the measured data in this study. Some modifications were made to the previous model based on recent literature studies and on the comparison of the simulated and measured results. The effect of the fuel type on the formation of PAHs and OPAHs was investigated through kinetic analysis using the model to discuss their reaction pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Fuel-rich oxidation of gasoline surrogate components in an atmospheric flow reactor

Fuel-rich oxidation of three typical gasoline surrogate components, toluene, isooctane, and n-heptane, was investigated in an atmospheric-pressure flow reactor at mean gas temperatures from 1050 to 1350 K, equivalence ratio of 9.0, and residence times of 0.45 and 1.2 s. Not only polycyclic aromatic hydrocarbons (PAHs) up to 3 ring structure but also small intermediate products from C 1 to C 5 were quantified by a gas chromatograph mass spectrometry coupled with photon ionization and gas chromatograph with flame ionization detector, respectively. The kinetic model recently developed by Lawrence Livermore National Laboratory was revised to reflect the results of many recent investigations. Basically, the updated model could satisfactorily reproduce the experimental mole fractions of many species. The experimental and simulated results showed that PAH mole fractions produced were in the order of toluene, isooctane, and n-heptane. The kinetic analysis using the model was carried out to explore PAH formation pathways, especially focusing on naphthalene, acenaphthylene, and phenanthrene. Through rate of production analysis, it was found that the main formation pathways of many PAHs were affected by the fuels. Although resonantly stabilized radicals, such as benzyl and fulvenallenyl radicals, played a crucial role in the formation pathways of many PAHs in every fuel, they were produced through hydrogen elimination of toluene in toluene fuel, while they were formed from small intermediate products in isooctane and n-heptane fuels. Sensitivity analysis revealed the difference and similarity of the reactions with large positive coefficients according to the fuels studied here. In conclusion, the molecular growth reactions of aromatic species were influential in PAH production in every fuel, whereas the ring formation reactions from small species and the reactions involving toluene had large positive sensitivity coefficients in isooctane/n-heptane and toluene fuels, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2019

Experimental data based combinatorial kinetic simulations for predictions of synergistic catalyst mixtures

Physical mixtures of catalysts can sometimes show synergistic activity which exceeds the sum of the activities of the same amount of catalysts employed separately. However, identifying such synergistic mixtures is non-trivial. Not all mixtures display synergy, and the number of combinations that are possible (even for binary mixtures of equal portions) scales very rapidly, and thus they would be costly to screen experimentally. In this work, we show that it is possible to predict synergistic mixtures using combinatorial kinetic simulations based on experimental data collected on individual catalysts. The data was first collected for conditions-of-interest under low conversions (also called near differential conditions, such that each condition approximates a small volume in a non-differential conditions reactor) to build a library of kinetic models (one model for each catalyst). This data was then used for combinatorial kinetic simulations of non-differential reactor scale conversions, and successfully predicted the qualitative behavior of two synergistic physical mixtures. The capability was utilized in the present work in the context of converting CO, C 3 H 6 , and NO species in simulated car exhaust, but is a general approach. We provide equations to estimate the costs of the screening method as well as the combinatorial kinetics method and show that the costs of identifying synergistic physical mixtures become much lower with the combinatorial kinetics method even before ten catalysts of interest. The cost of using prediction by the combinatorial kinetics method continue to become lower (relative to screening) when more catalysts are added -- and even more so if mixtures beyond two components, or consideration of multiple catalyst ratios, is of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

The Argon Power Cycle, a zero emission, flexible and reliable power generation technology

The proposed project was aimed at assisting in the development of the Argon Power Cycle (APC), a zero emission, flexible, and efficient power generation technology. In the effort of bringing the APC to the market, NTS proposed in the scope of this project several partial but key activities to advance the development of the APC. These tasks encompassed 1) the analysis of the market, developing a techno-economic model with integrated market inputs, 2) further optimization of our chemical kinetics models for more accurate simulation of the process and 3) scale up of our original prototype to a larger, commercially viable platform. NTS expanded on this scope of work to secure additional R&D funding in amounts north of $7M, carry out customer discovery with the sponsorship of the NSF and grow the team’s engineering talent.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS2002

Origins of Solar Systems Program

O'D. Alexander, Conel↗

NASA NTRS2003

Final Summary of Research Report to the National Aeronautics and Space Administration Origins of Solar Systems Program

The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments. Cosmic spherules present an opportunity to test our evaporation models under flash heating conditions that would be difficult to simulate experimentally. However, there is surprisingly little known about the isotopic compositions of silicate cosmic spherules, and a number of questions need to be addressed. Is the range of compositions they exhibit due to evaporation? If they are, are the relative volatilities consistent with the models/experiments and are the isotopic fractionations consistent with Rayleigh conditions? For instance, do the alkalis and S evaporate prior to significant melting so that conditions did not meet the Rayleigh criteria of rapid diffusion? If so, their isotopic fractionation might be considerably suppressed. Could this mechanism of K loss apply to chondrule formation? The Fe isotopic fractionation during evaporation of silicates has not been measured, so cosmic spherules might provide a clue to whether FeO diffusion is fast enough to maintain Rayleigh conditions during evaporation. And so on.

O'D. Alexander, Conel↗

NASA NTRS2013

NASA Models of Space Radiation Induced Cancer, Circulatory Disease, and Central Nervous System Effects

The risks of late effects from galactic cosmic rays (GCR) and solar particle events (SPE) are potentially a limitation to long-term space travel. The late effects of highest concern have significant lethality including cancer, effects to the central nervous system (CNS), and circulatory diseases (CD). For cancer and CD the use of age and gender specific models with uncertainty assessments based on human epidemiology data for low LET radiation combined with relative biological effectiveness factors (RBEs) and dose- and dose-rate reduction effectiveness factors (DDREF) to extrapolate these results to space radiation exposures is considered the current "state-of-the-art". The revised NASA Space Risk Model (NSRM-2014) is based on recent radio-epidemiology data for cancer and CD, however a key feature of the NSRM-2014 is the formulation of particle fluence and track structure based radiation quality factors for solid cancer and leukemia risk estimates, which are distinct from the ICRP quality factors, and shown to lead to smaller uncertainties in risk estimates. Many persons exposed to radiation on earth as well as astronauts are life-time never-smokers, which is estimated to significantly modify radiation cancer and CD risk estimates. A key feature of the NASA radiation protection model is the classification of radiation workers by smoking history in setting dose limits. Possible qualitative differences between GCR and low LET radiation increase uncertainties and are not included in previous risk estimates. Two important qualitative differences are emerging from research studies. The first is the increased lethality of tumors observed in animal models compared to low LET radiation or background tumors. The second are Non- Targeted Effects (NTE), which include bystander effects and genomic instability, which has been observed in cell and animal models of cancer risks. NTE's could lead to significant changes in RBE and DDREF estimates for GCR particles, and the potential effectiveness of radiation mitigator's. The NSRM- 2014 approaches to model radiation quality dependent lethality and NTE's will be described. CNS effects include both early changes that may occur during long space missions and late effects such as Alzheimer's disease (AD). AD effects 50% of the population above age 80-yr, is a degenerative disease that worsens with time after initial onset leading to death, and has no known cure. AD is difficult to detect at early stages and the small number of low LET epidemiology studies undertaken have not identified an association with low dose radiation. However experimental studies in mice suggest GCR may lead to early onset AD. We discuss modeling approaches to consider mechanisms whereby radiation would lead to earlier onset of occurrence of AD. Biomarkers of AD include amyloid beta (A(Beta)) plaques, and neurofibrillary tangles (NFT) made up of aggregates of the hyperphosphorylated form of the micro-tubule associated, tau protein. Related markers include synaptic degeneration, dentritic spine loss, and neuronal cell loss through apoptosis. Radiation may affect these processes by causing oxidative stress, aberrant signaling following DNA damage, and chronic neuroinflammation. Cell types to be considered in multi-scale models are neurons, astrocytes, and microglia. We developed biochemical and cell kinetics models of DNA damage signaling related to glycogen synthase kinase-3(Beta) (GSK3(Beta)) and neuroinflammation, and considered multi-scale modeling approaches to develop computer simulations of cell interactions and their relationships to A(Beta) plaques and NFTs. Comparison of model results to experimental data for the age specific development of A(Beta) plaques in transgenic mice will be discussed.

Cucinotta, Francis A.↗

DOE OSTI2021

A priori examination of reduced chemistry models derived from canonical stirred reactors using three-dimensional direct numerical simulation datasets

Data-driven approaches to construct reduced chemical kinetic models, that rely heavily on thermo-chemical datasets with full chemical kinetics, have been gaining popularity. Datasets from direct numerical simulations (DNS) under three-dimensional (3-D) realistic turbulent flow conditions are desirable but limited to carefully designed parametric conditions due to the computational cost. Constructing datasets from a large ensemble of zero-dimensional stirred reactors like perfectly stirred reactor (PSR) and partially stirred reactor (PaSR) is a computationally efficient solution to consider the turbulence-chemistry interactions and cover a broad range of parametric conditions. In this paper, we derive reduced chemistry models from solutions of a large number of PSR and PaSR reactors using autoencoder (AE) neural networks and principal component analysis (PCA), and conduct a priori examination of the reduced models in three temporally evolving 3-D DNS jet flames featuring local extinction and re-ignition. The results show that the reduced models derived from PaSR datasets, i.e., AE-PaSR and PCA-PaSR, generally show significant improvement over the ones derived from PSR datasets. Among all the reduced models, AE-PaSR shows the best agreement with DNS results on the reconstruction accuracy and the representation of temporally evolving local extinction and re-ignition events.

Zhang, Pei↗

DOE OSTI2024

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Experimental and computational investigations of ethane and ethylene kinetics with copper oxide particles for Chemical Looping Combustion

In this work, reaction pathways for the oxidation of methane, ethane, and ethylene with CuO was obtained by ReaxFF Molecular Dynamics (MD) simulations between temperatures of 1000 K and 2000 K. Experiments in a fixed-bed flow reactor were preformed with methane, ethane, and ethylene at temperatures ranging from 500 K to 1000 K with time-dependent species measurements from an Electron-Ionization Molecular Beam Mass Spectrometer (MBMS), and species validation with Gas Chromatography (GC) for detection of complete and intermediate combustion products. The MBMS and GC allow for the detection of oxygenated species and larger species produced from radical reformation. The simulation and experiment agree on the production of such species as CH 3 CHO, CH 2 O, CO, and H 2 O, which allow for the creation of simple C1 and C2 reaction pathways, which can be used in kinetic models of C2 species and larger fuels such as biofuels, which inherently depend on C1 and C2 kinetics and reaction pathway. The simulation and experiment disagree on the formation of C 2 H 2 , CH 3 OH, and CO 2 with large amounts of C 2 H 2 being measured in the ethylene oxidation simulations and CH 3 OH being formed in methane oxidation simulations, while neither species were experimentally found. In the case of CO 2 large amounts of CO 2 are rapidly produced in experiments with C2 fuels at 800 K, while little-to-no CO 2 was observed in simulations. This is believed to be resulting from the extremely short timescale of the simulations, preventing total oxidation of the fuel. Here, the differences in products produced between simulation and experiment allow for the potential to modify the ReaxFF potential functions to more accurately model the experimental products of Cu–H–O–C reaction kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Mechanistic modelling of SARS-CoV-2 and other infectious diseases and the effects of therapeutics

Modern viral kinetic modeling and its application to therapeutics is a field that attracted the attention of the medical, pharmaceutical and modeling communities during the early days of the AIDS epidemic. Its successes lead to applications of modeling methods not only to HIV but a plethora of other viruses, such as hepatitis C virus (HCV), hepatitis B virus (HBV) and cytomegalovirus (CMV), which along with HIV cause chronic diseases, and viruses such as influenza, respiratory syncytial virus (RSV), West Nile virus (WNV), Zika virus, and SAR-CoV-2, which generally cause acute infections. Here we first review the historical development of mathematical models to understand HIV and HCV infections and the effects of treatment by fitting the models to clinical data. Furthermore, we then focus on recent efforts and contributions of applying these models towards understanding SARS-CoV-2 infection and highlight outstanding questions where modeling can provide crucial insights and help to optimize non-pharmaceutical and pharmaceutical interventions of the COVID-19 pandemic. The review is written from our personal perspective emphasizing the power of simple target cell limited models that provided important insights and then their evolution into more complex models that captured more of the virology and immunology. To quote Albert Einstein “Everything should be made as simple as possible, but not simpler”, and this idea underlies the modeling we describe below.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗