Electron correlation and relativistic effects in atomic structure calculations of the thorium atom

Relativistic two-component ab initio calculations have been performed for the Th atom. The spin free low lying states have been calculated at state-averaged complete active space self-consistent field (SA-CASSCF) and multi-state complete active space second-order perturbation (MS-CASPT2) level of theories using different sets of active orbitals. The spin-orbit states have been computed using Douglas-Kroll type of atomic mean-field integral approach. The effects of dynamic electron correlation have been studied at the MS-CASPT2 level. The energy levels of spin-orbit states below 30,000 cm(-1) obtained by the inclusion of dynamic electron correlation are in very good agreement with the experimental values. The radiative properties such as weighted transition probabilities (gA) and oscillator strengths (gf) among several spin-orbit states have been calculated at the SA-CASSCF and MS-CASPT2 levels and are expected to be very helpful for future experiments.     

Conformational analysis of singlet-triplet state mixing in Paternò-Büchi diradicals

Conformational dependence of spin-orbit coupling (SOC) in flexible Paternò-Büchi (PB) diradicals has been studied with high-level ab initio methods using both (i) one-electron spin-orbit Hamiltonian with parametrized (effective) nuclear charges in conjunction with a state-averaged MCSCF wave function as implemented by Robb in Gaussian 98 and (ii) complete one- and two-electron SOC with a fully optimized MCSCF triplet wave function and frozen core singlet as implemented by Furlani in the GAMESS computational package. The ab initio results revealed two distinct areas of elevated SOC values, one corresponding to the region whereby a cisoid conformation in the C-C-O-C fragment brings the two odd-electron orbitals closer to each other, and the other area corresponding to the partially eclipsed conformation lacking direct overlap between the spin centers. In this second region the 1,4-electronic communication is mediated by the oxygen's 2p-lone pair, which is suitably oriented to play the role of a "relay-antenna". The other critical factor affecting the rate of intersystem crossing (ISC)--singlet-triplet energy separation--was computed utilizing a multireference CASSCF-MP2 method to include dynamic correlation effects. The largest singlet-triplet energy gap, approximately 2 kcal/mol, was found for a gauche conformer (also a minimum SOC conformation). Rotation about the central C-O bond either toward the fully eclipsed (0 degrees ) or the partially eclipsed (120 degrees ) conformations decreases the singlet-triplet gap while increasing the value of the SOC matrix element. These computational findings support the Griesbeck model for stereochemistry of triplet PB reactions and provide a rigorous basis for predicting the probability of ISC in diradicals separated by a partially conjugated spacer.     

Accurate calculation of core-electron binding energies: multireference perturbation treatment

Multireference perturbation theory (MRPT) with multiconfigurational self-consistent field (MCSCF) reference functions is applied to the calculations of core-electron binding energies (CEBEs) of atoms and molecules. Orbital relaxations in a core-ionized state and electron correlation are both taken into account in a conventional MCSCF-MRPT procedure. In the MCSCF calculation, the target core ionized state is directly optimized as an excited state and this treatment can completely prevent a variational collapse. Multireference Moller-Plesset perturbation theory and multiconfigurational self-consistent field reference quasidegenerated perturbation theory were used to treat electron correlation. The present method quite accurately reproduced the 1s CEBEs of CH4, NH3, H2O, and FH; the average deviation from the experimental data is 0.11 eV using Ahlrichs' VTZ basis set. The C 1s and O 1s CEBEs of formic acid and acetic acid were calculated and the results are consistent with the bonding characters of the atoms in these molecules. The present procedure can also be applied to CEBEs of higher angular momentum orbitals by including spin-orbit coupling. The calculated CEBEs of Ar 2p, HCl 2p, Kr 3d, and HBr 3d are in reasonable agreement with the available experimental values. In the calculation of the 3d CEBEs, a relativistic correction significantly improves the agreements. The effect of polarization functions is also discussed.     

Generalized-molecular-orbital theory: Simple multiconfiguration self-consistent-field method

Even after completing a multiconfiguration self-consistent-field (MCSCF ) calculation, one must often include additional configuration interaction (CI ) to obtain quantitative or semiquantitative results. There is some question of whether the prior MCSCF calculation is worthwhile, if additional CI is needed later. We have developed a new MCSCF computational method, which, because of our assumptions about the nature of the configurations, yields one Fock-like operator for all the “filled” orbitals (high occupation numbers) and a second Fock-like operator for all the “virtual” orbitals (low occupation numbers). Since there are only two matrices to build, our method is considerably faster than other MCSCF approaches. Because of these similarities to standard molecular-orbital (MO ) calculations, we have termed our approach generalized-molecular-orbital (GMO ) theory. However, the “virtual” orbitals, unlike those of standard MO theory, are optimized to correlate the “filled” ones and can he used in a subsequent CI calculation. Results are presented for the correlation energy of H₂O, the spectroscopic constants of N₂, the singlet–triplet energy separations in CH₂, and the nature of the chromium–chromium quadruple bond. Although these results are at a very low level of CI , the GMO approach appears to correct for the gross deficiencies of the single-determinant SCF procedure.     

Theoretical investigation of the vibronic spectrum in the X(2)Pi(u) electronic state of C(6)(+)

In this study we employ the recently developed model for handling the Renner-Teller effect in Pi electronic states of six-atomic molecules with linear equilibrium geometry to calculate the vibronic spectrum in the X(2)Pi(u) electronic state of the C(6)(+) ion. The applied model Hamiltonian excludes the stretching vibrations and end-over-end rotations. On the other hand, it considers the interplay between the vibronic and spin-orbit couplings. The parameters determining the shape of the bending potential energy surfaces are computed by means of a Density functional theory, and the spin-orbit coupling constant by the Multireference CI program using state-averaged complete active space self-consistent field (SA-CASSCF) wavefunctions. The results of the present study are expected to motivate and help future experimental investigations on C(6)(+).     

Theoretical study of some bis-verdazyl diradicals: singlet–triplet energy gap

The bis-verdazyl diradical (BVD) system is closely examined by using the multiconfiguration wavefunctions as well as the density functional theory (DFT). The totally symmetric singlet ground state turns out to have strong multiconfiguration character at all levels of theory. The singlet ground state takes on the planar structure while the most stable triplet state corresponds to the twisted form. The MCSCF+MCQDPT2 calculations are shown to be sufficient to predict the singlet–triplet energy gap which is insensitive to the electronic characters of the ring substituents.     

Gauge-origin-independent magnetizabilities of solvated molecules using the polarizable continuum model

We present an implementation of the polarizable continuum model in its integral equation formulation for the calculation of the magnetizabilities of solvated molecules. The gauge-origin independence of the calculated magnetizabilities and the fast basis set convergence are ensured through the use of London atomic orbitals. Our implementation can use Hartree-Fock and multiconfigurational self-consistent-field (MCSCF) wave functions as well as density-functional theory including hybrid functionals such as B3LYP. We present the results of dielectric continuum effects on water and pyridine using MCSCF wave functions, as well as dielectric medium effects on the magnetizability of the aromatic amino acids as a model for how a surrounding protein environment affects the magnetizability of these molecules. It is demonstrated that the dielectric medium effects on the magnetizability anisotropies of the aromatic amino acids may be substantial, being as large as 25% in the case of tyrosine.     

Superior performance of Mukherjee’s state-specific multi-reference coupled-cluster theory at the singles and doubles truncation scheme with localized active orbitals

The state-specific multi-reference coupled-cluster (SS-MRCC) theory of Mukherjee et al., in its singles and doubles truncation scheme (SS-MRCCSD), misses important couplings between the virtual functions reached by single and double excitations from different model functions. Since the SS-MRCC theory is not invariant with respect to the transformations among the active orbitals, the results are dependent on the active orbitals chosen. We demonstrate in this paper with results for potential energy curves for several example molecules involving single and multiple bond dissociation that the performance of SS-MRCCSD is significantly improved if localized active orbitals are used. The improvement is remarkable both in terms of the non-parallelity error and the magnitude of correlation energy recovered vis-a-vis the full configuration interaction results with the same basis set. The results bolster our claim that SS-MRCCSD with localized orbitals is an accurate general theory for potential energy surfaces.     

Unitary group adapted state specific multireference perturbation theory: Formulation and pilot applications

We present here a comprehensive account of the formulation and pilot applications of the second‐order perturbative analogue of the recently proposed unitary group adapted state‐specific multireference coupled cluster theory (UGA‐SSMRCC), which we call as the UGA‐SSMRPT2. We also discuss the essential similarities and differences between the UGA‐SSMRPT2 and the allied SA‐SSMRPT2. Our theory, like its parent UGA‐SSMRCC formalism, is size‐extensive. However, because of the noninvariance of the theory with respect to the transformation among the active orbitals, it requires the use of localized orbitals to ensure size‐consistency. We have demonstrated the performance of the formalism with a set of pilot applications, exploring (a) the accuracy of the potential energy surface (PES) of a set of small prototypical difficult molecules in their various low‐lying states, using natural, pseudocanonical and localized orbitals and compared the respective nonparallelity errors (NPE) and the mean average deviations (MAD) vis‐a‐vis the full CI results with the same basis; (b) the efficacy of localized active orbitals to ensure and demonstrate manifest size‐consistency with respect to fragmentation. We found that natural orbitals lead to the best overall PES, as evidenced by the NPE and MAD values. The MRMP2 results for individual states and of the MCQDPT2 for multiple states displaying avoided curve crossings are uniformly poorer as compared with the UGA‐SSMRPT2 results. The striking aspect of the size‐consistency check is the complete insensitivity of the sum of fragment energies with given fragment spin‐multiplicities, which are obtained as the asymptotic limit of super‐molecules with different coupled spins. © 2015 Wiley Periodicals, Inc.     

Dynamically weighted multiconfiguration self-consistent field: multistate calculations for F+H2O-->HF+OH reaction paths

A novel method of dynamically adjusted weighting factors in state-averaged multiconfiguration self-consistent-field calculations (SA-MCSCF) is described that is applicable to systems of arbitrary dimensionality. The proposed dynamically weighted approach automatically weights the relevant electronic states in each region of the potential energy surface, smoothly adjusting between these regions with an energy dependent functional. This method is tested on the F(2P)+H2O-->HF+OH(2Pi) reaction, which otherwise proves challenging to describe with traditional SA-MCSCF methods due to (i) different asymptotic degeneracies of reactant (threefold) and product (twofold) channels, and (ii) presence of low-lying charge transfer configurations near the transition state region. The smoothly varying wave functions obtained by dynamically weighted multiconfigurational self-consistent field represent excellent reference states for high-level multireference configuration interaction calculations and offer an ideal starting point for construction of multiple state potential energy surfaces.     

Towards the Development and Applications of Manifestly Spin-free Multi-reference Coupled Electron-pair Approximation-like Methods: A State Specific Approach

As a practical tool of being applicable to bigger molecules, a full-blown state-specific multi-reference coupled cluster formalism developed by us (Mahapatra et al. in J Chem Phys 110:6171, 1999) would be rather demanding computationally, and it is worthwhile to look for physically motivated approximation schemes which capture a substantial portion of the correlation of the full-blown theory. In this spirit, we have recently proposed coupled electron-pair approximation (CEPA)-like various approximants to the parent spin-adapted state-specific multi-reference coupled cluster (SS-MRCC) theory which depend on the inclusion of EPV terms to various degree. Here, the space of excitations is confined to the first order interactive virtual space generated by the cluster operator, but the EPV terms are included exactly. We call them spin-free state specific multi-reference CERA (SS-MRCEPA) theories. They work within the complete active space (CAS) and have been found to be very effective in bypassing the intruders, similar in performance to that of the parent SS-MRCC theory. The spin-adaptation of the working equations of both the SS-MRCC and the CEPA-like approximants is a non-trivial exercise. In this paper, we delineate briefly the essentials of a spin-free formulation of the SS-MRCC and SS-MRCEPA theories. This allows us to include open-shell configuration state functions (CSF) in the CAS. We consider three variants of SS-MRCEPA method. Two are explicitly orbital invariant: (1) SS-MRCEPA(0), a purely lineralized version of the SS-MRCC theory, (2) SS-MRCEPA(I), which includes all the EPV terms explicitly and exactly in an orbital invariant manner and (3) the SS-MRCEPA(D), which emerges when we keep only the diagonal terms of a set of dressed operators in the working equations. Unlike the first two, the third version is not invariant under the orbital transformation within the set of doubly occupied core, valence and virtual orbitals. The SS-MRCEPA methods produce very encouraging results as was evidenced in the applications on the computation of potential energy surfaces for the ground states of LiH and HF molecules.     

Pathway analysis of super-exchange electronic couplings in electron transfer reactions using a multi-configuration self-consistent field method

We present a novel pathway analysis of super-exchange electronic couplings in electron transfer reactions using localized molecular orbitals from multi-configuration self-consistent field (MCSCF) calculations. In our analysis, the electronic coupling and the tunneling pathways can be calculated in terms of the configuration interaction (CI) Hamiltonian matrix obtained from the localized MCSCF wave function. Making use of the occupation restricted multiple active spaces (ORMAS) method can effectively produce the donor, acceptor, and intermediate configuration state functions (CSFs) and CIs among these CSFs. In order to express the electronic coupling as a sum of individual tunneling pathways contributions, we employed two perturbative methods: L?wdin projection-iteration method and higher-order super-exchange method. We applied them to anion couplings of butane-1,4-diyl and pentane-1,5-diyl. The results were (1) the electronic couplings calculated from the two perturbative methods were in reasonable agreement with those from a non-perturbative method (one-half value of the energy difference between the ground and first excited states), (2) the main tunneling pathways consisted of a small number of lower-order super-exchange pathways where bonding, anti-bonding, or extra-valence-shell orbitals were used once or twice, and (3) the interference among a huge number of higher-order super-exchange pathways significantly contributed to the overall electronic coupling, whereas each of them contributed only fractionally. Our method can adequately take into account both effects of non-dynamical electron correlation and orbital relaxation. Comparing with the analyses based on the Koopmans' theorem (ignoring both effects) and the ORMAS-CIs from frozen localized reference orbitals (ignoring the effect of orbital relaxation), we discuss these effects.     

Analyzing Pt chemical shifts calculated from relativistic density functional theory using localized orbitals: the role of Pt 5d lone pairs

Pt chemical shifts were calculated from two-component relativistic density functional theory (DFT). The shielding tensors were analyzed by using a recently developed method to decompose the spin-orbit DFT results into contributions from spin-free localized orbitals (here: natural localized molecular orbitals (NLMOs) and natural bond orbitals (NBOs)). Seven chemical shifts in six Pt complexes with Pt oxidation states II, III, and IV; and halide, amino, and amidate ligands were analyzed, with particular focus on the role of nonbonding Pt 5d orbitals. A simple d-orbital 'rotation' model has been used to rationalize some of the observed trends such as the main difference between Pt(II) and Pt(IV) chemical shifts. The localized orbital analysis data showed that most of this difference as well as trends among different Pt complexes with similar coordination can be rationalized by comparing properties of the nonbonding Pt 5d orbitals. We have also analyzed the spin-orbit effects on the chemical shifts of [PtCl4](2-) compared to [PtBr4](2-).     

Theoretical study of low-lying electronic states of Mn2 using a newly developed relativistic model core potential

We investigated the electronic structure of low-lying electronic states of Mn₂ using a newly developed relativistic model core potential (spdsMCP). Calculations were performed at complete active space self-consistent field (CASSCF) and second-order multiconfiguration quasidegenerate perturbation theory (MCQDPT2) levels. The MCQDPT2 calculations reveal that the 1Σg⁺ state is the ground state. Calculated spectroscopic constants are very similar to the results of recent all-electron calculations and experimental values, indicating that the spdsMCP works well for Mn₂, which requires a highly correlated calculation. The wave functions of low-lying states are also analyzed at the CASSCF level.     

Ab initio study of bond stretching: Implications in force-field parametrization for molecular mechanics and dynamics

We report a theoretical study of the stretching of chemical bonds and its implications on the force-field parametrization. Computations were performed at the SCF and MCSCF levels by using minimal, split-valence, and large extended and polarized basis sets. The stretching energy profiles were determined considering up to 25 perturbed geometries of 11 different bonds (6 singles, 2 doubles, and 3 triples). The energy profiles and stretching parameters are compared with the experimental data compiled in the most popular force fields. MCSCF stretching energy profiles are mainly anharmonic and can be only roughly reproduced by quadratic equations. The use of Allinger's MM2 quasiharmonic expression appears as the best choice because it fits with reasonable accuracy a large percentage of the stretching profile without increasing the complexity of the formalism and of the parametrization procedure. MCSCF computations are needed to obtain reliable stretching force parameters. In this respect, MCSCF calculations considering as active space only the bonded and nonbonded orbitals of the perturbed bond seems to be the best strategy to obtain good results at a minimum computational cost, especially if small split-valence basis sets like the 3-21G are used. Results obtained at this level of sophistication are completely comparable to stretching parameters compiled on empirical force fields. © 1993 John Wiley & Sons, Inc.     

大体系多电子相关研究中应用群对称定域轨道的构想

大体系多电子相关研究中应用群对称定域轨道的构想周泰锦，刘爱民（厦门大学化学系，厦门３６１００５）关键词：组态相关，多构型自治叠代，多中心积分，群对称定域轨道，对称约化有关原子簇化合物及化学吸附、过渡态、激发态、催化反应等大体系的量子化学研究，对于探讨...     

A simplified relativistic time-dependent density-functional theory formalism for the calculations of excitation energies including spin-orbit coupling effect

In the present work we have proposed an approximate time-dependent density-functional theory (TDDFT) formalism to deal with the influence of spin-orbit coupling effect on the excitation energies for closed-shell systems. In this formalism scalar relativistic TDDFT calculations are first performed to determine the lowest single-group excited states and the spin-orbit coupling operator is applied to these single-group excited states to obtain the excitation energies with spin-orbit coupling effects included. The computational effort of the present method is much smaller than that of the two-component TDDFT formalism and this method can be applied to medium-size systems containing heavy elements. The compositions of the double-group excited states in terms of single-group singlet and triplet excited states are obtained automatically from the calculations. The calculated excitation energies based on the present formalism show that this formalism affords reasonable excitation energies for transitions not involving 5p and 6p orbitals. For transitions involving 5p orbitals, one can still obtain acceptable results for excitations with a small truncation error, while the formalism will fail for transitions involving 6p orbitals, especially 6p1/2 spinors.     

Approximate MCSCF optimization for multiconfigurational spin‐tensor electron propagator method (MCSTEP): The vertical ionization potentials of CO,HCN, HNC,H2CO,and O3

The multiconfigurational spin tensor electron propagator method (MCSTEP) was developed as an implementation of electron propagator/single particle Green's function methods. MCSTEP was specifically designed for open shell and highly correlated (nondynamically correlated) initial states. The initial state used in MCSTEP is typically a small complete active space (CAS) with multiconfigurational self‐consistent field (MCSCF) state. In some cases, because of our use of a small CAS in MCSTEP, the Lagrangian eigenvalues of the MCSCF reference state are in an undesired order (u). The desired order (d) can usually be obtained by excluding one or more orbital rotations in MCSCF optimization between the doubly occupied and partially occupied orbitals. We systematically examine several cases where the undesired order occurs for the low‐lying vertical MCSTEP ionization potentials (IPs) of the molecules CO, HCN, HNC, H₂CO, and O₃ with our recently established CAS choices for MCSCF/MCSTEP. By excluding one or more orbital rotations between the partially and doubly occupied orbitals, an approximate MCSCF reference state with the same CAS choice is obtained for use in standard MCSTEP calculations that, in general, gives more reliable vertical MCSTEP IPs. © 2007 Wiley Periodicals, Inc. J Quantum Chem, 2008     

On the choice of active space orbitals in MCSCF calculations

We describe a procedure which may be used to aid selection of the active space in multiconfigurational self-consistent field (MCSCF) calculations for general chemical systems. Starting from a restricted Hartree-Fock calculation, we define a hierarchy of interacting virtual orbitals for every occupied orbital. The most strongly interacting orbitals are then taken to constitute the active space in a configuration interaction (CI) calculation. The natural orbital occupation numbers obtained from the CI calculation are then used to choose the active space to be used in a subsequent MCSCF calculation. We illustrate our method on a number of systems (Li₂, B₂, C₂, carbonyl oxide and the transition state for oxidation of H₂S by dioxirane). In all these cases, ‘intuitive’ active spaces are inadequate, as are active spaces derived from the natural orbitals of unrestricted Hartree-Fock calculations.