“…Successive layers of this combination of polyhedra create channels parallel to the c axis of the hexagonal crystal structure. The channels have a maximum width of 5.1 Å (0.51 nm) (Kolesov and Geiger, 2000;Krambrock et al, 2002;Mashkovtsev et al, 2016;Fukuda and Shinoda, 2008;Blak et al, 1982) alternating with bottlenecks of 2.8 Å (0.28 nm) (Kolesov and Geiger, 2000;Mashkovtsev et al, 2016;Blak et al, 1982;Fukuda et al, 2009), thus creating two possible crystallographic positions, referred to as 2a and 2b (Hawthorne and Černý, 1977;Aurisicchio et al, 1988;Mashkovtsev and Lebedev, 1993;Artioli et al, 1995;Fukuda and Shinoda, 2008;Mashkovtsev et al, 2016;Andersson, 2019) at the coordinates (0, 0, 0.25) for the 2a site and (0, 0, 0) for the 2b site, respectively (Hawthorne and Černý, 1977;Kolesov and Geiger, 2000;Aurisicchio et al, 1988;Fukuda and Shinoda, 2008;Artioli et al, 1995;Gatta et al, 2006;Fridrichová et al, 2018). In addition to the replacement of constituting ions by ions of similar valence state and radii as is possible in many minerals, the large dimension of the Si tetrahedra channels in beryl provides space for various large molecules and ions, such as H 2 O or Na, K, Ca, and Cs, to be incorporated (Krambrock et al, 2002;Kolesov and Geiger, 2000).…”