The simulation of open molecular systems requires explicit or implicit reservoirs of energy and particles. Whereas full atomistic resolution is desired in the region of interest, there is some freedom in the implementation of the reservoirs. Here, a combined, explicit reservoir is constructed by interfacing the atomistic region with regions of point-like, non-interacting particles (tracers) embedded in a thermodynamic mean field. The tracer molecules acquire atomistic resolution upon entering the atomistic region and equilibrate with this environment, while atomistic molecules become tracers governed by an effective mean-field potential after crossing the atomistic boundary. The approach is extensively tested on thermodynamic, structural, and dynamic properties of liquid water. Conceptual and numerical advantages of the procedure as well as new perspectives are highlighted and discussed.Molecular dynamics (MD) simulations have proven as a powerful means of investigation of molecular systems for half a century now. [1] The steady increase of computing resources has brought larger system sizes and longer time scales into the scope of direct simulations, which may soon be run in parallel to experiments. Diverse application fields know about situations where the region of interest exchanges particles with an environment, for example, droplet evaporation, [2] condensation in porous hosts, [3,4] nanoflows, [5] and chemical reactions in