A Synergistic Approach towards Optimization of Coupled Cluster Amplitudes by Exploiting Dynamical Hierarchy

The coupled cluster iteration scheme for determining the cluster amplitudes involves a set of nonlinearly coupled difference equations. In the space spanned by the amplitudes, the set of equations are analyzed as a multivariate time-discrete map where the concept of time appears in an implicit manner. With the observation that the cluster amplitudes have difference in their relaxation timescales with respect to the distributions of their magnitudes, the coupled cluster iteration dynamics are considered as a synergistic motion of coexisting slow and fast relaxing modes, manifesting a dynamical hierarchical structure. With the identification of the highly damped auxiliary amplitudes, their time variation can be neglected compared to the principal amplitudes which take much longer time to reach the fixed points. We analytically establish the adiabatic approximation where each of these auxiliary amplitudes are expressed as unique parametric functions of the collective principal amplitudes, allowing us to study the optimization with the latter taken as the independent degrees of freedom. Such decoupling of the amplitudes significantly reduces the computational scaling without sacrificing the accuracy in the ground state energy as demonstrated by a number of challenging molecular applications. A road-map to treat higher order post-adiabatic effects is also discussed.