Low-temperature heat capacities,entropies and high-pressure phase relations of MgSiO<Subscript>3</Subscript> ilmenite and perovskite

Low-temperature isobaric heat capacities (C _p ) of MgSiO₃ ilmenite and perovskite were measured in the temperature range of 1.9–302.4 K with a thermal relaxation method using the Physical Properties Measurement System. The measured C _p of perovskite was higher than that of ilmenite in the whole temperature range studied. From the measured C _p , standard entropies at 298.15 K of MgSiO₃ ilmenite and perovskite were determined to be 53.7 ± 0.4 and 57.9 ± 0.3 J/mol K, respectively. The positive entropy change (4.2 ± 0.5 J/mol K) of the ilmenite–perovskite transition in MgSiO₃ is compatible with structural change across the transition in which coordination of Mg atoms is changed from sixfold to eightfold. Calculation of the ilmenite–perovskite transition boundary using the measured entropies and published enthalpy data gives an equilibrium transition boundary at about 20–23 GPa at 1,000–2,000 K with a Clapeyron slope of −2.4 ± 0.4 MPa/K at 1,600 K. The calculated boundary is almost consistent within the errors with those determined by high-pressure high-temperature in situ X-ray diffraction experiments.     

Peierls dislocation modelling in perovskite (CaTiO<Subscript>3</Subscript>): comparison with tausonite (SrTiO<Subscript>3</Subscript>) and MgSiO<Subscript>3</Subscript> perovskite

We present here a numerical modelling study of dislocations in perovskite CaTiO₃. The dislocation core structures and properties are calculated through the Peierls–Nabarro model using the generalized stacking fault (GSF) results as a starting model. The GSF are determined from first-principles calculations using the VASP code. The dislocation properties such as collinear, planar core spreading and Peierls stresses are determined for the following slip systems: [100](010), [100](001), [010](100), [010](001), [001](100), [001](010), and All dislocations exhibit lattice friction, but glide appears to be easier for [100](010) and [010](100). [001](010) and [001](100) exhibit collinear dissociation. Comparing Peierls stresses among tausonite (SrTiO₃), perovskite (CaTiO₃) and MgSiO₃ perovskite demonstrates the strong influence of orthorhombic distortions on lattice friction. However, and despite some quantitative differences, CaTiO₃ appears to be a satisfactory analogue material for MgSiO₃ perovskite as far as dislocation glide is concerned.     

High-pressure phase transitions of CaRhO<Subscript>3</Subscript> perovskite

High-pressure phase transitions of CaRhO₃ perovskite were examined at pressures of 6–27 GPa and temperatures of 1,000–1,930°C, using a multi-anvil apparatus. The results indicate that CaRhO₃ perovskite successively transforms to two new high-pressure phases with increasing pressure. Rietveld analysis of powder X-ray diffraction data indicated that, in the two new phases, the phase stable at higher pressure possesses the CaIrO₃-type post-perovskite structure (space group Cmcm) with lattice parameters: a = 3.1013(1) Å, b = 9.8555(2) Å, c = 7.2643(1) Å, V _m  = 33.43(1) cm³/mol. The Rietveld analysis also indicated that CaRhO₃ perovskite has the GdFeO₃-type structure (space group Pnma) with lattice parameters: a = 5.5631(1) Å, b = 7.6308(1) Å, c = 5.3267(1) Å, V _m  = 34.04(1) cm³/mol. The third phase stable in the intermediate P, T conditions between perovskite and post-perovskite has monoclinic symmetry with the cell parameters: a = 12.490(3) Å, b = 3.1233(3) Å, c = 8.8630(7) Å, β = 103.96(1)°, V _m  = 33.66(1) cm³/mol (Z = 6). Molar volume changes from perovskite to the intermediate phase and from the intermediate phase to post-perovskite are –1.1 and –0.7%, respectively. The equilibrium phase relations determined indicate that the boundary slopes are large positive values: 29 ± 2 MPa/K for the perovskite—intermediate phase transition and 62 ± 6 MPa/K for the intermediate phase—post-perovskite transition. The structural features of the CaRhO₃ intermediate phase suggest that the phase has edge-sharing RhO₆ octahedra and may have an intermediate structure between perovskite and post-perovskite.     

Elastic behavior of vanadinite,Pb<Subscript>10</Subscript>(VO<Subscript>4</Subscript>)<Subscript>6</Subscript>Cl<Subscript>2</Subscript>, a microporous non-zeolitic mineral

The high-pressure behavior of a vanadinite (Pb₁₀(VO₄)₆Cl₂, a = b = 10.3254(5), c = 7.3450(4) Å, space group P6₃/m), a natural microporous mineral, has been investigated using in-situ HP-synchrotron X-ray powder diffraction up to 7.67 GPa with a diamond anvil cell under hydrostatic conditions. No phase transition has been observed within the pressure range investigated. Axial and volume isothermal Equations of State (EoS) of vanadinite were determined. Fitting the P–V data with a third-order Birch-Murnaghan (BM) EoS, using the data weighted by the uncertainties in P and V, we obtained: V ₀ = 681(1) Å³, K ₀ = 41(5) GPa, and K′ = 12.5(2.5). The evolution of the lattice constants with P shows a strong anisotropic compression pattern. The axial bulk moduli were calculated with a third-order “linearized” BM-EoS. The EoS parameters are: a ₀ = 10.3302(2) Å, K ₀(a) = 35(2) GPa and K′(a) = 10(1) for the a-axis; c ₀ = 7.3520(3) Å, K ₀(c) = 98(4) GPa, and K′(c) = 9(2) for the c-axis (K ₀(a):K ₀(c) = 1:2.80). Axial and volume Eulerian-finite strain (fe) at different normalized stress (Fe) were calculated. The weighted linear regression through the data points yields the following intercept values: Fe _a (0) = 35(2) GPa for the a-axis, Fe _c (0) = 98(4) GPa for the c-axis and Fe _V (0) = 45(2) GPa for the unit-cell volume. The slope of the regression lines gives rise to K′ values of 10(1) for the a-axis, 9(2) for the c-axis and 11(1) for the unit cell-volume. A comparison between the HP-elastic response of vanadinite and the iso-structural apatite is carried out. The possible reasons of the elastic anisotropy are discussed.     

First-principles simulation of high-pressure polymorphs in MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript>

We have used density functional theory to investigate the stability of MgAl₂O₄ polymorphs under pressure. Our results can reasonably explain the transition sequence of MgAl₂O₄ polymorphs observed in previous experiments. The spinel phase (stable at ambient conditions) dissociates into periclase and corundum at 14 GPa. With increasing pressure, a phase change from the two oxides to a calcium-ferrite phase occurs, and finally transforms to a calcium-titanate phase at 68 GPa. The calcium-titanate phase is stable up to at least 150 GPa, and we did not observe a stability field for a hexagonal phase or periclase + Rh₂O₃(II)-type Al₂O₃. The bulk moduli of the phases calculated in this study are in good agreement with those measured in high-pressure experiments. Our results differ from those of a previous study using similar methods. We attribute this inconsistency to an incomplete optimization of a cell shape and ionic positions at high pressures in the previous calculations.     

Lattice preferred orientation in CaIrO<Subscript>3</Subscript> perovskite and post-perovskite formed by plastic deformation under pressure

Lattice preferred orientations (LPO) developed in perovskite and post-perovskite structured CaIrO₃ were studied using the radial X-ray diffraction technique combined with a diamond anvil cell. Starting materials of each phase were deformed from 0.1 MPa to 6 GPa at room temperature. Only weak LPO was formed in the perovskite phase, whereas strong LPO was formed in the post-perovskite phase with an alignment of the (010) plane perpendicular to the compression axis. The present result suggests that the (010) is a dominant slip plane in the post-perovskite phase and it is in good agreement with the crystallographic prediction, dislocation observations via transmission electron microscopy, and a recent result of simple shear deformation experiment at 1 GPa–1,173 K. However, the present result contrasts markedly from the results on MgGeO₃ and (Mg,Fe)SiO₃, which suggested that the (100) or (110) is a dominant slip plane with respect to the post-perovskite structure. Therefore it is difficult to discuss the behavior of the post-perovskite phase in the Earth’s deep interior based on existing data of MgGeO₃, (Mg,Fe)SiO₃ and CaIrO₃. The possible sources of the differences between MgGeO₃, (Mg,Fe)SiO₃ and CaIrO₃ are discussed.     

Low-temperature heat capacities of MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> and spinels of the MgCr<Subscript>2</Subscript>O<Subscript>4</Subscript>–MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> solid solution

We present results from low-temperature heat capacity measurements of spinels along the solid solution between MgAl₂O₄ and MgCr₂O₄. The data also include new low-temperature heat capacity measurements for MgAl₂O₄ spinel. Heat capacities were measured between 1.5 and 300 K, and thermochemical functions were derived from the results. No heat capacity anomaly was observed for MgAl₂O₄ spinel; however, we observe a low-temperature heat capacity anomaly for Cr-bearing spinels at temperatures below 15 K. From our data we calculate standard entropies (298.15 K) for Mg(Cr,Al)₂O₄ spinels. We suggest a standard entropy for MgAl₂O₄ of 80.9 ± 0.6 J mol⁻¹ K⁻¹. For the solid solution between MgAl₂O₄ and MgCr₂O₄, we observe a linear increase of the standard entropies from 80.9 J mol⁻¹ K⁻¹ for MgAl₂O₄ to 118.3 J mol⁻¹ K⁻¹ for MgCr₂O₄.     

Low-temperature heat capacities,entropies and enthalpies of Mg<Subscript>2</Subscript>SiO<Subscript>4</Subscript> polymorphs,and α−β−γ and post-spinel phase relations at high pressure

The low-temperature isobaric heat capacities (C _p) of β- and γ-Mg₂SiO₄ were measured at the range of 1.8–304.7 K with a thermal relaxation method using the Physical Property Measurement System. The obtained standard entropies (S°₂₉₈) of β- and γ-Mg₂SiO₄ are 86.4 ± 0.4 and 82.7 ± 0.5 J/mol K, respectively. Enthalpies of transitions among α-, β- and γ-Mg₂SiO₄ were measured by high-temperature drop-solution calorimetry with gas-bubbling technique. The enthalpies of the α−β and β−γ transitions at 298 K (ΔH°₂₉₈) in Mg₂SiO₄ are 27.2 ± 3.6 and 12.9 ± 3.3 kJ/mol, respectively. Calculated α−β and β−γ transition boundaries were generally consistent with those determined by high-pressure experiments within the errors. Combining the measured ΔH°₂₉₈ and ΔS°₂₉₈ with selected data of in situ X-ray diffraction experiments at high pressure, the ΔH°₂₉₈ and ΔS°₂₉₈ of the α−β and β−γ transitions were optimized. Calculation using the optimized data tightly constrained the α−β and β−γ transition boundaries in the P, T space. The slope of α−β transition boundary is 3.1 MPa/K at 13.4 GPa and 1,400 K, and that of β−γ boundary 5.2 MPa/K at 18.7 GPa and 1,600 K. The post-spinel transition boundary of γ-Mg₂SiO₄ to MgSiO₃ perovskite plus MgO was also calculated, using the optimized data on γ-Mg₂SiO₄ and available enthalpy and entropy data on MgSiO₃ perovskite and MgO. The calculated post-spinel boundary with a Clapeyron slope of −2.6 ± 0.2 MPa/K is located at pressure consistent with the 660 km discontinuity, considering the error of the thermodynamic data.     

The nature of blue luminescence from natural benitoite BaTiSi<Subscript>3</Subscript>O<Subscript>9</Subscript>

It is proved that blue luminescence from benitoite is connected with intrinsic luminescence centers, namely isolated TiO₆ octahedra. The metastable level ³T_1u is the emitting level at low temperatures with a long decay time of 1.1 ms. At higher temperatures an energy level with higher radiation probability must be involved in the emission process, and this level is situated at 0.06 eV higher than the lowest level. These two levels may be connected with ³T_1u level splitting or with closely spaced ³T_1u and ³T_2u levels. Decay time shortening and thermal quenching are connected with nonradiative decay within the TiO₆ luminescence center, while energy migration does not take place at least up to room temperature.     

Tilting and distortion of CaSnO<Subscript>3</Subscript> perovskite to 7 GPa determined from single-crystal X-ray diffraction

The structural changes of CaSnO₃, a GdFeO₃-type perovskite, have been investigated to 7 GPa in a diamond-anvil cell at room temperature using single-crystal X-ray diffraction. Significant changes are observed in both the octahedral Sn–O bond lengths and tilt angles between the SnO6 octahedra. The octahedral (SnO6) site shows anisotropic compression and consequently the distortion of SnO6 increases with pressure. Increased pressure also results in a decrease of both of the inter-octahedral angles, Sn–O1–Sn and Sn–O2–Sn, indicating that octahedral tilting increases with increasing pressure, chiefly equivalent to rotation of the SnO6 octahedra about the pseudocubic <001>_p axis. The distortion in the CaO12 and SnO6 sites, along with the octahedral SnO6 tilting, is attributed to the SnO6 site being less compressible than the CaO12 site.Acknowledgments The authors acknowledge with gratitude the financial support for this work from NSF grant EAR-0105864. Ruby pressure measurements were conducted with the Raman system in the Vibrational Spectroscopy Laboratory in the Department of Geosciences at Virginia Tech with the help of Mr. Charles Farley.     

New high-pressure phases in BaCO<Subscript>3</Subscript>

Barium carbonate (BaCO₃) was examined in a diamond anvil cell up to a pressure of 73 GPa using an in situ angle-dispersive X-ray diffraction technique. Three new phases of BaCO₃ were observed at pressures >10 GPa. From 10 to 24 GPa, BaCO₃-IV had a post-aragonite structure with space group Pmmn. There are two molecules in a single unit cell (Z = 2) of the orthorhombic phase, which is same as the high-pressure phases of CaCO₃ and SrCO₃. The isothermal bulk modulus of BaCO₃-IV is K ₀ = 84(4) GPa, with V ₀ = 129.0(7) Å³ when K ₀′ = 4. The c axis of the unit cell parameter is less compressible than the a and b axes. The relative change in volume that accompanies the transformation between BaCO₃-III and BaCO₃-IV is ～6%. BaCO₃-V, which has an orthorhombic symmetry, was synthesized at 50 GPa. As the pressure increases, BaCO₃-V is transformed into tetragonal BaCO₃-VI. This transformation is likely to be second order, because the diffraction pattern of BaCO₃-V is similar to that of BaCO₃-VI, and some single peaks in BaCO₃-VI become doublets in BaCO₃-V. After decompression, the new high-pressure phases transform into BaCO₃-II. Our findings resolve a dispute regarding the stable high-pressure phases of BaCO₃.     

Equation of state and crystal structure of Sb<Subscript>2</Subscript>S<Subscript>3</Subscript> between 0 and 10 GPa

High-precision unit-cell volume data of stibnite, collected in the pressure range of 0–10 GPa, was used for fitting a third-order Birch–Murnaghan equation of state. The zero-pressure volume, bulk modulus and its pressure derivative were found to be 487.73(6) ų, 26.91(14) GPa and 7.9(1), respectively. A series of X-ray intensity data was collected in the same pressure range using a CCD-equipped Bruker diffractometer. The high-pressure structures were all refined to R1(|F₀|>4) values of approximately 0.03. Crystal-chemical parameters as polyhedron volume, centroid and eccentricity were calculated for the seven coordinated cation positions using the software IVTON. The cation eccentricity appears to be a very useful tool for quantification of the lone electron pair activity. U₂S₃, Dy₂S₃ and Nd₂Te₃ are all isostructural with stibnite, but the cations in these materials have no lone electron pair. Their eccentricity is much smaller than that of Sb, and close to zero. This confirms that the stibnite structure type alone does not force eccentricity upon the cations involved and it is the lone electron pairs of Sb that generate the eccentricity of cation positions in the structures of stibnite. At increasing pressure the eccentricity of Sb is decreasing. It is therefore reasonable to conclude that the lone electron pair activity is decreasing with increasing pressure.     

Equation of state,structural behaviour and phase diagram of synthetic MgFe<Subscript>2</Subscript>O<Subscript>4</Subscript>, as a function of pressure and temperature

The behaviour of synthetic Mg-ferrite (MgFe₂O₄) has been investigated at high pressure (in situ high-pressure synchrotron radiation powder diffraction at ESRF) and at high temperature (in situ high-temperature X-ray powder diffraction) conditions. The elastic properties determined by the third-order Birch–Murnaghan equation of state result in K₀=181.5(± 1.3) GPa, K=6.32(± 0.14) and K= –0.0638 GPa^–1. The symmetry-independent coordinate of oxygen does not show significant sensitivity to pressure, and the structure shrinking is mainly attributable to the shortening of the cell edge (homogeneous strain). The lattice parameter thermal expansion is described by _a0+_a1*(T–298)+_a2/(T–298)², where _a0=9.1(1) 10^–6 K^–1, _a1=4.9(2) 10^–9 K^–2 and _a2= 5.1(5) 10^–2 K. The high-temperature cation-ordering reaction which MgFe-spinel undergoes has been interpreted by the ONeill model, whose parameters are = 22.2(± 1.8) kJ mol^–1 and =–17.6(± 1.2) kJ mol^–1. The elastic and thermal properties measured have then been used to model the phase diagram of MgFe₂O₄, which shows that the high-pressure transition from spinel to orthorombic CaMn₂O₄-like structure at T < 1700 K is preceded by a decomposition into MgO and Fe₂O₃.     

Fine structure in photoluminescence spectrum of S<Subscript>2</Subscript><Superscript>−</Superscript> center in sodalite

The photoluminescence and excitation spectra of sodalites from Greenland, Canada and Xinjiang (China) are observed at 300 and 10 K in detail. The features of the emission and excitation spectra of the orange-yellow fluorescence of these sodalites are independent of the locality. The emission spectra at 300 and 10 K consist of a broad band with a series of peaks and a maximum peak at 648 and 645.9 nm, respectively. The excitation spectra obtained by monitoring the orange-yellow fluorescence at 300 and 10 K consist of a main band with a peak at 392 nm. The luminescence efficiency of the heat-treated sodalite from Xinjiang is about seven times as high as that of untreated natural sodalite. The emission spectrum of the S₂ ⁻ center in sodalite at 10 K consists of a band with a clearly resolved structure with a series of maxima spaced about 560 cm⁻¹ (20–25 nm) apart. Each narrow band at 10 K shows a fine structure consisting of a small peak due to the stretching vibration of the isotopic species of ³²S³⁴S⁻, a main peak due to that of the isotopic species of ³²S₂ ⁻ and five peaks due to phonon sidebands of the main peak.     

Minerals at high pressure. Mechanics of compression from quantum mechanical calculations in a case study: the beryl (Al<Subscript>4</Subscript>Be<Subscript>6</Subscript>Si<Subscript>12</Subscript>O<Subscript>36</Subscript>)

The Bader topological analysis has been applied to ab initio computed electron densities of beryl, in order to clarify its mechanism of compression. Full structural optimization and total energy (E) calculations were performed at different cell volumes (V _c). The pressure at each volume and the equation of state were estimated from the first and second derivatives of the resultant E(V _c) curve. The total (negative) potential energy of the crystal, sum of both attractive and repulsive electrostatic terms, was found to systematically decrease (i.e., it moved to more negative values) up to the highest pressure considered (28.4 GPa), indicating that interelectronic and internuclear repulsions are not the only terms controlling the compressibility, at least in the pressure range investigated. Electronic kinetic energy increases as the cell volume is reduced, leading to a parallel increase of the total energy. Both structure at equilibrium and compressibility are therefore due to the balance between the opposing kinetic and potential energy terms. The Bader theory has been used to identify the topological atoms within the structure and to calculate their properties, with particular attention to the forces driving the structural relaxation at high pressure. On a qualitative basis, the obtained results are expected to be transferable to the discussion of compressibility of other mineral phases.     

Ab initio study of the high-pressure behavior of CaSiO3 perovskite

Using density functional simulations, within the generalized gradient approximation and projector-augmented wave method, we study structures and energetics of CaSiO₃ perovskite in the pressure range of the Earths lower mantle (0–150 GPa). At zero Kelvin temperature the cubic CaSiO₃ perovskite structure is unstable in the whole pressure range, at low pressures the orthorhombic (Pnam) structure is preferred. At 14.2 GPa there is a phase transition to the tetragonal (I4/mcm) phase. The CaIrO₃-type structure is not stable for CaSiO₃. Our results also rule out the possibility of decomposition into oxides.

First-principles-based calculations of the CaCO<Subscript>3</Subscript>–MgCO<Subscript>3</Subscript> and CdCO<Subscript>3</Subscript>–MgCO<Subscript>3</Subscript> subsolidus phase diagrams

 Planewave pseudopotential calculations of supercell total energies were used as bases for first-principles calculations of the CaCO₃–MgCO₃ and CdCO₃–MgCO₃ phase diagrams. Calculated phase diagrams are in qualitative to semiquantitative agreement with experiment. Two unobserved phases, Cd₃Mg (CO₃)₄ and CdMg₃(CO₃)₄, are predicted. No new phases are predicted in the CaCO₃–MgCO₃ system, but a low-lying metastable Ca₃Mg(CO₃)₄ state, analogous to the Cd₃Mg(CO₃)₄ phase is predicted. All of the predicted lowest-lying metastable states, except for huntite CaMg₃(CO₃)₄, have dolomite-related structures, i.e. they are layer structures in which A _m B _n cation layers lie perpendicular to the rhombohedral [111] vector. Received: 6 May 2002 / Accepted: 23 October 2002 Acknowledgements This work was partially supported by NSF contract DMR-0080766 and NIST.     

Stability at high-pressure,elastic behaviour and pressure-induced structural evolution of CsAlSi<Subscript>5</Subscript>O<Subscript>12</Subscript>, a potential host for nuclear waste

The elastic and structural behaviour of the synthetic zeolite CsAlSi₅O₁₂ (a = 16.753(4), b = 13.797(3) and c = 5.0235(17) Å, space group Ama2, Z = 2) were investigated up to 8.5 GPa by in situ single-crystal X-ray diffraction with a diamond anvil cell under hydrostatic conditions. No phase-transition occurs within the P-range investigated. Fitting the volume data with a third-order Birch–Murnaghan equation-of-state gives: V ₀ = 1,155(4) ų, K _T0 = 20(1) GPa and K′ = 6.5(7). The “axial moduli” were calculated with a third-order “linearized” BM-EoS, substituting the cube of the individual lattice parameter (a ³, b ³, c ³) for the volume. The refined axial-EoS parameters are: a ₀ = 16.701(44) Å, K _T0a = 14(2) GPa (β_a = 0.024(3) GPa^?1), K′ _a = 6.2(8) for the a-axis; b ₀ = 13.778(20) Å, K _T0b = 21(3) GPa (β_b = 0.016(2) GPa^?1), K′ _b = 10(2) for the b-axis; c ₀ = 5.018(7) Å, K _T0c = 33(3) GPa (β_c = 0.010(1) GPa^?1), K′ _c = 3.2(8) for the c-axis (K _T0a:K _T0b:K _T0c = 1:1.50:2.36). The HP-crystal structure evolution was studied on the basis of several structural refinements at different pressures: 0.0001 GPa (with crystal in DAC without any pressure medium), 1.58(3), 1.75(4), 1.94(6), 3.25(4), 4.69(5), 7.36(6), 8.45(5) and 0.0001 GPa (after decompression). The main deformation mechanisms at high-pressure are basically driven by tetrahedral tilting, the tetrahedra behaving as rigid-units. A change in the compressional mechanisms was observed at P ≤ 2 GPa. The P-induced structural rearrangement up to 8.5 GPa is completely reversible. The high thermo-elastic stability of CsAlSi₅O_12, the immobility of Cs at HT/HP-conditions, the preservation of crystallinity at least up to 8.5 GPa and 1,000°C in elastic regime and the extremely low leaching rate of Cs from CsAlSi₅O₁₂ allow to consider this open-framework silicate as functional material potentially usable for fixation and deposition of Cs radioisotopes.     

Effects of nitrogen on the ecosystem respiration,CH<Subscript>4</Subscript> and N<Subscript>2</Subscript>O emissions to the atmosphere from the freshwater marshes in northeast China

Freshwater marshes could be a source of greenhouse gases emission because they contain large amounts of soil carbon and nitrogen. These emissions are strongly influenced by exogenous nitrogen. We investigate the effects of exogenous nitrogen on ecosystem respiration (CO₂), CH₄ and N₂O emissions from freshwater marshes in situ in the Sanjiang Plain Northeast of China during the growing seasons of 2004 and 2005, using a field fertilizer experiment and the static opaque chamber/GC techniques. The results show that there were no significant differences in patterns of seasonal variations of CO₂ and CH₄ among the fertilizer and non-fertilizer treatments, but the seasonal patterns of N₂O emission were significantly influenced by the exogenous nitrogen. Seasonal averages of the CO₂ flux from non-fertilizer and fertilizer were 987.74 and 1,344.35 mg m ⁻² h ⁻¹, respectively, in 2004, and 898.59 and 2,154.17 mg m ⁻² h ⁻¹, respectively, in 2005. And the CH₄ from the control and fertilizer treatments were 6.05 and 13.56 mg m ⁻² h ⁻¹ and 0.72 and 1.88 mg m ⁻² h ⁻¹, respectively, in 2004 and 2005. The difference of N₂O flux between the fertilizer and non-fertilizer treatments is also significant either in 2004 and 2005. On the time scale of 20-, 100-, and 500-year periods, the integrated global warming potential (GWP) of CO₂ + CH₄ + N₂O released during the two growing seasons for the treatment of fertilizer was 97, 94 and 89%, respectively, higher than that for the control, which suggested that the nitrogen fertilizer can enhance the GWP of the CH₄ and N₂O either in long time or short time scale.     

Hydrostatic compression of galenobismutite (PbBi<Subscript>2</Subscript>S<Subscript>4</Subscript>): elastic properties and high-pressure crystal chemistry

A single-crystal sample of galenobismutite was subjected to hydrostatic pressures in the range of 0.0001 and 9 GPa at room temperature using the diamond-anvil cell technique. A series of X-ray diffraction intensities were collected at ten distinct pressures using a CCD equipped 4-circle diffractometer. The crystal structure was refined to R1(|F₀| > 4σ) values of approximately 0.05 at all pressures. By fitting a third-order Birch-Murnaghan equation of state to the unit-cell volumes V ₀ = 700.6(2) ų, K ₀ = 43.9(7) GPa and dK/dP = 6.9(3) could be determined for the lattice compression. Both types of cations in galenobismutite have stereochemically active lone electron pairs, which distort the cation polyhedra at room pressure. The cation eccentricities decrease at higher pressure but are still pronounced at 9 GPa. Galenobismutite is isotypic with CaFe₂O₄ (CF) but moves away from the idealised CF-type structure during compression. Instead of the two octahedral cation sites and one bi-capped trigonal-prismatic site, PbBi₂S₄ attains a new high-pressure structure characterised by one octahedral site and two mono-capped trigonal-prismatic sites. Analyses of the crystal structure at high pressure confirm the preference of Bi for the octahedral site and the smaller one of the two trigonal-prismatic sites.