In stable isotope geochemistry, the Urey–Bigeleisen–Mayer equation, also known as the Bigeleisen–Mayer equation or the Urey model,[1] is a model describing the approximate equilibrium isotope fractionation in an isotope exchange reaction.[2][3][4][5][6] While the equation itself can be written in numerous forms, it is generally presented as a ratio of partition functions of the isotopic molecules involved in a given reaction.[7][8] The Urey–Bigeleisen–Mayer equation is widely applied in the fields of quantum chemistry and geochemistry and is often modified or paired with other quantum chemical modelling methods (such as density functional theory) to improve accuracy and precision and reduce the computational cost of calculations.[1][6][9]
The equation was first introduced by Harold Urey and, independently, by Jacob Bigeleisen and Maria Goeppert Mayer in 1947.[2][7][8]
Since its original descriptions, the Urey–Bigeleisen–Mayer equation has taken many forms. Given an isotopic exchange reaction A + B ∗ = A ∗ + B {\displaystyle A+B^{*}=A^{*}+B} , such that ∗ {\displaystyle ^{*}} designates a molecule containing an isotope of interest, the equation can be expressed by relating the equilibrium constant, K e q {\displaystyle K_{eq}} , to the product of partition function ratios, namely the translational, rotational, vibrational, and sometimes electronic partition functions.[10][11][12] Thus the equation can be written as: K e q = [ A ∗ ] [ B ] [ A ] [ B ∗ ] {\displaystyle K_{eq}={\frac {[A^{*}][B]}{[A][B^{*}]}}} where [ A ] = ∏ n Q n , A {\displaystyle [A]=\prod ^{n}Q_{n,A}} and Q n {\displaystyle Q_{n}} is each respective partition function of molecule or atom A {\displaystyle A} .[12][13] It is typical to approximate the rotational partition function ratio as quantized rotational energies in a rigid rotor system.[11][14] The Urey model also treats molecular vibrations as simplified harmonic oscillators and follows the Born–Oppenheimer approximation.[11][14][15]
Isotope partitioning behavior is often reported as a reduced partition function ratio, a simplified form of the Bigeleisen–Mayer equation notated mathematically as s s ′ f {\displaystyle {\frac {s}{s'}}f} or ( Q ∗ Q ) r {\displaystyle ({\frac {Q^{*}}{Q}})_{r}} .[16][17] The reduced partition function ratio can be derived from power series expansion of the function and allows the partition functions to be expressed in terms of frequency.[16][18][19] It can be used to relate molecular vibrations and intermolecular forces to equilibrium isotope effects.[20]
As the model is an approximation, many applications append corrections for improved accuracy.[15] Some common, significant modifications to the equation include accounting for pressure effects,[21] nuclear geometry,[22] and corrections for anharmonicity and quantum mechanical effects.[1][2][23][24] For example, hydrogen isotope exchange reactions have been shown to disagree with the requisite assumptions for the model but correction techniques using path integral methods have been suggested.[1][8][25]
One aim of the Manhattan Project was increasing the availability of concentrated radioactive and stable isotopes, in particular 14C, 35S, 32P, and deuterium for heavy water.[26] Harold Urey, Nobel laureate physical chemist known for his discovery of deuterium,[27] became its head of isotope separation research while a professor at Columbia University.[28][29]: 45 In 1945, he joined The Institute for Nuclear Studies at the University of Chicago, where he continued to work with chemist Jacob Bigeleisen and physicist Maria Mayer, both also veterans of isotopic research in the Manhattan Project.[11][28][30][31] In 1946, Urey delivered the Liversidge lecture at the then-Royal Institute of Chemistry, where he outlined his proposed model of stable isotope fractionation.[2][7][11] Bigeleisen and Mayer had been working on similar work since at least 1944 and, in 1947, published their model independently from Urey.[2][8][11] Their calculations were mathematically equivalent to a 1943 derivation of the reduced partition function by German physicist Ludwig Waldmann.[8][11][a]
Initially used to approximate chemical reaction rates,[7][8] models of isotope fractionation are used throughout the physical sciences. In chemistry, the Urey–Bigeleisen–Mayer equation has been used to predict equilibrium isotope effects and interpret the distributions of isotopes and isotopologues within systems, especially as deviations from their natural abundance.[35][36] The model is also used to explain isotopic shifts in spectroscopy, such as those from nuclear field effects or mass independent effects.[1][22][35] In biochemistry, it is used to model enzymatic kinetic isotope effects.[37][38] Simulation testing in computational systems biology often uses the Bigeleisen–Mayer model as a baseline in the development of more complex models of biological systems.[39][40] Isotope fractionation modeling is a critical component of isotope geochemistry and can be used to reconstruct past Earth environments as well as examine surface processes.[41][42][43][44]
After this paper had been completed, Professor W.F. Libby kindly called a paper by L. Waldmann[32] to our attention. In this paper, Waldmann discusses briefly the fact that the chemical separation of isotopes is a quantum effect. He gives formulae which are equivalent to our (11') and (11a) and discusses qualitatively their application to two acid base exchange equilibria. These are the exchange between NH3 and NH4+ and HCN and CN- studies by Urey[33][34] and co-workers.