Highly Accurate Expectation Values Using High-Order Relativistic Coupled Cluster Theory

This work presents the automatic generation of analytic first derivatives of the energy for general coupled-cluster models using the \text{tenpi} toolchain. We report the first implementation of expectation values for CCSDT and CCSDTQ methods within the DIRAC program package for relativistic molecular calculations. As pivotal calculations, we focus on the electric field gradient (EFG) evaluated at the lithium nucleus in LiX (X = H, F, Cl) compounds, enabling the extraction of the nuclear electric quadrupole moment $Q({}^{7}Li)$, and at the aluminum nucleus in AlY (Y =H, F, Cl, Br) compounds, for the determination of $Q({}^{27}Al)$. These high-order methods are applied to compute corrections for triple and quadruple excitations for the EFG, a crucial quantity for determining nuclear quadrupole moments. We obtain $Q({}^{27}Al)$ = 0.1466 b, in excellent agreement with the recommended value, and $Q({}^{7}Li)$ = -0.0386 b, which is smaller than the currently recommended value, that indicates the need for further investigation.