Investigation of benzene-hexafluorobenzene dynamics in liquid binary mixtures

The structure and microscopic dynamics of liquid mixtures of benzene and hexafluorobenzene at room temperature and several compositions have been studied by molecular-dynamics simulations. In this implementation we have rescaled the intermolecular H-F cross potential parameters obtained from the Lorentz-Berthelot combining rules, in order to avoid the substantial overestimation of the energy of mixing predicted by the model when the usual rules are employed. We found that a reduction in the strength of cross H-F interactions by 50% relative to the geometric mean is required in order to get a good agreement with experiments. Radial-angular pair-correlation functions between like and unlike species have been computed and analyzed, by comparing them with the correlations in the corresponding neat liquids. We have also studied the microscopic intermolecular momentum transfer, by computing the time correlation function between the initial velocity of a central molecule and later velocities of neighboring molecules. Structural and dynamical information extracted from the mentioned functions seem to be consistent with the picture of relatively long-lived benzene-hexafluorobenzene (Bz-Hf) complexes present in the mixtures, which would be responsible for the considerable perturbation of the structure in the first shell of like species, and would be moving within the liquid in a parallel face-to-face configuration. Using the tools developed originally to estimate hydrogen-bond lifetimes in liquids, we have computed the lifetimes of the Bz-Hf complexes as a function of the mixture composition, by two different methods: the direct time-averaging scheme and from the autocorrelation function of bond occupation numbers. The obtained lifetimes are strongly dependent on the scheme chosen to compute the characteristic times. We have obtained for the Bz-Hf dimer in solution, at room temperature, lifetimes in the range of 30-40 ps from averaging schemes and around 60-120 ps from autocorrelation function methods. In the latter case, the longest times correspond to the equimolar mixture.