Yadong Lia,
Xiaopeng Jiaa,
Bingmin Yanb,
Ning Chena,
Chao Fanga,
Yong Lic,
Shishuai Sund and
Hongan Ma*a
aState Key Lab of Superhard Materials, Jilin University, Changchun 130012, China. E-mail: jlu_li1989@163.com; Fax: +86-431-85168858; Tel: +86-431-85168858
bCenter for High Pressure Science & Technology Advanced Research, Changchun 130012, China
cPhysical and Applied Engineering Department, Tongren University, China
dCollege of Science, Tianjin University of Technology, Tianjin 300384, China
First published on 19th April 2016
The temperature and convection fields of a catalyst with three different heights were simulated in a temperature gradient growth (TGG) system under high pressure and high temperature (HPHT) conditions. Temperature fields were simulated to rule out the influence of temperature on the crystal morphology. The features of the calculated convection field could predict the particular situations of diamond growth systems very well, and could explain the change of diamond morphology accompanying the growth process. According to the calculated results, we predict that the morphology of a diamond crystal changes from cubic crystal to cub-octahedral. A good agreement has been obtained between the calculated results and the observed experimental data. The morphology and structural properties of the synthesized samples are characterized by optical microscopy and a Raman spectrum. The results illustrate that the synthesized diamond crystals have less lattice distortion and with high quality.
Generally, the district for diamond growth is a V-shape region bounded by diamond–graphite equilibrium line and solvent–carbon eutectic melting line in the metal solvent–carbon system.10 The synthesized diamond usually present the {100} face at lower growth temperatures, while the morphology of diamond is dominantly by {111} face at higher temperatures. Therefore, changing the shape of diamond crystal is mainly by altering the synthetic temperature in traditional methods. Besides, during the crystal growth process, the carbon convection field around the seed diamond also has significant effect on the morphology of the diamond.17 In order to understand and optimize the growth process of diamond single crystals using the TGG method, it is necessary to investigate the detailed information about the nature of the convection field of particular growth systems. However, it is very difficult to characterize the convection field experimentally, such as the flow tendency and intensity of carbon in the molten catalyst/solvent under HPHT conditions. Until now, there have been very few reports about the effect of carbon convection on the change of the morphology of diamond crystals.
In this paper, we have simulated the temperature field, convection field and growth process of the diamond crystal by the application of finite element method (FEM). The simulation successfully explained the influence of the carbon convection field on diamond morphology in diamond synthesis process. The simulation results indicated that the morphology of diamond crystals could be adjusted by changing the carbon convection field. At the same time, the experimental data obtained from the corresponding synthesis experiments were consistent with the simulation results. The synthesized diamonds were characterized by optical microscope (OM) and Raman spectrum.
As shown in Fig. 1(a), the growth cell was heated when voltage V1 was applied on the base of the up steel ring and V2 was set on the base of the low steel ring which contacted a graphite heater and the WC anvils. The bottoms of the six steel rings were kept at a uniform temperature. The velocity constraints on the symmetric surface were symmetry constraints and the other areas of 1/4 catalyst were set as VX = VY = VZ = 0, so that only convection in the catalyst region was allowed. The material parameters and boundary conditions used in the finite element simulation can be found in Tables 1 and 2.18–20
Boundary parameters | Value |
---|---|
The temperature of circulation water (°C) | 45 |
The temperature of environment (°C) | 35 |
Sample assembly initial temperature (°C) | 60 |
Terminal surface temperature of the steel ring (°C) | 50 |
Material parameters | Value |
---|---|
Density (kg m−3) | |
Diamond | 3510 |
Pyrophyllite | 2950 |
Steel | 7900 |
WC | 14![]() |
Nickel | 9331–1.046T |
Graphite | 2200 |
Insulator | 5600 |
![]() |
|
Specific heat (J kg−1 °C−1) | |
Diamond | 471 |
Pyrophyllite | 900 |
Steel | 450 |
WC | 200 |
Nickel | 460 |
Graphite | 710 |
Insulator | 450 |
![]() |
|
Electric resistivity (Ω m) | |
Diamond | 1.0 × 1014 |
Pyrophyllite | 5.0 × 109 |
Steel | 9.71 × 10−8 |
WC | 5.83 × 10−7 |
Nickel | 1.15 × 10−8 |
Graphite | 2.46 × 10−6 |
Insulator | 8.30 × 1013 |
![]() |
|
Thermal conductivity (W m−1 K−1) | |
Diamond | 470 |
Pyrophyllite | 9.04 |
Steel | 79 |
WC | 160 |
Nickel | 20 |
Graphite | 200 |
Insulator | 2.09 |
![]() |
|
Viscosity (Pa s) | |
Nickel | 0.041–3.478 × 10−5T + 7.941 × 10−9T2 |
Temperature has significant effect on diamond morphology. The synthesized diamond usually present the {100} face at lower growth temperatures, while the morphology of diamond is dominantly by {111} face at higher temperatures. Fig. 2 shows the distributions of temperature fields on different heights of catalyst and carbon source. From the simulation results in Fig. 2, we can see that the temperatures of the three seed locations are all 1263 °C, which means that the crystal morphology should be same. Besides, through the color of the cloud of Fig. 2(a)–(c), we can infer that the temperature at the top of the catalyst abutting the carbon source become higher and higher. The solubility of carbon in molten metal also increased with the temperature rising. The carbon concentration modes varied with different growth cells because of the change of solubility gradient in the solvent fluid21 and the difference among the increased temperature may have influence on the carbon convection field in molten metallic solvent.
The calculated distributions of the carbon convection are shown in Fig. 3. In the simulated section figure, the arrowheads indicate the flow direction of carbon, the length and the color of the arrows represent the magnitude of carbon flow velocity. Therefore, the relationship between carbon convection distributions and colors (or lengths) of arrowheads are gained. First of all, it is clear that there is a region near the periphery of the solvent where the carbon convection rate is the highest with red arrows. Additionally, it is important to find that the vortex center (black circle with direction) of convection field move closer to the solvent center gradually as the solvent height increases. It is obviously to see that the colors of the arrowheads become deepen and the lengths become larger in the center area of the catalyst/solvent around the seed crystal and the arrow density and length are symmetrically distributed along the vertical direction (inside the blue dotted line).
![]() | ||
Fig. 3 Distributions of convection field of the catalyst with three different heights: (a) 2.0 mm (b) 2.4 mm (c) 2.8 mm. |
According to the melt–solvent theory, the temperature gradient between the carbon source and diamond seed is the driving force of diamond growth, the carbon element can transport from the carbon source to the crystal surface by temperature gradient. The convection is caused under the action of temperature gradient in the gravitational field. Thus, the convection is also an important factor for diamond growth by HPHT. Convection in crystal growth cell provide carbon transport from carbon source to seed crystal face, so the bigger convection velocity of the carbon, the faster transport of the carbon atoms from source to the seed. From the Fig. 3(a)–(c), we can see that the convection of carbon is deflect and concentrate on the seed crystal, especially in the region near the seed bed and the color and length of arrowheads increase in the centre of the solvent along Y-axis. That is to say, the driving force of diamond crystal growth also increased gradually in the centre of the solvent along Y-axis. Therefore, in the process of diamond growth, the faster carbon transport from source to the seed surface can lead to the quick growth of {100} face. In order to further study the situations of the convection field in the catalyst, we compared the X-axis path and Y-axis path convection velocity. Fig. 4 plots the convection velocity on the X-axis path and Y-axis path. It shows that the convection velocity in X-axis and Y-axis of different heights catalyst both increased (change from 2.0 mm to 2.8 mm). Besides, compared to the X-axis, the convection velocity of Y-axis path grows faster than X-axis. In the diamond crystal growth cell, convection provides a driving force for carbon transport from the carbon source to the crystal surface. Thus, the increased velocity of convection can enhance the carbon transport capacity. As we have mentioned before, the solubility and concentration of carbon in 2.8 mm catalyst are higher than that of in 2.4 mm and 2.0 mm catalyst near the top surface of the catalyst. Thus, there are more carbon can be transported to seed crystal in the growth cell of 2.8 mm catalyst. According to the crystal growth theory, the face growing rapidly disappears firstly and leaves the face which grows slowly.22 As we know, {111} and {100} appear most dominantly on synthetic diamond. So we can predict from the simulated results that convection field of Fig. 3(a) is suitable for growing diamond dominantly with {100} faces, Fig. 3(b) is suitable for {100} and {111} faces morphology diamond and Fig. 3(c) is suitable for diamond with dominated {111} faces and less {100} faces.
![]() | ||
Fig. 4 Convection velocity of different heights of catalyst. (a) The convection velocity on the X-axis path. (b) The convection velocity on the Y-axis path. |
Since the diamond itself is a good heat conductor, and with the growth of diamond, it will alters the temperature distribution in the cell. To further study the influence of temperature and convection field on crystal growth process, we have simulated the temperature and convection field of the initial stage for 5 h. Fig. 5 shows the temperature field of diamond growth of 5 h. It can be seen that the distribution of temperature field is similar to Fig. 2. The temperature in seed location is the same and the temperature at the top of the catalyst in Fig. 5(c) is highest. Fig. 6 shows the convection field of diamond growth of 5 h. The tendency of convection distribution is consistent with the previous results. Therefore, we believe that the prediction that the convection fields have effect on the crystal morphology is credible.
![]() | ||
Fig. 5 The distribution of temperature field in initial stage of diamond growth for 5 h. (a) 2.0 mm (b) 2.4 mm (c) 2.8 mm. |
![]() | ||
Fig. 6 The distribution of convection field in initial stage of diamond growth for 5 h. (a) 2.0 mm (b) 2.4 mm (c) 2.8 mm. |
![]() | ||
Fig. 7 Synthesized diamond crystals in different heights of Ni-based catalyst under 5.7 GPa and 1260 °C. (a) 2.0 mm (b) 2.4 mm (c) 2.8 mm. |
Raman spectroscopy was used to characterize crystal lattice distortion or crystalline quality of diamond crystals. The Raman peak values and the full width at half maximum (FWHM) indicated the relative quality of crystals. Fig. 8 shows Raman spectra recorded from diamond crystals of our sample obtained. The diamond Raman peaks values of those diamond crystals are near 1332 cm−1 and the FWHM of the diamonds are 3.8 cm−1. It is supposed that the observed features of the Raman spectrum of diamond are related to the increasing concentrations of micro-crystallites in boundary, impurities, residual stress, or defects, etc. Therefore, our results illustrate that the change of the convection field only has influence on the crystal morphology and no effect on the crystal quality of diamond.
![]() | ||
Fig. 8 Raman spectra of diamond crystals for the samples obtained. (a) 2.0 mm (b) 2.4 mm (c) 2.8 mm. |
This journal is © The Royal Society of Chemistry 2016 |