Predicting and rationalizing the Soret coefficient of binary Lennard-Jones mixtures in the liquid state

05 July 2022, Version 2
This content is a preprint and has not undergone peer review at the time of posting.


The thermodiffusion behavior of binary Lennard-Jones mixtures in the liquid state was investigated by combining the individual strengths of non-equilibrium molecular dynamics (NEMD) and equilibrium molecular dynamics (EMD) simulations. On the one hand, boundary-driven NEMD simulations are useful to quickly predict Soret coefficients because they are easy to set up and straightforward to analyze. However, careful interpolation is required because the mean temperature in the measurement region does not exactly reach the target temperature. On the other hand, EMD simulations attain the target temperature precisely and yield a multitude of properties that clarify the microscopic origins of Soret coefficient trends. An analysis of the Soret coefficient suggests a straightforward dependence on thermodynamic properties, whereas its dependence on dynamic properties is far more complex. Furthermore, a limit of applicability of a popular theoretical model, which mainly relies on thermodynamic data, was identified by virtue of an uncertainty analysis in conjunction with efficient empirical Soret coefficient predictions, which rely on model parameters instead of simulation output. Finally, the present study underscores that a combination of predictive models and EMD and NEMD simulations is a powerful approach to shed light onto the thermodiffusion behavior of liquid mixtures.


Soret coefficient
molecular simulation
non-equilibrium molecular dynamics
Green-Kubo formalism
Lennard-Jones potential
thermodynamic models

Supplementary materials

Supporting Information (SI) accompanying the research paper
Details of the methodology are presented in this SI. Furthermore, additional validation results and tests are included. Numerical results of EMD simulations are found here as well as an uncertainty analysis of a predictive model. Finally, the list of symbols is given.


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.