We study the remanent magnetization in antiferromagnetic, many-body localized quantum spin chains, initialized in a fully magnetized state. Its long time limit is an order parameter for the localization transition, which is readily accessible by standard experimental probes in magnets. We analytically calculate its value in the strong-disorder regime exploiting the explicit construction of quasi-local conserved quantities of the localized phase. We discuss analogies in cold atomic systems.Introduction. The non-equilibrium dynamics in disordered, isolated quantum systems have been subject to theoretical investigations ever since the notion of localization was introduced in [1]. Spin systems in random fields are prototypical models to analyze the disorderinduced breakdown of thermalization: a large number of studies on disordered spin chains [2][3][4][5][6][7][8][9][10][11][12][13] has provided evidence for a dynamical phase transition between a weakdisorder phase which thermalizes, and a Many-Body Localized (MBL) phase in which excitations do not diffuse, ergodicity is broken and local memory of the initial conditions persists for infinite time [14][15][16].Signatures of MBL are found in the properties of individual many-body eigenstates. Even highly excited eigenstates exhibit area-law scaling of the bipartite entanglement entropy [4,7,8,17] and Poissonian level statistics [2,18,19], both being incompatible with thermalization [20][21][22]. Novel dynamical properties such as the logarithmic spreading of entanglement have been observed in direct simulations of the time evolution [3,23]. The non-equilibrium physics of MBL systems has been probed experimentally in artificial quantum systems made of cold atomic gases [24,25] and trapped ion systems [26], while an indirect signature in the from of strongly suppressed absorption of radiation was found in electron-glasses [27]. However, direct observations of MBL in solid-state materials are still lacking.It has been argued [28][29][30] that the properties of MBL systems are related to the existence of extensively many quasi-local conserved operators that strongly constrain the quantum dynamics, preventing both transport and thermalization. Their existence also follows as a corollary from Imbrie's rigorous arguments in favor of MBL [31,32].In this work, we propose a experimentally readily observable consequence of MBL in quantum magnets: the out-of-equilibrium remanent magnetization that persists after ferromagnetically polarizing an antiferromagnet whose total magnetization is not a conserved. The remanence implies non-ergodicity, since ergodic dynamics would relax the magnetization completely (cf. Fig. 1 for a schematic sketch of the protocol). As an example, we consider an antiferromagnetic, anisotropic Heisenberg spin-1/2 chain(1) subject to random fields h k along the Ising axis. We assume J z < 0, as well as J x = J y to ensure the nonconservation of the total magnetization. Such Hamiltonians can be realized, e.g., in Ising compounds with both exchange and dipola...