A theoretical study of CO adsorption on gold by Hückel theory and density functional theory calculations

It is crucial to understand the nature of CO adsorption on gold so as to elucidate the mechanism of low-temperature CO oxidation on nanogold catalysts. We performed theoretical analysis of CO adsorption on gold by using Hückel theory and density functional theory (DFT) calculations. Hückel theory indicates that CO adsorption on gold is dominated by the electron distribution at the Au atom, which is greatly affected by neighboring Au atoms, coadsorbed or doping species. The increase of σ-bonding electrons should weaken the CO adsorption, while the increase of π-electrons should strengthen the adsorption. DFT calculations proved this prediction quantitatively for various systems, including CO adsorption on a Au(100)-hex surface with locally varying subsurface configurations and CO coadsorption with acceptor or donor species.     

Intermolecular potentials of the silane dimer calculated with Hartree-Fock theory, Møller-Plesset perturbation theory, and density functional theory

We have calculated the intermolecular interaction potentials of the silane dimer at the D3d conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with 108 functionals chosen from the combinations of 9 exchange and 12 correlation functionals. Single-point coupled cluster [CCSD(T)] calculations have also been carried out to calibrate the correlation effect. The HF calculations yield unbound potentials largely because of the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater type orbitals fitted with Gaussian functions (STO-nG, n = 3 approximately 6), Pople's medium size basis sets [up to 6-311++G(3df,3pd)], to Dunning's correlation consistent basis sets (cc-pVXZ and aug-cc-pVXZ, X = D, T, Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(3d,3p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy ( approximately 0.05 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the expected C6 value from molecular polarizability calculations. We attribute the slow convergence partly to the inefficacy of using the MP2 calculations with Gaussian type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the expected potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energies calculated using the OPTXHCTH147, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals and the equilibrium bond lengths calculated using the MPWHCTH93, TPSSHCTH, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals are close to the MP2 results using the 6-311++G(3df,3pd) basis set. A correlation between the calculated DFT potentials and the exchange and correlation enhancement factors at the low-density region has been elucidated. The asymptotic behaviors of the DFT potentials are also analyzed.     

Surface-enhanced Raman scattering of 4,4′-bipyridine on silver by density functional theory calculations

Surface-enhanced Raman scattering (SERS) of 4,4′-bipyridine (BPy) on silver foil substrate was measured using the 488, 514.5, and 1064 nm excitation lines. Density functional theory (DFT) methods were used to calculate the structure and vibrational spectra of Ag–BPy, Ag₃–BPy and Ag₄–BPy complexes with B3LYP/6-31++G(d,p)(C,H,N)/Lanl2dz(Ag) basis set. The Raman bands of BPy were assigned on the basis of the calculation of potential energy distribution. The calculated spectra of Ag–BPy and Ag₄–BPy complexes were much closer to the experimental results of BPy adsorbed on silver surface than that of Ag₃–BPy complexes. The vibrational frequencies that are sensitive to the planar or non-planar structure of BPy and to the dihedral angle of two pyridyl rings were discussed. The DFT results showed that the angles between two pyridyl rings for Ag–BPy and Ag₄–BPy were skewed by about 38.44° and 37.1°, respectively. The energy gaps of the HOMO and LUMO from DFT were 415–912 nm for BPy–Ag complexes. The relative intensities of SERS bands changed with different excitation laser lines. Thus, a chemical enhancement mechanism should play an important role in the SERS of BPy on silver substrate.     

Influence of impurity elements on the corrosion of α-uranium surface: a density functional theory study

Journal of Radioanalytical and Nuclear Chemistry - Influence of substitutional Fe, Al and Mg atoms on the corrosion of α-U(110) surface was studied by DFT?+?U method. The...     

Intermolecular potentials of the methane dimer calculated with Møller-Plesset perturbation theory and density functional theory

We have calculated the intermolecular interaction potentials of the methane dimer at the minimum-energy D(3d) conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with the Perdew-Wang (PW91) functional as the exchange or the correlation part. The HF calculations yield unbound potentials largely due to the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater-type orbitals fitted with Gaussian functions (STO-nG) (n=3-6) [Quantum Theory of Molecular and Solids: The Self-Consistent Field for Molecular and Solids (McGraw-Hill, New York, 1974), Vol. 4], Pople's medium size basis sets of Krishnan et al. [J. Chem. Phys. 72, 650 (1980)] [up to 6-311++G(3df,3pd)] to Dunning's correlation consistent basis sets [J. Chem. Phys. 90, 1007 (1989)] (cc-pVXZ and aug-cc-pVXZ) (X=D, T, and Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(2d,2p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy (approximately 0.01 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the destined C6 value from molecular polarizability calculations. The slow convergence could indicate the inefficacy of using the MP2 calculations with Gaussian-type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the destined potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energy calculated using the PW91PW91 functional and the equilibrium bond length calculated using the PW91VP86 functional are close to the MP2 results at the basis set limit.     

Electronic excitation of sulfur-organic compounds – performance of time-dependent density functional theory

 To define the scope and limitations of the time-dependent density functional theory (TDDFT) method, spectral absorption data of a series of about 100 neutral or charged sulfur-organic compounds with up to 24 non-hydrogen atoms and up to four sulfur atoms were calculated in the near-UV, visible and IR regions. Although the theoretical vertical transition energies correspond only approximately to experimental absorption band maxima, the mean absolute deviation was calculated to be 0.21 eV (1600 cm⁻¹). The main absorption features of various compounds with monocoordinated or dicoordinated sulfur atoms are well reproduced. As far as possible TDDFT results were compared with those of semiempirical Zerner's intermediate neglect of differential overlap (ZINDO/S) and of Pariser–Parr–Pople (PPP) calculations. TDDFT also works well in cases where the semiempirical methods fail. Limitations of TDDFT were encountered with calculations of spectral absorptions of dye molecules. The “vinylene shift” of polymethine dyes is not reproduced by TDDFT. Whereas electronic excitation energies delocalized polar and betainic chromophores are reasonably well reproduced, excitation energies of charge-transfer-type and charge-resonance-type transitions of weakly interacting composite chromophores are significantly underestimated. Received: 30 October 2000 / Accepted: 29 November 2000 / Published online: 22 May 2001     

Vibrational spectroscopy investigation using ab initio and density functional theory on 3′-chloropropiophenone and 3′-nitropropiophenone

The FTIR and FT Raman spectra of 3′-chloropropiophenone and 3′-nitropropiophenone have been recorded in the regions 4000–400 and 3500–100 cm^?1 respectively. The optimized geometry, frequency and intensity of the vibrational bands of 3′-chloropropiophenone and 3′-nitropropiophenone were obtained by ab initio and DFT levels of theory with complete relaxation in the potential energy surface using 6-31G (d,p) basis set. A complete vibrational assignment aided by the theoretical harmonic frequency analysis has been proposed. The harmonic vibrational frequencies calculated have been compared with experimental FTIR and FT Raman spectra. The observed and the calculated frequencies are found to be in good agreement. The experimental spectra also coincide satisfactorily with those of theoretically constructed simulated spectrograms.     

The importance of Colle–Salvetti for computational density functional theory

The purpose of this presentation is to show the importance of the Colle–Salvetti (Theor Chim Acta 37:329, 1975) paper in the development of modern computational density functional theory. To do this we cover the following topics (1) the Bright Wilson understanding (2) the Kohn–Sham equations (3) local density exchange (4) the exchange-hole (5) generalised gradient approximation for exchange (Becke and Cohen) (6) left–right correlation and dynamic correlation (7) the development of the Lee–Yang–Parr dynamic correlation functional from the Colle–Salvetti paper (8) the early success of GGA DFT. Finally we observe that the the BLYP and OLYP exchange-correlation functionals are not semi-empirical; this may explain their great success.     

Concurrent cyclopropanation by carbenes and carbanions? A density functional theory study on the reaction pathways

The mechanisms for the addition reactions of phenylhalocarbenes and phenyldihalomethide carbanions with acrylonitrile (ACN) and trimethylethylene (TME) have been investigated using an ab initio BH and HLYP/6-31G (d, p) level of theory. Solvent effects on these reactions have been explored by calculations that included a polarizable continuum model (PCM) for the solvent (THF). These model calculations show that for the addition of phenylhalocarbenes, a TME species may readily undergo addition reactions with carbenes while ACN has a high-energy barrier to overcome. It was also found that phenyldihalomethide carbanions do not readily add to the electron-rich TME. The cyclopropane yields only appear to occur via addition of PhCBr to TME. However, the cyclopropanation proceeds not only via slow addition of phenylhalocarbenes to ACN but also forms through the stepwise reaction of phenyldihalomethide carbanions with ACN. Our calculation results are in good agreement with experimental observations (Moss, R.A.; Tian, J.-Z. J. Am. Chem. Soc. 2005, 127, 8960) that indicate that the cyclopropanation of phenylhalocarbenes and phenyldihalomethide carbanions with ACN are concurrent in THF.     

A density functional theory study on mechanisms of [4 + 2] annulation of enal with α-methylene cycloalkanone catalyzed by N-heterocyclic carbene

Density functional theory was used to study the reaction mechanisms of the N-heterocyclic carbene (NHC)-catalyzed [4 + 2] annulation reaction between enal and α-methylene cycloalkanone for the formation of tricyclic benzopyran-2-one. The simulations suggest that the energy-favorable catalytic cycle includes five steps: (a) the nucleophilic addition of enal by NHC catalyst; (b) [1, 2]-proton transfer for the formation of Breslow intermediate; (c) [1, 4]-proton transfer for the formation of enolate intermediate; (d) the [4 + 2] cycloaddition process, product formation; and (e) catalyst regeneration. The proton transfer process was particularly designed in two independent ways, the direct proton transfer and the Brønsted acid DBU∙H⁺-mediated proton transfer. Our study reveals that [1, 2]-proton transfer process is the rate-determining step with the accessible energy barrier, which agrees with the experimental observation. Further analysis of the global reaction index confirmed that NHC was primary used as a Lewis base during the reaction processes. The frontier molecular orbital (FMO) analysis indicates that the introduction of the NHC catalyst facilitates the reaction to occur due to the narrower energy gap of FMO.     

Synthesis of α-allyl-α-amino acids by means of a palladium catalyzed intramolecular rearrangement

Allylic esters of imines of amino acids provide α-alkylated-α-amino acids in good yields, on rearrangement by catalytic quantities of Pd(O) complexes.     

Are electrophilicity and electrofugality related concepts? A density functional theory study

It is proposed that the electrofugality of a fragment within a molecule is determined by its group nucleophilicity. The variation of electrofugality should be tightly related to the electron releasing ability of the substituent attached to the electrofuge moiety. This contribution closes the set of relationships between philicity and fugality quantities: while nucleofugality appears related to the group electrophilicity of the leaving group, electrofugality is related to the group nucleophilicity of the permanent group.     

Near-infrared-emitting phthalocyanines. A combined experimental and density functional theory study of the structural, optical, and photophysical properties of Pd(II) and Pt(II) α-butoxyphthalocyanines

The structural, optical, and photophysical properties of 1,4,8,11,15,18,22,25-octabutoxyphthalocyaninato-palladium(II), PdPc(OBu)(8), and the newly synthesized platinum analogue PtPc(OBu)(8) are investigated combining X-ray crystallography, static and transient absorption spectroscopy, and relativistic zeroth-order regular approximation (ZORA) Density Functional Theory (DFT)/Time Dependent DFT (TDDFT) calculations where spin-orbit coupling (SOC) effects are explicitly considered. The results are compared to those previously reported for NiPc(OBu)(8) (J. Phys. Chem. A 2005, 109, 2078) in an effort to highlight the effect of the central metal on the structural and photophysical properties of the group 10 transition metal octabutoxyphthalocyanines. Different from the nickel analogue, PdPc(OBu)(8) and PtPc(OBu)(8) show a modest and irregular saddling distortion of the macrocycle, but share with the first member of the group similar UV-vis spectra, with the deep red and intense Q-band absorption experiencing a blue shift down the group, as observed in virtually all tetrapyrrolic complexes of this triad. The blue shift of the Q-band along the MPc(OBu)(8) (M = Ni, Pd, Pt) series is interpreted on the basis of the metal-induced electronic structure changes. Besides the intense deep red absorption, the title complexes exhibit a distinct near-infrared (NIR) absorption due to a transition to the double-group 1E (π,π*) state, which is dominated by the lowest single-group (3)E (π,π*) state. Unlike NiPc(OBu)(8), which is nonluminescent, PdPc(OBu)(8) and PtPc(OBu)(8) show both deep red fluorescence emission and NIR phosphorescence emission. Transient absorption experiments and relativistic spin-orbit TDDFT calculations consistently indicate that fluorescence and phosphorescence emissions occur from the S(1)(π,π*) and T(1)(π,π*) states, respectively, the latter being directly populated from the former, and the triplet state decays directly to the S(0) surface (the triplet lifetime in deaerated benzene solution was 3.04 μs for Pd and 0.55 μs for Pt). Owing to their triplet properties, PdPc(OBu)(8) and PtPc(OBu)(8) have potential for use in photodynamic therapy (PDT) and are potential candidates for NIR light emitting diodes or NIR emitting probes.     

Effects of H-bonding and structural constituents on the acidity and potential "anticancer activity" of D-mandelonitrile-β-D-glucuronic acid by density functional theory

Research on Chemical Intermediates - Laetrile (D-mandelonitrile-β-D-glucuronic acid, 1COOH) is being promoted, more and more, as an alternate cancer treatment, in many countries. Here we show...     

Application of capillary affinity electrophoresis and density functional theory to the investigation of valinomycin–lithium complex

Capillary affinity electrophoresis (CAE) and quantum mechanical density functional theory (DFT) have been applied to the investigation of interactions of valinomycin (Val), a macrocyclic dodecadepsipeptide antibiotic ionophore, with lithium cation Li⁺. Firstly, from the dependence of effective electrophoretic mobility of Val on the Li⁺ ion concentration in the background electrolyte (BGE) (methanolic solution of 50 mM chloroacetic acid, 25 mM Tris, pH_MeOH 7.8, 0–40 mM LiCl), the apparent binding (stability) constant (K_b) of Val–Li⁺ complex in methanol was evaluated as log K_b = 1.50 ± 0.24. The employed CAE method include correction of the effective mobilities measured at ambient temperature, at different input power (Joule heating) and at variable ionic strength of the BGEs to the mobilities related to the reference temperature 25 °C and to the constant ionic strength 25 mM. Secondly, using DFT calculations, the most probable structures of the non-hydrated Val–Li⁺ and hydrated Val–Li⁺·3H₂O complex species were predicted.     

How are the ready and unready states of nickel-iron hydrogenase activated by H2? A density functional theory study

We have explored possible mechanisms for the formation of the catalytically active Ni(a)-S state of the enzyme, nickel iron hydrogenase, from the Ni*(r) (ready) or Ni*(u) (unready) state, by reaction with H(2), using density functional theory calculations with the BP86 functional in conjunction with a DZVP basis set. We find that for the reaction of the ready state, which is taken to have an -OH bridge, the rate determining step is the cleavage of H(2) at the Ni(3+) centre with a barrier of approximately 15 kcal mol(-1). We take the unready state to have a -OOH bridge, and find that reaction with H(2) to form the Ni(r)-S state can proceed by two possible routes. One such path has a number of steps involving electron transfer, which is consistent with experiment, as is the calculated barrier of approximately 19 kcal mol(-1). The alternative pathway, with a lower barrier, may not be rate determining. Overall, our predictions give barriers in line with experiment, and allow details of the mechanism to be explored which are inaccessible from experiment.