Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Formation and emission of large furans and oxygenated hydrocarbons from flames

K. Olof Johansson, Tyler Dillstrom, Matteo Monti, Farid El Gabaly, Matthew F. Campbell, Paul E. Schrader, Denisia M. Popolan-Vaida, Nicole K. Richards-Henderson, Kevin R. Wilson, Angela Violi, and Hope A. Michelsen
PNAS July 26, 2016 113 (30) 8374-8379; first published July 7, 2016; https://doi.org/10.1073/pnas.1604772113
K. Olof Johansson
aCombustion Research Facility, Sandia National Laboratories, Livermore, CA 94550;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tyler Dillstrom
bDepartment of Mechanical Engineering, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Matteo Monti
cDepartment of Materials Science & Engineering, Stanford University, Stanford, CA 94305;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Farid El Gabaly
dMaterials Physics, Sandia National Laboratories, Livermore, CA 94550;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Matthew F. Campbell
aCombustion Research Facility, Sandia National Laboratories, Livermore, CA 94550;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul E. Schrader
aCombustion Research Facility, Sandia National Laboratories, Livermore, CA 94550;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Denisia M. Popolan-Vaida
eChemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720;
fDepartment of Chemistry, University of California, Berkeley, CA 94720;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nicole K. Richards-Henderson
eChemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kevin R. Wilson
eChemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Angela Violi
bDepartment of Mechanical Engineering, University of Michigan, Ann Arbor, MI 48109;
gDepartment of Chemical Engineering, University of Michigan, Ann Arbor, MI 48109;
hDepartment of Macromolecular Science and Engineering, University of Michigan, Ann Arbor, MI 48109;
iBiophysics Program, University of Michigan, Ann Arbor, MI 48109
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: avioli@umich.edu hamiche@sandia.gov
Hope A. Michelsen
aCombustion Research Facility, Sandia National Laboratories, Livermore, CA 94550;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: avioli@umich.edu hamiche@sandia.gov
  1. Edited by Barbara J. Finlayson-Pitts, University of California, Irvine, CA, and approved June 7, 2016 (received for review March 22, 2016)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Furans and related large oxygenated organic carbon species (OC) are highly toxic pollutants. Their integration into soot particles may greatly enhance soot’s hygroscopicity, leading to regional and global climate change. We show that furans are the primary oxygenated functional group on soot formed in hydrocarbon combustion and report a reaction scheme that elucidates the interplay between nonoxygenated and oxygenated hydrocarbons. We expect this reaction pathway to be important in many hydrocarbon oxidation systems spanning geosciences, astrophysics, and energy research. We discovered ∼100 oxygenated species previously unaccounted for in hydrocarbon models. This study advances the understanding of the oxidation chemistry of OC, which is critical to many processes, from controlling emissions of toxic combustion by-products to reducing anthropogenic climate change.

Abstract

Many oxygenated hydrocarbon species formed during combustion, such as furans, are highly toxic and detrimental to human health and the environment. These species may also increase the hygroscopicity of soot and strongly influence the effects of soot on regional and global climate. However, large furans and associated oxygenated species have not previously been observed in flames, and their formation mechanism and interplay with polycyclic aromatic hydrocarbons (PAHs) are poorly understood. We report on a synergistic computational and experimental effort that elucidates the formation of oxygen-embedded compounds, such as furans and other oxygenated hydrocarbons, during the combustion of hydrocarbon fuels. We used ab initio and probabilistic computational techniques to identify low-barrier reaction mechanisms for the formation of large furans and other oxygenated hydrocarbons. We used vacuum-UV photoionization aerosol mass spectrometry and X-ray photoelectron spectroscopy to confirm these predictions. We show that furans are produced in the high-temperature regions of hydrocarbon flames, where they remarkably survive and become the main functional group of oxygenates that incorporate into incipient soot. In controlled flame studies, we discovered ∼100 oxygenated species previously unaccounted for. We found that large alcohols and enols act as precursors to furans, leading to incorporation of oxygen into the carbon skeletons of PAHs. Our results depart dramatically from the crude chemistry of carbon- and oxygen-containing molecules previously considered in hydrocarbon formation and oxidation models and spearhead the emerging understanding of the oxidation chemistry that is critical, for example, to control emissions of toxic and carcinogenic combustion by-products, which also greatly affect global warming.

  • furans
  • oxygenated hydrocarbons
  • soot
  • organic carbon
  • black carbon

Oxygenated hydrocarbons produced during combustion can have a wide range of detrimental effects on human health, air quality, and regional and global climate. Furans, for example, are compounds that contain five-membered rings with four carbon atoms and one oxygen atom. They are frequently observed in the exhaust plumes and nearby environment of combustion sources. Many studies have shown that they are toxic, whether ingested or inhaled, and thus pose a considerable threat to human health (1⇓⇓–4). The simplest of these compounds (i.e., unsubstituted furan, C4H4O) is a cyclic, dienic ether with a low molecular weight, high volatility, and high lipophilicity. Studies on rats and mice have shown a dose-dependent increase in hepatocellular adenomas and carcinomas, indicating that furan is carcinogenic (4), and furan is marked as a high-priority substance and a carcinogenic risk by the World Health Organization (5).

Combustion sources of furans include biomass burning (6⇓⇓–9), cigarette and pipe smoke (10, 11), waste incineration (12), electronic waste recycling (13, 14), and volcanic activity (15). The polychlorinated dibenzofurans (PCDFs) are among the most notorious environmental pollutants, and the main source of PCDFs is biomass burning (6⇓⇓–9, 16). Nonchlorinated furans and PCDFs have been shown to be kinetically linked (17, 18). Furans released during combustion are often partitioned into particles and are found in ash from peat (9) and wood (6) burning, in primary organic aerosols from meat cooking (19), and in secondary organic aerosols from hydrocarbon oxidation (20, 21). Wood burning for heating and cooking constitute a major human exposure to airborne particulate PCDFs in some parts of the world (22, 23).

Previous work has suggested that oxygenated species can be attached to surfaces of soot particles of varying maturity emitted from flames and diesel engines, even before atmospheric processing (24⇓⇓⇓⇓⇓⇓⇓–32). Functional groups that have been identified include alcohols/enols, carbonyls, peroxies, and ethers. Oxygen atoms bound to organic species on the particle surface have been shown to greatly affect soot hygroscopicity (28) and the ability of soot particles to adsorb atmospheric water vapor and act as cloud-condensation or ice nuclei. Soot particles emitted from combustors, such as diesel engines, are generally hydrophobic, and enhancements in hygroscopic particle emissions could have substantial indirect climate effects via their influence on cloud formation (33). The effect of soot emissions on cloud-nucleation properties is a major uncertainty in climate predictions (34⇓–36).

Despite the impact of large oxygenated hydrocarbons on combustion chemistry, the environment, and human health, very little is known about their formation mechanisms and emissions. In this paper we present evidence of the formation of oxygenated compounds, including furans, during the combustion of hydrocarbon fuels. Via a synergistic approach that includes ab initio methods and a stochastic model in conjunction with experimental measurements, we identify reaction pathways leading to formation of oxygenated compounds during the combustion of ethylene. We recorded aerosol mass spectra sampled from premixed and diffusion flames, using synchrotron-generated vacuum-UV (VUV) radiation for ionization, for comparison with masses of the predicted chemical compositions. The mass spectra show masses of oxygenated species that agree with the atomic compositions predicted by the simulations. Both experiments and simulations demonstrate that ∼50% of the mass peaks observed at some flame heights in the mass range 140–350 u (unified atomic mass units) contain signal from oxygenated species. We also recorded X-ray photoelectron spectroscopy (XPS) spectra of soot samples extracted from these flames for further validation of these mechanisms by comparison with functional groups of the predicted oxygenated species incorporated into particles. The XPS measurements confirm formation of furan precursors, hydroxyl groups, early in the soot-formation process and evolution of furan signatures, ether groups, as the combustion and particles evolve.

The present study represents an important step toward the development of predictive models for the oxidation of hydrocarbons, which will require that the presence and reactivity of these oxygenated compounds are taken into account. Understanding the chemistry related to high-temperature hydrocarbon oxidation may provide a key to controlling emissions of harmful combustion byproducts, such as soot, nonchlorinated furans, and PCDFs, leading to multiple environmental and health benefits.

Results and Discussion

Simulations.

We investigated new reaction pathways for the formation and decomposition of oxygenated species using ab initio electronic-structure calculations to elucidate how oxygen becomes incorporated into the carbon framework of organic carbon (OC) species and soot. We found a generic furan-formation route with low reaction barriers involving species that are abundant in hydrocarbon flames, such as acetylene. The reaction pathway is likely to end with a unimolecular ring closure, suggesting that there should be a high fraction of furans present among oxygenated OC formed during combustion. Fig. 1 illustrates the barriers (TS1 and TS2) of the unimolecular furan-ring closure by showing the formation of benzofuran from phenoxy-acetylene. The energy barriers are relatively low compared with the average temperature fluctuations in combustion environments.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Potential-energy diagram for the formation of benzofuran. The energies are in kilocalories per mole and are referenced to the reactant species. The ring-closing reaction is exothermic, and the closed-ring structure is thus favored over the open-ring structure. The energy barriers are low compared with the average temperature fluctuations in flames. For example, at 1,500 K, TS1 is 3kBT and TS2 is −1kBT with respect to the reactant species.

Reaction rates for different oxygen chemistry pathways were calculated and incorporated into the Stochastic NAnoParticle Simulator (SNAPS) (37) along with kinetic rates from the literature that accounted for different oxygen-adsorption reactions and ring-closure reactions. SNAPS is a probabilistic code that predicts the formation of different soot-precursor molecules and other OC species in flames based on individual reaction trajectories of initial seed molecules as they evolve in the flame.

A premixed C2H4/O2 flame, in which fuel and oxygen are mixed before ignition, was chosen for the SNAPS simulations. SNAPS yielded six main classes of oxygenated groups: alcohol/enol, peroxy acid/radical, ketene, pyran, noncyclic ether, and furan. Alcohols/enols and peroxy acids and radicals are formed via reactions of hydrocarbons with small oxygen-containing molecules, mainly OH and O2, and serve as precursors to the other four oxygenated groups. Enols containing 6–14 carbon atoms (Fig. S1) account for ∼20% of the oxygenated structures close to the burner, where oxygen-containing radicals are present, demonstrating that large enols are important for the chemistry of oxygenated hydrocarbons, and extending the studies of small enols by Taatjes et al. (38) (Important Enol Intermediates Predicted by SNAPS). We found that, in addition to the isomerization reaction, removal of enols may proceed via reactions with flame radicals, particularly H and OH, producing oxyradicals. Production of large enols can proceed through radical attacks on the carbon–carbon double bond of enols, followed by alkyl radical addition.

Fig. S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S1.

Enols predicted by the SNAPS simulations. Enols predicted to be important to the OC chemistry at small DFFOs in premixed combustion are (A) 94 u (phenol), (B) 108 u (C7H8O), (C) 118 u (C8H6O), (D) 168 u (C12H8O), and (E) 192 u (C14H8O).

The SNAPS simulations showed that ketenes are formed as the result of oxyradicals present on the terminal carbon of aliphatic chains. Noncyclic ethers are produced when hydrocarbons add to oxyradical sites, following three evolution pathways: (i) decomposition, (ii) continued hydrocarbon addition, or (iii) ring-closure reactions to form furan/pyran groups. Pyran groups are rare, forming when propargyl or methyl and acetylene add to oxyradicals on free-edge sites of aromatic rings, or when acetylene adds to oxyradicals on zigzag sites (see Fig. S2 for site definition). The results also showed potential precursor molecules for formation of dioxins, species with six-membered rings including two oxygen atoms. However, the code does not include any elementary reactions that close a ring containing two oxygen atoms.

Fig. S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S2.

Definition of a zigzag site and a free-edge site.

Our simulation results demonstrate that oxygen addition to polycyclic aromatic hydrocarbons (PAHs) occurs mainly through consecutive bimolecular reactions resulting in H-abstraction followed by addition of OH or O2 onto PAH edge sites (Comments on the SNAPS Modeling). The most probable pathway for embedding oxygen into the hydrocarbon molecules is via ethers formed when H is abstracted from hydroxyl groups or OH/O is abstracted from peroxy acid/radical groups, followed by hydrocarbon addition to the oxyradical and furan-ring closure. Acetylene is the most frequently added hydrocarbon, and Fig. 2 shows the most probable reaction sequence leading to formation of a furan group. We expect this scheme to be important to a wide range of hydrocarbon oxidation processes and virtually any hydrocarbon-combustion system because of high acetylene concentrations and low reaction barriers. Thus, apart from the toxic potential of furan species, furans may also play a dominant role in determining growth and oxidation sites, radiative forcing potential, and hygroscopic properties of the OC formed during combustion. As these oxygenated OC species evolve in the flame, they can become large enough to condense onto incipient soot particles, leading to incorporation of oxygen onto the surface of soot. Oxygen embedded in the particle surface is expected to contribute to gaseous oxidation products and affect the further growth and oxidation of the particle (17).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Most probable reaction sequence leading to formation of a furan group. (Left to Right) H-abstraction followed by OH addition to the radical free-edge site on an aromatic ring; H-abstraction from the OH group, followed by acetylene addition, forming an ether group; H-elimination during ring closure to form a furan group. The left side of the molecule has been left attached to an indeterminate PAH backbone to illustrate an arbitrary molecular size (see Comments on the SNAPS Modeling for details on OH addition).

Fig. 3 shows some of the most common structures obtained from SNAPS at 160, 168, 194, and 220 u low in the premixed flame. Many of the ether structures formed are furan precursors. Competing reactions to furan-ring closure in the final step in Fig. 2 includes addition of a second hydrocarbon species (Fig. 3 A and B) or H-addition forming an R–O–CH=CH2 group (Fig. 3 D and F) or acetylene removal. The most common oxygenated structure at 192 u is the 194-u structure shown in Fig. 3D with a furan ring instead of the –OCHCH2 chain; the most common structure at 158 u is the furan formed with the oxygen attached to the phenyl ring in Fig. 3B (Information Related to Fig. 3). The species in Fig. 3B has a ketene group at the terminal carbon of the aliphatic branch, showing that, if acetylene adds to a radical site and OH or O2 adds to the acetylene chain, a ketene may form instead of a furan and that aliphatic side chains provide sites for the formation of carbonyl bonds.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Frequently predicted oxygen-containing structures of selected masses low in a premixed flame. Red atoms, oxygen; gray, carbon; white, hydrogen. (A and B) Ether and ether/ketene, 160 u. (C) Furan, 168 u. (D and E) Ether and furan, 194 u. (F) Ether, 220 u.

The evolution of the structure in Fig. 3F frequently proceeds via acetylene loss from the oxygen atom or formation of a pyran ring with the zigzag-site carbon, yielding a structure with mass 218 u. Many of the furans formed are potential PCDF precursors (e.g., the substituted dibenzofuran in Fig. 3E), although naphtho[2,1-b]furan is more common in the simulations than dibenzofuran at 168 u (Fig. 3C). The majority of the species predicted by SNAPS experienced molecular growth as they evolved in the flame, and furans constituted the largest group of oxygenated species at large distance from the fuel outlet (DFFO), where the flame temperature was around 1,750 K. An important furan-destruction pathway identified in the simulations is CO reactions: CO can open furan rings and abstract the oxygen atom to form CO2, leaving the original furan ring as an aliphatic radical side chain. This mechanism may explain the observed aliphatics on the surface of soot particles under some conditions (25, 27, 31, 39).

Aerosol Mass Spectrometry.

To support and validate our modeling results, we studied two McKenna-type premixed C2H4/O2 flames with equivalence ratios matching the flame studied using SNAPS and a counterflow diffusion (CF) flame (Burners and Gas Flows). In premixed combustion the oxidizer is mixed with the fuel before reaching the flame front, and the concentration of volatile oxygenated hydrocarbons is expected to peak at small DFFOs where OH and O2 are present. In the CF flame, fuel and oxidizer are fed separately via counterpropagating flows; mixing of fuel and oxidizer takes place across the gas-stagnation plane formed by the counterpropagating flows. The CF flame is an excellent flame for comparison with premixed combustion because there is no oxygen available at small DFFOs, and large oxygenated structures cannot be present at small DFFOs. Hydrocarbon growth initially takes place under pyrolytic conditions, and particles are horizontally convected away from the flame parallel to the gas-stagnation plane, where the vertical velocity component vanishes.

We performed synchrotron-coupled VUV aerosol mass spectrometry (AMS) of soot particles drawn from the different flames and injected directly into the mass spectrometer using an aerodynamic lens system (ADLS) (see AMS for more details). Large gas-phase species can condense onto particles in the sampling line between the flame and the ADLS, which leads to physical growth of the particles while the chemistry is frozen, allowing small particles to reach sizes that can be efficiently focused by the ADLS. The particle growth also permits us to detect large gas-phase OC species that are prone to condense onto soot particles at atmospheric temperatures but not at flame temperatures. Species condensed onto soot particles are vaporized in the ionization region of the mass spectrometer without changing species internal composition and mass. Fig. 4 shows example mass spectra of species drawn from three DFFOs in a premixed flame. The results are consistent with the SNAPS simulations and demonstrate that a remarkably large fraction (∼50%) of the mass peaks is associated with oxygen inclusion (Determination of Precise Masses and Atomic Constituents, Table S1, Nominal Masses of Oxygen-Containing Species Determined from the Aerosol Mass Spectra, and Table S2). Mass peaks identified as being associated with oxygenated species are highlighted in red in Fig. 4. The masses of these species (Table S2) agree with the masses of oxygenated OC species predicted by SNAPS, that is, the experimentally observed and predicted species have the same atomic compositions (Table S1). Our AMS studies revealed that these oxygenated species are present under a very wide range of combustion conditions, suggesting generic formation mechanisms.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

AMS spectra from a premixed flame. Mass spectra are shown for particles extracted from selected heights in the flame; that is, DFFOs of (A) 3.5, (B) 5.0, and (C) 7.5 mm. Red peaks contain signal from oxygenated species. The arrows indicate the peaks at 160, 194, and 220 u for comparison with Fig. 5A; see Fig. 3 for main predicted structures at a DFFO of ∼3.5 mm.

View this table:
  • View inline
  • View popup
Table S1.

Analysis of some mass peaks

View this table:
  • View inline
  • View popup
Table S2.

Nominal masses of oxygen-containing molecules observed in the flames

Fig. 4 shows that the signal from oxygenated OC gets weaker at larger DFFOs in premixed flames, which indicates that they are only formed at low DFFOs in premixed combustion (i.e., where small oxygen-containing species, such as OH and O2, are present). The SNAPS simulations showed that oxygen addition to PAHs occurred predominantly below a DFFO of 2 mm (Fig. S3), where the concentrations of reactive oxygen (e.g., OH and O2) are significant in these flames. The soot particles are, however, too small below a DFFO of 3.4 mm to yield signal in the mass spectrometer (see AMS for details). Hence, the mass spectrum in Fig. 4A was recorded at a DFFO of ∼3.5 mm. The SNAPS simulations predicted that the oxygenated species continue to grow as they evolve in the flame, which is consistent with the signal decrease for the red peaks in Fig. 4 with increasing DFFO. As their sizes increase and they become less volatile, they are likely to condense onto incipient soot particles and become integrated into the soot, which further reduces their volatility and prevents them from being vaporized in the ionization region of the AMS instrument. Hence, large oxygenated species are expected to avoid detection at large DFFOs.

Fig. S3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. S3.

Predicted H and OH mole fractions at low flame heights in the Ar-diluted premixed flame using the Appel–Bockhorn–Frenklach (47) mechanism. The ratio between the H and OH mole fractions is also shown. The vertical dashed lines are the edges of the region where the majority of the oxygenation occurs, according to the SNAPS simulations.

Fig. 5A shows the ion-signal DFFO dependence for 160, 194, and 220 u for a premixed flame (Fig. 3 displays major predicted structures for these masses). These three masses are indicated by arrows in Fig. 4. Their spatial profiles peak at the lowest DFFO sampled, and we are unable to experimentally establish the location of maximum concentrations because of lack of AMS signal at lower DFFOs. Profiles inferred from mass peaks identified by the SNAPS simulations to stem from pure (nonoxygenated) hydrocarbons, for example C18H10 isomers (226 u) and C20H10 isomers (250 u), have low intensities at small DFFOs for this flame. Fig. 5B shows that the same trend with reactive oxygen was obtained from the CF flame. In the CF flame the growth of hydrocarbon species initially takes place on the fuel side of the flame under pyrolytic conditions, and the reaction scheme in Fig. 2 cannot occur at small DFFOs. We identified many of the same oxygenated masses in the CF flame as in the premixed flames (e.g., 160, 194, and 220 u) (Table S2), and, as expected, they all appear skewed toward the oxygen side of the flame compared with masses identified as pure hydrocarbon masses (e.g., 226 and 250 u), as shown in Fig. 5B.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Normalized signal profiles at selected masses for (A) a premixed flame and (B) the CF flame. Solid lines refer to masses that include oxygenated species; dashed lines show profiles of nonoxygenated hydrocarbons. For clarity, every third symbol is shown in B.

The profiles at 194 and 220 u reach minima at ∼6.5 mm DFFO in Fig. 5A and then increase with DFFO. The intensity of the 160-u mass peak does not increase at large DFFOs. Analysis of the exact masses of the 194- and 220-u peaks shows that these two peaks shift slightly toward higher masses with increasing DFFO. These shifts in mass-peak locations indicate that pure hydrocarbons are formed at these two nominal masses at large DFFOs, because the combined mass of one carbon and four hydrogen atoms is slightly larger than the mass of one oxygen atom. These nonoxygenated hydrocarbons increase with increasing DFFO; this observation is consistent with previous gas-phase PAH measurements in sooting premixed flames (e.g., ref. 40). A similar result is obtained in the CF flame (Fig. 5B), where the profile at 160 u has a larger shift toward the oxygen side than the profiles at 194 and 220 u, indicating that the signal at 160 u stems entirely from oxygenated species whereas pure hydrocarbon structures contribute to the signals at 194 and 220 u.

The mass spectra in Fig. 4 are composed of clusters of peaks. The width of each cluster is a crude measure of the spread in hydrogenation among species containing the same number of carbon atoms, because each cluster is largely composed of peaks corresponding to OC species containing a specific number of carbon atoms. The low-mass side of the clusters contain peaks stemming from less-saturated (more aromatic) species without oxygen content. Highly saturated species and species with oxygen atoms appear on the high-mass side of the clusters because an oxygen atom weighs more than the most saturated hydrocarbon group for which it can substitute (i.e., CH3). The high-mass side of a cluster may also contain signal from species of very low saturation (e.g., odd-carbon-numbered clusters may include signal from the polyyne species containing one additional carbon atom). Hence, on the high-mass side of each cluster are peaks stemming from species that may contain a different number of carbon atoms than the species yielding peaks in the center and on the low-mass side of the cluster. Some oxygenated OC species have nominal masses identical to those of large polyynes and are likely to be misidentified in soot-sample aerosol mass spectra.

XPS.

The AMS study confirmed the atomic constituents of the oxygenated OC species predicted by the SNAPS simulations. To validate the predicted OC structures, information on how the oxygen is bonded in the OC molecules is needed. We thus performed XPS studies on soot samples extracted from the premixed flames, as shown in Fig. 6 (XPS). The XPS data verified the SNAPS prediction that alcohol/enol groups are formed low in the flame, and their concentrations decrease with increasing DFFO. The abundance of ethers included in soot particles increases with increasing DFFO, suggesting that alcohols/enols are precursors to ethers. The fractions of C-OH, C-O-C, and C=O species at different DFFOs were determined by analysis of the O 1s spectra (30, 41) and are shown in Fig. 6D. The low fraction of C=O species suggests that our soot samples had not been significantly affected by atmospheric oxygen between flame sampling and XPS measurements (42). The strong XPS signal from oxygen-containing functional groups at large DFFOs suggests that the reduced signal from the mass peaks correlated with oxygen inclusion in the mass spectra recorded at DFFOs of 5.0 mm and 7.5 mm (Fig. 4 B and C) compared with the mass spectrum recorded at a DFFO of 3.5 mm (Fig. 4A) is not due to lack of oxygenated OC in the soot. This observation is consistent with incorporation of oxygenated OC in the soot at larger DFFOs, as predicted based on the SNAPS simulations and the AMS studies.

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

XPS O 1s spectra of soot sampled from a premixed flame. Spectra are shown for particles extracted from (A) a low flame height (DFFO of 3.5 mm), (B) an intermediate flame height (DFFO of 5.0 mm), and (C) high in the flame (DFFO of 7.5 mm). (D) Fractional contributions of oxygenated functional groups were inferred from fits to these and similar XPS O 1s spectra at other DFFOs and plotted as a function of DFFO. Error bars represent ± 1 SD of the uncertainties associated with the peak fits.

Summary and Conclusions

Despite the distinctly different flame conditions for premixed and diffusion-controlled combustion, we observed consistent masses and behavior regarding formation of oxygenated OC, supporting the generality of large-oxygenate formation, as described in Fig. 2. These reaction pathways are expected to have substantial impact on soot properties, such as hygroscopicity and toxicity. The apparent generality of these chemical mechanisms also suggests that they are important to hydrocarbon oxidation systems other than combustion and may be important in geophysics, astrophysics, and energy research.

We discovered ∼100 oxygenated species previously unaccounted for in hydrocarbon-mass-growth and oxidation models. Under the present configurations, oxygen addition to PAHs occurs mainly via reactions with OH and O2. Enols and alcohols that result from these reactions act as precursors to ethers, such as furans and pyrans. Furans are produced in the high-temperature regions of hydrocarbon flames, and they remarkably survive and become the main functional group of oxygenates in the OC on soot particles. Cyclic ethers therefore play a more important role than carbonyl and hydroxyl groups during some phases of organic carbon oxidation and may require inclusion in many hydrocarbon oxidation mechanisms and in models aimed at predicting toxicity and hygroscopicity of soot particles and related OC.

The oxygenated OC species discovered in the present study seem to readily become partitioned into incipient soot particles, the lighter ones at atmospheric temperatures and the heavier ones also at flame temperatures. These species should therefore be integral to the molecular-growth pathways leading to soot formation as well as to OC coating of soot particles in the combustor and in the atmosphere. The emission of oxygenated OC and soot can be minimized by an oxygen-rich postcombustion zone at an elevated temperature that should be kept below the threshold for NOx formation. Furans can be removed through CO reactions as identified by the SNAPS simulations.

Further studies are needed to improve the knowledge of carbon–oxygen interactions in the oxidation and growth chemistry discussed here. Understanding how soot forms and interacts with OC can lead to more efficient combustion technologies that reduce soot emissions, which would have many benefits (e.g., reduced global warming and anthropogenic hydrological impact and reduced concentrations of toxic and carcinogenic environmental pollutants). Flame studies of soot and related OC are, however, also important because they can provide benchmarks for hydrocarbon growth and oxidation models important in other research fields, such as astrophysics (43) and fuel cells (44).

Methods

Computational Approach.

Chemical kinetics simulations were performed using SNAPS (37, 45) (Comments on the SNAPS Modeling). Briefly, SNAPS is Monte Carlo-based and simulates trajectories of reversible reactions for given seed molecules. General insights into the growth/consumption processes require analysis of a large number of trajectories. SNAPS necessitates information on temperature and gas-phase species concentrations. These inputs were computed by solving the gas-energy equation using the PREMIX program in CHEMKIN (46) with the gas-phase mechanism from Appel et al. (47) because it best represented an average of the speciation profiles generated when using four well-validated gas-phase mechanisms [Appel et al. (47), Miller and Melius (48), Richter et al. (49), and Raj et al. (50)]. Benzene and toluene were used as seed molecules because formation of the first aromatic ring is considered to be the first step in soot formation (51). Molecules formed in the simulations were classified by carbon configuration using social permutation invariant topological coordinates (52); see ref. 37 for more details.

Electronic-structure calculations were carried out using the CBS-QB3 method (53) as implemented in Gaussian 09 (54). Singlet, doublet, and triplet spin multiplicities were tested for the reactant, and the lowest energy (doublet) was determined to be the ground state and thus used for all subsequent calculations. Intrinsic reaction-coordinate calculations were carried out at the B3LYP/CBS level to ensure that each optimized transition state connected the expected reactant and product.

Experimental Approach.

Information about the burners and the gas flows are available in Burners and Gas Flows. Briefly, the premixed burner is a McKenna design with a 38.1-mm sintered-brass plug surrounded by a shroud-gas sintered-brass ring. The counterflow diffusion burner consists of two vertically mounted, central flow tubes, facing each other, 12 mm apart. They have a 12.7-mm inner diameter and are surrounded by outer flow tubes that help shield the flame from the surroundings. Fuel and oxidizer were fed separately via counterpropagating flows; mixing of fuel and oxidizer took place across the gas-stagnation plane formed by the counterpropagating flows. Hydrocarbon growth initially takes place under pyrolytic conditions, and particles are horizontally convected away from the flame parallel to the gas-stagnation plane, where the vertical velocity component vanishes.

Soot particles were sampled along the vertical centerline of the flames using a quartz probe with a tapered tip. The probe assembly is fixed, and, to sample from different DFFOs, the burner assembly is translated vertically. In the AMS measurements, the sample enters an ADLS (55, 56), which focuses particles larger than ∼50 nm diameter into a beam. Condensation in the sampling line leads to particle growth, which enables some particles that are originally smaller than 50 nm to be focused by the ADLS onto a target heated to ∼570 K inside a vacuum chamber (∼7 × 10−7 Torr). Species that vaporize from the heated target are photoionized using VUV radiation at 9.6 eV, generated at the Advanced Light Source (ALS) synchrotron facility at Lawrence Berkeley National Laboratory. The molecular ions generated are pulse-extracted into a time-of-flight mass spectrometer at a rate of 15 kHz, and mass spectra are recorded using a multichannel scaler. Additional information about the aerosol mass spectrometry investigation is available in AMS.

Soot samples for the XPS measurements were collected by extracting particles from the N2-diluted premixed flame using quartz microprobes from the same batch of probes that was used for the AMS measurements. The particles were drawn from different points along the vertical centerline of the flame at the same flow rate as was used during the AMS studies. The particles were passed through an aerosol neutralizer (TSI 3087) and collected on clean disks of pure aluminum (99.5 pure, 12.5-mm diameter; Goodfellow) using a nanometer aerosol sampler (TSI 3089).

The XPS measurements were performed under ultra-high vacuum conditions (residual pressure <1 × 10−9 Torr) using an Omicron DAR400 Mg K-alpha X-ray source and Physical Electronics 10–360 electron energy analyzer. The spectra shown have not been corrected or processed apart from a background subtraction. The calibration of our instrument was confirmed by measuring a clean sample of highly oriented pyrolytic graphite‎, which gives a sharp C 1s peak centered at 284.4 eV and FWHM of 1.2 eV, in good agreement with the literature (30, 41). The fitting process, quantitative analysis, and error estimation were performed using the CasaXPS software. All XPS peaks were baseline-corrected using a Shirley background, and a mixed Gaussian–Lorentzian (70%/30%) line shape was used to fit the spectra. The O 1s photoemission spectra were deconvolved using three components corresponding to the oxygen functional groups: C-OH (533.2 eV), C-O-C (531.9 eV), and C=O (530.2 eV), according to Müller et al. (30).

Details About the Approaches

Burners and Gas Flows.

The premixed burner was a McKenna design with a 38.1-mm sintered-brass plug surrounded by a shroud-gas sintered-brass ring. Two different water-cooled stabilization plates were used: a square aluminum plate with 100-mm side length and a stainless-steel cylinder with 38.1-mm diameter. The aluminum plate was located 12 mm above the sintered-brass plug and was maintained at ∼335 K. The stainless-steel cylinder was mounted 13.3 mm above the plug, and the surface temperature of the cylinder was kept at ∼400 K. The gas-flow rates through the sintered-brass plug used in conjunction with the aluminum plate were 0.75 L/min (lpm) C2H4, 1.10 lpm O2, and 3.40 lpm Ar. All gas flows are referenced to 273 K and 1.01325 × 105 Pa. The shroud-gas flow rate was 14 lpm Ar. The corresponding gas-flow rates used in combination with the stainless-steel cylinder were 0.67 lpm C2H4, 0.96 lpm O2, and 3.62 lpm N2 through the porous plug. The shroud-gas flow rate was 19 lpm N2.

The counterflow diffusion burner consists of two vertically mounted, central flow tubes, facing each other, 12 mm apart. They have a 12.7-mm inner diameter and a maximum outer diameter of 15.875 mm. The outer walls have a 9-mm taper to the edge at the outlets. Fuel and oxidizer were fed separately via counterpropagating flows; mixing of fuel and oxidizer took place across the gas-stagnation plane formed by the counterpropagating flows. Hydrocarbon growth initially took place under pyrolytic conditions, and particles were horizontally convected away from the flame parallel to the gas-stagnation plane, where the vertical velocity component vanished. The fuel-gas mixture was supplied through the lower flow tube and consisted of 0.23 lpm C2H4 and 1.10 lpm Ar. Thus, oxygenated structures could not form at small DFFOs because there was no oxygen available, which was verified experimentally. The flow through the upper tube was for the oxidizer, and it consisted of 0.25 lpm O2 diluted in 1.20 lpm Ar. Outer flow tubes with inner diameter 18.364 mm surround the central flow tubes. The outer flow tubes were used for Ar flows, 2.30 lpm through the lower one and 3.00 lpm through the upper one.

AMS.

Sample gas and soot particles were extracted from the flames using a quartz microprobe manufactured from a straight quartz tube with an o.d. of 2.9 mm and an i.d. of 2 mm. The last ∼5 mm of the probe end facing the flame was tapered down to a tip opening of ∼300 μm, and the probe drew sample from the flame at a rate of ∼70 standard cubic centimeters per minute (sccm). The relatively large opening diameter of the probe was necessary to avoid the soot particles generated in the present flames from clogging the probe during the measurements. The probe was slightly longer than 40 mm, because the outer radius of the premixed burner is 38.1 mm. The quartz probe was connected to an aluminum probe block using a Swagelok connection. The temperature of the probe block was close to ambient. The probe block is ∼25 mm long, and the i.d. of the sample channel inside the block is ∼8.5 mm. The probe block was connected to a ∼50-cm-long, 6.35-mm o.d., 4.57-mm i.d. stainless steel tube that was connected to the inlet nozzle of the ADLS on the AMS. The total residence time in the sampling system leading to the ADLS was ∼8 s. The ADLS has been characterized by Zhang et al. (55) and Headrick et al. (56). The pressure behind the inlet nozzle feeding the ADLS is about 2 Torr. The outlet of the ADLS is in a vacuum chamber (4.5-inch spherical cube from Kimball Physics) with an operating pressure that is lower than 1 mTorr.

The ADLS acts as a filter against gas-phase molecules and small particles; it focuses particles larger than ∼50 nm into a beam whereas very small particles and gas-phase molecules follow Brownian motion and are pumped away by the turbomolecular vacuum pumps. Heat transfer is a diffusive process and the sample temperature drops to room temperature within 3 cm in the sampling line. The chemistry can be considered to be frozen following extraction from the flame. Large hydrocarbon species that were in the gas phase in the flame condense onto soot particles in the sampling line, making the soot particles grow. This growth does not alter the atomic compositions of the different species, because no covalent bonds are formed or broken.

Because the ADLS acts like a particle filter, particle growth through condensation in the sampling line is favorable; it allows small particles to grow to sizes that are more efficiently focused by the ADLS, and it allows the detection of gas-phase species that are prone to condense onto particles at temperatures lower than in the flame. For example, molecules as small as naphthalene (C10H8, mass 128 u) are often solids at room temperature, but they are in gas phase at flame temperatures. Naphthalene is, however, too volatile to remain condensed on soot particles at the low pressures of the AMS and vaporizes and is pumped away before reaching the ionization region. Nevertheless, for several species with masses as low as ∼135 u (Fig. 4), a fraction of the molecules remain condensed all of the way to the ionization region and are detected by the instrument.

A majority of the present AMS spectra was recorded using an ionizing-photon energy of 9.6 eV. There are several gas-phase species with ionization thresholds below 9.6 eV (e.g., benzene, toluene, and naphthalene), and these should completely dominate the signal from the AMS if gas-phase molecular concentrations were not suppressed by several orders of magnitude. The fact that there is virtually no signal at these masses in the AMS supports the conclusion that volatile species were pumped out of the aerosol mass spectrometer after being guided through the ADLS.

The first vacuum chamber, which houses the ADLS, leads to the second vacuum chamber (4.5-inch spherical cube from Kimball Physics), which has an operating pressure in the 10−5 Torr range. The third chamber is the main chamber (custom 10-inch-diameter cylindrical vessel with multiple conflat ports), and it has an operating pressure of ∼7 × 10−7 Torr. The particles that were focused by the ADLS hit an aerosol target heated to ∼570 K and located in the center of the main chamber. The VUV beam passed ∼1.9 mm in front of the aerosol target and ionized species that had vaporized from the soot on the aerosol target. We have investigated the effect on the signal as a function of different aerosol target temperatures. Increasing the temperature beyond 570 K did not yield any notable new signal peaks on the high-mass side of the spectra but could potentially increase fragmentation.

The O2 concentrations were very low throughout the regions sampled in the three flames, which was verified experimentally and by our gas-phase calculations for the SNAPS simulations. The combination of low O2 concentrations in the flame locations probed and efficient removal of gas-phase molecules in the different pumping stages of the AMS meant that the O2 concentration in the ionization region was negligibly small and insufficient for any surface-catalyzed hydrocarbon–O2 reactions to influence the detected signals. Gas-phase OH was not a concern at the heated target because it was consumed in and around the quartz probe used for sampling and was therefore not even detected in later gas-phase molecular beam mass spectrometry studies performed in the CF flame.

The energy distribution of the ionizing photons used for the AMS studies was relatively broad (∼0.25 eV FWHM), which was not a problem for the present investigation. It is virtually impossible to perform definitive isomer assignments based on photoionization efficiency (PIE) scans (i.e., how the signals of individual mass peaks respond to the tuning of the photoionization energy) for species as large as those under consideration in the present study (57). The reasons are many: the number of isomers grows with the molecular size, the isomers of potential importance have small differences in heats of formation, and similar chemical structures lead to similar ionization energies and PIE curves (58). We have, however, performed PIE scans in the present flames in the range 7.6–10.5 eV. We recorded the signal response to the photon energy for all masses simultaneously, because the instrument used was a time-of-flight mass spectrometer. The PIE curves of the signals at different masses are generally very similar and the ionization thresholds are often below 7.6 eV, which is the lower photon-energy limit at the beamline where we performed the present investigation.

XPS.

In addition to recording highly averaged signals in the O 1s energy range, as shown in Fig. 6 A–C, we also acquired complete XPS scans between 0 eV and 1,200 eV. These broad scans were performed to verify that the signals that we obtained stemmed only from the hydrocarbon soot sample and that the samples were clean from impurities of other elements.

The XPS results showed that OH and ether groups are almost equally common in the nascent soot particles drawn at small DFFOs. As these particles evolved through the flame, oxygen atoms that were initially bound in OH groups became embedded into ether structures and the fraction of ether functional groups increased and the fraction of OH functional groups decreased with increasing DFFO. It is conceivable from the observed trend that OH functional groups are more abundant than ether groups at very low flame heights, strongly supporting the SNAPS model, which predicted that small oxygen-containing radicals such as OH attached to hydrocarbon structures at low flame heights and then evolved into ether groups (mainly furans).

Comments on the SNAPS Modeling.

The presently used version of SNAPS includes hydrocarbon growth through reactions of methyl, methane, acetylene, vinyl, ethylene, ethane, ethyl, propargyl, butatriene, cyclopentadienyl, cyclopentadiene, phenyl, and benzene. Oxygen chemistry is included through reactions with O, O2, OH, HO2, H2O2, H2O, CO, and CO2. For details concerning earlier versions of the model and methods used in the SNAPS code, consult ref. 37. New implementations to SNAPS were focused on the inclusion of pathways and kinetic rates describing oxygenation of PAHs, including but not limited to addition of hydrocarbons to oxyradicals and furan- and pyran-ring formations. We accounted for key kinetic pathways leading to oxygen incorporation into hydrocarbons. These pathways were oxygen adsorption reactions and ring-closure reactions, together accounting for the plurality of oxygen-embedment pathways.

The new reaction pathways for furan formation were incorporated into SNAPS via elementary steps. In the code, reaction sites are specified in terms of properties including the atom type (e.g., aromatic, aliphatic, sp3 hybridized), the bond type, and molecular connectivity [e.g., membership in ring(s) or a specific aliphatic group]. Reactions are defined in terms of these reaction sites and corresponding kinetic-rate constants. SNAPS assumes identical reaction-rate constants for furan-ring closure regardless of the overall molecular size and structure attached to the reactant aromatic ring, as long as the O‒CH=CH· chain is attached to a free-edge site (Fig. S2) of the aromatic ring.

Approximately 70 oxygenation/deoxygenation reactions were implemented in the SNAPS mechanism. Most of these reactions (80%) fall into three categories: bimolecular addition/removal of oxygenated groups to the target molecule, ether formation/destruction, and oxygenated ring closures/openings. The majority of the kinetic reaction rates were sourced from existing literature, and the remaining were estimated based on similarities to analogous reactions with known rates. The high frequency with which we observed furans in the initial simulations demonstrated the significance of the furan formation pathways. Therefore, we used ab initio techniques to validate the kinetic reaction rates of the furan-group formation. The CBS-QB3 method (53) used for the electronic structure calculations involves a sequence of five calculations: (i) geometry optimization at the B3LYP level and (ii) frequency calculation followed by single-point energy calculations at (iii) CCSD(T), (iv) MP4SDQ, and (v) MP2 levels. The extrapolation method unique to CBS then gives the final zero-point corrected energies (59). Intrinsic reaction coordinate calculations were carried out at the B3LYP/CBS level to ensure that each optimized transition state connected the expected reactant and product. Subsequent analysis revealed the importance of the intermediate enol and ether chemistry to oxygenated compounds observed at a DFFO of ∼3.5 mm.

The SNAPS code does not simulate nucleation, particle formation, and particle chemistry. The species suggested by SNAPS can therefore be compared with large gas-phase species in the flames and species that have condensed onto particle surfaces but not participated in surface chemistry on the particles. In addition to precursors to the PCDFs we also observed potential precursor molecules for formation of dioxins, species with six-membered rings including two oxygen atoms. Presently the SNAPS mechanism does not include any elementary reactions that close a ring containing two oxygen atoms. We ran several SNAPS simulations to assess the effect of gas-phase species profiles that had been shifted downstream by a few microprobe orifices. The qualitative SNAPS results remained the same when using nonshifted and shifted gas-phase species profiles.

In Fig. 2, an H atom is abstracted/eliminated from the aromatic ring before OH addition, because, at high temperatures, the prereaction complex formed from direct OH addition to an aromatic ring is unstable; above 350 K the equilibrium highly favors decomposition back to the reactants (60, 61). In addition, electronic structure calculations have shown that benzene + OH reactions do not lead primarily to oxygen addition to benzene at high temperatures. The potential energy barrier for the benzene + OH → C6H5OH + H reaction proceeding via the C6H6OH prereaction complex is 9–82 times greater than the reaction barrier for benzene + OH → phenyl + H2O (60, 62). Hence, H-abstraction reactions by OH are much more favorable than direct OH addition to the benzene ring (60, 62, 63). The substitution reaction (PAH + OH → PAH-OH + H) accounts for less than 10% of the H-atom abstraction rate, and the formation of C6H6OH is negligible (62). Electronic structure calculations have also shown that the products following H abstraction (phenyl + H2O) are 24–40 kJ/mol more stable than the products following oxygenation (phenol + H) (60, 62).

In addition, the mole fraction of H is roughly 8–15 times higher than the OH mole fraction at the flame height (DFFO) where most of the oxygenation occurs, that is, between 0.5 mm and 2 mm, as shown in Fig. S3. Hence, PAHs are predominantly activated by H abstraction. An OH molecule can then add to the radical site on the aromatic ring. This OH group then readily and rapidly decomposes to an oxyradical + H at high temperatures (64). This rapid decomposition leads to an abundance of active oxygen sites attached to aromatic rings, which is why the frequency of acetylene addition to the oxyradical is so high in our SNAPS simulations. At lower temperatures, however, the alcohols are more stable in the SNAPS simulations. We established a step-by-step elementary approach to the oxygenation pathways, but, because the rates are very fast for some reactions (e.g., H elimination from oxygen at high temperatures), the results are akin to a “global” reaction, which proceeds as PAH + OH → PAH-O + 2H (or H2), where the carbon ring maintains its aromaticity.

Important Enol Intermediates Predicted by SNAPS.

The SNAPS modeling results showed the importance of enols in the growth mechanisms of furans. We found that, in addition to the isomerization reaction, removal of enols may proceed via reactions with flame radicals, particularly H and OH, producing oxyradicals. Production of large enols can proceed through radical attacks on the carbon–carbon double bond of enols, followed by alkyl radical addition. Masses identified by SNAPS to contain enol contributions at small DFFOs in the C2H4/O2/Ar flame are 94 u (phenol), 108 u (C7H8O), 118 u (C8H6O), 168 u (C12H8O), and 192 u (C14H8O). The structures are displayed in Fig. S1.

Information Related to Fig. 3.

The SNAPS simulations revealed that, between a DFFO of ∼3.5 mm and a DFFO of ∼8.2 mm, the structures in Fig. 3 A and B evolved into larger furan compounds; the structure in Fig. 3D either evolved into furan compounds (70% probability) or became deoxygenated (30%), and the structure in Fig. 3E maintained its dibenzofuran substructure while following traditional HACA growth pathways on the other sites of the molecule. These sites are where chlorine atoms may add to form PCDFs. Hence, the structure in Fig. 3E is a likely PCDF precursor.

Although ethers and ketenes are common in Fig. 3, the majority of the oxygenated compounds predicted by SNAPS at a DFFO of ∼3.5 mm were furans. The fraction of ethers, and furans in particular, among the oxygenated structures increased with increasing DFFO in the simulations. This result was validated by the XPS measurements.

Determination of Precise Masses and Atomic Constituents.

The aerosol mass spectrometer yielded different flight times for species of different masses. The flight times were converted into molecular masses through precise calibrations. The resulting masses were sufficiently precise to distinguish between species of different atomic compositions of carbon, hydrogen, and oxygen based on the differences in mass between these elements. The mass peaks were bell-shaped and we fitted Gaussians to the mass peaks to optimize the accuracy of the actual peak locations. We then compared the measured masses with the masses of potential atomic combinations of carbon, hydrogen, and oxygen. A maximum of four oxygen atoms was allowed in a structure because of the high C/O ratios required for soot formation. No constraints were placed on the number of C and H atoms allowed. The atomic masses (65) used were 12C, 12.000000 u; 1H, 1.007825 u; and 16O, 15.994915 u. Some of the results are summarized in Table S1 along with the results from the SNAPS simulations. The mass peak at 141.0705 u included in Table S1 is not influenced by 13C interference, as can be seen in the mass spectra in Fig. 4 (the signals at 139 and 140 u are weak compared with the signal at 141 u).

Nominal Masses of Oxygen-Containing Species Determined from the Aerosol Mass Spectra.

Table S2 shows the nominal masses of peaks identified as containing signal from oxygenated structures in the AMS investigation. Ions of different masses appear at the detector with decreasing time separations for increasing molecular masses because of the quadratic relation between molecular mass and molecular-ion-flight time (m = a2t2 + a1t + a0) in the aerosol mass spectrometer. As a consequence, the analysis of the precise mass-peak locations becomes less reliable for increasing molecular masses, and, beyond a mass of ∼250 u, we could not completely rely on the measured masses to determine the atomic compositions. Hence, for species heavier than ∼250 u, we also considered the spatial profiles of the species when we determined whether they contained oxygen or not. Profiles related to heavy species, however, tend to be more difficult to interpret than profiles of lighter molecules, indicating a larger degree of coexistence of pure hydrocarbons and oxygenated hydrocarbons for heavy species. Therefore, we limited our analysis to species weighing less than 350 u.

Acknowledgments

We thank Prof. Barbara Finlayson-Pitts for providing valuable input on an earlier version of our manuscript and to Dr. Paolo Elvati for insightful discussions. This work was funded by the US Department of Energy (DOE) Office of Basic Energy Sciences (BES), Single Investigator Small Group Research Grant DE-SC0002619 (to A.V., T.D., and K.O.J.), and an Alexander von Humboldt Foundation Feodor Lynen Fellowship (to D.M.P.-V.). Experimental expenses, including burner design and construction, were funded under DOE BES, the Division of Chemical Sciences, Geosciences, and Biosciences (M.F.C., P.E.S., and H.A.M.). Aerosol mass spectrometry measurements were performed at the Advanced Light Source (ALS) at Lawrence Berkeley National Laboratory. The ALS, N.K.R.-H., and K.R.W. were supported by the Director, DOE BES, under Contract DE-AC02-05CH11231. Experimental preparations and XPS measurements were performed at Sandia National Laboratories, which is a multiprogram laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Company, for the DOE’s National Nuclear Security Administration under Contract DE-AC04-94-AL85000.

Footnotes

  • ↵1To whom correspondence may be addressed. Email: avioli{at}umich.edu or hamiche{at}sandia.gov.
  • Author contributions: A.V. and H.A.M. designed research; K.O.J., T.D., M.M., F.E.G., M.F.C., P.E.S., D.M.P.-V., N.K.R.-H., K.R.W., A.V., and H.A.M. performed research; K.O.J., T.D., M.M., and F.E.G. analyzed data; and K.O.J., T.D., A.V., and H.A.M. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1604772113/-/DCSupplemental.

Freely available online through the PNAS open access option.

View Abstract

References

  1. ↵
    1. Peterson LA
    (2006) Electrophilic intermediates produced by bioactivation of furan. Drug Metab Rev 38(4):615–626.
    .
    OpenUrlCrossRefPubMed
  2. ↵
    1. Ravindranath V,
    2. Burka LT,
    3. Boyd MR
    (1984) Reactive metabolites from the bioactivation of toxic methylfurans. Science 224(4651):884–886.
    .
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Monien BH,
    2. Herrmann K,
    3. Florian S,
    4. Glatt H
    (2011) Metabolic activation of furfuryl alcohol: Formation of 2-methylfuranyl DNA adducts in Salmonella typhimurium strains expressing human sulfotransferase 1A1 and in FVB/N mice. Carcinogenesis 32(10):1533–1539.
    .
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. National Institutes of Health
    (1993) Toxicology and carcinogenesis studies of furan (CAS No. 110-00-9) in F344 rats and B6C3F1 mice (gavage studies). Natl Toxicol Program Tech Rep Ser 402:1–286.
    .
    OpenUrlPubMed
  5. ↵
    1. World Health Organization
    (2014) IARC monographs on the evaluation of carcinogenic risks to humans, Internal report 14/002 (World Health Organization, Lyon, France).
    .
  6. ↵
    1. Lavric ED,
    2. Konnov AA,
    3. De Ruyck J
    (2004) Dioxin levels in wood combustion—A review. Biomass Bioenergy 26(2):115–145.
    .
    OpenUrlCrossRef
  7. ↵
    1. Tame NW,
    2. Dlugogorski BZ,
    3. Kennedy EM
    (2007) Formation of dioxins and furans during combustion of treated wood. Prog Energy Combust Sci 33(4):384–408.
    .
    OpenUrlCrossRef
  8. ↵
    1. Czuczwa JM,
    2. McVeety BD,
    3. Hites RA
    (1984) Polychlorinated dibenzo-p-dioxins and dibenzofurans in sediments from Siskiwit Lake, Isle Royale. Science 226(4674):568–569.
    .
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Meharg AA,
    2. Killham K
    (2003) Environment: A pre-industrial source of dioxins and furans. Nature 421(6926):909–910.
    .
    OpenUrlCrossRefPubMed
  10. ↵
    1. Sanders EB,
    2. Goldsmith AI,
    3. Seeman JI
    (2003) A model that distinguishes the pyrolysis of D-glucose, D-fructose, and sucrose from that of cellulose. Application to the understanding of cigarette smoke formation. J Anal Appl Pyrolysis 66(1):29–50.
    .
    OpenUrlCrossRef
  11. ↵
    1. Schubert J,
    2. Bewersdorff J,
    3. Luch A,
    4. Schulz TG
    (2012) Waterpipe smoke: A considerable source of human exposure against furanic compounds. Anal Chim Acta 709:105–112.
    .
    OpenUrlCrossRefPubMed
  12. ↵
    1. Valberg PA,
    2. Drivas PJ,
    3. McCarthy S,
    4. Watson AY
    (1996) Evaluating the health impacts of incinerator emissions. J Hazard Mater 47(1):205–227.
    .
    OpenUrlCrossRef
  13. ↵
    1. Man M,
    2. Naidu R,
    3. Wong MH
    (2013) Persistent toxic substances released from uncontrolled e-waste recycling and actions for the future. Sci Total Environ 463-464:1133–1137.
    .
    OpenUrl
  14. ↵
    1. Grant K, et al.
    (2013) Health consequences of exposure to e-waste: A systematic review. Lancet Glob Health 1(6):e350–e361.
    .
    OpenUrlCrossRef
  15. ↵
    1. Tassi F,
    2. Montegrossi G,
    3. Capecchiacci F,
    4. Vaselli O
    (2010) Origin and distribution of thiophenes and furans in gas discharges from active volcanoes and geothermal systems. Int J Mol Sci 11(4):1434–1457.
    .
    OpenUrlCrossRefPubMed
  16. ↵
    1. Meharg AA,
    2. Osborn D
    (1995) Dioxins released from chemical accidents. Nature 375(6530):353–354.
    .
    OpenUrlCrossRefPubMed
  17. ↵
    1. Wilhelm J,
    2. Stieglitz L,
    3. Dinjus E,
    4. Will R
    (2001) Mechanistic studies on the role of PAHs and related compounds in PCDD/F formation on model fly ashes. Chemosphere 42(5-7):797–802.
    .
    OpenUrlCrossRefPubMed
  18. ↵
    1. Iino F,
    2. Imagawa T,
    3. Takeuchi M,
    4. Sadakata M
    (1999) De novo synthesis mechanism of polychlorinated dibenzofurans from polycyclic aromatic hydrocarbons and the characteristic isomers of polychlorinated naphthalenes. Environ Sci Technol 33(7):1038–1043.
    .
    OpenUrlCrossRef
  19. ↵
    1. Seinfeld JH,
    2. Pandis SN
    (1998) Atmospheric Chemistry and Physics: From Air Pollution to Climate Change (Wiley, New York).
    .
  20. ↵
    1. Atkinson R,
    2. Arey J,
    3. Aschmann SM
    (2008) Atmospheric chemistry of alkanes: Review and recent developments. Atmos Environ 42(23):5859–5871.
    .
    OpenUrlCrossRef
  21. ↵
    1. Schilling Fahnestock KA, et al.
    (2015) Secondary organic aerosol composition from C12 alkanes. J Phys Chem A 119(19):4281–4297.
    .
    OpenUrlCrossRefPubMed
  22. ↵
    1. Northcross AL,
    2. Hammond SK,
    3. Canuz E,
    4. Smith KR
    (2012) Dioxin inhalation doses from wood combustion in indoor cookfires. Atmos Environ 49:415–418.
    .
    OpenUrlCrossRef
  23. ↵
    1. Pieters R,
    2. Focant J-F
    (2014) Dioxin, furan and PCB serum levels in a South African Tswana population: Comparing the polluting effects of using different cooking and heating fuels. Environ Int 66:71–78.
    .
    OpenUrlCrossRefPubMed
  24. ↵
    1. Cain J,
    2. Laskin A,
    3. Kholghy MR,
    4. Thomson MJ,
    5. Wang H
    (2014) Molecular characterization of organic content of soot along the centerline of a coflow diffusion flame. Phys Chem Chem Phys 16(47):25862–25875.
    .
    OpenUrlCrossRefPubMed
  25. ↵
    1. McKinnon JT,
    2. Meyer E,
    3. Howard JB
    (1996) Infrared analysis of flame generated PAH samples. Combust Flame 105:161–166.
    .
    OpenUrlCrossRef
  26. ↵
    1. Santamaria A, et al.
    (2006) FT-IR and 1H NMR characterization of the products of an ethylene inverse diffusion flame. Combust Flame 146:52–62.
    .
    OpenUrlCrossRef
  27. ↵
    1. di Stasio S,
    2. Braun A
    (2006) Comparative NEXAFS study on soot obtained from an ethylene/air flame, a Diesel engine, and graphite. Energy Fuels 20:187–194.
    .
    OpenUrlCrossRef
  28. ↵
    1. Commodo M,
    2. Falco GD,
    3. Larciprete R,
    4. D’Anna A,
    5. Minutolo P
    (2016) On the hydrophilic/hydrophobic character of carbonaceous nanoparticles formed in laminar premixed flames. Exp Therm Fluid Sci 73:56–63.
    .
    OpenUrlCrossRef
  29. ↵
    1. Gaddam CK,
    2. Vander Wal RL
    (2013) Physical and chemical characterization of SIDI engine particulates. Combust Flame 160:2517–2528.
    .
    OpenUrlCrossRef
  30. ↵
    1. Müller J-O,
    2. Su DS,
    3. Wild U,
    4. Schlögl R
    (2007) Bulk and surface structural investigations of diesel engine soot and carbon black. Phys Chem Chem Phys 9(30):4018–4025.
    .
    OpenUrlCrossRefPubMed
  31. ↵
    1. Cain JP,
    2. Gassman PL,
    3. Wang H,
    4. Laskin A
    (2010) Micro-FTIR study of soot chemical composition-evidence of aliphatic hydrocarbons on nascent soot surfaces. Phys Chem Chem Phys 12(20):5206–5218.
    .
    OpenUrlCrossRefPubMed
  32. ↵
    1. Mustafi N,
    2. Raine RR,
    3. James B
    (2010) Characterization of exhaust particulates from a dual fuel engine by TGA, XPS, and Raman techniques. Aerosol Sci Technol 44(11):954–963.
    .
    OpenUrlCrossRef
  33. ↵
    1. Bond TC, et al.
    (2013) Bounding the role of black carbon in the climate system: A scientific assessment. J Geophys Res 118(11):5380–5552.
    .
    OpenUrl
  34. ↵
    1. Zhang R, et al.
    (2008) Variability in morphology, hygroscopicity, and optical properties of soot aerosols during atmospheric processing. Proc Natl Acad Sci USA 105(30):10291–10296.
    .
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Kaufman YJ,
    2. Koren I,
    3. Remer LA,
    4. Rosenfeld D,
    5. Rudich Y
    (2005) The effect of smoke, dust, and pollution aerosol on shallow cloud development over the Atlantic Ocean. Proc Natl Acad Sci USA 102(32):11207–11212.
    .
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Ramanathan V, et al.
    (2005) Atmospheric brown clouds: Impacts on South Asian climate and hydrological cycle. Proc Natl Acad Sci USA 102(15):5326–5333.
    .
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Lai JY,
    2. Elvati P,
    3. Violi A
    (2014) Stochastic atomistic simulation of polycyclic aromatic hydrocarbon growth in combustion. Phys Chem Chem Phys 16(17):7969–7979.
    .
    OpenUrlCrossRefPubMed
  38. ↵
    1. Taatjes CA, et al.
    (2005) Enols are common intermediates in hydrocarbon oxidation. Science 308(5730):1887–1889.
    .
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Öktem B,
    2. Tolocka MP,
    3. Zhao B,
    4. Wang H,
    5. Johnston MV
    (2005) Chemical species associated with the early stage of soot growth in a laminar premixed ethylene–oxygen–argon flame. Combust Flame 142(4):364–373.
    .
    OpenUrlCrossRef
  40. ↵
    1. Wartel M,
    2. Pauwels J-F,
    3. Desgroux P,
    4. Mercier X
    (2011) Pyrene measurements in sooting low pressure methane flames by jet-cooled laser-induced fluorescence. J Phys Chem A 115(49):14153–14162.
    .
    OpenUrlCrossRefPubMed
  41. ↵
    1. McFeely F, et al.
    (1974) X-ray photoemission studies of diamond, graphite, and glassy carbon valence bands. Phys Rev B 9(12):5268–5278.
    .
    OpenUrlCrossRef
  42. ↵
    1. Han C,
    2. Liu Y,
    3. Ma J,
    4. He H
    (2012) Key role of organic carbon in the sunlight-enhanced atmospheric aging of soot by O2. Proc Natl Acad Sci USA 109(52):21250–21255.
    .
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Gudipati MS,
    2. Yang R
    (2012) In-situ probing of radiation-induced processing of organics in astrophysical ice analogs—Novel laser desorption laser ionization time-of-flight mass spectroscopic studies. Astrophys J 756(1):L24.
    .
    OpenUrlCrossRef
  44. ↵
    1. Park S,
    2. Vohs JM,
    3. Gorte RJ
    (2000) Direct oxidation of hydrocarbons in a solid-oxide fuel cell. Nature 404(6775):265–267.
    .
    OpenUrlCrossRefPubMed
  45. ↵
    1. Johansson K, et al.
    (2015) Soot precursor formation and limitations of the stabilomer grid. Proc Combust Inst 35(2):1819–1826.
    .
    OpenUrlCrossRef
  46. ↵
    1. Kee R, et al.
    (2006) CHEMKIN Release 4.1 (Reaction Design, San Diego).
    .
  47. ↵
    1. Appel J,
    2. Bockhorn H,
    3. Frenklach M
    (2000) Kinetic modeling of soot formation with detailed chemistry and physics: Laminar premixed flames of C2 hydrocarbons. Combust Flame 121(1):122–136.
    .
    OpenUrlCrossRef
  48. ↵
    1. Miller JA,
    2. Melius CF
    (1992) Kinetic and thermodynamic issues in the formation of aromatic compounds in flames of aliphatic fuels. Combust Flame 91(1):21–39.
    .
    OpenUrl
  49. ↵
    1. Richter H,
    2. Benish TG,
    3. Mazyar OA,
    4. Green WH,
    5. Howard JB
    (2000) Formation of polycyclic aromatic hydrocarbons and their radicals in a nearly sooting premixed benzene flame. Proc Combust Inst 28(2):2609–2618.
    .
    OpenUrlCrossRef
  50. ↵
    1. Raj A,
    2. Prada IDC,
    3. Amer AA,
    4. Chung SH
    (2012) A reaction mechanism for gasoline surrogate fuels for large polycyclic aromatic hydrocarbons. Combust Flame 159(2):500–515.
    .
    OpenUrlCrossRef
  51. ↵
    1. Richter H,
    2. Howard JB
    (2000) Formation of polycyclic aromatic hydrocarbons and their growth to soot: A review of chemical reaction pathways. Pror Energy Combust Sci 26:565–608.
    .
    OpenUrlCrossRef
  52. ↵
    1. Pietrucci F,
    2. Andreoni W
    (2011) Graph theory meets ab initio molecular dynamics: Atomic structures and transformations at the nanoscale. Phys Rev Lett 107(8):085504.
    .
    OpenUrlCrossRefPubMed
  53. ↵
    1. Montgomery JA Jr,
    2. Frisch MJ,
    3. Ochterski JW,
    4. Petersson GA
    (1999) A complete basis set model chemistry. VI. Use of density functional geometries and frequencies. J Chem Phys 110(6):2822–2827.
    .
    OpenUrlCrossRef
  54. ↵
    1. Frisch MJ, et al.
    (2009) Gaussian 09, Revision A.02 (Gaussian, Inc., Wallingford, CT).
    .
  55. ↵
    1. Zhang X, et al.
    (2004) Numerical characterization of particle beam collimation: Part II integrated aerodynamic-lens–nozzle system. Aerosol Sci Technol 38(6):619–638.
    .
    OpenUrlCrossRef
  56. ↵
    1. Headrick JM,
    2. Schrader PE,
    3. Michelsen HA
    (2013) Radial-profile and divergence measurements of combustion-generated soot focused by an aerodynamic-lens system. J Aerosol Sci 58:158–170.
    .
    OpenUrlCrossRef
  57. ↵
    1. Moshammer K,
    2. Lucassen A,
    3. Togbé C,
    4. Kohse-Höinghaus K,
    5. Hansen N
    (2015) Formation of oxygenated and hydrocarbon intermediates in premixed combustion of 2-methylfuran. Z Phys Chem 229(4):507–528.
    .
    OpenUrl
  58. ↵
    1. Hansen N,
    2. Cool TA,
    3. Westmoreland PR,
    4. Kohse-Höinghaus K
    (2009) Recent contributions of flame-sampling molecular-beam mass spectrometry to a fundamental understanding of combustion chemistry. Pror Energy Combust Sci 35(2):168–191.
    .
    OpenUrlCrossRef
  59. ↵
    1. Guner V, et al.
    (2003) A standard set of pericyclic reactions of hydrocarbons for the benchmarking of computational methods: The performance of ab initio, density functional, CASSCF, CASPT2, and CBS-QB3 methods for the prediction of activation barriers, reaction energetics, and transition state geometries. J Phys Chem A 107(51):11445–11459.
    .
    OpenUrlCrossRef
  60. ↵
    1. Tokmakov I,
    2. Lin M
    (2002) Kinetics and mechanism of the OH+ C6H6 reaction: A detailed analysis with first-principles calculations. J Phys Chem A 106(46):11309–11326.
    .
    OpenUrlCrossRef
  61. ↵
    1. Bohn B,
    2. Zetzsch C
    (1999) Gas-phase reaction of the OH–benzene adduct with O2: Reversibility and secondary formation of HO2. Phys Chem Chem Phys 1(22):5097–5107.
    .
    OpenUrlCrossRef
  62. ↵
    1. Seta T,
    2. Nakajima M,
    3. Miyoshi A
    (2006) High-temperature reactions of OH radicals with benzene and toluene. J Phys Chem A 110(15):5081–5090.
    .
    OpenUrlCrossRefPubMed
  63. ↵
    1. Blanquart G,
    2. Pepiot-Desjardins P,
    3. Pitsch H
    (2009) Chemical mechanism for high temperature combustion of engine relevant fuels with emphasis on soot precursors. Combust Flame 156(3):588–607.
    .
    OpenUrlCrossRef
  64. ↵
    1. Xu ZF,
    2. Lin MC
    (2006) Ab initio kinetics for the unimolecular reaction C6H5OH–> CO + C5H6. J Phys Chem A 110(4):1672–1677.
    .
    OpenUrlCrossRefPubMed
  65. ↵
    1. Audi G,
    2. Wapstra A
    (1995) The 1995 update to the atomic mass evaluation. Nucl Phys A 595(4):409–480.
    .
    OpenUrlCrossRef
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Formation and emission of large furans and oxygenated hydrocarbons from flames
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Large furan and oxygenate formation in flames
K. Olof Johansson, Tyler Dillstrom, Matteo Monti, Farid El Gabaly, Matthew F. Campbell, Paul E. Schrader, Denisia M. Popolan-Vaida, Nicole K. Richards-Henderson, Kevin R. Wilson, Angela Violi, Hope A. Michelsen
Proceedings of the National Academy of Sciences Jul 2016, 113 (30) 8374-8379; DOI: 10.1073/pnas.1604772113

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Large furan and oxygenate formation in flames
K. Olof Johansson, Tyler Dillstrom, Matteo Monti, Farid El Gabaly, Matthew F. Campbell, Paul E. Schrader, Denisia M. Popolan-Vaida, Nicole K. Richards-Henderson, Kevin R. Wilson, Angela Violi, Hope A. Michelsen
Proceedings of the National Academy of Sciences Jul 2016, 113 (30) 8374-8379; DOI: 10.1073/pnas.1604772113
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 113 (30)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Physical Sciences
  • Chemistry

Jump to section

  • Article
    • Abstract
    • Results and Discussion
    • Summary and Conclusions
    • Methods
    • Details About the Approaches
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Abstract depiction of a guitar and musical note
Science & Culture: At the nexus of music and medicine, some see disease treatments
Although the evidence is still limited, a growing body of research suggests music may have beneficial effects for diseases such as Parkinson’s.
Image credit: Shutterstock/agsandrew.
Scientist looking at an electronic tablet
Opinion: Standardizing gene product nomenclature—a call to action
Biomedical communities and journals need to standardize nomenclature of gene products to enhance accuracy in scientific and public communication.
Image credit: Shutterstock/greenbutterfly.
One red and one yellow modeled protein structures
Journal Club: Study reveals evolutionary origins of fold-switching protein
Shapeshifting designs could have wide-ranging pharmaceutical and biomedical applications in coming years.
Image credit: Acacia Dishman/Medical College of Wisconsin.
White and blue bird
Hazards of ozone pollution to birds
Amanda Rodewald, Ivan Rudik, and Catherine Kling talk about the hazards of ozone pollution to birds.
Listen
Past PodcastsSubscribe
Goats standing in a pin
Transplantation of sperm-producing stem cells
CRISPR-Cas9 gene editing can improve the effectiveness of spermatogonial stem cell transplantation in mice and livestock, a study finds.
Image credit: Jon M. Oatley.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Latest Articles
  • Archive

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490