PubMed Health⌕ Search

Biomedical subjects

Lingchun Song

Publications and source records attributed to Lingchun Song.

7 recordsLinked to original sources

XMVB: a program for ab initio nonorthogonal valence bond computations.

An ab initio nonorthogonal valence bond program, called XMVB, is described in this article. The XMVB package uses Heitler-London-Slater-Pauling (HLSP) functions as state functions, and calculations can be performed with either all independent state functions for a molecule or preferably a few selected important state functions. Both our proposed paired-permanent-determinant approach and conventional Slater determinant expansion algorithm are implemented for the evaluation of the Hamiltonian and overlap matrix elements among VB functions. XMVB contains the capabilities of valence bond self-consistent field (VBSCF), breathing orbital valence bond (BOVB), and valence bond configuration interaction (VBCI) computations. The VB orbitals, used to construct VB functions, can be defined flexibly in the calculations depending on particular applications and focused problems, and they may be strictly localized, delocalized, or bonded-distorted (semidelocalized). The parallel version of XMVB based on MPI (Message Passing Interface) is also available.

Journal Article↗

Valence bond calculations of hydrogen transfer reactions: a general predictive pattern derived from theory.

Hydrogen abstraction reactions of the type X(*) + H-H' --> X-H + H'(*) (X = F, Cl, Br, I) are studied by ab initio valence bond methods and the VB state correlation diagram (VBSCD) model. The reaction barriers and VB parameters of the VBSCD are computed by using the breathing orbital valence bond and valence bond configuration interaction methods. The combination of the VBSCD model and semiempirical VB theory leads to analytical expressions for the barriers and other VB quantities that match the ab initio VB calculations fairly well. The barriers are influenced by the endo- or exothermicity of the reaction, but the fundamental factor of the barrier is the average singlet-triplet gap of the bonds that are broken or formed in the reactions. Some further approximations lead to a simple formula that expresses the barrier for nonidentity and identity hydrogen abstraction reactions as a function of the bond strengths of reactants and products. The semiempirical expressions are shown to be useful not only for the model reactions that are studied in this work, but also for other nonidentity and identity hydrogen abstraction reactions that have been studied in previous articles.

Journal Article↗

The ground and excited states of polyenyl radicals C2n-1H2n + 1 (n = 2-13): a valence bond study.

The semiempirical valence bond (VB) method, VBDFT(s), is applied to the ground states and the covalent excited states of polyenyl radicals C2n - 1H2n + 1 (n = 2-13). The method uses a single scalable parameter with a value that carries over from the study of the covalent excited states of polyenes (W. Wu, D. Danovich, A. Shurki, S. Shaik, J. Phys. Chem. A, 2000, 104, 8744). Whenever comparison is possible, the VB excitation energies are found to be in good accord with sophisticated molecular orbital (MO)-based methods like CASPT2. The symmetry-adapted Rumer structures are used to discuss the state-symmetry and VB constitution of the ground and excited states, and the expansion to VB determinants is used to gain insight on spin density patterns. The theory helps to understand in a coherent and lucid manner the properties of polyenyl radicals, such as the makeup of the various states, their geometries and energies, and the distribution of the unpaired electrons (the neutral solitons).

Journal Article↗

Charge transfer in the electron donor-acceptor complex BH3NH3.

As a simple yet strongly binding electron donor-acceptor (EDA) complex, BH(3)NH(3) serves as a good example to study the electron pair donor-acceptor complexes. We employed both the ab initio valence bond (VB) and block-localized wave function (BLW) methods to explore the electron transfer from NH(3) to BH(3). Conventionally, EDA complexes have been described by two diabatic states: one neutral state and one ionic charge-transferred state. Ab initio VB self-consistent field (VBSCF) computations generate the energy profiles of the two diabatic states together with the adiabatic (ground) state. Our calculations evidently demonstrated that the electron transfer between NH(3) and BH(3) falls in the abnormal regime where the reorganization energy is less than the exoergicity of the reaction. The nature of the NH(3)-BH(3) interaction is probed by an energy decomposition scheme based on the BLW method. We found that the variation of the charge-transfer energy with the donor-acceptor distance is insensitive to the computation levels and basis sets, but the estimation of the amount of electron transferred heavily depends on the population analysis procedures. The recent resurgence of interest in the nature of the rotation barrier in ethane prompted us to analyze the conformational change of BH(3)NH(3), which is an isoelectronic system with ethane. We found that the preference of the staggered structure over the eclipsed structure of BH(3)NH(3) is dominated by the Pauli exchange repulsion.

Journal Article↗

A practical valence bond method: a configuration interaction method approach with perturbation theoretic facility.

The previously developed valence bond configuration interaction (VBCI) method (Wu, W.; Song, L.; Cao, Z.; Zhang, Q.; Shaik, S., J. Phys. Chem. A, 2002, 105, 2721) that borrows the general CI philosophy of the MO theory, is further extended in this article, and its methodological features are improved, resulting in three accurate and cost-effective procedures: (a) the effect of quadruplet excitation is incorporated using the Davidson correction, such that the new procedure reduces size consistency problems, with due improvement in the quality of the computational results. (b) A cost-effective procedure, named VBCI(D, S), is introduced. It includes doubly excited structures for active electrons and singly excited structures for inactive pairs. The computational results of VBCI(D, S) match those of VBCISD with much less computational effort than VBCISD. (c) Finally, a second-order perturbation theory is utilized as a means of configuration selection, and lead to considerable reduction of the computational cost, with little or no loss in accuracy. Applications of the new procedures to bond energies and barriers of chemical reactions are presented and discussed.

Journal Article↗

An accurate barrier for the hydrogen exchange reaction from valence bond theory: is this theory coming of age?

One of the landmark achievements of quantum chemistry, specifically of MO-based methods that include electron correlation, was the precise calculation of the barrier for the hydrogen-exchange reaction (B. Liu, J. Chem. Phys. 1973, 58, 1925; P. Siegbahn, B. Liu, J. Chem. Phys. 1978, 68, 2457). This paper reports an accurate calculation of this barrier by two recently developed VB methods that use only the eight classical VB structures. To our knowledge, the present work is the first accurate ab initio VB barrier that matches an experimental value. Along with the accurate barrier, the VB method provides accurate bond energies and diabatic quantities that enable the barrier height to be analyzed by the VB state correlation diagram approach, VBSCD (S. Shaik, A. Shurki, Angew. Chem. 1999, 111, 616; Angew. Chem. Int. Ed. Engl. 1999, 38, 586). This is a proof of principal that VB theory with appropriate account of dynamic electron correlation can achieve quantitative accuracy of reaction barriers, and still retain a compact and interpretable wave function. A sample of S(N)2 barriers and dihalogen bonding energies, which are close to CCSD(T) and G2(+) values, show that the H(3) problem is not an isolated case, and while it is premature to conclude that VB theory has come of age, the occurrence of this event is clearly within sight.

Journal Article↗