Heat and matter transport in binary liquid mixtures

Charge Transport in Water-NaCl Electrolytes with Molecular Dynamics Simulations

A systematic description of microscopic mechanisms is necessary to understand mass transport in solid and liquid electrolytes. From Molecular Dynamics (MD) simulations, transport properties can be computed and provide a detailed view of the molecular and ionic motions. In this work, ionic conductivity and transport numbers in electrolyte systems are computed from equilibrium and nonequilibrium MD simulations. Results from the two methods are compared with experimental results, and we discuss the significance of the frame of reference when determining and comparing transport numbers. Two ways of computing ionic conductivity from equilibrium simulations are presented: the Nernst-Einstein approximation or the Onsager coefficients. The Onsager coefficients take ionic correlations into account and are found to be more suitable for concentrated electrolytes. Main features and differences between equilibrium and nonequilibrium simulations are discussed, and some potential anomalies and critical pitfalls of using nonequilibrium molecular dynamics to determine transport properties are highlighted.

Non-equilibrium molecular dynamics simulations of the transient Ludwig-Soret effect in a binary Lennard-Jones/spline mixture

A binary isotope mixture of Lennard-Jones/spline particles at equilibrium was perturbed by a sudden change in the system's boundary temperatures. The system's response was determined by non-equilibrium molecular dynamics (NEMD). Three transient processes were studied: 1) The propagation of a pressure (shock) wave, 2) heat diffusivity and conduction, and 3) thermal diffusion (the Ludwig-Soret effect). These three processes occur at different time scales, which makes it possible to separate them in one single NEMD run. The system was studied in liquid, supercritical, and dense gas states with various forms and strengths of the thermal perturbation. The results show that heat was initially transported by two separate mechanisms: 1) heat diffusion as described by the transient heat equation and 2) as a consequence of a pressure wave. The pressure wave travelled faster than the speed of sound, generating a shock wave in the system. Local equilibrium was found in the transient phase, even with very strong perturbations and in the shock front. Although the mass separation due to the Ludwig-Soret effect developed much slower than the pressure and temperature fields in the system at large, it was found that the Soret coefficient could be accurately determined from the initial phase of the transient and close to the heat source. This opens the possibility of a new way to analyse results from transient experiments and thereby minimize effects of gravity and convection due to buoyancy.

The rich phase behavior of the thermopolarization of water: from a reversal in the polarization, to enhancement near criticality conditions

Non-equilibrium molecular dynamics simulations show that the polarization of water induced by thermal gradients depends strongly on the thermodynamic conditions, with a large enhancement near the critical point.

Modeling thermophoretic effects in solid-state nanopores

Local modulation of temperature has emerged as a new mechanism for regulation of molecular transport through nanopores. Predicting the effect of such modulations on nanopore transport requires simulation protocols capable of reproducing non-uniform temperature gradients observed in experiment. Conventional molecular dynamics (MD) method typically employs a single thermostat for maintaining a uniform distribution of temperature in the entire simulation domain, and, therefore, can not model local temperature variations. In this article, we describe a set of simulation protocols that enable modeling of nanopore systems featuring non-uniform distributions of temperature. First, we describe a method to impose a temperature gradient in all-atom MD simulations based on a boundary-driven non-equilibrium MD protocol. Then, we use this method to study the effect of temperature gradient on the distribution of ions in bulk solution (the thermophoretic effect). We show that DNA nucleotides exhibit differential response to the same temperature gradient. Next, we describe a method to directly compute the effective force of a thermal gradient on a prototypical biomolecule-a fragment of double-stranded DNA. Following that, we demonstrate an all-atom MD protocol for modeling thermophoretic effects in solid-state nanopores. We show that local heating of a nanopore volume can be used to regulate the nanopore ionic current. Finally, we show how continuum calculations can be coupled to a coarse-grained model of DNA to study the effect of local temperature modulation on electrophoretic motion of DNA through plasmonic nanopores. The computational methods described in this article are expected to find applications in rational design of temperature-responsive nanopore systems.

Do thermal diffusion and Dufour coefficients satisfy Onsager's reciprocity relation?

It is commonly admitted that in liquids the thermal diffusion and Dufour coefficients D(T) and D(F) satisfy Onsager's reciprocity. From their relation to the cross-coefficients of the phenomenological equations, we are led to the conclusion that this is not the case in general. As illustrative and physically relevant examples, we discuss micellar solutions and colloidal suspensions, where D(T) arises from chemical reactions or viscous effects but is not related to the Dufour coefficient D(F). The situation is less clear for binary molecular mixtures; available experimental and simulation data do not settle the question whether D(T) and DF are reciprocal coefficients.

Computation of thermodynamic and transport properties to predict thermophoretic effects in an argon-krypton mixture

Thermophoresis is the movement of molecules caused by a temperature gradient. Here we report the results of a study of thermophoresis using non-equilibrium molecular dynamics simulations of a confined argon-krypton fluid subject to two different temperatures at thermostated walls. The resulting temperature profile between the walls is used along with the Soret coefficient to predict the concentration profile that develops across the channel. We obtain the Soret coefficient by calculating the mutual diffusion and thermal diffusion coefficients. We report an appropriate method for calculating the transport coefficients for binary systems, using the Green-Kubo integrals and radial distribution functions obtained from equilibrium molecular dynamics simulations of the bulk fluid. Our method has the unique advantage of separating the mutual diffusion and thermal diffusion coefficients, and calculating the sign and magnitude of their individual contributions to thermophoresis in binary mixtures.

Influence of the interaction potential and of the temperature on the thermodiffusion (Soret) coefficient in a model system

In this paper the thermodiffusive behavior of an equimolar binary mixture subject to repulsive potentials of the form (sigma/r)(n) is investigated by using nonequilibrium molecular dynamics (NEMD) and the thermodiffusion (Soret) coefficient, S(T), is computed in a wide range of temperatures. With the aim to contribute to the study of the dependence of the Soret coefficient on the interaction potential, the exponent n of the potential is varied from 1 to 12, that is from a pseudocoulombian to a pseudohard-sphere interaction. The steady state equation is integrated for the composition function under reasonable assumptions and it is shown that in some cases the request for it to be linear cannot be satisfied. For this reason nonlinear functions are used to fit the NEMD composition data. The simulations indicate a negligible dependence of S(T) on the composition (in the composition range here considered) while the dependence on the temperature is more marked. The computed values of S(T) as a function of the temperature are fitted with analytical functions. It is found that with n> or =3 (medium and short range interaction) the model system behaves like a dilute gas mixture with the Soret coefficient varying with the temperature almost like 1/T. In the case of n=1 (long range interaction), S(T) has a more complex dependence on T: in particular it shows a change of sign. The analytical fitting functions, S(T)(T), are used in the integrated steady state equation thus obtaining the steady state composition profile and its comparison with the NEMD results indicates the grounding of the approach here proposed.

Equilibrium and nonequilibrium molecular dynamics simulations of the thermal conductivity of molten alkali halides

The thermal conductivity of molten NaCl and KCl was calculated through the Evans-Gillan nonequilibrium molecular dynamics (NEMD) algorithm and Green-Kubo equilibrium molecular dynamics (EMD) simulations. The EMD simulations were performed for a "binary" ionic mixture and the NEMD simulations assumed a pure system for reasons discussed in this work. The cross thermoelectric coefficient obtained from Green-Kubo EMD simulations is discussed in terms of the homogeneous thermoelectric power or Seebeck coefficient of these materials. The thermal conductivity obtained from NEMD simulations is found to be in very good agreement with that obtained through Green-Kubo EMD simulations for a binary ionic mixture. This result points to a possible cancellation between the neglected "partial enthalpy" contribution to the heat flux associated with the interdiffusion of one species through the other and that part of the thermal conductivity related to the coupled fluxes of charge and heat in "binary" ionic mixtures.

Microscopic interpretation of a pure chemical contribution to the Soret effect

In a recent Letter by Köhler [Phys. Rev. Lett. 87, 055901 (2001)], it has been shown that the Soret effect or thermal diffusion can be split into three different contributions: mass, moment of inertia, and a so-called chemical effect, but only the chemical effect gives rise to a composition dependent contribution. As it is experimentally difficult to deal with the chemical contribution without changing the two others, it has not been studied accurately yet. Our Letter presents both equilibrium and nonequilibrium Molecular Dynamics in simple Lennard-Jones mixtures. By thoroughly changing the strength of direct and cross interaction energies between particles, we show that the composition dependence and the change of sign of the Soret coefficient is driven only by the nature of interactions between unlike particles and propose a microscopic interpretation of the Soret effect.

Computing the Soret coefficient in aqueous mixtures using boundary driven nonequilibrium molecular dynamics

We have computed the Soret coefficient in aqueous mixtures using a boundary driven nonequilibrium molecular dynamics algorithm and standard molecular force fields. The choice of this specific approach is justified by the nature of the mixtures studied here. Four aqueous solutions, including methanol, ethanol, acetone, and dimethyl-sulfoxide (DMSO) have been studied at ambient conditions for different compositions. The experimental behavior of water-alcohol mixtures was reproduced, including the change of sign of the Soret coefficient with composition, in excellent agreement with existing experimental data. The methodology has been applied to obtain pure predictions for water-acetone and water-DMSO where no experimental data are accessible. A change of sign is also observed in the same range of composition as in water-alcohol mixtures. It is suggested that the nature and strength of the molecular interactions, rather than the mass or shape ratio of the components, dominates the behavior of the Soret coefficient versus composition for the aqueous associating mixtures studied here.

Soret coefficient for liquid argon-krypton mixtures via equilibrium and nonequilibrium molecular dynamics: a comparison with experiments

The Soret and the other transport coefficients, characterizing the heat and mass transport in binary mixtures, have been obtained by equilibrium and nonequilibrium molecular dynamics (EMD and NEMD, respectively). Two state points of the argon-krypton mixture are considered, for which experimental values of the Soret coefficient are available. To attempt a comparison between simulations and experiments the common enthalpy-diffusion-free expression for the heat flux has been chosen. The comparison of the simulations with the experiments shows a remarkable agreement, for all the several utilized EMD and NEMD techniques (dynamical and stationary). The techniques, used over 0.3 micros of total simulation time span, are slow convergent but have comparable performances.

Coupled normal heat and matter transport in a simple model system

We introduce the first simple mechanical system that shows fully realistic transport behavior while still being exactly solvable at the level of equilibrium statistical mechanics. The system is a Lorentz gas with fixed freely rotating circular scatterers which scatter point particles via perfectly rough collisions. Upon imposing either a temperature gradient and/or a chemical potential gradient, a stationary state is attained for which local thermal equilibrium holds. Transport in this system is normal in the sense that the transport coefficients which characterize the flow of heat and matter are finite in the thermodynamic limit. Moreover, the two flows are nontrivially coupled, satisfying Onsager's reciprocity relations.