Insight into Relationship between Crystal Structure and Crystal-Field

The results will show that there exists no contradiction. 2. Analysis of garnet crystal structure. Composition change of the above-mentioned garnets h...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCC

Cite This: J. Phys. Chem. C 2018, 122, 3567−3574

Insight into the Relationship between Crystal Structure and CrystalField Splitting of Ce3+ Doped Garnet Compounds Zhen Song,†,‡ Zhiguo Xia,† and Quanlin Liu*,† †

Beijing Key Laboratory for New Energy Materials and Technologies, School of Materials Science and Engineering, University of Science and Technology Beijing, Beijing 100083, China ‡ Departamento de Química, Universidad Autónoma de Madrid, 28049 Madrid, Spain S Supporting Information *

ABSTRACT: The common understanding of the negative relationship between bond lengths and crystal-field splitting (CFS) is renewed by Ce3+ doped garnets in this work. We represent the contradictory relationship between structure data and spectroscopic crystal-field splitting in detail. A satisfactory explanation is given by expressing crystal-field splitting in terms of crystal-field parameters, on the basis of structural data. The results show that not only the bond length, but also the geometrical configuration have influence on the magnitude of crystal-field splitting. Also it is found that the ligand oxygen behaves differently with regard to multiple site substitution in garnet structure. parity-allowed d−f transitions of Ce3+ are first recognized in 1960s in the pioneer work by Blasse et al.,8,9 in which it is used as scintillator. Later, researchers began to pay attention to the absorption and emission spectra of Ce3+ in the garnet series, namely, Y3Al5O12(YAG), Y3Ga5O12(YGG), Gd3Al5O12(GdAG), Gd3Ga5O12(GdGG), Lu3Al5O12 (LuAG), Lu3Ga5O12(LuGG), and their solid solutions.10−13 As will be shown in this work, the cell parameters, volumes of dodecahedron (A) and mean bondlengths are consistent in trend with ionic radius change for different compositions. However, when the crystal-field splitting is taken into account, it gives self-contradictory relationships. For rare-earth elements doped phosphors involving d−f transitions, it is commonly recognized that the crystal-field splitting is closely related to the mean bond lengths between central lanthanides and ligands. In most cases this statement is true, as can be seen from the widely known d5 relation14,15 deduced from the octahedral ligands (Oh symmetry), which is

1. INTRODUCTION Lanthanide aluminum/gallium garnet and its compositional derivatives with Ce3+ as the luminescence activator form a series of highly efficient luminescent materials. Besides the conventional application in cathode-ray display, fluorescent tube, and scintillator, nowadays they are most widely used in the field of white light emitting diodes (WLED). Those phosphors are fabricated together with GaN chips, with the function of converting blue light (from GaN chips) to greenyellow light. The whole device could produce white light due to the combination of blue and green-yellow light. In order to improve the lighting quality, much effort has been paid to tune the spectroscopy. For example, with Ga replacing Al in Y3Al5O12:Ce3+, the emission changes gradually from yellow to green, which could be fabricated together with narrow-band red-emitting phosphors to realize wide gamut of liquid crystal display (LCD).1 Meanwhile, the red shift by substituting Y by Gd could reduce CCT (correlated color temperature) and increase CRI (Color Rendering Index).2 In addition, the garnet phosphors are relatively cheap, chemically stable and easy for mass production. Until now, the garnet series phosphor is the dominant yellow phosphor used in WLED fabrication industry. The garnet series compound has the space group Ia3̅d, and the crystal chemical formula can be expressed as (A)3{B}2[C]3O12. Generally, the larger lanthanides occupy the dodecahedral (A) site, while Al and Ga are accommodated in the octahedral {B} and tetrahedral [C] sites. In some cases, the occupation is rather complex, like the antisite defect caused by Y, Lu entering {B} sites3−5 and the preferential occupation of Ga to [C] sites.6,7 The importance of garnet phosphors arises from the spectroscopy of Ce3+ embedded in garnet lattice. The © 2018 American Chemical Society

Δ = 10Dq ,

D=

35Ze 2 , 4d5

q=

2 4 ⟨r ⟩ 105

In this equation, Δ stands for the crystal-splitting between eg and t2g levels of an d electron. Z is the ligand charge, and d is the distance between central and ligand atoms. This relation serves as a bridge between crystal structure data of the host and luminescent spectroscopy of the doped lanthanides, such as Eu2+, Ce3+. It reveals a simple but intuitive rule: the shorter the Received: December 30, 2017 Revised: January 26, 2018 Published: January 27, 2018 3567

DOI: 10.1021/acs.jpcc.7b12826 J. Phys. Chem. C 2018, 122, 3567−3574

Article

The Journal of Physical Chemistry C Table 1. Crystal Structure Informaiton of Garnet Compounds dodecahedron (A) compd lattice parameter (Å) bond length (Å) mean bond length (Å) polyhedron volume (Å3) ICSD #

YAG 12.002 2.439 2.306 2.372 22.895 4114524

YGG 12.280 2.428 2.338 2.383 23.052 2385225

GdAG 12.114 2.522 2.353 2.437 24.895 19218426

bond length, the larger the crystal-field splitting. Researchers have applied this rule to explain the photoluminescent excitation (PLE) peak shifts in many solid solutions. For example, in M2MgSi2O7:Eu2+ (M = Ca,Sr,Ba) series,16 a red shift toward longer wavelength of PLE peaks is observed when the composition varies from Ba to Ca, with Ca−O bond lengths shorter than those of Ba−O. However, the simple relationship between bond length and crystal-field splitting is challenged by the lanthanide aluminum/ gallium garnet, as mentioned above. For example, from YAG to YGG, the cell is expanding due to the larger Ga, and the crystalfield splitting is decreasing,17 which is in accordance with expectation. On the other hand, the shrinkage of cell is observed from YAG to LuAG, while the crystal-field splitting is surprisingly reduced.18,19 Some researchers name those compounds as reverse garnets.14,20 However, what is confusing is that we can observe the normal and reverse phenomena in the same compound, just as the example of YAG. This selfcontradictory phenomenon greatly challenges the relationship between bond lengths and crystal-field splitting. The main interest of this work is focused on the explanation of normal/reverse garnet phenomena. One important attempt has been done by Seijo and Barandiarán using AIMPs (Ab Initio Model Potentials) methods,21 which could yield all five nondegenerate 5d levels. They find the results in good accordance with the normal/reverse phenomena. In addition, as pointed out by Dorenbos,22 the type of polyhedron surrounding Ce3+ also influences the crystal-field splitting, and Wu et al. explore the correlation between the dodecahedral distortion and crystal-field splitting.20 In this paper, we provide a more direct approach linking the garnet structure and crystalfield splitting by point-charge electrostatic model (PCEM) and simple perturbation calculation. The results will show that there exists no contradiction.

GdGG 12.379 2.527 2.354 2.440 25.026 19218126

LuAG 11.906 2.383 2.276 2.330 21.621 2384625

LuGG 12.188 2.393 2.303 2.348 22.027 2385025

Figure 1. Mean bond lengths and dodecahedral volumes of garnet compounds as a function of lattice parameter.

composition varies from Al to Ga. The above analysis show that it is the ionic size of different elements constituting garnets that determines the bond lengths between (A) site and ligands. In the next part, we will check whether this statement holds for crystal-field splitting.

3. ENERGY LEVEL SCHEME OF Ce3+ IN GARNET STRUCTURES AND CRYSTAL-FIELD SPLITTING For rare-earth elements doped phosphors, the dopant lanthanides generally occupy some specific sites. In garnet structure, Ce3+ is accommodated in the dodecahedral (A) site, with point symmetry D2. Within the eight ligand-coordinated crystal-field potential, the 5-fold degenerate 5d levels will be totally split. The energy level scheme has been fully discussed,27−29 and it is widely accepted that the degenerated energy levels are first split by Td crystal-field potential, with 2E (2-fold, low-lying) and 2T2 (3-fold, higher) levels, and then further split by D2 crystal-field potential. Researches are particularly interested in the two lowest 5d levels, because absorption transitions from the 4f ground levels to the two lowest 5d levels result in two PLE (photoluminescent excitation) peaks located in the visible-blue range. In this work, we treat the energy difference between those two peaks which is abbreviated as E12, as the representation of crystal-field splitting based on two reasons. First, the higher 5d levels may be immersed in the conduction band or located beyond the ultraviolet region,9 which makes them difficult to detect. Second, as pointed by Dorenbos,30 the location of PLE peaks are determined by both centroid shift and crystal-field splitting. A large survey of literature on Eu2+ and Ce3+ doped sulfide,

2. ANALYSIS OF GARNET CRYSTAL STRUCTURE Composition change of the above-mentioned garnets has no effect on the space group. Therefore, the structure change induced by dopant incorporation could be examined by checking lattice parameters and bond lengths. The detailed structural data is listed in Table 1. For the (A) site, there exist two kinds of bond lengths, therefore the mean bond lengths are taken into account. The relationship between composition, lattice parameter, bond length, and volume of dodecahedra is displayed in Figure 1. It can be easily seen that the trends are in accordance with the ionic radius under 8 coordinated ligands, which is 0.977 Å for Lu3+, 1.019 Å for Y3+, and 0.938 Å for Gd3+.23 The solid lines connecting the whole Al or Ga series show that the cell expands together with the bond length and dodecahedra as the sequence Lu < Y < Gd. Meanwhile, those garnets belonging to the same rare-earth element are connected by dashed lines, which also give increasing trend as the 3568

DOI: 10.1021/acs.jpcc.7b12826 J. Phys. Chem. C 2018, 122, 3567−3574

Article

The Journal of Physical Chemistry C oxide, fluoride and nitride compounds shows that the centroid shift is connected to composition, while the crystal-field splitting is mainly affected by geometrical configuration.22,31 Therefore, the energy difference, other than the energy location is discussed in the following sections. The E12 value could be easily determined by picking up the two peaks from PLE spectra in energy scale. The data is compiled in Table 2. Table 2. Data of Crystal-Field Splitting of Ce Garnets wavenumber (cm−1) YAG

YGG GdAG GdGG LuAG

29239 29239

21739 21929

29400 28571 29762 28735 28700 28233

22000 23810 21253 23364 23500 21138

3+

GdGG. However, the normal-garnet phenomenon also exists, as the cases of Figure 2c and d, in which E12 decreases with larger Ga substituting Al. This contradiction puts up great challenge to the simple relationship between bond length and crystal-field splitting, and we will show in the next section that the origin arises from the limitation of d5 relation.

4. CRYSTAL-FIELD ANALYSIS OF Ce3+ IN GARNET STRUCTURE Crystal-field theory has been widely applied to rare-earth doped crystals.38,39 By fitting crystal-field parameters to spectroscopic measurements, the energy levels of doped rare-earth elements could be deduced.40 The simple d5 relation is deduced under Oh symmetry, i.e., center of 8-coordinated octahedron. The symmetry elements in Oh point group will eliminate all the second rank crystal-field parameters in the crystal-field potential.41,42 As a result, for d electron the Oh crystal-field potential could be expressed using only one-fourth rank crystalfield parameter, i.e., B40. If expressed by the point-charge electrostatic model (PCEM), B40 has the form of d5 in the denominator. However, when the site symmetry degrades, the symmetry limitations on the second crystal-field parameters are broken, and they become nonzero, which will introduce d3 together with d5. As mentioned above, for garnet structure, the dopant Ce3+ is accommodated in the 8-coordinated (A) site, with point symmetry D2. The crystal-field potential is expressed as41

Doped

crystal-field splitting

ref

7500 7310 8360 7400 4762 8509 5371 5200 7100

Bachmann et al.33 Chen et al.2 Ogieglo et al.34 Blasse et al.9 Ueda35 Chen et al.2 Ogieglo et al.36 Kaminska et al.37 Ogieglo et al.34

Ogieglo et al. reported a crystal-field splitting larger than 8000 cm−1, which is far away from values of other references and is beyond consideration. The collection lacks data on YGG, with the possible reason that the luminescence of Ce3+ is quenched due to the down-shift of conduction band.32 Investigation into the relationship between bond lengths and crystal-field splitting of Ce3+ doped garnets, as shown in Figure 2, produces self-contradictory results. For the aluminum-garnet series, E12 is increasing together with the cell expansion in the sequence of Lu < Y < Gd, as shown in Figure 2a. This reversegarnet phenomenon is also observed in Figure 2b for YGG and

V (D2 ) = [B20 U02 + B22 (U22 + U −22)]⟨l C 2 l⟩ + [B40 U04 + B42 (U24 + U −42) + B44 (U44 + U −44)]⟨l C 4 l⟩

Therefore, it is expected that the second rank crystal-field parameters will be included in the analytical expression of split

Figure 2. Relationship between crystal-field splitting and bond lengths compiled from references. 3569

DOI: 10.1021/acs.jpcc.7b12826 J. Phys. Chem. C 2018, 122, 3567−3574

Article

3570

⎡ ⎛ 3.01⟨r 2⟩ ⎞⎤2 =⎢Ze 2⎜ 4 ⎟⎥ ⎣ ⎝− 1.33⟨r ⟩⎠⎦ ⎡ ⎛1.24⟨r 2⟩ ⎞⎤2 =⎢Ze 2⎜ 4 ⎟⎥ ⎣ ⎝− 1.50⟨r ⟩⎠⎦

⎛ 239.86B20 ⎞2 ⎟ ≈⎜+ 78.23B ⎝ 40 ⎠ ⎛ 240.00B20 ⎞2 ⎟ ≈⎜+ 74.59B ⎝ 40 ⎠

⎡ ⎛ 3.11⟨r 2⟩ ⎞⎤2 =⎢Ze 2⎜ 4 ⎟⎥ ⎣ ⎝− 1.41⟨r ⟩⎠⎦ ⎡ ⎛ 2.64⟨r 2⟩ ⎞⎤2 =⎢Ze 2⎜ 4 ⎟⎥ ⎣ ⎝− 1.57⟨r ⟩⎠⎦

B20 B22 B40 B42 B44 (E12)2

⎡ ⎛ 4.38⟨r 2⟩ ⎞⎤2 =⎢Ze 2⎜ 4 ⎟⎥ ⎣ ⎝− 1.20⟨r ⟩⎠⎦

⎛ 240.00B20 ⎞2 ⎟ ≈⎜+ 64.73B ⎝ 40 ⎠ ⎛ 240.04B20 ⎞2 ⎟ ≈⎜+ 70.68B ⎝ 40 ⎠ ⎛ 239.97B20 ⎞2 ⎟ ≈⎜+ 71.34B ⎝ 40 ⎠

⎡ ⎛ 0.996⟨r 2⟩ ⎞⎤2 =⎢Ze 2⎜ 4 ⎟⎥ ⎣ ⎝− 1.67⟨r ⟩⎠⎦

−0.004151e ⟨r ⟩Z 0.005085e2⟨r2⟩Z 0.02140e2⟨r4⟩Z 0.01312e2⟨r4⟩Z 0.002870e2⟨r4⟩Z 57600.00B220 + 37499.39B20B40 + 6104.26B240 −0.01824e ⟨r ⟩Z 0.02234e2⟨r2⟩Z 0.01728e2⟨r4⟩Z 0.00872e2⟨r4⟩Z 0.003737e2⟨r4⟩Z 57560.00B220 + 31065.01B20B40 + 4190.97B240 −0.01294e ⟨r ⟩Z 0.01585e2⟨r2⟩Z 0.02000e2⟨r4⟩Z 0.01106e2⟨r4⟩Z 0.003600e2⟨r4⟩Z 57618.0B220 + 33926.3B20B40 + 4995.75B240 −0.01100e ⟨r ⟩Z 0.01347e2⟨r2⟩Z 0.02201e2⟨r4⟩Z 0.01228e2⟨r4⟩Z 0.003865e2⟨r4⟩Z 57584.7B220 + 34232.5B20B40 + 5089.22B240

⎛ 240.00B20 ⎞2 ⎟ ≈⎜+ 78.13B ⎝ 40 ⎠

−0.005184e ⟨r ⟩Z 0.006349e2⟨r2⟩Z 0.02005e2⟨r4⟩Z 0.01172e2⟨r4⟩Z 0.003125e2⟨r4⟩Z 57600.0B220 + 35797.0B20B40 + 5563.02B240 2 2 2

2

2

2

GdAG YAG LuAG

2 2

Table 3. Expressons of Crystal-Field Parameters and Crystal-Field Splittings of Garnet Compounds

2

LuGG

2

YGG

2

GdGG

energy levels. In this case, the simple d5 relationship takes no more effect. We will follow a reported routine17 to treat the crystal-field splitting of Ce3+ doped garnets. In brief, the perturbation matrix is constructed according to the D2 crystal-field potential, and then the eigenvalues are obtained by matrix diagonalization. The important point lies in how to treat crystal-field potential and spin−orbit coupling. Unlike transition-metal elements, for rare-earth elements the spin−orbit coupling is significant. If spin−orbit coupling is treated prior to crystal-field potential, the energy levels of 5d1 electron are first split into 2D5/2 and 2 D3/2 terms, which will be further split under the crystal-field potential. However, this treatment is inconsistent with the real case, because the energy difference between the two 2D terms is 5ξ/2 (ξ, the spin−orbit coupling constant, has a value of 991 cm−1 for 5d electron of Ce3+),43 much smaller than the energy level separation. If the energy levels are first split by crystal-field potential and after that spin−orbit coupling is considered, the final split energy will be only related to ξ, which means the energy level separation keeps unchanged among different compounds. This is also inconsistent with the real case. In fact, in the correct method the spin−orbit coupling and crystal-field potential should be considered simultaneously, but this treatment will lead to solving a 5-power secular equation during the diagonalization process, which seldom gives analytical results. A more realistic routine is to assume cubic field ≫ the spin−orbit interaction ≫ D2 crystal-field,28 and by this method we have shown the decrease of crystal-field splitting from YAG to YGG, as explained in Appendix of ref 17. In this work, the other garnets are treated in the same way, and the results are listed in Table 3. Interpretation of those results provide convincing explanation for the above-mentioned selfcontradictory phenomena. For the aluminum and gallium series, it is clear that E12 increases as Lu < Y < Gd. For LnAG/ LnGG paris, i.e., LuAG/LuGG, YAG/YGG and GdAG/GdGG, E12 drops as Ga substitutes Al. Those results are in accordance with Figure 2 and the contradiction disappears. It should be noted that the sequence is related to the values of ⟨r2⟩ and ⟨r4⟩, for which a coarse estimation is provided in the Supporting Information. A close examination of Table 3 shows that B20 changes magnitude sharply among different compounds, while the change of B40 is relatively small. For both Al and Ga garnet series, the magnitude of B20 increases as the sequence Lu < Y < Gd. The same trend also holds for every Al/Ga pair, just in accordance with E12. Meanwhile, when Ga completely substitutes Al, the magnitude of B20 drops nearly half. To get a more vivid understanding, the contribution of every ligand to B20 of all the compounds is summarized in the Supporting Information. It is easily seen that for every compound only ligand 4 and 5 have positive contribution. These two ligands are related by one 2-fold axis. They share the same bond lengths and have the same contribution to B20. In addition, ligand 4 and 5 are more close to the polar axis (z axis) than the other ligands. Based on these facts, we consider the sum of contributions of ligand 4 and 5 as a whole, whereas the sum of contributions from the other ligands is considered as the counterpart for comparison. They are listed in Table 4, and plotted in Figure 3. For both positive and negative contributions, the magnitudes are in accordance with the bond lengths. Figure 3a and b shows that for both Al and Ga series, Lu garnet has the largest magnitudes of the two items, and Gd garnet has the smallest. This is reasonable because B20

−0.01253e ⟨r2⟩Z 0.01534e2⟨r2⟩Z 0.01694e2⟨r4⟩Z 0.01040e2⟨r4⟩Z 0.002259e2⟨r4⟩Z 57533.9B220 + 37525.9B20B40 + 6119.87B240

The Journal of Physical Chemistry C

DOI: 10.1021/acs.jpcc.7b12826 J. Phys. Chem. C 2018, 122, 3567−3574

Article

The Journal of Physical Chemistry C Table 4. Polar Angle (θ), Positive and Negative B20 Contributions of Al and Ga Series Garnet garnet compound

B20/e2⟨r2⟩Z

polar angle of ligand 5 (deg) sum of 4 and 5 (positive contribution) sum of the other (negative contribution) total

LuAG

YAG

GdAG

LuGG

YGG

GdGG

11.29 0.139 −0.150 −0.011

11.24 0.131 −0.144 −0.013

11.59 0.117 −0.135 −0.018

10.54 0.139 −0.143 −0.0041

10.79 0.132 −0.138 −0.0052

10.24 0.118 −0.130 −0.0125

Figure 3. Absolute value of B20 coefficient for: sum of ligand 4 and 5, blue triangle; sum of the other ligands, red square; total sum of all the ligands, red circle; polar angle of ligand 5, pink diamond. Panels (a) and (b) represent the Al and Ga garnet series, respectively, while panels (c) and (d) represent LnAG/LnGG paris.

Figure 4. (a) Coordinated ligands of dodecahedral (A) site. The crystallographically equivalent eight oxygen atoms could be grouped into two categories according to the bond lenths. They are rendered by different colors: (b) Perspective view of the distorted cube of the dedecahedral ligands. (c) Schematic of the adjacent dodecahedron, octahedron and tetrahedron with respect to the oxygen ligands. (d) Top view of the distorted cube of the dodecahedral ligands. Solid lines indicate shared edges with another dodecahedron, dashed lines octahedron and dotted lines tetrahedron. (e) Ligand movement pattern for LuAG to YAG. (f) Ligand movement pattern for LuAG to LuGG. The arrows indicate the moving direction of ligands when the garnet compound change. The magnitude of the arrow has been magnified 30 times for a better visualization. Plotted by VESTA.44

3571

DOI: 10.1021/acs.jpcc.7b12826 J. Phys. Chem. C 2018, 122, 3567−3574

Article

The Journal of Physical Chemistry C is reverse proportional to r3, which could be seen from the PCEM expression B20 =

e 2⟨r 2⟩Z(1 + 3 cos 2θ ) 4r 3

However, the difference between the positive and negative parts has reverse trend in magnitude. For Gd garnet, although the magnitudes of positive and negative contributions are the smallest, the difference is the largest, as shown by the blue triangle line in Figure 3a and b. As a result, Gd garnet has the largest B20 among Al or Ga series. For the LnAG/LnGG pairs, Figure 3c shows that the magnitude of the negative contribution also decreases as cell expands. However, the positive contribution has a rising trend for all the Al/Ga pairs, despite the fact that bond lengths are stretching from Al to Ga. This arises from the effect of polar angle (θ) which appears in the nominator of B20 expression. When θ < 45°, B20 is positive and has a larger magnitude at smaller polar angles. From Figure 3d, it is known that all the Ga garnet series have smaller polar angles than Al series. In the case of LnAG/LnGG pairs, the polar angle dominates the magnitude of positive contribution to B20. As a result, with smaller negative contribution and larger positive contribtuion, the magnitude of B20 shrinks from Al to Ga, which leads to the decrease of E12. It should be noted that the results of crystal-field contributions of every ligand depend on the axis choice. Under different coordinate-systems the values may change, but the overall result remains the same due to physical equality.

Figure 5. Cube distortion (d23/d21) for garnet compounds. Red squares belong to the Al series, while blue circles belong to Ga series. Data is extracted by VESTA.

larger Ga substituting Al results to a decrease in d23/d21. It is evident that the ratio between d23 and d21 has the same trend with crystal-field splitting (Figure 2) and crystal-field parameter (Figure 3). However, to the best of our knowledge, it is difficult to connect crystal-field splitting with cube distortion. As will be shown in this section, the similarity in the varying trend of crystal-field splitting and cube distortion arises from the movements of dodecahedral ligands through garnet compounds. It is known from the above section that both incorporation of larger size lanthanides and substituting Al for Ga have the effect of expanding the cell and dodecahedron volume. However, there exist subtle differences between how the oxygen ligands respond to different site substitutions. A detailed study of the dodecahedral ligand coordinates of garnet compounds reveals that the ligand movement has two distinct patterns with regards to site substitutions. In Figure 4e and f, the evolvements from LuAG to YAG and LuAG to LuGG are selected to elucidate the two ligand movement patterns. The arrows indicate the ligand movement between different garnet compounds. For the (A) site substitution, the ligand movement shown in Figure 4e indicates a much larger elongation of distance between O2 and O3 than that between O1 and O2. This pattern is in accordance with the dodecahedron inflation from LuAG to YAG, considering the fact that edge 2−1 is shared by tetrahedron, while edge 2−3 is shared by another dodecahedron. Therefore, d23/d21 is increasing through Lu−Y−Gd among Al and Ga garnet series. However, for the {B} and [C] sites substitutions, which are referred to Ga substitution for Al from LuAG to LuGG in Figure 4e, they share a quite different mode. The ligand movement pattern indicates a further separation between O1 and O2, but an approach between O2 and O3. This means inflation of octahedron and tetrahedron simultaneously squeezes dodecahedron. As a result, we expect a shorter d23 and longer d12, then finally a smaller d23/d21. Meanwhile, it is noted that ligand 4 and 5 approach closer to z axis in the second movement pattern, which accounts for the smaller polar angles of Ga series garnet in Figure 3d. From this point of view, it is known that the ligand oxygen atoms respond in different ways regarding to different site substitutions. Although the whole effect is the same, i.e., larger cell and longer bond lengths, the two ligand movement pattern

5. STRUCTURE−PROPERTY CORRELATIONS IN GARNET STRUCTURE The eight coordinated ligands of the dodecahedral (A) site could be categorized into two groups, as colored in red and blue in Figure 4a. Although they are crystallographically equivalent, ligands 1, 2, 6, 7 share the same bond length, while ligands 3, 4, 5, 8 share the other bond length. The dodecahedron could be thought as a distorted cube, as shown in Figure 4b, with ligand 1−4−7−8 constituting the bottom face, and ligand 2−3−6−5 the top face. Under the point symmetry D2, some of the edges share the same distances, which are d23 = d56 = d18 = d47, d14 = d78 = d25 = d36, d21 = d67 and d58 = d34. In addition, the equivalent edges have the same connectivity with polyhedrons, as illustrated in Figure 4c. In the top and bottom face, the four edges are shared by octahedron and another dodecahedron successively. For the four edges connecting the top and bottom faces, edge 2−1 and 6−7 are shared by tetrahedrons, while edge 5−8 and 3−4 have no share with other polyhedrons. In Figure 4d, the distortion could be visualized more clearly from the top view of the cube. In addition, distortion of the cube has been adopted by Euler and Bruce to investigate the crystal chemistry of garnets.25 Wu et al. discovered the ratio between d88 and d18 could be used to characterize the crystal-field splitting of Ce3+ in garnets,20 while d88 and d18 represent the O−O distances shared by adjacent dodecahedron and tetrahedron, respectively. In this work, d88 and d18 are equivalent to the distances between ligand 2, 3 and 2, 1 as shown in Figure 4c. Consequently, in Figure 4d, the distortion of the cube could be represented by d23/d21, which are plotted for both Al and Ga garnet series in Figure 5. For Al/ Ga series, incorporation of larger cations into dodecahedral (A) site increases the ratio. On the contrary, for LnAG/LnGG pairs 3572

DOI: 10.1021/acs.jpcc.7b12826 J. Phys. Chem. C 2018, 122, 3567−3574

Article

The Journal of Physical Chemistry C Notes

leads to the reverse trend of crystal-field splitting and cube distortions. The normal and reverse phenomena on crystal-field splitting phenomena of garnet structures could be understood from the above treatment. However, it should be noted that PCEM is a coarse model, because the ligands are regarded as negative point-charges in this model, which deviates from the real case. Moreover, the bond lengths are extracted from pure garnet structures without dopant Ce3+, and it is known that incorporation of larger Ce3+ will cause cell expansion, which means the real bond lengths should be larger than those used in this work. Nevertheless, PCEM could be regarded as the first approximation, and it could be applied to explain the crystalfield splitting in Ce3+ doped garnets.

The authors declare no competing financial interest.



REFERENCES

(1) Cui, D. P.; Song, Z.; Xia, Z. G.; Liu, Q. L. Luminescence Tuning, Thermal Quenching, and Electronic Structure of Narrow-Band RedEmitting Nitride Phosphors. Inorg. Chem. 2017, 56, 11837−11844. (2) Chen, L.; Chen, X.; Liu, F.; Chen, H.; Wang, H.; Zhao, E.; Jiang, Y.; Chan, T. S.; Wang, C. H.; Zhang, W.; et al. Charge Deformation and Orbital Hybridization: Intrinsic Mechanisms on Tunable Chromaticity of Y3Al5O12:Ce3+ Luminescence by Doping Gd3+ for Warm White LEDs. Sci. Rep. 2015, 5, 11514. (3) Kanai, T.; Satoh, M.; Miura, I. Characteristics of a Nonstoichiometric Gd3+Δ(Al,Ga)5−ΔO12:Ce Garnet Scintillator. J. Am. Ceram. Soc. 2008, 91, 456−462. (4) Babin, V.; Blazek, K.; Krasnikov, A.; Nejezchleb, K.; Nikl, M.; Savikhina, T.; Zazubovich, S. Luminescence of Undoped LuAG and YAG Crystals. Phys. Status Solidi C 2005, 2, 97−100. (5) Muñoz-García, A. B.; Barandiarán, Z.; Seijo, L. Antisite Defects in Ce-Doped YAG (Y3Al5O12): First-Principles Study on Structures and 4f−5d Transitions. J. Mater. Chem. 2012, 22, 19888−19897. (6) Marezio, M.; Remeika, J. P.; Dernier, P. D. Cation Distribution in Y3Al5−CGacO12 Garnet. Acta Crystallogr., Sect. B: Struct. Crystallogr. Cryst. Chem. 1968, 24, 1670−1674. (7) Laguta, V.; Zorenko, Y.; Gorbenko, V.; Iskaliyeva, A.; Zagorodniy, Y.; Sidletskiy, O.; Bilski, P.; Twardak, A.; Nikl, M. Aluminum and Gallium Substitution in Yttrium and Lutetium Aluminum−Gallium Garnets: Investigation by Single-Crystal NMR and TSL Methods. J. Phys. Chem. C 2016, 120, 24400−24408. (8) Blasse, G.; Bril, A. A New Phosphor for Flying-Spot Cathode-Ray Tubes for Color Television: Yellow-Emitting Y3Al5O12-Ce3+. Appl. Phys. Lett. 1967, 11, 53−55. (9) Blasse, G.; Bril, A. Investigation of Some Ce3+ Activated Phosphors. J. Chem. Phys. 1967, 47, 5139−5145. (10) Robbins, D. J. The Effects of Crystal Field and Temperature on the Photoluminescence Excitation Efficiency of Ce3+ in YAG. J. Electrochem. Soc. 1979, 126, 1550−1555. (11) Robbins, D. J.; Cockayne, B.; Glasper, J. L.; Lent, B. The Temperature Dependence of Rare-Earth Activated Garnet Phosphors. J. Electrochem. Soc. 1979, 126, 1213−1220. (12) Liu, X. R.; Wang, X. J.; Shun, W. F. Luminescence Properties of the Ce3+ Ion in Yttrium Gallium Garnet. Phys. Status. Solidi. A 1987, 101, K161−K165. (13) Tien, T. Y.; Gibbons, E. F.; DeLosh, R. G.; Zacmanidis, P. J.; Smith, D. E.; Stadler, H. L. Ce3+ Activated Y3Al5O12 and Some of Its Solid Solutions. J. Electrochem. Soc. 1973, 120, 278−281. (14) George, N. C.; Denault, K. A.; Seshadri, R. Phosphors for SolidState White Lighting. Annu. Rev. Mater. Res. 2013, 43, 481−501. (15) Dunn, T. M.; McClure, D. S.; Pearson, R. G. Some Aspects of Crystal Field Theory; Harper & Row: New York, 1965. (16) Kim, T.; Kim, Y.; Kang, S. Luminescence Properties of Eu2+ in M2MgSi2O7 (M= Ca, Sr, and Ba). Appl. Phys. B: Lasers Opt. 2012, 106, 1009−1013. (17) Song, Z.; Wang, Z. Z.; He, L. Z.; Zhang, R. J.; Liu, X. L.; Xia, Z. G.; Geng, W. T.; Liu, Q. L. After-Glow, Luminescent Thermal Quenching, and Energy Band Structure of Ce-Doped Yttrium Aluminum-Gallium Garnets. J. Lumin. 2017, 192, 1278−1287. (18) Ahn, W.; Kim, Y. J. Substitutional Solubility Limit for Ce3+ Ions in Lu3Al5O12:xCe3+ and Its Effect on Photoluminescence. Ceram. Int. 2017, 43, S412−S416.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.7b12826. Evaluation of ⟨r2⟩ and ⟨r4⟩ and contribution of every ligand to B20 of all the garnet compounds (PDF)



ACKNOWLEDGMENTS

This work is supported by National Natural Science Foundation of China (No. 51602019 and 51672027) and Fundamental Research Funds for the Central Universities (FRF-TP-17-005A2). Z.S. thanks China Scholarship Council (CSC) for a visiting scholarship support in UAM.

6. CONCLUSIONS This work is devoted to explain the crystal-field splitting in Ce3+ doped garnet compounds. Generally, incorporation of larger size atoms will expand the cell and prolong the bond lengths. This includes the successive substitution of Lu, Y, Gd for dodecahedral (A) sites, and Ga substituting Al. The crystal-field splittings of garnet compounds show reverse trends. For Al or Ga garnet series, E12 increases in the sequence LuAG < YAG < GdAG, and YGG < GdGG. For LnAG/LnGG pairs, E12 decreases in the sequence YAG > YGG and GdAG > GdGG. The reverse trends are contradictory to the common understanding that crystal-field splitting has a negative correlation with bond lengths. With the help of PCEM and energy level scheme of Ce3+ in garnet structure, we obtain the expressions for E12, which is in accordance with the reverse trends. More detailed investigation reveals that E12 is related to the difference between the positive and negative contributions of the eight ligands around (A) site. Although the magnitudes of both contributions are mainly determined by bond lengths, their differences exhibit different correlations with bond lengths. In addition, the geometrical configuration has a profound effect on the magnitudes of difference in the cases of LnAG/LnGG pairs. This arises from the two different ligand movement patterns among garnet compounds, which also accounts for the cube distortion of dodecahedral ligands. The method provided in this work could be applied to other luminescent materials to elucidate the relationship between crystal-field splitting and crystal structure.





AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Phone/Fax: 86-10-62334705. ORCID

Zhiguo Xia: 0000-0002-9670-3223 Quanlin Liu: 0000-0003-3533-7140 Author Contributions

The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. 3573

DOI: 10.1021/acs.jpcc.7b12826 J. Phys. Chem. C 2018, 122, 3567−3574

Article

The Journal of Physical Chemistry C

cence Properties for NdxY1‑XAl3(BO3)4 Nonlinear Laser Crystal. J. Mater. Chem. C 2017, 5, 7174−7181. (40) Ju, M.; Lu, C.; Yeung, Y.; Kuang, X.; Wang, J.; Zhu, Y. Structural Evolutions and Crystal Field Characterizations of Tm-Doped YAlO3: New Theoretical Insights. ACS Appl. Mater. Interfaces 2016, 8, 30422− 30429. (41) Görller-Walrand, C.; Binnemans, K. Rationalization of CrystalField Parametrization. In Handbook on the Physics and Chemistry of Rare Earths; Elsevier: 1996; Vol. 23, pp 121−283. (42) Liu, G.; Jacquier, B. Spectroscopic Properties of Rare Earths in Optical Materials; Springer Science & Business Media, 2006; Vol. 83. (43) Tanner, P. A.; Fu, L.; Ning, L.; Cheng, B.-M.; Brik, M. G. Soft Synthesis and Vacuum Ultraviolet Spectra of YAG: Ce3+ Nanocrystals: Reassignment of Ce3+ Energy Levels. J. Phys.: Condens. Matter 2007, 19, 216213. (44) Momma, K.; Izumi, F. Vesta 3 for Three-Dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 1272−1276.

(19) Dorenbos, P. Electronic Structure and Optical Properties of the Lanthanide Activated Re3(Al1−XGax)5O12 (Re = Gd, Y, Lu) Garnet Compounds. J. Lumin. 2013, 134, 310−318. (20) Wu, J. L.; Gundiah, G.; Cheetham, A. K. Structure−Property Correlations in Ce-Doped Garnet Phosphors for Use in Solid State Lighting. Chem. Phys. Lett. 2007, 441, 250−254. (21) Seijo, L.; Barandiarán, Z. Large Splittings of the 4f Shell of Ce3+ in Garnets. Phys. Chem. Chem. Phys. 2014, 16, 3830−3834. (22) Dorenbos, P. Crystal Field Splitting of Lanthanide 4fn−15dLevels in Inorganic Compounds. J. Alloys Compd. 2002, 341, 156−159. (23) Shannon, R. Revised Effective Ionic Radii and Systematic Studies of Interatomic Distances in Halides and Chalcogenides. Acta Crystallogr., Sect. A: Cryst. Phys., Diffr., Theor. Gen. Crystallogr. 1976, 32, 751−767. (24) Bagdasarov, K. S.; Bolotina, N.; Kalinin, V.; Karyagin, V.; Kuz’min, B.; Muradyan, L.; Ryadnov, S.; Uyukin, E.; Chernaya, T.; Fedorov, E. Photoinduced Effects and Real Structure of Crystals of Yttrium Aluminum Garnet. Soviet physics. Crystallography 1991, 36, 398−405. (25) Euler, F.; Bruce, J. A. Oxygen Coordinates of Compounds with Garnet Structure. Acta Crystallogr. 1965, 19, 971−978. (26) Hamilton, A. S.; Lampronti, G.; Rowley, S.; Dutton, S. Enhancement of the Magnetocaloric Effect Driven by Changes in the Crystal Structure of Al-Doped GGG, Gd3Ga5−XAlxO12 (0 ≤ X ≤ 5). J. Phys.: Condens. Matter 2014, 26, 116001. (27) Li, J. G.; Sakka, Y. Recent Progress in Advanced Optical Materials Based on Gadolinium Aluminate Garnet (Gd3Al5O12). Sci. Technol. Adv. Mater. 2015, 16, 014902. (28) Hoshina, T.; Kuboniwa, S. 4f-5d Transition of Tb3+ and Ce3+ in MPO4(M = Sc, Y and Lu). J. Phys. Soc. Jpn. 1971, 31, 828−840. (29) Song, Z.; Liu, X. L.; He, L. Z.; Liu, Q. L. Correlation between the Energy Level Structure of Cerium-Doped Yttrium Aluminum Garnet and Luminescent Behavior at Varying Temperatures. Mater. Res. Express 2016, 3, 055501. (30) Dorenbos, P. The 5d Level Positions of the Trivalent Lanthanides in Inorganic Compounds. J. Lumin. 2000, 91, 155−176. (31) Wang, S. X.; Song, Z.; Kong, Y. W.; Xia, Z. G.; Liu, Q. L. Crystal Field Splitting of 4fn−15d-Levels of Ce3+ and Eu2+ in Nitride Compounds. J. Lumin. 2018, 194, 461−466. (32) Fasoli, M.; Vedda, A.; Nikl, M.; Jiang, C.; Uberuaga, B. P.; Andersson, D. A.; McClellan, K. J.; Stanek, C. R. Band-Gap Engineering for Removing Shallow Traps in Rare-Earth Lu3Al5O12 Garnet Scintillators Using Ga3+Doping. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, 081102. (33) Bachmann, V.; Ronda, C.; Meijerink, A. Temperature Quenching of Yellow Ce3+ Luminescence in YAG: Ce. Chem. Mater. 2009, 21, 2077−2084. (34) Ogieglo, J. M.; Zych, A.; Ivanovskikh, K. V.; Justel, T.; Ronda, C. R.; Meijerink, A. Luminescence and Energy Transfer in Lu3Al5O12 Scintillators Co-Doped with Ce3+ and Tb3+. J. Phys. Chem. A 2012, 116, 8464−8474. (35) Ueda, J. Analysis of Optoelectronic Properties and Development of New Persistent Phosphor in Ce3+-Doped Garnet Ceramics. J. Ceram. Soc. Jpn. 2015, 123, 1059−1064. (36) Ogieglo, J. M.; Katelnikovas, A.; Zych, A.; Justel, T.; Meijerink, A.; Ronda, C. R. Luminescence and Luminescence Quenching in Gd3(Ga,Al)5O12 Scintillators Doped with Ce3+. J. Phys. Chem. A 2013, 117, 2479−84. (37) Kaminska, A.; Duzynska, A.; Berkowski, M.; Trushkin, S.; Suchocki, A. Pressure-Induced Luminescence of Cerium-Doped Gadolinium Gallium Garnet Crystal. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, 155111. (38) Ju, M.; Kuang, X.; Lu, C.; Li, H.; Wang, J.; Zhang, C.; Zhu, Y.; Yeung, Y. Determination of the Microstructure, Energy Levels and Magnetic Dipole Transition Mechanism for Tm3+ Doped Yttrium Aluminum Borate. J. Mater. Chem. C 2016, 4, 1988−1995. (39) Ju, M.; Sun, G.-l.; Kuang, X.-y.; Lu, C.; Zhu, Y.-s.; Yeung, Y.-y. Theoretical Investigation of the Electronic Structure and Lumines3574

DOI: 10.1021/acs.jpcc.7b12826 J. Phys. Chem. C 2018, 122, 3567−3574