where, O N is the normal-ordered property operator, i.e., the derivative of normal-ordered perturbed Hamiltonian with respect to an external perturbation. On the other hand, a desired property can be simply computed as the expectation-value of the property operator of interest. Once the cluster amplitudes are solved, the expectation-value of O N can be calculated using the CC wavefunction as given below (Ref ):