The Soret effect and isotopic fractionation in high-temperature silicate melts

The role of thermodiffusion in transpiration

Plant leaf temperatures can differ from ambient air temperatures. A temperature gradient in a gas mixture gives rise to a phenomenon known as thermodiffusion, which operates in addition to ordinary diffusion. Whilst transpiration is generally understood to be driven solely by the ordinary diffusion of water vapour along a concentration gradient, we consider the implications of thermodiffusion for transpiration. We develop a new modelling framework that introduces the effects of thermodiffusion on the transpiration rate, E. By applying this framework, we quantify the proportion of E attributable to thermodiffusion for a set of physiological and environmental conditions, varied over a wide range. Thermodiffusion is found to be most significant (in some cases > 30% of E) when a leaf-to-air temperature difference coincides with a relatively small water vapour concentration difference across the boundary layer; a boundary layer conductance that is large as compared to the stomatal conductance; or a relatively low transpiration rate. Thermodiffusion also alters the conditions required for the onset of reverse transpiration, and the rate at which this water vapour uptake occurs.

Composition-dependent sign inversion of the Soret coefficient of SiO2 in binary borosilicate melts

Using a laser-induced local-heating experiment combined with temperature analysis, we observed the composition-dependent sign inversion of the Soret coefficient of SiO2 in binary silicate melts, which was successfully explained by a modified Kempers model used for describing the Soret effect in oxide melts. In particular, the diffusion of SiO2 to the cold side under a temperature gradient, which is an anomaly in silicate melts, was observed in the SiO2-poor compositions. The theoretical model indicates that the thermodynamic mixing properties of oxides, partial molar enthalpy of mixing, and partial molar volume are the dominant factors for determining the migration direction of the SiO2 component under a temperature gradient.

Molecular Dynamics Simulation of the Soret Effect on Two Binary Liquid Solutions with Equimolar n-Alkane Mixtures

Molecular dynamics is employed to simulate the Soret effect on two binary liquid solutions with equimolar mixtures: normal pentane (n-pentane, nC-5) and normal heptane (n-heptane, nC-7) molecules plus normal decane (n-decane, nC-10) and normal pentane molecules. Moreover, two coarse-grained force field (the CG-FF) potentials, which may depict inter-/intramolecular interactions fairly well among n-alkane molecules, are developed to fulfill such investigations. In addition, thermal diffusion for the mass fraction of each of these n-alkane molecules is simulated under an effect of a weak thermal gradient (temperature difference) exerting on solution systems from their hot to cold boundary sides. Finally, quantities of the Soret coefficient (SC) for two binary solutions are calculated by means of the developed CG-FF potentials, so as to improve the calculation rationality. As a result, first, it is found that molecules with light molar masses will migrate toward the hot boundary side, while those with heavy molar masses will migrate toward the cold boundary one ; second, the SC quantities indicate that they match relevant experimental determinations fairly well, i.e., trends of these SC quantities show inverse proportionality to the thermal gradient on the systems.

Determination of thermodynamic and microscopic origins of the Soret effect in sodium silicate melts: Prediction of sign change of the Soret coefficient

The Soret effect in silicate melts has attracted attention in earth and material sciences, particularly in glass science and engineering, because a compositional change caused by the Soret effect modifies the material properties of silicate melts. We investigated the Soret effect in an Na₂O-SiO₂ system, which is the most common representative of silicate melts. Our theoretical approach based on the modified Kempers model and non-equilibrium molecular dynamics simulation was validated for 30Na₂O-70SiO₂(mol. %). The sign and order of the absolute values of the calculated Soret coefficients were consistent with the experimental values. The positive Soret coefficient of SiO₂ in the SiO₂-poor composition range was accurately predicted. Previous experimental studies have focused on SiO₂-rich compositions, and only the negative sign, indicating SiO₂ migration to the hot side, has been observed. In the SiO₂-poor composition range, the Q₀ structure was dominant and had four Si-O-Na bonds around an SiO₄ unit. The Si-O-Na bond had high enthalpic stability and contributed to the large negative enthalpy of SiO₂ mixing. According to our model, components with a large negative partial molar enthalpy of mixing will concentrate in the cold region. The microscopic and thermodynamic origins of the sign change in the Soret effect were determined.

Role of partial molar enthalpy of oxides on Soret effect in high-temperature CaO-SiO2 melts

The Soret effect or thermodiffusion is the temperature-gradient driven diffusion in a multicomponent system. Two important conclusions have been obtained for the Soret effect in multicomponent silicate melts: first, the SiO₂ component concentrates in the hot region; and second, heavier isotopes concentrate in the cold region more than lighter isotopes. For the second point, the isotope fractionation can be explained by the classical mechanical collisions between pairs of particles. However, as for the first point, no physical model has been reported to answer why the SiO₂ component concentrates in the hot region. We try to address this issue by simulating the composition dependence of the Soret effect in CaO–SiO₂ melts with nonequilibrium molecular dynamics and determining through a comparison of the results with those calculated from the Kempers model that partial molar enthalpy is one of the dominant factors in this phenomenon.

Microscopic mechanism of thermomolecular orientation and polarization

Recent molecular dynamics simulations show that thermal gradients can induce electric fields in water that are comparable in magnitude to electric fields seen in ionic thin films and biomembranes. This surprising non-equilibrium phenomenon of thermomolecular orientation is also observed more generally in simulations of polar and non-polar size-asymmetric dumbbell fluids. However, a microscopic theory linking thermomolecular orientation and polarization to molecular properties is yet unknown. Here, we formulate an analytically solvable microscopic model of size-asymmetric dumbbell molecules in a temperature gradient using a mean-field, local equilibrium approach. Our theory reveals the relationship between the extent of thermomolecular orientation and polarization, and molecular volume, size anisotropy and dipole moment. Predictions of the theory agree quantitatively with molecular dynamics simulations. Crucially, our framework shows how thermomolecular orientation can be controlled and maximized by tuning microscopic molecular properties.

Why charged molecules move across a temperature gradient: the role of electric fields

Methods to move solvated molecules are rare. Apart from electric fields, only thermal gradients are effective enough to move molecules inside a fluid. This effect is termed thermophoresis, and the underlying mechanisms are still poorly understood. Nevertheless, it is successfully used to quantify biomolecule binding in complex liquids. Here we show experiments that reveal that thermophoresis in water is dominated by two electric fields, both established by the salt ions of the solution. A local field around the molecule drives molecules along an energy gradient, whereas a global field moves the molecules by a combined thermoelectrophoresis mechanism known as the Seebeck effect. Both mechanisms combined predict the thermophoresis of DNA and RNA polymers for a wide range of experimental parameters. For example, we correctly predict a complex, nonlinear size transition, a salt-species-dependent offset, a maximum of thermophoresis over temperature, and the dependence of thermophoresis on the molecule concentration.

Unidirectional thermal diffusion in bimetallic Cu@Au nanoparticles

Understanding the atomic diffusions at the nanoscale is important for controlling the synthesis and utilization of nanomaterials. Here, using in situ X-ray absorption spectroscopy coupled with theoretical calculations, we demonstrate a so far unexplored unidirectional diffusion from the Au shell to the Cu core in thermally alloying Cu@Au core@shell architecture of ca. 7.1 nm. The initial diffusion step at 423 K is found to be characterized by the formation of a diffusion layer composed of a Au-dilute substitutional CuAu-like intermetallic compound with short Cu-Au bond length (2.61 Å). The diffusion further happens by the migration of the Au atoms with large disorder into the interior Cu matrix at higher temperatures (453 and 553 K). These results suggest that the structural preference of a CuAu-like compound, along with the nanosized effect, plays a critical role in determining the atomic diffusion dynamics.

Thermodiffusion in binary liquids: the role of irreversibility

We study thermal diffusion in binary mixtures in the framework of non-equilibrium thermodynamics. Our formal result displays the role of partial enthalpies hi and Onsager's generalized mobilities Ai. The ratio A1/A2 provides a measure for the irreversible character of thermal diffusion. Comparison with experimental data on benzene, cyclohexane, toluene and n-alkanes shows that irreversibility is essential for thermal diffusion, and in particular for the isotope effect.

The Physical Chemistry of Mass-Independent Isotope Effects and Their Observation in Nature

Historically, the physical chemistry of isotope effects and precise measurements in samples from nature have provided information on processes that could not have been obtained otherwise. With the discovery of a mass-independent isotopic fractionation during the formation of ozone, a new physical chemical basis for isotope effects required development. Combined theoretical and experimental developments have broadened this understanding and extended the range of chemical systems where these unique effects occur. Simultaneously, the application of mass-independent isotopic measurements to an extensive range of both terrestrial and extraterrestrial systems has furthered the understanding of events such as solar system origin and evolution and planetary atmospheric chemistry, present and past.

Isotope fractionation by thermal diffusion in silicate melts

Isotopes fractionate in thermal gradients, but there is little quantitative understanding of this effect in complex fluids. Here we present results of experiments and molecular dynamics simulations on silicate melts. We show that isotope fractionation arises from classical mechanical effects, and that a scaling relation based on Chapman-Enskog theory predicts the behavior seen in complex fluids without arbitrary fitting parameters. The scaling analysis reveals that network forming elements (Si and O) fractionate significantly less than network modifiers (e.g., Mg, Ca, Fe, Sr, Hf, and U).