LAMMPS is a very customizable molecular dynamics simulation software, which can be used to simulate a large diversity of systems. We introduce a new package for simulation of polarizable systems with LAMMPS using thermalized Drude oscillators. The implemented functionalities are described and are illustrated by examples. The implementation was validated by comparing simulation results with published data and using a reference software. Computational performance is also analyzed.

Over the past twenty years experiments performed on thin polymer films deposited on substrates have shown that the glass transition temperature T(g) can either decrease or increase depending on the strength of the interactions. Over the same period, experiments have also demonstrated that the dynamics in liquids close to the glass transition temperature is strongly heterogeneous, on the scale of a few nanometers. A model for the dynamics of non-polar polymers, based on percolation of slow subunits, has been proposed and developed over the past ten years. It proposes a unified mechanism regarding these two features. By extending this model, we have developed a 3D model, solved by numerical simulations, in order to describe and calculate the mechanical properties of polymers close to the glass transition in the linear regime of deformation, with a spatial resolution corresponding to the subunit size. We focus on the case of polymers confined between two substrates with non-negligible interactions between the polymer and the substrates, a situation which may be compared to filled elastomers. We calculate the evolution of the elastic modulus as a function of temperature, for different film thicknesses and polymer-substrate interactions. In particular, this allows to calculate the corresponding increase of glass transition temperature, up to 20 K in the considered situations. Moreover, between the bulk T(g) and T(g) + 50 K the modulus of the confined layers is found to decrease very slowly in some cases, with moduli more than ten times larger than that of the pure matrix at temperatures up to T(g) + 50 K. This is consistent with what is observed in reinforced elastomers. This slow decrease of the modulus is accompanied by huge fluctuations of the stress at the scale of a few tens of nanometers that may even be negative as compared to the solicitation, in a way that may be analogous to mechanical heterogeneities observed recently in molecular dynamics simulations. As a consequence, confinement may result not only in an increase of the glass transition temperature, but in a huge broadening of the glass transition.