The memory kernel for a tagged particle in a fluid, computed from molecular dynamics simulations, decays algebraically as t −3/2 . We show how the hydrodynamic Basset-Boussinesq force naturally emerges from this long-time tail and generalize the concept of hydrodynamic added mass. This mass term is negative in the present case of a molecular solute, at odds with incompressible hydrodynamics predictions. We finally discuss the various contributions to the friction, the associated time scales and the cross-over between the molecular and hydrodynamic regimes upon increasing the solute radius.The Brownian motion of a particle in a fluid finds its origin in the fluctuating force exerted by the solvent molecules on the solute. It has long been known that the canonical description of this random force by a Gaussian Markov process is only valid in limiting cases. Even in the limit where the solute is much heavier than the solvent particles, for which multiple time-scale analysis allows to recover the Smoluchowski equation for diffusion [1], nonMarkovian effects are expected when the mass density ratio is close to unity [2] -a situation which is rather the rule than the exception e.g. in colloidal suspensions. These non-Markovian effects arise because of momentum conservation, leading to slow hydrodynamic modes that manifest themselves as long-time tails in the velocity autocorrelation function (VACF) [3][4][5][6]. Recent experiments have demonstrated that the force exerted by the bath includes a deterministic component [7], well described for large colloidal spheres by the Basset-Boussinesq (BB) hydrodynamic force [8,9]:where R is the sphere radius, η the solvent viscosity and ρ 0 its mass density. The first term is the usual Stokes friction. The other two account for the inertia of the displaced fluid and involve a finite added mass m BB 0 = 2 3 πR 3 ρ 0 and a viscosity-dependent retarded component describing the transient effects of momentum diffusion in the solvent. While continuous descriptions of steady-state flows appear to hold down to the nanoscale [10][11][12], possibly at the price of adapting the hydrodynamic radius R or the boundary conditions [13], their validity for the transient regimes should be questioned. The implicit assumption of a separation of time scales between the solvent and solute dynamics, which holds a priori for colloidal particles [14], is expected to break down with smaller solutes such as nanoparticles or biomolecules.Here we address the fundamental questions that arise when approaching the regime of molecular solutes by computing directly from Molecular Dynamics (MD) simulations the memory kernel and the random noise of the Generalized Langevin Equation (GLE). A novel algorithm based on the Mori-Zwanzig formalism with high numerical stability allows us to explore long time scales for the first time. We consider the extreme case of a tagged particle (identical masses and sizes) in a pure supercritical fluid.By examining the long-time behaviour of the memory kernel, we demonstrate the gen...