Photoluminescence Mapping of Optical Defects in HPHT Synthetic Diamond
ABSTRACT
Photoluminescence (PL) mapping provides a means to identify the distribution of optical centers in diamond. To demonstrate the impact of this method on the field of gemology and the study of diamonds, photoluminescence maps were acquired from a laboratory-irradiated brownish orange HPHT synthetic diamond of mixed diamond type. Acquisition time for each PL map was less than four minutes. Analysis of the maps confirmed that optical centers are incorporated in diamond growth sectors, such as {111} octahedral, {100} cubic, {110} dodecahedral, and {113} trapezohedral. The remarkable correlation between optical defects and growth sectors in the sample demonstrates that high-speed photoluminescence mapping is ideal for rapidly determining the distribution of optical defects resulting from both diamond growth and treatments. We anticipate that this technique will allow gemological labs to continue to identify increasingly sophisticated synthetic diamonds and color treatments, helping to ensure consumer confidence in the diamond industry.
INTRODUCTION
Optical centers in a diamond’s crystal lattice, also known as point defects, absorb visible light to produce color (Collins, 1982). In this way, optical centers give rise to the variety of fancy-color diamonds available in the trade. Optical centers can form as a result of impurities in the diamond (such as nitrogen or boron), deformation of the crystal lattice, missing carbon atoms (vacancies), or a combination of these (Collins, 2003; Breeding and Shigley, 2009). Photoluminescence occurs when an optical center is excited to a higher energy state by the absorption of photons and then returns to its ground energy state, emitting light at a specific wavelength or band of wavelengths (Collins, 1992; Eaton-Magaña and Breeding, 2016). There are currently several hundred known optical centers in diamond that produce photoluminescence (e.g., Collins, 1999; Zaitsev, 2001).
Photoluminescence (PL) spectroscopy is a powerful tool used in gemological laboratories to verify whether a specimen is natural or synthetic and determine whether its color origin is natural or due to treatment (Eaton-Magaña and Breeding, 2016). One drawback is that spectra are often collected from a single spot, or a few spots, on the sample surface and therefore represent only a tiny fraction of the entire sample (figure 1). In general, the spectra collected are sufficient for a trained gemologist to determine the diamond’s origin, especially when the PL spectra are interpreted in concert with DiamondView images, Fourier-transform infrared (FTIR) spectra, and other absorption and luminescence data. However, as synthetic diamond technology improves for both HPHT- and CVD-grown diamonds (e.g., Moe et al., 2016; Wang and Poon, 2016; Law and Wang, 2016; Eaton-Magaña and Shigley, 2016), and color treatments advance (e.g., Moe et al., 2015; Ardon and Loudin, 2016), gemological laboratories must introduce new techniques to continue to identify all synthetic and treated diamonds.
One useful technique is PL mapping (Eaton-Magaña and Breeding, 2016), which can be used to plot the distribution of an optical defect’s intensity. A PL mapping microscope can collect thousands of spectra in a raster pattern by either continuously scanning the sample’s surface or by moving from one point to another (pixel to pixel) and collecting a spectrum at each point. The spectral data are then interpolated and combined to produce a map of a defect’s distribution (again, see figure 1). PL mapping was first used for gemological purposes in the summer of 2015; the first public demonstration of the technique was by Johnson et al. (2015) at the annual Geological Society of America conference in Baltimore. The technique was further developed by Dieck et al. (2015) to show the distribution of the silicon vacancy (SiV–) in a CVD synthetic diamond. Ardon and Loudin (2016) used PL mapping to show the influence of focused beam irradiation on the distribution of nitrogen vacancies. Most recently, Johnson and Myagkaya (2017) used PL mapping to show the relationship between SiV– and Ni+ in an HPHT synthetic diamond.
One major challenge to implementing a PL mapping technique is that PL emission is temperature dependent. Therefore, the diamond must be cooled to liquid nitrogen temperature (approximately 77 K, or –321°F) to produce sharp PL bands (Collins, 1982; Eaton-Magaña and Breeding, 2016). Overcoming this challenge requires that the diamond be cooled in either a cryogenic stage or an open bath of liquid nitrogen. If the open-bath cooling method is used, then a high-speed mapping technique is needed because the liquid nitrogen boils off and vaporizes at atmospheric temperature and pressure. Using a Raman imaging microscope equipped with an electron-multiplying charge-coupled device (EMCCD), we were able to efficiently produce PL maps in less than four minutes, allowing us to overcome the temperature/time constraints and evaluate the distribution of optical centers in our synthetic diamond sample. Here we present an introduction to PL mapping and demonstrate the technique’s effectiveness by reporting the distribution of optical defects in an HPHT synthetic diamond.
MATERIALS AND METHODS
GIA’s New York laboratory received an irradiated 1.06 ct Fancy Deep brownish orange HPHT synthetic diamond of mixed type (figure 2) for a synthetic colored diamond grading report. Absorption features typical of both type IIb and type Ib diamonds were observed in the specimen’s infrared absorption spectrum (figure 3).
The type IIb component, attributed to boron in the diamond lattice, was identified by absorptions at 2803 and 2930 cm–1. The type Ib component, composed of isolated nitrogen in the form of C centers, was identified by infrared absorptions at 1130 and 1344 cm–1. DiamondView images showed growth sectors corresponding to the {100} cubic, {111} octahedral, {110} dodecahedral, and {113} trapezohedral sectors (figure 4). PL spectra collected using a Renishaw inVia Raman microscope with 488, 514, and 830 nm laser excitation wavelengths revealed several PL features, including well-known defects such as NV0 at 575 nm and NV– at 637 nm, as well as several lesser-known defects such as 550.5 nm (figure 5). The combination of well-defined growth sectors and abundant optical defects made this sample well suited for demonstrating the effectiveness of PL imaging microscopy in mapping the distribution of optical defects in diamond.
PL mapping was performed at GIA’s New York laboratory. Prior to mapping, the pavilion facets were coated with Pelco colloidal graphite to reduce the amount of light reflecting off the pavilion facet junctions and returning to the table. Additionally, the sample was mounted with the table level. Mapping was conducted using a Thermo Scientific DXRxi Raman imaging microscope with 455, 532, 633, and 780 nm laser excitation wavelengths. The DXRxi is equipped with an Olympus optical microscope and an Andor Technology Newton 970 EMCCD. An Olympus 10×, 0.25 numerical aperture objective lens was used.
The DXRxi uses a continuously moving, variable-speed sample stage driven by a combination of linear magnetic motors and controlled by optical encoders. The stage movement is synchronized with the EMCCD detector exposures so that the spectra are collected at the desired positions with repeatability within 100 nm (A. Rzhevskii, pers. comm., 2016). For each laser excitation, the sample was analyzed at liquid nitrogen temperature to produce sharp PL peaks. The combination of an EMCCD detector and a variable-speed, exposure-synchronized stage allowed the DXRxi to produce high-quality PL maps while maintaining the appropriate temperature and level of liquid nitrogen. Table 1 summarizes the analytical conditions used for each laser excitation. The data were processed using Thermo Scientific’s OMNICxi analysis software package, and baseline-corrected peak area profiles were used to produce the observed PL maps.
TABLE 1. Analytical conditions and Raman imaging microscope settings for each laser excitation wavelength. | ||||
Laser wavelength (nm) | ||||
455 | 532 | 633 | 780 | |
Grating (lines/mm) | 1200 | 400 | 600 | 400 |
Laser power (mW) | 5.9 | 0.3 | 8 | 24 |
Exposure time (ms) | 1.67 | 1.67 | 1.67 | 8.33 |
Pinhole aperture size (μm) | 50 | 25 | 50 | 50 |
Pixel size (μm) | 15 | 15 | 15 | 15 |
EM gain | Off | Off | Off | Off |
Total spectra | 61,005 | 60,270 | 59,535 | 21,756 |
Collection time (min) | 3.4 | 3.4 | 3.3 | 3.8 |
RESULTS
For each laser excitation wavelength, representative PL maps of the baseline-corrected peak area intensity are shown in figure 6. The maps are interpreted using a color scale. Warmer colors indicate a region that contains a greater intensity of a defect, while cooler colors represent lower intensities of that defect. The scale for each map depends on the defect’s minimum and maximum intensity. To fully interpret the map, the spatial distribution of the defect must be taken in combination with the PL spectra. Figure 7 shows the average spectra from the {100} and {111} growth sectors around points 1 and 2 in figure 6A. Table 2 summarizes the peak area distribution per growth sector for a majority of the peaks observed in the DXRxi spectra. The data presented here are qualitative and are used for the purpose of demonstrating the technique.
TABLE 2. Relative peak distribution per growth sectormeasured by the Raman imaging microscope usingvarious laser excitation wavelengths. | ||
Laser | Peak position (nm) | Peak intensitya per growth sector |
455 nm | 470.1 489.4 490.2 490.9 492.0 503.2 506.1 508.8 517.6 537.8 |
{100}>{113}>{111}>>{110} {111}≥{100}>>{113} {111}>{100}>{113} {111}≥{100}>{113} {111}≥{100}>{113} {111}≥{100}>{113}>>{110} {113}≥{110} {100}>{113}>{111}>>{110} (fig. 6A) {111}>>{110} {111} |
532 nm | 550.7 574.9 636.9 741.0 |
{111}>{100} {111}>{100}>>{113}>{110} (fig. 6B) {111}>{100}>>{113}>{110} {110}>{113} |
633 nm | 637.0 648.1 662.5 676.2 681.1 684.0 718.2 741.1 |
{111}>{100}>{113}≈{110} (fig. 6C) {111} {100}>{111}>>{113} {100}≈{111}>>{113} {100}>{111}>{113} {110}>{113} {111}>{100}>{113} {110}>{113}>{100}>>{111} |
780 nm | 795.7 815.1 816.9 823.8 828.8 833.0 843.2 845.0 849.4 853.5 855.6 859.6 861.0 864.7 866.7 873.4 878.6 884.6 886.0 889.1 897.1 898.6 921.4 931.7 936.3 939.4 970.1 |
{111} {110}>{113}>{100}>{111} (fig. 6D) {111} {110}>{113}>{100}≈{111} {110}>>{113} {110}>>{113} {110}>{100}>{113}>{111} {110}>{113}>{111}>{100} {111} {110}>{113}>>{100}≈{111} {110}>{113}>{111}>{100} {110}>>{113}>{100} {110}>>{111}>{113}>{100} {110}>{113}>{100}≈{111} {110}>{113}>{111} {110}>{113}>{111}>{100} {111}>{100}≈{110}≈{113} {111} {111} {110}>{113} {110}>{111}>{113}>{100} {110}>{113}>{100} {111} {110} {110} {111} {111} |
aBaseline-corrected peak area used for comparison
|
The PL map in figure 6A, collected using the 455 nm laser, shows the peak area distribution of a band at 508.8 nm. The four growth sectors in the map—{100}, {111}, {110}, and {113}—correspond to those observed in the DiamondView image in figure 4. The peak area for the 508.8 nm band is greatest in the {100} sector, followed by {113} and {111}, respectively. The 508.8 band is undetectable in {110}. The same distribution is observed for the 470.1 nm peak. A band located at 506.1 nm displays a peak area distribution greatest in {113}; the band is also present in {110}, at approximately half the intensity in {113}. The remaining peaks listed in table 2 for the 455 nm laser were most intense in the {111} growth sector.
The relative peak area distributions for the 532 nm laser are summarized in table 2. Figure 6B shows the peak area distribution for the 575 nm PL band caused by NV0. For NV0, the {100}, {111}, {110}, and {113} growth sectors are all visible and correspond with the DiamondView image (again, see figure 4). The NV0 band has the greatest peak area in the {111} growth sector, with the second-greatest intensity found in {100}. The NV0 peak area is an order of magnitude lower in the {113} and {110} sectors than in {100}. The peak area distribution for NV– at 637 nm follows a similar distribution. The peak area intensity of the NV– center is approximately two times greater in the {111} sector than in {100}. The peak area for the GR1 defect at 741 nm was greatest in the {110} and {113} growth sectors; it was not present in the {111} and {100} sectors. Similar results for both NV– and GR1 were observed in PL maps obtained using the 633 nm laser.
The peak area distribution for the NV– center at 637 nm collected using the 633 nm laser is shown in figure 6C. Peak area distributions for the NV– center correlate with the four growth sectors visible in the DiamondView image (figure 4). The peak area intensity for the NV– center is greatest in the {111} growth sector, nearly double its intensity in {100} and two orders of magnitude greater than in the {110} and {113} sectors. At 633 nm excitation, the GR1 peak at 741 nm was present in all four visible growth sectors. However, the GR1 intensity in {111} and {100} was less than half the GR1 intensity found in the {110} sector.
Meanwhile, the 780 nm excitation revealed the most peaks of any laser used for this study (table 2). Figure 6D shows the peak area distribution for the 815.1 nm peak. The four growth sectors seen in the DiamondView image (again, see figure 4) are visible in the PL image for the 815.1 nm peak area (figure 6D). For the 815.1 nm peak, the {110} sector had the greatest peak area intensity, an order of magnitude greater than in {113}. The peak area intensity in both {111} and {100} was an order of magnitude less than in the {113} sector. Many of the peaks observed using the 780 nm laser and reported in table 2 were most intense in the {110} sector.
DISCUSSION
The data in this study and in previous studies (e.g., Johnson et al., 2015; Dieck et al., 2015; Ardon and Loudin, 2016; Johnson and Myagkaya, 2017) show a clear correlation between the growth features observed in PL maps with DiamondView images. The strong correlation demonstrates the effectiveness of PL mapping in resolving the spatial distribution of optical defects in diamond. The visibility of growth features in the PL maps is a function of both the peak’s intensity and its distribution. For example, in this study the 849.4 nm peak is only visible in the {111} growth sector (table 2). Accordingly, the remaining growth sectors cannot be identified using the 849.4 nm peak intensity. Therefore, it is likely that the 849.4 nm peak is preferentially incorporated in the {111} octahedral growth sector. Using similar reasoning, at least 10 peaks from this study are preferentially incorporated in {111} and two peaks in {110} (table 2).
The distribution of the 849.4 nm band also demonstrates the advantage of PL mapping over standard PL. If in a standard PL measurement the instrument operator unknowingly focused on the part of the table corresponding to the {110} sector, then the 849.4 nm peak would not be observed in the spectra. In the same instance, the intensity of the NV– (637 nm) and NV0 (575 nm) centers would be lower and the intensity of GR1 (741 nm) would be higher than if the spectra were collected in the {111} sector. This is important to note because, as stated in Eaton-Magaña and Breeding (2016), determining some diamond color treatments relies on the “presence and absence of particular PL peaks.” Therefore, PL mapping provides scientists with a three-dimensional understanding of the diamond’s optical defects, its PL intensity, and its spatial distribution.
Finally, to assess the effectiveness of the PL mapping technique introduced in this study, we can compare our data to previous studies that have assigned optical defects to particular growth sectors. The concept of optical defects being incorporated in a specific growth sector is well documented (van Enckevort and Lochs, 1988; Collins et al., 1990; Lawson et al., 1996; Welbourn et al., 1996; Zaitsev, 2001). For instance, the NV– center is preferentially incorporated in the {111} octahedral sector in HPHT synthetic diamonds (van Enckevort and Lochs, 1988; Zaitsev, 2001). In agreement with the study presented by van Enckevort and Lochs (1988), our results show that the NV– center is most intense in the {111} octahedral growth sector (table 2). Furthermore, the GR1 intensity was greatest in the {110} growth sector (also in table 2). Comparing the distribution of GR1 to that of the NV– center, we see a nearly opposite distribution of the two defects. This is expected because a vacancy introduced in the {111} growth sector can readily combine with the available N to form the nitrogen vacancies (Burns et al., 1990). The correlation between our PL maps and the results from previous studies validates the effectiveness of using an EMCCD Raman imaging microscope for rapidly collecting PL maps.
CONCLUSIONS
Using a unique HPHT synthetic diamond of mixed type, our high-speed photoluminescence mapping technique produced images that reveal the intensity variations of PL features between growth sectors. The collected PL images, obtained in less than four minutes each, establish the usefulness of this method in understanding the preferential uptake of impurities and optical defects in diamond growth. We conclude that using an EMCCD Raman imaging microscope is an effective and efficient method for determining the distribution of optical centers in diamond and producing PL maps. This technique will allow scientists to rapidly collect and analyze PL maps in order to better understand the nature and formation of optical defects in diamond and evaluate new criteria for identifying treated and synthetic diamonds.