Multi-reference coupled-cluster linear

Debashis Mukherjee is a Professor of Physical Chemistry and the Director of the Indian Association for the Cultivation of Science, Calcutta, India. He has been one of the earliest developers of a class of multi-reference coupled cluster theories and also of the coupled cluster based linear response theory. Other contributions by him are in the resolution of the size-extensivity problem for multi-reference theories using an incomplete model space and in the size-extensive intermediate Hamiltonian formalism. His research interests focus on the development and applications of non-relativistic and relativistic theories of many-body molecular electronic structure and theoretical spectroscopy, quantum many-body dynamics and statistical held theory of many-body systems. He is a member of the International Academy of the Quantum Molecular Science, a Fellow of the Third World Academy of Science, the Indian National Science Academy and the Indian Academy of Sciences. He is the recipient of the Shantiswarup Bhatnagar Prize of the Council of Scientihc and Industrial Research of the Government of India. 

Although the past 20 years have witnessed a great progress in the Hilbert space multi-reference coupled cluster methods (see, for example, the work of Mukherjee and Pal [99],Paldus [101], Jeziorski and Paldus [102], Jankowski et al. [103],Paldus et al. [104], Paldus et al. [105], Meissner et al. [106], Kucharski and Bartlett [107], Balkovd et al. [108], Baikova and Bartlett [109], Balkovd et al. [110], Baikova et al. [Ill], Berkovic and Kaldor [112]) only a few applications of this approach have been reported, mostly oriented to the simple model systems exploiting a lowdimensional model space. Among the reasons for this paucity of applications are the choice of an appropriate model space, convergence difficulties arising from intruder state problems and from multiple solutions of non-linear coupled cluster equations. 

Another category of approaches that avoids the symmetry breaking problem of the orbitals for the target state is based on using a related, well-behaved HF reference instead. Examples of such techniques include quasi-restricted Hartree-Fock coupled-cluster (QRHF CC) (11,19), symmetry adapted cluster configuration interaction (SAC-CI) (22,23), coupled-cluster linear response theory (CCLRT) (24-26) or equivalently equation-of-motion coupled-cluster (EOM-CC) (27-32), Fock space multi-reference coupled-cluster (FSMRCC) (33-37), and similarity transformed equation-of-motion coupled-cluster (STEOM-CC) (38-40). In the case of NO3 and N03, the restricted Hartree-Fock (RHF) orbitals of the closed-shell N03 anion system can be used as the reference. The anion orbitals are stable with respect to symmetry perturbations, and the system does not suffer from the artifactual symmetry breaking of the neutral and cation. 

Abstract The modified equation-of-motion coupled cluster approach of Sekino and Bartlett is extended to computations of the mixed electric-dipole/magnetic-dipole polarizability tensor associated with optical rotation in chiral systems. The approach - referred to here as a linearized equation-of-motion coupled cluster (EOM-CCl) method - is a compromise between the standard EOM method and its linear response counterpart, which avoids the evaluation of computationally expensive terms that are quadratic in the field-perturbed wave functions, but still yields properties that are size-extensive/intensive. Benchmark computations on five representative chiral molecules, including (P)-hydrogen peroxide, (5)-methyloxirane, (5 )-2-chloropropioniuile, (/ )-epichlorohydrin, and (75,45)-norbornenone, demonstrate typically small deviations between the EOM-CCl results and those from coupled cluster linear response theory, and no variation in the signs of the predicted rotations. In addition, the EOM-CCl approach is found to reduce the overall computing time for multi-wavelength-specific rotation computations by up to 34%. 

AIMD = ab initio molecular dynamics B-LYP = Becke-Lee-Yang-Parr CCSD = coupled cluster single double excitations CVC = core-valence correlation ECP = effective core potential DF = density functional GDA = gradient corrected density approximation MCLR = multiconfigurational linear response MP2 = M0ller-Plesset second-order (MRD)CI = multi-reference double-excitation configuration interaction RPA = random phase approximation TD-MCSCF = time-dependent multiconfigurational self-consistent field TD-SCF = time-dependent self-consistent field. 

PHF methods can, in turn, be classified as the variational and nonvariational ones. In the former gronp of methods the coefficients in linear combination of Slater determinants and in some cases LCAO coefficients in HF MOs are optimized in the PHF calculations, in the latter such an optimization is absent. To the former group of PHF methods one refers different versions of the configuration interaction (Cl) method, the multi-configuration self-consistent field (MCSCF) method, the variational coupled cluster (CC) approach and the rarely used valence bond (VB) and generaUzed VB methods. The nonvariational PHF methods inclnde the majority of CC reaUza-tions and many-body perturbation theory (MBPT), called in its molecular realization the MoUer-Plessett (MP) method. In MP calculations not only RHF but UHF MOs are also used [107]. 

---