The Anisotropic Thermal Expansion of Non-linear Optical Crystal BaAlBO3F2 Below Room Temperature

Introduction

Thermal expansion behavior is a crucial performance factor for materials used in lasers, due to its affinity with the ability of the optical outputting (Wynne et al., 1999; Wang et al., 2007; Mangin et al., 2011; Ito et al., 2017; Fang et al., 2018). Non-linear optical (NLO) crystal, a type of crystals for laser frequency conversion, has played a key role in the broadening of laser spectrum and become one of the most prevailing branch of optical materials (Chen et al., 1999, 2012; Cyranoski, 2009; Meng et al., 2009; Lin et al., 2014). Among the commercial NLO crystals, BaAlBO₃F₂(BABF) crystal, is a very excellent member, and exhibit many superiority over its congeners in the output of laser in 532 and 355 nm (Zhou et al., 2009; Hu et al., 2011; Yue et al., 2011; Yang et al., 2016). As for the practical application of BABF, the optical spot with high-power is usually focused on one point of the crystal, and thus “tearing” stress resulted from anisotropic thermal expansion would severely affect the beam quality and outputting power. The previous study about the thermal property of BABF revealed that above room temperature, its thermal expansion coefficient along c-axis is about eight times of that along a(b)-axis, and the strongly anisotropic thermal expansion has been the major disadvantage to restrict its practical application (Yue et al., 2011). However, the microscopic mechanism of the thermal expansion anisotropy in BABF has not been investigated. Recently, it is proposed that the thermal effect on the performance of some optical crystal used in laser can be eliminated to some extent below room temperature (Marrazzo et al., 2016; Veselov et al., 2016; Roitero et al., 2018), which provide a feasible method to overcome the problem about the thermal effect. However, the thermal expansion behavior of BABF under low temperature has not been investigated yet, and its performance index below room temperature remain unclear. Therefore, it is desirable to perform a study about the thermal expansion to give a comprehensive evaluation on the application prospect below room temperature and elaborate the structure-property relationship of the anisotropic thermal expansion of BABF.

In this work, using variable-temperature X-ray diffraction, the thermal expansion behavior of BABF below room temperature is studied, and the mechanism of the thermal expansion behavior is elaborated by Raman spectrum and first-principles phonon analysis. The optical property below room temperature is also investigated by first-principles calculation. It is elucidated that BABF still exhibit a relatively strong anisotropic thermal expansion below room temperature, but its anisotropy is less prominent than that above room temperature. The optical performance is also slightly improved under low temperature. These result indicates that the performance below room temperature for the laser generation is superior than that above room temperature, and BABF is a potential NLO crystals in cryogenic system.

Experimental and Computational Method

Sample Preparation

Polycrystalline BABF was synthesized through solid-state reaction. Analytically pure BaF₂, Al₂O₃, and B₂O₃ in stoichiometric ratio as the starting materials were mixed homogeneously by agate mortar. The well-grinded reactant then was placed into muffle furnace and was heated gradually up to 800°C at the rate of 0.5°C/min with several intermediate careful grinding at 300, 600, and 800K. After cooling to room temperature, the white powder of target compound was obtained.

Variable Temperature X-Ray Diffraction

Variable temperature X-ray diffraction patterns were recorded from 13 to 300K with the internal of 20K. Each pattern was recorded with Bruker D8 advanced X-ray diffractometer Cu Kα radiation (Kα1 = 1.5406Å and Kα2 = 1.5443Å) on the finely grounded powder samples. The angular scanning range were set to 10–90° with a step of 0.01° and scanning rate 0.5 s/step. The crystal structures under different temperature were refined by Rietveld method (Rietveld, 1969) using TOPAS 4.2 program (BAXS, 2008). Based on the refined cell parameters. The thermal expansion coefficient was fitted by the PASCal software (Cliffe and Goodwin, 2012).

Raman Spectrum

The Raman pattern was recorded from 100 to 1,500 cm⁻¹ at room temperature, using in Via-Reflex, equipped with a solid state laser with a wavelength of 532 nm. In order to improve the signal to noise ratio of the spectra, 10 integrations were carried out with an integration time of 10 s at a nominal resolution of 1 cm⁻¹ and a precision of 1 cm⁻¹.

First-Principles Calculation

The first-principles calculation were performed by CASTEP (Clark et al., 2005), a plane-wave pseudopotential total energy package based on density functional theory (DFT) (Kohn and Sham, 1965; Payne et al., 1992). The functionals developed by Perdew, Burke, and Ernzerhof (PBE) (Perdew et al., 1996) in generalized gradient approximation (GGA) (Perdew and Wang, 1992) form were adopted to describe the exchange-correlation energy. The optimized norm-conserving pseudopential (Hamann et al., 1979) in (Kleinman and Bylander, 1982) form were used to model the effective interaction between the valcence electrons and atom cores, which allow us to use a small plane basis set without compromising the accuracy required by the calculation. High kinetic energy cutoff of 900 eV and Monkhorst-pack (Monkhorst and Pack, 1976) k-point mesh spanning less than 0.04 Å⁻¹ in the Brillouin zone were chosen. The vibrational property was calculated by linear response formalism (Baroni et al., 2001), in which the phonon frequencies were obtained by the second derivative of the total energy with respect to a given perturbation. The band gaps at variable temperature were predicted by hybridized functionals PBE0 based on the refined structure at respective temperature, and the refractive index and SHG coefficients were calculated by Kramers–Kronig transform based on the electronic transition matrix (PED, 1985) and the software developed by our group based on length-gauge formalism (Lin et al., 1999, 2014).

Result and Discussion

The crystal structure of BABF was determined by Hu et al. (2011). BABF possesses a hexagonal lattice with layered-tendency structure, as displayed in Figure 1. One boron atom is bonded with three oxygen atoms to form a planar BO₃ triangles, and one aluminum atom coordinates with three oxygen and two fluorine atoms to generate the AlO₃F₂ trigonal bipyramid. The BO₃ triangles and AlO₃F₂ trigonal bipyramids are aligned alternatively in the ratio 1:1 by sharing the vertical oxygen atoms, giving rise to the infinite two-dimensional [AlBO₃F₂]_∞ layer, in which the BO₃ triangle and the AlO₃ base of AlO₃F₂ trigonal bipyramids are aligned absolutely parallel to (a, b) plane. [AlBO₃F₂]_∞ layer are further connected with each other via the Coulomb interaction between the dangling fluorine atoms and interstitial barium atoms to generate the layered-tendency structure.

FIGURE 1

Figure 1. Crystal structure of BaAlBO₃F₂ viewed along (A) (a,b) plane and (B) c-axis. The BO₃ and AlO₃F₂ groups are represented by pink triangles and light blue trigonal bipyramids, respectively.

The variable-temperature X-ray diffraction (VT-XRD) patterns are plotted in Figure 2A. It is observed that no new peaks appear as temperature elevate from 13 to 300K, indicating the absence of structural phase-transition and manifesting its high thermodynamical stability. The temperature-deduced diffraction peak-shifting of the lattice plane corresponding to c-axis is much more prominent than that of a(b)-axis (Figure 2B), suggesting the high anisotropic thermal expansion in BABF. According to the cell parameters extracted from the VT-XRD by Rietveld refinement (Figure S1, Figure 2C and Table S1), the average thermal expansion coefficient are 4.42(39)/MK and 20.08(71)/MK along a(b)-axis and c-axis respectively, which confirmed that the thermal expansion along c-axis is much more prominent than that along a(b)-axis. More importantly, the ratio (4.5:1) between the thermal expansion coefficients along c- and a(b)-axis below room temperature is much lower than that above room temperature (7.6:1), and the low anisotropy of the thermal expansion below room temperature is more favorable for elimination of the stress induced by temperature gradient when irradiated by laser. Additionally, It should be emphasized that at low temperature (13–100K) the a(b)-axis exhibit very low thermal expansion with the fitted linear thermal expansion coefficient 0.11(25)/MK, which is lower than that of c-axis [2.44(98)/MK] by one order and can be categorized to typical zero thermal expansion.

FIGURE 2

Figure 2. Result of VT-XRD on BABF: (A) the VT-XRD patterns, (B) 2θ-temperature contour map of the diffraction peak of (002) and (010) plane and (C) evolution of cell parameters with respect to temperature. The insert in (C) displays the special distribution of thermal expansion coefficient plotted by PASCal software (Cliffe and Goodwin, 2012).

The most intuitive way to investigate the mechanism of the thermal expansion is to trace the modification of bond length and angles with respect to temperature. The a(b)-axis is exclusively determined by B-O, Al-O bonds, and ∠ O-O-O angle, while c-axis is dominantly determined by the Ba-F and Al-F bond. According to the temperature-dependent bond length and angles (Figure 3 and Table S2), despite the clutter of the point due to the difficulty to exactly determine the position of the constituted light atoms in BABF, a roughly changing tendency can be observed. Below 140K the bond length and angles accounting for the change of a(b)-axis almost keep constant, resulting in the rigidity of a(b)-axis under low temperature. And dramatic modification only occur above 140K: B-O bond, and ∠ O-O-O angle increase from 1.301 to 1.340Å (by 3%) and 60.892–62.437° (by 2.5%) respectively as temperature increase from 140K to 300K, which both positively contribute to the thermotropic expansion. Extruded by the increase of ∠ O-O-O angle, Al-O bond decrease from 1.919 to 1.865Å (by 2.8%), which slightly cancel out the expansion effect originated from the increase of B-O bond and ∠ O-O-O angle, leading to the relative low thermal expansion along a(b)-axis. On the contrary, both Al-F and Ba-F bond increase as temperature increase, and thus c-axis exhibit a normal thermal expansion ~10/MK of inorganic crystals.

FIGURE 3

Figure 3. Function of refined bond length and angles with respect to temperature in BABF.

The thermal expansion of materials is closely related with the temperature-induced stimulation of lattice vibration, and it is anticipated that first-principles phonon mode assignment would compensate the disadvantage of the littery temperature-dependent bond length and angles and make the mechanism of the thermal expansion in BABF more clear. The irreducible representation of P-62c space group at Γ-point yields a sum of 30E+7A1+8A2 phonon modes (see Table S3). According to the calculated Raman spectrum (Figure 4), six principal peaks are observed at around 87, 266, 338, 473, 1,030, and 1,331cm⁻¹ respectively. The calculated Raman spectrum is in rather agreement with the measured ones, verifying the accuracy of the computational method. According to the vibrational vector projected onto the real space of respective phonon modes, it is revealed that as frequency increase the phonon roughly experience a evolution from the type-I (mode 1–7), -II (mode 8–29) to -III (mode 30–45). In type-I, the configuration of 2D [AlBO₃F₂]_∞ layer keep constant and it vibrate as a rigid unit along c-axis with the reverse phase with that of interstitial barium atoms. This indicates that at low temperature the vibration modes related with the stretch between [AlBO₃F₂]_∞ and barium would be firstly stimulated, in which the size of [AlBO₃F₂]_∞ layer along a(b)-axis remain constant. Therefore, below 100K, BABF exhibit slight thermal expansion along c-axis and rigidity along a(b)-axis. In type-II, the oxygen atoms within [AlBO₃F₂]_∞ layer vibrate along the direction almost vertical to B-O bonds. This implies that in these modes, the BO₃ triangles rotate within (a,b) plane as the rigid units. This effect corresponds to the slight contraction of Al-O bond in AlO₃F₂ trigonal bipyramids and also would lead to the expansion of the interspace within [AlBO₃F₂]_∞ layer, which eventually lead to the weak thermal expansion along a(b)-axis. Meanwhile, the stretch vibration of Al-F bonds also partially account for the modes of type-II, which enhanced the thermal expansion along c-axis. In type-III modes, the phonon vibration mainly originate from the stretch of B-O bonds, which is related with the elongation of B-O bonds. It should be emphasized that high temperature is required to stimulate these vibration modes (such as ~ 1300K for the mode at 1030 cm⁻¹), and these modes can only afford the thermal expansion at high temperature and almost contribute nothing to that below 300K. Moreover, as temperature increase, all the Raman peaks are red-shifted, affording the normal general thermal expansion. Besides, it should be emphasized that the weight of modes related with the thermal expansion along a(b)-axis at 300K is much more prominent than that at 13K, which also result in the enhanced thermal expansion along a-(b) axis as temperature increase. Therefore, the anisotropic thermal expansion behavior in BABF is mainly attributed to the progressive stimulation of the respective vibration modes related with the expansion of c-axis to that of a(b)-axis. Moreover, it should be emphasized that as temperature increase, the amplitude of atomic vibration would increase, and the phonon anharmonicity would be enhanced. Since the thermal expansion is originated from phonon anharmonicity, the thermal expansion of both a-axis [4.42(39)/MK] and c-axis [20.08(71)/MK] between 13 and 300 are less prominent than that those between 303 and 1,073K [6.3 and 48.1/MK, respectively (Yue et al., 2011)].

FIGURE 4

Figure 4. Measured and calculated Raman spectrum of BABF. The modes with the highest peaks at 87, 473, and 1,030 cm⁻¹ are used to schematically describe the atomic vibration of type-I, II, and III phonon modes respectively. The peaks at around 720cm⁻¹ labeled by green triangle in the experimental spectrum is attributed to the impurity.

Additionally, considering that the optical performance is also crucial to the usage in apparatus operated at low temperature, the temperature-dependent optical property was also studied based on the refined structure by first-principles calculation (Figure 5 and Table S4). Accordingly, both the SHG coefficient and birefringence present ultra-stability under temperature fluctuation, and the band gaps are broadened as temperature decrease, which guarantee the optical transmittance of BABF below room temperature. All these observation elucidated that the good optical performance of BABF can be well-kept, which is favorable to its practical application in the apparatus operated below room temperature.

FIGURE 5

Figure 5. Function of the calculated band gap, SHG coefficient and birefringence (@1,064 nm) of BABF with respect to temperature.

Conclusion

The thermal expansion behavior below room temperature of BABF, an important ultraviolet NLO crystal, is investigated. It is revealed that the BABF exhibit relatively high anisotropic thermal expansion, and the near-zero thermal expansion behavior along a(b)-axis below 100K was observed. Based on the refined temperature-dependent crystal structure, Raman spectrum and first-principles vibration phonon analysis, it is elucidated that the anisotropic thermal expansion behavior in BABF mainly stem from the progressive stimulation of the phonon modes related with the thermal expansion along a(b)- and c-axis. Moreover, it is revealed that the optical performance of BABF can also be well-kept under low temperature, which is very favorable for its practical application. Our study exhibit a comprehensive investigation about the thermal and optical property below room temperature and shed light on the structural origin of the thermal expansion of BABF. We believed that the conclusion deduced from our study further deepens the understanding about the performance of BABF and would be beneficial for its practical application under complex environment.

Author Contributions

XJ, NW, and MM performed the experiment. XJ performed the first-principles simulation. WW, SG, RH, and LL help in the variable X-ray diffraction. ZL, ZH, and XJ written the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by the National Scientific Foundations of China (Grants 11474292, 51702330, 11611530680, 91622118, and 91622124), Russian Foundation for Basic Research (Grant 17-52-53031), the Special Foundation of the Director of Technical Institute of Physics and Chemistry (TIPC) and the Youth Innovation Promotion Association, CAS (outstanding member for ZL and Grant 2017035 for XJ).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2018.00252/full#supplementary-material

References

Baroni, S., de Gironcoli, S., Dal Corso, A., and Giannozzi, P. (2001). Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 73, 515–562. doi: 10.1103/RevModPhys73.515

CrossRef Full Text | Google Scholar

BAXS (2008). Bruker AXS TOPAS V4: General Profile and Structure Analysis Software for Powder Diffraction Data. Karlsruhe: User's Manual.

Chen, C. T., Sasaki, T., Li, R. K., Wu, Y. C., Lin, Z. S., Yoshimura, M., et al. (2012). Nonlinear Optical Borate Crystals-Principles and Applications. Germany: Wiley-VCH.

Chen, C. T., Ye, N., Lin, J., Jiang, J., Zeng, W. R., and Wu, B. C. (1999). Computer-assisted search for nonlinear optical crystals. Adv. Mater. 11, 1071–1078. doi: 10.1002/(SICI)1521-4095(199909)11:13<1071::AID-ADMA1071>3.0.CO;2-G

CrossRef Full Text | Google Scholar

Clark, S. J., Segall, M. D., Pickard, C. J., Hasnip, P. J., Probert, M. J., and Payne, M. C. (2005). First principles methods using CASTEP. Z. Kristallogr. 220, 567–570. doi: 10.1524/zkri.220.5.567.65075

CrossRef Full Text | Google Scholar

Cliffe, M. J., and Goodwin, A. L. (2012). PASCal: a Principal axis strain calculator for thermal expansion and compressibility determination. J. Appl. Crystallogr. 45, 1321–1329. doi: 10.1107/S0021889812043026

CrossRef Full Text | Google Scholar

Cyranoski, D. (2009). Materials science: china's crystal cache. Nature 457, 953–955. doi: 10.1038/457953a

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, Z., Liu, L., Wang, X., and Chen, C. (2018). Thermo-physical properties of a new UV nonlinear optical crystal: NaSr3Be3B3O9FF4. Appl. J. Crystallogr. 51, 357–360. doi: 10.1107/S1600576718001218

CrossRef Full Text | Google Scholar

Hamann, D. R., Schluter, M., and Chiang, C. (1979). Norm-conserving pseudopotentials. Phys. Rev. Lett. 43, 1494–1497. doi: 10.1103/PhysRevLett.43.1494

CrossRef Full Text | Google Scholar

Hu, Z., Yue, Y., Chen, X., Yao, J., Wang, J., and Lin, Z. (2011). Growth and structure redetermination of a nonlinear BaAlBO3F2 crystal. Solid State Sci. 13, 875–878. doi: 10.1016/j.solidstatesciences.2011.03.002

CrossRef Full Text | Google Scholar

Ito, I., Silva, A., Nakamura, T., and Kobayashi, Y. (2017). Laser based on low thermal expansion ceramic cavity with 4.9 MHz/s frequency drift. Opt. Express 25, 26020–26028. doi: 10.1364/OE.25.026020

CrossRef Full Text | Google Scholar

Kleinman, L., and Bylander, D. M. (1982). Efficacious form for model pseudopotentials. Phys. Rev. Lett. 48, 1425–1428. doi: 10.1103/PhysRevLett.48.1425

CrossRef Full Text | Google Scholar

Kohn, W., and Sham, L. J. (1965). Self-consistent equations including exchange and correlation effects. Phys. Rev. 140:1133. doi: 10.1103/PhysRev.140.A1133

CrossRef Full Text | Google Scholar

Lin, J., Lee, M. H., Liu, Z. P., Chen, C. T., and Pickard, C. J. (1999). Mechanism for linear and nonlinear optical effects in beta-BaB2O4 crystals. Phys. Rev. B 60, 13380–13389. doi: 10.1103/PhysRevB.60.13380

CrossRef Full Text | Google Scholar

Lin, Z., Jiang, X., Kang, L., Gong, P., Luo, S., and Lee, M.-H. (2014). First-principles materials applications and design of nonlinear optical crystals. J. Phys. D 47:253001. doi: 10.1088/0022-3727/47/25/253001

CrossRef Full Text | Google Scholar

Mangin, J., Mennerat, G., and Villeval, P. (2011). Thermal expansion, normalized thermo-optic coefficients, and condition for second harmonic generation of a Nd:YAG laser with wide temperature bandwidth in RbTiOPO4. J. Opt. Soc. Am. B Opt. Phys. 28, 873–881. doi: 10.1364/JOSAB.28.000873

CrossRef Full Text | Google Scholar

Marrazzo, S., Goncalves-Novo, T., Millet, F., and Chanteloup, J.-C. (2016). Low temperature diode pumped active mirror Yb³⁺: YAG disk laser amplifier studies. Opt. Express 24, 12651–12660. doi: 10.1364/OE.24.012651

CrossRef Full Text | Google Scholar

Meng, J., Liu, G., Zhang, W., Zhao, L., Liu, H., Jia, X., et al. (2009). Coexistence of fermi arcs and fermi pockets in a high-Tc copper oxide superconductor Nature 462, 335–338. doi: 10.1038/nature08521

CrossRef Full Text | Google Scholar

Monkhorst, H. J., and Pack, J. D. (1976). Special points for brillouin zone integrations. Phys. Rev. B 13, 5188–5192. doi: 10.1103/PhysRevB.13.5188

CrossRef Full Text | Google Scholar

Payne, M. C., Teter, M. P., Allan, D. C., Arias, T. A., and Joannopoulos, J. D. (1992). Iterative minimization techniques for abinitio total energy calculations molecular dynamics and conjugate gradients. Rev. Mod. Phys. 64, 1045–1097. doi: 10.1103/RevModPhys.64.1045

CrossRef Full Text | Google Scholar

PED (1985). Handbook of Optical Constants of Solids. New York, NY: Academic Press.

Perdew, J. P., and Wang, Y. (1992). Pair-distribution function and its coupling-constant average for the spin-polarized electron-gas. Phys. Rev. B 46, 12947–12954. doi: 10.1103/PhysRevB.46.12947

PubMed Abstract | CrossRef Full Text | Google Scholar

Perdew, J. P., Burke, K., and Ernzerhof, M. (1996). Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868. doi: 10.1103/PhysRevLett.77.3865

PubMed Abstract | CrossRef Full Text | Google Scholar

Rietveld, H. M. (1969). A profile refinement method for nuclear and magnetic structures. J. Appl. Crystallogr. 2, 65–71. doi: 10.1107/S0021889869006558

CrossRef Full Text | Google Scholar

Roitero, E., Ochoa, M., Anglada, M., Muecklich, F., and Jimenez-Pique, E. (2018). Low temperature degradation of laser patterned 3y-tzp: enhancement of resistance after thermal treatment. J. Eur. Ceram. Soc. 38, 1742–1749. doi: 10.1016/j.jeurceramsoc.2017.10.044

CrossRef Full Text | Google Scholar

Veselov, D. A., Shashkin, I. S., Bobretsova, Y. K., Bakhvalov, K. V., Lutetskiy, A. V., and Tarasov, I. S. (2016). Study of the pulse characteristics of semiconductor lasers with a broadened waveguide at low temperatures (110-120 K). Semiconductors 50, 1396–1402. doi: 10.1134/S1063782616100249

CrossRef Full Text | Google Scholar

Wang, Z., Zhang, Q., Sun, D., and Yin, S. (2007). Study on thermal expansion of Nd3+:Gd3Ga5O12 laser crystal. J. Rare Earths 25, 244–246. doi: 10.1016/S1002-0721(07)60480-3

CrossRef Full Text | Google Scholar

Wynne, R., Daneu, J. L., and Fan, T. Y. (1999). Thermal coefficients of the expansion and refractive index in YAG. Appl. Opt. 38, 3282–3284. doi: 10.1364/AO.38.003282

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Yue, Y., Yang, F., Hu, Z., and Xu, Z. (2016). 266 nm ultraviolet light generation in Ga-doped BaAlBO3F2 crystals. Opt. Lett. 41, 1598–1600. doi: 10.1364/OL.41.001598

CrossRef Full Text | Google Scholar

Yue, Y., Hu, Z., Zhou, Y., Wang, J., Zhang, X., Chen, C., et al. (2011). Growth and nonlinear optical properties of BaAlBO3F2 crystal. J. Opt. Soc. Am. B: Opt. Phys. 28, 861–866. doi: 10.1364/JOSAB.28.000861

CrossRef Full Text | Google Scholar

Zhou, Y., Wang, G., Yue, Y., Li, C., Lu, Y., Cui, D., et al. (2009). High-efficiency 355 nm generation in barium aluminum borate diflouride BaAlBO3FF2. Opt. Lett. 34, 746–748. doi: 10.1364/OL.34.000746

CrossRef Full Text | Google Scholar