改进的相对论量子化学计算ZORA方法

提出一种改进的ZORA(Zeroth-OrderRegularApproximationtotheDiracEquation)方法,其单电子方程为：[σ·p(c^2/2c^2-V~0)σ·p+V(r)]ψ=εψ。式中V~0为空间限域的势能函数：V~0(r=~A∑V~0^A(r~A),r~A=|r-R~A|,V~0^A(r~A)=V^A(r~A){1+exp[α(r~A-r~0^A]}^-^1。其中A表示分子的某个组成原子,R~A为A原子的位置矢量,V^A(r~A)为自由A原子的势函数,α和r~0^A为参数。改进的ZORA方法具有原来方法的所有优点,避免了原有ZORA方法因不满足标度变换不变性带来的缺陷,而且计算过程简单。具体计算表明,通过适当选择参数α和r~0^A,用本研究提出来的方法,在计算分子几何结构和键合能时,基本上消除了ZORA方法由于标度变换依赖性产生的误差。     

Fully relativistic correlated benchmark results for uranyl and a critical look at relativistic effective core potentials for uranium

 A complete set of fully-relativistic benchmark results for the bond lengths and vibrational frequencies of uranyl at various levels of correlation treatment are presented. It is shown that the relativistic and correlation effects are of the same magnitude and should be treated on an equal footing. Results of uranyl calculations using various relativistic effective core potentials (RECP) are presented. Appropriate choices for RECPs for use in density functional theory (DFT) calculations in the local density approximation (LDA) and with the gradient-corrected B3LYP exchange-correlation functional are discussed. The conclusion is reached that small-core RECPs need to be used and that the best results as compared to the benchmark values are obtained by using a DFT functional that includes some fraction of Hartree-Fock Exchange. Received: 18 May 2001 / Accepted: 25 July 2001 / Published online: 30 October 2001     

A singularity excluded approximate expansion scheme in relativistic density functional theory

 A singularity excluded approximate expansion (SEAX) scheme, which can be considered as one between Breit-Pauli expansion and RA expansion schemes, is proposed to expand the total energy of 4-component relativistic density functional theory. The one-electron equation can be derived variationally from the approximate total energy expression. The Hamiltonian of the one-electron equation is bounded from below and can be dealt with variationally, and the gauge dependency error in the ZORA method is essentially eliminated. It is easier to solve the SEAX equation than the IORA equation. The results related to the valence orbitals by solving the scalar SEAX equation agree very well with those by the scalar ZORA ESA method, and the results related to the inner-shell electrons of heavy elements by the two component SEAX calculations agree quite well with those by the 4-component relativistic density functional calculations. Received: 20 February 2002 / Accepted: 29 April 2002 / Published online: 8 July 2002     

Realistic sampling of amino acid geometries for a multipolar polarizable force field

The Quantum Chemical Topological Force Field (QCTFF) uses the machine learning method kriging to map atomic multipole moments to the coordinates of all atoms in the molecular system. It is important that kriging operates on relevant and realistic training sets of molecular geometries. Therefore, we sampled single amino acid geometries directly from protein crystal structures stored in the Protein Databank (PDB). This sampling enhances the conformational realism (in terms of dihedral angles) of the training geometries. However, these geometries can be fraught with inaccurate bond lengths and valence angles due to artefacts of the refinement process of the X‐ray diffraction patterns, combined with experimentally invisible hydrogen atoms. This is why we developed a hybrid PDB/nonstationary normal modes (NM) sampling approach called PDB/NM. This method is superior over standard NM sampling, which captures only geometries optimized from the stationary points of single amino acids in the gas phase. Indeed, PDB/NM combines the sampling of relevant dihedral angles with chemically correct local geometries. Geometries sampled using PDB/NM were used to build kriging models for alanine and lysine, and their prediction accuracy was compared to models built from geometries sampled from three other sampling approaches. Bond length variation, as opposed to variation in dihedral angles, puts pressure on prediction accuracy, potentially lowering it. Hence, the larger coverage of dihedral angles of the PDB/NM method does not deteriorate the predictive accuracy of kriging models, compared to the NM sampling around local energetic minima used so far in the development of QCTFF. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

A comparison of quantum chemical models for calculating NMR shielding parameters in peptides: mixed basis set and ONIOM methods combined with a complete basis set extrapolation

This article compares several quantum mechanical approaches to the computation of chemical shielding tensors in peptide fragments. First, we describe the effects of basis set quality up to the complete basis set (CBS) limit and level of theory (HF, MP2, and DFT) for four different atoms in trans N-methylacetamide. For both isotropic shielding and shielding anisotropy, the MP2 results in the CBS limit show the best agreement with experiment. The HF values show quite a different tendency to MP2, and even in the CBS limit they are far from experiment for not only the isotropic shielding of carbonyl carbon but also most shielding anisotropies. In most cases, the DFT values differ systematically from MP2, and small basis-set (double- or triple-zeta) results are often fortuitously in better agreement with the experiment than the CBS ones. Second, we compare the mixed basis set and ONIOM methods, combined with CBS extrapolation, for chemical shielding calculations at a DFT level using various model peptides. From the results, it is shown that the mixed basis set method provides better results than ONIOM, compared to CBS calculations using the nonpartitioned full systems. The information studied here will be useful in guiding the selection of proper quantum chemical models, which are in a tradeoff between accuracy and cost, for shielding studies of peptides and proteins.     

Four‐component relativistic DFT calculations of 77Se NMR chemical shifts: A gateway to a reliable computational scheme for the medium‐sized organoselenium molecules

A versatile high‐accuracy computational scheme for the ⁷⁷Se nuclear magnetic resonance (NMR) chemical shifts of the medium‐sized organoselenium compounds is suggested within a framework of a full four‐component relativistic density functional theory (DFT). The main accuracy factors (DFT functionals, relativistic geometry, vibrational corrections, and solvent effects) are addressed. The best result is achieved with NMR‐oriented KT2 functional of Keal–Tozer characterized by a fairly small error of only 30 ppm for the span of about 1700 ppm (<2%). © 2015 Wiley Periodicals, Inc.     

DFT/NMR integrated approach: a valid support to the total synthesis of chiral molecules

The analysis of the configuration of kedarcidin chromophore and palau'amine through quantum chemical calculation of Js and chemical shifts suggests a fast and convenient quantum chemical approach that can be applied prior to proceeding to the total synthesis of complex natural compounds in order to avoid loss of time and resources employed in the total synthesis of wrong diastereoisomers.     

WebMTA: a web-interface for ab initio geometry optimization of large molecules using molecular tailoring approach

A web-interface for geometry optimization of large molecules using a linear scaling method, i.e., cardinality guided molecular tailoring approach (CG-MTA), is presented. CG-MTA is a cut-and-stitch, fragmentation-based method developed in our laboratory, for linear scaling of conventional ab initio techniques. This interface provides limited access to CG-MTA-enabled GAMESS. It can be used to obtain fragmentation schemes for a given spatially extended molecule depending on the maximum allowed fragment size and minimum cut radius values provided by the user. Currently, we support submission of single point or geometry optimization jobs at Hartree-Fock and density functional theory levels of theory for systems containing between 80 to 200 first row atoms and comprising up to 1000 basis functions. The graphical user interface is built using HTML and Python at the back end. The back end farms out the jobs on an in-house Linux-based cluster running on Pentium-4 Class or higher machines using an @Home-based parallelization scheme (http://chem.unipune.ernet.in/ approximately tcg/mtaweb/).