Growth, Structure, Thermal Properties and Spectroscopic Characteristics of Nd3+-Doped KGdP4O12 Crystal

A single crystal of Nd3+-doped KGdP4O12 was successfully grown with the top-seeded solution growth and slow cooling (TSSG−SC) technique. It crystallizes in space group C2/c with cell parameters a = 7.812(2) Å, b = 12.307(3) Å, c = 10.474(2) Å, β = 110.84(3)° and Z = 4. The IR and Raman spectra also indicated that the phosphoric polyhedra of Nd:KGdP4O12 has a cyclic symmetry. The chemical composition of the crystal was analyzed and the distribution coefficient of Nd3+ was calculated. The crystal morphology of KGdP4O12 was identified using X-ray diffraction. The compound has good thermal stability to 920°C. Its specific heat and thermal conductivity were determined for potential applications. The spectral properties of Nd:KGdP4O12 indicates that it exhibits broad absorption and emission bands, which are attributed to low symmetry of the crystal. The broad absorption band around 798 nm has a full-width at half-maximum (FWHM) of 14.8 nm and is suitable for AlGaAs laser diode pumping. Moreover, 5 at% Nd3+-doped KGdP4O12 crystal has a long luminescence lifetime of 300 μs and a high quantum efficiency of 96%.


Introduction
With the development of the diode-pumped solid-state (DPSS) lasers based on neodymium-doped crystals, research on new laser host materials has gained much interest [1]. Nd:YAG crystal is still the most common Nd 3+ laser material because of its good physical and laser properties, but it is limited to low Nd 3+ -doped concentration and narrow absorption bandwidth near the AlGaAs diode emission wavelength of 808 nm. Because the emission wavelength of the laser diode increases at 0.220.3 nm?K 21 of the laser device, the temperature stability of the laser diode needs to be tightly controlled. Therefore, it is necessary to explore new and more efficient crystals for DPSS lasers. Furthermore, the absorption band of laser crystals close to the laser output of AlGaAs (l = 808 nm) needs to have a large full-width at halfmaximum (FWHM).
The title compound, KGdP 4 O 12 , belongs to the broader family of condensed phosphates including those double phosphates of alkali and lanthanide ions with the general formula M I Ln III (PO 3 ) 4 . The structure of M I Ln III (PO 3 ) 4 is highly dependent not only on the sizes of the alkali and the lanthanide ions but also on the crystallization environment. Double tetra-metaphosphates of potassium and gadolinium have low symmetry (monoclinic system) including three structural types with different space groups: P2 1 , P2 1 /n and C2/c. The first two types are written as KGd(PO 3 ) 4 in general and are double polyphosphates. Here, the phosphoric anions have a long-chain geometry of (PO 3 ) n{ n . The C2/c is a cyclotetraphosphate, and these phosphoric anions have a cyclic geometry of [P 4 O 12 ] 42 ; hence its formula is usually written as KGdP 4 O 12 . The growth and structure of KGd(PO 3 ) 4 (P2 1 ) crystal was reported by Parreu et al [20,21] and it is isostructural to KNd(PO 3 ) 4 [22]-a famous Nd stoichiometric laser material. The KGd(PO 3 ) 4 (P2 1 ) can be used as a self-frequency doubling (SFD) laser host material because of its noncentrosymmetrical structure and has attracted more attention [23226]. KGdP 4 O 12 crystallized in space group C2/c and KGd(PO 3 ) 4 in P2 1 /n were successively synthesized in phosphoric acid by Naïli et al [27,28].
Compared to other double phosphates of alkali and lanthanide ions with the general formula M I Ln III (PO 3 ) 4 , KGdP 4 O 12 has the following characteristics as a new laser host material. First, Gd 3+ is much closer to Nd 3+ in size; hence KGdP 4 O 12 is easily substituted with Nd 3+ . Second, raw materials for preparing KGdP 4 O 12 are easily available and very cheap in comparison with double metaphosphates of cesium or rubidium. Third, the spectroscopic spectrum of doped active ions may be broadened because of its low symmetry in the space group C2/c. Finally, the KGdP 4 O 12 crystal has a cyclic geometry that is difficult to cleave in contrast with compounds having long chain geometry.
To the best of our knowledge, the crystal growth of KGdP 4 O 12 has not yet been reported. Here, we have grown a single crystal of Nd 3+ -doped KGdP 4 O 12 with high quality and determined its structural character to be a cyclic geometry via single crystal structural analysis and IR and Raman spectroscopy. The spectral characteristics of the Nd:KGdP 4 O 12 crystal were studied through its absorption spectrum, emission and excitation fluorescent spectra, and fluorescence decay curve. Finally, we measured its thermal stability, specific heat and thermal conductivity.

Crystal growth
The crystal growth used a single-zone vertical tubular electric furnace with a Kantal AF heater. The temperature of the furnace was measured near the heater with a B-type thermocouple and was regulated with a SHIMADEN FP21 type temperature controller.
We There was marked bubbling of NH 3 , CO 2 and H 2 O, and the mixture was heated slowly to a melt. The melt was homogenized at 900uC with a stirrer until the melt was clear and there were no bubbles. Meanwhile, the crucible was rotated alternately clockwise and counterclockwise at 15 revolutions per minute for complete homogenization. The oscillating crucible rotation was also necessary during the crystal growth for effective mass transport. The coldest position in the solution was 0.5 cm below the melt surface and the temperature gradient in the axial direction was 1.3uC?cm 21 . After several attempts, the melt temperature decreased to 812uC. Meanwhile, a thick wire of platinum attached to a corundum rod was introduced into the furnace and was placed above the melt. About 10 min later, the platinum wire was immersed below the melt surface as a seed. While maintaining the oscillating crucible rotation, the melt temperature was reduced at 0.4uC?d 21 until it reached 794uC. During this process, the solute first crystallized on the platinum wire, which served as the seed for larger crystals.
After the growth was completed, the crystals were removed from the solution. The furnace was cooled to room temperature slowly and then some single crystals clustering around the platinum wire were taken out of the furnace. Figure 1(a) shows that the biggest one of these had good transparency with no cracks or inclusions. The colour of the crystal shows that the neodymium ion had already entered the crystal lattice. The stoichiometry of the crystal was analyzed by inductively coupled plasma2atomic emission spectrometer (ICP2AES) (Thermo IRIS Advantage). The density of the Nd:KGdP 4 O 12 crystal at room temperature (23uC) was measured using a density unit of Sartorius balance (YDK012C) by the buoyancy method.

X-ray diffraction analysis
The X-ray diffraction (XRD) technique was used to analyze several aspects of the as-grown Nd:KGdP 4 O 12 single crystal such as determining its structure, identifying the crystalline phase, studying the crystal morphology as well as crystallinity.
A single crystal of Nd:KGdP 4 O 12 cut from the as-grown crystal was used for indexing and intensity data collections on a Rigaku Saturn CCD area detector diffractometer using Mo Ka graphite monochromated radiation (l = 0.71073 Å ). The structure was resolved with the SHELXS-97 computer utility, and refined by a full-matrix least-squares method with SHELXL-97 software.
Polycrystalline powder of Nd:KGdP 4 O 12 was prepared from the as-grown single crystal. The powder X-ray diffraction (PXRD) pattern of Nd:KGdP 4 O 12 was recorded at room temperature with a Fangyuan DX2700 (Dandong, China) powder diffractometer using graphite monochromatized Cu Ka radiation in the 2h range 10u270u with a step of 0.02u.
An XRD goniometer was used to identify and orientate the crystalline forms that comprised the morphological habit for the as-grown single crystal had well-developed facets. XRD was also used to check the quality of the as-grown Nd:KGdP 4 O 12 crystal. The crystal in Figure 1(a) was polished and used as a sample. The rocking curve of the (004) diffraction crystal face was obtained on an XRD diffractometer (PANalytical X'Pert PRO) equipped with a four-crystal Ge (022) monochromator. The setting of the X-ray tube (PW3373/10) was 40 kV and 10 mA. The step time and step size were 0.1 s and 0.001u, respectively.

IR and Raman spectroscopy
The molecular spectroscopy and lattice vibrations of Nd:KGdP 4 O 12 were studied with Fourier transform infrared (FTIR) and Raman spectroscopy. The IR spectrum of a KBr pressed pellet of the powered sample was recorded from 4000 to 400 cm 21 on a Nicolet Magna-IR 560 ESP FTIR spectrometer. The as-grown crystal shown in Figure 1(a) was also studied on a Renishaw inVia Raman microscope. The output wavelength of the excited CW argon laser was 514 nm.

Thermal analysis
The thermal stability of Nd:KGdP 4 O 12 was measured with both thermogravimetry and differential scanning calorimetry (TG2 DSC) using a Netzsch STA 449C from 40 to 1000uC at 10uC?min 21 . The specific heat of the Nd:KGdP 4 O 12 crystal was also measured by DSC using a simultaneous thermal analyzer (Netzsch DSC 200 F3) in an atmosphere of N 2 .
The thermal diffusivity was measured by the laser flash method using a laser flash apparatus (Netzsch LFA 457 Nanoflash) in the 25-300uC range at intervals of 50uC. A 0.95 mm64 mm64 mm sample was used to carry out the measurements whose {001} faces were polished and coated with graphite on both sides.

Spectroscopic characterization
A 4.30-mm thick sample was polished and used to measure the spectroscopic properties. The polished facet was {001}. The unpolarized absorption spectrum was measured on a Hitachi UV2Vis2NIR spectrophotometer (U4100) from 175 to 2000 nm. The fluorescence spectra were recorded using an Edinburgh Instruments FLS920 spectrophotometer with a xenon lamp as the excitation source. The incident light and fluorescence were dispersed with two M300 monochromators using ruling gratings from Bentham Instruments. A cooled Hamamatsu R5509-72 photomultiplier was used for detection. The decay curve was measured using a pulsed xenon lamp as the pump source. The pulse duration was 10 ns, and the excitation wavelength was 808 nm. All the measurements were performed at room temperature (300 K). The splitting of the Nd 3+ emission band was also studied at 10 K using the same FLS920 spectrophotometer equipped with a close cycle helium cryostat (Advanced Research Systems DE202).

Crystal growth and chemical composition
As mentioned above, KGdP 4 O 12 is a polymorphous compound. In addition to the C2/c phase, it has two isomers, i.e. KGd(PO 3 ) 4 in P2 1 and P2 1 /n. Both KGdP 4 O 12 in C2/c and KGd(PO 3 ) 4 in P2 1 /n can be prepared by crystallizing KH 2 PO 4 and Gd 2 O 3 from phosphoric acid below 550uC [27,28]. Because of the lower crystallization temperature and the resulting higher viscosity, it is difficult to grow bulk crystals of Nd:KGdP 4 O 12 . Parreu grew a KGd(PO 3 ) 4 crystal in P2 1 from a melt made from NH 4 H 2 PO 4 , K 2 CO 3 and Gd 2 O 3 at a higher temperature [20]. In studying the crystallization region of KGd(PO 3 ) 4 (P2 1 ) in the ternary system Gd 2 O 3 2K 2 O2P 2 O 5 , six neighbouring phases were found, and KGdP 4 O 12 in C2/c was one of them. By choosing the ratio in the melt composition, we successfully grew bulk crystals of Nd:KGdP 4 O 12 (Figure 1(a)). The composition of a melt influences the formation of phase. A similar situation occurred to a BaTeMo 2 O 9 polymorphous compound recently reported-the a and b crystal phases can be grown from a melt with different compositions [29].
The chemical composition of the as-grown crystal was analyzed by ICP2AES (Table 1). The chemical formula of the crystal corresponds to KGd 0.95 Nd 0.05 P 4 O 12 , which is consistent with results of single crystal diffraction. The distribution coefficient of the neodymium ion at the structural sites of gadolinium ion is defined as tion . This can be calculated from the data in Table 1 and the stoichiometry of the starting materials. The value (0.984) is very close to 1 because the ionic radius of Nd 3+ is only slightly bigger than that of Gd 3+ . This means that the active Nd 3+ ion can be easily doped in the KGdP 4 O 12 host crystal, and that the composition inside a single crystal would be homogeneous in general.   Table S1. The atomic coordinate information is given in Table S2, and the anisotropic displacement parameter is given in  [27]. The Nd 3+ -active ion occupies the lattice site of Gd 3+ . The density of a single crystal of Nd:KGdP 4 O 12 at 23uC is 3.54 g?cm 23 . This agrees well with the calculated value from the crystallographic data.
The basic structural unit of Nd:KGdP 4 O 12 is a centrosymmetric cyclotetraphosphate ring anion of [P 4 O 12 ] 42 that is composed of four PO 4 tetrahedra. Each PO 4 tetrahedron shares its two corners (i.e. bridging oxygen) with the others. The four phosphorus atoms of each [P 4 O 12 ] 42 ring lie in the same plane due to symmetry, and all the oxygen atoms are either above or below the P 4 plane. As shown in projection into the ac plane (Figure 2), the [P 4 O 12 ] 42 rings form layers parallel to (001) at z = 0 and 1/2. A previous description [27] suggesting that these layers were perpendicular to the [001] direction is erroneous because the normal of the (001) faces is different from the [001] direction. We note that the [P 4 O 12 ] 42 rings stack in the structure similar to fish scales. The P 4 plane, i.e. (0.2908, 20.0415, 0.9559), is not parallel to (001), and the interfacial angle between them is 17.08(2)u.
The Gd(Nd)O 8 polyhedra and KO 10 polyhedra link the [P 4 O 12 ] 42 rings to form a three-dimensional framework. The eight vertexes of the Gd(Nd)O 8 polyhedra come from terminal oxygens of the PO 4 polyhedra, and they belong to six [P 4 O 12 ] 42 rings that are equal in the two adjacent layers. The Gd(Nd)O 8 polyhedra do not share any oxygen atoms and are isolated by the PO 4 and KO 10 polyhedra. Ettis et al [27] erroneously reported the shortest distance between the two rare earth ions to be 5.269 Å . Thus, they thought that it was small compared to the other alkali metal and rare earth tetraphosphates. Actually, the shortest Gd(Nd)2Gd(Nd) distance in the Nd:KGdP 4 O 12 crystal is 6.0059 Å , which is actually relatively long for tetraphosphates. Figure S1 shows the experimental XRD pattern of the pulverized Nd:KGdP 4 O 12 crystal as well as our simulated pattern. The peak positions and diffraction intensities are consistent between the experimental and simulated XRD patterns. This confirms that our proposed structure is accurate.
For a laser crystal, the mechanical strength is due to the host material. Cleavage behaviour of the host crystal is related not only to the bond strength but also the geometry. There are three kinds of chemical bonds in potassium gadolinium phosphates: P2O, Gd2O and K2O. According to their covalent character, P2O bonds are much stronger than Gd2O and K2O bonds. The K2 O bonds are the weakest of all. This is confirmed by the calculated results of bond valences shown in

Crystal morphology and rocking curve
The as-grown crystal has some well-developed facets (Figure 1(a)) oriented using an XRD goniometer. The habit contains crystalline forms {001}, {010}, {110} and {021}. The measured interfacial angles among the facets were in good agreement with the calculated angles.
The morphology of a crystal reflects its structure. Taking into account only the point group and the crystal cell parameters, we established the theoretical morphological scheme of the Nd:KGdP 4 O 12 crystal via the WinXMorph software [31] using the Bravais-Friedel and Donnay-Harker (BFDH) law. This morphological scheme is shown in Figure 1(b). Table 2  While the as-grown crystal did not show a completely ideal morphology due to the restriction of the crystal growth method, all the predicted crystalline forms can be found on the as-grown single crystal. It should be noted that the f 1 111g facet was smaller and limited by the growth conditions. Thus, it was quite pristine and we did not identify f 1 111g for the as-grown crystal. The narrow tetragonal facet (021) on the as-grown crystal shows that {021} has a faster growth rate, which is consistent with the predicted scheme. The ideal morphology of KGdP 4 O 12 is like a drum and hexagonal {001} facets reveal much more than other facets similar to drumheads. This relates to the structural feature of KGdP 4 O 12 in which [P 4 O 12 ] 42 rings regularly array to form layers parallel to (001). The distance between the two [P 4 O 12 ] 42 rings adjacent to the two layers is larger than that from the same layer.
The rocking curve of the (004) diffraction plane of the as-grown Nd:KGdP 4 O 12 crystal is presented in Figure S2. The diffraction peak is intense and with good symmetry without splitting, and its FWHM value is 0.006u. These results demonstrate that the asgrown Nd:KGdP 4 O 12 crystal has a nearly perfect lattice structure.

IR and Raman spectra
The IR and Raman spectra of Nd:KGdP 4 O 12 at room temperature are shown in Figure 3. The vibration frequencies above 600 cm 21 and their corresponding assignment are listed in Table S5 and are based on published findings [27,32]. In the lowfrequency region below 600 cm 21 , it is very difficult to distinguish the symmetric and antisymmetric bending modes of the O2P2O and P2O2P groups.
Compared to other cyclotetraphosphate compounds, we notice an absence of bands in the 75021000 cm 21 region of the IR spectrum. This result is in agreement with the cyclic structure of Nd:KGdP 4 O 12 . The non-coincidence of the observed bands in the IR and Raman spectra indicates that the [P 4 O 12 ] 42 anion is centrosymmetric, which is in agreement with the above structural results.
A distinguishing characteristic also exits in the Raman spectrum. The symmetric stretching vibration of the P2O2P linkage-n s (P2O2P)-has a single peak at 684 cm 21 . This is the strongest of all the Raman vibration peaks. However, for n s (P2 O2P) of KGd(PO 3 ) 4 (P2 1 ), two peaks are observed in the range of 660-730 cm 21 and have the same intensity [24]. These are less than one half of the intensity of the n s (O2P2O). That is because of the symmetrical difference in the crystal structures and the different positions of the lanthanide and alkali ions. The obvious differences between the IR and Raman spectra of KGdP 4 O 12 and KGd(PO 3 ) 4 highlight how IR and Raman spectroscopy can identify the structure of alkali-metal lanthanide metaphosphates.   Figure S3. A single sharp endothermic peak is observed at 920uC, which exhibits the characteristics of a first-order phase transition. The sample weight of Nd:KGdP 4 O 12 does not show representative variation in the measurement temperature up to 1000uC. Apart from a crystalline phase, an amorphous phase was formed in the sample cell after the DSC measurement. The M I Ln III P 4 O 12 compounds often decompose irreversibly into lanthanide trimetaphosphates and alkali metal metaphosphates [33]. We thus conclude that the exothermic peak at 920uC may be related to a decomposition of Nd:KGdP 4 O 12 in accordance with the reaction: The decomposition temperature of the Nd:KGdP 4 O 12 crystal is much higher than that of KGd(PO 3 ) 4 (P2 1 ) (878uC) [20]. The difference in stability between these two potassium gadolinium metaphosphates can be explained by comparing their crystal data. . Thus, the thermal stability of the cyclic C2/c space group should be better than that of the chain structure in the P2 1 space group. For laser crystal materials, the damage threshold and possible laser applications can be greatly influenced by the specific heat. The specific heat is also an important value used to calculate thermal conductivity. Figure S4 shows that the constant pressure specific heat of the Nd:KGdP 4 O 12 crystal varies as a function of the temperature. The specific heat of the Nd:KGdP 4 O 12 crystal at room temperature (25uC) is 0.521 J?g 21 ?K 21 . This is comparable to Nd:YAG (0.59 J?g 21 ?K 21 ) and Nd:YVO 4 (0.51 J?g 21 ?K 21 ) [34]. Because Nd:YAG and Nd:YVO 4 both have a high optical damage threshold, it follows that Nd:KGdP 4 O 12 should also have a high damage threshold [5]. The specific heat increases almost linearly from 0.485 to 0.799 J?g 21 ?K 21 with temperature increases from 210 to 510uC. This suggests that the Nd:KGdP 4 O 12 crystal can tolerate even more thermal energy at a high temperature. Figure 4 shows the thermal diffusivity and thermal conductivity of the Nd:KGdP 4 O 12 crystal. The thermal conductivity (k) was calculated using the measured thermal diffusivity and specific heat according to k = lrC p , where l, r and C p denote the thermal diffusivity, density and specific heat of the crystal at the same temperature, respectively. The thermal diffusivity of the crystal is 0.906 mm 2 ?s 21 along the direction perpendicular to (001), i.e. the c* direction, at 27.8uC. The calculated thermal conductivity of the crystal is 1.66 W?m 21 ?K 21 along the c* direction correspondingly. This value is similar to that of Nd:SrLaGa 3 O 7 [4] and larger than those of some laser crystals with broad absorption such as NaLn(WO 4 ) 2 [2,3] and Ca 3 La 2 (BO 3 ) 4 [5]. In laser designs, the thermal loading causes a temperature gradient in the crystal and leads to thermal expansion that results in thermal lensing and other thermo-optic effects. All these effects would cause a decline in the quality of the laser beams and even crack the active crystal. The deposited heat would be easily transferred to the environment if the laser crystal possesses high thermal conductivity-this would minimize the thermal loading effects. It is clearly shown in Figure 4 that the thermal diffusivity and thermal conductivity component of the Nd:KGdP 4 O 12 crystal decreases with increasing temperature. Thermal conductivity is dominated primarily by its phonon thermal conductivity for the dielectric. This has two important determinants: the heat capacity and the phonon mean free path. Though the heat capacity of the Nd:KGdP 4 O 12 crystal increases with temperature, the phononphonon scattering becomes much stronger and the phonon mean free path markedly decreases. Therefore, the final result is that the thermal conductivity decreases with temperature increases. Figure 5 shows the unpolarized absorption spectrum of the Nd:KGdP 4 O 12 crystal at room temperature. The UV cut off is below 200 nm. This means that a large band gap would increase the damage resistivity of the crystal. All absorption lines are due to the 4f 3 -4f 3 transition of the Nd 3+ ions. Very strong absorption lines occur near 581, 747 and 798 nm, corresponding to the transitions of 4 I 9/2 R 2 G 7/2 + 4 G 5/2 , 4 I 9/2 R 4 S 3/2 + 4 F 7/2 and 4 I 9/ 2 R 2 H 9/2 + 4 F 5/2 , respectively. The interesting feature of this spectrum relative to the potential application in the diode-pumped laser, is the strong absorption band with a FWHM of 14.8 nm near 798 nm. This is close to the laser output of the AlGaAs laser diode. It is well known that temperature control of the laser diode is crucial, because the emission wavelength increases with temperature. Large line-width in the Nd:KGdP 4 O 12 crystal is very suitable for laser diode pumping because it is not temperature dependent.

Absorption spectrum
The Judd-Ofelt theory was used to calculate the optical parameters of the electric dipole transition within the 4f electronic configuration of the Nd 3+ ion. Important parameters included the oscillator intensity parameters (V l ), spontaneous emission probability (A), fluorescence branching ratio (b) and radiative lifetime (t r ) [35,36]. A detailed calculation procedure was used similar to literature [37]. Nine absorption bands were used to obtain the intensity parameters (V 2,4,6 ), which are V 2 = 1.79610 -20 cm 2 , V 4 = 3.48610 -20 cm 2 and V 6 = 6.05610 -20 cm 2 , respectively. The values of the experimental (f exp ) and calculated (f cal ) oscillator strengths are listed with the absorption cross sections (s) in Table 3. The root mean square error (d) between f exp and f cal is 6.97610 28 . The properties of 4 F 3/2 R 4 I J' (J' = 9/2, 11/2, 13/2, 15/2) transitions are well known. These are the main channels of a Nd 3+ laser and depend only on the values of V 4 and V 6 because the reduced matrix elements |,||U (2) ||.| for these transitions are zero. The spectroscopic quality parameter (X = V 4 /V 6 ) is 0.57, which is less than 1 and indicates that the emission to the 4 I 11/2 manifold is more feasible than that to 4 I 9/2 manifold. This result can be confirmed by the fluorescence emission spectrum (see below). In this regard, it is similar to the most common laser materials based on the Nd 3+ ion including Nd:YVO 4 and Nd:YAG. Table 4 presents the spontaneous emission probability (A), radiative branching ratio (b) and radiative lifetime (t r ) of the multiplet 4 F 3/2 . It can be seen that the calculated branching ratios are in good agreement with the experimental ones obtained from the measured fluorescence spectrum later.

Fluorescence spectra analysis
The room temperature fluorescence spectrum of the Nd:KGdP 4 O 12 crystal at 808 nm excitation was recorded from 830 to 1600 nm ( Figure 6). Three emission bands corresponding to the 4 F 3/2 R 4 I 9/2 , 4 I 11/2 and 4 I 13/2 transitions are observed at 860-920 nm, 1040-1080 nm and 1300-1380 nm, respectively. The room temperature emission FWHM of the transition 4 F 3/ 2 R 4 I 11/2 is 14 nm. This is evidence of the dominant contribution of inhomogeneous broadening to the spectra line-width in the low symmetry structure. The emission cross sections for the laser channel from 4 F 3/2 to 4 I 9/2 , 4 I 11/2 and 4 I 13/2 were calculated through the Füchtbauer-Ladenburg (FL) formula [12,18]. These values are 1.11610 220 , 6.25610 220 and 2.13610 220 cm 2 , respectively.
The emission spectra of the Nd:KGdP 4 O 12 at low temperature were studied to determine the Stark sublevels of the 4 I 13/2 and 4 I 11/2 excited multiplets and the ground state multiplet 4 I 9/2 . Figure S5 shows the emission spectra of Nd:KGdP 4 O 12 at 10K and 300K. These correspond to the transitions from 4 F 3/2 to 4 I 9/2 , 4 I 11/2 and 4 I 13/2 levels, respectively. Figure 7 shows a schematic diagram of the Stark sublevels of Nd 3+ in KGdP 4 O 12 crystals obtained from the low temperature emission spectra. The 4 F 3/2 term is split into two components with D = 70 cm 21 . The low value of the crystal field splitting indicates that the crystal field of the KGdP 4 O 12 host is weak. This is due to weak interactions between rare earth ions and oxygen atoms in contrast with the strong covalent bonds between P and O. In this regard, the Nd:KGdP 4 O 12 crystal is similar to other double phosphates of alkali and lanthanide ions [38,39].
The fluorescence lifetime (t f ) of the 4 F 3/2 energy level of Nd 3+ in the Nd:KGdP 4 O 12 crystal was determined to be 300 ms by fitting the decay curve exponentially ( Figure S6). The long fluorescence lifetime would be beneficial to high energy storage during laser operation. Though the interatomic distance (6.0059 Å ) between the nearest rare earth ions for KGdP 4 O 12 is smaller than that for KGd(PO 3 ) 4 (P2 1 ) (6.5865 Å ), it does not obviously affect the fluorescence lifetime. This is a function of the isolated GdO 8 polyhedra in KGdP 4 O 12 . The t f of the Nd:KGd(PO 3 ) 4 crystal (P2 1 ) is 246 ms with a very low doping concentration (0.6 at%) [38]. The values of Nd:LaP 5 O 14 and Nd:GdP 5 O 14 crystals are around 310 ms when doped with ,1 at% Nd 3+ [40].
It is well known that the t f of Nd 3+ ions generally decreases with increasing Nd 3+ concentrations due to the interactions between the Nd 3+ ions, i.e. the concentration quenching effect. However, The calculated t r of the 4 F 3/2 energy level is 312 ms (Table 4). Thus, the fluorescent quantum efficiency (g = t f/ t r ) of the 4 F 3/2 level is 96%. For solid-state laser materials, low crystal field strength leads to weak electron-phonon interactions and further leads to high quantum efficiency. The previous data describing weak Stark splitting of the Nd 3+ energy levels further supports the calculated result showing that Nd:KGdP 4 O 12 has a high quantum efficiency.
Some spectroscopic properties of the Nd:KGdP 4 O 12 crystal are listed in Table 5 as well as those of other Nd 3+ -doped crystals including YAG and YVO 4 crystals with ordered structures and SrLaGa 3 O 7 and KBaGd(MoO 4 ) 3 crystals with disordered structures. The FWHM value of Nd:KGdP 4 O 12 at l a is not only much larger than those of Nd:YAG and Nd:YVO 4 crystals, but also larger than those of Nd:SrLaGa 3 O 7 and Nd:KBaGd(MoO 4 ) 3 crystals. This suggests that the Nd:KGdP 4 O 12 crystal can be pumped more effectively by AlGaAs laser diode, but not be restricted to the temperature stability of the output wavelength in the laser diode. The quality factor (M) of a laser material is proportional to its doping concentration (N) and fluorescence lifetime (t f ). A high value of M generally means a low oscillation threshold in subsequent laser operations. Therefore, the high doping concentration and long fluorescence lifetime of Nd:KGdP 4 O 12 makes it possible to achieve a continuous wave laser action. In addition, Nd:KGdP 4 O 12 can be used as a tunable laser material with a very broad emission band near 1060 nm.

Conclusions
We successfully grew a macro-defect-free single crystal of 5 at% Nd 3+ -doped KGdP 4 O 12 with very good crystallinity using TSSG2 SC from self-flux. To the best of our knowledge, both the growth of the bulk crystals and doping with neodymium are shown here for the first time in KGdP 4 O 12 . The Nd:KGdP 4 O 12 crystallizes in space group C2/c, and the phosphoric anions have a cyclic geometry of [P 4 O 12 ] 42 . The crystal has good chemical stability and does not show cleavage like its isomer KGd(PO 3 ) 4 . The distribution coefficient of neodymium ion is very close to 1, so gadolinium ion can be evenly substituted. The absorption and fluorescence spectra of the Nd:KGdP 4 O 12 crystal were investigated at room temperature. The peak absorption cross section at 798 nm is 2.59610 220 cm 2 with a FWHM of 14.8 nm. Such a broad FWHM in the absorption band is suitable for InGaAs laser diode pumping. The radiation and fluorescence lifetimes of the excited state 4 F 3/2 are 312 and 300 ms, respectively. These result in a high luminescent quantum efficiency of 96%. The emission from the 4 F 3/2 energy level to the 4 I 11/2 manifold is more feasible than that of the 4 I 9/2 manifold. The peak emission cross section at 1057 nm (corresponding to 4 F 3/2 R 4 I 11/2 ) is 6.25610 220 cm 2 . In addition, the KGdP 4 O 12 crystal presents positive thermal characteristics. It has a good thermal stability with decomposition at 920uC. At room temperature, the specific heat and the thermal conductivity along the c* direction are 0.521 J?g 21 ?K 21 and 1.66 W?m 21 ?K 21 , respectively. In summary, the Nd:KGdP 4 O 12 crystal may be regarded as a potential solid-state laser material for laser diode pumping.