Carbon Core Electron Spectra of Polycyclic Aromatic Hydrocarbons.

Aromaticity profoundly affects molecular orbitals in polycyclic aromatic hydrocarbons. X-ray core electron spectroscopy has observed that carbon 1s-π* transitions can be broadened or even split in some polycyclic systems, although the origin of the effect has remained obscure. The π electrons in polycyclic systems are typically classified in the Clar model as belonging to either true aromatic sextets (similar to benzene) or isolated double bonds (similar to olefins). Here, bulk-sensitive carbon core excitation spectra are presented for a series of polycyclic systems and show that the magnitude of the 1s-π* splitting is determined primarily by the ratio of true aromatic sextets to isolated double bonds. The observed splitting can be rationalized in terms of ground state energetics as described by Hückel, driven by the π electron structure described by Clar. This simple model including only ground state energetics is shown to explain the basics physics behind the spectral evolution for a broad set of polycyclic aromatic hydrocarbons, although some residual deviations between this model and experiment can likely be improved by including a more detailed electronic structure and the core hole effect.

Realistic molecular model of kerogen's nanostructure.

Despite kerogen's importance as the organic backbone for hydrocarbon production from source rocks such as gas shale, the interplay between kerogen's chemistry, morphology and mechanics remains unexplored. As the environmental impact of shale gas rises, identifying functional relations between its geochemical, transport, elastic and fracture properties from realistic molecular models of kerogens becomes all the more important. Here, by using a hybrid experimental-simulation method, we propose a panel of realistic molecular models of mature and immature kerogens that provide a detailed picture of kerogen's nanostructure without considering the presence of clays and other minerals in shales. We probe the models' strengths and limitations, and show that they predict essential features amenable to experimental validation, including pore distribution, vibrational density of states and stiffness. We also show that kerogen's maturation, which manifests itself as an increase in the sp(2)/sp(3) hybridization ratio, entails a crossover from plastic-to-brittle rupture mechanisms.

Two-step laser mass spectrometry of asphaltenes.

Defined by their solubility in toluene and insolubility in n-heptane, asphaltenes are a highly aromatic, polydisperse mixture consisting of the heaviest and most polar fraction of crude oil. Although asphaltenes are critically important to the exploitation of conventional oil and are poised to rise in significance along with the exploitation of heavy oil, even as fundamental a quantity as their molecular weight distribution is unknown to within an order of magnitude. Laser desorption/ionization (LDI) mass spectra vary greatly with experimental parameters so are difficult to interpret: some groups favor high laser pulse energy measurements (yielding heavy molecular weights), arguing that high pulse energy is required to detect the heaviest components of this mixture; other groups favor low pulse energy measurements (yielding light molecular weights), arguing that low pulse energy is required to avoid aggregation in the plasma plume. Here we report asphaltene mass spectra recorded with two-step laser mass spectrometry (L2MS), in which desorption and ionization are decoupled and no plasma is produced. L2MS mass spectra of asphaltenes are insensitive to laser pulse energy and other parameters, demonstrating that the asphaltene molecular weight distribution can be measured without limitation from insufficient laser pulse energy or plasma-phase aggregation. These data resolve the controversy from LDI, showing that the asphaltene molecular weight distribution peaks near 600 Da and previous measurements reporting much heavier species suffered from aggregation effects.

Nanoscale geochemical and geomechanical characterization of organic matter in shale.

Solid organic matter (OM) plays an essential role in the generation, migration, storage, and production of hydrocarbons from economically important shale rock formations. Electron microscopy images have documented spatial heterogeneity in the porosity of OM at nanoscale, and bulk spectroscopy measurements have documented large variation in the chemical composition of OM during petroleum generation. However, information regarding the heterogeneity of OM chemical composition at the nanoscale has been lacking. Here we demonstrate the first application of atomic force microscopy-based infrared spectroscopy (AFM-IR) to measure the chemical and mechanical heterogeneity of OM in shale at the nanoscale, orders of magnitude finer than achievable by traditional chemical imaging tools such as infrared microscopy. We present a combination of optical microscopy and AFM-IR imaging to characterize OM heterogeneity in an artificially matured series of New Albany Shales. The results document the evolution of individual organic macerals with maturation, providing a microscopic picture of the heterogeneous process of petroleum generation.

Regioselectivity control in a ruthenium-catalyzed cycloisomerization of diyne-ols.

The ruthenium-catalyzed cycloisomerization of diynes containing one silyl alkyne and one propargyl alcohol yields 2-silyl-[6H]-pyrans instead of the expected unsaturated acylsilanes except when additional conjugation of a aromatic ring is present at the delta-position. Under certain conditions, a facile ruthenium-catalyzed isomerization of the product takes place as well. This regioselectivity of the cyclization can be controlled by the choice of solvent system. DFT calculations confirm the expected greater stability of the silyl-pyrans relative to the acylsilanes.

Downhole fluid analysis and asphaltene science for petroleum reservoir evaluation.

Petroleum reservoirs are enshrouded in mysteries associated with all manner of geologic and fluid complexities that Mother Nature can inspire. Efficient exploitation of petroleum reservoirs mandates elucidation of these complexities; downhole fluid analysis (DFA) has proven to be indispensable for understanding both fluids and reservoir architecture. Crude oil consists of dissolved gases, liquids, and dissolved solids, known as the asphaltenes. These different fluid components exhibit fluid gradients vertically and laterally, which are best revealed by DFA, with its excellent precision and accuracy. Compositional gradient analysis falls within the purview of thermodynamics. Gas-liquid equilibria can be treated with a cubic equation of state (EoS), such as the Peng-Robinson EoS, a modified van der Waals EoS. In contrast, the first EoS for asphaltene gradients, the Flory-Huggins-Zuo (FHZ) EoS, was developed only recently. The resolution of the asphaltene molecular and nanocolloidal species in crude oil, which is codified in the Yen-Mullins model of asphaltenes, enabled the development of this EoS. The combination of DFA characterization of gradients of reservoir crude oil with the cubic EoS and FHZ EoS analyses brings into view wide-ranging reservoir concerns, such as reservoir connectivity, fault-block migration, heavy oil gradients, tar mat formation, huge disequilibrium fluid gradients, and even stochastic variations of reservoir fluids. New petroleum science and DFA technology are helping to offset the increasing costs and technical difficulties of exploiting ever-more-remote petroleum reservoirs.

Minimization of fragmentation and aggregation by laser desorption laser ionization mass spectrometry.

Measuring average quantities in complex mixtures can be challenging for mass spectrometry, as it requires ionization and detection with nearly equivalent cross-section for all components, minimal matrix effect, and suppressed signal from fragments and aggregates. Fragments and aggregates are particularly troublesome for complex mixtures, where they can be incorrectly assigned as parent ions. Here we study fragmentation and aggregation in six aromatic model compounds as well as petroleum asphaltenes (a naturally occurring complex mixture) using two laser-based ionization techniques: surface assisted laser desorption ionization (SALDI), in which a single laser desorbs and ionizes solid analytes; and laser ionization laser desorption mass spectrometry (L(2)MS), in which desorption and ionization are separated spatially and temporally with independent lasers. Model compounds studied include molecules commonly used as matrices in single laser ionization techniques such as matrix assisted laser desorption ionization (MALDI). We find significant fragmentation and aggregation in SALDI, such that individual fragment and aggregate peaks are typically more intense than the parent peak. These fragment and aggregate peaks are expected in MALDI experiments employing these compounds as matrices. On the other hand, we observe no aggregation and only minimal fragmentation in L(2)MS. These results highlight some advantages of L(2)MS for analysis of complex mixtures such as asphaltenes.

Analysis of petroleum compositional similarity using multiway principal components analysis (MPCA) with comprehensive two-dimensional gas chromatographic data.

The accurate establishment of oil similarity is a longstanding problem in petroleum geochemistry and a necessary component for resolving the architecture of an oil reservoir. Past limitations have included the excessive reliance on a relatively small number of biomarkers to characterize such complex fluids as crude oils. Here we use multiway principal components analysis (MPCA) on large numbers of specific chemical components resolved with comprehensive two-dimensional gas chromatography-flame ionization detection (GC×GC-FID) to determine the molecular relatedness of eight different maltene fractions of crude oils. MPCA works such that every compound eluting within the same first and second dimension retention time is quantitatively compared with what elutes at that same retention times within the other maltene fractions. Each maltene fraction and corresponding MPCA analysis contains upwards of 3500 quantified components. Reservoir analysis included crude oil sample pairs from around the world that were collected sequentially at depth within a single well, collected from multiple depths in the same well, and from different depths and different wells but thought to be intersected by the same permeable strata. Furthermore, three different regions of each GC×GC-FID chromatograms were analysed to evaluate the effectiveness of MPCA to resolve compositional changes related to the source of the oil generating sediments and its exposure to biological and/or physical weathering processes. Compositional and instrumental artefacts introduced during sampling and processing were also quantitatively evaluated. We demonstrate that MPCA can resolve multi-molecular differences between oil samples as well as provide insight into the overall molecular relatedness between various crude oils.

Inverting MRI measurements to heterogeneity spectra.

Spatially resolved MRI measurements of porosity and relaxation time have been performed on a series of sandstone and carbonate rock cores in order to assess spatial heterogeneity in these samples. Geostatistical techniques such as the construction of experimental variograms provide a quantitative measure of heterogeneity, although the interpretation of standard techniques is at times ambiguous. Here, we attempt to resolve some of that ambiguity by addressing the influence of regularization (spatial averaging over the volume of a voxel) on the variogram. Modeling the influence of regularization allows measured variograms to be inverted, yielding a heterogeneity spectrum that shows the extent of spatial heterogeneity as a function of length scale. The current experiment is sensitive to heterogeneity on the 0.3-100mm length scale, and heterogeneity spectra of carbonates are shown to vary widely from sample to sample over this range. Thus, this analysis is shown to provide a more detailed description of these porous media than the variogram or the first two moments of the porosity distribution provide. The magnetic resonance aspects of this technique are described here, while details of the geostatistical methodology are presented in a companion paper [A.E. Pomerantz, P.G. Tilke, Y.-Q. Song, Math. Geosci., submitted for publication].

Magnetic resonance in porous media: recent progress.

Recent years have seen significant progress in the NMR study of porous media from natural and industrial sources and of cultural significance such as paintings. This paper provides a brief outline of the recent technical development of NMR in this area. These advances are relevant for broad NMR applications in material characterization.

Construction and calibration of an instrument for three-dimensional ion imaging.

We describe a new instrument based on a delay-line detector for imaging the complete three-dimensional velocity distribution of photoionized products from photoinitiated reactions. Doppler-free [2+1] resonantly enhanced multiphoton ionization (REMPI) of H and D atoms formed upon photolysis of HBr and DBr in the range 203 nm < or = lambda photolysis < or = 243 nm yields radial speeds measured to be accurate within 1% of those calculated. The relative speed resolution is about 5% and limited by photoionization recoil broadening. A relative speed resolution of 3.4% is obtained for [3+1] REMPI, which minimizes the ionization recoil. We also determine the branching ratio between ground-state and spin-orbit-excited product channels and their associated anisotropies. We find that DBr photolysis dynamics differs slightly from its HBr counterpart.

Collision-energy dependence of HD(nu'=1,j') product rotational distributions for the H + D2 reaction.

Product rotational distributions for the reaction H + D2 --> HD(nu'=1,j') + D have been measured for 16 collision energies in the range of 1.43 < or = E(coll) < or = 2.55 eV. Time-dependent quantum-mechanical calculations agree well in general with the experimental results, but they consistently yield slightly colder distributions. In terms of the average energy channeled into rotation, the differences between experiment and theory amount to approximately 10% for all collision energies sampled. No peculiarity is found for E(coll)=2.55 eV at which the system has sufficient energy to access the first HD2 electronically excited state.

Reaction products with internal energy beyond the kinematic limit result from trajectories far from the minimum energy path: an example from H + HBr --> H2 + Br.

The importance of reactive trajectories straying far from the minimum energy path is demonstrated for the bimolecular reaction H + HBr --> H2(v', j') + Br at 53 kcal/mol collision energy. Product quantum state distributions are measured and calculated using the quasi-classical trajectory technique, and the calculations indicate that highly internally excited H2 products result from indirect reactive trajectories with bent transition states. A general argument is made suggesting that reaction products with internal energy exceeding a kinematic constraint can, in general, be attributed to reactive collisions straying far from the minimum energy path.

Observation and interpretation of a time-delayed mechanism in the hydrogen exchange reaction.

Extensive theoretical and experimental studies have shown the hydrogen exchange reaction H+H2 --> H2+H to occur predominantly through a 'direct recoil' mechanism: the H--H bonds break and form concertedly while the system passes straight over a collinear transition state, with recoil from the collision causing the H2 product molecules to scatter backward. Theoretical predictions agree well with experimental observations of this scattering process. Indirect exchange mechanisms involving H3 intermediates have been suggested to occur as well, but these are difficult to test because bimolecular reactions cannot be studied by the femtosecond spectroscopies used to monitor unimolecular reactions. Moreover, full quantum simulations of the time evolution of bimolecular reactions have not been performed. For the isotopic variant of the hydrogen exchange reaction, H+D2 --> HD+D, forward scattering features observed in the product angular distribution have been attributed to possible scattering resonances associated with a quasibound collision complex. Here we extend these measurements to a wide range of collision energies and interpret the results using a full time-dependent quantum simulation of the reaction, thus showing that two different reaction mechanisms modulate the measured product angular distribution features. One of the mechanisms is direct and leads to backward scattering, the other is indirect and leads to forward scattering after a delay of about 25 femtoseconds.

Disagreement between theory and experiment in the simplest chemical reaction: collision energy dependent rotational distributions for H + D2 --> HD(nu' = 3,j') + D.

We present experimental rotational distributions for the reaction H + D2 --> HD(nu' = 3,j') + D at eight different collision energies between 1.49 and 1.85 eV. We combine a previous measurement of the state-resolved excitation function for this reaction [Ayers et al., J. Chem. Phys. 119, 4662 (2003)] with the current data to produce a map of the relative reactive cross section as a function of both collision energy and rotational quantum number (an E-j' plot). To compare with the experimental data, we also present E-j' plots resulting from both time-dependent and time-independent quantum mechanical calculations carried out on the BKMP2 surface. The two calculations agree well with each other, but they produce rotational distributions significantly colder than the experiment, with the difference being more pronounced at higher collision energies. Disagreement between theory and experiment might be regarded as surprising considering the simplicity of this system; potential causes of this discrepancy are discussed.

Collision energy dependence of the HD(nu' = 2) product rotational distribution of the H + D2 reaction in the range 1.30-1.89 eV.

An experimental and theoretical investigation of the collision energy dependence of the HD(nu' = 2,j') rotational product state distribution for the H + D2 reaction in the collision energy range of Ecol = 1.30-1.89 eV has been carried out. Theoretical results based on time-dependent and time-independent quantum mechanical methods agree nearly perfectly with each other, and the agreement with the experiment is good at low collision energies and very good at high collision energies. This behavior is in marked contrast to a previous report on the HD(nu' = 3,j') product state rotational distribution [Pomerantz et al., J. Chem. Phys. 120, 3244 (2004)] where a systematic difference between experiment and theory was observed, especially at the highest collision energies. The reason for this different behavior is not yet understood. In addition, this study employs Doppler-free spectroscopy to resolve an ambiguity in the E, F-X resonantly enhanced multiphoton ionization transition originating from the HD(nu' = 2,j' = 1) state, which is found to be caused by an accidental blending with the transition coming from the HD(nu' = 1,j' = 14) state.

Rovibrational product state distribution for inelastic H+D2 collisions.

Experimental measurements of rovibrational product state distributions for the inelastic scattering process H + D2(nu=0,j)-->H + D2(nu' = 1,2,j') are presented and compared with the results of quasiclassical and quantum mechanical calculations. Agreement between theory and experiment is almost quantitative. Two subtle trends are found: the relative amount of energy in product rotational excitation decreases slightly with increasing collision energy and increases slightly with increasing product vibrational excitation. These trends are the reverse of what has been found for reactive scattering in which the opposite trends are much more pronounced.