The initial source fluctuates event-by-event (ebe), leading to the ebe fluctuation of the harmonic flow coefficients in the final state [1]}, [2]}. These flow fluctuations can be probed by studying the correlation coefficient between the harmonic flow vectors in two different transverse momentum bins [3]}, the factorization-breaking coefficients. The fluctuations result in the decorrelation (deviation of the correlation coefficients from 1) between these flow vectors, which manifests the combined effect of flow magnitude decorrelation as well as the flow angle decorrelation[4]}, [5]}. These two effects can be separated by constructing four particle correlators between two transverse momentum bins. However, to avoid the experimental difficulty due to the low statistics at higher momenta, the coefficients are constructed keeping two particles fixed at transverse momentum bins (differential flow) and the other two are momentum averaged (integrated flow)[6]}. We compare our model results with the preliminary results by ALICE collaboration. We also study the similar factorization coefficients between the mixed flow harmonics e.g. \(V_2^2\)  - \(V_4(p)\)  or \(V_2 V_3\)  - \(V_5(p)\)  , which serve as measures of the non-linear hydrodynamic response of the medium.
