Experimental determination of the Raman CH4 symmetric stretching (ν1) band position from 1-650 bar and 0.3-22 °C: Application to fluid inclusion studies

The position of the Raman methane (CH₄) symmetric stretching band (ν₁) over the range 1-650 bar and 0.3-22 °C has been determined using a high-pressure optical cell mounted on a Raman microprobe. Two neon emission lines that closely bracket the CH₄ band were collected simultaneously with each CH₄ spectrum. The peak position was determined after least squares fitting using a summed Gaussian-Lorentzian method, resulting in a precision of ≈±0.02 cm⁻¹ in peak position determination. The CH₄ν₁ band position shifts to lower wave number with increasing pressure. At a given pressure, the band shifts to lower wave number with decreasing temperature, and the magnitude of the temperature shift increases with increasing pressure. The relationship between the Raman CH₄ν₁ band position and temperature and pressure determined here may be used to estimate the internal pressure in natural or synthetic CH₄-bearing fluid inclusions. This information, in turn, may be used to determine the density of pure CH₄ fluid inclusions and the salinity of CH₄-bearing aqueous inclusions.     

Thermodynamic modeling of binary CH4-H2O fluid inclusions

This work reports the application of thermodynamic models, including equations of state, to binary (salt-free) CH₄-H₂O fluid inclusions. A general method is presented to calculate the compositions of CH₄-H₂O inclusions using the phase volume fractions and dissolution temperatures of CH₄ hydrate. To calculate the homogenization pressures and isolines of the CH₄-H₂O inclusions, an improved activity-fugacity model is developed to predict the vapor-liquid phase equilibrium. The phase equilibrium model can predict methane solubility in the liquid phase and water content in the vapor phase from 273 to 623 K and from 1 to 1000 bar (up to 2000 bar for the liquid phase), within or close to experimental uncertainties. Compared to reliable experimental phase equilibrium data, the average deviation of the water content in the vapor phase and methane solubility in the liquid phase is 4.29% and 3.63%, respectively. In the near-critical region, the predicted composition deviations increase to over 10%. The vapor-liquid phase equilibrium model together with the updated volumetric model of homogenous (single-phase) CH₄-H₂O fluid mixtures (Mao S., Duan Z., Hu J. and Zhang D. (2010) A model for single-phase PVTx properties of CO₂-CH₄-C₂H₆-N₂-H₂O-NaCl fluid mixtures from 273 to 1273 K and from 1 to 5000 bar. Chem. Geol.275, 148-160), is applied to calculate the isolines, homogenization pressures, homogenization volumes, and isochores at specified homogenization temperatures and compositions. Online calculation is on the website: http://www.geochem-model.org/.     

Determination of methane concentrations in water in equilibrium with sI methane hydrate in the absence of a vapor phase by in situ Raman spectroscopy

Most submarine gas hydrates are located within the two-phase equilibrium region of hydrate and interstitial water with pressures (P) ranging from 8 to 60 MPa and temperatures (T) from 275 to 293 K. However, current measurements of solubilities of methane in equilibrium with hydrate in the absence of a vapor phase are limited below 20 MPa and 283.15 K, and the differences among these data are up to 30%. When these data were extrapolated to other P-T conditions, it leads to large and poorly known uncertainties. In this study, in situ Raman spectroscopy was used to measure methane concentrations in pure water in equilibrium with sI (structure one) methane hydrate, in the absence of a vapor phase, at temperatures from 276.6 to 294.6 (±0.3) K and pressures at 10, 20, 30 and 40 (±0.4%) MPa. The relationship among concentration of methane in water in equilibrium with hydrate, in mole fraction [X(CH₄)], the temperature in K, and pressure in MPa was derived as: X(CH₄) = exp [11.0464 + 0.023267 P − (4886.0 + 8.0158 P)/T]. Both the standard enthalpy and entropy of hydrate dissolution at the studied T-P conditions increase slightly with increasing pressure, ranging from 41.29 to 43.29 kJ/mol and from 0.1272 to 0.1330 kJ/K · mol, respectively. When compared with traditional sampling and analytical methods, the advantages of our method include: (1) the use of in situ Raman signals for methane concentration measurements eliminates possible uncertainty caused by sampling and ex situ analysis, (2) it is simple and efficient, and (3) high-pressure data can be obtained safely.     

Raman spectroscopic measurements of CO2 density: Experimental calibration with high-pressure optical cell (HPOC) and fused silica capillary capsule (FSCC) with application to fluid inclusion observations

Raman spectroscopy is a powerful method for the determination of CO₂ densities in fluid inclusions, especially for those with small size and/or low fluid density. The relationship between CO₂ Fermi diad split (Δ, cm⁻¹) and CO₂ density (ρ, g/cm³) has been documented by several previous studies. However, significant discrepancies exist among these studies mainly because of inconsistent calibration procedures and lack of measurements for CO₂ fluids having densities between 0.21 and 0.75 g/cm³, where liquid and vapor phases coexist near room temperature.In this study, a high-pressure optical cell and fused silica capillary capsules were used to prepare pure CO₂ samples with densities between 0.0472 and 1.0060 g/cm³. The measured CO₂ Fermi diad splits were calibrated with two well established Raman bands of benzonitrile at 1192.6 and 1598.9 cm⁻¹. The relationship between the CO₂ Fermi diad split and density can be represented by: ρ = 47513.64243 − 1374.824414 × Δ + 13.25586152 × Δ² − 0.04258891551 × Δ³ (r² = 0.99835, σ = 0.0253 g/cm³), and this relationship was tested by synthetic fluid inclusions and natural CO₂-rich fluid inclusions. The effects of temperature and the presence of H₂O and CH₄ on this relationship were also examined.     

The system H2O-NaCl. Part I: Correlation formulae for phase relations in temperature-pressure-composition space from 0 to 1000 °C, 0 to 5000 bar, and 0 to 1 XNaCl

Realistic simulations of fluid flow in geologic systems have severely been hampered by the lack of a consistent formulation for fluid properties for binary salt-water fluids over the temperature-pressure-composition ranges encountered in the Earth’s crust. As the first of two companion studies, a set of correlations describing the phase stability relations in the system H₂O-NaCl is developed. Pure water is described by the IAPS-84 equation of state. New correlations comprise the vapor pressure of halite and molten NaCl, the NaCl melting curve, the composition of halite-saturated liquid and vapor, the pressure of vapor + liquid + halite coexistence, the temperature-pressure and temperature-composition relations for the critical curve, and the compositions of liquid and vapor on the vapor + liquid coexistence surface. The correlations yield accurate values for temperatures from 0 to 1000 °C, pressures from 0 to 5000 bar, and compositions from 0 to 1 X_NaCl (mole fraction of NaCl). To facilitate their use in fluid flow simulations, the correlations are entirely formulated as functions of temperature, pressure and composition.     

Strain response and re-equilibration of CH4-rich synthetic aqueous fluid inclusions in calcite during pressure drops

Aqueous fluids in sedimentary basins often contain dissolved methane, particularly in petroleum environments. PVTX (Pressure-Volume-Temperature-Composition) reconstructions performed using fluid inclusion data are largely based on the assumption that inclusions do not change from the time of trapping until the present. Many authors, however, consider that fluid inclusions can re-equilibrate, particularly in fragile minerals like calcite. In order to understand this re-equilibration phenomenon in the metamorphic domain, previous experiments have been performed under high PT conditions, but few have been performed at low to medium PT conditions such as those associated with sedimentary burial diagenesis, and no previous studies have examined CH₄-bearing aqueous inclusions in calcite.An experimental study of the preservation/modification of CH₄-rich synthetic fluid inclusions in calcite during isothermal decompression was conducted. An autoclave was used for accurate PTX control allowing equilibrium between liquid and vapour in the CH₄-H₂O system. PTX conditions were maintained at four stages of decreasing pressure, with each stage held for 7 days to simulate an isothermal pressure drop. In order of decreasing pressure, the pressure-temperature conditions monitored were 276 ± 10 bar at 180 ± 7 °C, 176 ± 10 bar at 180 ± 7 °C, 76 ± 10 bar at 180 ± 7 °C and 10 ± 3 bar at 180 ± 15 °C. At the end of the experiment, the calcite was recovered and analyzed by microthermometry and Raman microspectroscopy for PTX reconstruction. A careful procedure was adopted to limit re-equilibration of inclusions during analytical procedures. Four types of inclusion shapes and four types of strain patterns were differentiated. Classification of the petrographic strain patterns was carried out. These strain patterns were associated with inclusion stretching and/or leakage regarding CH₄, T_h and P_h compared to experimental conditions. Factors controlling the preservation or acquisition of strain patterns included the initial shape and size of the inclusion, and the pressure differential (ΔP) between the confining pressure (P_cf) and the internal pressure (P_i) within the inclusion. Most fluid inclusions seemed to be trapped during the first 7 days of the experiment, although few (4%) of these preserved the initial PT conditions of 276 ± 10 bar, whereas 8% preserved the second and third run of PT conditions. Overall, the majority of inclusions (88%) did not reflect accurately the PTX trapping conditions. A petrographic guide to the inclusions is presented here that allows strain identification for PVT reconstructions. Re-equilibration patterns and evidence for preferential methane leakage from aqueous inclusions in calcite are important findings revealed by this study, and may be useful for the reconstruction of post-trapping events in investigations of natural samples, and in other experiments using synthetic inclusions in calcite.     

Formation and abundance of doubly-substituted methane isotopologues (CH3D) in natural gas systems

Formation of the Carbon-13 (¹³C) and deuterium (D) doubly-substituted methane isotopologues (¹³CH₃D) in natural gases is studied utilizing both first-principle quantum mechanism molecular calculation and direct FTIR laboratorial measurements of specifically synthesized high isotope concentration methane gas. For ¹³CH₃D, the symmetrically breathing mode A₀ emerges as IR-detectable attributed to the molecular symmetry lowering to C_3v from T_d of the non-isotopic methane (CH₄), along with a large vibrational frequency shift from ∼3000 to ∼2250 cm⁻¹. Our studies also indicate that the concentration of ¹³CH₃D is dependent on the environmental temperature through isotope exchanges among methane isotopologues; and the Gibbs’ Free Energy difference due to Quantum Mechanics Zero-Point vibrational motions has the major contribution to this temperature dependency. Potential geologic applications of the ¹³CH₃D measurement to natural gas exploration and assessments are also discussed. In order to detect the ¹³CH₃D concentration change of each 50 °C in the natural gas system, a 10⁻⁹ resolution is desirable. Such a measurement could provide important add-on information to distinguish natural gas origin and distribution.     

Molecular dynamics simulation of the CH4 and CH4-H2O systems up to 10 GPa and 2573 K

Based on our previous study of the intermolecular potential for pure H₂O and the strict evaluation of the competitive potential models for pure CH₄ and the ab initio fitting potential surface across CH₄-H₂O molecules in this study, we carried out more than two thousand molecular dynamics simulations for the PVTx properties of pure CH₄ and the CH₄-H₂O mixtures up to 2573 K and 10 GPa. Comparison of 1941 simulations with experimental PVT data for pure CH₄ shows an average deviation of 0.96% and a maximum deviation of 2.82%. The comparison of the results of 519 simulations of the mixtures with the experimental measurements reveals that the PVTx properties of the CH₄-H₂O mixtures generally agree with the extensive experimental data with an average deviation of 0.83% and 4% in maximum, which is equivalent to the experimental uncertainty. Moreover, the maximum deviation between the experimental data and the simulation results decreases to about 2% as temperature and pressure increase, indicating that the high accuracy of the simulation is well retained in the high temperature and pressure region.After the validation of the simulation method and the intermolecular potential models, we systematically simulated the PVTx properties of this binary system from 673 K and 0.05 GPa to 2573 K and 10 GPa. In order to integrate all the simulation results and the experimental data for the calculation of thermodynamic properties, an equation of state (EOS) is developed for the CH₄-H₂O system covering 673-2573 K and 0.01-10 GPa. Isochores for compositions <4 mol% CH₄ up to 773 K and 600 MPa are also determined in this paper. The program for the EOS can be downloaded from www.geochem-model.org/programs.htm.     

Raman spectroscopic study of CO2-NaCl-H2O mixtures in synthetic fluid inclusions at high temperatures

Mixtures of CO₂-NaCl-H₂O contained in synthetic fluid inclusions are studied by laser Raman spectroscopy at high temperatures. With increasing temperature, the band splitting (X) of υ₁-2υ₂ diad of spectrum of CO₂ presents more variations, and the intensity ratios of the hot bands to the υ₁-2υ₂ diad increase. For mixtures of gas phase rich in CO₂ and liquid phase rich in H₂O before homogenization, the strength of hydrogen bonding of water in the liquid phase decreases almost linearly with increasing temperature. When mixtures become homogeneous liquid phases, carbon dioxide content increases significantly as a result of the abrupt decrease in hydrogen bonds. Our results show that the hydrogen bonds change only slightly at higher temperatures above the homogeneous point, and a certain extent of the hydrogen bonds still remains at the highest temperature of 550°C of this work. The study is helpful to Raman spectroscopic analysis of natural fluid inclusions at high temperatures.     

Experimental study of quartz inclusions in garnet at pressures up to 3.0 GPa: evaluating validity of the quartz-in-garnet inclusion elastic thermobarometer

Garnet crystals with quartz inclusions were hydrothermally crystallized from oxide starting materials in piston–cylinder apparatuses at pressures from 0.5 to 3 GPa and temperatures ranging from 700 to 800 °C to study how entrapment conditions affect remnant pressures of quartz inclusions used for quartz-in-garnet (QuiG) elastic thermobarometry. Systematic changes of the 128, 206 and 464 cm^?1 Raman band frequencies of quartz were used to determine pressures of quartz inclusions in garnet using Raman spectroscopy calibrations that describe the P–T dependencies of Raman band shifts for quartz under hydrostatic pressure. Within analytical uncertainties, inclusion pressures calculated for each of the three Raman band frequencies are equivalent, which suggests that non-hydrostatic stress effects caused by elastic anisotropy in quartz are smaller than measurement errors. The experimental quartz inclusions have pressures ranging from ??0.351 to 1.247 GPa that span the range of values observed for quartz inclusions in garnets from natural rocks. Quartz inclusion pressures were used to model P–T conditions at which the inclusions could have been trapped. The accuracy of QuiG thermobarometry was evaluated by considering the differences between pressures measured during experiments and pressures calculated using published equation of state parameters for quartz and garnet. Our experimental results demonstrate that Raman measurements performed at room temperature can be used without corrections to estimate garnet crystallization pressures. Calculated entrapment pressures for quartz inclusions in garnet are less than ~?10% different from pressures measured during the experiments. Because the method is simple to apply with reasonable accuracy, we expect widespread usage of QuiG thermobarometry to estimate crystallization conditions for garnet-bearing silicic rocks.     

Hydroxyl stretching in phyllosilicates at high pressures and temperatures: an infrared spectroscopic study

The hydroxyl stretching frequencies of four phyllosilicates have been measured at high pressures and temperatures using an externally heated diamond-anvil cell and synchrotron infrared spectroscopy. Spectra were measured up to 26, 31, 21 and 8 GPa at room temperature for samples of talc, pyrophyllite, muscovite and 10-Å phase, respectively. Spectra were also measured in the range 273–500 K at ambient pressure for all samples and at 8–9 GPa for talc and pyrophyllite. The frequency of the Mg₃OH band in talc increases with pressure due to the absence of hydrogen bonding. The different orientation of the hydroxyl group in pyrophyllite and muscovite leads to hydrogen bonding and a decrease in the frequency of the Al₂OH band with pressure. 10-Å phase is approximately equivalent to talc with the addition of interlayer H₂O. In a spectrum of a sample synthesised for 143 h, two hydroxyl stretching bands are clearly resolved on compression. One is the same as the Mg₃OH band in talc, indicating the presence of intra-layer hydroxyl in a talc-like environment with no hydrogen bonding. The other, which separates from the talc-like band at 1 GPa, is associated with intra-layer hydroxyl that is hydrogen bonded to interlayer H₂O. There are equivalent bands in high-pressure spectra of a sample of deuterated 10-Å phase, synthesised for 400 h. This sample shows a greater extent of hydrogen bonding at ambient pressure than the 143 h sample. For all of the phases studied, increasing temperature leads to a decrease in frequency for every hydroxyl stretching vibration, both at low and high pressures. The shifts in frequency with temperature are an order of magnitude greater than the shifts with pressure when normalised to previously measured structural parameters.     

Combined microthermometric and Raman spectroscopic technique to determine the salinity of H2O-CO2-NaCl fluid inclusions based on clathrate melting

Fluid inclusions approximated by the system H₂O-CO₂-NaCl are common in many geologic environments. In order to apply microthermometric data from these inclusions to infer P-T (pressure-temperature) trapping conditions, the composition of the inclusions, including the salinity, must be known. Normally, salinities of aqueous inclusions are determined from ice-melting temperatures obtained during microthermometry. However, when CO₂-bearing aqueous fluid inclusions are cooled they often form a hydrate that incorporates H₂O into the structure, and salinities estimated from ice-melting temperatures are therefore higher than the actual salinity. A technique that combines data from Raman spectroscopic and microthermometric analyses of individual inclusions was developed to determine the salinity of CO₂-bearing aqueous inclusions based on measured clathrate melting temperatures and CO₂ pressures obtained from Raman analyses. In this study, the pressure within inclusions was determined using Raman spectroscopy based on the splitting of the Fermi diad of CO₂, measured at the clathrate melting temperature. The CO₂ densities (and pressures) predicted by the equation developed in this study are in relatively good agreement with previously published equations, except for very low densities and correspondingly low pressures. The combined Raman spectroscopy - microthermometry technique thus provides both the temperature and the pressure in the inclusion at clathrate melting. For inclusions in which the clathrate melts in the presence of CO₂ liquid, the salinity can be determined with a precision of a few tenths of a wt% NaCl, whereas for inclusions in which clathrate melts in the presence of CO₂ vapor the salinity error may be a few wt% NaCl. Applying the method to synthetic fluid inclusions with known salinity suggests that the technique is valid for determining salinity of H₂O-CO₂-NaCl fluid inclusions in which clathrate melts in the presence of liquid CO₂ only or vapor CO₂ only.     

Synthesis of Fluid Inclusion Reference Materials in the H2O–NaCl–CH4 System

In this study, we detail procedures to synthesise inclusions of fluids of known composition under diagenetic conditions – low temperatures (50–180 °C) and pressures (2–25 MPa) – hosted by quartz. Pre‐existing microcavities in polished quartz chips were used, decrepitated and rinsed before being healed in an autoclave. Each synthesis required around 4 weeks. A range of biphase liquid–vapour inclusions in the H₂O–CH₄–NaCl system, similar to natural ones, were synthesised in quartz, varying the content of all three components. The quartz chips were checked before and after synthesis by optical microscopy, microthermometry and Raman spectroscopy. The method was used to build a series of calibration curves relating methane‐to‐water count ratios to the methane content of aqueous fluid inclusions. Good correlation was obtained between Raman band area ratios and the expected CH₄ solubilities from equations of state. The chosen procedure's main advantages, compared with other published methods, are as follows: quicker healing of the microfractures at lower temperatures, larger inclusions, simple quartz sample preparation, control of fluid inclusion composition and the use of commercial reactors. Such reference material inclusions allow the calibration of Raman microspectrometers with the objective of PTX reconstruction of palaeofluid migration in petroleum reservoirs.     

Bromine cycle in subduction zones through in situ Br monitoring in diamond anvil cells

The geochemical partitioning of bromine between hydrous haplogranitic melts, initially enriched with respect to Br and aqueous fluids, has been continuously monitored in situ during decompression. Experiments were carried out in diamond anvil cells from 890 °C to room temperature and from 1.7 GPa to room pressure, typically from high P, T conditions corresponding to total miscibility (presence of a supercritical fluid). Br contents were measured in aqueous fluids, hydrous melts and supercritical fluids. Partition coefficients of bromine were characterized at pressure and temperature between fluids, hydrous melts and/or glasses, as appropriate: D^Br_fluid/melt = (Br)_fluid/(Br)_melt, ranges from 2.18 to 9.2 ± 0.5 for conditions within the ranges 0.66-1.7 GPa, 590-890 °C; and D^Br_fluid/glass = (Br)_fluid/(Br)_glass ranges from 60 to 375 at room conditions. The results suggest that because high pressure melts and fluids are capable of accepting high concentrations of bromine, this element may be efficiently removed from the slab to the mantle source of arc magmas. We show that Br may be highly concentrated in subduction zone magmas and strongly enriched in subduction-related volcanic gases, because its mobility is strongly correlated with that of water during magma degassing. Furthermore, our experimental results suggest that a non negligible part of Br present in the subducted slab may remain in the down-going slab, being transported toward the transition zone. This indicates that the Br cycle in subduction zones is in fact divided in two related but independent parts: (1) a shallower one where recycled Br may leave the slab with a water and silica-bearing “fluid” leading to enriched arc magmas that return Br to the atmosphere. (2) A deeper cycle where Br may be recycled back to the mantle maybe to the transition zone, where it may be present in high pressure water-rich metasomatic fluids.     

Thermodynamics, structure, dynamics, and freezing of Mg2SiO4 liquid at high pressure

We perform first principles molecular dynamics simulations of Mg₂SiO₄ liquid and crystalline forsterite. On compression by a factor of two, we find that the Grüneisen parameter of the liquid increases linearly from 0.6 to 1.2. Comparison of liquid and forsterite equations of state reveals a temperature-dependent density crossover at pressures of ∼12-17 GPa. Along the melting curve, which we calculate by integration of the Clapeyron equation, the density crossover occurs within the forsterite stability field at P = 13 GPa and T = 2550 K. The melting curve obtained from the root mean-square atomic displacement in forsterite using the Lindemann law fails to match experimental or calculated melting curves. We attribute this failure to the liquid structure that differs significantly from that of forsterite, and which changes markedly upon compression, with increases in the degree of polymerization and coordination. The mean Si coordination increases from 4 in the uncompressed system to 6 upon twofold compression. The self-diffusion coefficients increase with temperature and decrease monotonically with pressure, and are well described by the Arrhenian relation. We compare our equation of state to the available highpressure shock wave data for forsterite and wadsleyite. Our theoretical liquid Hugoniot is consistent with partial melting along the forsterite Hugoniot at pressures 150-170 GPa, and complete melting at 170 GPa. The wadsleyite Hugoniot is likely sub-liquidus at the highest experimental pressure to date (200 GPa).     

Speciation of reduced C-O-H volatiles in coexisting fluids and silicate melts determined in-situ to ∼1.4 GPa and 800 °C

The structure of silicate melts in the system Na₂O·4SiO₂ saturated with reduced C-O-H volatile components and of coexisting silicate-saturated C-O-H solutions has been determined in a hydrothermal diamond anvil cell (HDAC) by using confocal microRaman and FTIR spectroscopy as structural probes. The experiments were conducted in-situ with the melt and fluid at high temperature (up to 800 °C) and pressure (up to 1435 MPa). Redox conditions in the HDAC were controlled with the reaction, Mo + H₂O = MoO₂ + H₂, which is slightly more reducing than the Fe + H₂O = FeO + H₂ buffer at 800 °C and less.The dominant species in the fluid are CH₄ + H₂O together with minor amounts of molecular H₂ and an undersaturated hydrocarbon species. In coexisting melt, CH₃ - groups linked to the silicate melt structure via Si-O-CH₃ bonding may dominate and possibly coexists with molecular CH₄. The abundance ratio of CH₃ - groups in melts relative to CH₄ in fluids increases from 0.01 to 0.07 between 500 and 800 °C. Carbon-bearing species in melts were not detected at temperatures and pressures below 400 °C and 730 MPa, respectively. A schematic solution mechanism is, Si-O-Si + CH₄?Si-O-CH₃+H-O-Si. This mechanism causes depolymerization of silicate melts. Solution of reduced (C-O-H) components will, therefore, affect melt properties in a manner resembling dissolved H₂O.     

Structure, thermodynamics, and diffusion in CaAl2Si2O8 liquid from first-principles molecular dynamics

We have performed first-principles molecular dynamics simulations of CaAl₂Si₂O₈ (anorthite) liquid at pressures up to 120 GPa and temperatures of 3000, 4000 and 6000 K. At the lowest degrees of compression the liquid is seen to accommodate changes in density through decreasing the abundance of 3- and 4-membered rings, while increases in coordination of network forming cations take effect at somewhat higher degrees of compression. Results are fit to a fundamental thermodynamic relation with 4th order finite strain and 1st order thermal variable expansions. Upon compression by a factor of two, the Grüneisen parameter (γ) is found to increase continuously from 0.35 to 1.10. Weak temperature dependence in γ is thermodynamically consistent with a slight decrease in isochoric heat capacity (C_V), for which values of between 4.4 and 5.2 Nk_B are obtained, depending on the temperature. Pressure and temperature dependence of self-diffusivities is found to be well represented by an Arrhenius relation, except at 3000 K and pressures lower than 5 GPa, where self-diffusivities of Si, Al, and O increase with pressure. Analysis of the lifetimes of individual coordination species reveals that this phenomenon arises due to the disproportionately high stability of 4-fold coordinated Si, and to a lesser extent 4-fold coordinated Al. Our results represent a marked improvement in accuracy and reliability in describing the physics of CaAl₂Si₂O₈ liquid at deep mantle pressures, pointing the way to a general thermodynamic model of melts at extreme pressures and temperatures relevant to planetary-scale magma oceans and deep mantle partial melting.     

Variable carbon isotope fractionation expressed by aerobic CH4-oxidizing bacteria

Carbon isotope fractionation factors reported for aerobic bacterial oxidation of CH₄(α_CH4-CO2) range from 1.003 to 1.039. In a series of experiments designed to monitor changes in the carbon isotopic fractionation of CH₄ by Type I and Type II methanotrophic bacteria, we found that the magnitude of fractionation was largely due to the first oxidation step catalyzed by methane monooxygenase (MMO). The most important factor that modulates the (α_CH4-CH3OH) is the fraction of the total CH₄ oxidized per unit time, which strongly correlates to the cell density of the growth cultures under constant flow conditions. At cell densities of less than 0.1 g/L, fractionation factors greater than 1.03 were observed, whereas at cell densities greater than 0.5 g/L the fractionation factors decreased to as low as 1.002. At low cell densities, low concentrations of MMO limit the amount of CH₄ oxidized, while at higher cell densities, the overall rates of CH₄ oxidation increase sufficiently that diffusion of CH₄ from the gaseous to dissolved state and into the cells is likely the rate-determining step. Thus, the residual CH₄ is more fractionated at low cell densities, when only a small fraction of the total CH₄ has been oxidized, than at high cell densities, when up to 40% of the influent CH₄ has been utilized. Therefore, since Rayleigh distillation behavior is not observed, δ¹³C values of the residual CH₄ cannot be used to infer the amount oxidized in either laboratory or field-studies. The measured (α_CH4-CH3OH) was the same for both Type I and Type II methanotrophs expressing particulate or soluble MMO. However, large differences in the δ¹³C values of biomass produced by the two types of methanotrophs were observed. Methylosinus trichosporium OB3b (Type II) produced biomass with δ¹³C values about 15‰ higher than the dissimilated CO₂, whereas Methylomonas methanica (Type I) produced biomass with δ¹³C values only about 6‰ higher than the CO₂. These effects were independent of the magnitude of the initial carbon isotope fractionation caused by MMO and were relatively constant despite changing ratios of assimilatory to dissimilatory carbon transformation by the organisms. This suggests that the difference in biomass carbon isotopes is primarily due to differences in the fractionation effect at the formaldehyde branch point in the metabolic pathway, rather than assimilation of CO₂ by Type II methanotrophs.     

An algorithm for finding composition,molar volume and isochors of CO2-CH4 fluid inclusions from Th and Tfm (for Th < Tfm)

A modified Redlich-Kwong equation of state is used to calculate the solubility of CO₂ in methane at various temperatures and pressures. From the solubility of CO₂ in CH₄ at the triple point and at final melting (T_h < T_fm), and the molar volume of solid CO₂, the volume of solid at the triple point, and the molar volume of the inclusion can be calculated using a mass balance. The pressure at the melting point is calculated from the equation of state.The algorithm predicts composition, molar volume, pressure at final melting and the isochor pressure (for a given temperature of trapping) for CO₂-CH₄ fluid inclusions for the case T_h < T_fm, given T_h, T_fm and experimental data on P_h and d_co2 (solid) at T_h.     

The system H2O-NaCl. Part II: Correlations for molar volume, enthalpy, and isobaric heat capacity from 0 to 1000 °C, 1 to 5000 bar, and 0 to 1 XNaCl

A set of correlations for the volumetric properties and enthalpies of phases in the system H₂O-NaCl as a function of temperature, pressure, and composition has been developed that yields accurate values from 0 to 1000 °C, 1 to 5000 bar, and 0 to 1 X_NaCl. The volumetric properties of all fluid phases from low-density vapor to hydrous salt melts and single-phase binary fluids at high pressures and temperatures, can be described by a simple equation