PubMed Health⌕ Search

Biomedical subjects

Yixiang Cao

Publications and source records attributed to Yixiang Cao.

4 recordsLinked to original sources

A localized orbital analysis of the thermochemical errors in hybrid density functional theory: achieving chemical accuracy via a simple empirical correction scheme.

This paper describes an empirical localized orbital correction model which improves the accuracy of density functional theory (DFT) methods for the prediction of thermochemical properties for molecules of first and second row elements. The B3LYP localized orbital correction version of the model improves B3LYP DFT atomization energy calculations on the G3 data set of 222 molecules from a mean absolute deviation (MAD) from experiment of 4.8 to 0.8 kcal/mol. The almost complete elimination of large outliers and the substantial reduction in MAD yield overall results comparable to the G3 wave-function-based method; furthermore, the new model has zero additional computational cost beyond standard DFT calculations. The following four classes of correction parameters are applied to a molecule based on standard valence bond assignments: corrections to atoms, corrections to individual bonds, corrections for neighboring bonds of a given bond, and radical environmental corrections. Although the model is heuristic and is based on a 22 parameter multiple linear regression to experimental errors, each of the parameters is justified on physical grounds, and each provides insight into the fundamental limitations of DFT, most importantly the failure of current DFT methods to accurately account for nondynamical electron correlation.

Fluorine↗

Nuclear-magnetic-resonance shielding constants calculated by pseudospectral methods.

We have developed an algorithm based upon pseudospectral (PS) ab initio electronic structure methods for evaluating nuclear magnetic shielding constants using gauge-including atomic orbitals (GIAOs) in the spin-restricted and spin-unrestricted formalisms of Hartree-Fock (HF) theory and density-functional theory (DFT). The nuclear magnetic shielding constants for both 1H and 13C calculated using PS methodology for 21 small molecules have absolute mean errors of less than 0.3 ppm in comparison with analytic integral results. CPU timing comparisons between PS methods and conventional methods carried out for seven large molecules ranging from 510 to 1285 basis functions demonstrate that the PS methods are an order of magnitude more efficient than the conventional methods. PS-HF was between 9 and 26 times faster than conventional integral technology, and PS-DFT (Becke three-parameter Lee-Yang-Parr) was between 6 and 21 times faster.

Journal Article↗

Molecular (hyper)polarizabilities computed by pseudospectral methods.

We have developed algorithms based on pseudospectral (PS) ab initio electronic structure methods for solving the first- and second-order Hartree-Fock/Kohn-Sham equations and evaluating molecular polarizabilities and first- and second-order hyperpolarizabilities in the spin-restricted and spin-unrestricted formalisms at the Hartree-Fock (HF) and density functional theory (DFT) levels. We carry out calculations on 50 small molecules to test the accuracy of the PS approach. Our results demonstrate that the molecular polarizability alpha computed by the PS method is essentially identical to the value obtained from conventional methods for both HF and DFT calculations, while the first-order hyperpolarizability beta and second-order hyperpolarizability gamma have mean unsigned percentage differences of 1.26% and 0.62% (HF) and 0.78% and 0.65% (DFT), respectively. We also present CPU timing comparisons between the PS and conventional methods at the 6-31 G(**) level for 14 molecules having 185 to 1185 basis functions. The timing results show that the PS method is 25 (PS-HF) and 13 (PS-DFT) times faster than the conventional method for a system with 500 basis functions. The PS methods are found scale as N(2.70) (PS-HF) and N(2.40) (PS-DFT), while the conventional methods scale as N(2.93) (PRISM-HF) and N(2.87) (PRISM-DFT), where N is the number of basis functions.

Journal Article↗

Integrated Modeling Program, Applied Chemical Theory (IMPACT).

We provide an overview of the IMPACT molecular mechanics program with an emphasis on recent developments and a description of its current functionality. With respect to core molecular mechanics technologies we include a status report for the fixed charge and polarizable force fields that can be used with the program and illustrate how the force fields, when used together with new atom typing and parameter assignment modules, have greatly expanded the coverage of organic compounds and medicinally relevant ligands. As we discuss in this review, explicit solvent simulations have been used to guide our design of implicit solvent models based on the generalized Born framework and a novel nonpolar estimator that have recently been incorporated into the program. With IMPACT it is possible to use several different advanced conformational sampling algorithms based on combining features of molecular dynamics and Monte Carlo simulations. The program includes two specialized molecular mechanics modules: Glide, a high-throughput docking program, and QSite, a mixed quantum mechanics/molecular mechanics module. These modules employ the IMPACT infrastructure as a starting point for the construction of the protein model and assignment of molecular mechanics parameters, but have then been developed to meet specialized objectives with respect to sampling and the energy function.

Computer Simulation↗