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

Abstract

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. File list (2)download file view on ChemRxiv Thermodiffusion-main-final.pdf (1.09 MiB) download file view on ChemRxiv Thermodiffusion-SI-final.pdf (495.88 KiB)

Thermodiffusion of CO2 in Water by Nonequilibrium Molecular Dynamics Simulations

“…The computations are based on the trajectory of the particles, obtained by integrating Newton’s equation of motion, and the use of a force field to represent the interaction potential between particles. Molecular dynamics simulations have been employed to compute thermodiffusion in systems containing hydrocarbons, − Lennard-Jones fluids, ,,− confined fluids, ,, associating liquids, ,− electrolytes in aqueous solution, − and silicate liquids …”

Mass Transport in “Water-in-Polymer Salt” Electrolytes

“…41 Such a force field combination was recently recommended in a benchmark study about the thermodiffusion of ions in aqueous solution. 42 The initial geometry of the polymer (see the SI) was obtained by geometry optimization at the DFT/M06/6-31G(d,p) theory level, and its atomic charges were determined by fitting the DFT electrostatic potential using the Merz−Singh−Kollman scheme 43 combined with the polarizable continuum model (PCM) to account for the solvent polarization effects. Both the geometry and the atomic charges of the polymer were calculated in the presence of the counterions, for the sake of the charge neutralization of the system.…”

Optical Manipulation Heats up: Present and Future of Optothermal Manipulation

“…To address such challenges, molecular dynamics (MD) simulations are commonly used to study the thermal drifts of ions and estimate the Soret coefficients that determine their drift direction and velocity. 64 MD simulations involve modeling the interaction potential between atoms and molecules in the system and can give insights into system characteristics at the subnanometer scale. However, such simulations are limited by computational cost and can only be applied to nanometer-sized domains with high pressure, high temperature gradient, and high concentrations of salt to observe a concentration gradient that helps visualize a thermoelectric potential difference (temperature gradients of the order of 10 9 K/m are required to observe a concentration change Δc of 2−5% within the nanometer scale simulation domains).…”

“…However, in biocompatible solvents such as phosphate-buffered saline containing 90% NaCl, the thermophoretic response of living entities still remains difficult to estimate due to the response of several other ions in different concentrations in the solution. To address such challenges, molecular dynamics (MD) simulations are commonly used to study the thermal drifts of ions and estimate the Soret coefficients that determine their drift direction and velocity . MD simulations involve modeling the interaction potential between atoms and molecules in the system and can give insights into system characteristics at the subnanometer scale.…”