Coupled‐cluster reaction barriers of : An application of the coupled‐cluster//Kohn–Sham density functional theory model chemistry

In this work, we report a theoretical investigation concerning the use of the popular coupled‐cluster//Kohn‐Sham density functional theory (CC//KS‐DFT) model chemistry, here applied to study the entrance channel of the reaction, namely by comparing CC//KS‐DFT calculations with KS‐DFT, MRPT2//CASSCF, and CC//CASSCF results from our previous investigations. This was done by performing single point energy calculations employing several coupled cluster methods and using KS‐DFT geometries optimized with six different functionals, while conducting a detailed analysis of the barrier heights and topological features of the curves and surfaces here obtained. The quality of this model chemistry is critically discussed in the context of the title reaction and also in a wider context. © 2013 Wiley Periodicals, Inc.     

Interaction of copernicium with gold: Assessment of applicability of simple density functional theories

Interactions of Cn (element 112) atom with small Au clusters are studied using accurate ab initio scalar relativistic coupled cluster method for correlation treatment and two‐component relativistic density functional theory (RDFT) to take account of spin‐dependent relativistic effects. The results demonstrate the failure of RDFT with simple generalized‐gradient and hybrid functionals in describing Cn–Au bonds in complex systems. © 2013 Wiley Periodicals, Inc.     

Energy landscapes of nucleophilic substitution reactions: a comparison of density functional theory and coupled cluster methods

We have carried out a detailed evaluation of the performance of all classes of density functional theory (DFT) for describing the potential energy surface (PES) of a wide range of nucleophilic substitution (SN2) reactions involving, amongst others, nucleophilic attack at carbon, nitrogen, silicon, and sulfur. In particular, we investigate the ability of the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA as well as hybrid DFT to reproduce high-level coupled cluster (CCSD(T)) benchmarks that are close to the basis set limit. The most accurate GGA, meta-GGA, and hybrid functionals yield mean absolute deviations of about 2 kcal/mol relative to the coupled cluster data, for reactant complexation, central barriers, overall barriers as well as reaction energies. For the three nonlocal DFT classes, the best functionals are found to be OPBE (GGA), OLAP3 (meta-GGA), and mPBE0KCIS (hybrid DFT). The popular B3LYP functional is not bad but performs significantly worse than the best GGA functionals. Furthermore, we have compared the geometries from several density functionals with the reference CCSD(T) data. The same GGA functionals that perform best for the energies (OPBE, OLYP), also perform best for the geometries with average absolute deviations in bond lengths of 0.06 A and 0.6 degrees, even better than the best meta-GGA and hybrid functionals. In view of the reduced computational effort of GGAs with respect to meta-GGAs and hybrid functionals, let alone coupled cluster, we recommend the use of accurate GGAs such as OPBE or OLYP for the study of SN2 reactions.     

Reciprocal adjustment of approximate coupled cluster and configuration interaction approaches

The linear version of the externally corrected coupled cluster method with singles and doubles (ecLCCSD), the recently proposed coupled cluster corrections to the multireference configuration interaction (ccMRCI) energies, and the so‐called self‐consistent, size‐consistent [(SC)²] approaches, which are designed to correct for the dynamic correlation effects and the size inconsistency of the MRCI energies, are analyzed and compared using several illustrative examples, including the dissociation of a triple‐zeta (TZ) model of the N₂ molecule. It is emphasized that the exponential cluster ansatz for the wave function is the basis of all these approaches, and appropriate cluster analysis of the MRCI wave function is the key step for both ecLCCSD and ccMRCI. The contributions from the orthogonal complement of the MRCI space, which can be generated by relying on such a cluster analysis, are responsible for a substantial part of the missing correlation energy. The ecLCCSD approach seems to represent a particularly attractive alternative to other highly accurate methods for the calculation of the ground‐state energy in the presence of quasidegeneracy, both due to its efficiency and affordability. It may in fact be regarded as a simple alternative to the iterative reduced multireference (RMR) CCSD method. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 693–703, 2000     

Cyanide adsorbed on coinage metal electrodes: A relativistic density functional investigation

The adsorptive properties of cyanide (CN⁻) on coinage metal (M) electrodes (M=Cu, Ag, Au) have been investigated using a relativistic density functional method. The way to model the electrochemical potential applied to the electrodes is to consider the systems in the presence of a perturbative external field F. The field-perturbative approach is proven to be a suitable method in interpreting the observed spectral shifts with electrode potential. The calculated potential-dependent shifts of ω_{M(SINGLE BOND)CN} and ω_{C(SINGLE BOND)M} are similar for the three metals, in agreement with experiment observations. The relativistic effects are required to account for the similarity in the frequency shifts of ω_{M(SINGLE BOND)CN}. The calculated vibrational tuning rates dω_{C(SINGLE BOND)N}/dF are 6.61×10⁻⁷, 6.61×10⁻⁷, and 5.64×10⁻⁷ cm⁻¹/(V/cm) for M=Cu, Ag, and Au, respectively. The coupling of the M(SINGLE BOND)CN and C(SINGLE BOND)N internal modes contributes significantly (about 25%) to the size of the frequency shifts Δω_{C(SINGLE BOND)N} of the ligand. The effect of electric fields on the metal(SINGLE BOND)CN⁻ bonding is also investigated. It is shown that changes in the magnitude of CN⁻ to the metal donation and M(SINGLE BOND)CN bond strength occur under the influence of the electric field. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 175–185, 1998     

Spin‐orbit ab initio and density functional theory investigation of bismuth monoboronyl,BiBO

The molecular properties of bismuth monoboronyl, BiBO, were investigated using high‐level ab initio and density functional theory calculations by including the effect of spin‐orbit coupling (SOC). SOC does not cause any change in the Bi? B bond length of BiBO, by contrast it causes significant elongation of the Bi? B bond of BiBO^?, by ～0.03 Å. The Bi? B bond length of BiBO^? that is calculated by considering SOC is almost identical to that of BiBO; this result is consistent with a recent experimental study. The term values of excited states of BiBO calculated by including SOC are in good agreement with the experimental results. One excited state which was not assigned in the previous experimental study is the Ω = 0⁺ state generated by strong SOC. In the theoretical calculations on molecules containing 6p‐block elements, including SOC is crucial for obtaining results that are consistent with the corresponding experimental results.     

Novel method of self‐interaction corrections in density functional calculations

It is demonstrated that the commonly applied self‐interaction correction (SIC) used in density functional theory does not remove all self‐interaction. We present as an alternative a novel method that, by construction, is totally free from self‐interaction. The method has the correct asymptotic 1/r dependence. We apply the new theory to localized f electrons in praseodymium and compare with the old version of SIC, the local density approximation (LDA) and with an atomic Hartree–Fock calculation. The results show a lowering of the f level, a contraction of the f electron cloud and a lowering of the total energy by 13 eV per 4 f electron compared to LDA. The equilibrium volume of the new SIC method is close to the ones given by LDA and the older SIC method and is in good agreement with experiment. The experimental cohesive energy is in better agreement using the new SIC method, both compared to LDA and another SIC method. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 247–252, 2001     

Reactivation pathway of the hydrogenase H‐cluster: Density functional theory study

This work puts forth a reaction pathway for the reactivation of exogenous ligand inhibited H‐cluster, the active site of Fe‐only hydrogenases. The H‐cluster is a dimetal complex, Fe–Fe, with the metal centers bridged by di(thiomethyl)amine. Exogenous ligands, H₂O, and OH^?, are bound to the distal iron (Fe_d). Density functional theory (DFT) calculations on the native and ruthenium‐modified H‐cluster have been performed using the B3LYP functional with 6‐31+G** and 6‐311+G** basis sets. We have ascertained that there is a thermodynamically favorable pathway for the reactivation of the OH^? inhibited H‐cluster, which proceeds by an initial protonation of the Fe_d–OH^? complex. The proposed reaction pathway has all its intermediate reactions ensue exothermically. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

金属Ir4簇催化乙烯加氢反应势能面的理论研究

应用密度泛函理论(DFT)对金属Ir4簇催化乙烯加氢反应的反应机理进行了详尽的理论研究．在B3LYP／ECP[C,H:6-311G(d)和6-31G(d);Ir:LANL2DZ]理论水平下优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型,并且用组态相互作用CCSD／ECP[C,H:6-311G(d,p);Ir:LANL2DZ]计算了各驻点的单点能,构建了该反应的基态势能面．为了验证过渡态的真实性,在B3LYP／ECP理论水平下做了内禀反应坐标(IRC)计算和频率分析．计算结果表明:金属Ir4簇催化乙烯加氢反应为双通道(a和b)反应,经过多个反应步骤完成;通道a:R→TSR-1→I1→TS1-2→I2→TS2-3→I3→TS3-P→P为较为可行的反应通道．     

Theoretical determination of lowest structures of all‐metal aromatic clusters M4L2 (M = Al,Ga, In,Tl; L = Li,Na, K,Rb, Cs)

The researches of all‐metal aromatic clusters have been a thermic theme in inorganic aromaticity domain both experimentally and theoretically since the Al₄L^? (L = Li, Na, Cu) clusters were created by laser vaporization. In systemic determination of the lowest structures of 20 gaseous all‐metal aromatic clusters M₄L₂ (M = Al, Ga, In, Tl; L = Li, Na, K, Rb, Cs), the isomer energy differences of four low‐lying structures of each cluster were evaluated at high‐quality quantum chemistry levels. Single point calculations at the coupled cluster level were performed at geometries optimized at the MP2, B3LYP, and B3PW91 levels, and harmonic frequency calculations and zero point energy corrections were implemented following optimizations at the B3LYP and B3PW91 levels. In addition to Li‐ and Na‐containing species, theoretical investigations came down to those new clusters including K, Rb, and Cs. For many clusters, the most convincing theoretical evidences indicate that the lowest structures are a square bipyramidal isomer rather than an edge‐caped square pyramidal species. A few discrepancies were addressed at the MP2, B3LYP, and B3PW91 levels in comparison with the coupled cluster results. These findings are significant because some clusters were generated by laser vaporization and served as theoretical prototypes to test the new means for assessing aromaticity. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

IR and Raman spectra,density functional computations of vibrational spectrum,molecular geometry,atomic charges,and some molecular properties of 3‐aminobenzonitrile molecule

Geometry, vibrational wavenumbers, and several thermodynamic parameters have been calculated using ab initio quantum chemical methods for the 3‐aminobenzonitrile molecule for the first time. The results were compared with experimental values. With the help of specific scaling procedures, the observed vibrational wavenumbers were analyzed and assigned to different normal modes of the molecule. In general, the error obtained was very low. Using potential energy distribution (PED), the contributions of the different modes to each wavenumber were determined. Other general conclusions were also deduced. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

DFT study of the carbon‐ and nitrogen‐pivot lariat crown ethers and their complexes with alkali metal cations: Na+, K+

In this work, a quantum mechanical research of five lariat crown ethers(LCEs), 2‐methoxy‐15‐crown‐5( A ), N‐methoxy‐4‐aza‐15‐crown‐5( B ), N‐methoxy‐4‐aza‐18‐crown‐6( C ), N‐methoxyethyl‐4‐aza‐18‐crown‐6( D ), N,N′‐bis(2‐metho xyethyl)‐4,13‐diaza‐18‐crown‐6( E ), which are based on either 15‐crown‐5 or 18‐crown‐6 frameworks and contain various pendant arms extending from either carbon or nitrogen atoms on the crown frameworks, had been done using density functional theory with B3LYP/6‐31G* method to obtain the electronic and geometrical structures of the LCEs and their complexes with alkali metal ions: Na⁺ and K⁺. The nucleophilicity of LCEs had been investigated by the Fukui functions. For complexes, the match between the cation and cavity size, the status of interaction between alkali metal ions and donor atoms in the LCEs, and the sidearm effect of the LCEs had been analyzed through the other calculated parameters, such as, highest occupied molecular orbital energy, lowest unoccupied molecular orbital energy, and energy gaps. In addition, the enthalpies of complexation reaction had been studied by the calculated thermodynamic data (298 K). The calculated results are all in a good agreement with the experimental data for the complexes. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

An open‐source framework for analyzing N‐electron dynamics. II. Hybrid density functional theory/configuration interaction methodology

In this contribution, we extend our framework for analyzing and visualizing correlated many‐electron dynamics to non‐variational, highly scalable electronic structure method. Specifically, an explicitly time‐dependent electronic wave packet is written as a linear combination of N‐electron wave functions at the configuration interaction singles (CIS) level, which are obtained from a reference time‐dependent density functional theory (TDDFT) calculation. The procedure is implemented in the open‐source Python program det CI@ORBKIT, which extends the capabilities of our recently published post‐processing toolbox (Hermann et al., J. Comput. Chem. 2016, 37, 1511). From the output of standard quantum chemistry packages using atom‐centered Gaussian‐type basis functions, the framework exploits the multideterminental structure of the hybrid TDDFT/CIS wave packet to compute fundamental one‐electron quantities such as difference electronic densities, transient electronic flux densities, and transition dipole moments. The hybrid scheme is benchmarked against wave function data for the laser‐driven state selective excitation in LiH. It is shown that all features of the electron dynamics are in good quantitative agreement with the higher‐level method provided a judicious choice of functional is made. Broadband excitation of a medium‐sized organic chromophore further demonstrates the scalability of the method. In addition, the time‐dependent flux densities unravel the mechanistic details of the simulated charge migration process at a glance. © 2017 Wiley Periodicals, Inc.     

A solvation‐free‐energy functional: A reference‐modified density functional formulation

The three‐dimensional reference interaction site model (3D‐RISM) theory, which is one of the most applicable integral equation theories for molecular liquids, overestimates the absolute values of solvation‐free‐energy (SFE) for large solute molecules in water. To improve the free‐energy density functional for the SFE of solute molecules, we propose a reference‐modified density functional theory (RMDFT) that is a general theoretical approach to construct the free‐energy density functional systematically. In the RMDFT formulation, hard‐sphere (HS) fluids are introduced as the reference system instead of an ideal polyatomic molecular gas, which has been regarded as the appropriate reference system of the interaction‐site‐model density functional theory for polyatomic molecular fluids. We show that using RMDFT with a reference HS system can significantly improve the absolute values of the SFE for a set of neutral amino acid side‐chain analogues as well as for 504 small organic molecules. © 2015 Wiley Periodicals, Inc.     

Torsional effects on excitation energies of thiophene derivatives induced by beta-substituents: comparison between time-dependent density functional theory and approximated coupled cluster approaches

The influence of methyl or phenyl substitution in beta-position of dioxygenated terthiophene and diphenylthiophene on the optical properties is investigated by first-principles calculations. We compare the approximated singles and doubles coupled cluster (CC2) approach with time-dependent density functional theory methods. CC2 reproduces experimental excitation energies with an accuracy of 0.1 eV. We find that the different substituents modify the inter-ring torsional angle which in turn strongly influences the excitation energies. The steric contribution to the excitation energies have been separated from the total substituent effects.     

A density functional theory study on diazotization and nitration of 3,5‐diamino‐1,2,4‐triazole

The reaction mechanism, thermodynamic and kinetic properties for diazotization and nitration of 3,5‐diamino‐1,2,4‐triazole were studied by a density functional theory. The geometries of the reactants, transition states, and intermediates were optimized at the B3LYP/6‐31G (d, p) level. Vibrational analysis was carried out to confirm the transition state structures, and the intrinsic reaction coordinate (IRC) method was used to explore the minimum energy path. The single‐point energies of all stagnation points were further calculated at the B3LYP (MP2)/6‐311+G (2d, p) level. The statistical thermodynamic method and Eyring transition state theory with Wigner correction were used to study the thermodynamic and kinetic characters of all reactions within 0–25°C. Two reaction channels are computed, including the diazotization and nitration of 3‐NH₂ or 5‐NH₂, and there are six steps in each channel. The reaction rate in each step is increased with temperature. The last step in each channel is the slowest step. The first, second, and fifth steps are exothermic reactions, and are favored at lower temperature in the thermodynamics. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

A screened hybrid density functional study on energetic complexes: Metal carbohydrazide nitrates

The molecular geometry, electronic structure and thermochemistry of a series of metal carbohydrazide nitrates were investigated using the Heyd–Scuseria–Ernzerhof (HSE) screened hybrid density functional. The results show that Ca, Sr, and Ba complexes have additional coordinated oxygen atoms from the nitrate ion, which differ obviously from Cu, Ni, Co, and Mg complexes in terms of the geometric structure. Detailed NBO analyses clearly indicate that the metal–ligand interactions in Cu, Ni, and Co complexes are covalent, whereas those of Mg, Ca, Sr, and Ba complexes are ionic in nature. Furthermore, the donor–acceptor interactions result in a reduction of occupancies of σ_{C? O} and σ_{N? H} orbitals. Consequently, the bond lengths increase and the bond orders decrease. Finally, the calculated heats of formation predict that the ionic alkaline‐earth metal carbohydrazide nitrates are more stable than the covalent transition metal carbohydrazide nitrates. It agrees well with the available experimental thermal stabilities, indicating that the metal–ligand bonding character plays an important role in the stabilities of these energetic complexes. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

An assessment of pure,hybrid, meta,and hybrid‐meta GGA density functional theory methods for open‐shell systems: The case of the nonheme iron enzyme 8R–LOX

The performance of a range density functional theory functionals combined in a quantum mechanical (QM)/molecular mechanical (MM) approach was investigated in their ability to reliably provide geometries, electronic distributions, and relative energies of a multicentered open‐shell mechanistic intermediate in the mechanism 8R–Lipoxygenase. With the use of large QM/MM active site chemical models, the smallest average differences in geometries between the catalytically relevant quartet and sextet complexes were obtained with the B3LYP^* functional. Moreover, in the case of the relative energies between ⁴II and ⁶II , the use of the B3LYP^* functional provided a difference of 0.0 kcal mol^–1. However, B3LYP^± and B3LYP also predicted differences in energies of less than 1 kcal mol^–1. In the case of describing the electronic distribution (i.e., spin density), the B3LYP^*, B3LYP, or M06‐L functionals appeared to be the most suitable. Overall, the results obtained suggest that for systems with multiple centers having unpaired electrons, the B3LYP^* appears most well rounded to provide reliable geometries, electronic structures, and relative energies. © 2012 Wiley Periodicals, Inc.     

Interaction between RNA segment (adenine‐uracil) and model of protein unit (formamide): A density‐functional theory study

Various possible structures of adenine‐uracil‐formamide hydrogen‐bond complexes were optimized at 6‐311++G(d,p) level, and the binding energies of these complexes were also calculated at DFT B3LYP/6‐311++G(d,p) level. Eight stable cyclic structures being involved in the interaction are found on the potential energy surface. By analyzing the structure, NPA charge and interaction energy of complexes, we obtain the most stable geometry structure. The results show that the interactions between formamide and adenine‐uracil (A‐U) base pair affect the stabilities of the base pairs. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010