Mechanical Properties of Single-Crystal Calcite and Their Temperature and Strain-Rate Effects

Calcite is the most stable crystalline phase of calcium carbonate. It is applied or found in composite products, the food industry, biomineralization, archaeology, and geology, and its mechanical properties have attracted more and more attention. In this paper, the mechanical behaviors of single-crystal calcite under uniaxial tension in different directions were simulated with the molecular dynamics method. The obtained elastic moduli are in good agreement with the experimental results. It has been found from further research that single-crystal calcite has typical quasi-brittle failure characteristics, and its elastic modulus, fracture strength, and fracture strain are all strongly anisotropic. The tensile failure is caused by dislocation emission, void formation, and phase transition along the [010] and [421] directions, but by continuous dislocation glide and multiplication along the [421¯] direction. The fracture strength, fracture strain, and elastic modulus are all sensitive to temperature, but only elastic modulus is not sensitive to strain rate. The effects of temperature and logarithmic strain rate on fracture strength are in good agreement with the predictions of fracture dynamics.


Introduction
Calcium carbonate (CaCO 3 ) is one of the most common of the mineral materials that provide important support for human life and industry. It occurs naturally in a wide variety of organisms, such as coccolithophores, foraminifers, and corals, as well as pearl oysters and animal bones [1,2]. Due to its biocompatibility and ease of use, calcium carbonate plays a very important role in biomedical aspects [3][4][5], such as regeneration [6,7]. It is also a biological mineral. When copolymer is added, its friction properties and other mechanical properties can be greatly improved [8,9]. In addition, it plays a key role in many anthropological and geological studies [10], materials protection [11], and accurate dating [12][13][14][15].
Calcium carbonate has three anhydrous crystal types, including calcite, aragonite, and vaterite. As a building material, calcite is found throughout nature. According to the statistics from European Calcium Carbonate Association (CCA-Europe), the global construction industry consumes hundreds of millions of tons of calcium carbonate every year. In addition, calcite accounts for more than 60% of the content of polymer fillers and is used to reinforce raw materials [16]. The addition of calcite not only improves the characteristics of raw materials, but also greatly reduces the product cost. These applications are inevitably affected by temperature and strain changes. Therefore, it is of great significance to understand the mechanical properties of calcite at different temperatures and strains.
At present, the research on calcite, from Bridgeman's pioneering work [17] to many static and dynamic experimental studies [18][19][20][21][22][23][24], and then to the recent atomic scale simulations of phase-transition processes [25][26][27][28], mainly focuses on its high-pressure phase-transition behavior. In addition, Sokolá et al. [29] studied the effects of calcite content on the mechanical properties of calcium carbonate at different temperatures. They found that calcite can reduce the shrinkage of samples during sintering, which can improve the characteristics of samples. Kunitake et al. [30] obtained biological calcite from the prismatic layer of mollusks and compared the hardness and modulus of biological calcite and geological calcite through an indentation experiment. They found that biological calcite has higher hardness, and pointed out that the hardening mechanism in which the dislocation motion is hindered causes the failure of calcite. Subsequently, Schuster et al. [31] twisted the calcite aggregates under high pressure and characterized the products by Xray diffraction. They found that the dominant slip plane 1014 has obvious dislocation characteristics.
Although some mechanical properties of calcite have been tested and some deformation mechanisms have been revealed, the research and understanding of these mechanisms are not complete. In this study, we used the molecular dynamics (MD) method to simulate the mechanical behavior of calcite under uniaxial tensile loading, studied the effects of temperature and strain rate on the mechanical properties of calcite, and tried to reveal its internal deformation mechanism at the atomic scale.

Interatomic Interaction Potential
An interatomic potential function describes the interaction between atoms. In MD simulations, the potential function plays a decisive role in the simulation results. For ionic materials, there are generally two interatomic potential function models: the rigid ion model (RIM) and the shell model (SM). In the RIM, the potential energy only depends on the distance of each atomic pair. In order to more accurately describe the polarization mechanism of negative ions, especially O 2− , the SM divides each negative ion into two particles, namely a core and a shell [32]. The polarization-coordinate-coupling method is used to shorten the simulation time and improve the convergence speed [33].
For calcite, a lot of RIMs [34][35][36][37][38][39][40] and SMs [35,37,[41][42][43] have been developed. They have different descriptions of C-O interaction. Although the SMs show the volume characteristics [41,42] of calcite crystals through the polarization of CO 3 2− , the addition of core and shell particles in O 2− significantly reduces the efficiency of calculation. Therefore, an available simulation model has a relatively small size.
The SM developed for calcite by Dove et al. [34] can give more consistent predictions of the mechanical properties with the experimental results. In this model, the total potential energy includes bonded and non-bonded parts, namely: In the bonded part, the intermolecular harmonic function is used to deal with the bond-bending between O-C and C-O in O-C-O, and the out-of-plane potential of the CO 3 2− group is realized by a four-body torsional term. Moreover: where k θ is the bond-bending constant, θ 0 is the equilibrium bond angle between O-C and C-O in the O-C-O, taken as 120 • here, θ is the actual bond angle, k t is the torsional constant, ϕ is the angle between two O-C-O planes in CO 3 2− group, taken as 180 • here, and subscripts i, j, k, and l represent different atoms. In the non-bonded part, Born-Mayer interaction potential is used to characterize the short-range interactions of Ca-O, C-O, and O-O, and the Coulomb potential of their long-range interactions. See: where A and ρ are Born-Mayer potential parameters, ε 0 is the dielectric constant, q is the charge of atom, and r is the distance between two atoms with cut-off radius of 10 Å. All the potential function parameters are listed in Table 1.

The Computer Model
Calcite belongs to the trigonal system, and the space group is R3c. In a hexagonal calcite unit cell (a = b = 4.990 Å, c = 17.061 Å, α = β = 90 • , and γ = 120 • ), there are 64 atoms ( Figure 1). A computer model of single-crystal calcite with a size of 50 nm × 50 nm × 25 nm and 4.8 million atoms in total was constructed by using the unit cell structure proposed by Graf [44], as shown in Figure 2 In the non-bonded part, Born-Mayer interaction potential is used to characterize the short-range interactions of Ca-O, C-O, and O-O, and the Coulomb potential of their longrange interactions. See: where A and ρ are Born-Mayer potential parameters, 0 ε is the dielectric constant, q is the charge of atom, and r is the distance between two atoms with cut-off radius of 10 Å. All the potential function parameters are listed in Table 1.

The Computer Model
Calcite belongs to the trigonal system, and the space group is 3 R c . In a hexagonal calcite unit cell (a = b = 4.990 Å, c = 17.061 Å, α = β = 90°, and γ = 120°), there are 64 atoms ( Figure 1). A computer model of single-crystal calcite with a size of 50 nm × 50 nm × 25 nm and 4.8 million atoms in total was constructed by using the unit cell structure proposed by Graf [44], as shown in Figure 2    ] directions, respectively. The NPT ensemble was also used in the tension process simulations. The pressure was controlled by air pressure and the temperature by a Nose-Hoover thermostat. The particle-particle-particle-mesh (PPPM) algorithm [46], with accuracy of 10 −6 , was adopted to solve the dynamic equations. hough the three stress-strain curves have very similar shapes, they obviously have different slopes of the elastic stage, different ultimate stresses or fracture strength, and different strains corresponding to the ultimate stresses. This reflects the remarkable anisotropy of single-crystal calcite in the mechanical properties. At the stage of elastic deformation, the stress increases with the strain in an approximately linear trend. When the fracture strength is reached, the curves enter into the fracture stage. The results indicate that the tensile failure behavior of single-crystal calcite exhibits typical quasi-brittleness characteristics along the three directions [47]. By fitting the elastic deformation stage with the least square method, the elastic modulus along the [010] direction is about 73.2 GPa, which is very close to the experimental result of 70. 23   The calculation was performed using the general MD simulation package LAMMPS [45]. The NPT ensemble was also used in the tension process simulations. The pressure was controlled by air pressure and the temperature by a Nose-Hoover thermostat. The particleparticle-particle-mesh (PPPM) algorithm [46], with accuracy of 10 −6 , was adopted to solve the dynamic equations. directions at 300 K temperature and 1 × 10 8 s −1 strain rate. It can be seen that although the three stress-strain curves have very similar shapes, they obviously have different slopes of the elastic stage, different ultimate stresses or fracture strength, and different strains corresponding to the ultimate stresses. This reflects the remarkable anisotropy of singlecrystal calcite in the mechanical properties. At the stage of elastic deformation, the stress increases with the strain in an approximately linear trend. When the fracture strength is reached, the curves enter into the fracture stage. The results indicate that the tensile failure behavior of single-crystal calcite exhibits typical quasi-brittleness characteristics along the three directions [47]. By fitting the elastic deformation stage with the least square method, the elastic modulus along the [010] direction is about 73.2 GPa, which is very close to the experimental result of 70.23 ± 1.35 GPa [48]. Moreover, the elastic moduli along the [421] and [421] directions are 58.5 GPa and 60.4 GPa respectively. They are about 20% smaller than that along the [010] direction.

Uniaxial Tensile Behaviors in Different Directions
In order to understand the atomistic mechanism of tensile deformation and failure behavior of single-crystal calcite, we extracted the atomic configurations at different stages of tensile deformation along the [010] direction, as shown in Figure 4. In this study, the radial distribution function (RDF) as well as the atomic coordination number was used to identify a new phase, as shown in Figure 5. The RDF analysis shows that Ca-Ca has the nearest-neighbor distance of 3.2 Å in the new phase, while it is 3.8 Å in calcite, which corresponds to the first peaks of RDFs in Figure 5. As an intermediate value between 3.2 Å and 3.8 Å, therefore, 3.5 Å was chosen as the cut-off radius to calculate the coordination number of calcium atoms. Figure 4a-c shows the atomic configurations when the strain reaches 0.084 in the elastic deformation stage, 0.121 at the rapid stress reduction stage, and 0.135 at the failure stage, respectively. It can be seen that at the elastic deformation stage, only some dislocations initiate from the middle of both sides along the [421] direction. After the curve enters into a sharp stress-declining stage with the increase of strain, the dislocations in both the sides develop inwards and an amorphous phase forms gradually. Meanwhile, dislocations are also initiating, as shown in Figure 4b. At the failure stage, the dislocations develop rapidly along the slip plane 1014 . As the most important cleavage plane of calcite crystals [49][50][51], it is orthogonal to the [421] direction. As the strain increases, the further development of dislocation leads to the occurrence of voids in the middle of both sides along the [421] direction, while an ordered new phase is constantly formed in the center of the model, as shown in Figure 4c. Although the new phase is different from calcite, it has a regular atomic arrangement and a fixed radial distribution function, as shown in Figure 5. For observation, the atoms are colored according to the coordination number of calcium atoms in Figure 4, with green for calcite, blue for the new phase, and cyan for the amorphous phase, respectively. The development of the new phase and voids is the fundamental reason for the rapid stress reduction at the second stage [47]. The failure process of calcite in tension along the [010] direction originates from initial dislocation emission and ends with the rapid development of the new phase and voids.  In order to understand the atomistic mechanism of tensile deformation and failure behavior of single-crystal calcite, we extracted the atomic configurations at different stages of tensile deformation along the [010] direction, as shown in Figure 4. In this study, the radial distribution function (RDF) as well as the atomic coordination number was used to identify a new phase, as shown in Figure 5. The RDF analysis shows that Ca-Ca has the nearest-neighbor distance of 3.2 Å in the new phase, while it is 3.8 Å in calcite, which corresponds to the first peaks of RDFs in Figure 5. As an intermediate value between 3.2 Å and 3.8 Å, therefore, 3.5 Å was chosen as the cut-off radius to calculate the coordination number of calcium atoms. Figure 4a-c shows the atomic configurations when the strain reaches 0.084 in the elastic deformation stage, 0.121 at the rapid stress reduction stage, and 0.135 at the failure stage, respectively. It can be seen that at the elastic deformation stage, only some dislocations initiate from the middle of both sides along the [421] direction. After the curve enters into a sharp stress-declining stage with the increase of strain, the dislocations in both the sides develop inwards and an amorphous phase forms gradually. Meanwhile, dislocations are also initiating, as shown in Figure 4b. At the failure stage, the dislocations develop rapidly along the slip plane ( ) 10 14 . As the most important cleavage plane of calcite crystals [49][50][51], it is orthogonal to the [421] direction. As the strain increases, the further development of dislocation leads to the occurrence of voids in the middle of both sides along the [421] direction, while an ordered new phase is constantly formed in the center of the model, as shown in Figure 4c. Although the new phase is different from calcite, it has a regular atomic arrangement and a fixed radial distribution function, as shown in Figure 5. For observation, the atoms are colored according to the coordination number of calcium atoms in Figure 4, with green for calcite, blue for the new phase, and cyan for the amorphous phase, respectively. The development of the        In order to understand the internal mechanism of dislocation slip induced by the phase transition in uniaxially stretched calcite, Figure 7 shows the process from the dislocation initiation and development to the formation of the new ordered phase. Compared with the initial perfect calcite structure shown in Figure 7a, when the stress reaches the fracture strength, a dislocation line is formed in the center of the model, as shown in Figure 7b, which corresponds to Figure 4b. In Figure 7c, the rapid propagation and development of dislocations causes the new phase to form. The Ca 2+ and CO3 2− interaction inside calcium carbonate molecule is much stronger than the interaction between calcium carbonate molecules, so dislocations easily occur between two adjacent calcium carbonate molecules. As reported in reference [52], they tend to break weaker bonds in the cleavage plane. In order to understand the internal mechanism of dislocation slip induced by the phase transition in uniaxially stretched calcite, Figure 7 shows the process from the dislocation initiation and development to the formation of the new ordered phase. Compared with the initial perfect calcite structure shown in Figure 7a, when the stress reaches the fracture strength, a dislocation line is formed in the center of the model, as shown in Figure 7b, which corresponds to Figure 4b. In Figure 7c, the rapid propagation and development of dislocations causes the new phase to form. The Ca 2+ and CO 3 2− interaction inside calcium carbonate molecule is much stronger than the interaction between calcium carbonate molecules, so dislocations easily occur between two adjacent calcium carbonate molecules. As reported in reference [52], they tend to break weaker bonds in the cleavage plane. fracture strength, a dislocation line is formed in the center of the model, as shown in Figure 7b, which corresponds to Figure 4b. In Figure 7c, the rapid propagation and development of dislocations causes the new phase to form. The Ca 2+ and CO3 2− interaction inside calcium carbonate molecule is much stronger than the interaction between calcium carbonate molecules, so dislocations easily occur between two adjacent calcium carbonate molecules. As reported in reference [52], they tend to break weaker bonds in the cleavage plane.   Figure 8a, the calcium carbonate molecules in the dotted blue lines, such as molecules 1, 4, 7, and 10 or 2, 5, 8, and 11, act as a close-packed plane. There is a high molecular density on these planes. Due to large distance between two neighboring close-packed planes, there is a weak atomic binding force and a small slip resistance between them, so that the slip preferentially occurs.   Figure 8a, the calcium carbonate molecules in the dotted blue lines, such as molecules 1, 4, 7, and 10 or 2, 5, 8, and 11, act as a close-packed plane. There is a high molecular density on these planes. Due to large distance between two neighboring closepacked planes, there is a weak atomic binding force and a small slip resistance between them, so that the slip preferentially occurs.  Figure 9b. This phenomenon can also be explained by the arrangement of calcium carbonate molecules. As shown in Figure 8a, the calcium carbonate molecules in the dotted red lines, such as molecules 4, 5, and 6 or 7, 8, and 9, act as a close-packed plane, and are stretched in a direction just perpendicular to the closepacked plane, resulting in a direct separation between the two close-packed planes and the formation of voids. It is more similar to the tensile failure of single-crystal aragonite along the [010] direction [53]. As the strain increases, the voids continue to grow up and phase transition takes place in the center of the model, as shown in Figure 9c. In the new phase, Ca-Ca has the nearest neighbor distance of 3.4 Å, which is different from the phase with the Ca-Ca nearest neighbor distance of 3.2 Å shown in Figure 4c.  Figure 9b. This phenomenon can also be explained by the arrangement of calcium carbonate molecules. As shown in Figure 8a, the calcium carbonate molecules in the dotted red lines, such as molecules 4, 5, and 6 or 7, 8, and 9, act as a close-packed plane, and are stretched in a direction just perpendicular to the close-packed plane, resulting in a direct separation between the two close-packed planes and the formation of voids. It is more similar to the tensile failure of single-crystal aragonite along the [010] direction [53]. As the strain increases, the voids continue to grow up and phase transition takes place in the center of the model, as shown in Figure 9c. In the new phase, Ca-Ca has the nearest neighbor distance of 3.4 Å, which is different from the phase with the Ca-Ca nearest neighbor distance of 3.2 Å shown in Figure 4c.
A similar quasi-brittle fracture can also be observed when the model is stretched along the [421] direction. Compared with the atomic configuration at the elastic deformation stage shown in Figure 10a, some dislocations occur on both sides along the [421] direction at the rapid stress-reduction stage, as shown in Figure 10b, and then rapidly multiply and are emitted to the center of the model at the failure stage, as shown in Figure 10c. As shown in Figure 8b, the calcium carbonate molecules in the dotted yellow lines, such as molecules 1, 2, 3, and 4; or 5, 6, 7, and 8, act as a close-packed plane. Different 010] direction, there are a lot of voids that form in advance at the rapid stress reduction stage, as shown in Figure 9b. This phenomenon can also be explained by the arrangement of calcium carbonate molecules. As shown in Figure 8a, the calcium carbonate molecules in the dotted red lines, such as molecules 4, 5, and 6 or 7, 8, and 9, act as a close-packed plane, and are stretched in a direction just perpendicular to the closepacked plane, resulting in a direct separation between the two close-packed planes and the formation of voids. It is more similar to the tensile failure of single-crystal aragonite along the [010] direction [53]. As the strain increases, the voids continue to grow up and phase transition takes place in the center of the model, as shown in Figure 9c. In the new phase, Ca-Ca has the nearest neighbor distance of 3.4 Å, which is different from the phase with the Ca-Ca nearest neighbor distance of 3.2 Å shown in Figure 4c. direction at the rapid stress-reduction stage, as shown in Figure 10b, and then rapidly multiply and are emitted to the center of the model at the failure stage, as shown in Figure  10c. As shown in Figure 8b, the calcium carbonate molecules in the dotted yellow lines, such as molecules 1, 2, 3, and 4; or 5, 6, 7, and 8, act as a close-packed plane. Different from

Influence of Temperature and Strain Rate
As with other inorganic materials, the mechanical properties of calcite are affected by temperature and strain rate [54,55]. In order to study the effect of temperature on the deformation behavior, the strain rate was fixed at 1 × 10 8 s −1 but the temperature was changed from 10 K to 100 K, 200 K, 300 K, 400 K, and 500 K. Figure 11 shows

Influence of Temperature and Strain Rate
As with other inorganic materials, the mechanical properties of calcite are affected by temperature and strain rate [54,55]. In order to study the effect of temperature on the deformation behavior, the strain rate was fixed at 1 × 10 8 s −1 but the temperature was changed from 10 K to 100 K, 200 K, 300 K, 400 K, and 500 K. Figure 11 shows Figure 11. Variations of (a) fracture strength, (b) fracture strain, and (c) elastic modulus with temperature.
In order to investigate the internal mechanisms of temperature effect, the atomic configurations at 100 K, 300 K, and 500 K when the tensile strain reaches 0.125 along the [010] direction are given in Figure 12. It can be seen that when the temperature elevates from 100 K to 300 K and 500 K, the volume of the new phase expands rapidly, which should be the main reason for the decline of fracture strength, fracture strain, and elastic modulus caused by the elevated temperature. In order to investigate the internal mechanisms of temperature effect, the atomic configurations at 100 K, 300 K, and 500 K when the tensile strain reaches 0.125 along the [010] direction are given in Figure 12. It can be seen that when the temperature elevates from 100 K to 300 K and 500 K, the volume of the new phase expands rapidly, which should be the main reason for the decline of fracture strength, fracture strain, and elastic modulus caused by the elevated temperature. In order to study the effect of strain rate on the mechanical properties of calcite, the temperature was fixed at 300 K, but the strain rate was changed from 0.00001 ps −1 to 0.0001 ps −1 , 0.001 ps −1 , 0.01 ps −1 , and 0.1 ps −1 . Figure 13 shows the variations of fracture strength, fracture strain, and elastic modulus with the strain rate along the .4% respectively, and the fracture strain along the directions increases by 16.6%, 15.6%, and 17.2% respectively. It can be seen that both the fracture strength and fracture strain increase monotonously with the increase of strain rate. However, the strain rate hardly has any effect on the elastic modulus. In order to study the effect of strain rate on the mechanical properties of calcite, the temperature was fixed at 300 K, but the strain rate was changed from 0.00001 ps −1 to 0.0001 ps −1 , 0.001 ps −1 , 0.01 ps −1 , and 0.1 ps −1 . Figure 13 14.4%, and 13.4% respectively, and the fracture strain along the directions increases by 16.6%, 15.6%, and 17.2% respectively. It can be seen that both the fracture strength and fracture strain increase monotonously with the increase of strain rate. However, the strain rate hardly has any effect on the elastic modulus. In order to study the effect of strain rate on the mechanical properties of calcite, the temperature was fixed at 300 K, but the strain rate was changed from 0.00001 ps −1 to 0.0001 ps −1 , 0.001 ps −1 , 0.01 ps −1 , and 0.1 ps −1 . Figure 13 15.6%, and 17.2% respectively. It can be seen that both the fracture strength and fracture strain increase monotonously with the increase of strain rate. However, the strain rate hardly has any effect on the elastic modulus. In order to investigate the internal mechanisms of strain rate effect, the atomic configurations with tensile strain of 0.125 along the [010] direction are given in Figure 14 for strain rates of 0.00001 ps −1 , 0.001 ps −1 , and 0.1 ps −1 . It can be seen that with the increase of strain rate, the volume of the new phase is noticeably reduced, which should be the main reason for the increase of the fracture strength and fracture strain caused by increasing strain rate. In general, the increase of strain rate leads to a shorter time for the atoms in calcite to respond to the load, so that they have less chance to overcome the energy barrier and destroy chemical bonds.

Fracture Dynamics Analysis
Fracture dynamics can explain the temperature and strain rate effects of bond failure through thermal activation. Arrhenius regarded the life cycle τ as the following function of tensile stress σ and temperature T [56]: where 0 τ is the average vibration period of atoms, In order to investigate the internal mechanisms of strain rate effect, the atomic configurations with tensile strain of 0.125 along the [010] direction are given in Figure 14 for strain rates of 0.00001 ps −1 , 0.001 ps −1 , and 0.1 ps −1 . It can be seen that with the increase of strain rate, the volume of the new phase is noticeably reduced, which should be the main reason for the increase of the fracture strength and fracture strain caused by increasing strain rate. In general, the increase of strain rate leads to a shorter time for the atoms in calcite to respond to the load, so that they have less chance to overcome the energy barrier and destroy chemical bonds. In order to investigate the internal mechanisms of strain rate effect, the atomic configurations with tensile strain of 0.125 along the [010] direction are given in Figure 14 for strain rates of 0.00001 ps −1 , 0.001 ps −1 , and 0.1 ps −1 . It can be seen that with the increase of strain rate, the volume of the new phase is noticeably reduced, which should be the main reason for the increase of the fracture strength and fracture strain caused by increasing strain rate. In general, the increase of strain rate leads to a shorter time for the atoms in calcite to respond to the load, so that they have less chance to overcome the energy barrier and destroy chemical bonds.

Fracture Dynamics Analysis
Fracture dynamics can explain the temperature and strain rate effects of bond failure through thermal activation. Arrhenius regarded the life cycle τ as the following function of tensile stress σ and temperature T [56]: where 0 τ is the average vibration period of atoms, U is the dissociation energy of interatomic chemical bonds, V is the activation volume, and c is the local overstress coefficient. In a uniaxial tensile process, the stress changes with time, that is According to Bailey's principle [57], the fracture time r t can be obtained by:

Fracture Dynamics Analysis
Fracture dynamics can explain the temperature and strain rate effects of bond failure through thermal activation. Arrhenius regarded the life cycle τ as the following function of tensile stress σ and temperature T [56]: where τ 0 is the average vibration period of atoms, k = 1.381 × 10 −23 J/K is the Boltzmann constant, U 0 is the dissociation energy of interatomic chemical bonds, V is the activation volume, and c is the local overstress coefficient. In a uniaxial tensile process, the stress changes with time, that is σ = σ(t). According to Bailey's principle [57], the fracture time t r can be obtained by: For linear elastic materials, we have: By substituting Equations (4), (6), and (7) into Equation (5), the fracture strength can be given by: where γ = cV. This shows that there is a linear relationship between the fracture strength and the logarithmic strain rate when the temperature is fixed. Equation (8) can be rewritten as:  (9) can be reduced to: For calcium carbonate, γ = 5.67 × 10 −29 m 3 . When T = 300 K, kT γ = 0.073 GPa in Equation (8), which lies in the range of 0.0671 to 0.1137 GPa, given by the above MD simulations. In addition, −26.3 k γ = −0.0064 GPa in Equation (10), which is between −0.0073 GPa and −0.0048 GPa, given by the above MD simulations. Figures 15 and 16 (8) and (10). Obviously, the MD simulations are in good agreement with the fracture dynamics analysis.

Conclusions
In this paper, the mechanical properties of single-crystal calcite, as well as their temperature and strain-rate effects, were studied by MD simulations. Single-crystal calcite's tensile failure behavior exhibits typical quasi-brittleness characteristics. The elastic moduli obtained by the simulations are very close to the experimental value, which confirms the reliability of MD simulations. The fracture strength, fracture strain, and elastic modulus are all highly anisotropic. The fracture strength and elastic modulus along the [010] direction is much higher than those along the [421] and [42 1] directions. High anisotropy is caused by different failure and deformation mechanisms in different directions, including dislocation initiation and development, phase transition, and void formation. It is found that the fracture strength, fracture strain, and elastic modulus are all very sensitive to temperature, but only elastic modulus is not sensitive to strain rate. The variations of fracture strength with temperature and logarithmic strain rate are in good agreement with the predictions of fracture dynamics.

Conclusions
In this paper, the mechanical properties of single-crystal calcite, as well as their temperature and strain-rate effects, were studied by MD simulations. Single-crystal calcite's tensile failure behavior exhibits typical quasi-brittleness characteristics. The elastic moduli obtained by the simulations are very close to the experimental value, which confirms the reliability of MD simulations. The fracture strength, fracture strain, and elastic modulus are all highly anisotropic. The fracture strength and elastic modulus along the [010] direction is much higher than those along the [421] and [421] directions. High anisotropy is caused by different failure and deformation mechanisms in different directions, including dislocation initiation and development, phase transition, and void formation. It is found that the fracture strength, fracture strain, and elastic modulus are all very sensitive to temperature, but only elastic modulus is not sensitive to strain rate. The variations of fracture strength with temperature and logarithmic strain rate are in good agreement with the predictions of fracture dynamics.
Funding: This work was supported by the National Natural Science Foundation of China (Grant Nos 12072120 and 11832013).