A systematic improvement to UGA-SSMRCCSD equations and its implication for potential energy curves.

The Unitary Group Adaptation (UGA) offers a very compact and efficient spin adaptation strategy for any spin-free Hamiltonian in a many body framework. Our use of UGA in the context of state-specific (SS) Jeziorski-Monkhorst Ansatz based multireference coupled cluster (MRCC) theory obviates the non-commutativity between the spin-free cluster operators via a normal ordered exponential parametrization in the wave operator. A previous formulation of UGA-SSMRCC by us [R. Maitra, D. Sinha, and D. Mukherjee, J. Chem. Phys. 137, 024105 (2012)], using the same ansatz, employed certain sufficiency conditions to reach the final working equations, which cannot be improved systematically. In this article, we will present a more rigorous formulation that follows from an exact factorization of the unlinked terms of the Bloch equation, resulting in equations on which a hierarchy of approximations can be systematically performed on the emergent additional terms. This derivation was shown in our recent article [D. Chakravarti, S. Sen, and D. Mukherjee, Mol. Phys. 119, e1979676 (2021)] in the context of a single open shell CC formalism and was applied to spectroscopic energy differences where the contribution of the new terms was found to be of the order of ∼0.001 eV for ionization potential, electron affinity, and excitation energy. In the current work, we will present a comparison between the earlier and current formulations via both a theoretical analysis and a numerical demonstration of the dramatic effect of the additional terms brought in by the factorization on potential energy curves. The contribution of such terms was found to gain importance with an increase in the number of singly occupied active orbitals in the model space functions.

Exploration of interlacing and avoided crossings in a manifold of potential energy curves by a unitary group adapted state specific multi-reference perturbation theory (UGA-SSMRPT).

The Unitary Group Adapted State-Specific Multi-Reference Perturbation Theory (UGA-SSMRPT2) developed by Mukherjee et al. [J. Comput. Chem. 36, 670 (2015)] has successfully realized the goal of studying bond dissociation in a numerically stable, spin-preserving, and size-consistent manner. We explore and analyze here the efficacy of the UGA-SSMRPT2 theory in the description of the avoided crossings and interlacings between a manifold of potential energy curves for states belonging to the same space-spin symmetry. Three different aspects of UGA-SSMRPT2 have been studied: (a) We introduce and develop the most rigorous version of UGA-SSMRPT2 that emerges from the rigorous version of UGA-SSMRCC utilizing a linearly independent virtual manifold; we call this the "projection" version of UGA-SSMRPT2 (UGA-SSMRPT2 scheme P). We compare and contrast this approach with our earlier formulation that used extra sufficiency conditions via amplitude equations (UGA-SSMRPT2 scheme A). (b) We present the results for a variety of electronic states of a set of molecules, which display the striking accuracy of both the two versions of UGA-SSMRPT2 with respect to three different situations involving weakly avoided crossings, moderate/strongly avoided crossings, and interlacing in a manifold of potential energy curves (PECs) of the same symmetry. Accuracy of our results has been benchmarked against IC-MRCISD + Q.

Unitary coupled-cluster based self-consistent polarization propagator theory: A third-order formulation and pilot applications.

In this article, the development of a third-order self-consistent polarization propagator method based on unitary coupled-cluster (UCC) parametrization of the ground-state wavefunction and the excitation manifold comprising unitary-transformed excitation operators, hereafter referred to as UCC3, is reported. The UCC3 method is designed to provide excitation energies correct up to the third order for excited states dominated by single excitations. An expansion for the UCC transformed Hamiltonian involving Bernoulli numbers as expansion coefficients is adopted in the derivation of UCC3 working equations. Interestingly, UCC-based polarization propagator theory offers an alternative derivation for the strict version of the third-order algebraic diagrammatic construction [ADC(3)-s] method. The UCC3 results for the excitation energies of excited states in H2O, HF, N2, Ne, CH2, BH, and C2 molecules are compared with benchmark full configuration interaction values as well as ADC(3) and equation-of-motion coupled-cluster singles and doubles results to demonstrate the accuracy of the UCC3 method. UCC-based self-consistent polarization propagator theory appears to be a promising framework for developing non-perturbative hermitian formulations for treating electronically excited states.

Inclusion of orbital relaxation and correlation through the unitary group adapted open shell coupled cluster theory using non-relativistic and scalar relativistic Hamiltonians to study the core ionization potential of molecules containing light to medium-heavy elements.

The orbital relaxation attendant on ionization is particularly important for the core electron ionization potential (core IP) of molecules. The Unitary Group Adapted State Universal Coupled Cluster (UGA-SUMRCC) theory, recently formulated and implemented by Sen et al. [J. Chem. Phys. 137, 074104 (2012)], is very effective in capturing orbital relaxation accompanying ionization or excitation of both the core and the valence electrons [S. Sen et al., Mol. Phys. 111, 2625 (2013); A. Shee et al., J. Chem. Theory Comput. 9, 2573 (2013)] while preserving the spin-symmetry of the target states and using the neutral closed-shell spatial orbitals of the ground state. Our Ansatz invokes a normal-ordered exponential representation of spin-free cluster-operators. The orbital relaxation induced by a specific set of cluster operators in our Ansatz is good enough to eliminate the need for different sets of orbitals for the ground and the core-ionized states. We call the single configuration state function (CSF) limit of this theory the Unitary Group Adapted Open-Shell Coupled Cluster (UGA-OSCC) theory. The aim of this paper is to comprehensively explore the efficacy of our Ansatz to describe orbital relaxation, using both theoretical analysis and numerical performance. Whenever warranted, we also make appropriate comparisons with other coupled-cluster theories. A physically motivated truncation of the chains of spin-free T-operators is also made possible by the normal-ordering, and the operational resemblance to single reference coupled-cluster theory allows easy implementation. Our test case is the prediction of the 1s core IP of molecules containing a single light- to medium-heavy nucleus and thus, in addition to demonstrating the orbital relaxation, we have addressed the scalar relativistic effects on the accuracy of the IPs by using a hierarchy of spin-free Hamiltonians in conjunction with our theory. Additionally, the contribution of the spin-free component of the two-electron Gaunt term, not usually taken into consideration, has been estimated at the Self-Consistent Field (ΔSCF) level and is found to become increasingly important and eventually quite prominent for molecules with third period atoms and below. The accuracies of the IPs computed using UGA-OSCC are found to be of the same order as the Coupled Cluster Singles Doubles (ΔCCSD) values while being free from spin contamination. Since the UGA-OSCC uses a common set of orbitals for the ground state and the ion, it obviates the need of two N5 AO to MO transformation in contrast to the ΔCCSD method.

4-Component relativistic calculations of L3 ionization and excitations for the isoelectronic species UO2(2+), OUN(+) and UN2.

We present a 4-component relativistic study of uranium 2p3/2 ionization and excitation in the isoelectronic series UO2(2+), OUN(+) and UN2. We calculate ionization energies by ΔSCF at the Hartree-Fock (HF) and Kohn-Sham (KS) level of theory. At the ΔHF level we observe a perfectly linear chemical shift of ionization energies with respect to uranium atomic charges obtained from projection analysis. We have also developed a non-canonical 2nd-order Møller-Plesset code for wave function based correlation studies. We observe the well-known failure of Koopmans' theorem for core ionization due to the dominance of orbital relaxation over electron correlation effects. More unexpectedly, we find that the correlation contribution has the same sign as the relaxation contribution and show that this is due to a strong coupling of relaxation and correlation. We simulate uranium L3 XANES spectra, dominated by 2p3/2 → U6d transitions, by restricted excitation window time-dependent density functional theory (REW-TDDFT) and the complex polarization propagator (CPP) approach and demonstrate that they give identical spectra when the same Lorentz broadening is chosen. We also simulate XANES spectra by the Hartree-Fock based static exchange (STEX) method and show how STEX excitation energies can be reproduced by time-dependent Hartree-Fock calculations within the Tamm-Dancoff approximation. We furthermore show that Koopmans' theorem provide a correct approximation of ionization energies in the linear response regime and use this observation to align REW-TDDFT and CPP spectra with STEX ones. We point out that the STEX method affords the most detailed assignment of spectra since it employs virtual orbitals optimized for the selected core ionization. The calculated XANES spectra reflect the loss of bound virtual orbitals as the molecular charge is reduced along the isoelectronic series.

Aspects of Size-Consistency of Orbitally Noninvariant Size-Extensive Multireference Perturbation Theories: A Case Study Using UGA-SSMRPT2 as a Prototype.

Profiling a potential energy surface (PES), all the way to dissociate a molecular state into particular fragments and to display real or avoided crossings, requires a multireference description and the maintenance of size-consistency. The many body methods, which suit this purpose, should thus be size-extensive. Size-extensive theories, which are invariant with respect to transformation among active orbitals are, in principle, size-consistent. Relatively cheaper size-extensive theories, which do not possess this invariance, can still be size-consistent if the active orbitals are localized on the asymptotic fragments. Such methods, if perturbative in nature, require the use of an unperturbed Hamiltonian, which has orbital invariance with respect to the transformation within active, core, and virtual orbitals. The principal focus of this paper is to numerically realize size-consistency with localized active orbitals using our recently developed orbitally noninvariant Unitary Group Adapted State Specific Multireference second order Perturbation Theory (UGA-SSMRPT2) as a prototype method. Our findings expose certain generic potential pitfalls of size-extensive but orbitally noninvariant MRPT theories, which are mostly related to the inability of reaching proper localized active orbitals in the fragments due to the artifacts of the orbital generation procedure. Despite the invariance of the zeroth order CAS function, lack of invariance of the MRPT itself then leads to size-inconsistency. In particular, reaching symmetry broken fragment active orbitals is an issue of concern where suitable state-averaging might ameliorate the problem, but then one has to abandon full orbital optimization. Additionally, there can be situations where the orbitals of the fragment reached as an asymptote of the supermolecule are not the same as those obtained from the optimization of the fragments individually and will require additional transformation. Moreover, for a certain PES, one may either abandon the use of optimized orbitals for that state to preserve proper symmetry and degeneracy in the fragment orbitals or be satisfied with the use of optimized orbitals, which generate broken symmetric orbitals in the fragmentation limit. All these pathologies are illustrated using the PES of various electronic states of multiply bonded systems like N2, C2H2, HCN, C2, and O2. Subject to such proviso, the UGA-SSMRPT2 turns out to be an excellent theory for studying the PES leading to fragmentation of strongly correlated systems satisfying the requirements of size-consistency with localized active orbitals. An unexpected spin-off of our studies is the realization that the size-inextensive MRMP2, which bears a close structural similarity with our theory, might under certain situations display size-consistency. We analyze this feature concretely in our paper. Our studies may serve as a benchmark for monitoring numerically the size-consistency of any state specific multireference theory which is size-extensive but not invariant with respect to transformation of active orbitals.

Pattern of Drug Use and Associated Behaviors Among Female Injecting Drug Users From Northeast India: A Multi-Centric, Cross-Sectional, Comparative Study.

BACKGROUND: Studies from developed countries document the presence of injecting drug use among females and significantly higher vulnerabilities and risks as compared with male injecting drug users (IDUs). Studies comparing vulnerabilities and drug use patterns between female and male IDUs are not available for developing countries. OBJECTIVES: The aim of the study was to assess the drug use pattern and related HIV vulnerabilities among female IDUs and compare these findings with those from male IDUs from four states of Northeast India. METHOD: The study used data collected as part of a nationwide study of drug use pattern and related HIV vulnerabilities among IDUs. Ninety-eight female and 202 male IDUs accessing services from harm reduction sites across the four states of Northeast region of India were chosen through random sampling methodology. Drug use pattern, injecting practices, and knowledge of HIV were assessed using a structured questionnaire. RESULTS: Significantly higher proportion of female IDUs was uneducated, unemployed, reported their occupation as sex workers, and switched to injecting drug use faster as compared with male IDUs. Female IDUs practicing sex work differed significantly from those who did not with respect to frequency of daily injections, choice of drugs injected, and concomitant use of non-injecting drugs. More than half of female IDUs initiated sharing within the first month of injecting. CONCLUSIONS: The study demonstrates that female IDUs differ from male IDUs in their drug use pattern, initiation into injection as well as injecting behavior, which would be an important consideration during designing of female-specific interventions.

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.

Excited states with internally contracted multireference coupled-cluster linear response theory.

In this paper, the linear response (LR) theory for the variant of internally contracted multireference coupled cluster (ic-MRCC) theory described by Hanauer and Köhn [J. Chem. Phys. 134, 204211 (2011)] has been formulated and implemented for the computation of the excitation energies relative to a ground state of pronounced multireference character. We find that straightforward application of the linear-response formalism to the time-averaged ic-MRCC Lagrangian leads to unphysical second-order poles. However, the coupling matrix elements that cause this behavior are shown to be negligible whenever the internally contracted approximation as such is justified. Hence, for the numerical implementation of the method, we adopt a Tamm-Dancoff-type approximation and neglect these couplings. This approximation is also consistent with an equation-of-motion based derivation, which neglects these couplings right from the start. We have implemented the linear-response approach in the ic-MRCC singles-and-doubles framework and applied our method to calculate excitation energies for a number of molecules ranging from CH2 to p-benzyne and conjugated polyenes (up to octatetraene). The computed excitation energies are found to be very accurate, even for the notoriously difficult case of doubly excited states. The ic-MRCC-LR theory is also applicable to systems with open-shell ground-state wavefunctions and is by construction not biased towards a particular reference determinant. We have also compared the linear-response approach to the computation of energy differences by direct state-specific ic-MRCC calculations. We finally compare to Mk-MRCC-LR theory for which spurious roots have been reported [T.-C. Jagau and J. Gauss, J. Chem. Phys. 137, 044116 (2012)], being due to the use of sufficiency conditions to solve the Mk-MRCC equations. No such problem is present in ic-MRCC-LR theory.

Communication: Extension of a universal explicit electron correlation correction to general complete active spaces.

We present the extension of a recently proposed universal explicit electron correlation (F12) correction for multi-reference perturbation theories to general complete active spaces and arbitrary choices of complete active space self-consistent field (CASSCF) orbitals. This F12 correction is applied to Mukherjee's multi-reference second-order perturbation theory (Mk-MRPT2). Pilot examples show the expected reduction of the basis sets incompleteness error of about two cardinal numbers.

Exploration of Various Aspects of UGA-SUMRCC: Size Extensivity, Possible Use of Sufficiency Conditions, and an Extension for Direct Determination of Energy Differences.

The Unitary Group Adapted State Universal Multireference Coupled Cluster (UGA-SUMRCC) theory, recently developed by us (J. Chem. Phys.2012, 137, 074104), contains exactly the right number of linearly independent cluster operators. This avoids any redundancy of the excitation manifold in a way exactly paralleling the traditional spin-orbital based SUMRCC. The choice of the linearly independent cluster operators inducing the same change of orbital occupancy becomes increasingly cumbersome if we go over to the cases of active CSFs with more than two active quasiparticles. In the present development, we explore several aspects of the UGA-SUMRCC theory: (a) The first is a variant where we have deliberately incorporated redundancy of the cluster amplitudes to simplify the working equations and have shown that it can serve as a very good approximation to the parent UGA-SUMRCC theory for states with more than two valence occupancies. This in turn suggests that it could be a useful avenue to pursue for arbitrary mh-np situation since the working equations assume simpler algebraic structure in such cases. (b) The analyses of the aspects of size extensivity are known to involve greater complexity if they involve various reduced density matrices (RDMs), since the RDMs are not size-extensive quantities. We have presented the proof for UGA-SUMRCC starting from equations containing h-p RDMs via a decomposition involving products of size-extensive cumulants and argue that it has relevance for general cases beyond the h-p model spaces. (c) A useful extension of UGA-SUMRCC lies in formulating the theory for direct calculations of energy differences of spectroscopic interest such as excitation energies, ionization potentials, and electron affinities relative to a closed shell ground state, thus providing attractive alternatives to other allied methods such as SAC-CI, CC-LRT, EOM-CC, STEOM-CC, or ADC. This extension, called UGA-based Quasi-Fock MRCC by us, also leads to exact cancellation of common correlation terms between the initial and final states. Taking a cue from the hierarchical development in Fock-space theories but keeping in mind the advantages of a state-universal (equivalently called a valence specific) theory, our formulation proposes a spin-adapted, accurate, and compact scheme for studying such energy differences. Our results demonstrate superior performance of the method as compared to EOM-CC.

Development and applications of a unitary group adapted state specific multi-reference coupled cluster theory with internally contracted treatment of inactive double excitations.

Any multi-reference coupled cluster (MRCC) development based on the Jeziorski-Monkhorst (JM) multi-exponential ansatz for the wave-operator Ω suffers from spin-contamination problem for non-singlet states. We have very recently proposed a spin-free unitary group adapted (UGA) analogue of the JM ansatz, where the cluster operators are defined in terms of spin-free unitary generators and a normal ordered, rather than ordinary, exponential parametrization of Ω is used. A consequence of the latter choice is the emergence of the "direct term" of the MRCC equations that terminates at exactly the quartic power of the cluster amplitudes. Our UGA-MRCC ansatz has been utilized to generate both the spin-free state specific (SS) and the state universal MRCC formalisms. It is well-known that the SSMRCC theory requires suitable sufficiency conditions to resolve the redundancy of the cluster amplitudes. In this paper, we propose an alternative variant of the UGA-SSMRCC theory, where the sufficiency conditions are used for all cluster operators containing active orbitals and the single excitations with inactive orbitals, while the inactive double excitations are assumed to be independent of the model functions they act upon. The working equations for the inactive double excitations are thus derived in an internally contracted (IC) manner in the sense that the matrix elements entering the MRCC equations involve excitations from an entire combination of the model functions. We call this theory as UGA-ICID-MRCC, where ICID is the acronym for "Internally Contracted treatment of Inactive Double excitations." Since the number of such excitations are the most numerous, choosing them to be independent of the model functions will lead to very significant reduction in the number of cluster amplitudes for large active spaces, and is worth exploring. Moreover, unlike for the excitations involving active orbitals, where there is inadequate coupling between the model and the virtual functions in the SSMRCC equations generated from sufficiency conditions, our internally contracted treatment of inactive double excitations involves much more complete couplings. Numerical implementation of our formalism amply demonstrates the efficacy of the formalism.

Formulation and implementation of a unitary group adapted state universal multi-reference coupled cluster (UGA-SUMRCC) theory: excited and ionized state energies.

The traditional state universal multi-reference coupled cluster (SUMRCC) theory uses the Jeziorski-Monkhorst (JM) based Ansatz of the wave operator: Ω = Σ(µ)Ω(µ)|φ(µ)><φ(µ)|, where Ω(µ) = exp(T(µ)) is the cluster representation of the component of Ω inducing virtual excitations from the model function φ(µ). In the first formulations, φ(µ)s were chosen to be single determinants and T(µ)s were defined in terms of spinorbitals. This leads to spin-contamination for the non-singlet cases. In this paper, we propose and implement an explicitly spin-free realization of the SUMRCC theory. This method uses spin-free unitary generators in defining the cluster operators, {T(µ)}, which even at singles-doubles truncation, generates non-commuting cluster operators. We propose the use of normal-ordered exponential parameterization for Ω:Σ(µ){exp(T(µ))}|φ(µ)><φ(µ)|, where {} denotes the normal ordering with respect to a common closed shell vacuum which makes the "direct term" of the SUMRCC equations terminate at the quartic power. We choose our model functions {φ(µ)} as unitary group adapted (UGA) Gel'fand states which is why we call our theory UGA-SUMRCC. In the spirit of the original SUMRCC, we choose exactly the right number of linearly independent cluster operators in {T(µ)} such that no redundancies in the virtual functions {χ(µ)(l)} are involved. Using example applications for electron detached/attached and h-p excited states relative to a closed shell ground state we discuss how to choose the most compact and non-redundant cluster operators. Although there exists a more elaborate spin-adapted JM-like ansatz of Datta and Mukherjee (known as combinatoric open-shell CC (COS-CC), its working equations are more complex. Results are compared with those from COS-CC, equation of motion coupled cluster methods, restricted open-shell Hartree-Fock coupled cluster, and full configuration interaction. We observe that our results are more accurate with respect to most other theories as a result of the use of the cluster expansion structure for our wave operator. Our results are comparable to those from the more involved COS-CC, indicating that our theory captures the most important aspects of physics with a considerably simpler scheme.

Unitary group adapted state-specific multi-reference coupled cluster theory: formulation and pilot numerical applications.

We present the formulation and the implementation of a spin-free state-specific multi-reference coupled cluster (SSMRCC) theory, realized via the unitary group adapted (UGA) approach, using a multi-exponential type of cluster expansion of the wave-operator Ω. The cluster operators are defined in terms of spin-free unitary generators, and normal ordered exponential parametrization is utilized for cluster expansion instead of pure exponentials. Our Ansatz for Ω is a natural spin-free extension of the spinorbital based Jeziorski-Monkhorst (JM) Ansatz. The normal ordered cluster Ansatz for Ω results in a terminating series of the direct term of the MRCC equations, and it uses ordinary Wick algebra to generate the working equations in a straightforward manner. We call our formulation as UGA-SSMRCC theory. Just as in the case of the spinorbital based SSMRCC theory, there are redundancies in the cluster operators, which are exploited to ensure size-extensivity and avoidance of intruders via suitable sufficiency conditions. Although there already exists in the literature a spin-free JM-like Ansatz, introduced by Datta and Mukherjee, its structure is considerably more complex than ours. The UGA-SSMRCC offers an easier access to spin-free MRCC formulation as compared to the Datta-Mukherjee Ansatz, which at the same time provides with quite accurate description of electron correlation. We will demonstrate the efficacy of the UGA-SSMRCC formulation with a set of numerical results. For non-singlet cases, there is pronounced M(s) dependence of the energy for the spinorbital based SSMRCC results. Although M(s) = 1 results are closer to full configuration interaction (FCI), the extent of spin-contamination is more. In most of the cases, our UGA-SSMRCC results are closer to FCI than the spinorbital M(s) = 0 results.

Inactive excitations in Mukherjee's state-specific multireference coupled cluster theory treated with internal contraction: development and applications.

One generic difficulty of most state-specific many-body formalisms using the Jeziorski-Monkhorst ansatz: ψ = Σ(µ)exp(T(µ))|φ(µ)>c(µ) for the wave-operators is the large number of redundant cluster amplitudes. The number of cluster amplitudes up to a given rank is many more in number compared to the dimension of the Hilbert Space spanned by the virtual functions of up to the same rank of excitations. At the same time, all inactive excitations--though linearly independent--are far too numerous. It is well known from the success of the contracted multi-reference configuration interaction (MRCI(SD)) that, at least for the inactive double excitations, their model space dependence (µ-dependence) is weak. Considerable simplifications can thus be obtained by using a partially internally contracted description, which uses the physically appealing approximation of taking the inactive excitations T(i) to be independent of the model space labels (µ-independent). We propose and implement in this paper such a formalism with internal contractions for inactive excitations (ICI) within Mukherjee's state-specific multi-reference coupled cluster theory (SS-MRCC) framework (referred to from now on as the ICI-SS-MRCC). To the extent the µ-independence of T(i) is valid, we expect the ICI-SS-MRCC to retain the conceptual advantages of size-extensivity yet using a drastically reduced number of cluster amplitudes without sacrificing accuracy. Moreover, greater coupling is achieved between the virtual functions reached by inactive excitations as a result of the internal contraction while retaining the original coupling term for the µ-dependent excitations akin to the parent theory. Another major advantage of the ICI-SS-MRCC, unlike the other analogous internally contracted theories, such as IC-MRCISD, CASPT2, or MRMP2, is that it can use relaxed coefficients for the model functions. However, at the same time it employs projection manifolds for the virtuals obtained from inactive n hole-n particle (nh-np) excitations on the entire reference function containing relaxed model space coefficients. The performance of the method has been assessed by applying it to compute the potential energy surfaces of the prototypical H(4); to the torsional potential energy barrier for the cis-trans isomerism in C(2)H(4) as well as that of N(2)H(2), automerization of cyclobutadiene, single point energy calculation of CH(2), SiH(2), and comparing them against the SS-MRCC results, benchmark full CI results, wherever available and those from the allied MR formalisms. Our findings are very much reminiscent of the experience gained from the IC-MRCISD method.

A spin-adapted size-extensive state-specific multi-reference perturbation theory. I. Formal developments.

We present in this paper a comprehensive formulation of a spin-adapted size-extensive state-specific multi-reference second-order perturbation theory (SA-SSMRPT2) as a tool for applications to molecular states of arbitrary complexity and generality. The perturbative theory emerges in the development as a result of a physically appealing quasi-linearization of a rigorously size-extensive state-specific multi-reference coupled cluster (SSMRCC) formalism [U. S. Mahapatra, B. Datta, and D. Mukherjee, J. Chem. Phys. 110, 6171 (1999)]. The formulation is intruder-free as long as the state-energy is energetically well-separated from the virtual functions. SA-SSMRPT2 works with a complete active space (CAS), and treats each of the model space functions on the same footing. This thus has the twin advantages of being capable of handling varying degrees of quasi-degeneracy and of ensuring size-extensivity. This strategy is attractive in terms of the applicability to bigger systems. A very desirable property of the parent SSMRCC theory is the explicit maintenance of size-extensivity under a variety of approximations of the working equations. We show how to generate both the Rayleigh-Schrödinger (RS) and the Brillouin-Wigner (BW) versions of SA-SSMRPT2. Unlike the traditional naive formulations, both the RS and the BW variants are manifestly size-extensive and both share the avoidance of intruders in the same manner as the parent SSMRCC. We discuss the various features of the RS as well as the BW version using several partitioning strategies of the hamiltonian. Unlike the other CAS based MRPTs, the SA-SSMRPT2 is intrinsically flexible in the sense that it is constructed in a manner that it can relax the coefficients of the reference function, or keep the coefficients frozen if we so desire. We delineate the issues pertaining to the spin-adaptation of the working equations of the SA-SSMRPT2, starting from SSMRCC, which would allow us to incorporate essentially any type open-shell configuration-state functions (CSF) within the CAS. The formalisms presented here will be applied extensively in a companion paper to assess their efficacy.

A spin-adapted size-extensive state-specific multi-reference perturbation theory with various partitioning schemes. II. Molecular applications.

Following the theoretical development of a spin-adapted state-specific multi-reference second-order perturbation theory (SA-SSMRPT2) as expounded in the preceding publication, we discuss here its implementation and the results of its applications to potential energy curves (PECs) of various electronic states of small molecules. In particular, we illustrate its efficacy in states of various spin multiplicities and varying multi-reference character. Both Møller-Plesset (MP) and Epstein-Nesbet (EN) type of partitions have been explored. Also, a straightforward Rayleigh-Schrödinger (RS) and Brillouin-Wigner (BW) version of the SA-SSMRPT2 have been studied. Ground state PECs were computed for singlet states of HF, BH, and H(2)O molecules as well as the doublet state of NH(2) and BeH radicals and compared to corresponding full configuration interaction numbers, which serve as benchmark results. As an extensive application on a production level, the ground state PECs of N(2), a classic example of multiple-bond breaking, were calculated using cc-pVXZ (X = 3,4,5) basis and then extrapolated to obtain estimates of the complete basis set limit. Vibrational energy levels were extracted from these N(2) PECs, which compare favorably to the experimental values. In addition, extensive studies were also carried out on PECs of the seven low-lying excited states of the N(2) molecule. Finally, it is shown that the flexibility to relax configuration coefficients in SA-SSMRPT2 helps to provide good descriptions for the avoided crossing between the two lowest (1)Σ states of the LiF molecule. Our results indicate (1) that more studies are needed to draw firm conclusions about the relative efficacies of the MP and EN results and (2) that the RS version works so well as compared to the BW version that the extra computational expenses needed in the later formalism is not warranted.

The spin-free analogue of Mukherjee's state-specific multireference coupled cluster theory.

In this paper, we develop a rigorously spin-adapted version of Mukherjee's state-specific multireference coupled cluster theory (SS-MRCC, also known as Mk-MRCC) [U. S. Mahapatra, B. Datta, and D. Mukherjee, J. Chem. Phys. 110, 6171 (1999)] for reference spaces comprising open-shell configurations. The principal features of our approach are as follows: (1) The wave operator Ω is written as Ω = ∑(µ)Ω(µ)|φ(µ)>c(µ), where {φ(µ)} is the set of configuration state functions spanning a complete active space. (2) In contrast to the Jeziorski-Monkhorst Ansatz in spin-orbital basis, we write Ω(µ) as a power series expansion of cluster operators R(µ) defined in terms of spin-free unitary generators. (3) The operators R(µ) are either closed-shell-like n hole-n particle excitations (denoted as T(µ)) or they involve valence (active) destruction operators (denoted as S(µ)); these latter type of operators can have active-active scatterings, which can also carry the same active orbital labels (such S(µ)'s are called to have spectator excitations). (4) To simulate multiple excitations involving powers of cluster operators, we allow the S(µ)'s carrying the same active orbital labels to contract among themselves. (5) We exclude S(µ)'s with direct spectator scatterings. (6) Most crucially, the factors associated with contracted composites are chosen as the inverse of the number of ways the S(µ)'s can be joined among one another leading to the same excitation. The factors introduced in (6) have been called the automorphic factors by us. One principal thrust of this paper is to show that the use of the automorphic factors imparts a remarkable simplicity to the final amplitude equations: the equations consist of terms that are at most quartic in cluster amplitudes, barring only a few. In close analogy to the Mk-MRCC theory, the inherent linear dependence of the cluster amplitudes leading to redundancy is resolved by invoking sufficiency conditions, which are exact spin-free analogues of the spin-orbital based Mk-MRCC theory. This leads to manifest size-extensivity and an intruder-free formulation. Our formalism provides a relaxed description of the nondynamical correlation in presence of dynamical correlation. Pilot numerical applications to doublet systems, e.g., potential energy surfaces for the first two excited (2)A' states of asymmetric H(2)S(+) ion and the ground (2)Σ(+)state of BeH radical are presented to assess the viability of our formalism over an wide range of nuclear geometries and the manifest avoidance of intruder state problem.

Inclusion of selected higher excitations involving active orbitals in the state-specific multireference coupled-cluster theory.

The parent state-specific multireference coupled-cluster (SS-MRCC) theory proposed by Mukherjee et al. [J. Chem. Phys. 110, 6171 (1999)], though rigorously size-extensive and also size-consistent with localized orbitals, has some deficiencies in the minimal truncation scheme, viz. at the singles and doubles (SD) level (SS-MRCCSD). SS-MRCCSD does not involve the direct coupling of all the model functions with a given virtual function belonging to the uncontracted multiconfiguration CISD space. It also does not involve, even in the linear power of a cluster operator T(µ), the direct coupling of the virtual functions χ(l(µ)), which are up to doubly excited with respect to a model function φ(µ) to the other virtual functions of the MRCISD space which can be generated by triple and quadruple excitations from φ(µ). We argue that inclusion of a selection of triples and quadruples involving at most two inactive orbital excitations from every φ(µ) would ameliorate the shortcoming of the incomplete coupling of the triply and quadruply excited virtual functions which can couple with the singly and doubly excited ones. This extended ansatz for our SS-MRCC theory, to be called SS-MRCCSDtq by us, would still miss the direct coupling of the manifold of the model functions {φ(λ),λ ≠ µ} to singly and doubly excited virtual functions. However, this effect is expected to be less significant than the lack of the more complete virtual space couplings, these functions being many more numerous, suggesting the new methods to be significantly improved schemes. Excellent results on the potential energy surfaces of small molecules involving single, double, and triple bond dissociation bear out our expectations fully.

An algebraic proof of generalized Wick theorem.

The multireference normal order theory, introduced by Kutzelnigg and Mukherjee [J. Chem. Phys. 107, 432 (1997)], is defined explicitly, and an algebraic proof is given for the corresponding contraction rules for a product of any two normal ordered operators. The proof does not require that the contractions be cumulants, so it is less restricted. In addition, it follows from the proof that the normal order theory and corresponding contraction rules hold equally well if the contractions are only defined up to a certain level. These relaxations enable us to extend the original normal order theory. As a particular example, a quasi-normal-order theory is developed, in which only one-body contractions are present. These contractions are based on the one-particle reduced density matrix.