Dynamically emergent correlations in bosons via quantum resetting

We study the nonequilibrium stationary state (NESS) induced by quantum resetting of a system of N noninteracting bosons in a harmonic trap. Our protocol consists of preparing initially the system in the ground state of a harmonic oscillator centered at +a, followed by a rapid quench where the center is shifted to -a and the system is allowed to evolve unitarily up to a random Poissonian time distributed via r\, e^-r\,. Then the trap center is reset to +a again and the system is assumed to cool instantaneously to the initial ground state. The system is again allowed to evolve unitarily in the trap centered at -a up to a random time, and the procedure is repeated. Under repeated resetting, the system reaches a NESS where the positions of bosons get strongly correlated due to simultaneous resetting induced by the trap. We fully characterize the steady state by analytically computing several physical observables such as the average density, extreme value statistics, order and gap statistics, and also the distribution of the number of particles in a region -L, L, known as the full counting statistics (FCS). In particular, we show that in the large N limit, the scaling function describing the FCS exhibits a striking feature: it is supported over a nontrivial finite interval, and moreover is discontinuous at an interior point of the support. Our results are supported by numerical simulations. This is a rare example of a strongly correlated quantum many-body NESS where various observables can be exactly computed.