Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reaction”

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 685 records · Page 38

The Rate Constant for the Reaction H + C2H5 at T = 295 - 150K

The reaction between the hydrogen atom and the ethyl (C2H3) radical is predicted by photochemical modeling to be the most important loss process for C2H5 radicals in the atmospheres of Jupiter and Saturn. This reaction is also one of the major sources for the methyl radicals in these atmospheres. These two simplest hydrocarbon radicals are the initial species for the synthesis of larger hydrocarbons. Previous measurements of the rate constant for the H + C2H5 reaction varied by a factor of five at room temperature, and some studies showed a dependence upon temperature while others showed no such dependence. In addition, the previous studies were at higher temperatures and generally higher pressures than that needed for use in planetary atmospheric models. The rate constant for the reaction H + C2H5 has been measured directly at T = 150, 202 and 295 K and at P = 1.0 Torr He for all temperatures and additionally at P = 0.5 and 2.0 Torr He at T = 202 K. The measurements were performed in a discharge - fast flow system. The decay of the C2H5 radical in the presence of excess hydrogen was monitored by low-energy electron impact mass spectrometry under pseudo-first order conditions. H atoms and C2H5 radicals were generated rapidly and simultaneously by the reaction of fluorine atoms with H2 and C2H6, respectively. The total rate constant was found to be temperature and pressure independent. The measured total rate constant at each temperature are: k(sub 1)(295K) = (1.02+/-0.24)x10(exp -10), k(sub 1)(202K) = (1.02+/-0.22)x10(exp -10) and k(sub 1)(150K) = (0.93+/-0.21)x10(exp -10), all in units of cu cm/molecule/s. The total rate constant derived from all the combined measurements is k(sub 1) = (l.03+/-0.17)x10(exp -10) cu cm/molecule/s. At room temperature our results are about a factor of two higher than the recommended rate constant and a factor of three lower than the most recently published study.

Pimentel, Andre S.↗

The rate of the reaction between C2H and C2H2 at interstellar temperatures

The reaction between the radical C2H and the stable hydrocarbon C2H2 is one of the simplest neutral-neutral hydrocarbon reactions in chemical models of dense interstellar clouds and carbon-rich circumstellar shells. Although known to be rapid at temperatures > or = 300 K, the reaction has yet to be studied at lower temperatures. We present here ab initio calculations of the potential surface for this reaction and dynamical calculations to determine its rate at low temperature. Despite a small potential barrier in the exit channel, the calculated rate is large, showing that this reaction and, most probably, more complex analogs contribute to the formation of complex organic molecules in low-temperature sources.

Non-NASA Center↗

Air oxidation of hydrazine. 1. Reaction kinetics on natural kaolinites, halloysites, and model substituent layers with varying iron and titanium oxide and O- center contents

Air oxidation of hydrazine was studied by using a group of kaolinites, halloysites, and substituent oxides as models for the tetrahedral and octahedral sheets. The rate was found to be linear with oxygen. The stoichiometry showed that oxygen was the primary oxidant and that dinitrogen was the only important nitrogen-containing product. The rates on kaolinites were strongly inhibited by water. Those on three-dimensional silica and gibbsite appeared not to be. That on a supposedly layered silica formed from a natural kaolinite by acid leaching showed transitional behavior--slowed relative to that expected from a second-order reaction relative to that on the gibbsite and silica but faster than those on the kaolinites. The most striking result of the reaction was the marked increase in the rate of reaction of a constant amount of hydrazine as the amount of clay was increased. The increase was apparent (in spite of the water inhibition at high conversions) over a 2 order of magnitude variation of the clay weight. The weight dependence was taken to indicate that the role of the clay is very important, that the number of reactive centers is very small, or that they may be deactivated over the course of the reaction. In contrast to the strong dependence on overall amount of clay, the variation of amounts of putative oxidizing centers, such as structural Fe(III), admixed TiO2 or Fe2O3, or O- centers, did not result in alteration of the rate commensurate with the degree of variation of the entity in question. Surface iron does play some role, however, as samples that were pretreated with a reducing agent were less active as catalysts than the parent material. These results were taken to indicate either that the various centers interact to such a degree that they cannot be considered independently or that the reaction might proceed by way of surface complexation, rather than single electron transfers.

NASA Discipline Number 52-20↗

Role of microbes in the smectite-to-illite reaction

Temperature, pressure, and time have been thought to control the smectite-to-illite (S-I) reaction, an important diagenetic process used for petroleum exploration. We demonstrated that microorganisms can promote the S-I reaction by dissolving smectite through reduction of structural Fe(III) at room temperature and 1 atmosphere within 14 days. This reaction typically requires conditions of 300 degrees to 350 degrees C, 100 megapascals, and 4 to 5 months in the absence of microbial activity. These results challenge the conventional concept of the S-I reaction and of reaction kinetic models.

Minerals/chemistry/metabolism↗

Minimum-Time and Vibration Avoidance Attitude Maneuver for Spacecraft with Torque and Momentum Limit Constraints in Redundant Reaction Wheel Configuration

In this paper, we present an optimal open-loop slew trajectory algorithm developed at GSFC for the so-called "Yardstick design" of the James Webb Space Telescope (JWST). JWST is an orbiting infrared observatory featuring a lightweight, segmented primary mirror approximately 6 meters in diameter and a sunshield approximately the size of a tennis court. This large, flexible structure will have significant number of lightly damped, dominant flexible modes. With very stringent requirements on pointing accuracy and image quality, it is important that slewing be done within the required time constraint and with minimal induced vibration in order to maximize observing efficiency. With reaction wheels as control actuators, initial wheel speeds as well as individual wheel torque and momentum limits become dominant constraints in slew performance. These constraints must be taken into account when performing slews to ensure that unexpected reaction wheel saturation does not occur, since such saturation leads to control failure in accurately tracking commanded motion and produces high frequency torque components capable of exciting structural modes. A minimum-time constraint is also included and coupled with reaction wheel limit constraints in the optimization to minimize both the effect of the control torque on the flexible body motion and the maneuver time. The optimization is on slew command parameters, such as maximum slew velocity and acceleration, for a given redundant reaction wheel configuration and is based on the dynamic interaction between the spacecraft and reaction wheel motion. Analytical development of the slew algorithm to generate desired slew position, rate, and acceleration profiles to command a feedback/feed forward control system is described. High-fidelity simulation and experimental results are presented to show that the developed slew law achieves the objectives.

Ha, Kong Q.↗

Researches on Preliminary Chemical Reactions in Spark-Ignition Engines

Chemical reactions can demonstrably occur in a fuel-air mixture compressed in the working cylinder of an Otto-cycle (spark ignition) internal-combustion engine even before the charge is ignited by the flame proceeding from the sparking plug. These are the so-called "prelinminary reactions" ("pre-flame" combustion or oxidation), and an exact knowledge of their characteristic development is of great importance for a correct appreciation of the phenomena of engine-knock (detonation), and consequently for its avoidance. Such reactions can be studied either in a working engine cylinder or in a combustion bomb. The first method necessitates a complicated experimental technique, while the second has the disadvantage of enabling only a single reaction to be studied at one time. Consequently, a new series of experiments was inaugurated, conducted in a motored (externally-driven) experimental engine of mixture-compression type, without ignition, the resulting preliminary reactions being detectable and measurable thermometrically.

Muehlner, E.↗

Blow-Out Velocities of Solutions of Hydrocarbons and Boron Hydride - Hydrocarbon Reaction Products in a 1 7/8-Inch-Diameter Combustor

Blow-out velocities were determined for JP-4 solutions containing: (1) 10 % ethylene - decaborane reaction product, (2) 10% and 20% acetylene - diborane reaction product, and (3) 5.5%, 15.7%, and 30.7% methylacetylene - diborane reaction product. These were compared with blow-out velocities for JP-4, propylene oxide, and neohexane and previously reported data for JP-4 solutions of pentaborane. For those reaction products investigated, the blow-out velocities at a fixed equivalence ratio were higher for those materials containing higher boron concentrations; that is, blow-out velocity increased in the following order: (1) methylacetylene - diborane, (2) acetylene - diborane, and (3) ethylene - decaborane reaction products.

Morris, James F.↗

Modeling Chemical Growth Processes in Titan's Atmosphere: 1. Theoretical Rates for Reactions between Benzene and the Ethynyl (C2H) and Cyano (CN) Radicals at Low Temperature and Pressure

Density functional theory calculations at the B3LYP/6-31+G** level were employed to characterize the critical points for adducts, isomers, products, and intervening transition states for the reactions between benzene and the ethynyl (C2H) or cyano (CN) radicals. Both addition reactions were found to have no barriers in their entrance channels, making them efficient at the low temperature and pressure conditions that prevail in the haze-forming region of Titan's atmosphere as well as in the dense interstellar medium (ISM). The dominant products are ethynylbenzene (C6H5C2H) and cyanobenzene (C6H5CN). Hydrogen abstraction reactions were also characterized but found to be non-competitive. Trajectory calculations based on potentials fit to about 600 points calculated at the ROMP2/6-31+G** level for each interaction surface were used to determine reaction rates. The rates incorporated any necessary corrections for back reactions as ascertained from a multiwell treatment used to determine outcome distributions over the range of temperatures and pressures pertinent to Titan and the ISM and are in good agreement with the limited available experimental data.

Woon, David E.↗

NASA Ares I Launch Vehicle Roll and Reaction Control Systems Overview

NASA's Ares I launch vehicle, consisting of a five segment solid rocket booster first stage and a liquid bi-propellant J-2X engine upper stage, is the vehicle that's been chosen to return humans to the moon, mars, and beyond. This paper provides an overview of the work that has taken place on the Ares I launch vehicle roll and reaction control systems. Reaction control systems are found on many launch vehicles and provide a vehicle with a three degree of freedom stabilization during the mission. The Ares I baseline configuration currently consists of a first stage roll control system that will provide the vehicle with a method of counteracting the roll torque that is expected during launch. An upper stage reaction control system will allow the upper stage three degrees of freedom control as needed. Design assessments and trade studies are being conducted on the roll and reaction control systems including: propellant selection, thruster arrangement, pressurization system configuration, and system component trades. Other vehicle considerations and issues include thruster plume impingement, thruster module aerothermal and aerodynamic effects, and system integration. This paper concludes by summarizing the process of down selecting to the current baseline configuration for the Ares I roll and reaction control systems.

Popp, Chris↗

Reaction Wheel Disturbance Model Extraction Software - RWDMES

The RWDMES is a tool for modeling the disturbances imparted on spacecraft by spinning reaction wheels. Reaction wheels are usually the largest disturbance source on a precision pointing spacecraft, and can be the dominating source of pointing error. Accurate knowledge of the disturbance environment is critical to accurate prediction of the pointing performance. In the past, it has been difficult to extract an accurate wheel disturbance model since the forcing mechanisms are difficult to model physically, and the forcing amplitudes are filtered by the dynamics of the reaction wheel. RWDMES captures the wheel-induced disturbances using a hybrid physical/empirical model that is extracted directly from measured forcing data. The empirical models capture the tonal forces that occur at harmonics of the spin rate, and the broadband forces that arise from random effects. The empirical forcing functions are filtered by a physical model of the wheel structure that includes spin-rate-dependent moments (gyroscopic terms). The resulting hybrid model creates a highly accurate prediction of wheel-induced forces. It accounts for variation in disturbance frequency, as well as the shifts in structural amplification by the whirl modes, as the spin rate changes. This software provides a point-and-click environment for producing accurate models with minimal user effort. Where conventional approaches may take weeks to produce a model of variable quality, RWDMES can create a demonstrably high accuracy model in two hours. The software consists of a graphical user interface (GUI) that enables the user to specify all analysis parameters, to evaluate analysis results and to iteratively refine the model. Underlying algorithms automatically extract disturbance harmonics, initialize and tune harmonic models, and initialize and tune broadband noise models. The component steps are described in the RWDMES user s guide and include: converting time domain data to waterfall PSDs (power spectral densities); converting PSDs to order analysis data; extracting harmonics; initializing and simultaneously tuning a harmonic model and a wheel structural model; initializing and tuning a broadband model; and verifying the harmonic/broadband/structural model against the measurement data. Functional operation is through a MATLAB GUI that loads test data, performs the various analyses, plots evaluation data for assessment and refinement of analysis parameters, and exports the data to documentation or downstream analysis code. The harmonic models are defined as specified functions of frequency, typically speed-squared. The reaction wheel structural model is realized as mass, damping, and stiffness matrices (typically from a finite element analysis package) with the addition of a gyroscopic forcing matrix. The broadband noise model is realized as a set of speed-dependent filters. The tuning of the combined model is performed using nonlinear least squares techniques. RWDMES is implemented as a MATLAB toolbox comprising the Fit Manager for performing the model extraction, Data Manager for managing input data and output models, the Gyro Manager for modifying wheel structural models, and the Harmonic Editor for evaluating and tuning harmonic models. This software was validated using data from Goodrich E wheels, and from GSFC Lunar Reconnaissance Orbiter (LRO) wheels. The validation testing proved that RWDMES has the capability to extract accurate disturbance models from flight reaction wheels with minimal user effort.

Blaurock, Carl↗

NASA Ares I Launch Vehicle Roll and Reaction Control Systems Design Status

This paper provides an update of design status following the preliminary design review of NASA s Ares I first stage roll and upper stage reaction control systems. The Ares I launch vehicle has been chosen to return humans to the moon, mars, and beyond. It consists of a first stage five segment solid rocket booster and an upper stage liquid bi-propellant J-2X engine. Similar to many launch vehicles, the Ares I has reaction control systems used to provide the vehicle with three degrees of freedom stabilization during the mission. During launch, the first stage roll control system will provide the Ares I with the ability to counteract induced roll torque. After first stage booster separation, the upper stage reaction control system will provide the upper stage element with three degrees of freedom control as needed. Trade studies and design assessments conducted on the roll and reaction control systems include: propellant selection, thruster arrangement, pressurization system configuration, and system component trades. Since successful completion of the preliminary design review, work has progressed towards the critical design review with accomplishments made in the following areas: pressurant / propellant tank, thruster assembly, and other component configurations, as well as thruster module design, and waterhammer mitigation approach. Also, results from early development testing are discussed along with plans for upcoming system testing. This paper concludes by summarizing the process of down selecting to the current baseline configuration for the Ares I roll and reaction control systems.

Butt, Adam↗

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗

Determination of the Rate Coefficients of the SO2 plus O plus M yields SO3 plus M Reaction

Rate coefficients of the title reaction R(sub 31) (SO2 +O+M yields SO3 +M) and R(sub 56) (SO2 + HO2 yields SO3 +OH), important in the conversion of S(IV) to S(VI),were obtained at T =970-1150 K and rho (sub ave) = 16.2 micro mol/cubic cm behind reflected shock waves by a perturbation method. Shock-heated H2/ O2/Ar mixtures were perturbed by adding small amounts of SO2 (1%, 2%, and 3%) and the OH temporal profiles were then measured using laser absorption spectroscopy. Reaction rate coefficients were elucidated by matching the characteristic reaction times acquired from the individual experimental absorption profiles via simultaneous optimization of k(sub 31) and k(sub 56) values in the reaction modeling (for satisfactory matches to the observed characteristic times, it was necessary to take into account R(sub 56)). In the experimental conditions of this study, R(sub 31) is in the low-pressure limit. The rate coefficient expressions fitted using the combined data of this study and the previous experimental results are k(sub 31,0)/[Ar] = 2.9 10(exp 35) T(exp ?6.0) exp(?4780 K/T ) + 6.1 10(exp 24) T(exp ?3.0) exp(?1980 K/T ) cm(sup 6) mol(exp ?2)/ s at T = 300-2500 K; k(sub 56) = 1.36 10(exp 11) exp(?3420 K/T ) cm(exp 3)/mol/s at T = 970-1150 K. Computer simulations of typical aircraft engine environments, using the reaction mechanism with the above k(sub 31,0) and k(sub 56) expressions, gave the maximum S(IV) to S(VI) conversion yield of ca. 3.5% and 2.5% for the constant density and constant pressure flow condition, respectively. Moreover, maximum conversions occur at rather higher temperatures (?1200 K) than that where the maximum k(sub 31,0) value is located (approximately 800 K). This is because the conversion yield is dependent upon not only the k(sup 31,0) and k(sup 56) values (production flux) but also the availability of H, O, and HO2 in the system (consumption flux).

Hwang, S. M.↗

Step-by-Step Simulation of Radiation Chemistry Using Green Functions for Diffusion-Influenced Reactions

Radiolytic species are formed approximately 1 ps after the passage of ionizing radiation through matter. After their formation, they diffuse and chemically react with other radiolytic species and neighboring biological molecules, leading to various oxidative damage. Therefore, the simulation of radiation chemistry is of considerable importance to understand how radiolytic species damage biological molecules [1]. The step-by-step simulation of chemical reactions is difficult, because the radiolytic species are distributed non-homogeneously in the medium. Consequently, computational approaches based on Green functions for diffusion-influenced reactions should be used [2]. Recently, Green functions for more complex type of reactions have been published [3-4]. We have developed exact random variate generators of these Green functions [5], which will allow us to use them in radiation chemistry codes. Moreover, simulating chemistry using the Green functions is which is computationally very demanding, because the probabilities of reactions between each pair of particles should be evaluated at each timestep [2]. This kind of problem is well adapted for General Purpose Graphic Processing Units (GPGPU), which can handle a large number of similar calculations simultaneously. These new developments will allow us to include more complex reactions in chemistry codes, and to improve the calculation time. This code should be of importance to link radiation track structure simulations and DNA damage models.

Plante, Ianik↗

A Theoretical Investigation of the Plausibility of Reactions Between Ammonia and Carbonyl Species (Formaldehyde, Acetaldehyde, and Acetone) in Interstellar Ice Analogs at Ultracold Temperatures

We have reexamined the reaction between formaldehyde and ammonia, which was previously studied by us and other workers in modestly sized cluster calculations. Larger model systems with up to 12H2O were employed, and reactions of two more carbonyl species, acetaldehyde and acetone, were also carried out. Calculations were performed at the B3LYP/6-31+G** level with bulk solvent effects treated with a polarizable continuum model; limited MP2/6-31+G** calculations were also performed. We found that while the barrier for the concerted proton relay mechanism described in previous work remains modest, it is still prohibitively high for the reaction to occur under the ultracold conditions that prevail in dense interstellar clouds. However, a new pathway emerged in more realistic clusters that involves at least one barrierless step for two of the carbonyl species considered here: ammonia reacts with formaldehyde and acetaldehyde to form a partial charge transfer species in small clusters (4H2O) and a protonated hydroxyamino intermediate species in large clusters (9H2O, 12H2O); modest barriers that decrease sharply with cluster size are found for the analogous processes for the acetone-NH3 reaction. Furthermore, if a second ammonia replaces one of the water molecules in calculations in the 9H2O clusters, deprotonation can occur to yield the same neutral hydroxyamino species that is formed via the original concerted proton relay mechanism. In at least one position, deprotonation is barrierless when zero-point energy is included. In addition to describing the structures and energetics of the reactions between formaldehyde, acetaldehyde, and acetone with ammonia, we report spectroscopic predictions of the observable vibrational features that are expected to be present in ice mixtures of different composition.

Chen, Lina↗

Controlled Nitric Oxide Production via O(1D) + N2O Reactions for Use in Oxidation Flow Reactor Studies

Oxidation flow reactors that use low-pressure mercury lamps to produce hydroxyl (OH) radicals are an emerging technique for studying the oxidative aging of organic aerosols. Here, ozone (O3) is photolyzed at 254 nm to produce O(1D) radicals, which react with water vapor to produce OH. However, the need to use parts-per-million levels of O3 hinders the ability of oxidation flow reactors to simulate NOx-dependent secondary organic aerosol (SOA) formation pathways. Simple addition of nitric oxide (NO) results in fast conversion of NOx (NO+NO2) to nitric acid (HNO3), making it impossible to sustain NOx at levels that are sufficient to compete with hydroperoxy (HO2) radicals as a sink for organic peroxy (RO2) radicals. We developed a new method that is well suited to the characterization of NOx-dependent SOA formation pathways in oxidation flow reactors. NO and NO2 are produced via the reaction O(1D)+N2O->2NO, followed by the reaction NO+O3->NO2+O2. Laboratory measurements coupled with photochemical model simulations suggest that O(1D)+N2O reactions can be used to systematically vary the relative branching ratio of RO2 +NO reactions relative to RO2 +HO2 and/or RO2+RO2 reactions over a range of conditions relevant to atmospheric SOA formation. We demonstrate proof of concept using high-resolution time-of-flight chemical ionization mass spectrometer (HR-ToF-CIMS) measurements with nitrate (NO-3 ) reagent ion to detect gas-phase oxidation products of isoprene and -pinene previously observed in NOx-influenced environments and in laboratory chamber experiments.

Lambe, Andrew↗

Low Friction Reaction Wheel System and Containment Structure with Integrated Electromagnets

The disclosed subject matter relates to a reaction wheel system for controlling and stabilizing a satellite or other spacecraft. The reaction wheel system includes a reaction wheel rotor and a containment structure stator including a plurality of contactor assemblies for securing the reaction wheel rotor when not in rotation and for electromagnetically inducing substantially frictionless rotation of the reaction wheel rotor.

Yew, Alvin G.↗

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗