Entropy production and statistical relaxation of dipolar bosons and fermions in interaction quench dynamics

We study the out-of-equilibrium dynamics of dipolar bosons and fermions after a sudden change in the interaction strength from zero to a finite repulsive value. We simulate the interaction quench on the initial state which is the ground state of harmonic potential with noninteracting bosons and fermions. We solve the time-dependent many-boson Schrödinger equation exactly using numerical methods. To understand the many-body dynamics we analyze several measures of many-body information entropy, monitoring their time evolution and assessing their dependence on interaction strength. We establish that for weak interaction quench the dynamics is statistics independent, both dipolar bosons and fermions do not relax, whereas it is significantly different for dipolar bosons from that of dipolar fermions in the stronger interaction quench. While dipolar bosons exhibit concurrent signature of relaxation in all entropy measures, dipolar fermions fail to relax, exhibiting modulated oscillation in all entropy dynamics. For dipolar bosons and larger interaction quench, the many-body information entropy measures dynamically approach the value predicted for the Gaussian orthogonal ensemble of random matrices, implying statistical relaxation and ergodicity. The entire relaxation process and signature of ergodicity are established by three unique features: dynamical delocalization in Hilbert space, violent orbital fragmentation, and saturation of all entropy measures to the maximum entropy states in the long-time dynamics. It highlights the importance of many-body effects with a possible exploration in quantum simulation with ultracold atoms.