The molecular dynamics method was applied to aluminous perovskite solid solution in the system MgSiO
3-Al
2O
3. The two-body interatomic potential model CMAS94 (Matsui, 1994) was used for this purpose. The solid solution crystals (Mg
1−x,Al
x)(Al
x,Si
1−x)O
3 (0≤
x≤1) accommodating Al
3+ ions without the production of vacant site were realized by the following two methods:
(1) To replace Mg
2+ and Si
4+ ions with 2Al
3+ ions at random disregarding the arrangement of Mg
2+ and Si
4+ in the crystal (Random substitution).
(2) To replace a pair of (Mg
2++Si
4+) with 2Al
3+ ions at the adjacent sites, so as to satisfy the local charge balance.
Although the latter solid solution model gives slightly lower values of enthalpy, both models show virtually the same energetic and structural characteristics. Under atmospheric condition, both the unit-cell parameter
c and the cell volume increase remarkably with the increase in the Al
2O
3 content. At around 40 - 60 GPa, however, the volume is almost insensitive to the Al
2O
3 content. Al
2O
3-containing perovskite is more compressible than aluminum-free perovskite. It is consistent with the results of recent high-pressure
in-situ X-ray studies below 15 GPa. This compositional effect on compressibility becomes much smaller with increase in pressure.
A positive excess enthalpy for the MgSiO
3-Al
2O
3 binary system is observed.
View full abstract