密度泛函理论计算研究表面应力对钯催化乙炔选择性加氢活性与选择性的影响（英文）

在石油催化裂解过程中,除了生成乙烯、丙烯及丁烯等烯烃,还会产生部分炔烃.目前工业上通常采用炔烃选择性加氢转化为对应的单烯烃,以除去其中炔烃.由于产品烯烃中的炔烃等杂质含量需极低,这就对用于加氢催化剂的活性和选择性提出了很高的要求,即催化剂需要选择性吸附炔烃并加氢,而不损失其中的烯烃.经过前期大量的基础研究工作,目前工业中炔烃选择性加氢应用最广泛的催化剂是负载型钯基催化剂.然而,单独的钯金属选择性并不理想,因而对其选择性以及活性进行调控成为了当前关注的研究课题.本文采用密度泛函理论计算结合微观反应动力学模拟手段,研究了钯金属表面应力存在条件下的活性与选择性,以及形成次表层物种的可能性和形成后的活性与选择性.研究发现,改变钯金属的晶格参数与表面应力,反应物、表面反应中间体和产物的吸附能都会产生相应的变化,且吸附能与晶格参数的变化存在线性关系,晶格参数越大,吸附越强.利用表面反应过渡态能量与初始态能量之间的线性关系,相应的乙炔加氢生成乙烯的反应速率可以通过微观反应动力学模拟得到.结果显示,不同晶格参数的钯催化剂催化乙炔加氢生成乙烯的反应活性位于相应火山型曲线的强吸附侧,即减弱乙炔和氢的吸附强度可提高乙烯的生成速率.在此基础上,本文研究了不同表面应力的钯催化剂在次表面吸附不同覆盖度碳原子和氢原子的情况,发现晶格参数越大越有利于碳原子和氢原子在次表面的吸附.同时,研究发现在次表面碳掺杂的条件下,不同表面应力条件下的钯催化剂的活性均有所增强.此外,由于乙烯在所有研究的钯催化剂表面脱附比进一步加氢容易,因而乙烯都可以选择性生成.     

Myers–Saito and Schmittel cyclizations of enyne-(hetero)-1,2,3-trienes: A DFT study on the structure–reactivity relationship

The Myers–Saito (C²–C⁷) cyclization and the alternative Schmittel (C²–C⁶) cyclization of enyne-1,2,3-triene 1R and its hetero analogues 2R–5R were investigated by using pure density functional theory method (BPW91) in connection with the 6-311G(d, p) basis set. It has been shown that heteroatom-containing reactants lower significantly reaction barriers for both cyclization modes and reduce the difference between the barrier of the C²–C⁷ cyclization mode and that of the C²–C⁶ one. The Myers–Saito cyclization of 4R is associated with the smallest reaction barrier and the highest exothermicity. Whereas the Schmittel cyclization of 4R is exothermic, all the others are predicted to be endothermic.     

DFT study of the reactivity of methane and dioxygen with d-L2M complexes

A density functional theory analysis of the reactions of methane and O₂ with d¹⁰-L₂M complexes (M = Pd, Pt; L = N-heterocyclic carbene (NHC), PMe₃) is presented. Computations suggest that reaction of L₂M⁰ with O₂/CH₄ to form cis- (L)₂M(OOH)(CH₃) is only slightly uphill (ΔG ∼ 10-11 kcal/mol). Based on calculated thermodynamics, reaction of (L)₂Pt⁰ with CH₄ and O₂ is predicted to be more favorable by first addition of CH₄ and then reaction of O₂ with the resulting methyl-hydrido complex. However, oxidative addition for either the C-H bond of methane or of O₂ to d¹⁰-L₂M complexes is kinetically prohibitive. If barriers to oxidative addition to d¹⁰-L₂M systems could be reduced, conversion of L₂M(H)(CH₃) to L₂M(OOH)(CH₃) via hydrogen atom abstract/radical rebound is calculated to be thermodynamically and kinetically feasible, particularly for NHC and Pd. As (NHC)₂Pd also has a lower free energy to methane C-H oxidative addition than does (NHC)₂Pt, the former combination would appear to be a promising starting point in the search for catalysts for partial hydrocarbon oxidation.     

A DFT analysis of the vibrational spectra of nitrobenzene

Raman and FTIR, spectra of nitrobenzene, nb, and its isotopomers, nb-¹⁵N, nb-¹³C₆ and nb-d₅, were obtained and the fundamental vibrational modes assigned with the aid of a B3LYP/6-311+G** calculation, without the need for scaling of the force constants. The changes in vibrational coupling between the nitro and benzene groups upon certain isotopic substitutions are well modelled by the calculation, which is able to reproduce the isotopic shifts in frequencies for the nitro vibrations, as well as changes in IR intensities.     

Grand canonical ensemble approach to electrochemical thermodynamics,kinetics, and model Hamiltonians

The unique feature of electrochemistry is the ability to control reaction thermodynamics and kinetics by the application of electrode potential. Recently, theoretical methods and computational approaches within the grand canonical ensemble (GCE) have enabled to explicitly include and control the electrode potential in first principles calculations. In this review, recent advances and future promises of GCE density functional theory and rate theory are discussed. Particular focus is devoted to considering how the GCE methods either by themselves or combined with model Hamiltonians can be used to address intricate phenomena such as solvent/electrolyte effects and nuclear quantum effects to provide a detailed understanding of electrochemical reactions and interfaces.     

CO2 capture properties of M-C-O-H (M=Li, Na, K) systems: A combined density functional theory and lattice phonon dynamics study

We have computed the phase diagrams for multi-component M-C-O-H (M=Li, Na, K) systems using first-principles density functional theory complemented with lattice phonon calculations. We have identified all CO₂ capture reactions that lie on the Gibbs free energy convex hull as a function of temperature and the partial pressures of CO₂ and H₂O. Our predicted phase diagrams for CO₂ capture reactions are in qualitative and in some instances quantitative agreement with experimental data. The Na₂CO₃/NaHCO₃ and K₂CO₃/KHCO₃ systems were found to be the most promising candidates of all those we investigated for both pre- and post-combustion CO₂ capture. Overall, we show that our calculation approach can be used to screen promising materials for CO₂ capture under different conditions of temperature and pressure.     

A DFT variational approach to Hooke''s atom based on local-scaling transformations

The local-scaling transformation version of density functional theory, LS-DFT, is employed in order to construct energy functionals for Hooke's atom. The components of the energy are analyzed and the resulting exchange and correlation potentials are compared with the exact ones. In addition, the representation of the exact one-particle density in terms of the various components of the total energy density is discussed.     

Recent progress on first-principles simulations of voltammograms

In this review article, we provide a brief summary of recent models and applications using first-principles methods to simulate cyclic voltammograms and linear sweep voltammograms. Electrochemical processes considered include hydrogen/hydroxyl adsorption, higher order anion adsorption, and electrocatalytic reactions. Reasonable agreement between experiment and theory is achieved in most systems, and we discuss approximations used and routes for improving the predictive power and accuracy of simulated voltammograms.     

Gas-phase properties and Fukui indices of sulfine (CH2SO). Potential energy surface and maximum hardness principle for its protonation process. A density functional study

Gas-phase thermochemical properties of sulfine (CH₂SO) and the potential energy surface of its protonation process were studied by the density functional method employing different exchange-correlation potentials. All calculations showed that the most stable protonated isomer is planar with the proton bonded to the oxygen atom in a trans arrangement of the skeleton. Three transition states were located that allow interconversion between the different isomers. Hardnesses and Fukui indices were calculated to follow the reactivity trend along the protonation path and to explain the preference for a particular protonation site on neutral sulfine. Proton affinity, gas-phase basicity and heat of formation values, obtained for the first time fully quantum mechanically, agree well with those derived by a recent mass spectrometry experimental study. Good agreement between density functional theory and previous high-level theoretical and experimental data was also found for the heat of formation of sulfine and its most stable protonated form. Received: 12 October 1998 / Accepted: 24 November 1998 / Published online: 16 March 1999     

A theoretical study on the regioselectivity of 1,3-dipolar cycloadditions using DFT-based reactivity indexes

The regioselectivity for a series of four 1,3-dipolar cycloaddition reactions has been studied using global and local reactivity indexes. The results of the theoretical analysis suggest that for asynchronous cycloadditions associated to polar processes, the regioselectivity is consistently explained by the most favorable two-center interactions between the highest nucleophilic and electrophilic sites of the reagents.     

Geometry and tautomerism of sapphyrin—DFT studies

The structure and total energy have been investigated, applying the density functional theory (DFT), for idealized sapphyrin derivatives, created by replacements of meso-phenyl groups of 5,10,15,20-tetraphenylsapphyrin or β-alkyl groups of decaalkylsapphyrin with hydrogens or two methyl groups, respectively. The following species have been considered: the dication of plain sapphyrin (SapH₅²⁺) 12,13-dimethylsapphyrin (12,13-DMSapH₃), 12,13-dimethylsapphyrin dication (12,13-DMSapH₅²⁺), and the 10,15-dimethylsapphyrin dication (10,15-DMSapH₅²⁺). To address the issue of tautomerization the calculations have been carried out for the six feasible distribution patterns of NH protons among the five nitrogen centers of the planar sapphyrin P0_i (i=1–6). The B3LYP/6-31G optimized bond lengths and angles of sapphyrin or (dimethylated sapphyrin) skeletons are in satisfactory agreement with X-ray crystallographic values determined for decaalkylsapphyrin dications. The calculated total energies, using the B3LYP/6-31G**//3-21G approach demonstrate that the plain sapphyrin and 12,13-dimethylsapphyrin dications favor the planar geometry of the macrocycle. Localization of the methyl groups at the 10,15-meso positions of sapphyrin reversed the order of stability leading to the energy preference for the inverted structure. Thus, this steric factor related to the meso substitution seems to be instrumental in the C pyrrolic ring inversion. The DFT calculation also provided the relative stability frame for the complex tautomeric equilibria of SapH₃ which involve six different forms. Their relative stability decreases in the range: {25-NH, 26-N, 27-NH, 28-NH, 29-N} P0₃>({25-NH, 26-N, 27-NH, 28-N, 29-NH} P0₁>{25-NH, 26-NH, 27-N, 28-NH, 29-N} P0₅>{25-N, 26-NH, 27-NH, 28-NH, 29-N} P0₂>{25-NH, 26-N, 27-N, 28-NH, 29-NH} P0₆>{25-N, 26-N, 27-NH, 28-NH, 29-NH} P0₄. The stability order can be qualitatively related to the number and the nature of the cis NH interactions present in the sapphyrin core.     

Influence of Mn-doping on the chemical shielding tensors of Ga and N in the armchair GaN nanotubes: A DFT study

The chemical shielding (CS) tensors of Gallium-71 and nitrogen-15 are computed for the first time in order to investigate the influence of Mn-doping on the electronic properties of the (5, 5) Gallium nitride nanotube (GaNNT). A GaNNT consisting of 40 Ga and 40 N atoms and having a 1.2 nm length was considered. One portal of the nanotube was capped by ten hydrogen atoms and other-end was kept open. Additionally, two other forms of this model of Mn-doped GaNNT were considered where a Mn-atom was substituted for a Ga atom either in the first or in the second layer. The calculations reveal that in both models of Mn-doped GaNNTs, the N atoms that are directly connected to the Mn atom have the smallest isotropic chemical shielding among other N atoms. These calculations were performed at the level of the density functional theory (DFT) using GAUSSIAN 03 package. The basis sets for Ga and N atoms were chosen to be 6-31G (d) and those for Mn atom were chosen to be LanL2DZ.     

Reactivity of the carbon-carbon double bond towards nucleophilic additions. A DFT analysis

The global and local electrophilicity indexes have been used to characterize the reactivity pattern of the CC double bond towards nucleophilic addition reactions. A wide family of molecules including ketones, esters, anhydrides, nitriles and nitrocompounds containing appropriate substitution on the CC double bond have been classified within an unique scale of reactivity. The predictive capability of the theoretical model is tested against a series of benzylidenemalononitriles and substituted α-nitrostilbenes.     

Lithium 2,3-dihydro-1-benzothiophene-1,1-dioxide: synthesis, characterization, DFT calculations, and reactivity toward aldehydes and azomethines

The sulfonyl carbanion derived from 2,3-dihydro-1-benzothiophene-1,1-dioxide and its lithium salt has been investigated by DFT calculations. NMR and IR spectroscopic analyses showed that the lithium sulfonyl carbanion exists in solution as a monomer in equilibrium with a dimer. The lithium carbanion was treated with aldehydes and azomethynes to give chiral hydroxy and amino derivatives. The stereochemistry of the products and the diastereoselectivity of the reaction were investigated.     

Single configuration interaction study on conjugated betainic chromophores based on DFT optimized geometries

The structure of a series of heterocyclic betaines was calculated by methods of density functional theory (DFT). The charge distribution and bond characteristics of these compounds were analyzed by Weinhold’s natural bond orbital analysis (NBO) and by natural resonance theory (NRT). In order to probe the aromatic character of the ring fragments, Schleyer’s nucleus-independent chemical shifts (NICSs) were calculated by GIAO-RHF. Ab initio single configuration interaction calculations (SCI) correctly predict intense π→π* transitions at low energies, but the transition energies of the color bands are overestimated. Torsion around the interfragmental bond increases the charge separation between the molecular fragments and the dipole moment. The molecular fragments become increasingly aromatic. The absorption wavelengths increase on torsion while the oscillator strengths decrease.     

A NMR, X-ray, and DFT combined study on the regio-chemistry of nucleophilic addition to platinum(II) coordinated terminal olefins

The series of platinum complexes [PtCl(η²-CH₂CH-C₆H₄-X)(tmeda)](ClO₄) (X = H, 1b; 4-OMe, 1c; 3-OMe, 1d; 4-CF₃, 1e; 3-CF₃, 1f; 3-NO₂, 1g; tmeda = N,N,N′,N′-tetramethyl-1,2-ethanediamine) has been considered. In the styrene complex (1b) both solution (NMR) and solid state (X-ray) data indicate a significant difference in the Pt-C bond lengths (the longer bond being that involving the olefin carbon atom carrying the phenyl ring). Such a difference increases when X is an electron donor group (EDG, 1c) and decreases when X is an electron withdrawing group (EWG, 1d-g). The attack of a nucleophile (MeO⁻) to the substituted carbon (Markovnikov type, M) is by far the most favoured in the case of unsubstituted (1b) or EDG-substituted (1c) styrenes. The presence of an EWG (compounds 1d-g) levels off the probability of M and anti-M type of attack. DFT calculations on 1b,c and 1e were also performed. The NLMO analysis reveals the crucial role of the interaction between the filled π orbital of the olefin and the empty d orbital of platinum; the carbon with greater electron density becoming less susceptible of nucleophilic attack.     

Local and nonlocal density functional calculations of the molecular structure of isomeric thiadiazoles

The chemistry of thiadiazoles and their derivatives is of considerable interest in chemistry owing to their pharmacological and potential industrial applications. In this context, a detailed study of isomeric thiadiazole molecules has been done using local (SVWN; Slater, and Vosko, Wilk and Nusair) and nonlocal (BLYP; Becke, and Lee, Yang and Parr) density functionals and optimizing the molecular geometries by means of the gradient technique. A charge sensitivity analysis of the studied molecule has been performed by resorting to density functional theory, obtaining several sensitivity coefficients such as the molecular energy, net atomic charges, global and local hardness, global and local softness and Fukui functions. With these results and the analysis of the dipole moments, the molecular electrostatic potentials and the total electron density maps, several conclusions have been inferred about the preferred sites of chemical reaction of the studied compounds. The condensed Fukui functions are shown to be one of the best criteria for predicting chemical reactivity.     

The C−NO2 bond dissociation energies of some nitroaromatic compounds: DFT study

The C−NO₂ bond dissociation energies in nitrobenzene; 3-amino-nitrobenze; 4-amino-nitrobenze; 1,3-dinitrobenzene; 1,4-dinitrobenzene; 2-methyl-nitrobenzene; 4-methyl-nitrobenzene; and 1,3,5-trinitrobenzene nitroaromatic molecules, are computed using B3LYP, B3PW91, B3P86 three-parameter hybrid Density Functional Theory (DFT) methods in conjunction with 6-31G^** basis set. By comparing the computed energies and experimental ones, it is found that B3P86/6-31G^** is not capable of predicting the satisfactory bond dissociation energy (BDE). The BDEs computed with both B3LYP/6-31G^** and B3PW91/6-31G^** for the nitroaromatic molecules are closer to the experimental ones than those obtained with B3P86/6-31G^**. But, when compared with the experimental one, the BDE from the B3LYP/6-31G^** has the maximum deviation, which is completely outside our desired target accuracy for chemical predictions (less than 2.00 kcal mol⁻¹). Therefore, we suggest B3PW91/6-31G^** method as a reliable method of computing the BDE for removal of the nitrogen dioxide group in the nitroaromatic compounds. In addition, the C−NO₂ BDEs for 2,4,6-trinitrotoluene (TNT), triaminotrinitrobenzene (TATB), diaminotrinitrobenzene (DATB), and picramide are studied with B3PW91/6-31G^** method.     

Chlorofluorocarbons adsorption structures and energetic over faujasite type zeolites—a first principle study

Adsorption of ozone depleting chlorofluorocarbons (CFC) over zeolite is of major global environmental concern. The current communication describes first-principle calculation performed on faujasite models to investigate the nature of CFCs including fluoro, chlorofluoro and hydrofluoro/chloro carbons (CF₄, CF₃Cl, CF₂Cl₂, CFCl₃, CHF₃, CHCl₃) adsorption. Experimentally it is observed that separation of halocarbons are possible using Na–Y, though the cause is unknown. Reactivity index within the helm of Hard Soft Acid Base (HSAB) principle was used to monitor the activity of the interacting CFCs using Density Functional Theory (DFT) to propose a qualitative order. The importance of both H-bonding and cation–F/Cl interactions in determining the low-energy sorption sites were monitored and rationalized. The host guest interactions show a distinctive difference between the adsorption phenomenon between H–Y and Na–Y and as well for Cl and F. It is observed that Cl has more favorable interaction with hydrogen of H–Y compared to Na–Y and for F the situation is just reversed. To validate this trend periodic optimization calculations were performed. The interaction energy as obtained matches well with the reactivity index order resulted from cluster calculations. This study is a combination of DFT and periodic calculation to rationalize the electronic phenomenon of the interaction process.