Potassium titanyl phosphate (KTiOPO4; KTP) has long been recognized as an important non-linear optical (NLO) material (Zumsteg et al., 1976 ▸) due to its unique combination of desirable physical properties including ‘a large hyperpolarizability, excellent temperature window, wide wavelength for phase matching and outstanding crystal stability’ (Stucky et al., 1989 ▸). Work continues to improve the performance of KTP waveguides in optoelectronics (Kores et al., 2021 ▸) and it is finding new uses as a frequency doubler (to 532 nm green light) for 1064 nm Nd–YAG laser radiation in many areas of medicine (Shim & Kim, 2021 ▸; McGarey et al., 2021 ▸). So far as crystal chemistry is concerned, the KTiOPO4 structure type (space group Pna21, a ≃ 12.8, b ≃ 6.4, c ≃ 10.6 Å, Z = 8, Z′ = 2) is remarkably accommodating with respect to partial or complete isovalent or aleovalent substitution at the potassium (Na+, Rb+, Cs+, Tl+, NH4
+…), titanium (ZrIV, HfIV, VIV, SnIV, SbV, Ga3+, Fe3+, Al3+, Cr3+…), phosphorus (AsV, SiIV, GeIV) and even oxygen (OH−, F−) sites and comprehensive reviews on its substitution chemistry have appeared (Sorokina & Voronkova, 2007 ▸).In an attempt to grow single crystals of the possible new KTP analogues NaZrOAsO4 and NaHfOAsO4 by reacting Na2CO3, MO2 (M = Zr, Hf) and As2O5 in a low-melting flux of NaCl and CsCl, the isostructural title compounds CsNa10Zr4(AsO4)9 (I) and CsNa10Hf4(AsO4)9 (II) were the unexpected result and their crystal structures are now described.
Structural commentary
Compounds (I) and (II) are isostructural and crystallize in the rhombohedral space group R
c (No. 167) with an unusually long c unit-cell parameter of nearly 77 Å. This is of course partly a consequence of our choosing the hexagonal (R-centred) setting of the unit cell [the equivalent primitive rhombohedral lattice for (I) has a = b = c ≃ 26.21 Å and α = β = γ ≃ 20.3°] but even so, it is notable that the l index runs well into three figures for (I) in the R-centred setting. This description will focus on the structure of (I) and note significant differences for (II) where applicable.The asymmetric unit of (I), expanded to show the full coordination polyhedra of the zirconium and arsenic atoms, is shown in Fig. 1 ▸. It consists of two zirconium atoms (both with site symmetry 3 on Wyckoff site 12c), two arsenic atoms [As1 on a general position (36f) and As2 with site symmetry 2 (18e)] and six oxygen atoms, one of which is disordered over two adjacent sites (all lying on general positions, 36f), which leads to the unusual 4:9 stoichiometry for the ZrIV and AsO4
3– moieties with a net charge of −11. The structure of (I) is completed by a Cs+ ion (site symmetry
, 6b) and four partly occupied sodium cations [one on a general position (36f), one with site symmetry 2 (18e) and two with site symmetry 3 (12c)]. To maintain charge balance, the four sodium ions must have a total occupancy of 10 based on Z = 6 (full occupancy of the four sites would give 13 sodium ions per caesium ion).
Figure 1
The asymmetric unit of (I) expanded to include the full Zr and As coordination polyhedra showing 50% displacement ellipsoids. Only one disorder component about As2 is shown. Symmetry codes: (i) 1 − x, 1 − y, −z; (ii) y, 1 − x − y, −z; (iii) 1 − y, 1 + x − y, z; (iv)
− x,
− x + y,
− z; (v) x − y, x, −z; (vi) y − x, 1 − x, z.
Both zirconium atoms adopt almost regular ZrO6 octahedral geometries (Müller-Buschbaum, 2010 ▸) when crystal symmetry is taken into account: the mean Zr1—O separation (to 3 × O3 and 3 × O5) is 2.070 Å and the quadratic elongation and angular variance are 1.001 and 4.43°2, respectively (Robinson et al., 1971 ▸). Equivalent data for Zr2 (bonded to to 3 × O2 and 3 × O4) are 2.072 Å, 1.003 and 9.86°2, respectively. The ‘extrapolated’ (Brese & O’Keeffe, 1991 ▸) bond-valence sums (BVS) in valence units are 4.10 and 4.07 for Zr1 and Zr2, respectively, in acceptable agreement with the expected value of 4.00. The mean Hf—O distances in (II) are 2.062 Å for Hf1 (BVS = 4.13, quadratic elongation = 1.002, angular variance = 5.38°2) and 2.065 Å for Hf2 (4.10, 1.004, 13.20°2). It may be seen that the Hf—O bonds are slightly shorter than the Zr—O bonds, which is in accordance with ionic radii data (Shannon, 1976 ▸): r
6(ZrIV) = 0.72 (6 = six-coordinate) and r
6(HfIV) = 0.71 Å and is presumed to arise from the lanthanide contraction effect.The As1 atom in (I) is surrounded by four oxygen atoms (O1–O4) in the geometry of a slightly distorted tetrahedron [mean As—O = 1.677 Å, spread of O—As—O angles = 103.0 (2)–114.9 (2)°, τ4 (Yang et al., 2007 ▸) = 0.95]. Atom As2 is also tetrahedral (to 2 × O5 and 2 × O6), with the latter O atom disordered over two adjacent sites in almost equal occupancies of 0.45 (3):0.55 (3) [O6A⋯O6B = 0.909 (13) Å]. Of the six oxygen atoms in the structure of (I), four of them (O2–O5) bridge zirconium and arsenic atoms with a mean Zr—O—As bond angle of 141.5° [equivalent mean Hf—O—As bond angle in (II) = 140.4°] and two (O1 and O6) are ‘terminal’ and only bonded to arsenic: all of the O atoms also form one or more bonds to nearby caesium and/or sodium ions.The caesium ion in (I) adopts a grossly squashed octahedral coordination to six O1 atoms with Cs1—O1 = 3.235 (4) Å: the cis O—Cs—O bond angles are compressed to 62.30 (10) or expanded to 117.70 (10)°: the Cs1 BVS of 0.61 compared to an expected value of 1.00 suggests significant underbonding. The interpretation of the sodium-ion coordination polyhedra are complicated by the positional disorder of atom O6 but can be described as distorted trigonal bipyramidal (Na1), very distorted tetrahedral (Na2), square-based pyramidal (Na3) and squashed trigonal pyramidal (Na4). It is notable that Na4 is only three coordinate but similar NaO3 geometries have been observed in dehydrated sodium aluminosilicate zeolites (Adams et al., 1982 ▸).The extended structure of (I) (Fig. 2 ▸) can be conceptually broken down into two different types of layers lying parallel to (001). The first layer (type ‘A’) occurs at z ≃ 0, 1/6, 1/3, 1/2, 2/3 and 5/6 with adjacent A-layers laterally displaced by 1/3 in x and 2/3 in y and consists of the Zr2 and As1 centred polyhedra as well as the caesium ions. Fig. 3 ▸ shows that each Zr2O6 octahedron is connected by two As1O4 tetrahedra (via O2 and O4) to result in a ‘honeycomb’ array of polyhedral 12-rings (six octahedra and 12 tetrahedra) encapsulating the Cs+ ions. Atom O3 of the arsenate group provides the link to the type ‘B’ layers on either side of the A layer. This inter-octahedral connectivity via O3 leads to a distinctive ‘lantern’ motif (Fig. 4 ▸) in which three tetrahedra link two octahedra [Zr1⋯Zr2 = 4.886 (2); Hf1⋯Hf2 in (II) = 4.863 (2) Å]: similar ‘lanterns’ are a feature of the polyhedral connectivity in the scandium tungstate [M
2(XO4)3] (Abrahams & Bernstein, 1966 ▸), Nasicon [AM
2(XO4)3] (Anantharamulu et al., 2011 ▸) and langbeinite [A
2
M
2(XO4)3] (Norberg, 2002 ▸) structure types but they differ from (I) because all the vertices of the constituent tetrahedra in these structures link to adjacent octahedra, hence their 2:3 M:X ratios compared to the 4:9 ratio for (I).
Figure 2
The unit-cell of (I) in polyhedral representation viewed approximately down [110]. A single O atom at the average location of O6A and O6B in the asymmetric unit has been used to construct the As2 tetrahedron. Colour code: Zr1O6 octahedra blue, Zr2O6 octahedra green, As1O4 tetrahedra peach, As2O4 tetrahedra rose, Cs sky blue, Na yellow, O (polyhedral corners) red.
Figure 3
View down [001] of an ‘A’-type layer in the structure of (I) in polyhedral representation. Atom and polyhedron colours as in Fig. 2 ▸ except O3 is blue.
Figure 4
Detail of the extended structure of (I) showing a Zr2As3O18 ‘lantern’ motif of Zr1 and Zr2 octahedra linked by three As1 tetrahedra via atoms O2 and O3. In (I), this motif has crystallographically imposed threefold symmetry about a rotation axis passing through the zirconium atoms. Symmetry codes: (i) 1 − y, 1 + x − y, z; (ii) y − x, 1 − x, z.
The B layers in (I) (Fig. 5 ▸) lie at z ≃ 1/12, 1/4, 5/12, 7/12, 3/4 and 11/12 and are associated with the Zr1 and As2 species. These also feature polyhedral 12-rings (six octahedra and six tetrahedra) but only one As2 tetrahedron (with two terminal As2—O6 bonds) links adjacent Zr1 octahedra via atom O5. There are numerous sodium sites associated with the B layers. The disorder of the sodium ions in the vicinities of the B layers and possible small [110] channels (see Fig. 2 ▸) suggests the possibility of ionic conductivity (Norberg, 2002 ▸). An analysis of the stucture with PLATON (Spek, 2020 ▸) with the sodium ions removed indicated that there was 119.4 Å3 of free space per unit cell (∼2.1%).
Figure 5
View down [001] of a ‘B’-type layer in the structure of (I) in polyhedral representation. Atom and polyhedron colours as in Fig. 2 ▸ except O3 is blue.
Database survey
A survey of the Inorganic Crystal Structure Database (ICSD) (Belsky et al., 2002 ▸) revealed 11 matches for crystal structures containing Zr + As + O, the majority of these being Nasicon (Anantharamulu et al., 2011 ▸) derivatives such as NaZr2(AsO4)3 (Chakir et al., 2003 ▸) or KZr2(AsO4)3 (Elbrahimi & Durand, 1990 ▸) as well as one KTP analogue, viz. RbZrOAsO4 (Simpson & Harrison, 2004 ▸). There were no hits for the combination of Hf + As + O.
Synthesis and crystallization
Compound (I) was prepared by mixing 1.00 g of Na2CO3, 0.581 g of ZrO2 and 1.399 g of As2O5 (Na:Zr:As molar ratio ≃ 4:1:3) in an agate mortar: 1.00 g of this mixture was added to 3.0 g of a eutectic-melt mixture (T
melt ≃ 500°C) of NaCl/CsCl (∼0.35:0.65 mol) and placed in a flat-bottom alumina crucible. The crucible was rapidly heated to 500°C in a muffle furnace and then ramped at 12°C min−1 to 700°C and cooled at the same rate to 400°C and then removed from the furnace and left to cool. The gummy white product was washed with copious amounts of hot water followed by acetone to result in a mass of tiny colourless rods of (I). Compound (II) was made in the same way starting from a pre-mixture of 1.00 g Na2CO3, 1.12 g HfO2 and 1.57 g As2O5 and tiny colourless rods of (II) were the result.Caution! Arsenic compounds are highly toxic and carcinogenic. Take all appropriate safety precautions, especially with respect to dust contamination.
Refinement
Crystal data, data collection and structure refinement details are summarized in Table 1 ▸. The crystal chosen for data collection for (I) was found to be twinned over its rhombohedral obverse and reverse settings (Herbst-Irmer & Sheldrick, 2002 ▸) in a 0.797 (3):0.203 (3) ratio, which was processed as a SHELXL HKLF 5 refinement. To ensure charge balance, the occupancies of the four partially occupied sodium sites must sum to 10.0 Na per caesium ion and this was achieved by using a SUMP card (linear free variable restraint) in SHELXL, as unrestrained refinements tended to drift to a collective occupancy of above 10 (full occupancy of the four sodium sites would give 13 Na to 1 Cs). This needed cautious damped refinement cycles to begin with, but as the refinement converged, the damping could be removed to give refined fractional site occupancies of Na1 = 0.852 (5), Na2 = 0.860 (9), Na3 = 0.731 (12) and Na4 = 0.423 (11) for (I) and Na1 = 0.887 (7), Na2 = 0.846 (11), Na3 = 0.735 (16) and Na4 = 0.337 (14) for (II). The final difference map for (II) features electron density peaks of ∼2 e Å−3 near some of the sodium ions, perhaps suggesting that they are localizing over split multiple sites at low temperatures, but efforts to model this did not lead to satisfactory refinements. The value of U
eq for Na4 is small, which might indicate partial occupancy of caesium on this site (i.e., a formula of Cs1+xNa10-xHf4(AsO4)9, but attempts to model this were inconclusive.
Table 1
Experimental details
(I)
(II)
Crystal data
Chemical formula
CsNa10Zr4(AsO4)9
CsNa10Hf4(AsO4)9
Mr
1977.97
2327.05
Crystal system, space group
Trigonal, Rc:H
Trigonal, Rc:H
Temperature (K)
293
120
a, c (Å)
9.2218 (5), 76.982 (5)
9.1795 (2), 76.527 (8)
V (Å3)
5669.6 (7)
5584.5 (6)
Z
6
6
Radiation type
Mo Kα
Mo Kα
μ (mm−1)
10.07
20.25
Crystal size (mm)
0.10 × 0.10 × 0.10
0.08 × 0.08 × 0.08
Data collection
Diffractometer
Bruker SMART CCD
Nonius KappaCCD
Absorption correction
Multi-scan (SADABS; Bruker, 1999 ▸)
Multi-scan (SORTAV; Blessing, 1995 ▸)
Tmin, Tmax
0.350, 0.495
0.40, 0.50
No. of measured, independent and observed [I > 2σ(I)] reflections
2288, 2288, 1694
11660, 1434, 1164
Rint
–
0.070
(sin θ/λ)max (Å−1)
0.756
0.650
Refinement
R[F2 > 2σ(F2)], wR(F2), S
0.040, 0.093, 1.00
0.032, 0.078, 1.06
No. of reflections
2288
1434
No. of parameters
112
107
No. of restraints
1
7
Δρmax, Δρmin (e Å−3)
1.85, −1.63
2.57, −2.00
Computer programs: SMART and SAINT (Bruker, 1999 ▸), DENZO/SCALEPACK (Otwinowski & Minor, 1997 ▸), SHELXS and SHELXS97 (Sheldrick, 2008 ▸), SHELXL2018/3 (Sheldrick, 2015 ▸), ORTEP-3 for Windows (Farrugia, 2012 ▸), ATOMS (Dowty, 2005 ▸) and publCIF (Westrip, 2010 ▸).
Crystal structure: contains datablock(s) I, II, global. DOI: 10.1107/S2056989022006338/pk2665sup1.cifStructure factors: contains datablock(s) I. DOI: 10.1107/S2056989022006338/pk2665Isup2.hklStructure factors: contains datablock(s) II. DOI: 10.1107/S2056989022006338/pk2665IIsup3.hklCCDC references: 2179420, 2179419Additional supporting information: crystallographic information; 3D view; checkCIF report
Primary atom site location: structure-invariant direct methods
R[F2 > 2σ(F2)] = 0.040
w = 1/[σ2(Fo2) + (0.0452P)2] where P = (Fo2 + 2Fc2)/3
wR(F2) = 0.093
(Δ/σ)max < 0.001
S = 1.00
Δρmax = 1.85 e Å−3
2288 reflections
Δρmin = −1.63 e Å−3
112 parameters
Geometry. All esds (except the esd in the dihedral angle between two l.s. planes)
are estimated using the full covariance matrix. The cell esds are taken
into account individually in the estimation of esds in distances, angles
and torsion angles; correlations between esds in cell parameters are only
used when they are defined by crystal symmetry. An approximate (isotropic)
treatment of cell esds is used for estimating esds involving l.s. planes.
Refinement. Refined as a 2-component obverse/reverse twin.
Primary atom site location: structure-invariant direct methods
R[F2 > 2σ(F2)] = 0.032
w = 1/[σ2(Fo2) + (0.0359P)2 + 136.7342P] where P = (Fo2 + 2Fc2)/3
wR(F2) = 0.078
(Δ/σ)max = 0.008
S = 1.06
Δρmax = 2.57 e Å−3
1434 reflections
Δρmin = −2.00 e Å−3
107 parameters
Geometry. All esds (except the esd in the dihedral angle between two l.s. planes)
are estimated using the full covariance matrix. The cell esds are taken
into account individually in the estimation of esds in distances, angles
and torsion angles; correlations between esds in cell parameters are only
used when they are defined by crystal symmetry. An approximate (isotropic)
treatment of cell esds is used for estimating esds involving l.s. planes.
Authors: Patrick O McGarey; Alissa Collins; Laura M Dominguez; Gregory R Dion; Charles Blake Simpson Journal: J Voice Date: 2020-03-01 Impact factor: 2.009