SEARCH · Engineering Papers
Results for “chemical dynamics”
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.
Real-space observation of the dissociation of a transition metal complex and its concurrent energy redistribution
Mechanistic insights into photodissociation dynamics of transition metal carbonyls, like Fe(CO) 5 , are fundamental for understanding active catalytic intermediates. Although extensively studied, the structural dynamics of these systems remain elusive. Using ultrafast X-ray scattering, we uncover the photochemistry of Fe(CO) 5 in real space and time, observing synchronous oscillations in atomic pair distances, followed by a prompt rotating CO release preferentially in the axial direction. This behavior aligns with simulations, reflecting the interplay between the axial Fe-C distances’ potential energy landscape and non-adiabatic transitions between metal-to-ligand charge-transfer states. Additionally, we characterize a secondary delayed CO release associated with a reduction of Fe-C steady state distances and structural dynamics of the formed Fe(CO) 4 . Our results quantify energy redistribution across vibration, rotation, and translation degrees of freedom, offering a microscopic view of complex structural dynamics, enhancing our grasp on Fe(CO) 5 photodissociation, and advancing our understanding of transition metal catalytic systems.
Stochastic self-tuning hybrid algorithm for reaction-diffusion systems
Many biochemical phenomena involve reactants with vastly different concentrations, some of which are amenable to continuum-level descriptions, while the others are not. We present a hybrid self-tuning algorithm to model such systems. The method combines microscopic (Brownian) dynamics for diffusion with mesoscopic (Gillespie-type) methods for reactions and remains efficient in a wide range of regimes and scenarios with large variations of concentrations. Furthermore, its accuracy, robustness, and versatility are balanced by redefining propensities and optimizing the mesh size and time step. We use a bimolecular reaction to demonstrate the potential of our method in a broad spectrum of scenarios: from almost completely reaction-dominated systems to cases where reactions rarely occur or take place very slowly. The simulation results show that the number of particles present in the system does not degrade the performance of our method. This makes it an accurate and computationally efficient tool to model complex multireaction systems.
Is there a single origin of life?
The emergence of the first life on the earth is now established as an early event, closely related to the evolving earth. Lack of knowledge of the primitive terrestrial conditions contemporary with the evolving prebolic organic chemistry limits reconstruction techniques. The primitive earth's aqueous history is essential to unraveling this problem. Based on current knowledge of other planets of the solar system, no close analog to the early earth can be expected. It is still not known if there was a second origin or if only earth has life. This may depend upon the question of the survival of information-bearing chemical systems in a dynamic or chaotic environment and the chemical protection afforded within such a system. Water is the central molecule of controversy: the blessing and the curse of the chemist. New and novel chemical mechanisms and systems abound.
Fluid dynamics and low gravity effects of chemical vapor deposition
Based on the comparison between experimental data and numerical results for the growth of GaAs from TMGa, it was shown that 3D simulations are necessary to simulate rectangular CVD reactors even when operated under subcritical (Ra) conditions. The important points found are summarized in the three attached reprints. The experimental studies of mixed convection in horizontal channels have shown three regimes of high Ra (22,220) number flows. At Re = 18.5, the rolls develop very quickly, significantly modulating the axial velocity even before it reaches the beginning of the hot plate. A few centimeters downstream, the velocities become asymmetric about the vertical centerplane and at x = 12 cm, become unsteady. These asymmetries were predicted theoretically, but experimental evidence has not been published prior to this work. At Re = 36, the axial velocity is only slightly modified at x = 0. Although the flow remains steady and symmetric about the vertical centerplane, there is a small spatial oscillation in the velocities over the length of the channel. The period of this oscillation was around 5 cm. At Re = 54, the longitudinal rolls developed smoothly over a length of 30 cm, with no asymmetries, unsteadiness, or spatial oscillations. Comparison of numerical simulations of these flows to experiments has revealed the importance and difficulty of setting proper thermal boundary conditions on the sidewalls. Calculated flows and experimentally measured flows showed very similar profiles, but at different axial locations, with the rolls developing more rapidly in the experiments. This is directly attributable to partially conducting sidewalls of the apparatus being hotter in the entrance section than the adiabatic walls of the simulations. A thorough comparison of the experimental data and numerical results for a variety of sidewall boundary conditions is in preparation.
New instrumentation and tools to study dynamic (non-equilibrium) physico-chemical interfacial processes across large length, rates, and time scales (Final Report)
The hydrophobic interaction dictates the behavior of many systems relevant to biological and industrial processes. However, the origin of the hydrophobic interaction has been extensively debated. In this article, we demonstrate the use of nanostructured surfaces to control the range of the hydrophobic interaction. We measure the strength and range of the hydrophobic interaction using bubble probe force microscopy, and interpret our results using the Stokes–Reynolds–Young–Laplace model. Our findings give insights into the mechanisms of the hydrophobic interaction and will lead to better understanding and control of diverse colloidal and interfacial phenomena.
Quantum dynamical effects of vibrational strong coupling in chemical reactivity
Abstract Recent experiments suggest that ground state chemical reactivity can be modified when placing molecular systems inside infrared cavities where molecular vibrations are strongly coupled to electromagnetic radiation. This phenomenon lacks a firm theoretical explanation. Here, we employ an exact quantum dynamics approach to investigate a model of cavity-modified chemical reactions in the condensed phase. The model contains the coupling of the reaction coordinate to a generic solvent, cavity coupling to either the reaction coordinate or a non-reactive mode, and the coupling of the cavity to lossy modes. Thus, many of the most important features needed for realistic modeling of the cavity modification of chemical reactions are included. We find that when a molecule is coupled to an optical cavity it is essential to treat the problem quantum mechanically to obtain a quantitative account of alterations to reactivity. We find sizable and sharp changes in the rate constant that are associated with quantum mechanical state splittings and resonances. The features that emerge from our simulations are closer to those observed in experiments than are previous calculations, even for realistically small values of coupling and cavity loss. This work highlights the importance of a fully quantum treatment of vibrational polariton chemistry.
Autonomously Responsive Membranes for Chemical Warfare Protection
Stimuli-responsive materials offer new opportunities to resolve long-standing material challenges and are rapidly gaining pivotal roles in diverse applications. For example, smart protective garments that rapidly transport water vapor and autonomously block chemical threats are expected to enable an effective new paradigm of adaptive personal protection. However, the incorporation of these seemingly incompatible properties into a single responsive system remains elusive. Herein, a bistable membrane that can rapidly, selectively, and reversibly transition from a highly breathable state in a safe environment to a chemically protective state when exposed to organophosphate threats such as sarin is demonstrated. Dynamic response to chemical stimuli is achieved through the physical collapse of an ultrathin copolymer layer on the membrane surface, which efficiently gates transport through membrane pores composed of single-walled carbon nanotubes (SWNTs). The adoption of nanometer-wide SWNTs for ultrafast moisture conduction enables a simultaneous boost in size-sieving selectivity and water-vapor permeability by decreasing nanotube diameter, thereby overcoming the breathability/protection trade-off that limits conventional membrane materials. Adaptive multifunctional membranes based on this platform greatly extend the active use of a protective garment and present exciting opportunities in many other areas including separation processes, sensing, and smart delivery.
SABER Observations of the OH Meinel Airglow Variability Near the Mesopause
The Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument, one of four on board the TIMED satellite, observes the OH Meinel emission at 2.0 m that peaks near the mesopause. The emission results from reactions between members of the oxygen and hydrogen chemical families that can be significantly affected by mesopause dynamics. In this study we compare SABER measurements of OH Meinel emission rates and temperatures with predictions from a 3-dimensional chemical dynamical model. In general, the model is capable of reproducing both the observed diurnal and seasonal OH Meinel emission variability. The results indicate that the diurnal tide has a large effect on the overall magnitude and temporal variation of the emission in low latitudes. This tidal variability is so dominant that the seasonal cycle in the nighttime emission depends very strongly on the local time of the analysis. At higher latitudes, the emission has an annual cycle that is due mainly to transport of oxygen by the seasonally reversing mean circulation.
Non-equilibrium chemistry in the solar nebula and early solar system: Implications for the chemistry of comets
Theoretical models of solar nebula and early solar system chemistry which take into account the interplay between chemical, physical, and dynamical processes have great utility for deciphering the origin and evolution of the abundant chemically reactive volatiles (H, O, C, N, S) observed in comets. In particular, such models are essential for attempting to distinguish between presolar and solar nebula products and for quantifying the nature and duration of nebular and early solar system processing to which the volatile constituents of comets have been subjected. The diverse processes and energy sources responsible for chemical processing in the solar nebula and early solar system are discussed. The processes considered include homogeneous and heterogeneous thermochemical and photochemical reactions, and disequilibration resulting from fluid transport, condensation, and cooling whenever they occur on timescales shorter than those for chemical reactions.
A Multireference Approach to Electron and Electron–Nuclear Dynamics in Nanomaterials (Final Report)
Many important chemical and physical phenomena involve dynamics on large number of electronic states. Thus, there is a critical need to develop methods to simulate dynamics in dense manifolds of states. Towards this end, we have: a) developed the multiple cloning in dense manifolds of states (MCDMS) method, which is capable of accurately modeling the quantum mechanical coherence between populations on a large number of electronic states, b) implemented MCDMS into the free, open-source PySpawn software package, c) developed graphics processing unit-accelerated algorithms modeling electron dynamics in light fields via Floquet time-dependent configuration interaction (F-TDCI), and d) critically compared different orbital bases in order to achieve an accurate and efficient F-TDCI expansion. This grant ended in August 2020, when our group moved to from Michigan State University to Stony Brook University, where this project continues under grant number DE-SC0021643.
Global impact of the Antarctic ozone hole - Dynamical dilution with a three-dimensional chemical transport model
The hypothesis that, once the Antarctic ozone hole is photochemically generated in the spring, the ozone deficit remains in the stratosphere throughout the year was investigated using a three-dimensional chemical transport model (CTM) for the stratosphere with linearized ozone chemistry. An Antarctic ozone hole was imposed on the basis of observations, and the dynamical dilution of the ozone-depleted air was followed for the next 12 months. The CTM calculations show that the transport of ozone-poor air from the Antarctic vortex results in measurable decreases to column ozone, extending to 30 deg S during the austral summer. By the end of the year, about 70 percent of the initially prescribed ozone deficit is replenished through stratospheric chemistry.
Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials
The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.
Graph characterization of higher-order structure in atmospheric chemical reaction mechanisms
Atmospheric chemical reactions play an important role in air quality and climate change. While the structure and dynamics of individual chemical reactions are fairly well understood, the emergent properties of the entire atmospheric chemical system, which can involve many different species that participate in many different reactions, are not well described. In this work, we leverage graph-theoretic techniques to characterize patterns of interaction (“motifs”) in three different representations of gas-phase atmospheric chemistry, termed “chemical mechanisms.” These widely used mechanisms, the master chemical mechanism, the GEOS-Chem mechanism, and the Super-Fast mechanism, vary dramatically in scale and application, but they all generally aim to simulate the abundance and variability of chemical species in the atmosphere. This motif analysis quantifies the fundamental patterns of interaction within the mechanisms, which are directly related to their construction. For example, the gas-phase chemistry in the very small Super-Fast mechanism is entirely composed of bimolecular reactions, and its motif distribution matches that of an individual bimolecular reaction well. The larger and more complex mechanisms show emergent motif distributions that differ strongly from any specific reaction type, consistent with their complexity. The proposed motif analysis demonstrates that while these mechanisms all have a similar design goal, their higher-order structure of interactions differs strongly and thus provides a novel set of tools for exploring differences across chemical mechanisms.
Chemical Ecosystem Selection on Mineral Surfaces Reveals Long-Term Dynamics Consistent with the Spontaneous Emergence of Mutual Catalysis
How did chemicals first become organized into systems capable of self-propagation and adaptive evolution? One possibility is that the first evolvers were chemical ecosystems localized on mineral surfaces and composed of sets of molecular species that could catalyze each other’s formation. We used a bottom-up experimental framework; chemical ecosystem selection (CES) to evaluate this perspective and search for surface-associated; and mutually catalytic chemical systems based on the changes in chemistry they are expected to induce. Here, we report the results of preliminary CES experiments conducted using a synthetic “prebiotic soup” and pyrite grains, which yield dynamical patterns that are suggestive of the emergence of mutual catalysis. While more research is needed to better understand the specific patterns observed here and determine whether they are reflective of self-propagation, these results illustrate the potential power of CES to test competing hypotheses for the emergence of protobiological chemical systems.
Advanced Multi-phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis
A Navier-Stokes code, finite difference Navier-Stokes (FDNS), is used to analyze the complicated internal flowfield of the SRM (solid rocket motor) to explore the impacts due to the effects of chemical reaction, particle dynamics, and slag accumulation on the solid rocket motor (SRM). The particulate multi-phase flowfield with chemical reaction, particle evaporation, combustion, breakup, and agglomeration models are included in present study to obtain a better understanding of the SRM design. Finite rate chemistry model is applied to simulate the chemical reaction effects. Hermsen correlation model is used for the combustion simulation. The evaporation model introduced by Spalding is utilized to include the heat transfer from the particulate phase to the gase phase due to the evaporation of the particles. A correlation of the minimum particle size for breakup expressed in terms of the Al/Al2O3 surface tension and shear force was employed to simulate the breakup of particles. It is assumed that the breakup occurs when the Weber number exceeds 6. A simple L agglomeration model is used to investigate the particle agglomeration. However, due to the large computer memory requirements for the agglomeration model, only 2D cases are tested with the agglomeration model. The VOF (Volume of Fluid) method is employed to simulate the slag buildup in the aft-end cavity of the redesigned solid rocket motor (RSRM). Monte Carlo method is employed to calculate the turbulent dispersion effect of the particles. The flowfield analysis obtained using the FDNS code in the present research with finite rate chemical reaction, particle evaporation, combustion, breakup, agglomeration, and VOG models will provide a design guide for the potential improvement of the SRM including the use of materials and the shape of nozzle geometry such that a better performance of the SRM can be achieved. The simulation of the slag buildup in the aft-end cavity can assist the designer to improve the design of the RSRM geometry.
Models and Strategies for Optimal Demand Side Management in the Chemical Industries
Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.
PubChemLite Plus Collision Cross Section (CCS) Values for Enhanced Interpretation of Nontarget Environmental Data
Finding relevant chemicals in the vast (known) chemical space is a major challenge for environmental and exposomics studies leveraging nontarget high resolution mass spectrometry (NT-HRMS) methods. Chemical databases now contain hundreds of millions of chemicals, yet many are not relevant. This article details an extensive collaborative, open science effort to provide a dynamic collection of chemicals for environmental, metabolomics, and exposomics research, along with supporting information about their relevance to assist researchers in the interpretation of candidate hits. The PubChemLite for Exposomics collection is compiled from ten annotation categories within PubChem, enhanced with patent, literature and annotation counts, predicted partition coefficient (logP) values, as well as predicted collision cross section (CCS) values using CCSbase. Monthly versions are archived on Zenodo under a CC-BY license, supporting reproducible research, and a new interface has been developed, including historical trends of patent and literature data, for researchers to browse the collection. This article details how PubChemLite can support researchers in environmental and exposomics studies, describes efforts to increase the availability of experimental CCS values, and explores known limitations and potential for future developments. The data and code behind these efforts are openly available.