Transport phenomena are central to physics, and transport in the many-body and fully-quantum regime is attracting an increasing amount of attention. It has been recently revealed that some quantum spin chains support ballistic transport of excitations at all energies. However, when joining two semi-infinite ballistic parts, such as the XX and XXZ spin-1/2 models, our understanding suddenly becomes less established. Employing a matrix-product-state ansatz of the wavefunction, we study the relaxation dynamics in this latter case. Here we show that it takes place inside a light cone, within which two qualitatively different regions coexist: an inner one with a strong tendency towards thermalization, and an outer one supporting ballistic transport. We comment on the possibility that even at infinite time the system supports stationary currents and displays a non-zero Kapitza boundary resistance. Our study paves the way to the analysis of the interplay between transport, integrability, and local defects.In 1941, P. Kapitza reported on the first measurement of the temperature drop near the boundary between helium and a solid when heat flows across the boundary [1]. The phenomenon, unrelated to that of conventional contact resistances, is ascribed to the mismatch between the energy carriers of the two materials, and appears even if the interface is perfect at the atomic scale. In the years 1951-1956, Kalatchnikov first, and Mazo and Onsager later, independently developed the so-called acoustic mismatch model, that gives a mathematical formulation to such intuition [2,3]. The quantitative comparison between experimental data and theoretical predictions motivated extensive investigations even in other interfaces without helium, e.g., between two solids [2], where the phenomenon has been often recovered.Owing to the difficulty of performing numerical simulations in order to benchmark the theory, simpler treatable models have been considered. In particular, focusing on junctions of classical one-dimensional (1D) harmonic chains enabled to quantitatively investigate the phenomenon of thermal boundary resistances [4][5][6]. The effect has been widely reproduced, and comparisons with suitable adaptations of the theory have been presented. Yet, a numerical investigation of the phenomenon in the fully quantum case has not been performed so far [7].The recent advances in the analytical understanding of the non-equilibrium dynamics of quantum spin chains [8][9][10][11][12] stimulated the improvement of state-of-theart numerical simulations addressing the unitary timeevolution of strongly-correlated 1D quantum models [13- * alberto.biella@univ-paris-diderot.fr 15]. Groundbreaking cold-atom experiments have also investigated several aspects of the coherent dynamics of closed quantum many-body systems [16][17][18][19][20][21][22][23][24]. These research efforts allowed to faithfully address how the presence of impurities and inhomogeneities may affect the quantum transport phenomena, also far from the linearresponse regime [25][26][2...