Promoted cobalt metal catalysts suitable for the production of lower olefins from natural gas.

Due to the surge of natural gas production, feedstocks for chemicals shift towards lighter hydrocarbons, particularly methane. The success of a Gas-to-Chemicals process via synthesis gas (CO and H2) depends on the ability of catalysts to suppress methane and carbon dioxide formation. We designed a Co/Mn/Na/S catalyst, which gives rise to negligible Water-Gas-Shift activity and a hydrocarbon product spectrum deviating from the Anderson-Schulz-Flory distribution. At 240 °C and 1 bar, it shows a C2-C4 olefins selectivity of 54%. At 10 bar, it displays 30% and 59% selectivities towards lower olefins and fuels, respectively. The spent catalyst consists of 10 nm Co nanoparticles with hcp Co metal phase. We propose a synergistic effect of Na plus S, which act as electronic promoters on the Co surface, thus improving selectivities towards lower olefins and fuels while largely reducing methane and carbon dioxide formation.

Revisiting Hansen Solubility Parameters by Including Thermodynamics.

The Hansen solubility parameter approach is revisited by implementing the thermodynamics of dissolution and mixing. Hansen's pragmatic approach has earned its spurs in predicting solvents for polymer solutions, but for molecular solutes improvements are needed. By going into the details of entropy and enthalpy, several corrections are suggested that make the methodology thermodynamically sound without losing its ease of use. The most important corrections include accounting for the solvent molecules' size, the destruction of the solid's crystal structure, and the specificity of hydrogen-bonding interactions, as well as opportunities to predict the solubility at extrapolated temperatures. Testing the original and the improved methods on a large industrial dataset including solvent blends, fit qualities improved from 0.89 to 0.97 and the percentage of correct predictions rose from 54 % to 78 %. Full Matlab scripts are included in the Supporting Information, allowing readers to implement these improvements on their own datasets.

Predicting adsorption on metals: simple yet effective descriptors for surface catalysis.

We present a simple and efficient model for predicting the adsorption of molecules on metal surfaces. This heuristic model uses six descriptors for each metal (number of d-electrons, surface energy, first ionization potential and atomic radius, volume and mass) and three for each adsorptive (HOMO-LUMO energy gap, molecular volume and mass). Strikingly, despite its simplicity and low computational cost, this model predicts well the chemisorption of a range of adsorptives (H2, HO˙, N2, CO, NO, O2, H2O, CO2, NH3 and CH4) on a range of metals (Fe, Co, Ni, Cu, Mo, Ru, Rh, Pd, Ag, W, Ir, Pt and Au) as calculated with DFT and taken from the literature. Using only a third of the data for fitting, the rest of the data were predicted with Q(2) = 0.91-0.95 and RMSEP = 0.94-1.16 eV. Furthermore, we measured experimental adsorption data for CO, CO2, CH4, H2, N2 and O2 on Ni, Pt and Rh supported on TiO2. Using the same descriptors, we then constructed a model for this experimental data set. Once again, the model explained the data well, with R(2) = 0.95 and Q(2) = 0.86, respectively.

Oxidation of methanol by FeO2+ in water: DFT calculations in the gas phase and ab initio MD simulations in water solution.

We investigate the mechanism of methanol oxidation to formaldehyde by ironoxido ([Fe(IV)O]2+), the alleged active intermediate in the Fenton reaction. The most likely reaction mechanisms are explored with density functional theory (DFT) calculations on microsolvated clusters in the gas phase and, for a selected set of mechanisms, with constrained Car-Parrinello molecular dynamics (CPMD) simulations in water solution. Helmholtz free energy differences are calculated using thermodynamic integration in a simulation box with 31 water molecules at 300 K. The mechanism of the reaction is investigated with an emphasis on whether FeO2+ attacks methanol at a C-H bond or at the O-H bond. We conclude that the most likely mechanism is attack by the oxido oxygen at the C-H bond ("direct CH mechanism"). We calculate an upper bound for the reaction Helmholtz free energy barrier in solution of 50 kJ/mol for the C-H hydrogen transfer, after which transfer of the O-H hydrogen proceeds spontaneously. An alternative mechanism, starting with coordination of methanol directly to Fe ("coordination OH mechanism"), cannot be ruled out, as it involves a reaction Helmholtz free energy barrier in solution of 44 +/- 10 kJ/mol. However, this coordination mechanism has the disadvantage of requiring a prior ligand substitution reaction, to replace a water ligand by methanol. Because of the strong acidity of [FeO(H2O)5]2+, we also investigate the effect of deprotonation of a first-shell water molecule. However, this is found to increase the barriers for all mechanisms.

Oxidative properties of FeO2+: electronic structure and solvation effects.

An electronic structure analysis is provided of the action of solvated FeO(2+), [FeO(H(2)O)(5)](2+), as a hydroxylation catalyst. It is emphasized that the oxo end of FeO(2+) does not form hydrogen bonds (as electron donor and H-bond acceptor) with H-bond donors nor with aliphatic C-H bonds, but it activates C-H bonds as an electron acceptor. It is extremely electrophilic, to the extent that it can activate even such poor electron donors as aliphatic C-H bonds, the C-H bond orbital acting as electron donor in a charge transfer type of interaction. Lower lying O-H bonding orbitals are less easily activated. The primary electron accepting orbital in a water environment is the 3sigma*alpha orbital, an antibonding combination of Fe-3d(z(2)) and O-2p(z), which is very low-lying relative to the pi*alpha compared with, for example, the sigma* orbital in O(2) relative to its pi*. This is ascribed to relatively small Fe-3d(z(2)) with O-2p(z) overlap, due to the nodal structure of the 3d(z(2)).The H-abstraction barrier is very low in the gas phase, but it is considerably enhanced in water solvent. This is shown to be due to strong screening effects of the dielectric medium, leading to relative destabilization of the levels of the charged [FeO(H(2)O)(5)](2+) species compared to those of the neutral substrate molecules, making it a less effective electron acceptor. The solvent directly affects the orbital interactions responsible for the catalytic reaction.

Hydroxyl radical and hydroxide ion in liquid water: a comparative electron density functional theory study.

Ab initio density functional theory molecular dynamics simulations of the solvated states of the hydroxyl radical and hydroxide ion are performed using the Becke-Lee-Yang-Parr (BLYP) exchange-correlation functional (Becke, A. D. Phys. Rev. A 1988, 38, 3098. Lee, C.; Yang, W.; Parr, R. G. Phys. Rev. B 1988, 37, 785). The structures of the solvation shells of the two species are examined. It is found that the OH radical forms a relatively well-defined solvation complex with four neighboring water molecules. Three of these molecules are hydrogen bonded to the OH, while the fourth is hemibonded via a three-electron two-centered bond between the oxygen atoms of the OH and water. The activity and the diffusion mechanism of the OH radical in water is discussed in comparison with the OH- ion. Although the results are partially influenced by the tendency of the BLYP density functional to overestimate hemibonded structure, the present simulations suggest that the widely accepted picture of rapid diffusion of OH radical in water through hydrogen exchange reaction may need to be reconsidered.

The merits of the frozen-density embedding scheme to model solvatochromic shifts.

We investigate the usefulness of a frozen-density embedding scheme within density-functional theory [J. Phys. Chem. 97, 8050 (1993)] for the calculation of solvatochromic shifts. The frozen-density calculations, particularly of excitation energies have two clear advantages over the standard supermolecule calculations: (i) calculations for much larger systems are feasible, since the time-consuming time-dependent density functional theory (TDDFT) part is carried out in a limited molecular orbital space, while the effect of the surroundings is still included at a quantum mechanical level. This allows a large number of solvent molecules to be included and thus affords both specific and nonspecific solvent effects to be modeled. (ii) Only excitations of the system of interest, i.e., the selected embedded system, are calculated. This allows an easy analysis and interpretation of the results. In TDDFT calculations, it avoids unphysical results introduced by spurious mixings with the artificially too low charge-transfer excitations which are an artifact of the adiabatic local-density approximation or generalized gradient approximation exchange-correlation kernels currently used. The performance of the frozen-density embedding method is tested for the well-studied solvatochromic properties of the n-->pi(*) excitation of acetone. Further enhancement of the efficiency is studied by constructing approximate solvent densities, e.g., from a superposition of densities of individual solvent molecules. This is demonstrated for systems with up to 802 atoms. To obtain a realistic modeling of the absorption spectra of solvated molecules, including the effect of the solvent motions, we combine the embedding scheme with classical molecular dynamics (MD) and Car-Parrinello MD simulations to obtain snapshots of the solute and its solvent environment, for which then excitation energies are calculated. The frozen-density embedding yields estimated solvent shifts in the range of 0.20-0.26 eV, in good agreement with experimental values of between 0.19 and 0.21 eV.

Modeling solvent effects on electron-spin-resonance hyperfine couplings by frozen-density embedding.

In this study, we investigate the performance of the frozen-density embedding scheme within density-functional theory [J. Phys. Chem. 97, 8050 (1993)] to model the solvent effects on the electron-spin-resonance hyperfine coupling constants (hfcc's) of the H2NO molecule. The hfcc's for this molecule depend critically on the out-of-plane bending angle of the NO bond from the molecular plane. Therefore, solvent effects can have an influence on both the electronic structure for a given configuration of solute and solvent molecules and on the probability for different solute (plus solvent) structures compared to the gas phase. For an accurate modeling of dynamic effects in solution, we employ the Car-Parrinello molecular-dynamics (CPMD) approach. A first-principles-based Monte Carlo scheme is used for the gas-phase simulation, in order to avoid problems in the thermal equilibration for this small molecule. Calculations of small H2NO-water clusters show that microsolvation effects of water molecules due to hydrogen bonding can be reproduced by frozen-density embedding calculations. Even simple sum-of-molecular-densities approaches for the frozen density lead to good results. This allows us to include also bulk solvent effects by performing frozen-density calculations with many explicit water molecules for snapshots from the CPMD simulation. The electronic effect of the solvent at a given structure is reproduced by the frozen-density embedding. Dynamic structural effects in solution are found to be similar to the gas phase. But the small differences in the average structures still induce significant changes in the computed shifts due to the strong dependence of the hyperfine coupling constants on the out-of-plane bending angle.