To prevent the flocculation and phase separation of nanoparticles in solution, nanoparticles are often functionalized with short chain surfactants. Here we present fully-atomistic molecular dynamics simulations which characterize how these functional coatings affect the interactions between nanoparticles and with the surrounding solvent. For 5 nm diameter silica nanoparticles coated with poly(ethylene oxide) (PEO) oligomers in water, we determined the hydrodynamic drag on two approaching nanoparticles moving through solvent and on a single nanoparticle as it approaches a planar surface. In most circumstances, macroscale fluid theory accurately predicts the drag on these nano-scale particles. Good agreement is seen with Brenner's analytical solutions for wall separations larger than the soft nanoparticle radius. For two approaching coated nanoparticles, the solvent-mediated (velocity-independent) and lubrication (velocity-dependent) forces are purely repulsive and do not exhibit force oscillations that are typical of uncoated rigid spheres.There is increasing interest in using nanoparticles in commercial and industrial applications. However, effective processing of nanoparticles requires that they do not aggregate and often involves solvation in a fluid. Functionalizing the nanoparticles accomplishes both goals. The behavior of functionalized nanoparticles depends strongly on the attached groups. The behavior of bare, nonfunctionalized nanoparticles has been studied via experiments, theory, and simulation, as have the interactions of polymer-grafted surfaces [1,2,3,4]. The hydrodynamic and nanoparticle-nanoparticle interactions involving small functionalized nanoparticles, however, are harder to characterize. Experimentally, it is difficult to manipulate and measure forces on individual nanoparticles smaller than ∼ 100 nm. Theoretical treatments are challenging because the coatings are relatively short, while the particles themselves are outside the large radius of curvature limit. Numerical simulations of discrete solvent effects on nanoparticles have been impractical until now because of the large systems required to avoid significant finite-size effects due to the long-range hydrodynamic interactions. For single-particle diffusion, these corrections have been shown [5,6] to scale as R/L where R is the particle radius and L is the simulation cell length.Recent studies have computed the potential of mean force for bare silica nanoparticles in an aqueous medium, with and without electrolytes present [7]. Forces between bare colloidal nanoparticles have also been studied in LennardJones fluids and in n-decane [8]. Hydrodynamic drag influenced by approach to a plane surface has been studied theoretically [9] and compared to simulations for rigid spheres [10,11]. Alignment effects for amorphous nanoparticles have also been studied [12]. Kim et al.[13] used molecular dynamics (MD) simulations to study the relaxation of a fullerene molecule coated with poly(ethylene oxide) (PEO). Other simulations have either relied u...