research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Sub-micrometer focusing setup for high-pressure crystallography at the Extreme Conditions beamline at PETRA III

crossmark logo

aDeutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany, bBayerisches Geoinstitut, University of Bayreuth, Universitätsstrasse 30, 95440 Bayreuth, Germany, cMaterial Physics and Technology at Extreme Conditions, Laboratory of Crystallography, University of Bayreuth, Universitätsstrasse 30, 95440 Bayreuth, Germany, dCenter for X-ray and Nano Science CXNS, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany, eHelmholtz Imaging Platform, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany, fDepartment Physik, Universität Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany, gRuhr-Universität Bochum, Universitätsstrasse 150, 44801 Bochum, Germany, hLaboratory for Micro- and Nanotechnology, Paul Scherrer Institut, Forschungsstrasse 111, 5232 Villigen-PSI, Switzerland, and iDepartment of Physics, Chemistry and Biology (IFM), Linköping University, Campus Valla, Fysikhuset F310, SE-581 83 Linköping, Sweden
*Correspondence e-mail: konstantin.glazyrin@desy.de

Edited by S. Svensson, Uppsala University, Sweden (Received 13 December 2021; accepted 7 March 2022; online 4 April 2022)

Scientific tasks aimed at decoding and characterizing complex systems and processes at high pressures set new challenges for modern X-ray diffraction instrumentation in terms of X-ray flux, focal spot size and sample positioning. Presented here are new developments at the Extreme Conditions beamline (P02.2, PETRA III, DESY, Germany) that enable considerable improvements in data collection at very high pressures and small scattering volumes. In particular, the focusing of the X-ray beam to the sub-micrometer level is described, and control of the aberrations of the focusing compound refractive lenses is made possible with the implementation of a correcting phase plate. This device provides a significant enhancement of the signal-to-noise ratio by conditioning the beam shape profile at the focal spot. A new sample alignment system with a small sphere of confusion enables single-crystal data collection from grains of micrometer to sub-micrometer dimensions subjected to pressures as high as 200 GPa. The combination of the technical development of the optical path and the sample alignment system contributes to research and gives benefits on various levels, including rapid and accurate diffraction mapping of samples with sub-micrometer resolution at multimegabar pressures.

1. Introduction

The exploration of systems subjected to high-pressure conditions using a combination of various pressurizing devices and a great number of sensing techniques has enabled great progress in the fields of geological and planetary sciences, materials science, solid-state physics and chemistry. In particular, the diamond anvil cell (DAC), where two opposing anvils compress the sample, has been widely used for the characterization of samples in single-crystalline form up to 2 Mbar (1 bar = 100 000 Pa) (Khandarkhaeva et al., 2020[Khandarkhaeva, S., Fedotenko, T., Bykov, M., Bykova, E., Chariton, S., Sedmak, P., Glazyrin, K., Prakapenka, V., Dubrovinskaia, N. & Dubrovinsky, L. (2020). Eur. J. Inorg. Chem. 2020, 2186-2190.]) and as powders at multiple megabars (Dorfman et al., 2012[Dorfman, S. M., Prakapenka, V. B., Meng, Y. & Duffy, T. S. (2012). J. Geophys. Res. 117, B08210.]; Dubrovinskaia et al., 2016[Dubrovinskaia, N., Dubrovinsky, L., Solopova, N. A., Abakumov, A., Turner, S., Hanfland, M., Bykova, E., Bykov, M., Prescher, C., Prakapenka, V. B., Petitgirard, S., Chuvashova, I., Gasharova, B., Mathis, Y., Ershov, P., Snigireva, I. & Snigirev, A. (2016). Sci. Adv. 2, e1600341.]; Dubrovinsky et al., 2015[Dubrovinsky, L., Dubrovinskaia, N., Bykova, E., Bykov, M., Prakapenka, V., Prescher, C., Glazyrin, K., Liermann, H.-P., Hanfland, M., Ekholm, M., Feng, Q., Pourovskii, L. V., Katsnelson, M. I., Wills, J. M. & Abrikosov, I. A. (2015). Nature, 525, 226-229.]; Tateno et al., 2010[Tateno, S., Hirose, K., Ohishi, Y. & Tatsumi, Y. (2010). Science, 330, 359-361.]). More recently, newly developed DAC designs such as toroidal anvils (Dewaele et al., 2018[Dewaele, A., Loubeyre, P., Occelli, F., Marie, O. & Mezouar, M. (2018). Nat. Commun. 9, 2913.]; Jenei et al., 2018[Jenei, Z., O'Bannon, E. F., Weir, S. T., Cynn, H., Lipp, M. J. & Evans, W. J. (2018). Nat. Commun. 9, 3563.]; McMahon, 2018[McMahon, M. (2018). Nat. Mater. 17, 858-859.]) and the double-stage DAC technique (Dubrovinskaia et al., 2016[Dubrovinskaia, N., Dubrovinsky, L., Solopova, N. A., Abakumov, A., Turner, S., Hanfland, M., Bykova, E., Bykov, M., Prescher, C., Prakapenka, V. B., Petitgirard, S., Chuvashova, I., Gasharova, B., Mathis, Y., Ershov, P., Snigireva, I. & Snigirev, A. (2016). Sci. Adv. 2, e1600341.]; Dubrovinsky et al., 2015[Dubrovinsky, L., Dubrovinskaia, N., Bykova, E., Bykov, M., Prakapenka, V., Prescher, C., Glazyrin, K., Liermann, H.-P., Hanfland, M., Ekholm, M., Feng, Q., Pourovskii, L. V., Katsnelson, M. I., Wills, J. M. & Abrikosov, I. A. (2015). Nature, 525, 226-229.], 2012[Dubrovinsky, L., Dubrovinskaia, N., Prakapenka, V. B. & Abakumov, A. M. (2012). Nat. Commun. 3, 1163.]) have reached pressures up to 0.6 and 1 TPa, respectively. In addition, it has become apparent that the new synthesis products in high-pressure chemistry experiments can only be evaluated efficiently if one can map and isolate the grain of interest from the surroundings, e.g. unreacted precursor materials, the pressure medium or other products of the synthesis.

Thus, one of the major challenges of structural studies at conventional and ultra-high pressures is the reduction of sample sizes, which often requires probes with very small beam size and high flux. As X-ray diffraction (XRD) is the most critical and one of the most powerful probes, the demand for sub-micrometer X-ray diffraction capabilities has increased significantly over the last decade. Developments at large-scale facilities are at the forefront of state-of-art studies employing small beams and enable the determination of crystal structures and their evolution under extreme conditions, using either single crystals (single-crystal X-ray diffraction, SCXRD) [e.g. Bykov et al. (2019[Bykov, M., Chariton, S., Fei, H., Fedotenko, T., Aprilis, G., Ponomareva, A. V., Tasnádi, F., Abrikosov, I. A., Merle, B., Feldner, P., Vogel, S., Schnick, W., Prakapenka, V. B., Greenberg, E., Hanfland, M., Pakhomova, A., Liermann, H., Katsura, T., Dubrovinskaia, N. & Dubrovinsky, L. (2019). Nat. Commun. 10, 2994.]), Bykova et al. (2016[Bykova, E., Dubrovinsky, L., Dubrovinskaia, N., Bykov, M., McCammon, C., Ovsyannikov, S. V., Liermann, H., Kupenko, I., Chumakov, A. I., Rüffer, R., Hanfland, M. & Prakapenka, V. (2016). Nat. Commun. 7, 10661.]), Friedrich et al. (2020[Friedrich, A., Collings, I. E., Dziubek, K. F., Fanetti, S., Radacki, K., Ruiz-Fuertes, J., Pellicer-Porres, J., Hanfland, M., Sieh, D., Bini, R., Clark, S. J. & Marder, T. B. (2020). J. Am. Chem. Soc. 142, 18907-18923.]), Laniel et al. (2020[Laniel, D., Winkler, B., Bykova, E., Fedotenko, T., Chariton, S., Milman, V., Bykov, M., Prakapenka, V., Dubrovinsky, L. & Dubrovinskaia, N. (2020). Phys. Rev. B, 102, 134109.])] or powders (powder X-ray diffraction, PXRD, Anzellini et al., 2019[Anzellini, S., Monteseguro, V., Bandiello, E., Dewaele, A., Burakovsky, L. & Errandonea, D. (2019). Sci. Rep. 9, 13034.]; Errandonea et al., 2020[Errandonea, D., Burakovsky, L., Preston, D. L., MacLeod, S. G., Santamaría-Perez, D., Chen, S., Cynn, H., Simak, S. I., McMahon, M. I., Proctor, J. E. & Mezouar, M. (2020). Commun. Mater. 1, 60.]; Kawamura et al., 2002[Kawamura, H., Akahama, Y., Umemoto, S., Takemura, K., Ohishi, Y. & Shimomura, O. (2002). J. Phys. Condens. Matter, 14, 10407-10410.]).

Review of the available literature shows that, in comparison with powder diffraction techniques, high-pressure SCXRD often provides a more reliable pathway for an accurate and precise structure solution, although the process of data acquisition can be more complex. Indeed, sub-micrometer-sized samples, e.g. synthesized in a DAC after laser heating, have to stay illuminated by the X-ray beam upon rotation by ±38° or higher angles. At the same time, the signal of interest should not be spoiled by a possible overlap with other phases, or grains of the same phase but different orientations. The centering of micrometer-sized samples requires the implementation of more precise and accurate sample-positioning systems on dedicated high-pressure beamlines. Thus, only a combination of precise motorization with the sub-micrometer beam capability of the instrument enables successful experiments, significantly improving the quality and usability of the data.

Here, we report the development of a sub-micrometer focusing setup on the general-purpose experiment table of beamline P02.2 at PETRA III, DESY, Hamburg, Germany (Liermann et al., 2015[Liermann, H.-P., Konôpková, Z., Morgenroth, W., Glazyrin, K., Bednarčik, J., McBride, E. E., Petitgirard, S., Delitz, J. T., Wendt, M., Bican, Y., Ehnes, A., Schwark, I., Rothkirch, A., Tischer, M., Heuer, J., Schulte-Schrepping, H., Kracht, T. & Franz, H. (2015). J. Synchrotron Rad. 22, 908-924.]). The new setup has a sub-micrometer focusing achieved by compound refractive lenses (CRLs), supplemented with a phase plate which reduces the spherical aberrations of the CRLs at 25.6 keV. The implementation of a correcting phase plate reduces the beam tails at the focal spot (Seiboth et al., 2017[Seiboth, F., Schropp, A., Scholz, M., Wittwer, F., Rödel, C., Wünsche, M., Ullsperger, T., Nolte, S., Rahomäki, J., Parfeniukas, K., Giakoumidis, S., Vogt, U., Wagner, U., Rau, C., Boesenberg, U., Garrevoet, J., Falkenberg, G., Galtier, E. C., Ja Lee, H., Nagler, B. & Schroer, C. G. (2017). Nat. Commun. 8, 14623.]). This new development significantly improves the capability of the beamline, enabling users to collect data with a beam size at the focal spot of 0.9 µm × 0.9 µm [horizontal (H) × vertical (V) at FWHM] and a photon flux higher than that without the correcting phase plate.

In addition, we describe an improved sample-positioning system with the implementation of a modern air-bearing rotation stage operating in combination with piezo actuator XY stages, enabling a spherical confusion of the order of a micrometer or less and a position accuracy of 0.1–0.2 µm. We demonstrate the capability of the new setup as an instrument by performing a test using micrometer-sized single crystals of CoSb3 and (Mg1.93Fe0.06)(Si1.93,Al0.06)O6 (orthoenstatite) prepared by the focused-ion-beam technique (FIB). Our studies, conducted at pressures above 150 GPa, characterize Fe3O4 and Fe-bearing perovskite Mg0.91(2)Fe0.09(2)SiO3 measured in a multigrain matrix after laser-heating in diamond anvil cells. These state-of-the-art examples highlight the full potential of the setup with application to high-pressure crystallographic studies.

2. Implementation of CRLs and a phase plate

A detailed description of the Extreme Conditions Beamline, P02.2, at the third-generation light source PETRA III can be found in the article by Liermann et al. (2015[Liermann, H.-P., Konôpková, Z., Morgenroth, W., Glazyrin, K., Bednarčik, J., McBride, E. E., Petitgirard, S., Delitz, J. T., Wendt, M., Bican, Y., Ehnes, A., Schwark, I., Rothkirch, A., Tischer, M., Heuer, J., Schulte-Schrepping, H., Kracht, T. & Franz, H. (2015). J. Synchrotron Rad. 22, 908-924.]). In order to improve the focusing capabilities of the beamline, e.g. 2 µm × 2 µm obtained using a Kirkpatrick–Baez (KB) mirror system, 136 2D parabolic Be CRLs with a radius of curvature at their apex of 50 µm and a geometric aperture of 360 µm were installed at around 70.5 m from the source for X-ray diffraction experiments at 25.6 keV. Using the full X-ray acceptance of the CRLs, a focused beam of 3 µm × 1.2 µm (H × V) was achieved at a focal distance of 395 mm. Here and below, the values are reported for the FWHM. In order to achieve a smaller beam size at the focal spot, we reduced the incident beam to 50 µm × 50 µm at the position of the high heat load slits system (SLT) located in the middle of the beamline about 35 m from the source [just before the double-crystal monochromator (DCM) installed at 40.4 m]. Thus, a small central area of the incident X-ray beam was used for focusing, resulting in a focal spot of ∼0.9 µm × 0.9 µm as determined with sharp-edge scans at the sample position (scans of absorbing cylindrical rods of 6 mm diameter made from hardened steel, in the horizontal and vertical directions in transmission mode).

Unfortunately, the size reduction of the incident beam at the position of the SLT also resulted in a significant loss of the focused beam intensity. In order to improve the photon flux at the focal spot, but at the same time keep the beam size to sub-micrometer level, we reduced the spherical aberration of the lenses (Celestre et al., 2020[Celestre, R., Berujon, S., Roth, T., Sanchez del Rio, M. & Barrett, R. (2020). J. Synchrotron Rad. 27, 305-318.]; Seiboth et al., 2016[Seiboth, F., Kahnt, M., Scholz, M., Seyrich, M., Wittwer, F., Garrevoet, J., Falkenberg, G., Schropp, A. & Schroer, C. G. (2016). Proc. SPIE, 9963, 99630P.]) by implementing a customized phase plate downstream from the CRLs. Prior to the phase plate production, a stack of 136 lenses was pre-characterized on the Hard X-ray Micro/Nano-Probe Beamline P06 at PETRA III via ptychography in order to achieve optimal lateral coherence in the horizonal direction at 25.6 keV (Falkenberg et al., 2020[Falkenberg, G., Seiboth, F., Koch, F., Falch, K. V., Schropp, A., Brückner, D. & Garrevoet, J. (2020). Powder Diffr. 35, S34-S37.]; Seiboth et al., 2020[Seiboth, F., Brückner, D., Kahnt, M., Lyubomirskiy, M., Wittwer, F., Dzhigaev, D., Ullsperger, T., Nolte, S., Koch, F., David, C., Garrevoet, J., Falkenberg, G. & Schroer, C. G. (2020). J. Synchrotron Rad. 27, 1121-1130.]). The calculated radial phase shift compensating for spherical aberrations is shown in Fig. 1[link](a). A polymer structure with a height profile that matches the necessary phase correction was fabricated by a two-photon polymerization method. A cross section of the 3D model of the phase plate is shown in Fig. 1[link](b). The structure was created on a silicon nitride membrane with a Nanoscribe Photonic Professional GT using a `dip-in' lithography mode with a 25× objective (numerical aperture NA = 0.8). Fig. 1[link](c) shows a scanning electron microscopy (SEM) image of the phase plate tailored to reduce the characteristic spherical aberration of the Be CRL stack used on P02.2.

[Figure 1]
Figure 1
The correcting phase plate and its design. (a) The phase shift was quantified for the 136 CRLs on the P06 beamline by means of ptychography. (b) The corresponding cross section of a 3D model of the phase plate providing appropriate aberration corrections for this specific CRL stack. (c) An SEM image of the phase plate produced on a silicon nitride substrate.

Considering the setup of P02.2, the phase plate was placed 29 mm downstream from the edge of the lens stack installed in the housing of a V-grooved CRL holder (Lengeler et al., 2005[Lengeler, B., Schroer, C. G., Kuhlmann, M., Benner, B., Günzler, T. F., Kurapova, O., Zontone, F., Snigirev, A. & Snigireva, I. (2005). J. Phys. D Appl. Phys. 38, A218-A222.]). The CRL holder provides an inert He environment protecting the CRLs from oxidation. The phase plate was installed externally. To align the custom-made correcting phase plate to the optical axis of the CRL stack in the horizontal and vertical directions, a motorized positioning system consisting of 2MFA-PP micro-steppers from Newport was mounted on the CRL V-groove housing. A high-resolution scintillator-based Optique Peter microscope (PCO.edge 4.2 CLHS camera, 20× objective, 10 µm thick YAG:Ce scintillator) was used to observe and control the alignment of the phase plate (Fig. 2[link]).

[Figure 2]
Figure 2
The phase plate alignment procedure. The phase plate is located between the CRLs and the microscope installed downstream from the focal spot position. (a) An image of a defocused X-ray beam without the phase plate. (b) The phase plate is moving in the field of view, approaching the central position from the right-hand side. (c) The phase plate is centered and the image is produced by subtracting the background, e.g. panel (a), from the image with the phase plate centered. In the final steps of alignment and with the sharp edge placed at the position of optimal focus, we make a small displacement of the phase plate perpendicular to the incident beam and scan the sharp edge in order to achieve the best beam shape.

The lenses of the CRL stack were numbered during our measurement at P06 and were installed on the beamline in the same order. This step is important, since the phase shift of the correcting phase plate was calculated for a specific lens configuration. After alignment of the CRLs and the correcting phase plate, further optimization steps were carried out in search of a compromise between the photon flux and the beam size at the focal spot. This included iterative adjustments of the SLT opening and the tracking of the resulting beam size at the optimal focal spot. We measured the latter by means of sharp-edge scans at different focal distances. In the end, the incident beam size at the SLT had dimensions of ∼150 µm × 150 µm and the corresponding optimal focal distance was 400–402 mm (the center of the CRL stack to the focal point).

By scanning the X-ray beam at the focal spot using an absorbing edge we measured a transmission curve with a sigmoid-like shape which corresponds to an edge spread function (ESF) (Wang et al., 2016[Wang, Y., Li, Q., Chen, N., Cheng, J.-M., Li, C.-G., Li, H., Long, Q.-H., Shi, J.-S. & Deng, J.-J. (2016). Chin. Phys. C, 40, 086205.], and references within). By derivation of the ESF we produced a linear spread function (LSF) response enabling determination of the focal spot size. In order to assess the situation before and after the installation of the phase plate, and to take advantage of using the latter, we collected data which are presented in Fig. 3[link]. For case (i), the red line depicts the beam shape at the focal spot without the phase plate and with a smaller opening of the SLT (50 µm × 50 µm), while for case (ii) the blue line shows the data with the phase plate installed and a larger opening of the SLT (∼150 µm × 150 µm). Considering case (i), we add that further SLT opening without installation of the phase plate increases the focal spot beyond 1 µm. Raw data for the curves were measured using the same pin-diode acquisition time, so the installation of the phase plate has produced a gain with a factor of ∼2.2. Comparing cases (i) and (ii), we conclude that a similar focal spot size in the horizontal direction was achieved in both cases, including similar or lower beam tail contributions, but installation of the phase plate enabled a higher photon flux.

[Figure 3]
Figure 3
A comparison of the horizontal focal size with (blue data) and without (red data) the correcting phase plate. The curves correspond to the derivative signal of the sharp-edge scans, the linear spread function (LSF). Analysis of the curves confirms the similar size of the X-ray beam in the horizontal direction. Given the same acquisition time per point, the higher values for the blue data indicate higher flux values for the setup when corrected with the phase plate. With the phase plate (blue), one can employ a larger opening of the SLT installed in front of the DCM position. Implementation of the phase plate enables a higher intensity X-ray beam while maintaining a small beam size with similar, or even slightly smaller, beam tails at the focal spot.

In one of the final steps, the sample stack position was readjusted in order to bring the position of the sample stack rotation axis (controlling the sample-to-detector distance, SDD) into focus. Fig. 4[link] illustrates the horizontal and vertical beam profiles in the focal place used for the measurement of the single-crystal data reported below. The optimization of the X-ray optical path included the installation of a Pt-based pinhole, ranging from 15 to 40 µm in diameter depending on the application. The photon flux for the `sub-micrometer' setup at the beam position with a 40 µm pinhole was measured to be around 3.7 × 109 photons s−1 using a calibrated passivated implanted planar silicon (PIPS) diode.

[Figure 4]
Figure 4
First-order derivatives of sharp-edge scans (LSF) collected during (left) horizontal and (right) vertical motor movement. The FWHM is shown together with the corresponding error bar estimates.

The opening of the SLT does have several effects on the X-ray optical path. On one hand, it controls the total energy dissipating at the DCM crystals (e.g. heat bump, non-ideal surface of crystals). On the other hand, the closing of the SLT effectively introduces a secondary source. It is hard to give a quantitative description of the realistic and non-ideal system, an analysis of which goes beyond the current publication. However, we note that, after further opening of the SLT, an increase in flux by at least of a factor 3–5 was observed. At the same time, the beam size at the focal spot increased to 1.3 µm × 1 µm (H × V) FWHM. This option may be of interest to user groups working with low-Z materials, where the compromise between the illuminating beam intensity and the beam size at the focal spot is shifted towards the higher photon flux at the sample position.

3. Sample positioning system with sub-micrometer resolution

Accurate sample positioning is a crucial aspect of a high-pressure DAC experimental setup. Along with the improved focusing, it was necessary to update the existing sample positioning system (Liermann et al., 2015[Liermann, H.-P., Konôpková, Z., Morgenroth, W., Glazyrin, K., Bednarčik, J., McBride, E. E., Petitgirard, S., Delitz, J. T., Wendt, M., Bican, Y., Ehnes, A., Schwark, I., Rothkirch, A., Tischer, M., Heuer, J., Schulte-Schrepping, H., Kracht, T. & Franz, H. (2015). J. Synchrotron Rad. 22, 908-924.]). We indicate the new parts in Fig. 5[link], labelled (1) to (5). An XY PILine piezo positioning stage (1) and a PIglide RM air-bearing rotation stage (2) were acquired from Physik Instrumente (PI) GmbH & Co. They replace the top portion of the standard sample six-motor positioning system (Liermann et al., 2015[Liermann, H.-P., Konôpková, Z., Morgenroth, W., Glazyrin, K., Bednarčik, J., McBride, E. E., Petitgirard, S., Delitz, J. T., Wendt, M., Bican, Y., Ehnes, A., Schwark, I., Rothkirch, A., Tischer, M., Heuer, J., Schulte-Schrepping, H., Kracht, T. & Franz, H. (2015). J. Synchrotron Rad. 22, 908-924.]). In addition, we exchanged the motor (3) that controls the position of the rotation axis and moves it in the horizontal plane perpendicular to the X-ray beam. It was upgraded through an HPS-170 high-precision linear stage also manufactured by PI. The factory characteristics of these positioning units are listed in Table 1[link].

Table 1
Factory characteristics of the new positioning parts used for the sub-micrometer setup on the Extreme Conditions Beamline (ECB)

The numbers in brackets refer to components in Fig. 5[link].

(1) XY stage with PILine piezo motors 25 mm × 25 mm travel range
100 nm min. incremental motion
10 nm sensor resolution
7 N drive force, maximum load capacity 50 N
 
(2) PIglide RM air-bearing rotation stage 150 mm motion platform diameter
50 mm travel length
0.0015 µrad sensor resolution
Absolute angle-measuring system
Slot-less brushless three-phase torque motor
 
(3) HPS-170 high-precision linear stage 170 mm width, 52 mm travel range
50 nm minimum incremental motion
Stepper motor
Linear encoder with sin/cos signal transmission
Optical limit switches
350 N load capacity
[Figure 5]
Figure 5
A schematic illustration of the sub-micrometer motor sample stack. `DAC' indicates the position of a sample mounted in a DAC. For clarity, individual components are represented by different colors and are numbered. The corresponding coordinate system employed at the beamline is shown at the bottom of the figure. The x axis is collinear to the X-ray beam direction. The Perkin–Elmer XRD1621 detector is shown for reference.

As indicated in Fig. 5[link], the sample stack is set up on the kinematic mounting plate. It has three cavities matching the three ball point joints mounted on the base granite installed on the beamline. After the installation, the kinematic plate is fixed to the underlying granite by means of several M12 screws. The change from a standard setup (also based on the kinematic stage approach) to the sub-micrometer configuration takes a few hours, including motor reconnection and reconfiguration. Optimal operation parameters of the motors were determined during the commissioning, e.g. the XY cross-roller stage was optimized for 0.2 nm reproducibility of movement. Loading capacities given in Table 1[link] indicate that the sub-micrometer sample stack is not compatible with heavy DAC environments, e.g. a vacuum chamber, which is typical for the standard setup. However, the sub-micrometer sample stack can easily hold DACs of various popular designs (e.g. symmetric, BX90, Boehler–Almax etc.), including a membrane connected to a pressure controller. The positioning of the samples by means of the sub-micrometer motor stack is fully integrated into the beamline software. It is fully compatible with our standard setup (e.g. synchronization of sample positions by means of offline and online microscopes available for users).

Below we discuss the performance of the sub-micrometer instrument setup from a user's perspective, with several examples of single-crystal studies under ambient and high-pressure conditions.

4. Performance of the sub-micrometer experimental setup with samples under ambient conditions

Two single crystals with different absorption lengths (high and low Z) were used to test the accuracy and precision of the sample-positioning system: orthoenstatite [(Mg1.93,Fe0.06)(Si1.93,Al0.06)O6, Pbca, a = 18.2391 (3), b = 8.8117 (2), c = 5.18320 (10) Å] and CoSb3 [[Im {\overline 3}], a = 9.0357 (3) Å], measured at Bayerisches Geoinstitute using the eight-position centering procedure (Angel & Finger, 2011[Angel, R. J. & Finger, L. W. (2011). J. Appl. Cryst. 44, 247-251.]). Both materials are stable under ambient conditions and have large unit cells that produce a considerable number of reflections, which is ideal for testing.

Tiny single crystals of orthoenstantite and CoSb3 were prepared using an FIB instrument (FEI SCIOS dual beam) located at DESY NanoLab. The pre-cut single-crystalline lamellae were reduced to smaller sizes with lateral dimensions of 2–4 µm (Fig. 6[link]). These grains were transferred by means of a Microsupport Axis Pro micromanipulator into DACs used as sample holders, ensuring sample safety in transit. We employed a Re gasket of ∼10 µm thickness indented by 40 µm culet diamond anvils of Boehler–Almax design. These test crystals were measured under ambient conditions.

[Figure 6]
Figure 6
SEM images produced by the SCIOS dual-beam FIB at the NanoLab. (a) Prepared pieces of orthoenstatite. (b) CoSb3 single crystals.

Single-crystal samples were measured using a Perkin–Elmer XRD1621 amorphous silicon detector bonded to a CsI scintillator. During the data acquisition, the crystals were rotated by ±32° and ±38° with a step size of 0.5° for CoSb3 and orthoenstatite, respectively. The wavelength of the X-ray beam was tuned to 0.4830 Å (25.67 keV). Alignment of the sample to the center of the sub-micrometer motor stack rotation was performed with the help of X-ray transmission scans. We used the CoSb3 absorption signal as a reference, while in the case of orthoenstatite (separate DAC, low-Z material) we had to use the absorption signal from a thin gasket and then locate the orthoenstatite crystal using a visible light microscope on the beamline.

Data were processed using the CrysAlis PRO software package (Rigaku, 2019[Rigaku (2019). CrysAlis PRO. Version 171.40. Rigaku Corporation, Wroclaw, Poland.]), including the SCALE3 ABSPACK routine for empirical absorption correction and conventional single-crystal data treatment. The output of the program demonstrates the exceptional sample stability during the data collection of CoSb3 (Fig. 7[link]). A minor deviation was detected in the case of the orthoenstatite sample, and we attribute this to the difficulty of sample centering for this poor X-ray scattering material.

[Figure 7]
Figure 7
The output of the SCALE3 ABSPACK empirical absorption correction routine of CrysAlis PRO, indicating the amount of scaling which had to be applied to each frame in order to compensate for the intensity mismatch between different Friedel equivalents. Each frame corresponds to a physical angle with a step size of 0.5° per frame. (a) The flat red line with 2σ = 1.2% is attributed to an average value of 1.011 (gray dashed line). It indicates the excellent stability of the CoSb3 sample with respect to the X-ray beam during single-crystal data acquisition. (b) Our analysis shows that the orthoenstatite single crystal was moving slightly out of the X-ray beam in the angular range corresponding to frame Nos. 0–20. We consider that the centering procedure involving the gasket absorption profile combined with visible light observations was not perfect enough. The difference in point scatter between the red and blue lines can be attributed to the difference in scattering factors. For the data in panel (a) we had to use a 50 µm Pt absorber foil in order to reduce the intensity of the diffraction signal coming from the sample. The much stronger scattering from CoSb3 than from orthoenstatite is the reason for the noise difference between the red and blue curves.

Structure solution and refinement were conducted using OLEX2 (Dolomanov et al., 2009[Dolomanov, O. V., Bourhis, L. J., Gildea, R. J., Howard, J. A. K. & Puschmann, H. (2009). J. Appl. Cryst. 42, 339-341.]) with the SHELX backend (Sheldrick, 2008[Sheldrick, G. M. (2008). Acta Cryst. A64, 112-122.], 2015[Sheldrick, G. M. (2015). Acta Cryst. A71, 3-8.]) and the JANA2006 package (Petříček et al., 2006[Petříček, V., Dušek, M. & Palatinus, L. (2006). JANA2006. Institute of Physics, Czech Academy of Sciences, Prague, Czech Republic.]). The refinement parameters and the corresponding structural information are summarized in Table 2[link]. Additional information, including CIFs, is provided in the supporting information. Illustrations of the crystal structure in Fig. 8[link] were prepared with the CrystalMaker software (Palmer, 2014[Palmer, D. C. (2014). CrystalMaker. CrystalMaker Software Ltd, Begbroke, England.]).

Table 2
Details of the crystal structure refinements of the test samples orthoenstatite (Mg1.93,Fe0.06)(Si1.93,Al0.06)O6 and CoSb3

Atomic displacement parameters of the cations were refined with anisotropic approximation. For additional information on the crystal structures and refinement parameters we refer the reader to the CIFs in the supporting information.

Crystallographic data
Chemical formula (Mg1.94Fe0.067)(Si1.93Al0.067)O6 CoSb3
Mr 202.9 424.2
Crystal system, space group Orthorhombic, Pbca Cubic, Im3
Temperature (K) 293 293
a, b, c (Å) 5.1815 (2), 18.2321 (11), 8.8085 (5) 9.0360 (1)
V3) 832.14 (8) 737.78 (1)
Z 8 8
Radiation type Synchrotron, λ = 0.483 Å Synchrotron, λ = 0.483 Å
μ (mm−1) 0.44 8.83
Crystal size (µm3) 4 × 4 × 2 3 × 3 × 2
     
Data collection
Diffractometer Single-circle (ω) diffractometer
Absorption correction Multi-scan, CrysAlis PRO 1.171.40.67a (Rigaku, 2019[Rigaku (2019). CrysAlis PRO. Version 171.40. Rigaku Corporation, Wroclaw, Poland.]), empirical absorption correction using spherical harmonics, implemented in the SCALE3 ABSPACK scaling algorithm
Tmin, Tmax 0.824, 1 0.888, 1
No. of measured, independent and observed [I > 2σ(I)] reflections 1072, 454, 409 490, 133, 121
Rint 0.017 0.016
(sin θ/λ)max−1) 0.652 0.645
     
Refinement
R[F2 > 2σ(F2)], wR(F2) 0.025, 0.039 0.012, 0.043
No. of reflections 454 133
No. of parameters 61 9
Δρmax, Δρmin (e Å−3) 0.58, −0.62 1.28, −1.03
[Figure 8]
Figure 8
The crystal structure of Pnma Fe3O4 formed after laser heating at 3000–3500 K. The structure is composed of layers of distorted octahedra (Fe1 and Fe3 sites) which are interconnected through the capped trigonal prisms (Fe2 sites). The structure solution converged to R(F2) = 0.071, wR(F2) = 0.075 using reflections with an intensity I > 2σ(I). The individual building blocks and the corresponding Fe—O bonds are shown. On the right-hand side is a projection of the structure along the b axis, providing a clear indication of edge sharing of the corresponding building blocks and an overview of the Fe–O framework.

The data quality reported in Table 2[link] is exceptional and in very good agreement with the literature data, even though the energy of the X-ray beam was very low (25.67 keV) which significantly limits our access to reciprocal space due to the finite aperture of DACs. Future developments and enhancements expected for diffraction-limited storage rings (DLSRs), such as the ESRF EBS, APS-U and PETRA IV, will boost these types of studies. Considering the ECB at PETRA III, the upgrade to PETRA IV (Schroer et al., 2019[Schroer, C., Agapov, I., Roehlsberger, R., Wanzenberg, R., Brinkmann, R., Weckert, E. & Leemans, W. (2019). PETRA IV: Upgrade of PETRA III to the Ultimate 3D X-ray Microscope, Conceptual Design Report (CDR). https://bib-pubdb1.desy.de/record/426140.]) promises improved capabilities. The upgrade will result in photon flux improvement at all photon energies interesting for high-pressure research and, in combination with sub-micrometer focusing, it will exceed the capabilities of third-generation sources.

5. Applications to materials compressed above ∼150 GPa in laser-heated DACs

A successful crystallographic solution for the test cases described below would not be possible with our standard motor setup. The implementation of the new motorization enabled more accurate and reproducible sample positioning in 3D space, including sample positioning at the center of rotation of the highly accurate rotation stage. Due to the improvements in the motorization and the implementation of the correcting phase plate, which reduces the overload of the sample signal in reciprocal space, the resulting data quality was greatly improved. Each of the described cases represents an important challenge for various fields of science, such as crystal chemistry, mineral physics and planetary sciences.

Iron oxides are closely related to planetary evolution (Bykova et al., 2016[Bykova, E., Dubrovinsky, L., Dubrovinskaia, N., Bykov, M., McCammon, C., Ovsyannikov, S. V., Liermann, H., Kupenko, I., Chumakov, A. I., Rüffer, R., Hanfland, M. & Prakapenka, V. (2016). Nat. Commun. 7, 10661.]; Dobson & Brodholt, 2005[Dobson, D. P. & Brodholt, J. P. (2005). Nature, 434, 371-374.]), but also play a significant role in thermal energy storage, as indicated by various studies (Grosu et al., 2017[Grosu, Y., Faik, A., Ortega-Fernández, I. & D'Aguanno, B. (2017). Solar Energy Mater. Solar Cells, 161, 170-176.]; Huang & Xu, 2019[Huang, H. & Xu, J. Y. (2019). Environ. Sci. 267, 042019.]). High-pressure high-temperature treatment of the Fe–O system at elevated pressures induces various phases and structural transformations. Below we report on the observation of the Pnma phase of Fe3O4 at ∼200 GPa. The starting material, α-Fe2O3 (hematite), was loaded into the sample chamber of a BX90 DAC (Kantor et al., 2012[Kantor, I., Prakapenka, V., Kantor, A., Dera, P., Kurnosov, A., Sinogeikin, S., Dubrovinskaia, N. & Dubrovinsky, L. (2012). Rev. Sci. Instrum. 83, 125102.]) together with Ne as the pressure-transmitting medium (Kurnosov et al., 2008[Kurnosov, A., Kantor, I., Boffa-Ballaran, T., Lindhardt, S., Dubrovinsky, L., Kuznetsov, A. & Zehnder, B. H. (2008). Rev. Sci. Instrum. 79, 045110.]). The starting material was laser-heated at the Bayerisches Geo­institut using a laser-heating system recently developed by Fedotenko et al. (2019[Fedotenko, T., Dubrovinsky, L., Aprilis, G., Koemets, E., Snigirev, A., Snigireva, I., Barannikov, A., Ershov, P., Cova, F., Hanfland, M. & Dubrovinskaia, N. (2019). Rev. Sci. Instrum. 90, 104501.]). This is a pulsed laser-heating system operating in continuous mode (NIR laser beam, 1070 nm, 5 µm at FWHM in diameter). Pressure was measured using the equation of state of Ne (Fei et al., 2007[Fei, Y., Ricolleau, A., Frank, M., Mibe, K., Shen, G. & Prakapenka, V. (2007). Proc. Natl Acad. Sci. USA, 104, 9182-9186.]). Single-crystal datasets were collected using the Perkin–Elmer XRD1621 detector with a wavelength of 0.483 Å. A complete description of the structure of Pnma Fe3O4 observed at 200 GPa, including a CIF, is provided in the supporting information. The small X-ray beam size was crucial as it constrained the number of illuminated grains and the overlay of scattered intensity in reciprocal space originating from various sources in the high-pressure environment.

The laser heating of the compressed hematite produced a polycrystalline mixture. Analysis of 2D scanning map data of the laser-heated area indicated several larger single-crystal grains that are potential candidates for structural solution. The indexing of the reflections belonging to two strongly scattering Pnma Fe3O4 (oP-Fe3O4) domains was identified after careful inspection of reciprocal space using the CrysAlis PRO software. A more detailed description of the process is provided in the supporting information. The lateral size of the domains within the 2D map did not exceed a couple of micrometers. The integration of intensities for the most strongly scattering domain converged with Rint = 0.029. The structure was solved and refined based on 195 measured reflections for 22 refined parameters with R1 = 0.071 (Table 3[link]). The structure of oP-Fe3O4 contains two different octahedrally coordinated positions of iron atoms (Fe1 and Fe3) with the Fe—O bond distance varying between ∼1.70 and 1.80 Å (Fig. 8[link]). The distorted octahedra share edges and form ribbons, and two enantiomorphic ribbons are joined together, hinging along lines of the shared oxygen parallel to the c axis. The Fe2 atom is surrounded by seven oxygen atoms with an Fe—O distance variation of ∼1.86–1.91 Å. Its coordination can be considered as a capped trigonal prism. The structure of novel oP-Fe3O4 belongs to the Yb3S4 structure type (Chevallier et al., 1967[Chevallier, R., Laruelle, P. & Flahaut, J. (1967). Bull. Soc. Fr. Mineral. Cristallogr. 90, 564-574.]). It is common for warwickites, mixed metal borates of the general formula AxT4–xB2O8, where A and T are alkaline earths, lanthanides or transition metals (Matos et al., 1996[Matos, M., Hoffmann, R., Latgé, A. & Anda, E. V. (1996). Chem. Mater. 8, 2324-2330.]).

Table 3
Detailed information for the crystal structure refinements of oP-Fe3O4 and Mg0.91 (2)Fe0.09 (2)SiO3 performed on laser-heated samples at ultra-high pressures

Samples were measured at ambient temperature after a laser-heating procedure. For additional information on the crystal structures and refinement parameters we refer the reader to the CIFs in the supporting information.

Crystal data
Chemical formula Fe3O4 Mg0.91 (2)Fe0.09 (2)SiO3
Pressure (GPa) ∼200 ∼155
Mr 231.5 103.3
Crystal system, space group Orthorhombic, Pnma Orthorhombic, Pbnm
Temperature (K) 293 293
a, b, c (Å) 7.932 (10), 2.5881 (13), 8.321 (3) 4.194 (2), 4.525 (1), 6.1910 (1)
V3) 170.8 (2) 117.49 (6)
Z 4 4
Radiation type Synchrotron, λ = 0.483 Å Synchrotron, λ = 0.483 Å
μ (mm−1) 7.97 0.98
Crystal size (µm3) 1 × 2 × 1 2 × 2 × 2
     
Data collection
Diffractometer Single-circle (ω) diffractometer
Absorption correction Multi-scan, CrysAlis PRO 1.171.40.67a (Rigaku, 2019[Rigaku (2019). CrysAlis PRO. Version 171.40. Rigaku Corporation, Wroclaw, Poland.]), empirical absorption correction using spherical harmonics, implemented in the SCALE3 ABSPACK scaling algorithm
Tmin, Tmax 0.559, 1 0.664, 1
No. of measured, independent and observed [I > 2σ(I)] reflections 195, 105, 88 106, 56, 53
Rint 0.029 0.008
(sin θ/λ)max−1) 0.643 0.607
     
Refinement
R[F2 > 2σ(F2)], wR(F2) 0.071, 0.075 0.041, 0.114
No. of reflections 105 56
No. of parameters 22 13
Δρmax, Δρmin (e Å−3) 2.04, −1.8 0.51, −0.54
†Approximate sizes based on X-ray diffraction 2D mapping.

Charge-balance considerations suggest that the Fe1 and Fe3 sites are occupied by Fe3+. Our analysis of the Fe1—O and Fe3—O distances, supplemented by analysis of the octahedron site volumes, reveals that the Fe1 and Fe3 atoms are in a low-spin state (Vasiukov, 2018[Vasiukov, D. (2018). PhD thesis. University of Bayreuth, Germany.]). The Fe2 site is attributed to Fe2+. It is clearly larger than the Fe1 and Fe3 sites. We suggest that the spin configuration of Fe2 should be also attributed to low spin. However, due to the deformation of the crystallographic site we cannot fully exclude the possibility of an intermediate spin-state configuration. The precise characterization of the iron atoms' electronic state in oP-Fe3O4 and the exploration of its phase diagram require further investigation.

The Mg–Si–O system is equally if not more important than Fe–O for our understanding of the Earth and other planetary environments. In addition to the example of the Fe3O4 system, we show a full structure solution and refinement of excellent quality for Fe-bearing bridgmanite Mg0.91 (2)Fe0.09 (2)SiO3. The data were obtained after laser heating and collected at a pressure of ∼155 GPa (Ne pressure scale; Fei et al., 2007[Fei, Y., Ricolleau, A., Frank, M., Mibe, K., Shen, G. & Prakapenka, V. (2007). Proc. Natl Acad. Sci. USA, 104, 9182-9186.]). Accurate measurements of the physical properties of bridgmanite and post-perovskite, including density as a function of pressure, temperature and composition, have always been a focus of solid-earth investigations. The studies continue with a doubled effort today involving static and dynamic compression methods [e.g. Fei et al. (2021[Fei, Y., Seagle, C. T., Townsend, J. P., McCoy, C. A., Boujibar, A., Driscoll, P., Shulenburger, L. & Furnish, M. D. (2021). Nat. Commun. 12, 876.])]. However, information characterizing the density alone can rightfully be considered insufficient. It must be supplemented with a thorough understanding of the crystal structure and its evolution and perturbation. Here, ultra-high-pressure studies of the Mg–Si–O system by means of SCXRD represent a very specific but very important challenge. Most experimental difficulties could be tracked down to the low scattering volumes of the sub-micrometer-sized samples, but they become even more severe due to the low atomic scattering factor in the case of the Mg–Si–O system. Using the example of bridgmanite laser-heated to 2500–3000 K, we have demonstrated the power of the new sub-micrometer setup on P02.2, enabling structural studies of chemically complex systems at ultra-high pressures, even those containing mostly light elements in their composition (e.g. Mg, Si, O). As we show in Table 3[link], the final refinement of the structural model converges with a very low R factor of R1 = 0.0411 based on 106 measured reflections and 13 refined parameters.

6. Conclusions

By applying the most recent technical developments within the capabilities of PETRA III, we report an upgrade of the existing high-pressure micro-X-ray diffraction instrument at P02.2. The upgrade involves the modification of the X-ray optical path and the implementation of a phase plate correcting the spherical aberration of the CRLs implemented at P02.2 in combination with the new, more accurate and precise, motorization of the sample stack.

Within our test studies conducted under ambient conditions and at pressures exceeding 1 Mbar, we have demonstrated the most challenging issues faced in high-pressure science and presented a solution addressing the challenges of single-crystal refinement from grains of micrometer to sub-micrometer dimensions. We have described `synthetic' examples of orthoenstatite and CoSb3 measured inside the DAC under ambient conditions and extended our studies to scientific cases exploring the crystal structures of Fe3O4 and Mg0.91 (2)Fe0.09 (2)SiO3. The structures of the latter two materials were successfully solved at pressures exceeding 1 Mbar, even though the data acquisition and analysis processes were highly challenging, e.g. complicated by sample finding and by signal extraction from micrometer-sized polycrystalline aggregates. The real case studies of Fe3O4 and Mg0.91 (2)Fe0.09 (2)SiO3 demonstrate the capacity of the current instrument upgrade.

We hope that our study will inspire the high-pressure community and engender further developments on dedicated high-pressure beamlines, extending the capabilities of the latter and driving high-pressure science beyond the boundaries of the commonly explored PT space.

7. Related literature

For further literature related to the supporting information, see Bykova (2015[Bykova, E. (2015). PhD thesis. University of Bayreuth, Germany. https://epub.uni-bayreuth.de/2124/]), Hrubiak et al. (2019[Hrubiak, R., Smith, J. S. & Shen, G. (2019). Rev. Sci. Instrum. 90, 025109.]), Petříček et al. (2014[Petříček, V., Dušek, M. & Palatinus, L. (2014). Z. Kristallogr. Cryst. Mater. 229, 345-352.]), Prescher & Prakapenka (2015[Prescher, C. & Prakapenka, V. B. (2015). High Press. Res. 35, 223-230.]) and Rigaku Oxford Diffraction (2020[Rigaku Oxford Diffraction (2020). CrysAlis PRO. Rigaku Oxford Diffraction, The Woodlands, Texas, USA.]).

Supporting information


Computing details top

Data collection: CrysAlis PRO 1.171.40.67a (Rigaku OD, 2019); cell refinement: CrysAlis PRO 1.171.40.67a (Rigaku OD, 2019); data reduction: CrysAlis PRO 1.171.40.67a (Rigaku OD, 2019); program(s) used to solve structure: Jana2006 Version : 24/08/2020; program(s) used to refine structure: Jana2006 Version : 24/08/2020; molecular graphics: Vesta, Version 3.4.0 : 08/01/2017.

Magnesium iron aluminium catena-alumosilicate top
Crystal data top
Mg1.94Fe0.067Si1.93Al0.067O6F(000) = 807
Mr = 202.9Dx = 3.238 Mg m3
Orthorhombic, PbcaSynchrotron radiation, λ = 0.483 Å
Hall symbol: -P -2xab;-2ybc;-2zacCell parameters from 547 reflections
a = 5.1815 (2) Åθ = 2.2–18.4°
b = 18.2321 (11) ŵ = 0.44 mm1
c = 8.8085 (5) ÅT = 293 K
V = 832.14 (8) Å3Plate, grey
Z = 80.004 × 0.004 × 0.002 mm
Data collection top
Customized ω-circle
diffractometer
454 independent reflections
Radiation source: synchrotron409 reflections with I > 2σ(I)
Detector resolution: 5.0 pixels mm-1Rint = 0.017
ω scansθmax = 18.4°, θmin = 2.2°
Absorption correction: multi-scan
CrysAlisPro 1.171.40.84a (Rigaku Oxford Diffraction, 2020) Empirical absorption correction using spherical harmonics, implemented in SCALE3 ABSPACK scaling algorithm.
h = 66
Tmin = 0.824, Tmax = 1k = 1819
1072 measured reflectionsl = 89
Refinement top
Refinement on F0 restraints
R[F2 > 2σ(F2)] = 0.02518 constraints
wR(F2) = 0.039Weighting scheme based on measured s.u.'s w = 1/(σ2(F) + 0.0001F2)
S = 1.87(Δ/σ)max = 0.045
454 reflectionsΔρmax = 0.58 e Å3
61 parametersΔρmin = 0.62 e Å3
Special details top

Experimental. Dataset measured while the sample was loaded inside a Diamond anvil cell (DAC)

Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/UeqOcc. (<1)
Si10.04913 (17)0.22856 (6)0.15836 (13)0.0037 (5)
Si20.29930 (18)0.02637 (6)0.33723 (12)0.0037 (5)0.933
Al20.29930 (18)0.02637 (6)0.33723 (12)0.0037 (5)0.067
Mg30.3580 (2)0.12261 (7)0.01415 (15)0.0095 (6)0.933
Fe30.3580 (2)0.12261 (7)0.01415 (15)0.0095 (6)0.067
Mg40.3653 (2)0.37586 (7)0.15360 (14)0.0029 (6)
O10.1055 (4)0.05238 (15)0.1946 (3)0.0058 (6)*
O20.0359 (4)0.31637 (16)0.1600 (3)0.0048 (7)*
O30.1989 (4)0.43755 (16)0.3397 (3)0.0046 (7)*
O40.1898 (4)0.06751 (16)0.4831 (3)0.0067 (7)*
O50.3291 (4)0.19721 (15)0.2236 (3)0.0059 (7)*
O60.0432 (4)0.31113 (16)0.4976 (3)0.0058 (7)*
Atomic displacement parameters (Å2) top
U11U22U33U12U13U23
Si10.0038 (5)0.0063 (10)0.0011 (10)0.0000 (5)0.0005 (4)0.0003 (3)
Si20.0038 (5)0.0055 (10)0.0020 (10)0.0008 (5)0.0003 (4)0.0001 (4)
Al20.0038 (5)0.0055 (10)0.0020 (10)0.0008 (5)0.0003 (4)0.0001 (4)
Mg30.0085 (6)0.0116 (11)0.0084 (11)0.0015 (6)0.0007 (5)0.0009 (4)
Fe30.0085 (6)0.0116 (11)0.0084 (11)0.0015 (6)0.0007 (5)0.0009 (4)
Mg40.0027 (6)0.0033 (11)0.0028 (12)0.0001 (6)0.0002 (5)0.0000 (4)
Geometric parameters (Å, º) top
Si1—Mg32.8119 (16)Mg3—Mg4ii3.1762 (19)
Si1—O21.603 (3)Mg3—O2vi2.107 (3)
Si1—O51.662 (2)Mg3—O3ii2.060 (3)
Si1—O5i1.646 (3)Mg3—O4iv1.991 (3)
Si1—O6ii1.590 (3)Mg3—O6ii2.035 (3)
Si2—Al20Fe3—Mg4v2.9501 (16)
Si2—O3iii1.619 (3)Fe3—O2vi2.107 (3)
Si2—O41.592 (3)Fe3—O3ii2.060 (3)
Al2—Mg3i3.1655 (16)Fe3—O4iv1.991 (3)
Al2—O11.677 (3)Fe3—O6ii2.035 (3)
Al2—O1iv1.680 (2)Mg4—Mg4i3.0977 (16)
Al2—O3iii1.619 (3)Mg4—Mg4iv3.0977 (16)
Al2—O41.592 (3)Mg4—O22.023 (3)
Mg3—Fe30Mg4—O3iv2.063 (3)
Mg3—Mg4v2.9501 (16)Mg4—O4ii2.037 (3)
Mg3—Mg4vi3.0153 (16)Mg4—O6iv2.004 (3)
Mg3—Si1—O2135.62 (10)O2vi—Fe3—O6ii89.14 (11)
Mg3—Si1—O554.75 (10)O3ii—Fe3—O4iv94.92 (12)
Mg3—Si1—O5i116.18 (11)O3ii—Fe3—O6ii86.67 (11)
Mg3—Si1—O6ii45.27 (10)O4iv—Fe3—O6ii172.91 (13)
O2—Si1—O5112.21 (13)Mg3v—Mg4—Mg3vi120.59 (5)
O2—Si1—O5i108.16 (14)Mg3v—Mg4—Mg3vii118.66 (5)
O2—Si1—O6ii117.49 (15)Mg3v—Mg4—Fe3vi120.59 (5)
O5—Si1—O5i105.47 (14)Mg3v—Mg4—Mg4i62.59 (4)
O5—Si1—O6ii99.67 (14)Mg3v—Mg4—Mg4iv176.07 (6)
O5i—Si1—O6ii113.04 (14)Mg3v—Mg4—O244.21 (8)
Al2—Si2—O3iii0Mg3v—Mg4—O3iv137.67 (9)
Al2—Si2—O40Mg3v—Mg4—O4ii40.97 (8)
O3iii—Si2—O4117.52 (15)Mg3v—Mg4—O6iv94.63 (9)
Si2—Al2—Mg3i0Mg3vi—Mg4—Mg3vii120.73 (5)
Si2—Al2—O10Mg3vi—Mg4—Fe3vi0
Si2—Al2—O1iv0Mg3vi—Mg4—Mg4i176.77 (6)
Si2—Al2—O3iii0Mg3vi—Mg4—Mg4iv63.31 (4)
Si2—Al2—O40Mg3vi—Mg4—O2138.51 (10)
Mg3i—Al2—O173.78 (10)Mg3vi—Mg4—O3iv44.27 (8)
Mg3i—Al2—O1iv126.49 (11)Mg3vi—Mg4—O4ii90.70 (8)
Mg3i—Al2—O3iii123.56 (10)Mg3vi—Mg4—O6iv43.49 (8)
Mg3i—Al2—O431.56 (9)Mg3vii—Mg4—Fe3vi120.73 (5)
O1—Al2—O1iv111.15 (14)Mg3vii—Mg4—Mg4i56.08 (4)
O1—Al2—O3iii107.24 (14)Mg3vii—Mg4—Mg4iv57.43 (4)
O1—Al2—O4104.96 (14)Mg3vii—Mg4—O288.11 (9)
O1iv—Al2—O3iii106.22 (13)Mg3vii—Mg4—O3iv88.66 (8)
O1iv—Al2—O4109.76 (14)Mg3vii—Mg4—O4ii136.68 (9)
O3iii—Al2—O4117.52 (15)Mg3vii—Mg4—O6iv132.46 (10)
Si1—Mg3—Al2iv127.24 (5)Fe3vi—Mg4—Mg4i176.77 (6)
Si1—Mg3—Fe30Fe3vi—Mg4—Mg4iv63.31 (4)
Si1—Mg3—Mg4v74.17 (4)Fe3vi—Mg4—O2138.51 (10)
Si1—Mg3—Mg4vi135.33 (6)Fe3vi—Mg4—O3iv44.27 (8)
Si1—Mg3—Mg4ii116.90 (5)Fe3vi—Mg4—O4ii90.70 (8)
Si1—Mg3—O2vi102.42 (9)Fe3vi—Mg4—O6iv43.49 (8)
Si1—Mg3—O3ii118.36 (8)Mg4i—Mg4—Mg4iv113.51 (5)
Si1—Mg3—O4iv146.36 (10)Mg4i—Mg4—O243.88 (8)
Si1—Mg3—O6ii33.73 (8)Mg4i—Mg4—O3iv133.25 (9)
Al2iv—Mg3—Fe30Mg4i—Mg4—O4ii91.77 (8)
Al2iv—Mg3—Mg4v146.87 (6)Mg4i—Mg4—O6iv138.51 (9)
Al2iv—Mg3—Mg4vi65.05 (4)Mg4iv—Mg4—O2133.66 (10)
Al2iv—Mg3—Mg4ii114.19 (5)Mg4iv—Mg4—O3iv44.24 (8)
Al2iv—Mg3—O2vi106.13 (8)Mg4iv—Mg4—O4ii141.04 (9)
Al2iv—Mg3—O3ii107.60 (9)Mg4iv—Mg4—O6iv88.83 (8)
Al2iv—Mg3—O4iv24.74 (8)O2—Mg4—O3iv176.72 (12)
Al2iv—Mg3—O6ii159.66 (10)O2—Mg4—O4ii85.16 (11)
Fe3—Mg3—Mg4v0O2—Mg4—O6iv95.21 (12)
Fe3—Mg3—Mg4vi0O3iv—Mg4—O4ii96.82 (11)
Fe3—Mg3—Mg4ii0O3iv—Mg4—O6iv87.39 (11)
Fe3—Mg3—O2vi0O4ii—Mg4—O6iv90.80 (12)
Fe3—Mg3—O3ii0Al2—O1—Al2i127.64 (17)
Fe3—Mg3—O4iv0Si1—O2—Mg3v122.66 (15)
Fe3—Mg3—O6ii0Si1—O2—Fe3v122.66 (15)
Mg4v—Mg3—Mg4vi120.59 (6)Si1—O2—Mg4119.95 (14)
Mg4v—Mg3—Mg4ii60.61 (4)Mg3v—O2—Fe3v0
Mg4v—Mg3—O2vi90.51 (8)Mg3v—O2—Mg493.77 (12)
Mg4v—Mg3—O3ii44.36 (8)Fe3v—O2—Mg493.77 (12)
Mg4v—Mg3—O4iv139.29 (10)Si2viii—O3—Al2viii0
Mg4v—Mg3—O6ii42.68 (8)Si2viii—O3—Mg3vii122.70 (15)
Mg4vi—Mg3—Mg4ii59.97 (4)Si2viii—O3—Fe3vii122.70 (15)
Mg4vi—Mg3—O2vi42.03 (7)Si2viii—O3—Mg4i123.38 (14)
Mg4vi—Mg3—O3ii89.32 (8)Al2viii—O3—Mg3vii122.70 (15)
Mg4vi—Mg3—O4iv42.10 (8)Al2viii—O3—Fe3vii122.70 (15)
Mg4vi—Mg3—O6ii131.15 (10)Al2viii—O3—Mg4i123.38 (14)
Mg4ii—Mg3—O2vi42.46 (8)Mg3vii—O3—Fe3vii0
Mg4ii—Mg3—O3ii42.57 (8)Mg3vii—O3—Mg4i91.36 (12)
Mg4ii—Mg3—O4iv90.36 (9)Fe3vii—O3—Mg4i91.36 (12)
Mg4ii—Mg3—O6ii86.14 (9)Si2—O4—Al20
O2vi—Mg3—O3ii85.02 (11)Si2—O4—Mg3i123.70 (15)
O2vi—Mg3—O4iv84.12 (11)Si2—O4—Fe3i123.70 (15)
O2vi—Mg3—O6ii89.14 (11)Si2—O4—Mg4vii132.45 (14)
O3ii—Mg3—O4iv94.92 (12)Al2—O4—Mg3i123.70 (15)
O3ii—Mg3—O6ii86.67 (11)Al2—O4—Fe3i123.70 (15)
O4iv—Mg3—O6ii172.91 (13)Al2—O4—Mg4vii132.45 (14)
Mg3—Fe3—Mg4v0Mg3i—O4—Fe3i0
Mg3—Fe3—O2vi0Mg3i—O4—Mg4vii96.93 (13)
Mg3—Fe3—O3ii0Fe3i—O4—Mg4vii96.93 (13)
Mg3—Fe3—O4iv0Si1—O5—Si1iv134.73 (18)
Mg3—Fe3—O6ii0Si1vii—O6—Mg3vii101.01 (14)
Mg4v—Fe3—O2vi90.51 (8)Si1vii—O6—Fe3vii101.01 (14)
Mg4v—Fe3—O3ii44.36 (8)Si1vii—O6—Mg4i150.27 (16)
Mg4v—Fe3—O4iv139.29 (10)Mg3vii—O6—Fe3vii0
Mg4v—Fe3—O6ii42.68 (8)Mg3vii—O6—Mg4i93.83 (13)
O2vi—Fe3—O3ii85.02 (11)Fe3vii—O6—Mg4i93.83 (13)
O2vi—Fe3—O4iv84.12 (11)
Symmetry codes: (i) x1/2, y, z+1/2; (ii) x, y+1/2, z1/2; (iii) x+1/2, y1/2, z; (iv) x+1/2, y, z+1/2; (v) x1/2, y+1/2, z; (vi) x+1/2, y+1/2, z; (vii) x, y+1/2, z+1/2; (viii) x+1/2, y+1/2, z.
 

Footnotes

Current affiliation: XRnanotech GmbH, Forschungsstrasse 111, 5232 Villigen-PSI, Switzerland.

Acknowledgements

We acknowledge DESY (Hamburg, Germany), a member of the Helmholtz Association HGF, for the provision of experimental facilities. Parts of this research were carried out at beamlines P02.2 and P06 of PETRA III and we would additionally like to thank Andreas Stierle, Satishkumar Kulkarni and Thomas Keller for assistance in using the NanoLab facilities. We are grateful to XRnanotech GmbH, current employer of Adam Kubec, which can manufacture the phase plates presented in this work. We acknowledge Mirtha Pillaca and Peter Gille (Ludwig-Maximilians-Universität, München, Germany) for their great help with the preparation of the CoSb3 samples. We are grateful to Tiziana Boffa-Ballaran (Bayerisches Geoinstitut, Bayreuth, Germany) for the eight-circle position measurements of orthoenstatite and CoSb3.

Funding information

ND and LD thank the Federal Ministry of Education and Research, Germany (BMBF, grant No. 05K19WC1) and the Deutsche Forschungsgemeinschaft (DFG projects DU 954-11/1, DU 393-9/2 and DU 393-13/1) for financial support. ND thanks the Swedish Government Strategic Research Area in Materials Science on Functional Materials at Linköping University (Faculty Grant SFO-Mat-LiU No. 2009 00971). DL thanks the Alexander von Humboldt Foundation for financial support.

References

First citationAngel, R. J. & Finger, L. W. (2011). J. Appl. Cryst. 44, 247–251.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationAnzellini, S., Monteseguro, V., Bandiello, E., Dewaele, A., Burakovsky, L. & Errandonea, D. (2019). Sci. Rep. 9, 13034.  Web of Science CrossRef PubMed Google Scholar
First citationBykova, E. (2015). PhD thesis. University of Bayreuth, Germany. https://epub.uni-bayreuth.de/2124/  Google Scholar
First citationBykova, E., Dubrovinsky, L., Dubrovinskaia, N., Bykov, M., McCammon, C., Ovsyannikov, S. V., Liermann, H., Kupenko, I., Chumakov, A. I., Rüffer, R., Hanfland, M. & Prakapenka, V. (2016). Nat. Commun. 7, 10661.  Web of Science CrossRef ICSD PubMed Google Scholar
First citationBykov, M., Chariton, S., Fei, H., Fedotenko, T., Aprilis, G., Ponomareva, A. V., Tasnádi, F., Abrikosov, I. A., Merle, B., Feldner, P., Vogel, S., Schnick, W., Prakapenka, V. B., Greenberg, E., Hanfland, M., Pakhomova, A., Liermann, H., Katsura, T., Dubrovinskaia, N. & Dubrovinsky, L. (2019). Nat. Commun. 10, 2994.  Web of Science CrossRef ICSD PubMed Google Scholar
First citationCelestre, R., Berujon, S., Roth, T., Sanchez del Rio, M. & Barrett, R. (2020). J. Synchrotron Rad. 27, 305–318.  Web of Science CrossRef IUCr Journals Google Scholar
First citationChevallier, R., Laruelle, P. & Flahaut, J. (1967). Bull. Soc. Fr. Mineral. Cristallogr. 90, 564–574.  Google Scholar
First citationDewaele, A., Loubeyre, P., Occelli, F., Marie, O. & Mezouar, M. (2018). Nat. Commun. 9, 2913.  Web of Science CrossRef PubMed Google Scholar
First citationDobson, D. P. & Brodholt, J. P. (2005). Nature, 434, 371–374.  Web of Science CrossRef PubMed CAS Google Scholar
First citationDolomanov, O. V., Bourhis, L. J., Gildea, R. J., Howard, J. A. K. & Puschmann, H. (2009). J. Appl. Cryst. 42, 339–341.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationDorfman, S. M., Prakapenka, V. B., Meng, Y. & Duffy, T. S. (2012). J. Geophys. Res. 117, B08210.  Google Scholar
First citationDubrovinskaia, N., Dubrovinsky, L., Solopova, N. A., Abakumov, A., Turner, S., Hanfland, M., Bykova, E., Bykov, M., Prescher, C., Prakapenka, V. B., Petitgirard, S., Chuvashova, I., Gasharova, B., Mathis, Y., Ershov, P., Snigireva, I. & Snigirev, A. (2016). Sci. Adv. 2, e1600341.  Web of Science CrossRef PubMed Google Scholar
First citationDubrovinsky, L., Dubrovinskaia, N., Bykova, E., Bykov, M., Prakapenka, V., Prescher, C., Glazyrin, K., Liermann, H.-P., Hanfland, M., Ekholm, M., Feng, Q., Pourovskii, L. V., Katsnelson, M. I., Wills, J. M. & Abrikosov, I. A. (2015). Nature, 525, 226–229.  Web of Science CrossRef CAS PubMed Google Scholar
First citationDubrovinsky, L., Dubrovinskaia, N., Prakapenka, V. B. & Abakumov, A. M. (2012). Nat. Commun. 3, 1163.  Web of Science CrossRef PubMed Google Scholar
First citationErrandonea, D., Burakovsky, L., Preston, D. L., MacLeod, S. G., Santamaría-Perez, D., Chen, S., Cynn, H., Simak, S. I., McMahon, M. I., Proctor, J. E. & Mezouar, M. (2020). Commun. Mater. 1, 60.  Web of Science CrossRef Google Scholar
First citationFalkenberg, G., Seiboth, F., Koch, F., Falch, K. V., Schropp, A., Brückner, D. & Garrevoet, J. (2020). Powder Diffr. 35, S34–S37.  Web of Science CrossRef CAS Google Scholar
First citationFedotenko, T., Dubrovinsky, L., Aprilis, G., Koemets, E., Snigirev, A., Snigireva, I., Barannikov, A., Ershov, P., Cova, F., Hanfland, M. & Dubrovinskaia, N. (2019). Rev. Sci. Instrum. 90, 104501.  Web of Science CrossRef Google Scholar
First citationFei, Y., Ricolleau, A., Frank, M., Mibe, K., Shen, G. & Prakapenka, V. (2007). Proc. Natl Acad. Sci. USA, 104, 9182–9186.  Web of Science CrossRef PubMed CAS Google Scholar
First citationFei, Y., Seagle, C. T., Townsend, J. P., McCoy, C. A., Boujibar, A., Driscoll, P., Shulenburger, L. & Furnish, M. D. (2021). Nat. Commun. 12, 876.  Web of Science CrossRef PubMed Google Scholar
First citationFriedrich, A., Collings, I. E., Dziubek, K. F., Fanetti, S., Radacki, K., Ruiz-Fuertes, J., Pellicer-Porres, J., Hanfland, M., Sieh, D., Bini, R., Clark, S. J. & Marder, T. B. (2020). J. Am. Chem. Soc. 142, 18907–18923.  Web of Science CSD CrossRef CAS PubMed Google Scholar
First citationGrosu, Y., Faik, A., Ortega-Fernández, I. & D'Aguanno, B. (2017). Solar Energy Mater. Solar Cells, 161, 170–176.  Web of Science CrossRef CAS Google Scholar
First citationHrubiak, R., Smith, J. S. & Shen, G. (2019). Rev. Sci. Instrum. 90, 025109.  Web of Science CrossRef PubMed Google Scholar
First citationHuang, H. & Xu, J. Y. (2019). Environ. Sci. 267, 042019.  Google Scholar
First citationJenei, Z., O'Bannon, E. F., Weir, S. T., Cynn, H., Lipp, M. J. & Evans, W. J. (2018). Nat. Commun. 9, 3563.  Web of Science CrossRef PubMed Google Scholar
First citationKantor, I., Prakapenka, V., Kantor, A., Dera, P., Kurnosov, A., Sinogeikin, S., Dubrovinskaia, N. & Dubrovinsky, L. (2012). Rev. Sci. Instrum. 83, 125102.  Web of Science CrossRef PubMed Google Scholar
First citationKawamura, H., Akahama, Y., Umemoto, S., Takemura, K., Ohishi, Y. & Shimomura, O. (2002). J. Phys. Condens. Matter, 14, 10407–10410.  Web of Science CrossRef CAS Google Scholar
First citationKhandarkhaeva, S., Fedotenko, T., Bykov, M., Bykova, E., Chariton, S., Sedmak, P., Glazyrin, K., Prakapenka, V., Dubrovinskaia, N. & Dubrovinsky, L. (2020). Eur. J. Inorg. Chem. 2020, 2186–2190.  Web of Science CSD CrossRef CAS Google Scholar
First citationKurnosov, A., Kantor, I., Boffa-Ballaran, T., Lindhardt, S., Dubrovinsky, L., Kuznetsov, A. & Zehnder, B. H. (2008). Rev. Sci. Instrum. 79, 045110.  Web of Science CrossRef PubMed Google Scholar
First citationLaniel, D., Winkler, B., Bykova, E., Fedotenko, T., Chariton, S., Milman, V., Bykov, M., Prakapenka, V., Dubrovinsky, L. & Dubrovinskaia, N. (2020). Phys. Rev. B, 102, 134109.  Web of Science CrossRef ICSD Google Scholar
First citationLengeler, B., Schroer, C. G., Kuhlmann, M., Benner, B., Günzler, T. F., Kurapova, O., Zontone, F., Snigirev, A. & Snigireva, I. (2005). J. Phys. D Appl. Phys. 38, A218–A222.  Web of Science CrossRef CAS Google Scholar
First citationLiermann, H.-P., Konôpková, Z., Morgenroth, W., Glazyrin, K., Bednarčik, J., McBride, E. E., Petitgirard, S., Delitz, J. T., Wendt, M., Bican, Y., Ehnes, A., Schwark, I., Rothkirch, A., Tischer, M., Heuer, J., Schulte-Schrepping, H., Kracht, T. & Franz, H. (2015). J. Synchrotron Rad. 22, 908–924.  Web of Science CrossRef ICSD CAS IUCr Journals Google Scholar
First citationMatos, M., Hoffmann, R., Latgé, A. & Anda, E. V. (1996). Chem. Mater. 8, 2324–2330.  CrossRef CAS Web of Science Google Scholar
First citationMcMahon, M. (2018). Nat. Mater. 17, 858–859.  Web of Science CrossRef CAS PubMed Google Scholar
First citationPalmer, D. C. (2014). CrystalMaker. CrystalMaker Software Ltd, Begbroke, England.  Google Scholar
First citationPetříček, V., Dušek, M. & Palatinus, L. (2006). JANA2006. Institute of Physics, Czech Academy of Sciences, Prague, Czech Republic.  Google Scholar
First citationPetříček, V., Dušek, M. & Palatinus, L. (2014). Z. Kristallogr. Cryst. Mater. 229, 345–352.  Google Scholar
First citationPrescher, C. & Prakapenka, V. B. (2015). High Press. Res. 35, 223–230.  Web of Science CrossRef CAS Google Scholar
First citationRigaku (2019). CrysAlis PRO. Version 171.40. Rigaku Corporation, Wroclaw, Poland.  Google Scholar
First citationRigaku Oxford Diffraction (2020). CrysAlis PRO. Rigaku Oxford Diffraction, The Woodlands, Texas, USA.  Google Scholar
First citationSchroer, C., Agapov, I., Roehlsberger, R., Wanzenberg, R., Brinkmann, R., Weckert, E. & Leemans, W. (2019). PETRA IV: Upgrade of PETRA III to the Ultimate 3D X-ray Microscope, Conceptual Design Report (CDR). https://bib-pubdb1.desy.de/record/426140Google Scholar
First citationSeiboth, F., Brückner, D., Kahnt, M., Lyubomirskiy, M., Wittwer, F., Dzhigaev, D., Ullsperger, T., Nolte, S., Koch, F., David, C., Garrevoet, J., Falkenberg, G. & Schroer, C. G. (2020). J. Synchrotron Rad. 27, 1121–1130.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationSeiboth, F., Kahnt, M., Scholz, M., Seyrich, M., Wittwer, F., Garrevoet, J., Falkenberg, G., Schropp, A. & Schroer, C. G. (2016). Proc. SPIE, 9963, 99630P.  Google Scholar
First citationSeiboth, F., Schropp, A., Scholz, M., Wittwer, F., Rödel, C., Wünsche, M., Ullsperger, T., Nolte, S., Rahomäki, J., Parfeniukas, K., Giakoumidis, S., Vogt, U., Wagner, U., Rau, C., Boesenberg, U., Garrevoet, J., Falkenberg, G., Galtier, E. C., Ja Lee, H., Nagler, B. & Schroer, C. G. (2017). Nat. Commun. 8, 14623.  Web of Science CrossRef PubMed Google Scholar
First citationSheldrick, G. M. (2008). Acta Cryst. A64, 112–122.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationSheldrick, G. M. (2015). Acta Cryst. A71, 3–8.  Web of Science CrossRef IUCr Journals Google Scholar
First citationTateno, S., Hirose, K., Ohishi, Y. & Tatsumi, Y. (2010). Science, 330, 359–361.  Web of Science CrossRef CAS PubMed Google Scholar
First citationVasiukov, D. (2018). PhD thesis. University of Bayreuth, Germany.  Google Scholar
First citationWang, Y., Li, Q., Chen, N., Cheng, J.-M., Li, C.-G., Li, H., Long, Q.-H., Shi, J.-S. & Deng, J.-J. (2016). Chin. Phys. C, 40, 086205.  Web of Science CrossRef 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.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds