Modeling the structure and relaxation in glycerol-silica nanocomposites

The relationship between the dynamics and structure of amorphous thin films and nanocomposites near their glass transition is an important problem in soft-matter physics. Here, we develop a simple theoretical approach to describe the density profile and the α-relaxation time of a glycerol-silica nanocomposite (S. Cheng et al., J. Chem. Phys., 2015, 143, 194704). We begin by applying the Derjaguin approximation, where we replace the curved surface of the particle with the planar one; thus, modeling the nanocomposite is reduced to that of a confined thin film. Subsequently, by employing the molecular dynamics (MD) simulation data of Cheng et al., we approximate the density profile of a supported liquid thin film as a stationary solution of a fourth-order partial differential equation (PDE). We then construct an appropriate density functional, from which the density profile emerges through the minimization of free energy. Our final assumption is that of a consistent, temperature-independent scaled density profile, ensuring that the free volume throughout the entire nanocomposite increases with temperature in a smooth, monotonic fashion. Considering the established relationship between glycerol relaxation time and temperature, we can employ Doolittle-type analysis ("naïve" free-volume model), to calculate the relaxation time based on temperature and film thickness. We then convert the film thickness into the interparticle distance and subsequently the filler volume fraction for the nanocomposites and compare our model predictions with experimental data, resulting in a good agreement. The proposed approach can be easily extended to other nanocomposite and film systems.