337 research outputs found
Structure, Energetics, and Dynamics of Cs+ and H2O in Hectorite: Molecular Dynamics Simulations with an Unconstrained Substrate Surface
Classical molecular dynamics simulations were performed for the smectite clay hectorite with charge-balancing Cs+ cations using a newly developed structural model with a disordered distribution of Li/Mg substitutions in the octahedral sheet and the fully flexible CLAYFF force field. Calculations for systems with interlayer galleries containing 0–19 H2O/Cs+ suggest that the monolayer hydrate is the only stable state at all relative humidities at ambient pressure and temperature, in agreement with experimental results and previous molecular calculations. The basal spacing of this structure is also in good agreement with experimental values. In contrast to previous molecular modeling results, however, the new simulations show that interlayer Cs+ occurs on 2 different inner sphere adsorption sites: above the center of ditrigonal cavities and above Si tetrahedra. Unlike previous simulations, which employed a rigid clay model and fixed orientations of the structural −OH groups, the present results are obtained for an unconstrained clay substrate structure, where the structural −OH groups are able to assume various orientations, including being nearly parallel to the clay layers. This flexibility allows the Cs+ ions to approach the surface more closely above the centers of the hexagonal rings. In this structural arrangement, Cs+ ions are not hydrated by the H2O molecules which share the same interlayer plane, but rather by the H2O molecules coordinated to the opposite surface. In contrast, on the external basal surface, a significant fraction of H2O molecules are adsorbed above the centers of ditrigonal cavities adjacent to adsorbed Cs+ ions. For these H2O molecules, both HH2O atoms coordinate and H-bond to Ob surface oxygen atoms. The mean residence times for the Cs+–H2O, Cs+–Ob, and H2O–Ob coordination pairs show that Cs+ ions are more strongly coordinated with Ob atoms than H2O molecules. This result is the opposite of the behavior in Ca-hectorite, due to the much smaller hydration energy of Cs+ compared to that of Ca2+
Cation and Water Structure, Dynamics, and Energetics in Smectite Clays: A Molecular Dynamics Study of Ca-Hectorite
The incorporation of Ca2+ into smectite minerals is well-known to have a significant effect on the swelling behavior and mechanical properties of this environmentally and technologically important group of materials. Relative to common alkali cations such as Na+, K+, and Cs+, Ca2+ has a larger charge/ionic radius ratio and thus interacts very differently with interlayer water molecules and the oxygens of the clay basal surface. Recent 2H and 43Ca NMR studies of the smectite mineral, hectorite, show that the molecular scale interlayer dynamics is quite different with Ca2+ than with alkali cations. Classical molecular dynamics (MD) simulations presented here use a newly developed hectorite model with a disordered distribution of Li+/Mg2+ substitutions in the octahedral sheet and provide new insight into the origin of the effects of Ca2+ on the structure, dynamics, and energetics of smectite interlayers. The computed basal spacings and thermodynamic properties suggest the potential for formation of stable monolayer hydrates that have partial and complete water contents, a bilayer hydrate, and possible expansion to higher hydration states. The system hydration energies are comparable to those previously calculated for Ca–montmorillonite and are more negative than for Cs– and Na–hectorite due to the higher hydration energy of Ca2+. The coordination environments of Ca2+ change significantly with increasing interlayer hydration, with the extent of coordination to basal oxygens decreasing as the number of interlayer molecules increases. On external (001) surfaces, the H2O molecules closest to the surface are adsorbed at the centers of ditrigonal cavities and bridge Ca2+ to the surface. The Ca2+ ions on the external surface are all in outer-sphere coordination with the basal oxygens of the surface, and the proximity-restricted region with a significant number of Ca2+ is approximately 6 Å thick. Quantification of these interactions provides a basis for understanding intercalation of Ca2+ by organic species and smectite minerals
Diffusion Behavior of Methane in 3D Kerogen Models
As global energy demand increases, natural gas recovery from source rocks is attracting considerable attention since recent development in shale extraction techniques has made the recovery process economically viable. Kerogens are thought to play an important role in gas recovery; however, the interactions between trapped shale gas and kerogens remain poorly understood due to the complex, heterogeneous microporous structure of kerogens. This study examines the diffusive behavior of methane molecules in kerogen matrices of different types (Type I, II, and II) and maturity levels (A to D for Type II kerogens) on a molecular scale. Models of each kerogen type were developed using simulated annealing. We employed grand canonical Monte Carlo simulations to predict the methane loadings of the kerogen models and then used equilibrium molecular dynamics simulations to compute the mean square displacement of methane molecules within the kerogen matrices under reservoir-relevant conditions, that is, 365 K and 275 bar. Our results show that methane self-diffusivity exhibits some degree of anisotropy in all kerogen types examined here except for Type I-A kerogens, where diffusion is the fastest and isotropic diffusion is observed. Self-diffusivity appears to correlate positively with pore volume for Type II kerogens, where an increase in diffusivity is observed with increasing maturity. Swelling of the kerogen matrix up to a 3% volume change is also observed upon methane adsorption. The findings contribute to a better understanding of hydrocarbon transport mechanisms in shale and may lead to further development of extraction techniques, fracturing fluids, and recovery predictions
Role of Cations in the Methane/Carbon Dioxide Partitioning in Nano- and Mesopores of Illite Using Constant Reservoir Composition Molecular Dynamics Simulation
We performed constant reservoir composition molecular dynamics (CRC-MD) simulations at 323 K and 124 bar to quantitatively study the partitioning of fluid species between the nano- and mesopores of clay and a bulk reservoir containing an equimolar mixture of CO2 and CH4. The results show that the basal (001) and protonated edge (010) surfaces of illite both demonstrate a strong preference for CO2 over CH4 adsorption; that the (001) surfaces show a stronger preference for CO2 than the (010) surfaces, especially with K+ as the exchangeable cation; and that the structuring of the near-surface CO2 by K+ is stronger than that by Na+. The protonated (010) surfaces have a somewhat greater preference for CH4, with the concentration near them close to that in the bulk fluid. The effects of the surfaces on the fluid composition extend to approximately 2.0 nm from them, with the fluid composition at the center of the pore becoming essentially the same as the bulk composition at a pore thickness of ∼5.7 nm. The preference of nano- and mesopores bounded by clay minerals for CO2 over CH4 suggests that injection of CO2 into tight reservoirs is likely to displace CH4 into larger pores, thus enhancing its production
Molecular Dynamics Study of CO2 and H2O Intercalation in Smectite Clays: Effect of Temperature and Pressure on Interlayer Structure and Dynamics in Hectorite
Grand Canonical Molecular Dynamics (GCMD) simulations were performed to
investigate the intercalation of CO2 and H2O molecules in the interlayers of the smectite clay,
Na-hectorite, at temperatures and pressures relevant to petroleum reservoir and geological carbon
sequestration conditions and in equilibrium with H2O-saturated CO2. The computed adsorption
isotherms indicate that CO2 molecules enter the interlayer space of Na-hectorite only when it is
hydrated with approximately 3 H2O molecules per unit cell. The computed immersion energies
show that the bilayer hydrate structure (2WL) contains less CO2 than the monolayer structure
(1WL), but that the 2WL hydrate is the most thermodynamically stable state, consistent with
experimental results for a similar Na-montmorillonite smectite. At all T and P conditions
examined (323-368 K and 90-150 bar), the CO2 molecules are adsorbed at the midplane of clay
interlayers for the 1WL structure and closer to one of the basal surfaces for the 2WL structure.
Interlayer CO2 molecules are dynamically less restricted in the 2WL structures. The CO2
molecules are preferentially located near basal surface oxygen atoms and H2O molecules rather
than in coordination with Na+
ions. Accounting for the orientation and flexibility of the structural
-OH groups of the clay layer has a significant effect on the details of the computed structure and
dynamics of H2O and CO2 molecules but does not affect the overall trends with changing basal
spacing or the principal structural and dynamical conclusions. Temperature and pressure in the
ranges examined have little effect on the principal structural and energetic conclusions, but the
rates of dynamical processes increase with increasing temperature, as expected
Understanding methane/carbon dioxide partitioning in clay nano- and meso-pores with constant reservoir composition molecular dynamics modeling
The interactions among fluid species such as H2O, CO2, and CH4 confined in nano- and meso-pores in shales and other rocks is of central concern to understanding the chemical behavior and transport properties of these species in the earth's subsurface and is of special concern to geological C-sequestration and enhanced production of oil and natural gas. The behavior of CO2, and CH4 is less well understood than that of H2O. This paper presents the results of a computational modeling study of the partitioning of CO2 and CH4 between bulk fluid and nano- and meso-pores bounded by the common clay mineral montmorillonite. The calculations were done at 323 K and a total fluid pressure of 124 bars using a novel approach (constant reservoir composition molecular dynamics, CRC-MD) that uses bias forces to maintain a constant composition in the fluid external to the pore. This purely MD approach overcomes the difficulties in making stochastic particle insertion-deletion moves in dense fluids encountered in grand canonical Monte Carlo and related hybrid approaches. The results show that both the basal siloxane surfaces and protonated broken edge surfaces of montmorillonite both prefer CO2 relative to CH4 suggesting that methods of enhanced oil and gas production using CO2 will readily displace CH4 from such pores. This preference for CO2 is due to its preferred interaction with the surfaces and extends to approximately 20 Ã… from them
Stellar winds from Massive Stars
We review the various techniques through which wind properties of massive
stars - O stars, AB supergiants, Luminous Blue Variables (LBVs), Wolf-Rayet
(WR) stars and cool supergiants - are derived. The wind momentum-luminosity
relation (e.g. Kudritzki et al. 1999) provides a method of predicting mass-loss
rates of O stars and blue supergiants which is superior to previous
parameterizations. Assuming the theoretical sqrt(Z) metallicity dependence,
Magellanic Cloud O star mass-loss rates are typically matched to within a
factor of two for various calibrations. Stellar winds from LBVs are typically
denser and slower than equivalent B supergiants, with exceptional mass-loss
rates during giant eruptions Mdot=10^-3 .. 10^-1 Mo/yr (Drissen et al. 2001).
Recent mass-loss rates for Galactic WR stars indicate a downward revision of
2-4 relative to previous calibrations due to clumping (e.g. Schmutz 1997),
although evidence for a metallicity dependence remains inconclusive (Crowther
2000). Mass-loss properties of luminous (> 10^5 Lo) yellow and red supergiants
from alternative techniques remain highly contradictory. Recent Galactic and
LMC results for RSG reveal a large scatter such that typical mass-loss rates
lie in the range 10^-6 .. 10^-4 Mo/yr, with a few cases exhibiting 10^-3 Mo/yr.Comment: 16 pages, 2 figures, Review paper to appear in Proc `The influence of
binaries on stellar population studies', Brussels, Aug 2000 (D. Vanbeveren
ed.), Kluwe
Tipping Point for Expansion of Layered Aluminosilicates in Weakly Polar Solvents: Supercritical CO2
Layered aluminosilicates play a dominant role in the mechanical and gas storage properties of the subsurface, are used in diverse industrial applications, and serve as model materials for understanding solvent-ion-support systems. Although expansion in the presence of H2O is well-known to be systematically correlated with the hydration free energy of the interlayer cation, particularly in environments dominated by nonpolar solvents (i.e., CO2), uptake into the interlayer is not well-understood. Using novel high-pressure capabilities, we investigated the interaction of dry supercritical CO2 with Na-, NH4-, and Cs-saturated montmorillonite, comparing results with predictions from molecular dynamics simulations. Despite the known trend in H2O and that cation solvation energies in CO2 suggest a stronger interaction with Na, both the NH4- and Cs-clays readily absorbed CO2 and expanded, while the Na-clay did not. The apparent inertness of the Na-clay was not due to kinetics, as experiments seeking a stable expanded state showed that none exists. Molecular dynamics simulations revealed a large endothermicity to CO2 intercalation in the Na-clay but little or no energy barrier for the NH4- and Cs-clays. Indeed, the combination of experiment and theory clearly demonstrate that CO2 intercalation of Na-montmorillonite clays is prohibited in the absence of H2O. Consequently, we have shown for the first time that in the presence of a low dielectric constant, gas swelling depends more on the strength of the interaction between the interlayer cation and aluminosilicate sheets and less on that with solvent. The finding suggests a distinct regime in layered aluminosilicate swelling behavior triggered by low solvent polarizability, with important implications in geomechanics, storage, and retention of volatile gases, and across industrial uses in gelling, decoloring, heterogeneous catalysis, and semipermeable reactive barriers
Saccharomyces cerevisiae mutants affected in vacuole assembly or vacuolar H+-ATPase are hypersensitive to lead (Pb) toxicity
Lead is an important environmental pollutant. The role of vacuole, in Pb detoxification, was studied using a vacuolar protein sorting mutant strain (vps16Δ), belonging to class C mutants. Cells disrupted in VPS16 gene, did not display a detectable vacuolar-like structure. Based on the loss of cell proliferation capacity, it was found that cells from vps16Δ mutant exhibited a hypersensitivity to Pb-induced toxicity, compared to wild type (WT) strain. The function of vacuolar H+-ATPase (V-ATPase), in Pb detoxification, was evaluated using mutants with structurally normal vacuoles but defective in subunits of catalytic (vma1Δ or vma2Δ) or membrane domain (vph1Δ or vma3Δ) of V-ATPase. All mutants tested, lacking a functional V-ATPase, displayed an increased susceptibility to Pb, comparatively to cells from WT strain. Modification of vacuolar morphology, in Pb-exposed cells, was visualized using a Vma2p-GFP strain. The treatment of yeast cells with Pb originated the fusion of the medium size vacuolar lobes into one enlarged vacuole. In conclusion, it was found that vacuole plays an important role in the detoxification of Pb in Saccharomyces cerevisiae; in addition, a functional V-ATPase was required for Pb compartmentalization.The authors thank the Fundacao para a Ciencia e a Tecnologia (FCT) through the Portuguese Government for their financial support of this work through the grant PEST-OE/EQB/LA0023/2011 to IBB
Parameterising a public good: how experiments on predation can be used to predict cheat frequencies
- …