These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
2 files

A Molecular Dynamics Approach to Calculate the Thermodiffusion (Soret and Seebeck) Coefficients of Salts in Aqueous Solutions

submitted on 21.04.2021, 16:50 and posted on 26.04.2021, 04:21 by Leandro Rezende Franco, André Luiz Sehnem, Antônio Martins Figueiredo Neto, Kaline Coutinho

An approach to investigate the physical parameters related to the ions thermodiffusion in aqueous solution is proposed herein by calculating the equilibrium hydration free energy and the self-diffusion coefficient as a function of temperature, ranging from 293 to 353 K, using molecular dynamics simulations of infinitely diluted ions in aqueous solutions. Several ion force field parameters are used in the simulations and new parameters are proposed for some ions to better describe their hydration free energy. Such a theoretical framework enables the calculation of some single-ion properties, such as heat of transport, Soret coefficient and mass current density, as well as properties of salts, such as effective mass and thermal diffusion, Soret and Seebeck coefficients. These calculated properties are compared with experimental data available from optical measurements and showed good agreement revealing an excellent theoretical predictability of salt thermodiffusion properties. Differences in single-ion Soret and self-diffusion coefficients of anions and cations give rise to a thermoelectric field, which affects the system response that is quantified by the Seebeck coefficient. The fast and slow Seebeck coefficients are calculated and discussed, resulting in values with mV/K order-of-magnitude, as observed in experiments involving several salts, such as K+Cl−, Na+Cl−, H+Cl−, Na+OH−, TMA+OH− and TBA+OH−. The present approach can be adopted for any ion or charged particle dispersed in water with the aim of predicting the thermoelectric field induced through the fluid. It has potential applications in designing electrolytes for ionic thermoelectric devices in order to harvest energy and thermoelectricity in biological nanofluids.


National Institute of Science and Technology of Complex Fluids (INCT-FCx) for the CNPq grant 141260/2017-3 and FAPESP grant 2014/50983-3

CAPES concerning the BioMol project 23038.004630/2014-35

FAPESP grants 2017/11631-2, 2019/10433- 33 8 and 2016/24531-3

High Performance Computing of USP (HPC-USP) for computational resources


Email Address of Submitting Author


Universidade de Sao Paulo, Instituto de Fisica



ORCID For Submitting Author


Declaration of Conflict of Interest

no conflict of interest

Version Notes

This version was accepted for publication in Journal of Chemical Theory and Computation (20-Apr-2021).