Search NASA⌕ Search

SEARCH · Search NASA

Results for “biochemical simulations”

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

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

At least 181 records · Page 10

Proceedings of the 1973 Lyndon B. Johnson Space Center Endocrine Program Conference

Papers given at the conference are presented. Subjects covered include the following: biochemical changes during 28 days of space flight, modulating the pituitary-adrenal response to stress, the significance of biorhythms in space flight, the importance of the rein-angiotensin system in normal cardiovascular homeostasis, a progress report of stress-induced changes in corticosteroid metabolism, recent studies of physiological factors involved in the regulation of serotonin content and turnover in the brain, the role of brain biogenic amines in the control of pituitary-adrenocortical activity, application of the water immersion model to man by studies of acid-base homeostasis during simulated weightlessness, the present status of physiological studies and analysis of calcium homeostasis in the Apollo and Skylab programs, and endocrine considerations in the red-cell-mass and plasma-volume changes of Skylab 2 and 3 crews.

Source record↗

Elucidating interprotein energy transfer dynamics within the antenna network from purple bacteria

In photosynthesis, absorbed light energy transfers through a network of antenna proteins with near-unity quantum efficiency to reach the reaction center, which initiates the downstream biochemical reactions. While the energy transfer dynamics within individual antenna proteins have been extensively studied over the past decades, the dynamics between the proteins are poorly understood due to the heterogeneous organization of the network. Previously reported timescales averaged over such heterogeneity, obscuring individual interprotein energy transfer steps. Here, we isolated and interrogated interprotein energy transfer by embedding two variants of the primary antenna protein from purple bacteria, light-harvesting complex 2 (LH2), together into a near-native membrane disc, known as a nanodisc. We integrated ultrafast transient absorption spectroscopy, quantum dynamics simulations, and cryogenic electron microscopy to determine interprotein energy transfer timescales. By varying the diameter of the nanodiscs, we replicated a range of distances between the proteins. The closest distance possible between neighboring LH2, which is the most common in native membranes, is 25 Å and resulted in a timescale of 5.7 ps. Larger distances of 28 to 31 Å resulted in timescales of 10 to 14 ps. Corresponding simulations showed that the fast energy transfer steps between closely spaced LH2 increase transport distances by ∼15%. Overall, our results introduce a framework for well-controlled studies of interprotein energy transfer dynamics and suggest that protein pairs serve as the primary pathway for the efficient transport of solar energy.

59 BASIC BIOLOGICAL SCIENCES↗

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. Overall, this protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

59 BASIC BIOLOGICAL SCIENCES↗

Mass balances for a biological life support system simulation model

Design decisions to aid the development of future space based biological life support systems (BLSS) can be made with simulation models. The biochemistry stoichiometry was developed for: (1) protein, carbohydrate, fat, fiber, and lignin production in the edible and inedible parts of plants; (2) food consumption and production of organic solids in urine, feces, and wash water by the humans; and (3) operation of the waste processor. Flux values for all components are derived for a steady state system with wheat as the sole food source. The large scale dynamics of a materially closed (BLSS) computer model is described in a companion paper. An extension of this methodology can explore multifood systems and more complex biochemical dynamics while maintaining whole system closure as a focus.

Volk, Tyler↗

Behavioral and Physiological Effects of Hindlimb Unloading in Rats

The overarching objective of this project was to identify changes in neural and biochemical systems of the central and peripheral nervous systems (the CNS and PNS) that are related to disruptions of functional motor responses, or motor control. The identification of neural and biochemical changes that are related to sensory-motor adaptation elicited as animals react to changes in the gravitational field was of particular interest. Thus, the major objective of this work was to study disruptions of motor responses that arise after (sic. due to) chronic exposure to altered gravity (G). To do this, parallel studies investigating changes in neural, sensory, and neuromuscular systems were conducted after animals (rats) experienced chronic exposure to conditions of altered-G. Conditions of altered-G included hyper-G produced by centrifugation, micro-G produced by orbital flight, and simulated micro-G produced by hind limb suspension. A second major interest was to examine the contribution of putative changes in sensory systems to disruptions of motor responses. To do this, motor responses and reflexes of rats were studied following chronic treatment with streptomycin sulfate (STP, an ototoxic chemical) to damage the vestibular hair cells.

Fox, Robert A.↗

Neonatal rat heart cells cultured in simulated microgravity

In vitro characteristics of cardiac cells cultured in simulated microgravity are reported. Tissue culture methods performed at unit gravity constrain cells to propagate, differentiate, and interact in a two-dimensional (2D) plane. Neonatal rat cardiac cells in 2D culture organize predominantly as bundles of cardiomyocytes with the intervening areas filled by nonmyocyte cell types. Such cardiac cell cultures respond predictably to the addition of exogenous compounds, and in many ways they represent an excellent in vitro model system. The gravity-induced 2D organization of the cells, however, does not accurately reflect the distribution of cells in the intact tissue. We have begun characterizations of a three-dimensional (3D) culturing system designed to mimic microgravity. The NASA-designed High-Aspect Ratio Vessel (HARV) bioreactors provide a low shear environment that allows cells to be cultured in static suspension. HARV-3D cultures were prepared on microcarrier beads and compared to control-2D cultures using a combination of microscopic and biochemical techniques. Both systems were uniformly inoculated and medium exchanged at standard intervals. Cells in control cultures adhered to the polystyrene surface of the tissue culture dishes and exhibited typical 2D organization. Cells cultured in HARVs adhered to microcarrier beads, the beads aggregated into defined clusters containing 8 to 15 beads per cluster, and the clusters exhibited distinct 3D layers: myocytes and fibroblasts appeared attached to the surfaces of beads and were overlaid by an outer cell type. In addition, cultures prepared in HARVs using alternative support matrices also displayed morphological formations not seen in control cultures. Generally, the cells prepared in HARV and control cultures were similar; however, the dramatic alterations in 3D organization recommend the HARV as an ideal vessel for the generation of tissuelike organization of cardiac cells in vitro.

Non-NASA Center↗

Integrative teaching of metabolic modeling and flux analysis with interactive python modules

Abstract The modeling of rates of biochemical reactions—fluxes—in metabolic networks is widely used for both basic biological research and biotechnological applications. A number of different modeling methods have been developed to estimate and predict fluxes, including kinetic and constraint‐based (Metabolic Flux Analysis and flux balance analysis) approaches. Although different resources exist for teaching these methods individually, to‐date no resources have been developed to teach these approaches in an integrative way that equips learners with an understanding of each modeling paradigm, how they relate to one another, and the information that can be gleaned from each. We have developed a series of modeling simulations in Python to teach kinetic modeling, metabolic control analysis, 13C‐metabolic flux analysis, and flux balance analysis. These simulations are presented in a series of interactive notebooks with guided lesson plans and associated lecture notes. Learners assimilate key principles using models of simple metabolic networks by running simulations, generating and using data, and making and validating predictions about the effects of modifying model parameters. We used these simulations as the hands‐on computer laboratory component of a four‐day metabolic modeling workshop and participant survey results showed improvements in learners' self‐assessed competence and confidence in understanding and applying metabolic modeling techniques after having attended the workshop. The resources provided can be incorporated in their entirety or individually into courses and workshops on bioengineering and metabolic modeling at the undergraduate, graduate, or postgraduate level.

Kaste, Joshua A. M.↗

Technology Case Study: Economic, Sustainability, and Deployment Considerations for Sustainable Aviation Fuels Produced via Lignocellulosic Sugar Catalysis

This report presents a technology case study reflecting one exemplary representative pathway for the conversion of lignocellulosic sugars to sustainable aviation fuels (SAF) via aqueous phase reforming (APR) catalysis, considered within a broader integrated biorefinery framework based on biochemical processing operations. While far from the only option for converting sugars to SAF, this pathway was selected as a case study here based on its relatively high technology maturity and simplistic processing approach (avoiding complex separations or other equipment scalability challenges), coupled with the potential for high fuel yields and favorable costs/carbon intensities with opportunities for further near-term optimization. The report considers key process integration and engineering design considerations for a modeled hypothetical, nth-plant commercial biorefinery, reflecting a number of processing options and parameters envisioned to be achievable as future goals. Resultant outputs from Aspen Plus process simulations are evaluated through techno-economic and life cycle analyses (TEA and LCA), including implications for marginal cost of CO2 abatement and inclusion of currently-applicable policy incentives. Moving beyond base case configurations, a number of alternative scenarios are also evaluated for their ability to further improve economics, greenhouse gas (GHG) emissions, and marginal cost of abatement, highlighting a path to achieve deep decarbonization goals of more than 70% GHG reduction for SAF (with the potential to reach net-negative carbon intensities in some cases) under reasonable fuel production costs. The report also highlights future opportunities and gaps for further research on this technology pathway.

09 BIOMASS FUELS↗

Base pairing, structural and functional insights into N 4 -methylcytidine (m 4 C) and N 4 ,N 4 -dimethylcytidine (m 4 2 C) modified RNA

The N 4 -methylation of cytidine (m 4 C and m 4 2 C) in RNA plays important roles in both bacterial and eukaryotic cells. In this work, we synthesized a series of m 4 C and m 4 2 C modified RNA oligonucleotides, conducted their base pairing and bioactivity studies, and solved three new crystal structures of the RNA duplexes containing these two modifications. Our thermostability and X-ray crystallography studies, together with the molecular dynamic simulation studies, demonstrated that m 4 C retains a regular C:G base pairing pattern in RNA duplex and has a relatively small effect on its base pairing stability and specificity. By contrast, the m 4 2 C modification disrupts the C:G pair and significantly decreases the duplex stability through a conformational shift of native Watson-Crick pair to a wobble-like pattern with the formation of two hydrogen bonds. This double-methylated m 4 2 C also results in the loss of base pairing discrimination between C:G and other mismatched pairs like C:A, C:T and C:C. The biochemical investigation of these two modified residues in the reverse transcription model shows that both mono- or di-methylated cytosine bases could specify the C:T pair and induce the G to T mutation using HIV-1 RT. In the presence of other reverse transcriptases with higher fidelity like AMV-RT, the methylation could either retain the normal nucleotide incorporation or completely inhibit the DNA synthesis. These results indicate the methylation at N 4 -position of cytidine is a molecular mechanism to fine tune base pairing specificity and affect the coding efficiency and fidelity during gene replication.

59 BASIC BIOLOGICAL SCIENCES↗

Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism

Abstract The properties of the water network in concentrated HCl acid pools in nanometer‐sized reverse nonionic micelles were probed with TeraHertz absorption, dielectric relaxation spectroscopy, and reactive force field simulations capable of describing proton hopping mechanisms. We identify that only at a critical micelle size of W 0 =9 do solvated proton complexes form in the water pool, accompanied by a change in mechanism from Grotthuss forward shuttling to one that favors local oscillatory hopping. This is due to a preference for H + and Cl − ions to adsorb to the micelle interface, together with an acid concentration effect that causes a “traffic jam” in which the short‐circuiting of the hydrogen‐bonding motif of the hydronium ion decreases the forward hopping rate throughout the water interior even as the micelle size increases. These findings have implications for atmospheric chemistry, biochemical and biophysical environments, and energy materials, as transport of protons vital to these processes can be suppressed due to confinement, aggregation, and/or concentration.

Adams, Ellen M.↗

Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism

Abstract The properties of the water network in concentrated HCl acid pools in nanometer‐sized reverse nonionic micelles were probed with TeraHertz absorption, dielectric relaxation spectroscopy, and reactive force field simulations capable of describing proton hopping mechanisms. We identify that only at a critical micelle size of W 0 =9 do solvated proton complexes form in the water pool, accompanied by a change in mechanism from Grotthuss forward shuttling to one that favors local oscillatory hopping. This is due to a preference for H + and Cl − ions to adsorb to the micelle interface, together with an acid concentration effect that causes a “traffic jam” in which the short‐circuiting of the hydrogen‐bonding motif of the hydronium ion decreases the forward hopping rate throughout the water interior even as the micelle size increases. These findings have implications for atmospheric chemistry, biochemical and biophysical environments, and energy materials, as transport of protons vital to these processes can be suppressed due to confinement, aggregation, and/or concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computer Code for Nanostructure Simulation

Due to their small size, nanostructures can have stress and thermal gradients that are larger than any macroscopic analogue. These gradients can lead to specific regions that are susceptible to failure via processes such as plastic deformation by dislocation emission, chemical debonding, and interfacial alloying. A program has been developed that rigorously simulates and predicts optoelectronic properties of nanostructures of virtually any geometrical complexity and material composition. It can be used in simulations of energy level structure, wave functions, density of states of spatially configured phonon-coupled electrons, excitons in quantum dots, quantum rings, quantum ring complexes, and more. The code can be used to calculate stress distributions and thermal transport properties for a variety of nanostructures and interfaces, transport and scattering at nanoscale interfaces and surfaces under various stress states, and alloy compositional gradients. The code allows users to perform modeling of charge transport processes through quantum-dot (QD) arrays as functions of inter-dot distance, array order versus disorder, QD orientation, shape, size, and chemical composition for applications in photovoltaics and physical properties of QD-based biochemical sensors. The code can be used to study the hot exciton formation/relation dynamics in arrays of QDs of different shapes and sizes at different temperatures. It also can be used to understand the relation among the deposition parameters and inherent stresses, strain deformation, heat flow, and failure of nanostructures.

Filikhin, Igor↗

Quantifying Impacts of Biomass Pelletization on Fast Pyrolysis Using a Single-Particle Reactor, X-ray Computed Tomography, and Computational Modeling

The pore structure and density of lignocellulosic feedstocks dictate intraparticle transport phenomena and thereby play an important role in thermochemical conversion processes such as fast pyrolysis for biofuel and biochemical production. Variations in microstructure are inherent from different biomass species and can be introduced by preprocessing techniques such as cutting and pelletization. Morphological changes also occur during conversion and lead to vastly different pore structures and behavior during pyrolysis, which impact required conversion times and product distributions. The current work presents a comprehensive comparison of fast pyrolysis of neat and pelletized pine feedstocks, which includes single-particle experiments, modeling, and 3D imaging by X-ray computed tomography (XCT). The particle-scale model included anisotropic heat and mass transport in a shrinking particle with pyrolysis reactions based on the CRECK mechanism with boundary conditions informed by reactor-scale simulations of the single-particle reactor. The models were validated by measurements of the temperature and mass loss from single-particle pyrolysis experiments of neat and pelletized pine. Quantitative analysis of XCT geometries revealed that pyrolytic conversion yielded chars with increased porosity and permeability compared to the unpyrolyzed materials, along with decreased tortuosity and anisotropy. Pelletization of the pine feedstock resulted in a much denser, less permeable material, which converted slower and produced more residual char after pyrolysis compared to neat pine. The results from particle modeling revealed that accounting for the dynamic and anisotropic heat and mass transport caused by differences in pore structure is critical to achieving agreement with experimental results. Overall, this study highlights the dramatic differences in conversion behavior imparted by pelletization and the importance of capturing microstructural attributes in computational models to guide the design and optimization of pyrolysis processes for specific biomass feedstocks.

09 BIOMASS FUELS↗

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH 2 -capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young’s modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. Here, we construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a “Goldilocks” tube radius at which capillary wave fluctuations in the tube wall are minimized.

36 MATERIALS SCIENCE↗

Use of spectral analogy to evaluate canopy reflectance sensitivity to leaf optical property

The spectral variation of canopy reflectance is mostly governed by the absorption properties of the elements, hence the leaves, since their intrinsic scattering properties show very little spectral variation. The relationship between canopy reflectance and leaf reflectance measured at the red edge over sugar beet canopies was used to simulate canopy reflectance from leaf reflectance spectra measured over the whole spectral domain. The results show that the spectral analogies found allows accurate reconstruction of canopy reflectance spectra. Explicit assumptions about the very low spectral variation of leaf intrinsic scattering properties are thus indirectly justified. The sensitivity of canopy reflectance (rho(sub c)) to leaf optical properties can then be investigated from concurrent spectral variations of canopy (delta rho(sub c)/delta lambda) and leaf reflectance (delta rho(sub l)/delta lambda): (delta rho(sub c))/(delta rho(sub l)) = ((delta rho(sub c))/(delta lambda) ((delta rho( sub l))/(delta lambda))(sup -1)). This expression is strictly valid only when the optical properties of the soil background or the other vegetation elements such as bark are either spectrally flat or do not contribute significantly to canopy reflectance. Simulations using the SAIL and PROSPECT models demonstrate that the sensitivity of canopy reflectance to leaf reflectance is significant for large vegetation cover fractions in spectral domains where absorption is low. In these conditions, multiple, scattering enhances the leaf absorption features by a factor that can be greater than 2.0. To override the limitations of the SAIL model for the description of the canopy architecture, we tested the previous findings on experimental data. Concurrent canopy and leaf reflectance spectra were measured for a range of sugar beet canopies. The results show good agreement with the theoretical findings. Conclusions are drawn about the applicability of these findings, with particular attention to the potential detectability of leaf biochemical composition from canopy reflectance sensed from space.

Baret, Frederic↗

Visualization of Microbial Biomarkers by Scanning Electron Microscopy

We are developing tools to link the biochemical structure of selected biomarkers with putative biogenic structures observed in mineralized samples. The detection of evidence of life on Mars and other planets will rely on methods that can discriminate compounds formed exclusively by living organisms. While biogenic compounds, such as amino acids and nucleotides have been discovered in extraterrestrial sources, such as meteorites and comets, their formation can be explained by abiotic means. The formation of cellular structures, or more elaborate organic molecules, such as complex lipids, proteins or nucleic acids, however, is strongly correlated to the presence of even the most primitive life processes. Recent evidence lends support to the hypothesis that life may have once existed on Mars. Carbonate globules and ppm concentrations of polycyclic aromatic hydrocarbons (PAHs) have been described in ALH84001, a meteorite originating from Mars ejecta captured by Earth over 13,000 years ago. The localized high concentration of PAHs that follow an increasing gradient from the intact fusion crust towards the interior corresponds to microgram quantities of hydrocarbon. Even though ALH84001 and other similar meteorites have withstood the forces capable of ejecting rock through Mars' escape velocity, upon entering Earth's atmosphere, their core temperatures are likely not to have been raised significantly, as evidenced by the survival of remanent magnetic signatures. Ideal biomarkers of ancient or modern biological life would include molecules that are (or were) pervasive and highly resistant to degradation. Also, requisite methods of detection should be simple, extremely sensitive and broadly inclusive (NASA SP-530). Lipopolysaccharide (LPS), peptidoglycan or pseudopeptidoglycan and beta-glucan are microbial cell wall components which together cover the entire microbial spectrum of eubacteria, archea and fungi. They are all remarkably resistant to thermal degradation. Fortunately, many antimicrobial defense systems of higher organisms require sensitive detection to combat microbial pathogens. We employ here the primitive immune system of the evolutionarily ancient horseshoe crab, Limulus polyphemus. This species relies on multi-enzyme signal amplification detection of cell wall molecules and they can be applied to the development of useful detectors of life. An extension of this work includes the visualization of microbial signatures by labeling LAL components with chromogenic or electron dense markers. The protein Limulus Anti-LPS Factor (LALF) has an extremely high affinity for LPS. By coupling LALF binding with colloidal gold labels we demonstrate a correlation of the structures visible by electron microscopy with biochemical evidence of microbial cell wall materials. Pure silica particles were mixed with cultures of E. coli (10(exp 6) cfu/mL). Samples were washed sequentially with buffered saline, LALF, antibody to LALF and finally colloidal gold-labeled Protein A. Negative controls were not exposed to E. coli but received identical treatment otherwise. Samples were coated with carbon and imaged on a JEOL JSM-840 scanning electron microscope with LaB6 source in the back scatter mode with the JEOL annular back scatter detector. 20 nm-scale black spots in this contrast-reversed image originate from electrons back-scattered by gold atoms. Negative controls did not give any signal. Future work will expand application of this technique to soil simulants and mineralized rock samples.

Wainwright, Norman R.↗

CoarsenConf: Equivariant Coarsening with Aggregated Attention for Molecular Conformer Generation

Molecular conformer generation (MCG) is an important task in cheminformatics and drug discovery. The ability to efficiently generate low-energy 3D structures can avoid expensive quantum mechanical simulations, leading to accelerated virtual screenings and enhanced structural exploration. Several generative models have been developed for MCG, but many struggle to consistently produce high-quality conformers for meaningful downstream applications. To address these issues, we introduce CoarsenConf, which coarse-grains molecular graphs based on torsional angles and integrates them into an SE(3)-equivariant hierarchical variational autoencoder. Through equivariant coarse-graining, we aggregate the fine-grained atomic coordinates of subgraphs connected via rotatable bonds, creating a variable-length coarse-grained latent representation. Our model uses a novel aggregated attention mechanism to restore fine-grained coordinates from the coarse-grained latent representation, enabling efficient generation of accurate conformers. Furthermore, we evaluate the chemical and biochemical quality of our generated conformers on multiple downstream applications, including property prediction and large-scale oracle-based protein docking. Overall, CoarsenConf generates more accurate conformer ensembles compared to prior generative models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗