research communications
Redetermination of metarossite, CaV5+2O6·2H2O
aDépartement des Sciences de la Terre, École Normale Supérieure de Lyon, Site Monod, 15 parvis René Descartes, BP 7000, 69342 Lyon, France, bDépartement des Sciences de la Terre, Université Claude Bernard Lyon 1, 43 Bd du 11 Novembre 1918, 69622 Villeurbanne Cedex, France, cDepartment of Geosciences, University of Arizona, 1040 E. 4th Street, Tucson, AZ 85721-0077, USA, and dLunar and Planetary Laboratory, University of Arizona, 1629 E. University Blv, Tucson, AZ 85721-0092, USA
*Correspondence e-mail: anais.kobsch@ens-lyon.fr
The 2O6·2H2O [chemical name: calcium divanadium(V) hexaoxide dihydrate], was first determined using precession photographs, with fixed isotropic displacement parameters and without locating the positions of the H atoms, leading to a reliability factor R = 0.11 [Kelsey & Barnes (1960). Can. Mineral. 6, 448–466]. This communication reports a structure redetermination of this mineral on the basis of single-crystal X-ray diffraction data of a natural sample from the Blue Cap mine, San Juan County, Utah, USA (R1 = 0.036). Our study not only confirms the structural topology reported in the previous study, but also makes possible the of all non-H atoms with anisotropic displacement parameters and all H atoms located. The metarossite structure is characterized by chains of edge-sharing [CaO8] polyhedra parallel to [100] that are themselves connected by chains of alternating [VO5] trigonal bipyramids parallel to [010]. The two H2O molecules are bonded to Ca. Analysis of the displacement parameters show that the [VO5] chains librate around [010]. In addition, we measured the Raman spectrum of metarossite and compared it with IR and Raman data previously reported. Moreover, heating of metarossite led to a loss of water, which results in a transformation to the brannerite-type structure, CaV2O6, implying a possible dehydration pathway for the compounds M2+V2O6·xH2O, with M = Cu, Cd, Mg or Mn, and x = 2 or 4.
of metarossite, ideally CaVKeywords: crystal structure; redetermination; metarossite; hydrogen bonds; phase transformation; brannerite.
CCDC reference: 1497229
1. Mineralogical and crystal-chemical context
Metarossite was originally described from Bull Pen Canyon, San Miguel County, Colorado, by Foshag & Hess (1927) as a yellow, platy, soft and friable mineral with composition CaV2O6·2H2O. It is soluble in hot water and generally is formed as a dehydration product of rossite, CaV2O6·4H2O, which itself crystallizes from aqueous solutions (Ahmed & Barnes, 1963).
Barnes & Qurashi (1952) reported triclinic symmetry (P) and unit-cell parameters [a = 6.215 (5), b = 7.065 (5), c = 7.769 (5) Å, α = 92.97 (17), β = 96.65 (17), γ = 105.78 (17)°] of metarossite from a sample from an area near Thompson's, Utah. Later, by means of precession photographs, Kelsey & Barnes (1960) determined its from the material used by Barnes & Qurashi (1952). For structure (R = 0.11), fixed isotropic displacement parameters were introduced without locating the positions of the hydrogen atoms.
This study reports the ) from the Blue Cap mine, San Juan County, Utah, USA, with anisotropic displacement parameters for all non-hydrogen atoms, positions of hydrogen atoms determined, and improvement of the reliability factor to 0.036. Raman spectra were also recorded and compared with that reported in the two studies by Frost et al. (2004, 2005), on a sample from the Burro mine, San Miguel County of Colorado, USA.
of the structure of a metarossite sample (Fig. 12. Structural commentary
The structural topology of metarossite for all non-hydrogen atoms from this study is identical to that reported by Kelsey & Barnes (1960). Chains of edge-sharing distorted [VO5] trigonal bipyramids run parallel to [010], with [V1O5] and [V2O5] polyhedra alternating along the chains (Fig. 2a). These chains are linked by chains of edge-sharing [CaO8] polyhedra aligned parallel to [100] (Fig. 2b). The water molecules are located at three vertices of the [CaO8] polyhedra [OW3, OW8i and OW8ii; symmetry codes: (i) −x + 1, −y + 1, −z; (ii) x + 1, y, z].
It is interesting to note that there is a radial orientation of the displacement ellipsoids associated with the [VO5] chains when viewed along the chain direction (Fig. 3). The amplitude also slightly radially increases, as indicated by the black dashed circles in Fig. 3. We interpret this as the oscillation or libration of the [VO5] chains around [010]. A similar behavior was reported for brackebuschite Pb2Mn3+(VO4)2(OH) (Lafuente & Downs, 2016) where the [Mn3+(VO4)2OH] chains oscillating about an axis.
Numerical data of the hydrogen-bonding scheme in metarossite are presented in Table 1. The bond-valence calculations (Brown, 2002) with the parameters given by Brese & O'Keeffe (1991) confirm that OW3 and OW8 correspond to the two H2O molecules (Table 2). The low bond-valence sum for O5 is because it is an acceptor for three hydrogen atoms (H2, H3 and H4; Table 1). In fact, all acceptor O atoms involved in hydrogen bonding are from VO5 polyhedra, providing the additional linkage between the [CaO8] and [VO5] chains.
3. Raman spectrum
The Raman spectrum of metarossite (Fig. 4) is comparable with the data recorded by Frost et al. (2005) below 1000 cm−1, but is different in the O–H stretching region between 2800 and 3700 cm−1 (Frost et al., 2004). Indeed, they recorded only three Raman bands (at 3177, 3401 and 3473 cm−1), whereas with the present data, it is possible to distinguish four to five bands depending on the orientation (2904, 2954, 3189, 3240 and 3398 cm−1), along with a broad shoulder around 3415–3480 cm−1 (Fig. 4). According to Libowitzky (1999), the band at 3398 cm−1 can be attributed to the OW8–H4 vibration, and the broad shoulder around 3415–3480 cm−1 may correspond to the OW8–H3 and OW3–H2 vibrations (Table 3). The last vibration (OW3–H1) cannot be seen on Fig. 4, but since the frequency currently accepted for free OH− ion is 3560 cm−1 (Lutz, 1995), it can be associated with the IR band at 3526 cm−1 observed by Frost et al. (2004).
4. Synthesis and crystallization
The natural sample used in this study is from the Blue Cap mine, San Juan County of Utah, USA (Fig. 1) and belongs to the RRUFF project collection (https://rruff.info/R100065). Chemical analysis was performed with a CAMECA SX100 electron microprobe operated at 20 kV and 20 nA and a beam size <1 µm. Eight analysis points yielded an average composition (wt.%): CaO 19.2 (1), V2O5 66.6 (4), trace amount of Sr, and H2O 13.06 estimated to provide two H2O molecules per formula unit. The empirical chemical formula, based on eight oxygen atoms, is Ca0.94V5+2.02O6·2H2O. The Raman spectrum of metarossite was collected from a randomly oriented crystal at 50% power of 150 mW on a Thermo–Almega microRaman system, using a solid-state laser with a wavelength of 532 nm and a thermoelectrically cooled CCD detector. The laser was partially polarized with 4 cm−1 resolution and a spot size of 1 µm.
5. Transformation of metarossite
When a small piece of metarossite (edge length in all dimensions 0.1 mm) was placed under a full power laser (150 mW, 532 nm), a change in its Raman spectrum was observed (Fig. 5). In particular, all bands originating from O–H stretching vibrations disappeared, suggesting a complete dehydration of the sample. Moreover, the spectrum below 1200 cm−1 was found to match that of synthetic CaV2O6 (Baran et al., 1987). In addition, we observed similar Raman spectra collected from a metarossite fragment that was heated in air in an oven at 373 K for 12 h. Single crystal X-ray on the heated crystal revealed monoclinic symmetry with parameters a = 10.0 (1), b = 3.6 (2), c = 6.9 (6) Å, β = 105 (6)°, which match those reported for brannerite (Szymanski & Scott, 1982). However, we were unable to obtain more detailed structure information for the heated sample due to its poor crystallinity (caused probably by dehydration).
A number of synthetic metavanadates, such as those with formula M2+V2O6 where M = Cu, Cd, Mg or Mn, are found to be isostructural with brannerite (Baran et al., 1987; Müller-Buschbaum & Kobel, 1991). There are also many hydrated forms of these compounds, including synthetic CuV2O6·2H2O (Leblanc & Ferey, 1990), and CdV2O6·2H2O (Ulická, 1988), as well as natural dickthomssenite MgV2O6·7H2O (Hughes et al., 2001) or ansermetite MnV2O6·4H2O (Brugger et al., 2003). Because tetrahydrated or dihydrated forms of these materials have structures related to rossite or metarossite, it is likely, then, that natural equivalents of the synthetic metavanadates M2+V2O6·xH2O (M = Cu, Cd, Mg or Mn and x = 0, 2 or 4) may exist.
6. details
Crystal data, data collection and structure . The electron microprobe analysis revealed traces of Sr in our sample. The shows a little deficiency for Ca and excess for V. For simplicity, the ideal chemical formula CaV2O6·2H2O was assumed during the Kelsey & Barnes (1960) underline that {101} is often a twin-plane in metarossite, but the crystal used for this X-ray analysis did not show Atomic coordinates of the previous study were taken as starting parameters for The H atoms were located from difference Fourier syntheses and their positions refined with fixed isotropic displacement parameters (Uiso = 0.04 Å2). The maximum residual electron density in the difference Fourier maps was located at 0.86 Å from O7 and the minimum density at 1.39 Å from Ca.
details are summarized in Table 4
|
Supporting information
CCDC reference: 1497229
https://doi.org/10.1107/S2056989016012433/wm5311sup1.cif
contains datablock I. DOI:Structure factors: contains datablock I. DOI: https://doi.org/10.1107/S2056989016012433/wm5311Isup2.hkl
Data collection: APEX2 (Bruker, 2004); cell
SAINT (Bruker, 2004); data reduction: SAINT (Bruker, 2004); program(s) used to solve structure: SHELXT (Sheldrick, 2015a); program(s) used to refine structure: SHELXL2014 (Sheldrick, 2015b); molecular graphics: XtalDraw (Downs & Hall-Wallace, 2003); software used to prepare material for publication: publCIF (Westrip, 2010).CaV2O6·2H2O | Z = 2 |
Mr = 273.99 | F(000) = 268 |
Triclinic, P1 | Dx = 2.814 Mg m−3 |
a = 6.2059 (4) Å | Mo Kα radiation, λ = 0.71073 Å |
b = 7.0635 (4) Å | Cell parameters from 1255 reflections |
c = 7.7516 (5) Å | θ = 2.7–29.8° |
α = 93.166 (4)° | µ = 3.68 mm−1 |
β = 96.548 (4)° | T = 293 K |
γ = 105.883 (4)° | Platy, pale yellow |
V = 323.36 (4) Å3 | 0.07 × 0.07 × 0.06 mm |
Bruker APEXII CCD area-detector diffractometer | 1508 reflections with I > 2σ(I) |
φ and ω scan | Rint = 0.037 |
Absorption correction: multi-scan (SADABS; Bruker, 2004) | θmax = 31.5°, θmin = 2.7° |
Tmin = 0.669, Tmax = 0.746 | h = −9→8 |
5576 measured reflections | k = −10→10 |
2075 independent reflections | l = −11→11 |
Refinement on F2 | Hydrogen site location: difference Fourier map |
Least-squares matrix: full | Only H-atom coordinates refined |
R[F2 > 2σ(F2)] = 0.036 | w = 1/[σ2(Fo2) + (0.0318P)2 + 0.0696P] where P = (Fo2 + 2Fc2)/3 |
wR(F2) = 0.076 | (Δ/σ)max < 0.001 |
S = 1.01 | Δρmax = 0.77 e Å−3 |
2075 reflections | Δρmin = −0.58 e Å−3 |
113 parameters | Extinction correction: SHELXL2014 (Sheldrick, 2015b), Fc*=kFc[1+0.001xFc2λ3/sin(2θ)]-1/4 |
0 restraints | Extinction coefficient: 0.005 (2) |
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. |
x | y | z | Uiso*/Ueq | ||
Ca | 0.76199 (10) | 0.46286 (9) | 0.14933 (8) | 0.01285 (14) | |
V1 | 0.44824 (9) | 0.10219 (7) | 0.33438 (6) | 0.01068 (13) | |
V2 | 0.37329 (8) | 0.58264 (7) | 0.34629 (6) | 0.01004 (13) | |
O1 | 0.4051 (4) | 0.8480 (3) | 0.4170 (3) | 0.0137 (4) | |
O2 | 0.5250 (3) | 0.3869 (3) | 0.3867 (3) | 0.0106 (4) | |
OW3 | 0.8519 (5) | 0.7632 (4) | 0.3671 (4) | 0.0236 (6) | |
O4 | 0.1030 (3) | 0.4726 (3) | 0.3330 (3) | 0.0157 (5) | |
O5 | 0.1896 (4) | 0.0684 (3) | 0.2322 (3) | 0.0207 (5) | |
O6 | 0.6214 (4) | 0.1115 (3) | 0.1910 (3) | 0.0232 (5) | |
O7 | 0.4308 (4) | 0.5998 (3) | 0.1408 (3) | 0.0155 (5) | |
OW8 | 0.0045 (4) | 0.7225 (3) | 0.0059 (3) | 0.0164 (5) | |
H1 | 0.937 (8) | 0.747 (6) | 0.444 (6) | 0.040* | |
H2 | 0.900 (8) | 0.856 (7) | 0.332 (6) | 0.040* | |
H3 | 0.078 (8) | 0.808 (7) | 0.066 (6) | 0.040* | |
H4 | −0.069 (7) | 0.766 (6) | −0.073 (6) | 0.040* |
U11 | U22 | U33 | U12 | U13 | U23 | |
Ca | 0.0120 (3) | 0.0142 (3) | 0.0116 (3) | 0.0024 (2) | 0.0017 (2) | 0.0017 (2) |
V1 | 0.0140 (3) | 0.0077 (2) | 0.0094 (2) | 0.00209 (19) | 0.00032 (18) | 0.00023 (18) |
V2 | 0.0124 (2) | 0.0082 (2) | 0.0087 (2) | 0.0026 (2) | −0.00078 (18) | −0.00002 (18) |
O1 | 0.0216 (11) | 0.0084 (10) | 0.0098 (10) | 0.0030 (9) | −0.0002 (8) | 0.0005 (8) |
O2 | 0.0116 (10) | 0.0092 (10) | 0.0109 (10) | 0.0032 (8) | −0.0001 (8) | 0.0002 (8) |
OW3 | 0.0237 (13) | 0.0201 (13) | 0.0243 (14) | 0.0026 (11) | 0.0011 (10) | 0.0030 (11) |
O4 | 0.0131 (10) | 0.0160 (11) | 0.0162 (11) | 0.0021 (9) | −0.0006 (8) | 0.0013 (9) |
O5 | 0.0211 (12) | 0.0136 (11) | 0.0232 (13) | 0.0024 (9) | −0.0075 (9) | 0.0003 (9) |
O6 | 0.0347 (14) | 0.0147 (11) | 0.0204 (12) | 0.0027 (11) | 0.0151 (10) | 0.0004 (9) |
O7 | 0.0188 (11) | 0.0170 (11) | 0.0113 (10) | 0.0064 (9) | 0.0007 (8) | 0.0020 (9) |
OW8 | 0.0196 (13) | 0.0132 (12) | 0.0130 (11) | 0.0005 (10) | −0.0014 (9) | 0.0006 (9) |
Ca—O7i | 2.381 (2) | V1—O1iii | 1.899 (2) |
Ca—O4ii | 2.394 (2) | V1—O2 | 1.943 (2) |
Ca—OW8ii | 2.441 (2) | V1—O1iv | 2.004 (2) |
Ca—O6 | 2.448 (2) | V1—V1v | 3.1033 (10) |
Ca—O2 | 2.476 (2) | V1—V2iv | 3.1187 (7) |
Ca—O7 | 2.496 (2) | V2—O4 | 1.633 (2) |
Ca—OW3 | 2.530 (3) | V2—O7 | 1.675 (2) |
Ca—OW8i | 2.554 (2) | V2—O1 | 1.874 (2) |
V1—O6 | 1.623 (2) | V2—O2 | 1.895 (2) |
V1—O5 | 1.655 (2) | V2—O2iv | 2.075 (2) |
O7i—Ca—O4ii | 145.60 (8) | O6—Ca—OW8i | 71.80 (7) |
O7i—Ca—OW8ii | 79.03 (8) | O2—Ca—OW8i | 133.65 (7) |
O4ii—Ca—OW8ii | 84.88 (8) | O7—Ca—OW8i | 149.58 (7) |
O7i—Ca—O6 | 89.84 (8) | OW3—Ca—OW8i | 135.12 (8) |
O4ii—Ca—O6 | 88.49 (8) | O6—V1—O5 | 109.11 (13) |
OW8ii—Ca—O6 | 148.85 (8) | O6—V1—O1iii | 105.21 (11) |
O7i—Ca—O2 | 116.63 (7) | O5—V1—O1iii | 97.89 (11) |
O4ii—Ca—O2 | 93.49 (7) | O6—V1—O2 | 94.99 (10) |
OW8ii—Ca—O2 | 145.90 (8) | O5—V1—O2 | 97.89 (10) |
O6—Ca—O2 | 64.79 (7) | O1iii—V1—O2 | 148.57 (9) |
O7i—Ca—O7 | 72.07 (8) | O6—V1—O1iv | 114.83 (11) |
O4ii—Ca—O7 | 140.72 (7) | O5—V1—O1iv | 135.88 (11) |
OW8ii—Ca—O7 | 97.51 (8) | O1iii—V1—O1iv | 74.72 (9) |
O6—Ca—O7 | 106.65 (8) | O2—V1—O1iv | 75.01 (8) |
O2—Ca—O7 | 63.17 (7) | O4—V2—O7 | 106.24 (11) |
O7i—Ca—OW3 | 131.12 (9) | O4—V2—O1 | 105.17 (10) |
O4ii—Ca—OW3 | 72.51 (9) | O7—V2—O1 | 100.84 (10) |
OW8ii—Ca—OW3 | 76.53 (8) | O4—V2—O2 | 106.81 (10) |
O6—Ca—OW3 | 129.97 (8) | O7—V2—O2 | 93.48 (10) |
O2—Ca—OW3 | 70.58 (8) | O1—V2—O2 | 139.47 (9) |
O7—Ca—OW3 | 70.05 (8) | O4—V2—O2iv | 102.14 (10) |
O7i—Ca—OW8i | 77.53 (7) | O7—V2—O2iv | 151.38 (10) |
O4ii—Ca—OW8i | 69.34 (7) | O1—V2—O2iv | 74.79 (8) |
OW8ii—Ca—OW8i | 77.39 (8) | O2—V2—O2iv | 74.69 (9) |
Symmetry codes: (i) −x+1, −y+1, −z; (ii) x+1, y, z; (iii) x, y−1, z; (iv) −x+1, −y+1, −z+1; (v) −x+1, −y, −z+1. |
D—H···A | D—H | H···A | D···A | D—H···A |
OW3—H1···O4iv | 0.78 (5) | 2.37 (5) | 2.965 (4) | 133 (4) |
OW3—H2···O5vi | 0.72 (4) | 2.25 (5) | 2.900 (3) | 150 (5) |
OW8—H3···O5vii | 0.75 (4) | 2.09 (5) | 2.810 (3) | 162 (5) |
OW8—H4···O5viii | 0.84 (4) | 1.97 (5) | 2.794 (3) | 166 (4) |
Symmetry codes: (iv) −x+1, −y+1, −z+1; (vi) x+1, y+1, z; (vii) x, y+1, z; (viii) −x, −y+1, −z. |
O1 | O2 | OW3* | O4 | O5 | O6 | O7 | OW8* | ΣM | |
Ca | 0.253 | 0.218 | 0.315 | 0.272 | 0.240 | 0.205 | 2.107 | ||
0.326 | 0.278 | ||||||||
V1 | 0.581 | 0.687 | 1.491 | 1.624 | 5.153 | ||||
0.770 | |||||||||
V2 | 0.825 | 0.780 | 1.584 | 1.415 | 5.083 | ||||
0.479 | |||||||||
ΣO | 2.176 | 2.200 | 0.218 | 1.899 | 1.491 | 1.897 | 1.981 | 0.482 |
Note: (*) H atoms not considered for calculation. |
This study | Frost et al. (2004) | |||
O—H···O | O···O | ν | ν | O···O |
OW3—H1···O4 | 2.965 | 3504 | 3526 | 2.9393 |
OW3—H2···O5 | 2.900 | 3482 | 3387 | 2.7995 |
OW8—H3···O5 | 2.810 | 3421 | 3181 | 2.6977 |
OW8—H4···O5 | 2.794 | 3404 | 2867 | 2.6227 |
Acknowledgements
This project was the work of summer intern AK, and she thanks the Downs lab at the University of Arizona, the École Normale Supérieure de Lyon and the Université Claude Bernard Lyon 1 for the opportunity. Funding was provided by the École Normale Supérieure de Lyon and the region Auvergne-Rhône-Alpes (France), through the grant Explo'ra Sup.
References
Ahmed, F. R. & Barnes, W. H. (1963). Can. Mineral. 7, 713–726. CAS Google Scholar
Baran, E. J., Cabello, C. I. & Nord, A. G. (1987). J. Raman Spectrosc. 18, 405–407. CrossRef CAS Web of Science Google Scholar
Barnes, W. H. & Qurashi, M. M. (1952). Am. Mineral. 37, 407–422. CAS Google Scholar
Brese, N. E. & O'Keeffe, M. (1991). Acta Cryst. B47, 192–197. CrossRef CAS Web of Science IUCr Journals Google Scholar
Brown, I. D. (2002). In The Chemical Bond in Inorganic Chemistry: The Bond Valence Model. Oxford University Press. Google Scholar
Brugger, J., Berlepsch, P., Meisser, N. & Armbruster, T. (2003). Can. Mineral. 41, 1423–1431. Web of Science CrossRef CAS Google Scholar
Bruker (2004). APEX2, SAINT and SADABS. Bruker AXS Inc., Madison, Wisconsin, USA. Google Scholar
Downs, R. T. & Hall-Wallace, M. (2003). Am. Mineral. 88, 247–250. CrossRef CAS Google Scholar
Foshag, W. F. & Hess, F. L. (1927). Proc. US Natl Museum, 72, 1–12. CrossRef CAS Google Scholar
Frost, R. L., Erickson, K. L. & Weier, M. L. (2004). Spectrochim. Acta A Mol. Biomol. Spectrosc. 60, 2419–2423. Web of Science CrossRef PubMed Google Scholar
Frost, R. L., Erickson, K. L., Weier, M. L. & Carmody, O. (2005). Spectrochim. Acta A Mol. Biomol. Spectrosc. 61, 829–834. Web of Science CrossRef PubMed Google Scholar
Hughes, J. M., Cureton, F. E., Marty, J., Gault, R. A., Gunter, M. E., Campana, C. F., Rakovan, J., Sommer, A. & Brueseke, M. E. (2001). Can. Mineral. 39, 1691–1700. Web of Science CrossRef CAS Google Scholar
Kelsey, C. H. & Barnes, W. H. (1960). Can. Mineral. 6, 448–466. CAS Google Scholar
Lafuente, B. & Downs, R. T. (2016). Acta Cryst. E72, 293–296. Web of Science CSD CrossRef IUCr Journals Google Scholar
Leblanc, M. & Ferey, G. (1990). Acta Cryst. C46, 15–18. CrossRef CAS Web of Science IUCr Journals Google Scholar
Libowitzky, E. (1999). Monatsh. Chem. 130, 1047–1059. Web of Science CrossRef CAS Google Scholar
Lutz, H. D. (1995). Struct. Bond. 82, 85–103. CrossRef CAS Google Scholar
Müller-Buschbaum, H. & Kobel, M. (1991). J. Alloys Compd. 176, 39–46. Google Scholar
Sheldrick, G. M. (2015a). Acta Cryst. A71, 3–8. Web of Science CrossRef IUCr Journals Google Scholar
Sheldrick, G. M. (2015b). Acta Cryst. C71, 3–8. Web of Science CrossRef IUCr Journals Google Scholar
Szymanski, J. T. & Scott, J. D. (1982). Can. Mineral. 20, 271–279. CAS Google Scholar
Ulická, Ľ. (1988). Chem. Pap. 42, 11–19. Google Scholar
Westrip, S. P. (2010). J. Appl. Cryst. 43, 920–925. Web of Science CrossRef CAS IUCr Journals Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.