Performance of dynamically weighted multiconfiguration self-consistent field and spin-orbit coupling calculations of diatomic molecules of Group 14 elements

The efficacy of several multiconfiguration self-consistent field (MCSCF) methods in the subsequent spin-orbit coupling calculations was studied. Three MCSCF schemes to generate molecular orbitals were analyzed: state-specific, state-averaged, and dynamically weighted MCSCF. With Sn(2)(+) as the representative case, we show that the state-specific MCSCF orbitals lead to discontinuities in potential energy curves when avoided crossings of electronic states occur; this problem can be solved using the state-averaged or dynamically weighted MCSCF orbitals. The latter two schemes are found to give similar results when dynamic electron correlation is considered, which we calculated at the level of multiconfigurational quasidegenerate perturbation theory (MCQDPT). We employed the recently developed Douglas-Kroll spin-orbit adapted model core potential, ZFK3-DK3, and the dynamically weighted MCSCF scheme to calculate the spectroscopic constants of the mono-hydrides and compared them to the results obtained using the older set of potentials, MCP-TZP. We also showed that the MCQDPT tends to underestimate the dissociation energies of the hydrides and discussed to what extent coupled-cluster theory can be used to improve results.     

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     

Potential energy curve for isomerization of N2H2 and C2H4 using the improved virtual orbital multireference Møller-Plesset perturbation theory

Multireference M?ller-Plesset (MRMP) perturbation theory [K. Hirao, Chem. Phys. Lett. 190, 374 (1992)] is modified to use improved virtual orbitals (IVOs) and is applied to study ground state potential energy curves for isomerization and dissociation of the N2H2 and C2H4 molecules. In contrast to traditional MRMP or multistate multiconfiguration quasidegenerate perturbation theory where the reference functions are obtained from (often difficult to converge) state averaged multiconfiguration self-consistent field methods, our reference functions are represented in terms of computationally efficient IVOs. For convenience in comparisons with other methods, a first order complete active space configuration interaction (CASCI) calculation with the IVOs is followed by the use of the IVOs in MRMP to incorporate residual electron correlation effects. The potential energy curves calculated from the IVO-MRMP method are compared with computations using state-of-the-art coupled cluster singles and doubles (CCSD) methods and variants thereof to assess the efficacy of the IVO-MRMP scheme. The present study clearly demonstrates that unlike the CCSD and its variants, the IVO-MRMP approach provides smooth and reliable ground state potential energy curves for isomerization of these systems. Although the rigorously size-extensive completely renormalized CC theory with noniterative triples corrections (CR-CC(2,3)) likewise provides relatively smooth curves, the CR-CC(2,3) calculations overestimate the cis-trans barrier height for N2H2. The ground state spectroscopic constants predicted by the IVO-CASCI method agree well with experiment and with other highly correlated ab initio methods.     

Low-lying excited states of C120 and C151: a multireference perturbation theory study

Excited states of two 7-aminocoumarin derivatives, coumarin 120 (7-amino-4-methylcoumarin) and coumarin 151 (7-amino-4-trifluoromethylcoumarin), were investigated using generalized multiconfigurational quasidegenerate perturbation theory (GMC-QDPT), multiconfigurational quasidegenerate perturbation theory (MC-QDPT) and time-dependent density functional theory (TDDFT) with the B3LYP and CAM-B3LYP functionals. The absorption and fluorescence spectra of C120 and C151 were calculated. We elucidated the characters of the low-lying states of C120 and C151. The absorption spectra calculated with GMC-QDPT and TDDFT B3LYP agreed well with the experimental data, while for the fluorescence spectra, the TDDFT calculations overestimated the fluorescence spectra compared to GMC-QDPT calculations. Utilizing active spaces with large numbers of electrons and orbitals for reference functions, GMC-QDPT showed a better performance than MC-QDPT with a complete active space self-consistent field (CASSCF) reference of active space with smaller number of electrons and orbitals. In our gas phase calculation, we found that the optimized structures for the first excited states have a planar amino group with a CN single bond, while the amino group is pyramidal in the ground state.     

The problem of hole localization in inner-shell states of N2 and CO2 revisited with complete active space self-consistent field approach

Potential energy curves for inner-shell states of nitrogen and carbon dioxide molecules are calculated by inner-shell complete active space self-consistent field (CASSCF) method, which is a protocol, recently proposed, to obtain specifically converged inner-shell states at multiconfigurational level. This is possible since the collapse of the wave function to a low-lying state is avoided by a sequence of constrained optimization in the orbital mixing step. The problem of localization of K-shell states is revisited by calculating their energies at CASSCF level based on both localized and delocalized orbitals. The localized basis presents the best results at this level of calculation. Transition energies are also calculated by perturbation theory, by taking the above mentioned MCSCF function as zeroth order wave function. Values for transition energy are in fairly good agreement with experimental ones. Bond dissociation energies for N(2) are considerably high, which means that these states are strongly bound. Potential curves along ground state normal modes of CO(2) indicate the occurrence of Renner-Teller effect in inner-shell states.     

The electronic structure and energetics of V(+)-benzene half-sandwiches of different multiplicities: Comparative multireference and single-reference theoretical study

Multireference [complete active space self-consistent field (CASSCF) and multiconfigurational quasidegenerate perturbation theory (MCQDPT)] and single-reference ab initio (Moller-Plesset second order perturbation theory (MP2) and coupled clusters with singles, doubles and noniterative triples [CCSD(T)]) and density functional theory (PBE and B3LYP) electronic structure calculations of V(C(6)H(6))(+) half-sandwich in the states of different multiplicities are described and compared. Detailed analyses of the geometries and electronic structures of the all found states are given; adiabatic and diabatic dissociation energies are estimated. The lowest electronic state of V(C(6)H(6))(+) half-sandwich was found to be the quintet (5)B(2) state with a slightly deformed upside-down-boat-shaped benzene ring and d(4) configuration of V atom, followed by a triplet (3)A(2) state lying about 4 kcal/mol above. The lowest singlet state (1)A(1)(d(4)) lies much ( approximately 28 kcal/mol) higher. MCQDPT calculated adiabatic dissociation energy (53.6 kcal/mol) for the lowest (5)B(2)(d(4)) state agrees well with the current 56.4 (54.4) kcal/mol experimental estimate, giving a preference to the lower one. Compared to MCQDPT, B3LYP hybrid exchange-correlation functional provides the best results, while CCSD(T) performs usually worse. Gradient-corrected PBE calculations tend to systematically overestimate metal-benzene binding in the row quintet相似文献   

A pair-excitation MCSCF calculation based on CNDO/INDO approximation I

Second-order multiconfigurational self-consistent field (MCSCF ) calculation has been programmed on the basis of CNDO /INDO molecular orbital bases, in which the configuration space employed is restricted within pair-excitations. Test calculations have been carried out for 17 small molecules. All the MCSCF ground states of these molecules have been successfully converged to their respective optimal states by employing a simple weighting scheme. This procedure provides a great savings in computer time. The MCSCF calculation on azetidine required only 27 min on a HITAC M-680H. The MCSCF energies of HF, F₂, and BH show improved behaviors up to large atomic distances (～7au).     

OClO里德堡态激发能的准确预测及其阴离子低能激发态的从头算研究

采用全活化空间自洽场方法(CASSCF)研究了OClO阴离子7个低能电子态及其自由基的基态. 为了进一步考虑动态电子相关效应, 采用二级多组态微扰理论(CASPT2)获得更加可靠的能量值. 此外, 在ANO-L基组的基础上, 在OClO自由基的电荷中心增加了为研究里德堡态所建立的1s1p1d的波函数, 并应用多组态二级微扰理论(MS-CASPT2)方法获得了里德堡态的准确电子激发能.     

Comparison of perturbative and multiconfigurational electron propagator methods

Ionization energies below 20 eV of 10 molecules calculated with electron propagator techniques employing Hartree-Fock orbitals and multiconfigurational self-consistent field orbitals are compared. Diagonal and nondiagonal self-energy approximations are used in the perturbative formalism. Three diagonal methods based on second- and third-order self-energy terms, all known as the outer valence Green's function, are discussed. A procedure for selecting the most reliable of these three versions for a given calculation is tested. Results with a polarized, triple ζ basis produce root mean square errors with respect to experiment of approximately 0.3 eV. Use of the selection procedure has a slight influence on the quality of the results. A related, nondiagonal method, known as ADC (3), performs infinite-order summations on several types of self-energy contributions, is complete through third-order, and produces similar accuracy. These results are compared to ionization energies calculated with the multiconfigurational spin-tensor electron propagator method. Complete active space wave functions or close approximations constitute the reference states. Simple field operators and transfer operators pertaining to the active space define the operator manifold. With the same basis sets, these methods produce ionization energies with accuracy that is comparable to that of the perturbative techniques. © 1996 John Wiley & Sons, Inc.     

Influence of multi-atom bridging ligands on the electronic structure and magnetic properties of homodinuclear titanium molecules

The electronic structure and magnetic properties of homodinuclear titanium(III) molecules with bridging ligands from groups 14, 15, and 16 are examined. Single- and multireference methods with triple-zeta plus polarization basis sets are employed. Dynamic electron correlation effects are included via second-order multireference perturbation theory. Isotropic interaction parameters are calculated, and two of the complexes studied are predicted to be ferromagnetic based on multireference second-order perturbation (MRMP2) theory, using the TZVP(fg) basis set. Zero-field splitting parameters are determined using spin-orbit coupling obtained from complete active space (CAS) self-consistent field (SCF) and multiconfigurational quasi-degenerate perturbation theory (MCQDPT) wave functions. Three Breit-Pauli-based spin coupling methods were employed: full Breit-Pauli (HSO2), the partial two-electron method (P2E), and the semiempirical one-electron method (HSO1).     

Structure and magnetic properties of iron-silicon clusters in a multiconfigurational calculation

Journal of Structural Chemistry - The structure and spin moments of iron-silicon clusters FeSi n are calculated by a multiconfigurational self-consistent field (MCSCF) method in order to refine the...     

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.     

Computation of conical intersections by using perturbation techniques

Multiconfigurational second-order perturbation theory, both in its single-state multiconfigurational second-order perturbation theory (CASPT2) and multistate (MS-CASPT2) formulations, is used to search for minima on the crossing seams between different potential energy hypersurfaces of electronic states in several molecular systems. The performance of the procedures is tested and discussed, focusing on the problem of the nonorthogonality of the single-state perturbative solutions. In different cases the obtained structures and energy differences are compared with available complete active space self-consistent field and multireference configuration interaction solutions. Calculations on different state crossings in LiF, formaldehyde, the ethene dimer, and the penta-2,4-dieniminium cation illustrate the discussions. Practical procedures to validate the CASPT2 solutions in polyatomic systems are explored, while it is shown that the application of the MS-CASPT2 procedure is not straightforward and requires a careful analysis of the stability of the results with the quality of the reference wave functions, that is, the size of the active space.     

Electronic circular dichroism of disulphide bridge: ab initio quantum-chemical calculations

Electronic circular dichroism (ECD) parameters of the disulphide chromophore have been calculated for dihydrogen disulphide, dimethyl disulphide, and cystine using density-functional theory, coupled-cluster theory, and multiconfigurational self-consistent field theory. The objective is twofold: first, to examine the performance of the Coulomb-attenuated CAM-B3LYP functional for the calculation of ECD spectra; second, to investigate the dependence of the ECD parameters on the conformation around the disulphide bridge. The CAM-B3LYP functional improves considerably on the B3LYP functional, giving results comparable to CCSD theory and to MCSCF theory in an extended active space. The conformational dependence of the ECD parameters does not change much upon substitution, which is promising for the application of ECD in structural investigations of proteins containing disulphide bridges.     

Equivalent orbitals for multiconfigurational spin‐tensor electron propagator method (MCSTEP): The vertical ionization potentials of B,NO, CF,and OF

The multiconfigurational spin tensor electron propagator method (MCSTEP) was developed as an implementation of electron propagator/single particle Green's function methods for ionization potentials (IPs) and electron affinities (EAs). MCSTEP was specifically designed for open shell and highly correlated (nondynamically correlated) initial states. For computational efficiency the initial state used in MCSTEP is typically a small complete active space (CAS) multiconfigurational self‐consistent field (MCSCF) state. If in a molecule there are some degenerate orbitals which are not fully or half occupied, usual MCSCF calculations will make these orbitals inequivalent, i.e., the occupied ones will be different from the nonoccupied ones, so that the degeneracy is broken. In this article, we use a state averaged MCSCF method to get equivalent orbitals for the initial state and import the integrals into the subsequent MCSTEP calculations. This gives, in general, more reliable MCSTEP vertical IPs. © 2008 Wiley Periodicals, Inc., 2008     

Intramolecular charge-transfer state formation of 4-(N,N-dimethylamino)benzonitrile in acetonitrile solution: RISM-SCF study

Intramolecular charge-transfer (ICT) state formation of 4-(N,N-dimethylamino)benzonitrile in acetonitrile solution is studied by the reference interaction site model self-consistent field (RISM-SCF) method. Geometry optimizations are performed for each electronic state in solution with the complete-active-space SCF wave functions. Dynamic electron correlation effects are taken into account by using the multiconfigurational quasidegenerate perturbation theory. Two-dimensional free energy surfaces are constructed as the function of the twisting and wagging angles of the dimethylamino group for the ground and locally excited (LE) states. The calculated absorption and fluorescence energies are in good agreement with experiments. The validity of the twisted ICT (TICT) model is confirmed in explaining the dual fluorescence, and the possibility of the planar ICT model is ruled out. To examine the mechanism of the TICT state formation, a "crossing" seam between the LE and charge-transfer (CT) state surfaces is determined. The inversion of two electronic states occurs at a relatively small twisting angle. The effect of solvent reorganization is also examined. It is concluded that the intramolecular twisting coordinate is more important than the solvent fluctuation for the TICT state formation, because the energy difference between the two states is minimally dependent on the solvent configuration.     

Ab initio study of the ground and lower-lying excited electronic states of NiX2 and FeX2 (X=F,Cl, Br,I) molecules

The ground and lower-lying excited electronic states of FeX₂ and NiX₂ (X=F, Cl, Br, I) molecules are systematically investigated by ab initio method at the complete active space self-consistent field (CASSCF) and multiconfigurational quasi-degenerate second-order perturbation (MCQDPT2) levels of theory. It is concluded that the dynamic electron correlation has to be taken into account in the prediction of the properties for such kind of molecules. The equilibrium bond lengths r_e(M–X), force constants and harmonic vibrational frequencies are calculated for the ground and lower-lying excited electronic states. The spin-orbit coupling (SOC) effects are analysed.     

Modification of optimal complete active space choices for the multiconfigurational spin‐tensor electron propagator method for ionization potentials

The multiconfigurational spin tensor electron propagator (MCSTEP) method 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. Ionization or electron attachment is always from a state of pure spin symmetry to a state of pure spin symmetry even if the initial state is open shell. MCSTEP can be used as well for molecules with initial states that can be accurately described by a single determinant‐based theory. The initial state that is used in MCSTEP is typically a small complete active space (CAS) multiconfigurational self‐consistent field (MCSCF) state. We previously examined different small CAS choices for MCSTEP initial states and have developed a generally workable scheme. This article further examines some different ways to choose the CAS for MCSTEP. With several logical CAS choices, we have calculated the low‐lying vertical MCSTEP ionization potentials (IPs) of C₂, N₂, linear H₂O, O₂, CH₂, and NH₂, comparing them with large multireference configuration interaction (MRCI) calculations. We conclude that generally a small modification and extension of our previous schemes for choosing the MCSTEP CAS gives IPs that most effectively mimics the results of large scale MRCI IPs in general. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007