A perturbative triples correction to relativistic Quadratic Unitary Coupled Cluster Method: Theory, Implementation and Benchmarking

We present a perturbative triples correction to the relativistic quadratic unitary coupled cluster singles and doubles (qUCCSD) method, denoted as qUCCSD[T]. The method builds upon the Hermitian structure of the unitary ansatz and employs a many-body perturbation theory framework to consistently include the effects of triple excitations without the need for computationally intensive iterative procedures. Relativistic effects are incorporated using the exact two-component atomic mean-field (X2CAMF) Hamiltonian, and the computational cost is further reduced through the frozen natural spinor (FNS) and Cholesky decomposition (CD) approximations. Benchmark results demonstrate that qUCCSD[T] outperforms previously proposed triples corrections to the unitary coupled cluster method in the clasical computing regime and yields excellent agreement with experimental data and Full CI benchmarks. Specifically, the method shows high accuracy in computing bond dissociation enthalpies, molecular geometries, vibrational frequencies, ionization potentials, and electron affinities of heavy-element-containing systems.