Perturbative triples correction to domain-based local pair natural orbital variants of Mukherjee's state specific coupled cluster method

In this article we report an implementation of the perturbative triples correction to Mukherjee's state-specific multireference coupled cluster method based on the domain-based pair natural orbital approach (DLPNO-MkCC). We tested the performance of DLPNO-MkCCSD(T) in calculations involving tetramethyleneethane and isomers of naphthynes. These tests show that more than 97% of triples energy was recovered with respect to the canonical MkCCSD(T) method, which together with the DLPNO-MkCCSD part accounts for about 99.70-99.85% of the total correlation energy. The applicability of the method was demonstrated on calculations of singlet-triplet gaps for several large systems: triangulene, dynemicin A, and a beryllium complex.

Spatial and Spin Symmetry Breaking in Semidefinite-Programming-Based Hartree-Fock Theory

The Hartree-Fock problem was recently recast as a semidefinite optimization over the space of rank-constrained two-body reduced-density matrices (RDMs) [ Phys. Rev. A 2014 , 89 , 010502(R) ]. This formulation of the problem transfers the nonconvexity of the Hartree-Fock energy functional to the rank constraint on the two-body RDM. We consider an equivalent optimization over the space of positive semidefinite one-electron RDMs (1-RDMs) that retains the nonconvexity of the Hartree-Fock energy expression. The optimized 1-RDM satisfies ensemble N-representability conditions, and ensemble spin-state conditions may be imposed as well. The spin-state conditions place additional linear and nonlinear constraints on the 1-RDM. We apply this RDM-based approach to several molecular systems and explore its spatial (point group) and spin ( Ŝ² and Ŝ₃) symmetry breaking properties. When imposing Ŝ² and Ŝ₃ symmetry but relaxing point group symmetry, the procedure often locates spatial-symmetry-broken solutions that are difficult to identify using standard algorithms. For example, the RDM-based approach yields a smooth, spatial-symmetry-broken potential energy curve for the well-known Be-H₂ insertion pathway. We also demonstrate numerically that, upon relaxation of Ŝ² and Ŝ₃ symmetry constraints, the RDM-based approach is equivalent to real-valued generalized Hartree-Fock theory.

Relativistic state-specific multireference perturbation theory incorporating improved virtual orbitals: Application to the ground state single-bond dissociation

Using four-component (4c) relativistic spinors, we present a computationally economical relativistic ab initio method for molecular systems employing our recently proposed second-order state-specific multireference perturbation theory (SSMRPT) incorporating the improved virtual orbital-complete active space configuration interaction (IVO-CASCI) reference wavefunction. The resulting method, 4c-IVO-SSMRPT [calculate one state at a time] is tested in pilot calculations on the homonuclear dimers including Li(2), Na(2), K(2), Rb(2), F(2), Cl(2), and Br(2) through the computations of the ground state potential energy curves (PECs). As SSMRPT curbs intruder effects, 4c-IVO-SSMRPT is numerically stable. To our knowledge, the SSMRPT in the 4c relativistic framework has not been explored in the past. Selective spectroscopic constants that are closely related to the correct shape and accuracy of the energy surfaces have been extracted from the computed PECs. For the halogen molecules, a relativistic destabilization of the bond has been found. Relativistic and electron correlation effects need to be incorporated to get reliable estimates. Our results are in good accordance with reference theoretical and experimental data which manifests the computational accuracy and efficiency of the new 4c-IVO-SSMRPT method. The method opens for an improved description of MR systems containing heavy elements. The inexpensiveness of IVO-CASCI makes 4c-IVO-SSMRPT method promising for studies on large systems of heavy elements.

Iterative universal state selective correction for the Brillouin-Wigner multireference coupled-cluster theory

As a further development of the previously introduced a posteriori Universal State-Selective (USS) corrections [K. Kowalski, J. Chem. Phys. 134, 194107 (2011); J. Brabec et al., ibid. 136, 124102 (2012)], we suggest an iterative form of the USS correction by means of correcting effective Hamiltonian matrix elements. We also formulate USS corrections via the left Bloch equations. The convergence of the USS corrections with excitation level towards the full configuration interaction (FCI) limit is also investigated. Various forms of the USS and simplified diagonal USS corrections at the singles and doubles and perturbative triple levels are numerically assessed on several model systems and on the ozone and tetramethyleneethane molecules. It is shown that the iterative USS correction can successfully replace the previously developed a posteriori Brillouin-Wigner coupled cluster size-extensivity correction, while it is not sensitive to intruder states and performs well also in other cases when the a posteriori one fails, like, e.g., for the asymmetric vibration mode of ozone.