Frozen natural orbital coupled-cluster theory: Forces and application to decomposition of nitroethane

The frozen natural orbital (FNO) coupled-cluster method increases the speed of coupled-cluster (CC) calculations by an order of magnitude with no consequential error along a potential energy surface. This method allows the virtual space of a correlated calculation to be reduced by about half, significantly reducing the time spent performing the coupled-cluster (CC) calculation. This paper reports the derivation and implementation of analytical gradients for FNO-CC, including all orbital relaxation for both noncanonical and semicanonical perturbed orbitals. These derivatives introduce several new orbital relaxation contributions to the CC density matrices. FNO-CCSD(T) and FNO-LambdaCCSD(T) are applied to a test set of equilibrium structures, verifying that these methods are capable of reproducing geometries and vibrational frequencies accurately, as well as energies. Several decomposition pathways of nitroethane are investigated using CCSD(T) and LambdaCCSD(T) with 60% of the FNO virtual orbitals in a cc-pVTZ basis, and find differences on the order of 5 kcalmol with reordering of the transition state energies when compared to B3LYP 6-311 + G(3df, 2p).