Predictive chemical kinetics : enabling automatic mechanism generation and evaluation

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Chemical Engineering, 2013.Cataloged from PDF version of thesis.Includes bibliographical references.The use of petroleum-based fuels for transportation accounted for more than 25% of the total energy consumed in 2012, both in the United States and throughout the world. The finite nature of world oil reserves and the effects of burning petroleum-based fuels on the world's climate have motivated efforts to develop alternative, renewable fuels. A major category of alternative fuels is biofuels, which potentially include a wide variety of hydrocarbons, alcohols, aldehydes, ketones, ethers, esters, etc. To select the best species for use as fuel, we need to know if it burns cleanly, controllably, and efficiently. This is especially important when considering novel engine technologies, which are often very sensitive to fuel chemistry. The large number of candidate fuels and the high expense of experimental engine tests motivates the use of predictive theoretical methods to help quickly identify the most promising candidates. This thesis presents several contributions in the areas of predictive chemical kinetics and automatic mechanism generation, particularly in the area of reaction kinetics. First, the accuracy of several methods of automatic, high-throughput estimation of reaction rates are evaluated by comparison to a test set obtained from the NIST Chemical Kinetics Database. The methods considered, including the classic Evans-Polanyi correlation, the "rate rules" method currently used in the RMG software, and a new method based on group contribution theory, are shown to not yet obtain the order-of-magnitude accuracy desired for automatic mechanism generation. Second, a method of very accurate computation of bimolecular reaction rates using ring polymer molecular dynamics (RPMD) is presented. RPMD rate theory enables the incorporation of quantum effects (zero-point energy and tunneling) in reaction kinetics using classical molecular dynamics trajectories in an extended phase space. A general-purpose software package named RPMD-rate was developed for conducting such calculations, and the accuracy of this method was demonstrated by investigating the kinetics and kinetic isotope effect of the reaction OH + CH4 --> CH3 + H2O. Third, a general framework for incorporating pressure dependence in thermal unimolecular reactions, which require an inert third body to provide or remove the energy needed for reaction via bimolecular collisions, was developed. Within this framework, several methods of reducing the full, master equation-based model to a set of phenomenological rate coefficients k(T, P) are compared using the chemically-activated reaction of acetyl radical with oxygen as a case study, and recommendations are made as to when each method should be used. This also resulted in a general-purpose code for calculating pressure-dependent kinetics, which was applied to developing an ab initio model of the reaction of the Criegee biradical CH 200 with small carbonyls that reproduces recent experimental results. Finally, the ideas and techniques of estimating reaction kinetics are brought together for the development of a detailed kinetics model of the oxidation of diisopropyl ketone (DIPK), a candidate biofuel representative of species produced from cellulosic biomass conversion using endophytic fungi. The model is evaluated against three experiments covering a range of temperatures, pressures, and oxygen concentrations to show its strengths and weaknesses. Our ability to automatically generate this model and systematically improve its parameters without fitting to the experimental results demonstrates the validity and usefulness of the predictive chemical kinetics paradigm. These contributions are available as part of the Reaction Mechanism Generator (RMG) software package.by Joshua W. Allen.Ph.D

RPMDrate: Bimolecular chemical reaction rates from ring polymer molecular dynamics

We present RPMDrate, a computer program for the calculation of gas phase bimolecular reaction rate coefficients using the ring polymer molecular dynamics (RPMD) method. The RPMD rate coefficient is calculated using the Bennett–Chandler method as a product of a static (centroid density quantum transition state theory (QTST) rate) and a dynamic (ring polymer transmission coefficient) factor. The computational procedure is general and can be used to treat bimolecular polyatomic reactions of any complexity in their full dimensionality. The program has been tested for the H+H2, H+CH4, OH+CH4 and H+C2H6 reactions.United States. Dept. of Energy (Office of Basic Energy Sciences under the Energy Frontier Research Center for Combustion Science (Grant No. DE-SC0001198))United States. Dept. of Energy (Energy Frontier Research Center for Combustion Science, Combustion Energy Research Fellowship)King Abdullah University of Science and Technology (Award No. KUS-I1-010-01

Assessing mesoscale material response under shock & isentropic compression via high-resolution line-imaging VISAR.

Of special promise for providing dynamic mesoscale response data is the line-imaging VISAR, an instrument for providing spatially resolved velocity histories in dynamic experiments. We have prepared two line-imaging VISAR systems capable of spatial resolution in the 10-20 micron range, at the Z and STAR facilities. We have applied this instrument to selected experiments on a compressed gas gun, chosen to provide initial data for several problems of interest, including: (1) pore-collapse in copper (two variations: 70 micron diameter hole in single-crystal copper) and (2) response of a welded joint in dissimilar materials (Ta, Nb) to ramp loading relative to that of a compression joint. The instrument is capable of resolving details such as the volume and collapse history of a collapsing isolated pore

Missing OH reactivity in the global marine boundary layer

The hydroxyl radical (OH) reacts with thousands of chemical species in the atmosphere, initiating their removal and the chemical reaction sequences that produce ozone, secondary aerosols, and gas-phase acids. OH reactivity, which is the inverse of OH lifetime, influences the OH abundance and the ability of OH to cleanse the atmosphere. The NASA Atmospheric Tomography (ATom) campaign used instruments on the NASA DC-8 aircraft to measure OH reactivity and more than 100 trace chemical species. ATom presented a unique opportunity to test the completeness of the OH reactivity calculated from the chemical species measurements by comparing it to the measured OH reactivity over two oceans across four seasons. Although the calculated OH reactivity was below the limit of detection for the ATom instrument used to measure OH reactivity throughout much of the free troposphere, the instrument was able to measure the OH reactivity in and just above the marine boundary layer. The mean measured value of OH reactivity in the marine boundary layer across all latitudes and all ATom deployments was 1.9 s⁻¹, which is 0.5 s⁻¹ larger than the mean calculated OH reactivity. The missing OH reactivity, the difference between the measured and calculated OH reactivity, varied between 0 and 3.5 s⁻¹, with the highest values over the Northern Hemisphere Pacific Ocean. Correlations of missing OH reactivity with formaldehyde, dimethyl sulfide, butanal, and sea surface temperature suggest the presence of unmeasured or unknown volatile organic compounds or oxygenated volatile organic compounds associated with ocean emissions

Eclipses During the 2010 Eruption of the Recurrent Nova U Scorpii

The eruption of the recurrent nova U Scorpii on 28 January 2010 is now the all-time best observed nova event. We report 36,776 magnitudes throughout its 67 day eruption, for an average of one measure every 2.6 minutes. This unique and unprecedented coverage is the first time that a nova has any substantial amount of fast photometry. With this, two new phenomena have been discovered: the fast flares in the early light curve seen from days 9-15 (which have no proposed explanation) and the optical dips seen out of eclipse from days 41-61 (likely caused by raised rims of the accretion disk occulting the bright inner regions of the disk as seen over specific orbital phases). The expanding shell and wind cleared enough from days 12-15 so that the inner binary system became visible, resulting in the sudden onset of eclipses and the turn-on of the supersoft X-ray source. On day 15, a strong asymmetry in the out-of-eclipse light points to the existence of the accretion stream. The normal optical flickering restarts on day 24.5. For days 15-26, eclipse mapping shows that the optical source is spherically symmetric with a radius of 4.1 R_sun. For days 26-41, the optical light is coming from a rim-bright disk of radius 3.4 R_sun. For days 41-67, the optical source is a center-bright disk of radius 2.2 R_sun. Throughout the eruption, the colors remain essentially constant. We present 12 eclipse times during eruption plus five just after the eruption.Comment: ApJ in press. 60 pages, 17 figure

Whole exome sequencing of circulating tumor cells provides a window into metastatic prostate cancer

Comprehensive analyses of cancer genomes promise to inform prognoses and precise cancer treatments. A major barrier, however, is inaccessibility of metastatic tissue. A potential solution is to characterize circulating tumor cells (CTCs), but this requires overcoming the challenges of isolating rare cells and sequencing low-input material. Here we report an integrated process to isolate, qualify and sequence whole exomes of CTCs with high fidelity, using a census-based sequencing strategy. Power calculations suggest that mapping of >99.995% of the standard exome is possible in CTCs. We validated our process in two prostate cancer patients including one for whom we sequenced CTCs, a lymph node metastasis and nine cores of the primary tumor. Fifty-one of 73 CTC mutations (70%) were observed in matched tissue. Moreover, we identified 10 early-trunk and 56 metastatic-trunk mutations in the non-CTC tumor samples and found 90% and 73% of these, respectively, in CTC exomes. This study establishes a foundation for CTC genomics in the clinic