Structures and Transport Properties of CaCO3 Melts under Earth's

Subscriber access provided by READING UNIV

Article 3

Structures and transport properties of CaCO melts under Earth’s mantle conditions Xiangpo Du, Min Wu, John S Tse, and Yuanming Pan ACS Earth Space Chem., Just Accepted Manuscript • DOI: 10.1021/ acsearthspacechem.7b00100 • Publication Date (Web): 21 Nov 2017 Downloaded from http://pubs.acs.org on November 23, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Earth and Space Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

1

TOC:

2 3

ABSTRACT

4

Carbonatitic and carbonated silicate melts generated from melting in the mantle are

5

the chief agents for liberating carbon from the solid Earth and exert important controls

6

on Earth’s deep carbon cycle. However, significant gaps in our knowledge about the

7

carbonatitic melts under conditions pertinent to Earth’s deep interior remain due to

8

experimental challenges. Here we report on a first-principles molecular dynamics

9

(FPMD) calculation of calcium carbonate (CaCO3) melts at pressures up to 52.5 GPa.

10

Our FPMD calculations reproduce the ultralow viscosity measured by experiments

11

and confirm the ideal liquid behavior of calcium carbonate melts at pressures below

12

11.2 GPa. However, calcium carbonate melts are characterized by a pronounced

13

nonideal liquid behavior at pressures above 11.2 GPa, arising from 1) the temporal

14

formation of small carbonate clusters and 2) increased interactions between Ca2+ and

15

CO32- ions. It is found that the Stokes-Einstein equation relating the viscosity with the

16

diffusion coefficients still holds at high pressure provided that a suitable effective

17

particle size can be chosen. 2

ACS Paragon Plus Environment


1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

1. Introduction

2

Calcium carbonate is the main component of skeletons of living organisms and

3

also occurs as the dominant minerals in limestones and marbles as well as a common

4

phase in many rocks, constituting an enormous carbon reservoir.1 At ambient

5

conditions, the two common crystalline phases of calcium carbonate are calcite and

6

aragonite. Aragonite is metastable with respect to calcite and is the stable form above

7

~650 K.2 Under ambient pressure, calcite decomposes into CaO and CO2 above 1170

8

K.3,4 High pressure experiments in laser-heated diamond anvil cells (DAC) have

9

shown that calcium carbonate melts congruently at least up to 3400 K when the

10

pressure is lower than 43 GPa. Moreover, compressed calcium carbonate melts

11

decompose into CaO and CO2 at temperatures higher than 3900 K.4 In addition, the

12

CO2 produced from the decomposition of carbonate melts disassociates further into

13

molecular oxygen and perhaps to form diamond at 30-80 GPa.5,6 Carbonated melting

14

that generates carbonatitic and carbonated silicate melts in the mantle has been

15

proposed to be the chief agents for liberating carbon from the solid Earth and thus

16

exerting important controls on the Earth’s deep carbon cycle.7

17

Multiple lines of evidence from inclusions in superdeep diamond, high P-T

18

melting experiments and theoretical calculations suggest that carbonatitic melts are

19

present down to the lower mantle.8-10 Therefore, knowledge about the structures and

20

properties of carbonatitic melts under pressure not only has significant impact on the

21

geological processes research from mantle geodynamics to volcanic activities, but

22

also important for the understanding of the related deep carbon cycle and global 3


Page 2 of 29

Page 3 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1

climate change.11-15

2

Compared to the much studied silicate melts, the physical properties of carbonate

3

melts, in particular under extreme conditions (high pressure and high temperature),

4

are scarcely investigated.12,15-18 Recently, an ultralow viscosity of carbonate melts (i.e.,

5

about one order of magnitude less than calcium silicate melts) up to 6.2 GPa has been

6

measured experimentally.19 The higher viscosity of silicate melts is attributable to the

7

polymerized SiO4 units that impede diffusion. The observed ultralow viscosity of

8

calcium carbonate melts was supported by a molecular dynamics (MD) calculation

9

using an empirical potential fitted to first-principles calculations.17 However, due to

10

the difficulty of in-situ high temperature and high pressure experiments, the viscosity

11

of calcium carbonate melts under conditions relevant to Earth’s transition zone and

12

lower mantle has not been studied. To better understand the transport properties of

13

calcium carbonate melts in Earth’s deep interior, experimental investigations at higher

14

pressures are needed but remain technically challenging. In this paper, we extended

15

the investigation of the structures and transport properties of calcium carbonate melts

16

to Earth’s lower mantle conditions.

17

It is well known that high pressure can alter the electronic structures of a material

18

dramatically, and consequently, changes its geometrical structures and related physical

19

properties. For instance, elemental calcium is known to transform from face centered

20

cubic (FCC) to body centered cubic (BCC) structures at around 8 GPa.20 Theoretical

21

work has also shown that orbital hybridization between Ca and Al starts to emerge at

22

~12 GPa in compressed Ca-Al metallic glasses.21 Therefore, it is often not possible to 4



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

extrapolate the property information of minerals and melts obtained at low pressures

2

to higher pressures. To circumvent the uncertainty of MD calculations employing

3

force fields parameterized at ambient conditions, first-principles molecular dynamic

4

(FPMD) calculation has been proven to be able to correctly describe the interatomic

5

interaction at high temperature and high pressure. The goals of the contribution are to

6

unravel the mechanisms for the viscosity evolution and accompanied structural

7

changes of calcium carbonate melts at high pressures.

8

2. Methods

9

FPMD calculations and electronic structure calculations were performed using the

10

Vienna Ab initio Simulation Package (VASP) based on the density functional

11

theory.22,23 Effects of the core electrons were replaced by atomic projected augmented

12

wave (PAW) potentials.24 Valence electronic configurations of 2s22p2 and 2s22p4 were

13

used for C and O atomic potentials, respectively. To better describe the electronic

14

structure of Ca atoms under pressure, a hard Ca potential with valence electronic

15

configurations of 3s23p64s2 (Ca_sv) was used in the present study. Calculation using a

16

softer Ca potential neglecting the 3s orbital (Ca_pv) was also performed for

17

comparison. The Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional was

18

used.25 The wave functions were represented by a plane-wave basis set with an energy

19

cutoff of 400 eV. A molten structure of CaCO3 at ambient pressure was generated by

20

melting a 2×2×1 supercell constructed from the calcite unit cell containing 120 atoms

21

(24 CaCO3 formula units) by heating to 3000 K. The liquid structure at ambient

22

pressure and 2000K was obtained by cooling slowly from 3000 K in several 5


Page 4 of 29

Page 5 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1

temperature steps. Higher pressure liquid structures were obtained by isotropically

2

reducing the volume of the liquid model obtained at ambient pressure. The pressure

3

was then determined from the stress tensor. Shear viscosities were computed from the

4

Green-Kubo auto-correlation function of the off-diagonal components of the stress

5

tensor (SACF) in the canonical ensemble (NVT).26-28 The Nose-Hoover thermostat

6

was used to control the temperature27. A single  point was used for Brillouin zone

7

sampling. A very long trajectory is required to achieve convergence in the calculated

8

viscosity. Therefore, at each pressure, MD calculation ran for at least 100 ps. Since

9

viscosity is usually difficult to measure at high pressure. Most often the experimental

10

viscosity was estimated from the use of the Stokes-Einstein relationship from the

11

diffusion constant, which is easier to determine using x-ray tomography.29,30 This

12

simple relationship, however, has not been validated theoretically for high pressure.

13

Therefore, we compare the viscosities calculated from the Green-Kubo correlation

14

function31 with those estimated from the Stokes-Einstein equation from the

15

corresponding diffusion coefficient computed from the same MD trajectory.

16

Table 1. Comparison of the densities (g/cm3) of calcite and aragonite crystalline

17

structures by first-principles calculations and experiments. Exp.

Ca_sv

Ca_pv

BLYP+D2

B3LYP

PBE

PBE

2.72 a

2.75c

2.78c

2.67d

2.604e

2.603g

2.650i

Aragonite 2.93b

3.00c

3.04c

2.83d

2.821f

2.851h

Calcite

18 19 20

a

Ref. 32 Ref. 33

b

6



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1 2 3 4 5 6 7 8

c

This work. Ref. 17 e Ref. 34 f Ref. 35 g Ref. 36 h Ref. 37 i Ref. 38 d

9 10

Figure 1, Structure of calcium carbonate polymorphs. (a) Calcite; (b) A snapshot of

11

molten calcium carbonate at 2000 K and 0.6 GPa. The grey, brown and red balls

12

represent Ca, C and O atoms, respectively.

13 14

3. Results and Discussion

15

3.1 Equation of states of crystalline and molten calcium carbonates

16

In calcite (Fig.1a), triangular CO3 anions formed from sp2 hybridized C with three

17

oxygen atoms were arranged to form planes where a layer of Ca atoms were

18

intercalated between these planes. To evaluate the accuracy of the parameters used in

19

the first-principles calculations, the theoretical densities of calcite are compared to

20

those from experiments and previous calculations in Table 1. For both calcite and

21

aragonite the present theoretical calculations gave the best densities as compared to

22

the experimental values. The density of calcite calculated with the Ca_sv potential at 7


Page 6 of 29

Page 7 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1

ambient pressure is 2.75 g/cm3 and agrees well with the experimental observation of

2

2.72 g/cm3.32 In comparison, the calculated density of 2.78 g/cm3 from the Ca_pv

3

potential is slightly overestimated. The density of aragonite calculated with the Ca_sv

4

potential of 3.00 g/cm3 agrees well with the experiment observation of 2.93 g/cm3 and,

5

once again, the density of 3.04 g/cm3 obtained from the Ca_pv potential is slightly

6

overestimated.33 These results show the importance of including the inner-core 3s2

7

electrons in the calculations even at ambient pressure. It is noteworthy that all

8

previous theoretical calculations employing different GGA functionals underestimated

9

the densities of both calcite and aragonite (Table 1). For instance, the densities of

10

calcite and aragonite predicted by Vuilleumier et al.17 using GGA (BLYP) together

11

with Van der Waal dispersion at 0 K are ~2% smaller than those from experimental

12

results at ambient condition and ~3% smaller than our calculated result using the

13

Ca_sv potential at 0 K, respectively. The results demonstrate that the structures

14

predicted by the present choice of computational parameters are more reliable than

15

earlier works.

16

From Bader atomic charge analysis,39 each cationic Ca loses 1.7 electrons to the

17

CO3 anion. The electronic density of states shows calcite is an insulator with an

18

energy gap of 5 eV. The highest occupied crystal orbitals (HOCO) were composed

19

mainly of O lone pair electrons. This description is similar to that from previous

20

theoretical calculation.40 In the melts (Fig. 1b), the arrangement of CO3 group is more

21

disordered and the bandgap energy decreases dramatically to 1.2 eV. There is no

22

spectroscopic study on the of CaCO3 melt. The evaluation on the accuracy of the 8



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 29

1

computed density of molten calcium carbonate at high temperature and high pressure

2

is more difficult due to lack of experimental data. A measurement on the CaCO3–

3

Li2CO3–Na2CO3–K2CO3 quaternary system by Liu et al.41 suggested the density of

4

CaCO3 melts to be ~2.6 g/cm3 at 1073 K and 0.7 GPa. From thermodynamic data,

5

Treiman et al.42 estimated a density of 2.4 g/cm3 at 1596 K and 0.1 GPa. These values

6

are close to the present calculated density of ~2.52 g/cm3 at 2000 K and 0.7 GPa. On

7

the other hand, the FPMD calculation by Vuilleumier et al.

8

2.94 g/cm3 at 2073 K and 12.2 GPa, which is close to the present value of ~3.1 g/cm3

9

at 12 GPa and 2000 K. The agreement is even more encouraging taking into the

10

consideration that the densities of calcite and aragonite were underestimated from

11

previous theoretical studies as discussed above (vide infra). Moreover, the different

12

exchange-correlation functionals used in the calculations and the procedures to

13

generate the liquid structure and the length of the MD run may also contribute to the

14

discrepancies. We conclude that the calculated densities of the melts in the present

15

study are reliable.

16

3.2 Viscosity of molten calcium carbonate under pressure

17

3.2.1 Viscosity calculated from Green-Kubo relation

17

predicted a density of

18

The goals of this study are to investigate the CaCO3 melts and to test the validity of

19

the Stokes-Einstein relationship29,30 at high pressure. The viscosities can be calculated

20

from first-principles from the time correlation function of the temporal stress tensors

21

from the MD calculation via the Green-Kubo relation,31

22

=𝑘

𝑉 𝐵

∞

∫ 〈𝑃𝑥𝑦 (𝑡)𝑃𝑥𝑦 (0)〉𝑑𝑡 𝑇 0 9


(1)

Page 9 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1

where η is the viscosity, V is the volume of the system, T is the temperature, kB is

2

the Boltzmann constant,〈〉denotes the ensemble average, and Pxy is the off-diagonal

3

component of the instantaneous stress tensor.

4 5

Figure 2, (a) Pressure dependent viscosities of molten calcium carbonates at 2000 K

6

and 3000 K. The experiment result at 2063 K and 6 GPa was marked as a star. The

7

solid red and black circles represent the calculated viscosity from the Green-Kubo

8

relationship at 2000 K and 3000 K, respectively. The open pink and blue circles

9

represent the calculated viscosity from the Stokes-Einstein (SE) equation using the 10



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

theoretical diffusion coefficients (see text) at 2000 K and 3000K, respectively. The

2

dashed lines are guides for eyes. (b) The SACF (black solid line) and the integral

3

(blue solid line) giving the viscosity calculated at 3000 K and 32.6 GPa. The pink

4

dashed line is a guide for eyes.

5 6

In order to maintain comparable viscosity at high pressure in the experiment, it

7

was necessary to increase the temperature. For example, the experiment of Kono et

8

al.19 showed that a measured viscosity of ~6 mPa·s at 0.9 GPa and 1380 oC (1653 K)

9

required an increase of temperature to 1790 OC (2063 K) at 6.2 GPa. Compared to the

10

experiment, the calculated viscosity at 0.7 GPa and 2000 K is ~3.8 mPa·s. To

11

maintain a similar viscosity at a higher pressure of 7 GPa, it was necessary to increase

12

the temperature to ~3000 K. The predicted trend is consistent with experimental

13

observations although the absolute temperature difference is higher. The calculated

14

pressure dependent viscosities at 2000 K and 3000 K are shown in Fig. 2. To illustrate

15

the convergence of the viscosity calculation, the SACF and the intergral (viscosity) as

16

a function of the simulation time at 32.6 GPa were presented in Fig. 2b. The

17

calculated viscosity at 4.1 GPa and 2000 K of 5.8 mPa·s agrees with the experimental

18

value of ~6 mPa·s measured at 2063 K and 6.2 GPa.19 However, at 2000 K and above

19

8 GPa, calcium carbonate starts to solidify and becomes a viscous liquid with very

20

low self-diffusion coefficients and very long simulation time is needed. Therefore, we

21

can only predict reliably the viscosity of the melt below this pressure. Since the

22

experimental melting temperature of CaCO3 is close to 2000 K, the very slow 11


Page 10 of 29

Page 11 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1

relaxation of the CaCO3 at 2000 K is indicative the melting temperature has increased

2

with pressure. In comparison, the calcium carbonate prepared at 3000 K is more fluid

3

like and accurate viscosities can be computed up to 32.6 GPa. The viscosity at 3000 K

4

and 0.6 GPa is 2 mPa-s and incidentally, close to the water viscosity of 0.89 mPa·s at

5

ambient conditions. The result confirms the experimental finding. Furthermore, the

6

calculated viscosity at 3000 K increases slowly at low pressure ( 11.2 GPa) and may

10

signal a structural change (vide supra).

11 12

Figure 3, Variation of calculated diffusion coefficient D with pressure at 2000 K and

13

3000K.

14 15 16

3.2.2 Self-diffusion coefficient The self-diffusion coefficients D for the ions in calcium carbonate melts were 12



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

Page 12 of 29

calculated from the Einstein mean square displacement formula, 1

〈𝑟𝑖 (𝑡)−𝑟𝑖 (0)〉2

𝑁

6𝑡

𝐷(𝑡, 𝑡 → ∞) = ∑𝑁 𝑖=1

2

(2)

3

where ri is the position of the ion i and the bracket represents the ensemble average

4

over many time origins. The results at 2000 K and 3000 K are shown in Fig. 3a and

5

3b with the estimated errors. The diffusion coefficients of Ca2+ and CO32- are very

6

similar. The results suggest both Ca2+ and CO32- diffused in unison. The diffusion

7

coefficients doubled from ~0 GPa when the temperature was raised from 2000 K to

8

3000 K. However, at both temperatures, the diffusion coefficients decrease to ca.

9

0.1×108 m2/s at high pressure (8 GPa at 2000 K and 25 GPa at 3000 K). We have

10

followed the suggestion from a previous study17 that the effective diffusion

11

coefficients D of the system were computed from the arithmetic mean of DCa and

12

DCO3 as the CO3 diffuse as a unit and no decomposition was observed in the P-T range

13

studied here. It should be noted that an alternative relationship Deff = [DCa ×

14

DCO3]/[0.5DCa+0.5DCO3] has been suggested to be suitable for mixed charged species.

15

In this case, since DCa is close to DCO3, therefore both expressions gave similar

16

effective diffusion coefficients. As shown in Fig.3, the calculated diffusion coefficient

17

of the CO32- anion of ~10-9 m2/s at 2000 K at low pressure is similar to previous

18

theoretical studies.17,

19

significantly upon compression at low pressure, indicating that at higher pressure (>4

20

GPa at 2000 K, >15 GPa at 3000K), the ions are much more restrictive with

21

compression.

22

3.2.3 Viscosity calculated from the Stokes-Einstein equation

43

The D values at both 2000 K and 3000 K decrease

13


Page 13 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1

In high pressure experimental studies, the liquid viscosity is often estimated from

2

the diffusion coefficient using the Stokes-Einstein equation ηSE=(kBT)/(2πDd), where

3

ηSE is the viscosity, D is the self-diffusion coefficient, d is the particle diameter, T is

4

the temperature and kB is the Boltzman constant. However, the Stokes-Einstein

5

equation is valid only in an ideal fluid composed of spherical particles that have no

6

slippage against the layer of the fluid in contact with it.29 It is well-known that the

7

Stokes-Einstein equation has to be corrected for fluids with non-spherical particles

8

because the friction coefficient is expected to be different. The viscosity calculated

9

from the Stokes-Einstein equation also depends on how the “particle” diameter is

10

determined. In the calcium carbonate melts, the diffusing particles are charged Ca2+

11

and CO32- ions with different shapes and sizes. The condition is very different from an

12

ideal fluid. For instance, the Stokes–Einstein relation has been shown to fail in

13

describing the viscosity of silicate liquids.44 Therefore, it is a somewhat surprising

14

that the Stokes-Einstein relation seems to work quite well for the calcium carbonate

15

melts below 15 GPa in the study of Vuilleumier et al..17

14



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1 2

Figure 4, (a) Calculated pair distribution functions at 3000 K, 0.7 GPa. (b) Effective

3

particle diameters at 2000 K (red) and 3000 K (black) calculated from the equation

4

(dCC + 2dCCa + dCaCa)/4;

5 6

To examine the Stokes-Einstein relation in a broader pressure range, we

7

employed the same procedure described earlier to determine the particle diameters:

8

viz. the effective particle diameter d was calculated from the equation (dCC + 2dCCa +

9

dCaCa)/4, where dCC, dCCa and dCaCa are the first maximum of the respective pair

10

distribution functions (PDFs) obtained from MD calculations. For this case, the

11

results at 0.7 GPa at 3000 K are shown in Fig. 4a. Instead of using a constant particle

12

diameter for all the viscosity calculations as assumed in the previous work, we 15


Page 14 of 29

Page 15 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1

employed the pressure dependent particle diameters derived from the respective pair

2

distribution functions calculated from the MD trajectories. The results are shown in

3

Fig.4b.

4

A comparison of the estimated ηSE from the Stokes-Einstein equation and ηGK

5

computed directly from the Green-Kubo relationship is shown in Fig. 2. The

6

calculated ηSE values at 2000 K agree well with ηGK up to 7 GPa. This observation is

7

consistent with a previous theoretical study.17 Intuitively, it is not expected that the ηSE

8

values at 3000 K would also agree well with ηGK at pressures up to 32.6 GPa, even

9

though there is an obvious abrupt change in the viscosity at ~11.2 GPa which may

10

indicate a structural transformation (vide supra) has occurred. The present results

11

show that the Stokes-Einstein relationship is valid for the calcium carbonate melts

12

within the pressure range studied here as long as an appropriate method is used to

13

estimate the particle size.

14

In retrospect, the success of the Stokes-Einstein equation with an appropriate

15

variable “particle size” at different pressures is not too surprisingly. The FPMD

16

calculations show that below 11.2 GPa, the structure of the carbonate melts, in

17

particular at higher temperature (3000 K), does not changed fundamentally. As will be

18

discussed in greater detail below, at 3000 K and below 11.2 GPa, there is no

19

noticeable change in the C=O distances of the carbonate anions indicating a more or

20

less rigid group. Furthermore, the diffusions of the Ca2+ and CO32- ions are almost

21

independent of each other. Therefore, the assumption of a constant particle size is

22

appropriate in this pressure region. Above 11.2 GPa, an abrupt increase in viscosity is 16



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

predicted. As will be shown below, the higher viscosity is due to increase in the

2

correlation of the Ca2+ and CO32- diffusive motions and appearance of instantaneous

3

clustering of CO32-. In this case, pressure-dependent corrections to the “particle” size

4

are needed. We must caution that the Stokes-Einstein equation should not be regarded

5

a universal relationship.

6 7

Figure 5, Comparison of the pressure dependent pair distribution function (PDF)

8

computed at 2000 K and 3000 K.

9

17


Page 16 of 29

Page 17 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1 2

Figure 6, Pressure dependent partial coordination numbers of atoms in calcium

3

carbonate melts at 3000 K. (a) number of O atoms surrounding a selected C atom, (b)

4

number of C atoms surrounding a selected O atom.

5

6

3.3 Structures of the CaCO3 melts

7

A comparison of the pair distribution functions (PDF) for 2000 K and 3000 K

8

MD calculations are shown in Fig. 5. The Ca PDF at 2000 K and 0.7 GPa shows three

9

principle features at 2.5, 4.2 and 6.6 Å. In comparison, the first peak of the Ca PDF at

10

3000 K is much broader indicating a more disordered nearest environment at higher

11

temperature. The first peak of C PDF and O PDF at ~1.3 Å can be assigned

12

unambiguously to the C-O covalent bonds of the CO3 units.32,33 Interestingly, the 18



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

positions of the peak are not sensitive to pressure, indicating that the CO32- anions are

2

rigid because of the strong C-O covalent bonds. Compression only affects the

3

distributions of the second or third nearest neighbors. A shoulder at ~2.7 Å starts to

4

emerge at high pressure in the O PDF for the structures at both 2000 K and 3000 K.

5

Inspection of the structure revealed that the distance between neighboring CO3 units

6

decreases with increasing pressure. Upon closer examination, it is found that although

7

the first peak of the O PDF is still separated from the second peak at 2000 K, the two

8

peaks start to merge at higher pressures at 3000 K. At 2000 K and 40 GPa, the

9

calculated coordination numbers of carbon and oxygen are 3 and 1, respectively. At

10

3000 K, the coordination numbers of both carbon and oxygen hardly change at

11

pressure below 11.2 GPa (Fig. 6). However, above 11.2 GPa at 3000 K, the

12

coordination numbers of both carbon and oxygen atoms increase steadily. The

13

increase in the CNs is best illustrated by the snapshot on the formation of

14

instantaneous clusters at 3000 K and 32.6 GPa (Fig. 7), although there is no extended

15

network structure connected by the C-O covalent bonds, a few clusters connected

16

instantaneously to two or even three CO3 units are formed to give rise to a few 4

17

coordinated carbons. The observation of these temporal poly-carbonate clusters at

18

high pressures support the previous report of Corradini et al..18 The change from the

19

isolated and independent CO3 to correlated instantaneous clusters is responsible for

20

the abrupt increase in the viscosity above 11.2 GPa.

19


Page 18 of 29

Page 19 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1 2

Figure 7, Instantaneous formation of small CO3 clusters in the calcium carbonate

3

melt at 3000 K and 32.6 GPa.

4

5 6

Figure 8, Pressure dependent velocity cross correlation functions of (a) Ca-C and (b)

7

Ca-O in calcium carbonate melts, together with (c) the velocity density of states at

8

3000 K.

9

20



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

3.4 Ca and CO3 interactions

2

Increasing interaction between Ca2+ and CO32- ions at high pressures also has an

3

important impact on the viscosity. This interaction, as a function of pressure, can be

4

analyzed from the corresponding velocity cross-correlation between Ca and O and

5

between Ca and C atoms. The velocity cross-correlation functions (VCC) at 3000 K

6

are depicted in Fig. 8. It can be seen that the VCC profiles are liquid like at pressures

7

below 7.1 GPa. Weak coupling between the Ca and the CO3 motions is reflected in

8

the lack of variation in the viscosity of liquid CaCO3 in this pressure range. Above

9

11.2 GPa, however, the interactions between Ca and C as well as Ca with O increase

10

noticeably. This is related to the temporal formation of small clusters as described

11

above (vide infra). A band attributed to the Ca…CO3 vibrations starts to emerge in the

12

plot of the Ca projected density of states (Fig. 8c). This vibrational band corresponds

13

to the three oscillation peaks observed at 11.2 GPa in VCC with a period of ~71 fs, or

14

a frequency of 470 cm-1. This frequency matches well the emergence of the new band

15

in the Ca PDOS between 400 and 500 cm-1. We can conclude unambiguously from

16

the analysis of the atomic dynamics that at high pressure (>11.2 GPa at 3000 K), the

17

Ca2+ and CO32- ions are compressed closer to each other and their motions are

18

correlated. As a consequence of the decrease in the diffusion rates, the viscosity of the

19

melt increases by almost one order of magnitude. This observation may have

20

important impact to the transport of melts under the mantle conditions.10, 16 Moreover,

21

the increased viscosity may affect the heat transport (thermal conductivity) in the

22

mantle as well. 21


Page 20 of 29

Page 21 of 29

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


1

4. Conclusions

2

In summary, the structures and viscosities of calcium carbonate melts under pressures

3

up to 52.5 GPa have been investigated with the FPMD method. The present study

4

reproduced the ultralow viscosity of calcium carbonate melts measured at low

5

pressures (< 11.2 GPa), The calculated velocity correlation functions show that

6

calcium carbonate melts behave like an ideal liquid at

Structures and Transport Properties of CaCO3 Melts under Earth's

Recommend Documents