Implicit solvent model for effective molecular dynamics simulations of systems composed of colloid nanoparticles and carbon nanotubes

Abstract

In this paper we propose an implicit solvent model which can be used in molecular dynamics simulations of systems comprising colloid nanoparticles and carbon nanotubes. Such systems, due to finite nanometer sizes of both components cannot be accurately approximated by a smaller slab geometry and thus represent a particularly difficult case in terms of computer simulations. In particular, nanoparticle sizes of a few tens of nanometers lead to billions of solvent molecules in a simulation box and require very long cut-off distances which drastically increases computation time. To overcome this difficulty we develop an implicit solvent model based on Hamaker theory of dispersive interactions. The predictions of our model are verified by comparison with the exact model, involving all atoms and full description of pair interactions. The proposed model correctly predicts the work of adhesion and average configuration in colloid – carbon nanotube systems. Moreover, application of the Langevin dynamics reproduces the dynamic behaviour of the exact model either.