Atrazine is one of the most widely used herbicide molecules in the triazine family. Despite its interdiction in the European Union in 2004, atrazine and its main degradation products remain among the most frequently found molecules in freshwater reservoirs in many European Union countries. Our study aims in obtaining insight into the desorption process of atrazine from the main soil absorbent material: clay. Constrained Molecular Dynamics simulations within the Density Functional Theory framework allow us to obtain a free energy desorption profile of atrazine from a Ca2+-montmorillonite surface. The results are interpreted in terms of atrazine inclination to the clay surface and moreover, in terms of hydration states of the cations present in the clay interlayer as well as the hydration state of the atrazine. The desorption mechanism is driven by atrazine alkyl groups and their sizes because of dispersion stabilizing effects. The highest barrier corresponds to the loss of the isopropyl interaction with the surface.