Here, i G N is the number of trajectories in cluster i. Note that in multidimensional cases, we can use a direct product of Gaussian functions for all the nuclear DoFs and introduce the covariant matrix to calculate the variance. The nuclear density is reconstructed by summing all the Gaussians with the corresponding weights as ( ) ( )