In this way, provided that the CG-MD trajectory is well converged, we can directly calculate the free energy surface (FES) of monomer states within the fibers based on the probability distribution in feature space (F = -K B T log(P)), using the two largest principal components obtained from the combined set of fibers.