Importance of a Fully Anharmonic Treatment of Equilibrium Isotope Fractionation Properties of Dissolved Ionic Species As Evidenced by Li+(aq)

Romain Dupuis, Magali Benoit, Mark Tuckerman, Merlin Méheut

Research output: Contribution to journalArticle

Abstract

ConspectusEquilibrium fractionation of stable isotopes is critically important in fields ranging from chemistry, including medicinal chemistry, electrochemistry, geochemistry, and nuclear chemistry, to environmental science. The dearth of reliable estimates of equilibrium fractionation factors, from experiment or from natural observations, has created a need for accurate computational approaches. Because isotope fractionation is a purely quantum mechanical phenomenon, exact calculation of fractionation factors is nontrivial. Consequently, a severe approximation is often made, in which it is assumed that the system can be decomposed into a set of independent harmonic oscillators. Reliance on this often crude approximation is one of the primary reasons that theoretical prediction of isotope fractionation has lagged behind experiment. A class of problems for which one might expect the harmonic approximation to perform most poorly is the isotopic fractionation between solid and solution phases.In order to illustrate the errors associated with the harmonic approximation, we have considered the fractionation of Li isotopes between aqueous solution and phyllosilicate minerals, where we find that the harmonic approximation overestimates isotope fractionation factors by as much as 30% at 25 °C. Lithium is a particularly interesting species to examine, as natural lithium isotope signatures provide information about hydrothermal processes, carbon cycle, and regulation of the Earth's climate by continental alteration. Further, separation of lithium isotopes is of growing interest in the nuclear industry due to a need for pure 6Li and 7Li isotopes. Moving beyond the harmonic approximation entails performing exact quantum calculations, which can be achieved using the Feynman path integral formulation of quantum statistical mechanics. In the path integral approach, a system of quantum particles is represented as a set of classical-like ring-polymer chains, whose interparticle interactions are determined by the rules of quantum mechanics. Because a classical isomorphism exists between the true quantum system and the system of ring-polymers, classical-like methods can be applied. Recent developments of efficient path integral approaches for the exact calculation of isotope fractionation now allow the case of the aforementioned dissolved Li fractionation properties to be studied in detail. Applying this technique, we find that the calculations yield results that are in good agreement with both experimental data and natural observations. Importantly, path integral methods, being fully atomistic, allow us to identify the origins of anharmonic effects and to make reliable predictions at temperatures that are experimentally inaccessible yet are, nevertheless, relevant for natural phenomena.

Original languageEnglish (US)
Pages (from-to)1597-1605
Number of pages9
JournalAccounts of Chemical Research
Volume50
Issue number7
DOIs
StatePublished - Jul 18 2017

Fingerprint

Fractionation
Isotopes
Lithium
Polymers
Nuclear industry
Statistical mechanics
Geochemistry
Quantum theory
Electrochemistry
Minerals
Carbon
Earth (planet)
Experiments

ASJC Scopus subject areas

  • Chemistry(all)

Cite this

Importance of a Fully Anharmonic Treatment of Equilibrium Isotope Fractionation Properties of Dissolved Ionic Species As Evidenced by Li+(aq). / Dupuis, Romain; Benoit, Magali; Tuckerman, Mark; Méheut, Merlin.

In: Accounts of Chemical Research, Vol. 50, No. 7, 18.07.2017, p. 1597-1605.

Research output: Contribution to journalArticle

@article{568d49266e1446be9d2e4382b8c9be5e,
title = "Importance of a Fully Anharmonic Treatment of Equilibrium Isotope Fractionation Properties of Dissolved Ionic Species As Evidenced by Li+(aq)",
abstract = "ConspectusEquilibrium fractionation of stable isotopes is critically important in fields ranging from chemistry, including medicinal chemistry, electrochemistry, geochemistry, and nuclear chemistry, to environmental science. The dearth of reliable estimates of equilibrium fractionation factors, from experiment or from natural observations, has created a need for accurate computational approaches. Because isotope fractionation is a purely quantum mechanical phenomenon, exact calculation of fractionation factors is nontrivial. Consequently, a severe approximation is often made, in which it is assumed that the system can be decomposed into a set of independent harmonic oscillators. Reliance on this often crude approximation is one of the primary reasons that theoretical prediction of isotope fractionation has lagged behind experiment. A class of problems for which one might expect the harmonic approximation to perform most poorly is the isotopic fractionation between solid and solution phases.In order to illustrate the errors associated with the harmonic approximation, we have considered the fractionation of Li isotopes between aqueous solution and phyllosilicate minerals, where we find that the harmonic approximation overestimates isotope fractionation factors by as much as 30{\%} at 25 °C. Lithium is a particularly interesting species to examine, as natural lithium isotope signatures provide information about hydrothermal processes, carbon cycle, and regulation of the Earth's climate by continental alteration. Further, separation of lithium isotopes is of growing interest in the nuclear industry due to a need for pure 6Li and 7Li isotopes. Moving beyond the harmonic approximation entails performing exact quantum calculations, which can be achieved using the Feynman path integral formulation of quantum statistical mechanics. In the path integral approach, a system of quantum particles is represented as a set of classical-like ring-polymer chains, whose interparticle interactions are determined by the rules of quantum mechanics. Because a classical isomorphism exists between the true quantum system and the system of ring-polymers, classical-like methods can be applied. Recent developments of efficient path integral approaches for the exact calculation of isotope fractionation now allow the case of the aforementioned dissolved Li fractionation properties to be studied in detail. Applying this technique, we find that the calculations yield results that are in good agreement with both experimental data and natural observations. Importantly, path integral methods, being fully atomistic, allow us to identify the origins of anharmonic effects and to make reliable predictions at temperatures that are experimentally inaccessible yet are, nevertheless, relevant for natural phenomena.",
author = "Romain Dupuis and Magali Benoit and Mark Tuckerman and Merlin M{\'e}heut",
year = "2017",
month = "7",
day = "18",
doi = "10.1021/acs.accounts.6b00607",
language = "English (US)",
volume = "50",
pages = "1597--1605",
journal = "Accounts of Chemical Research",
issn = "0001-4842",
publisher = "American Chemical Society",
number = "7",

}

TY - JOUR

T1 - Importance of a Fully Anharmonic Treatment of Equilibrium Isotope Fractionation Properties of Dissolved Ionic Species As Evidenced by Li+(aq)

AU - Dupuis, Romain

AU - Benoit, Magali

AU - Tuckerman, Mark

AU - Méheut, Merlin

PY - 2017/7/18

Y1 - 2017/7/18

N2 - ConspectusEquilibrium fractionation of stable isotopes is critically important in fields ranging from chemistry, including medicinal chemistry, electrochemistry, geochemistry, and nuclear chemistry, to environmental science. The dearth of reliable estimates of equilibrium fractionation factors, from experiment or from natural observations, has created a need for accurate computational approaches. Because isotope fractionation is a purely quantum mechanical phenomenon, exact calculation of fractionation factors is nontrivial. Consequently, a severe approximation is often made, in which it is assumed that the system can be decomposed into a set of independent harmonic oscillators. Reliance on this often crude approximation is one of the primary reasons that theoretical prediction of isotope fractionation has lagged behind experiment. A class of problems for which one might expect the harmonic approximation to perform most poorly is the isotopic fractionation between solid and solution phases.In order to illustrate the errors associated with the harmonic approximation, we have considered the fractionation of Li isotopes between aqueous solution and phyllosilicate minerals, where we find that the harmonic approximation overestimates isotope fractionation factors by as much as 30% at 25 °C. Lithium is a particularly interesting species to examine, as natural lithium isotope signatures provide information about hydrothermal processes, carbon cycle, and regulation of the Earth's climate by continental alteration. Further, separation of lithium isotopes is of growing interest in the nuclear industry due to a need for pure 6Li and 7Li isotopes. Moving beyond the harmonic approximation entails performing exact quantum calculations, which can be achieved using the Feynman path integral formulation of quantum statistical mechanics. In the path integral approach, a system of quantum particles is represented as a set of classical-like ring-polymer chains, whose interparticle interactions are determined by the rules of quantum mechanics. Because a classical isomorphism exists between the true quantum system and the system of ring-polymers, classical-like methods can be applied. Recent developments of efficient path integral approaches for the exact calculation of isotope fractionation now allow the case of the aforementioned dissolved Li fractionation properties to be studied in detail. Applying this technique, we find that the calculations yield results that are in good agreement with both experimental data and natural observations. Importantly, path integral methods, being fully atomistic, allow us to identify the origins of anharmonic effects and to make reliable predictions at temperatures that are experimentally inaccessible yet are, nevertheless, relevant for natural phenomena.

AB - ConspectusEquilibrium fractionation of stable isotopes is critically important in fields ranging from chemistry, including medicinal chemistry, electrochemistry, geochemistry, and nuclear chemistry, to environmental science. The dearth of reliable estimates of equilibrium fractionation factors, from experiment or from natural observations, has created a need for accurate computational approaches. Because isotope fractionation is a purely quantum mechanical phenomenon, exact calculation of fractionation factors is nontrivial. Consequently, a severe approximation is often made, in which it is assumed that the system can be decomposed into a set of independent harmonic oscillators. Reliance on this often crude approximation is one of the primary reasons that theoretical prediction of isotope fractionation has lagged behind experiment. A class of problems for which one might expect the harmonic approximation to perform most poorly is the isotopic fractionation between solid and solution phases.In order to illustrate the errors associated with the harmonic approximation, we have considered the fractionation of Li isotopes between aqueous solution and phyllosilicate minerals, where we find that the harmonic approximation overestimates isotope fractionation factors by as much as 30% at 25 °C. Lithium is a particularly interesting species to examine, as natural lithium isotope signatures provide information about hydrothermal processes, carbon cycle, and regulation of the Earth's climate by continental alteration. Further, separation of lithium isotopes is of growing interest in the nuclear industry due to a need for pure 6Li and 7Li isotopes. Moving beyond the harmonic approximation entails performing exact quantum calculations, which can be achieved using the Feynman path integral formulation of quantum statistical mechanics. In the path integral approach, a system of quantum particles is represented as a set of classical-like ring-polymer chains, whose interparticle interactions are determined by the rules of quantum mechanics. Because a classical isomorphism exists between the true quantum system and the system of ring-polymers, classical-like methods can be applied. Recent developments of efficient path integral approaches for the exact calculation of isotope fractionation now allow the case of the aforementioned dissolved Li fractionation properties to be studied in detail. Applying this technique, we find that the calculations yield results that are in good agreement with both experimental data and natural observations. Importantly, path integral methods, being fully atomistic, allow us to identify the origins of anharmonic effects and to make reliable predictions at temperatures that are experimentally inaccessible yet are, nevertheless, relevant for natural phenomena.

UR - http://www.scopus.com/inward/record.url?scp=85025177567&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85025177567&partnerID=8YFLogxK

U2 - 10.1021/acs.accounts.6b00607

DO - 10.1021/acs.accounts.6b00607

M3 - Article

VL - 50

SP - 1597

EP - 1605

JO - Accounts of Chemical Research

JF - Accounts of Chemical Research

SN - 0001-4842

IS - 7

ER -