Activity and coke formation of nickel and nickel carbide in dry reforming: A deactivation scheme from density functional theory

Dry reforming is a promising reaction to utilise the greenhouse gases CO2 and CH4. Nickel-based catalysts are the most popular catalysts for the reaction, and the coke formation on the catalysts is the main obstacle to the commercialisation of dry reforming. In this study, the whole reaction network of dry reformation on both flat and stepped nickel catalysts (Ni(111) and Ni(211)) as well as nickel carbide (flat: Ni3C(001); stepped: Ni3C(111)) is investigated using density functional theory calculations. The overall reaction energy profiles in the free energy landscape are obtained, and kinetic analyses are utilised to evaluate the activity of the four surfaces. By careful examination of our results, we find the following regarding the activity: (i) flat surfaces are more active than stepped surfaces for the dry reforming and (ii) metallic nickel catalysts are more active than those of nickel carbide, and therefore, the phase transformation from nickel to nickel carbide will reduce the activity. With respect to the coke formation, the following is found: (i) the coke formation probability can be measured by the rate ratio of CH oxidation pathway to C oxidation pathway (r(CH)/r(C)) and the barrier of CO dissociation, (ii) on Ni(111), the coke is unlikely to form, and (iii) the coke formations on the stepped surfaces of both nickel and nickel carbide can readily occur. A deactivation scheme, using which experimental results can be rationalised, is proposed. 

A density functional theory study of hydrogen dissociation and diffusion at the perimeter sites of Au/TiO2

Reactivity of supported gold catalysts is a hot topic in catalysis for many years. This communication reports an investigation on the dissociation of molecular hydrogen at the perimeter sites of Au/TiO2 and the spillover of hydrogen atoms from the gold to the support using density functional theory calculations. It is found that the heterolytic dissociation is favoured in comparison with homolytic dissociation of molecular hydrogen at the perimeter sites. However, the surface oxygen of the rutile TiO2(110) surface at these sites can be readily passivated by the formed OH, suggesting that further dissociation of molecular hydrogen may occur at pure gold sites.

Some Understanding of Fischer-Tropsch Synthesis from Density Functional Theory Calculations

The combination of density functional theory (DFT) calculations and kinetic analyses is a very useful approach to study surface reactions in heterogeneous catalysis. The present paper reviews some recent work applying this approach to Fischer-Tropsch (FT) synthesis. Emphasis is placed on the following fundamental issues in FT synthesis: (i) reactive sites for both hydrogenation and C-C coupling reactions; (ii) reaction mechanisms including carbene mechanism, CO-insertion mechanism and hydroxyl-carbene mechanism; (iii) selectivity with a focus on CH(4) selectivity, alpha-olefin selectivity and chain growth probability; and (iv) activity.

An Energy Descriptor To Quantify Methane Selectivity in Fischer-Tropsch Synthesis:A Density Functional Theory Study

Selectivity is a fundamental issue in heterogeneous catalysis. In this study, the CH(4) selectivity in Fischer-Tropsch synthesis is chosen to be investigated: CH4 selectivity on Rh, Co, Ru, Fe, and Re surfaces is computed by first-principles methods. In conjunction with kinetic analyses, we are able to derive the effective barrier difference between methane formation and chain growth (Delta E(eff)) to quantify the CH(4) selectivity. By using this energy descriptor, the ranking of methane selectivity predicted from density functional theory (DFT) calculations is consistent with experimental work. Moreover, a linear correlation between Delta E(eff) and the chemisorption energy of C + 4H (Delta H) is found. This fundamental finding possesses the following significance: (i) it shows that the selectivity, which appears to have kinetic characteristics, is largely determined by thermodynamic properties; and (ii) it suggests that an increase of the binding strength of C + 4H will suppress methane selectivity.

A density functional theory study of CH2 and H adsorption on Ni(111)

Ab initio total energy calculations within the density functional theory framework have been used to study the adsorption of CH2 and H as well as the coadsorption of CH2 and H on Ni(111). H binds strongly at threefold hollow sites with calculated adsorption energies of 2.60 and 2.54 eV at the face-centered-cubic (fcc) and hexagonal-close-packed (hcp) hollow sites, respectively. Adsorption energies and H-Ni distances are found to agree well with both experimental and theoretical results. CH2 adsorbs strongly at all high symmetry sites with calculated adsorption energies of 3.26, 3.22, 3.14 and 2.36 eV at the fcc, hcp, bridge and top sites, respectively. Optimized structures are reported at all sites, and, in the most stable hollow sites there is considerable internal reorganization of the CH2 fragment. The CH2 molecule is tilted, the hydrogens are inequivalent and the C-H bonds are lengthened relative to the gas phase. In the CH2-H coadsorption systems the adsorbates have a tendency to move toward bridge sites. The bonding of all adsorbates to the surface is analyzed in detail. (C) 2000 American Institute of Physics. [S0021-9606(00)71213-X].

CO oxidation on Pt(111):An ab initio density functional theory study

CO oxidation on Pt(111) is studied with ab initio density functional theory. The low energy pathway and transition state for the reaction are identified. The key event is the breaking of an O-metal bond prior to the formation of a chemisorbed CO2 molecule. The pathway can be rationalized in terms of competition of the O and C atoms for bonding with the underlying surface, and the predominant energetic barrier is the strength of the O-metal bond.

GRADIENT CORRECTIONS IN DENSITY-FUNCTIONAL THEORY CALCULATIONS FOR SURFACES - CO ON PD(110)

Ab initio total energy calculations have been performed for CO chemisorption on Pd(110). Local density approximation (LDA) calculations yield chemisorption energies which are significantly higher than experimental values but inclusion of the generalised gradient approximation (GGA) gives better agreement. In general, sites with higher coordination of the adsorbate to surface atoms lead to a larger degree of overbinding with LDA, and give larger corrections with GGA. The reason is discussed using a first-order perturbation approximation. It is concluded that this may be a general failure of LDA for chemisorption energy calculations. This conclusion may be extended to many surface calculations, such as potential energy surfaces for diffusion.

Turbulent heat transfer and coherent vortical structures in gas turbine film cooling

Unsteady coherent structures and turbulent heat transfer in a film cooling flow is studied by using detached eddy simulation (DES). Detailed computations for an inclined jet in crossflow by a single row of 35 degree round holes on a flat plate were performed at blowing ratios of 0.5 and 1.0, and a density ratio of 2.0. The correlation between the coherent vortical structures and the unsteady heat transfer is carefully examined. The instantaneous flow fields and heat transfer distributions are found to be characterized by the formation of large coherent vortical structures. These structures enhance the thermal mixing process and turbulent heat transfer to the wall. From the inspection of both unsteady adiabatic film cooling effectiveness and heat transfer coefficient, these two are found to have substantial local fluctuations due to the large unsteadiness of coherent structures. The fluctuation of the adiabatic effectiveness and heat transfer coefficient, for example, can be as high as 15 and 50 percent of the time-mean value, respectively. It could result in the detrimental effect on film cooling performance.

First Detection of Gas-phase Methanol in a Protoplanetary Disk

The first detection of gas-phase methanol in a protoplanetary disk (TW Hya) is presented. In addition to being one of the largest molecules detected in disks to date, methanol is also the first disk organic molecule with an unambiguous ice chemistry origin. The stacked methanol emission, as observed with the Atacama Large Millimeter/submillimeter Array, is spectrally resolved and detected across six velocity channels (>3σ), reaching a peak signal-to-noise of 5.5σ, with the kinematic pattern expected for TW Hya. Using an appropriate disk model, a fractional abundance of 3 x 10^-12 – 4 x 10^-11 (with respect to H₂) reproduces the stacked line profile and channel maps, with the favored abundance dependent upon the assumed vertical location (midplane versus molecular layer). The peak emission is offset from the source position, suggesting that the methanol emission has a ring-like morphology: the analysis here suggests it peaks at ≈30 au, reaching a column density ≈3–6 x 10¹² cm⁻². In the case of TW Hya, the larger (up to millimeter-sized) grains, residing in the inner 50 au, may thus host the bulk of the disk ice reservoir. The successful detection of cold gas-phase methanol in a protoplanetary disk implies that the products of ice chemistry can be explored in disks, opening a window into studying complex organic chemistry during planetary system formation.

Gas Solubility and Vapor Pressure of Electrolytes Based on Alkylcarbonates Containing Litfsi, Lifap, and LiPF6 Lithium Salt: Measurement and Prediction

Most liquid electrolytes used in commercial lithium-ion batteries are composed by alkylcarbonate mixture containing lithium salt. The decomposition of these solvents by oxidation or reduction during cycling of the cell, induce generation of gases (CO2, CH4, C2H4, CO …) increasing of pressure in the sealed cell, which causes a safety problem [1]. The prior understanding of parameters, such as structure and nature of salt, temperature pressure, concentration, salting effects and solvation parameters, which influence gas solubility and vapor pressure of electrolytes is required to formulate safer and suitable electrolytes especially at high temperature.

We present in this work the CO2, CH4, C2H4, CO solubility in different pure alkyl-carbonate solvents (PC, DMC, EMC, DEC) and their binary or ternary mixtures as well as the effect of temperature and lithium salt LiX (X = LiPF6, LiTFSI or LiFAP) structure and concentration on these properties. Furthermore, in order to understand parameters that influence the choice of the structure of the solvents and their ability to dissolve gas through the addition of a salt, we firstly analyzed experimentally the transport properties (Self diffusion coefficient (D), fluidity (h-1), and conductivity (s) and lithium transport number (tLi) using the Stock-Einstein, and extended Jones-Dole equations [2]. Furthermore, measured data for the of CO2, C2H4, CH4 and CO solubility in pure alkylcarbonates and their mixtures containing LiPF6; LiFAP; LiTFSI salt, are reported as a function of temperature and concentration in salt. Based on experimental solubility data, the Henry’s law constant of gases in these solvents and electrolytes was then deduced and compared with values predicted by using COSMO-RS methodology within COSMOthermX software. From these results, the molar thermodynamic functions of dissolution such as the standard Gibbs energy, the enthalpy, and the entropy, as well as the mixing enthalpy of the solvents and electrolytes with the gases in its hypothetical liquid state were calculated and discussed [3]. Finally, the analysis of the CO2 solubility variations with the salt addition was then evaluated by determining specific ion parameters Hi by using the Setchenov coefficients in solution. This study showed that the gas solubility is entropy driven and can been influenced by the shape, charge density, and size of the anions in lithium salt.

References

[1] S.A. Freunberger, Y. Chen, Z. Peng, J.M. Griffin, L.J. Hardwick, F. Bardé, P. Novák, P.G. Bruce, Journal of the American Chemical Society 133 (2011) 8040-8047.

[2] P. Porion, Y.R. Dougassa, C. Tessier, L. El Ouatani, J. Jacquemin, M. Anouti, Electrochimica Acta 114 (2013) 95-104.

[3] Y.R. Dougassa, C. Tessier, L. El Ouatani, M. Anouti, J. Jacquemin, The Journal of Chemical Thermodynamics 61 (2013) 32-44.

Modeling of Turbulent Natural Gas and Biogas Flames of the Delft Jet-in-Hot-Coflow Burner: Effects of Coflow Temperature, Fuel Temperature and Fuel Composition on the Flame Lift-Off Height

The structure of a turbulent non-premixed flame of a biogas fuel in a hot and diluted coflow mimicking moderate and intense low dilution (MILD) combustion is studied numerically. Biogas fuel is obtained by dilution of Dutch natural gas (DNG) with CO2. The results of biogas combustion are compared with those of DNG combustion in the Delft Jet-in-Hot-Coflow (DJHC) burner. New experimental measurements of lift-off height and of velocity and temperature statistics have been made to provide a database for evaluating the capability of numerical methods in predicting the flame structure. Compared to the lift-off height of the DNG flame, addition of 30 % carbon dioxide to the fuel increases the lift-off height by less than 15 %. Numerical simulations are conducted by solving the RANS equations using Reynolds stress model (RSM) as turbulence model in combination with EDC (Eddy Dissipation Concept) and transported probability density function (PDF) as turbulence-chemistry interaction models. The DRM19 reduced mechanism is used as chemical kinetics with the EDC model. A tabulated chemistry model based on the Flamelet Generated Manifold (FGM) is adopted in the PDF method. The table describes a non-adiabatic three stream mixing problem between fuel, coflow and ambient air based on igniting counterflow diffusion flamelets. The results show that the EDC/DRM19 and PDF/FGM models predict the experimentally observed decreasing trend of lift-off height with increase of the coflow temperature. Although more detailed chemistry is used with EDC, the temperature fluctuations at the coflow inlet (approximately 100K) cannot be included resulting in a significant overprediction of the flame temperature. Only the PDF modeling results with temperature fluctuations predict the correct mean temperature profiles of the biogas case and compare well with the experimental temperature distributions.