Computational investigation of the pyrolysis product selectivity for α-hydroxy phenethyl phenyl ether and phenethyl phenyl ether: analysis of substituent effects and reactant conformer selection.

Using computational methods, we determine product selectivities for the pyrolysis of two model compounds for the ß-O-4 linkage in lignin: phenethyl phenyl ether (PPE) and α-hydroxy phenethyl phenyl ether (α-hydroxy PPE). We investigate the dependence of the product selectivities on the number of reactant conformers included. Utilizing density functional theory in combination with transition state theory, we obtain rate constants for hydrogen abstraction reactions by the key chain-carrying radicals, which determine the product selectivity within a steady-state kinetic model. The inclusion of the energetically second lowest reactant conformer of PPE and α-hydroxy PPE has a large effect on the product selectivity. The final product selectivity computed for PPE agrees well with experiment. We find that the α-hydroxy substituent affects energetic as well as entropic contributions to the rate constant differently for alternative pathways of hydrogen abstraction and, thereby, significantly alters product distributions.

Role of carbon-carbon phenyl migration in the pyrolysis mechanism of ß-O-4 lignin model compounds: phenethyl phenyl ether and α-hydroxy phenethyl phenyl ether.

We investigate phenyl shift and subsequent ß-scission reactions for PhCHXCH·OPh [X = H, OH], which are part of the pyrolysis mechanism of phenethyl phenyl ether (PPE) and α-hydroxy PPE. PPE and its derivatives are model compounds for the most common linkage in lignin, the ß-O-4 linkage. We use density functional theory to locate transition states and equilibrium structures and kinetic Monte Carlo in combination with transition-state theory for kinetic simulations. Oxygen-carbon and carbon-carbon phenyl shift reactions proceed through cyclic intermediates with similar barriers. However, while subsequent ß-scission of the oxygen-carbon shift products proceeds with virtually no barrier, the activation energy for ß-scission of the carbon-carbon shift products exceeds 15 kcal/mol. We found that about 15% of ß-radical conversion can be attributed to carbon-carbon shift for PPE and α-hydroxy PPE at 618 K. Whereas the oxygen-carbon shift reaction has been established as an integral part of the pyrolysis mechanism of PPE and its derivatives, participation of the carbon-carbon shift reaction has not been shown previously.

Computational study of bond dissociation enthalpies for substituted ß-O-4 lignin model compounds.

The biopolymer lignin is a potential source of valuable chemicals. Phenethyl phenyl ether (PPE) is representative of the dominant ß-O-4 ether linkage. DFT is used to calculate the Boltzmann-weighted carbon-oxygen and carbon-carbon bond dissociation enthalpies (BDEs) of substituted PPE. These values are important for understanding lignin decomposition. Exclusion of all conformers that have distributions of less than 5% at 298 K impacts the BDE by less than 1 kcal mol(-1). We find that aliphatic hydroxyl/methylhydroxyl substituents introduce only small changes to the BDEs (0-3 kcal mol(-1)). Substitution on the phenyl ring at the ortho position substantially lowers the C-O BDE, except in combination with the hydroxyl/methylhydroxyl substituents, for which the effect of methoxy substitution is reduced by hydrogen bonding. Hydrogen bonding between the aliphatic substituents and the ether oxygen in the PPE derivatives has a significant influence on the BDE. CCSD(T)-calculated BDEs and hydrogen-bond strengths of ortho-substituted anisoles, when compared with M06-2X values, confirm that the latter method is sufficient to describe the molecules studied and provide an important benchmark for lignin model compounds.

Pyrolysis of phenethyl phenyl ether tethered in mesoporous silica. Effects of confinement and surface spacer molecules on product selectivity.

There has been expanding interest in exploring porous metal oxides as a confining environment for organic molecules resulting in altered chemical and physical properties including chemical transformations. In this paper, we examine the pyrolysis behavior of phenethyl phenyl ether (PPE) confined in mesoporous silica by covalent tethers to the pore walls as a function of tether density and the presence of cotethered surface spacer molecules of varying structure (biphenyl, naphthyl, octyl, and hexadecyl). The PPE pyrolysis product selectivity, which is determined by two competitive free-radical pathways cycling through the two aliphatic radical intermediates (PhCH·CH(2)OPh and PhCH(2)CH·OPh), is shown to be significantly different from that measured in the liquid phase as well as for PPE tethered to the exterior surface of nonporous silica nanoparticles. Tailoring the pore surface with spacer molecules further alters the selectivity such that the PPE reaction channel involving a molecular rearrangement (O-C phenyl shift in PhCH(2)CH·OPh), which accounts for 25% of the products in the liquid phase, can be virtually eliminated under pore confinement conditions. The origin of this change in selectivity is discussed in the context of steric constraints on the rearrangement path inside the pores, surface and pore confinement effects, pore surface curvature, and hydrogen bonding of PPE with residual surface silanols supplemented by nitrogen physisorption data and molecular dynamics simulations.

Kinetic analysis of the phenyl-shift reaction in ß-O-4 lignin model compounds: a computational study.

The phenyl-shift reaction for the ß-radical of phenethyl phenyl ether (PhCH(2)CHOPh, ß-PPE) is an integral step in the pyrolysis of PPE, which is a model compound for the ß-O-4 linkage in lignin. We investigated the influence of natural occurring substituents (hydroxy, methoxy) on the reaction rate by calculating relative rate constants using density functional theory in combination with transition state theory, including anharmonic correction for low-frequency modes. The phenyl-shift reaction proceeds through an oxaspiro[2.5]octadienyl radical intermediate and the overall rate constants were computed invoking the steady-state approximation (its validity was confirmed). Substituents on the phenethyl ring have only little influence on the rate constants. If a methoxy substituent is located in the para position of the phenyl ring adjacent to the ether oxygen, the energies of the intermediate and second transition state are lowered, but the overall rate constant is not significantly altered. This is a consequence of the dominating first transition from reactant to intermediate in the overall rate constant. In contrast, o- and di-o-methoxy substituents significantly accelerate the phenyl-migration rate compared to ß-PPE.

Characterization of biochars produced from cornstovers for soil amendment.

Through cation exchange capacity assay, nitrogen adsorption-desorption surface area measurements, scanning electron microscopic imaging, infrared spectra and elemental analyses, we characterized biochar materials produced from cornstover under two different pyrolysis conditions, fast pyrolysis at 450 °C and gasification at 700 °C. Our experimental results showed that the cation exchange capacity (CEC) of the fast-pyrolytic char is about twice as high as that of the gasification char as well as that of a standard soil sample. The CEC values correlate well with the increase in the ratios of the oxygen atoms to the carbon atoms (O:C ratios) in the biochar materials. The higher O:C ratio was consistent with the presence of more hydroxyl, carboxylate, and carbonyl groups in the fast pyrolysis char. These results show how control of biomass pyrolysis conditions can improve biochar properties for soil amendment and carbon sequestration. Since the CEC of the fast-pyrolytic cornstover char can be about double that of a standard soil sample, this type of biochar products would be suitable for improvement of soil properties such as CEC, and at the same time, can serve as a carbon sequestration agent.

Computational study of bond dissociation enthalpies for lignin model compounds. Substituent effects in phenethyl phenyl ethers.

Lignin is an abundant natural resource that is a potential source of valuable chemicals. Improved understanding of the pyrolysis of lignin occurs through the study of model compounds for which phenethyl phenyl ether (PhCH(2)CH(2)OPh, PPE) is the simplest example representing the dominant beta-O-4 ether linkage. The initial step in the thermal decomposition of PPE is the homolytic cleavage of the oxygen-carbon bond. The rate of this key step will depend on the bond dissociation enthalpy, which in turn will depend on the nature and location of relevant substituents. We used modern density functional methods to calculate the oxygen-carbon bond dissociation enthalpies for PPE and several oxygen-substituted derivatives. Since carbon-carbon bond cleavage in PPE could be a competitive initial reaction under high-temperature pyrolysis conditions, we also calculated substituent effects on these bond dissociation enthalpies. We found that the oxygen-carbon bond dissociation enthalpy is substantially lowered by oxygen substituents situated at the phenyl ring adjacent to the ether oxygen. On the other hand, the carbon-carbon bond dissociation enthalpy shows little variation with different substitution patterns on either phenyl ring.

Computational prediction of alpha/beta selectivities in the pyrolysis of oxygen-substituted phenethyl phenyl ethers.

Phenethyl phenyl ether (PPE; PhCH 2CH 2OPh) is the simplest model for the most common beta-O-4 linkage in lignin. Previously, we developed a computational scheme to calculate the alpha/beta product selectivity in the pyrolysis of PPE by systematically exploiting error cancellation in the computation of relative rate constants. The alpha/beta selectivity is defined as the selectivity between the competitive hydrogen abstraction reaction paths on the alpha- and beta-carbons of PPE. We use density functional theory and employ transition state theory where we include diagonal anharmonic correction in the vibrational partition functions for low frequency modes for which a semiclassical expression is used. In this work we investigate the effect of oxygen substituents (hydroxy, methoxy) in the para position on the phenethyl ring of PPE on the alpha/beta selectivities. The total alpha/beta selectivity increases when substituents are introduced and is larger for the methoxy than the hydroxy substituent. The strongest effect of the substituents is observed for the alpha-pathway of the hydrogen abstraction by the phenoxyl chain carrying radical for which the rate increases. For the beta pathway and the abstraction by the R-benzyl radical (R = OH,OCH 3) the rate decreases with the introduction of the substituents. These findings are compared with results from recent experimental studies.

Kinetic analysis of the pyrolysis of phenethyl phenyl ether: computational prediction of alpha/beta-selectivities.

We calculated an overall alpha/beta-selectivity for the pyrolysis of phenethyl phenyl ether as a composite of the alpha/beta-selectivities in the hydrogen abstraction reactions by the phenoxyl and by the benzyl radical that is in excellent agreement with experiment. The difference between the individual selectivities for these radicals is explained by analyzing the electronic structure of the transition states. Spin delocalization of the single electron favors the alpha-pathways. An opposing effect occurs for polarized transition states, such as the transition states for the hydrogen abstraction by the electrophilic phenoxyl radical, where the adjacent ether oxygen in phenethyl phenyl ether stabilizes the beta-transition states. These results indicate that theory will be able to provide excellent predictions of alpha/beta-product selectivities for more complicated lignin model compounds bearing multiple substituents. We have developed a scheme to predict alpha/beta-product selectivities in the pyrolysis of model compounds for the beta-ether linkage in lignin. The approach is based on computation of the relative rate constant, which profits from error cancellation in the individual rate constants. The Arrhenius prefactors depend strongly on the description of the low-frequency modes for which anharmonic contributions are important. We use density functional theory in combination with transition-state theory in this analysis. Diagonal anharmonic effects for individual low-frequency modes are included by employing a second-order Wigner-Kirkwood expansion in a semiclassical expression for the vibrational partition function. The composite alpha/beta-product selectivity is obtained by applying quasi-steady-state kinetic analysis for the intermediate radicals.

Confinement effects on product selectivity in the pyrolysis of phenethyl phenyl ether in mesoporous silica.

Pyrolysis of phenethyl phenyl ether confined in mesoporous silicas by covalent grafting results in significantly increased product selectivity compared with fluid phases.

Pyrolysis of mesoporous silica-immobilized 1,3-diphenylpropane. Impact of pore confinement and size.

Mesoporous silicas such as SBA-15 and MCM-41 are being actively investigated for potential applications in catalysis, separations, and synthesis of nanostructured materials. A new method for functionalizing these mesoporous silicas with aromatic phenols is described. The resulting novel hybrid materials possess silyl aryl ether linkages to the silica surface that are thermally stable to ca. 550 degrees C, but can be easily cleaved at room temperature with aqueous base for quantitative recovery of the organic moieties. The materials have been characterized by nitrogen physisorption, FTIR, NMR, and quantitative analysis of surface coverages. The maximum densities of 1,3-diphenylpropane (DPP) molecules that could be grafted to the surface were less than those measured on a nonporous, fumed silica (Cabosil) and were also found to decrease as a function of decreasing pore size (5.6-1.7 nm). This is a consequence of steric congestion in the pores that is magnified at the smaller pore sizes, consistent with parallel studies conducted using a conventional silylating reagent, 1,1,3,3-tetramethyldisilazane. Pyrolysis of the silica-immobilized DPP revealed that pore confinement leads to enhanced rates and altered product selectivity for this free-radical reaction compared with the nonporous silica, and the rates and selectivities also depended on pore size. The influence of confinement is discussed in terms of enhanced encounter frequencies for bimolecular reaction steps and pore surface curvature that alters the accessibility and resultant selectivity for hydrogen transfer steps.

Extraction of cesium ions from aqueous solutions using calix[4]arene-bis(tert-octylbenzo-crown-6) in ionic liquids.

Solvent extraction of cesium ions from aqueous solution to hydrophobic ionic liquids without the introduction of an organophilic anion in the aqueous phase was demonstrated using calix[4]arene-bis(tert-octylbenzo-crown-6) (BOBCalixC6) as an extractant. The selectivity of this extraction process toward cesium ions and the use of a sacrificial cation exchanger (NaBPh(4)) to control loss of imidazolium cation to the aqueous solutions by ion exchange have been investigated.

Pore size effects in the pyrolysis of 1,3-diphenylpropane confined in mesoporous silicas.

A new method for derivatizing mesoporous silicas, SBA-15 and MCM-41, with a substituted phenol is described, and pore confinement and surface curvature are shown to impact the reaction rate and product selectivity for the pyrolysis of surface-immobilized 1,3-diphenylpropane.

Free-radical reactions under diffusional constraints: orientation does matter in hydrogen transfer.

Pyrolysis of silica-immobilized 1,3-diphenylpropane at 375 degrees C has been examined in the presence of a series of isomeric (by point of attachment) co-attached hydroaromatic spacer molecules. Under the diffusional constraints, the pyrolysis rate is sensitive to the orientation of the spacer molecule, which must transfer hydrogen to intermediate benzylic radicals on the surface. Spacer molecules possessing a meta-orientation of the benzylic hydrogens with respect to the surface linkage are able to attain a better geometry for the hydrogen transfer on the surface resulting in faster reaction rates.

Synthetic utilization of polynitroaromatic compounds. 1. S-derivatization of 1-substituted 2,4,6-trinitrobenzenes with thiols.

Reactions of 1-R-2,4,6-trinitrobenenes (R = alkyl, protected aldehyde, aminocarbonyl, cyano groups, or isoxazole ring) with thiol salts were investigated. In most cases, these reactions gave a mixture of minor para and major ortho substitution products. Reactions of N,N-disubstituted 2,4,6-trinitrobenzamides with S-,O-, and N-nucleophiles afforded products of substitution of the p-nitro group exclusively. 1-Cyano-2,4,6-trinitrobenzene was found to be the most reactive and the least selective: all three nitro groups can be substituted using an excess of thiol salts. 2-R-4, 6-dinitrobenzamides showed no regioselectivity under similar conditions to yield 1:1 mixtures of para and ortho isomers.

Flash vacuum pyrolysis of methoxy-substituted lignin model compounds.

The flash vacuum pyrolysis (FVP) of methoxy-substituted beta-O-4 lignin model compounds has been studied at 500 degrees C to provide mechanistic insight into the primary reaction pathways that occur under conditions of fast pyrolysis. FVP of PhCH(2)CH(2)OPh (PPE), a model of the dominant beta-O-4 linkage in lignin, proceeds by C-O and C-C cleavage, in a 37:1 ratio, to produce styrene plus phenol as the dominant products and minor amounts of toluene, bibenzyl, and benzaldehyde. From the deuterium isotope effect in the FVP of PhCD(2)CH(2)OPh, it was shown that C-O cleavage occurs by homolysis and by 1,2-elimination in a ratio of 1.4:1, respectively. Methoxy substituents enhance the homolysis of the beta-O-4 linkage, relative to PPE, in o-CH(3)O-C(6)H(4)OCH(2)CH(2)Ph (o-CH(3)O-PPE) and (o-CH(3)O)(2)-C(6)H(3)OCH(2)CH(2)Ph ((o-CH(3)O)(2)-PPE) by a factor of 7.4 and 21, respectively. The methoxy-substituted phenoxy radicals undergo a complex series of reactions, which are dominated by 1,5-, 1,6-, and 1,4-intramolecular hydrogen abstraction, rearrangement, and beta-scission reactions. In the FVP of o-CH(3)O-PPE, the dominant product, salicylaldehyde, forms from the methoxyphenoxy radical by a 1,5-hydrogen shift to form 2-hydroxyphenoxymethyl radical, 1,2-phenyl shift, and beta-scission of a hydrogen atom. The 2-hydroxyphenoxymethyl radical can also cleave to form formaldehyde and phenol in which the ratio of 1, 2-phenyl shift to beta-scission is ca. 4:1. In the FVP of o-CH(3)O-PPE and (o-CH(3)O)(2)-PPE, products (ca. 20 mol %) are also formed by C-O homolysis of the methoxy group. The resulting phenoxy radicals undergo 1,5- and 1,6-hydrogen shifts in a ratio of ca. 2:1 to the aliphatic or benzylic carbon, respectively, of the phenethyl chain. In the FVP of (o-CH(3)O)(2)-PPE, o-cresol was the dominant product. It was formed by decomposition of 2-hydroxy-3-hydroxymethylbenzaldehyde and 2-hydroxybenzyl alcohol, which are formed from a complex series of reactions from the 2, 6-dimethoxyphenoxy radical. The key step in this reaction sequence was the rapid 1,5-hydrogen shift from 2-hydroxy-3-methoxybenzyloxy radical to 2-hydroxymethyl-6-methoxyphenoxy radical before beta-scission of a hydrogen atom to give the substituted benzaldehyde. The 2-hydroxybenzyl alcohols rapidly decompose under the reaction conditions to o-benzoquinone methide and pick up hydrogen from the reactor walls to form o-cresol.

Preparation of trifluoromethylpyridine libraries.

S-Alkylation followed by heterocyclization of trifluoromethyl-3-cyano-2(1H)-pyridinethiones was used for preparation of libraries of S-alkyl trifluoromethylpyridines and thieno[2,3-b]pyridines. The S-alkylation (in water--DMF mixtures) was successful for all 18 alkylating agents employed (yields typically > 50%). S-Alkyl derivatives were further converted to corresponding thieno[2,3-b]pyridines via heterocyclization in base conditions (yields > 65%). Structures of new compounds were elucidated by a combination of IR and 1H NMR spectroscopy and elemental analysis and were confirmed by means of single-crystal X-ray diffraction analysis.