Thermodynamic analysis of mixed and dry reforming of methane for solar thermal applications

Thermodynamic analysis of the reforming of methane with carbon dioxide alone ("dry reforming") and with carbon dioxide and steam together ("mixed reforming") is performed as part of a project which investigates the suitability of these endothermic reactions for the storage of solar thermal energy. The Gibbs free energy minimization method was employed to identify thermodynamically optimal operating conditions for dry reforming as well as mixed reforming with a desired H2/CO molar ratio of 2. The non-stoichiometric equilibrium model was developed using FactSage software to conduct the thermodynamic calculations for carbon formation, H2/CO ratio, CH4 conversion and H2 yield as a function of reaction temperature, pressure and reactant molar ratios. Thermodynamic calculations demonstrate that in the mixed reforming process, optimal operating conditions in a carbon-free zone are under H2O/CH4 /CO2 =1.0/1.0/0.5, p = 1 to 10 bar and T = 800 to 850℃ for the production of syngas with a H2 /CO molar ratio of 2. Under the optimal conditions, the maximum H2 yield of 88.0% is achieved at 1 bar and 850℃ with a maximum CH4 conversion of 99.3%. In the dry reforming process, a carbon formation regime is always present at a CO2/CH4 molar ratio of 1 for T = 700 1000℃ and p = 1-30 bar, whereas a carbon-free regime can be obtained at a CO2/CH4 molar ratio greater than 1.5 and T≥800℃.     

Rotational bands in 167Hf

High-spin states in ¹⁶⁷Hf, populated in the ¹⁴¹Pr(³⁰Si, p3n)¹⁶⁷Hf reaction, have been studied using the nordball Ge detector array. Three rotational cascades have been observed for the first time and the previously-known level scheme has been extended to significantly higher spin. Band-crossing effects are discussed within the framework of Woods-Saxon cranking calculations and are found to be in good agreement. Received: 9 April 1999 / Revised version: 13 May 1999     

Binding affinities in the SAMPL3 trypsin and host–guest blind tests estimated with the MM/PBSA and LIE methods

We have estimated affinities for the binding of 34 ligands to trypsin and nine guest molecules to three different hosts in the SAMPL3 blind challenge, using the MM/PBSA, MM/GBSA, LIE, continuum LIE, and Glide score methods. For the trypsin challenge, none of the methods were able to accurately predict the experimental results. For the MM/GB(PB)SA and LIE methods, the rankings were essentially random and the mean absolute deviations were much worse than a null hypothesis giving the same affinity to all ligand. Glide scoring gave a Kendall's τ index better than random, but the ranking is still only mediocre, τ = 0.2. However, the range of affinities is small and most of the pairs of ligands have an experimental affinity difference that is not statistically significant. Removing those pairs improves the ranking metric to 0.4-1.0 for all methods except CLIE. Half of the trypsin ligands were non-binders according to the binding assay. The LIE methods could not separate the inactive ligands from the active ones better than a random guess, whereas MM/GBSA and MM/PBSA were slightly better than random (area under the receiver-operating-characteristic curve, AUC = 0.65-0.68), and Glide scoring was even better (AUC = 0.79). For the first host, MM/GBSA and MM/PBSA reproduce the experimental ranking fairly good, with τ = 0.6 and 0.5, respectively, whereas the Glide scoring was considerably worse, with a τ = 0.4, highlighting that the success of the methods is system-dependent.     

Fast generation of broken‐symmetry states in a large system including multiple iron–sulfur assemblies: Investigation of QM/MM energies,clusters charges,and spin populations

A density functional theory study is presented regarding the energetics and the Mulliken population analyses of a quantum mechanical/molecular mechanical (QM/MM) system including multiple iron–sulfur clusters in the QM region. The [FeFe]‐hydrogenase from Desulfovibrio desulfuricans was studied, and both the active site (an Fe₆S₆ assembly generally referred to as the H‐cluster) and an ancillary Fe₄S₄ site were treated at the BP86‐RI/TZVP level. The antiferromagnetic coupling that characterizes both sites was modeled using the broken‐symmetry (BS) approach. For such a QM system, 36 different BS couplings can be defined, depending on the localization of spin excess on the various spin centers. All the BS states were obtained by means of an effective and simple method for spin localization, that is here described and compared with more sophisticated approaches already available in literature. The variation of the QM/MM energy among the various geometry‐optimized protein models was found to be less than 25 kJ mol^–1. This energy variation almost doubles if no geometry optimization is performed. A detailed analysis of the additive nature of these variations in QM/MM energy is reported. The Mulliken charges show very small variations among the 36 BS states, whereas the Mulliken spin populations were found to be somewhat more variable. The relevance of such variations is discussed in light of the available Mössbauer and Electron Paramagnetic Resonance (EPR) spectroscopic data for the enzyme. Finally, the influence of the basis set on the spin populations, charges, and structural parameters of the models was investigated, by means of QM/MM computations on the same system at the BP86‐RI/SVP level. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Isocyanide in Biochemistry? A Theoretical Investigation of the Electronic Effects and Energetics of Cyanide Ligand Protonation in [FeFe]‐Hydrogenases

The presence of Fe‐bound cyanide ligands in the active site of the proton‐reducing enzymes [FeFe]‐hydrogenases has led to the hypothesis that such Brønsted–Lowry bases could be protonated during the catalytic cycle, thus implying that hydrogen isocyanide (HNC) might have a relevant role in such crucial microbial metabolic paths. We present a hybrid quantum mechanical/molecular mechanical (QM/MM) study of the energetics of CN^? protonation in the enzyme, and of the effects that cyanide protonation can have on [FeFe]‐hydrogenase active sites. A detailed analysis of the electronic properties of the models and of the energy profile associated with H₂ evolution clearly shows that such protonation is dysfunctional for the catalytic process. However, the inclusion of the protein matrix surrounding the active site in our QM/MM models allowed us to demonstrate that the amino acid environment was finely selected through evolution, specifically to lower the Brønsted–Lowry basicity of the cyanide ligands. In fact, the conserved hydrogen‐bonding network formed by these ligands and the neighboring amino acid residues is able to impede CN^? protonation, as shown by the fact that the isocyanide forms of [FeFe]‐hydrogenases do not correspond to stationary points on the enzyme QM/MM potential‐energy surface.     

Free-energy perturbation and quantum mechanical study of SAMPL4 octa-acid host–guest binding energies

We have estimated free energies for the binding of nine cyclic carboxylate guest molecules to the octa-acid host in the SAMPL4 blind-test challenge with four different approaches. First, we used standard free-energy perturbation calculations of relative binding affinities, performed at the molecular-mechanics (MM) level with TIP3P waters, the GAFF force field, and two different sets of charges for the host and the guest, obtained either with the restrained electrostatic potential or AM1-BCC methods. Both charge sets give good and nearly identical results, with a mean absolute deviation (MAD) of 4 kJ/mol and a correlation coefficient (R ²) of 0.8 compared to experimental results. Second, we tried to improve these predictions with 28,800 density-functional theory (DFT) calculations for selected snapshots and the non-Boltzmann Bennett acceptance-ratio method, but this led to much worse results, probably because of a too large difference between the MM and DFT potential-energy functions. Third, we tried to calculate absolute affinities using minimised DFT structures. This gave intermediate-quality results with MADs of 5–9 kJ/mol and R ² = 0.6–0.8, depending on how the structures were obtained. Finally, we tried to improve these results using local coupled-cluster calculations with single and double excitations, and non-iterative perturbative treatment of triple excitations (LCCSD(T0)), employing the polarisable multipole interactions with supermolecular pairs approach. Unfortunately, this only degraded the predictions, probably because of a mismatch between the solvation energies obtained at the DFT and LCCSD(T0) levels.