To achieve this, the simulation box undergoes deformation in one lateral direction by adjusting the angle of the x-y plane relative to the x direction. The deformation rate equals the selected shear rate multiplied by the length of the simulation box in the y direction. Additionally, sufficient time is allotted for the simulation to converge towards its viscosity value. In this model, we utilized 7×10 7 timesteps, each of size ∆t = 0.001 µs. The viscosity was derived by averaging the instantaneous viscosity values after the simulation convergence. Each simulation required approximately two days utilizing one node.