It is believed that an isolated and far-from-equilibrium quantum many-body system should try to attain equilibrium via a mechanism whereby any given subsystem acts as an open quantum system that is coupled to an environment, which is the complementary part of the full system, and undergoes a complicated equilibration process such that all the subsystems in the long-time limit attain equilibrium states compatible with the global equilibrium state. This picture begs the question whether the dynamics of any given subsystem is Markovian (monotonic loss of information and memory) or non-Markovian. In this work, by numerically probing the dynamical behaviour of the quantum distances between
temporally-separated states of small subsystems, we reveal the telltale signatures of (non-)Markovianity of the dynamics of subsystems of an isolated quantum spin system brought in the far-from-equilibrium regime, exemplified with the mixed-field Ising spin chain quenched between parameter regimes deep inside its magnetically ordered and disordered regimes. Additionally, remarkably systematic behaviour is seen in a measure of classical distances between the quantum states of the considered subsystems. These features strongly depend on the direction of quenching in the parameter space, with paramagnetic-to-ferromagnetic quenches offering considerably stronger signatures of subsystem non-Markovianity, for which we offer heuristic arguments.
Quantum systems are invariably open, evolving under surrounding influences rather than in isolation. Standard open quantum system methods eliminate all information on the environmental state to yield a tractable description of the system dynamics. By incorporating a collective coordinate of the environment into the system Hamiltonian, we circumvent this limitation. Our theory provides straightforward access to important environmental properties that would otherwise be obscured, allowing us to quantify the evolving system-environment correlations. As a direct result, we show that the generation of robust system-environment correlations that persist into equilibrium (heralded also by the emergence of non-Gaussian environmental states) renders the canonical system steady-state almost always incorrect. The resulting equilibrium states deviate markedly from those predicted by standard perturbative techniques and are instead fully characterised by thermal states of the mapped system-collective coordinate Hamiltonian. We outline how noncanonical system states could be investigated experimentally to study deviations from canonical thermodynamics, with direct relevance to molecular and solid-state nanosystems.