Optimization of single crystal growth of candidate quantum spin-ice Pr2Hf2O7 by optical floating-zone method
V. K. Anand, A. T. M. N. Islam, A. Samartzis, J. Xu, N. Casati, and B., Lake

TL;DR
This paper details the optimized growth conditions for high-quality single crystals of the quantum spin-ice candidate Pr2Hf2O7 using an optical floating-zone method, and characterizes their structural and magnetic properties.
Contribution
It introduces optimized growth parameters for Pr2Hf2O7 crystals and provides comprehensive structural and magnetic characterization of these crystals.
Findings
Optimal growth in flowing argon atmosphere yields highest quality crystals.
High-resolution SXRD confirms structural quality of grown crystals.
Magnetic measurements characterize the magnetic behavior of the crystals.
Abstract
We report the optimized conditions for growing the high quality single crystals of candidate quantum spin-ice Pr2Hf2O7 using the optical floating-zone method. Large single crystals of Pr2Hf2O7 have been grown under different growth conditions using a four-mirror type optical floating-zone furnace and their microscopic structural differences have been probed by high-resolution synchrotron x-ray diffraction (SXRD). The SXRD data reveal that the crystals grown under fowing argon (~ 2 L/h) atmosphere with slightly off-stoichiometric (optimized) starting composition yields the highest quality crystals. The magnetic susceptibility, isothermal magnetization and heat capacity data of optimally grown crystals are presented.
| crystal (# X2) | polycrystalline | |
|---|---|---|
| Lattice parameters | ||
| (Å) | 10.6727(1) | 10.6728(1) |
| (Å3) | 1215.70(1) | 1215.72(1) |
| Atomic coordinate | ||
| 0.3353(3) | 0.3351(4) | |
| Refinement quality | ||
| (%) | 8.49 | |
| (%) | 10.1 | |
| (%) | 4.57 |
| direction | H | Fit range | ||||
|---|---|---|---|---|---|---|
| (T) | (K) | (emu/mol Pr) | (emu K/mol Pr) | (K) | /Pr | |
| 0.1 | K | 0.79(1) | 2.51 | |||
| 0.1 | K | 0.79(1) | 2.51 | |||
| 0.1 | K | 0.79(1) | 2.51 | |||
| 1.0 | K | 1.22(1) | 3.12 | |||
| 1.0 | K | 1.22(1) | 3.12 | |||
| 1.0 | K | 1.22(1) | 3.12 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Optimization of single crystal growth of candidate quantum spin-ice by optical floating-zone method
V. K. Anand
A. T. M. N. Islam
Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Hahn-Meitner Platz 1, D-14109 Berlin, Germany
A. Samartzis
J. Xu
Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Hahn-Meitner Platz 1, D-14109 Berlin, Germany
Institut für Festkörperphysik, Technische Universität Berlin, Hardenbergstraße 36, D-10623 Berlin, Germany
N. Casati
Paul Scherrer Institute, 5232 Villigen PSI, Switzerland
B. Lake
Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Hahn-Meitner Platz 1, D-14109 Berlin, Germany
Institut für Festkörperphysik, Technische Universität Berlin, Hardenbergstraße 36, D-10623 Berlin, Germany
Abstract
We report the optimized conditions for growing the high quality single crystals of candidate quantum spin-ice Pr2Hf2O7 using the optical floating-zone method. Large single crystals of Pr2Hf2O7 have been grown under different growth conditions using a four-mirror type optical floating-zone furnace and their microscopic structural differences have been probed by high-resolution synchrotron x-ray diffraction (SXRD). The SXRD data reveal that the crystals grown under flowing argon ( L/h) atmosphere with slightly off-stoichiometric (optimized) starting composition yields the highest quality crystals. The magnetic susceptibility, isothermal magnetization and heat capacity data of optimally grown crystals are presented.
I INTRODUCTION
The 227 rare earth pyrochlores have been a topic of interest owing to their exotic properties, for example Dy2Ti2O7 and Ho2Ti2O7 present spin-ice behavior for which magnetic monopoles are believed to be the fundamental excitations Gardner2010 ; Castelnovo2012 ; Gingras2014 ; Ramirez1999 ; Siddharthan1999 ; Hertog2000 ; Bramwell2001 ; Castelnovo2008 . In these materials rare earth ions form a three dimensional network of corner-sharing tetrahedra where ferromagnetic interaction becomes frustrated under the action of the crystal electric field (CEF) that forces the spins to point along the local direction. The CEF induced Ising anisotropic ground state of the classical spin-ice materials consists of ‘two-in/two-out’ spin configuration referred as ‘ice rule’ Harris1997 . Depending upon the relative strengths of CEF anisotropy and the nearest neighbor dipolar and exchange interactions these materials exhibit a wide range of magnetic ground states. For example, spin-ice behavior when ferromagnetic dipolar interaction dominates over the antiferromagnetic exchange interaction as in Dy2Ti2O7 and Ho2Ti2O7 Ramirez1999 ; Siddharthan1999 ; Hertog2000 ; Bramwell2001 , and all-in/all-out antiferromagnetic order when antiferromagnetic exchange interaction dominates such as in Nd2Zr2O7 Lhotel2015 ; Xu2015 and Nd2Sn2O7 Bertin2015 .
In our efforts to extend work on rare earth pyrochlores, we recently started working on hafnate pyrochlores Anand2015 ; Anand2017 ; Anand2016 which remained poorly investigated. Our investigations revealed a long-range antiferromagnetic ordering in Nd2Hf2O7 ( K) with an all-in/all-out arrangement of Nd3+ moments determined by the neutron diffraction study Anand2015 . The neutron diffraction data showed a strongly reduced ordered moment of 0.62(1) /Nd at 0.1 K, much lower than the expected 2.5 /Nd for the Ising anisotropic ground state of Nd2Hf2O7. Persistent dynamic spin fluctuations deep inside the ordered state were found in the muon spin relaxation study which seems to be responsible for the strongly reduced ordered moment Anand2017 . The inelastic neutron scattering study revealed the Kramers doublet ground state of Nd2Hf2O7 to be of dipolar-octupolar type Anand2017 .
Very recently the hafnate pyrochlore Pr2Hf2O7 was identified to be a potential candidate quantum spin-ice system Anand2016 ; Sibille2016 . Pr2Hf2O7 was found to show no long-range magnetic ordering down to 90 mK, however the frequency dependent behavior of ac magnetic susceptibility clearly reflects slow spin dynamics and spin-freezing behavior Anand2016 ; Sibille2016 . While the magnetic ground state of Pr2Hf2O7 is Ising anisotropic in nature, the presence of a weak non-Ising contribution is also inferred from the isothermal magnetization data Anand2016 thus favoring dynamic spin-ice behavior in Pr2Hf2O7. Further, Pr2Hf2O7 has a non-Kramers doublet ground state which is considered favorable for the dynamic spin-ice behavior Lee2012 .
In order to probe and understand the intrinsic properties of Pr2Hf2O7 in more detail it is highly desired to have a high quality single crystal of this compound. Because of the higher melting point of hafnate pyrochlores (compared to titanate pyrochlores) the synthesis of large single crystals of Pr2Hf2O7 using optical floating-zone method turns out to be somewhat tricky. While one can grow titanate pyrochlores in air at ambient pressure, the same recipe does not work for Pr2Hf2O7 which tends to decompose to Pr6O11 and HfO2. As such one needs to grow Pr2Hf2O7 under a pressurized/argon atmosphere. Very recently Ciomaga Hatnean et al. Hatnean2017 also reported a successful growth of Pr2Hf2O7 single crystal under pressure. We find that the quality of crystals not only depends on the atmospheric conditions but also on the stoichiometry of the starting material. We have optimized the growth conditions to yield a high quality crystal as report here. We used high-resolution synchrotron x-ray diffraction (SXRD) measurements to distinguish the subtle structural differences of crystals grown under different growth conditions. We also present the magnetic susceptibility , isothermal magnetization and heat capacity data of single crystal grown under optimized condition. The Ising anisotropy is confirmed, and the bulk properties are found to be consistent with those of polycrystalline Pr2Hf2O7 Anand2016 .
II EXPERIMENTAL DETAILS
Single crystal growth of Pr2Hf2O7 was carried out by the Floating-zone technique using a four-mirror type high-temperature optical floating zone furnace (Crystal Systems corp. FZ-T-12000-X-VII-VPO-PC) equipped with four 3 kW Xenon arc lamps at the Core Lab for Quantum Materials, Helmholtz-Zentrum Berlin (HZB). Crystals were grown under high purity argon atmosphere with a growth rate of 1–5 mm/h. The average power required for the growth was 36 % of lamp power (four 3 kW lamps). During the growth the seed and feed rods were counter-rotated at 5 and 20 rpm, respectively. The seed and feed rods were prepared from high purity materials: Pr6O11 (99.99%, Alfa Aesar) and HfO2 (99.95%, Alfa Aesar) mixed and ground thoroughly and fired at 1300 ∘C for 15 h, reground and fired again at 1300 ∘C for 24 h, then powdered and pressed hydrostatically in the form of cylindrical rod. The cylindical rods were then sintered for 60 h at 1500 ∘C. The sintered rods were typically about 5 mm in diameter and 8-10 cm in length. All grindings and heat treatments were done in air.
For subsequent growths we used single crystals from previous growths as seed crystals to avoid random nucleation and to obtain one large single crystal. Laue diffraction was used to check the quality of the single crystals and orient along specific crystallographic direction. The oriented crystal was then cut into small pieces for magnetic measurements. The crystal quality was also examined by polarized optical microscopy. The chemical composition of crystals were checked by energy dispersive x-ray (EDX) analysis using a scanning electron microscope.
Further, the quality of the single crystals were probed by high-resolution synchrotron x-ray diffraction measurements on crushed crystals performed on the MS-beamline Willmott2013 at Paul Scherrer Institute (PSI), Switzerland. The crystals were very finely ground and filled into a thin glass capillary of 0.3 mm diameter. During the measurement the sample was rotated continuously to reduce the effect of preferred orientation and profile shape dependence. Measurements were performed with incident x-rays of 25 keV energy at a temperature of 300 K. In order to accurately determine the wavelength of the x-rays (0.4148 Å) and the instrument profile parameters, the standard LaB6 powder (NIST) was also measured under the identical experimental condition.
The dc magnetic susceptibility versus temperature and isothermal magnetization versus magnetic field measurements were performed using a superconducting quantum interference device vibrating sample magnetometer (SQUID-VSM, Quantum Design, Inc.) and VSM option of a physical properties measurement system (PPMS, Quantum Design, Inc.) at Core Lab for Quantum Materials, HZB. The heat capacity was measured by the adiabatic-relaxation technique using the heat capacity option of PPMS at Core Lab for Quantum Materials, HZB.
III Crystal Growth
The optimization of the crystal growth of Pr2Hf2O7 was done by varying the growth conditions in several subsequent growths. According to the phase diagram of Pr2O3–HfO2 system Kravchinskaya1978 ; Shevchenko1984 Pr2Hf2O7 melts incongruently around 2420 ∘C in vacuum. However, the melting behavior in other atmospheric environments remains unclear. To check the melting behavior of the material in air, in the first growth attempt we used a feed-rod having stoichiometric composition and an ambient air atmosphere. During the growth we observed evaporation of large quantity of black fumes. The black evaporated powders condensed on the inner surface of the quartz tube and made the growth impossible to continue after a while. The deposit inside the quartz tube was found to be Pr6O11, which suggests that Pr3+ is oxidizing to higher oxidation state. This oxidizing behavior suggests that air does not serve as a good atmosphere, thus urging for an inert atmosphere to be more suitable for the growth. From this growth attempt we could also assume that Pr2Hf2O7 melts incongruently in air atmosphere. Therefore crystal growth using a self-flux, that can reduce the melting temperature, was thought be more suitable.
In our second attempt the growth was done by the traveling solvent floating zone technique, using a solvent with excess Pr2O3 (60 mol% Pr2O3 and 40 mol% Hf2O4). With excess Pr2O3 the eutectic melting point is lowered which allows growth at a relatively lower temperature. In order to compensate for the loss due to evaporation, we used a feed rod with 3% excess Pr2O3 (3% excess Pr6O11 was added to the starting composition of feed rod). The growth atmosphere was changed to argon in order to prevent Pr2O3 from oxidizing to higher oxidized state. Crystal growth was done at a rate of 1.2 mm/h. Under these conditions we observed significant reduction in evaporation during growth. Single crystal of Pr2Hf2O7 could be grown by this solvent floating zone growth process, though the crystal was not of high quality. Further, the visual observation of the melting behavior of the feed-rod suggested that the compound possibly melts congruently in Ar atmosphere.
That the Pr2Hf2O7 may be congruently melting in Ar atmosphere, prompted us to try floating zone growths (without traveling solvent) so that the growth rate can be increased, and large crystals can be grown in a relatively shorter time. Indeed, in Ar atmosphere we could successfully grow large single crystals by the floating zone method using the feed rods prepared from both stoichiometric and off-stoichiometric starting compositions. Subsequently, the crystals were grown under flowing argon atmosphere (2 L/h) by the floating zone technique. Crystal growths were done at a relatively higher speed of 5 mm/h. Large single crystals of about 5 mm in diameter and 4–8 cm in lengths could be successfully grown by this process.
We noticed that the appearance of the crystal depends on the stoichiometry of the starting material. We will be focussing here mainly on two crystals, one grown with a stoichiometric feed rod (crystal X1), and the other grown (optimally) with an offstoichiometric feed rod (crystal X2). The crystal X1 grown with stoichiometric feed rod (growth rate 5 mm/h, Ar amosphere) was bottle green in color [Fig. 1(a)]. The crystal X2 that was grown with feed rod having 3% excess Pr2O3 (growth rate 5 mm/h, Ar amosphere) was fluorescent green in color [see Fig. 1(a)]. For crystal X2 the feed rod was prepared with the starting materials taken as: Pr6O11 + HfO2 = Pr2Hf2O7 + Pr2O3) + O2. Growth conditions for crystal X2 was found to be optimum in terms of crystal quality and stoichiometry.
The as-grown crystal X1 and a piece of crystal X2 are shown in Fig. 1(a). Despite a clear change in the appearance of crystals, we could not notice any difference in the quality of crystals from X-ray Laue photographs. The representative Laue photographs for the three principal crystallographic directions [100], [110] and [111] of optimally grown crystal X2 are shown in Fig. 1(b–d)). No difference was noticed in the laboratory based powder X-ray diffraction (data not shown) performed on the crushed single crystals X1 and X2. We therefore performed synchrotron X-ray diffraction measurements to probe the microscopic details of crystals. The cross sections of the single crystals were examined by polarized optical microscopy and found to be free of grain boundaries or inclusions. The EDX composition analysis of single crystals revealed the stoichiometric ratio for Pr and Hf for both the crystals. However, due to the experimental limitations of EDX analyzer it was not possible to obtain any reliable estimate for light element oxygen.
IV Crystallography
Pr2Hf2O7 forms in face-centered cubic (space group ) pyrochlore structure Karthik2012 ; Anand2016 . The crystal structure is illustrated in Fig. 2(a). In this pyrochlore structure, while Pr and Hf atoms occupy only one site: 16d (1/2,1/2,1/2) and 16c (0,0,0), respectively, the O atoms occupy two sites: 48f (,1/8,1/8) and 8b (3/8,3/8,3/8), designated as O1 and O2, respectively Subramanian1983 . The formula unit can therefore be viewed as Pr2Hf2O16O2. The Pr(16d) atoms are surrounded by 6 O1(48f) atoms and 2 O2(8b) atoms, and are thus eight-fold coordinated Subramanian1983 . The coordination unit is a scalenohedra PrO16O22 and the Pr atoms sit at the centre of a distorted cube formed by 6 O1 and 2 O2, the distance Pr–O2 being slightly smaller than Pr–O1. The Hf(16c) atoms are surrounded by 6 O1(48f) atoms, and are thus six-fold cordinated Subramanian1983 . The coordination unit is a trigonal antiprism HfO16 and the Hf atoms sit at the centre of a distorted octahedra formed by 6 equally distant O1 atoms. Both Pr(16d) and Hf(16c) separately form three-dimensional networks of corner-sharing tetrahedra. The corner-sharing tetrahedra formed by Pr atoms is illustrated in Fig. 2(b).
The room temperature high-resolution synchrotron XRD pattern for powdered Pr2Hf2O7 single crystal X2 is shown in Fig. 3. The structural Rietveld refinement profile using the software FullProf Rodriguez1993 is also shown in Fig. 3. The Rietveld refinement confirmed the -type face-centered cubic (space group ) pyrochlore structure. The refinement also reflects the single phase nature of the sample. The crystallographic and refinement quality parameters are listed in Table 1. The lattice parameters agree very well with that of the polycrystalline Pr2Hf2O7 Anand2016 ; Karthik2012 . No Pr/Hf site mixing and/or oxygen deficiency could be detected in our refinement, which if present, is expected to be less than 0.5%.
The microscopic differences in the SXRD patterns of the single crystals X1 and X2 are compared in Fig. 2(d,e). For comparison the SXRD pattern of polycrystalline Pr2Hf2O7 Anand2016 is also presented [Fig. 2(c)]. Peak splitting is clearly seen in the case of crystal X1 [see Fig. 2(d)] which is not there in the case of optimally grown crystal X2 [Fig. 2(e)] or in the case of the polycrystalline sample [Fig. 2(c)]. The peak splitting reflects a distribution of lattice parameters in crystal X1 due to the presence of sub-crystals of slightly different lattice parameters suggesting off-stoichiometry. On the other hand no such distribution of lattice parameter is present in the case of crystal X2 [Fig. 2(e)] for which the pattern is very similar to that of the polycrystalline sample [Fig. 2(c)]. Thus we see that the optimally grown crystal X2 is the best crystal.
V Magnetization and Heat Capacity
The zero-field-cooled magnetic susceptibility data of single crystal Pr2Hf2O7 are shown in Fig. 4 for and 1.0 T applied along the three crystallographic directions , and . As can be seen from the upper insets of Fig. 4 the for all , and are very similar in magnitude and do not show any anomaly at K. The follows modified Curie-Weiss behavior, . We fitted the data in K to get an estimate of and for the Ising anisotropic ground state of Pr2Hf2O7. The fits of the T data by the modified Curie-Weiss law over K are shown in the lower insets of Fig. 4. The parameters obtained from the fits are listed in Table 2. An effective moment of /Pr is obtained from the values of . The obtained is much lower than the free ion value of /Pr for Pr3+ ions, reflecting the presence of strong Ising anisotropy.
We also fitted the high- data with modified Curie-Weiss law yielding the values of and in the paramagnetic state. The fits of T data by the modified Curie-Weiss law in K are shown in Fig. 4. The fit parameters are listed in Table 2. The value gives an effective moment of /Pr which is again smaller than the expected value of /Pr for free Pr3+ ions. The smaller value of could be the result of strong crystal field effect. As the first excited CEF state in Pr2Hf2O7 lies only at 9.1 meV Anand2016 , both the ground state doublet and the excited singlet are thermally occupied in this temperature range. Thus the estimate of may be affected by the thermal population of this CEF level.
The isothermal magnetization data of single crystal Pr2Hf2O7 are shown in Fig. 5 for along the three crystallographic directions, i.e. , and at K, 5 K, 10 K and 20 K. It is seen that the isotherms for all three directions of become nonlinear at lower temperatures reflecting the saturation tendency of magnetization in a mean-field way. A comparison of the isotherms for , and at 1.8 K is shown in Fig. 5(d). The data show clear anisotropy with as expected for Ising anisotropic behavior in 227 pyrochlores. The data follow saturation values of for [100] (two-in/two-out), for [110] (one-in/one-out, two free) and for [111] (three-in/one-out) as expected for a system with Ising anisotropy Fukazawa2002 . As seen from the Fig. 5(d), at T the values of are /Pr, /Pr and /Pr which are consistent with the estimated Ising anisotropic values /Pr, /Pr and /Pr taking , the effective moment. The observed higher value of possibly indicates the presence of non-Ising contributions. Further we also notice that in contrast to the expected saturation behavior for Ising anisotropic spin-ice materials, the for all three directions keeps on increasing up to the measured field of 14 T and this weak field dependence is again a signature for the presence of non-Ising term in the Hamiltonian.
The heat capacity data of Pr2Hf2O7 single crystal for 1.8 K 300 K are shown in Fig. 6. As can be seen from the inset of Fig. 6 the low- data show an upturn at K which is a residual of peak near 1.8 K Sibille2016 associated with the slow spin dynamics. The high- data are well represented by the Debye and Einstein models of lattice heat capacity. The fit of data by Debye and Einstein models of heat capacity as detailed in Ref. 26 plus crystal field contribution deduced from the inelastic neutron scattering (INS) measurement on Pr2Hf2O7 Anand2016 is shown by the solid red curve in Fig. 6. The fit gives K and Einstein temperature K with 66% weight to Debye term and 34% to Einstein term.
We note that the magnetic susceptibility, isothermal magnetization and heat capacity measurements on crystal X1 and X2 did not reveal any noticeable difference within the accuracy of instruments used. The bulk properties of optimized crystal X2 are consistent with those of polycrystalline sample Anand2016 . Our magnetic susceptibility and isothermal magnetization data are also in good agreement with the recently reported data of single crystal Pr2Hf2O7 Hatnean2017 .
VI Summary and Conclusions
We have successfully grown the large single crystals of candidate quantum spin-ice Pr2Hf2O7 using optical floating-zone method and optimized the growth conditions. Further, the microscopic structural differences of crystals grown under different conditions were probed by the high-resolution synchrotron x-ray diffraction which revealed that the crystal grown under the optimized condition is free from sub-lattice formation or site mixing. The best crystal was obtained with starting material having 3% excess Pr6O11, grown at a growth rate of 5 mm/hour under flowing argon atmosphere ( L/h). As the presence of distribution of lattice parameters can affect the quantum spin ice physics, the optimally grown crystals should provide a scope for further study of the intrinsic properties of quantum spin-ice candidate Pr2Hf2O7. The optimally grown crystal was investigated by the magnetic susceptibility, isothermal magnetization and heat capacity measurements. The and data confirm the Ising anisotropy and presence of weak non-Ising contribution as previously observed for polycrystalline Pr2Hf2O7.
Acknowledgements
We acknowledge Helmholtz Gemeinschaft for funding via the Helmholtz Virtual Institute (Project No. VH-VI-521).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) J. S. Gardner, M. J. P. Gingras, and J. E. Greedan, Magnetic pyrochlore oxides , Rev. Mod. Phys. 82 , 53 (2010).
- 2(2) C. Castelnovo, R. Moessner, and S. L. Sondhi, Spin ice, fractionalization, and topological order , Annu. Rev. Condens. Matter Phys. 3 , 35 (2012).
- 3(3) M J P Gingras and P A Mc Clarty, Quantum spin ice: a search for gapless quantum spin liquids in pyrochlore magnets , Rep. Prog. Phys. 77 , 056501 (2014).
- 4(4) A. P. Ramirez, A. Hayashi, R. J. Cava, R. B. Siddharthan, and S. Shastry, Zero-point entropy in ‘spin ice’ , Nature 399 , 333 (1999).
- 5(5) R. Siddharthan, B. S. Shastry, A. P. Ramirez, A. Hayashi, R. J. Cava, and S. Rosenkranz, Ising pyrochlore magnets: Low-temperature properties, “ice rules,” and beyond , Phys. Rev. Lett. 83 , 1854 (1999).
- 6(6) B. C. Den Hertog and M. J. P. Gingras, Dipolar interactions and origin of spin ice in ising pyrochlore magnets , Phys. Rev. Lett. 84 , 3430 (2000).
- 7(7) S. T. Bramwell and M. J. P. Gingras, Spin ice state in frustrated magnetic pyrochlore materials , Science 294 , 1495 (2001).
- 8(8) C. Castelnovo, R. Moessner, and S. L. Sondhi, Magnetic monopoles in spin ice , Nature 451 , 42 (2008).
