Various molecular parameters in quantum chemistry could be computed as derivatives of energy over different arguments. Unfortunately, it is quite complicated to obtain analytical expression for characteristics that are of interest in the framework of methods that account electron correlation. Especially it relates to the coupled cluster (CC) theory. In such cases, numerical differentiation comes to rescue. This approach, like any other numerical method has empirical parameters and restrictions that require investigation. Current work is called to clarify the details of Finite-Field method usage for high-order derivatives calculation in CC approaches. General approach to the parameter choice and corresponding recommendations about numerical steadiness verification are proposed. As an example of Finite-Field approach implementation characterization of optical properties of fullerene passing process through the aperture of carbon nanotorus is given.
Accuracy and Stability of Computing High-order Derivatives of Analytic Functions by Cauchy Integrals