A Shack-Hartmann sensor for single-shot multi-contrast imaging with hard X-rays
Tomy dos Santos Rolo Stefan Reich Dmitry Karpov Sergey Gasilov Danays Kunka Edwin Fohtung Tilo Baumbach Anton Plech





dos Santos Rolo et al.]A Shack-Hartmann sensor for single-shot multi-contrast imaging with hard X-rays



^1Karlsruhe Institute of Technology (KIT), Institute for Photon Science and Synchrotron Radiation (IPS), Herrmann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany.
^2New Mexico State University (NMSU), Department of Physics, 1255 N Horseshoe, Las Cruces, NM 88003-8001, USA.
^3Physical-Technical Institute, Tomsk Polytechnic University (TPU), Lenin Avenue 30, 634050 Tomsk, Russia.
^4Current affiliation: Canadian Light Source, 44 Innovation Boulevard Saskatoon, S7N 2V3, Canada.
^5Karlsruhe Institute of Technology (KIT), Institute of Microstructure Technology, Hermann-von-Helmholtz Platz 1, 76344 Eggenstein-Leopoldshafen, Germany.
^6Los Alamos National Laboratory, Los Alamos, NM 87545, USA.
^7Karlsruhe Institute of Technology (KIT), Laboratory for Applications of Synchrotron Radiation (LAS), Engesserstrasse 15, 76131 Karlsruhe, Germany.

tomy.rolo@googlemail.com
anton.plech@kit.edu



An array of compound refractive X-ray lenses (CRL) with 20x20 lenslets, a focal distance of 20 cm and a visibility of 0.93 is presented. It can be used as a Shack-Hartmann sensor for hard X-rays (SHARX) for wavefront sensing and  permits for true single-shot multi-contrast imaging the dynamics of materials with a spatial resolution in the micrometer range, sensitivity on nanosized structures and temporal resolution on the microsecond scale. The object's absorption and its induced wavefront shift can be assessed simultaneously together with information from diffraction channels. This enables the imaging of hierarchical materials. In contrast to the established Hartmann sensors the SHARX has an increased flux efficiency through focusing of the beam rather than blocking parts of it. We investigated the spatiotemporal behavior of a cavitation bubble induced by laser pulses. Furthermore, we validated the SHARX by measuring refraction angles of a single diamond CRL, where we obtained an angular resolution better than 4 rad.











[
    Tomy dos Santos Rolo^1, Stefan Reich^1, Dmitry Karpov^2, Sergey Gasilov^3,4, Danays Kunka^1, Edwin Fohtung^2,6, Tilo Baumbach ^1,7 and Anton Plech^1
    [
========================================================================================================================================================





§ INTRODUCTION


Assessing distinct features on different length scales in hierarchical materials typically involves the sequential application of different imaging techniques, impeding the simultaneous investigation of more than a single hierarchy level at once. Examining the dynamics of the evolution of hierarchical materials on the microsecond time scale has therefore been a considerable challenge. Furthermore, when the interrelation of dynamics on length scales from the millimeter down to the nanometer range is of interest, simultaneous hierarchical imaging  is mandatory.

The broad availability of X-ray sources such as synchrotrons and free electron lasers with very high brightness paved the way to tackle the formidable demands of imaging the structural evolution of hierarchical materials by means of multi-contrast imaging. 
Phase contrast imaging modalities assess the relative variations in the real part of the refractive index (propagation-based phase contrast <cit.>, crystal-based <cit.> and grating-based interferometry <cit.>). 
Grating based interferometry and diffraction-enhanced imaging can additionally provide scattering contrast, thereby extending the accessible information to the full set of scalar wave-matter interactions in transmission geometry <cit.>.
However, the requirement to record several sub-images interrupted by mechanical motion (either with different sample positions, different relative positions of the gratings, or different angles of the analyzer crystal) complicate the investigation of dynamics on the microsecond scale <cit.>. Successful approaches have used a canted interferometer setup to retrieve the phase shifts by Fourier analysis of the Moiré pattern <cit.> and could use polychromatic X-rays <cit.>.

Hartmann sensors for X-rays provide single-shot differential phase-front measurements <cit.>. By using absorption gratings, however, a large fraction of the incident radiation is absorbed in the optical element and can't contribute to the signal and therefore reducing the achievable signal-to-noise ratio for a given integration time.

We show the application of the analogue of a Shack-Hartmann sensor in the hard X-ray spectral range (SHARX), which permits for single-shot imaging of absorption, phase and diffraction contrast. This is quantitatively equivalent to the Hartmann-mask approach <cit.>, or even conceptually similar to methods with structured X-ray beams <cit.>. Due to the use of focusing lenslets <cit.> the local X-ray flux density is increased strongly, leading to an improved signal-to-noise ratio even at shortest integration times. The availability of  diffraction contrast extends the accessible length scales <cit.>. With a spatial resolution of the pitch of the SHARX (here 50 m) structures in the nanometer up to sub-micrometer regime can be imaged.

Since a single image is sufficient for obtaining the three contrast types simultaneously, contrast variations on the microsecond time scale can be used to judge how the structural dynamics on different length scales influence each other and the overall evolution of the material properties <cit.>.





§ MATERIAL AND METHODS



 §.§ Fabrication of the 2D lens array



In Fig. <ref> we present the hard X-ray analogue to a visible-light array of lenslets in the shape of a 3D  hole structure fabricated from a low-absorbing polymer <cit.>. The 3D structure with crossed hollow cylinders of 40 m  in diameter, a periodicity of 50 m  and a total volume of 1 mm^3 has been fabricated by using a state-of-the-art three-dimensional direct laser writing system (Photonic Professional GT, Nanoscribe GmbH). A commercial photoresist (IP-S, Nanoscribe GmbH) and a 25x magnification objective (0.8 NA) were used. The structure was fabricated in dip-in mode, meaning that the objective is immersed in the photoresist. After the writing process, the structure was developed in mr-Dev 600 (micro resist technology GmbH).

Each cylinder acts as a one-dimensional lens and the serial stacking shortens the focal length of the compound refractive lens (CRL)<cit.>. The orthogonal combination of crossed cylinders in the beam propagation direction leads to a focusing of the beam in two directions. This results in well-defined spots that can be used analogously to Hartmann-mask beamlets. The benefits of this arrangement are manifold. The majority of photons from a lenslet aperture are concentrated on the beamlet area leading to a reduced exposure time due to the increased local photon flux density compared to the unfocused primary beam when using gratings or Hartmann-masks. Furthermore, the differential phase sensitivity can be adjusted to a certain extent by changing the sample-to-detector distance. The required lenslet pitches of tens of micrometers for hard X-ray diffraction experiments permit to vary cylinder diameters, grating materials and fabrication methods, and therefore tune the respective optical element for experimental demands. In particular, diffraction contrast is sensitive to length scales that can be tuned by the distance of the object to the detection plane and the array pitch <cit.>.



 §.§ Experimental setup




The experiments were performed at the bending magnet beamline TOPO-TOMO at the synchrotron at KIT (Karlsruhe, Germany). The scheme of the Shack-Hartmann setup is shown in Fig. <ref>. A white beam (filtered by 0.2 mm of aluminum and 0.25 mm of beryllium), which has a beam size defined by mechanical blades (motorized slits) illuminates the lens array homogeneously, which leads to the formation of the beamlets on an X-ray scintillator. The visible light from the scintillator is imaged by a lens system onto a 2D detector. The sample can either be placed before or after the SHARX, which essentially determines the sensitivity (and susceptibility to phase wrapping) on phase contrast and on the length scale in diffraction contrast. Several setups for the imaging and optical characterisation were employed:

(i) The internal microstructure of the SHARX was examined by means of micro-tomography acquired in white beam with a 7 cm distance between the SHARX and the scintillator, which is much smaller than the focal distance of the lens array. A PCO.dimax camera was used to image a 10 m thick LSO:Tb scintillator by a 10x objective (Mitutoyo LWD 10X 0.28) and a 180 mm tube lens, resulting in a ninefold optical magnification. The effective pixel size was 1.22 m. Projections were continuously recorded while rotating the sample with constant angular velocity around an axis perpendicular to the beam propagation directions <cit.>.

(ii) The focal distance of the SHARX was characterized in monochromatic beam (double-multilayer monochromator with 2% bandwidth at 9 keV). As detector an Andor Neo camera coupled to a 200 m LuAG:Ce scintillator via a single objective lens (Nikon Nikkor 85/1.4) was used. The optical magnification was 3.6x, resulting in an effective pixel size of 1.8 m.

(iii) The phase-front reconstruction of the diamond lens was performed in monochromatic beam (8.5 keV). The camera used was the Andor Neo using the same optics as in setup (i) resulting in an effective pixel size of 0.72 m. The distance between the detector and the SHARX was 16 cm and the diamond lens was placed before the SHARX.

(iv) The fast multi-contrast imaging during the ablation process was performed in white beam. The camera was a PCO.dimax, lens-coupled with twofold magnification to a 50 m  thick LuAg:Ce scintillator. The effective pixel size was 5.5 m. For each laser pulse an image sequence was recorded with a frame rate of 15 kHz and an exposure time of 33 s to limit motion blurring. The ablation process occurred in an in-situ chamber, which was placed behind the SHARX close  to the scintillator (7 cm distance). The ablation, in brief, was performed on a silver wire target (0.7 mm diameter) continuously transported through a sealed chamber of 0.4 cm^3 volume and flushed continuously by a water flow to avoid laser shielding by nanoparticles produced at predecessing laser pulses. Pulses from a nanosecond laser (Continuum Minilite I, 1064 nm, 10 mJ) entered the chamber via a lens (38 mm effective focal length in water) to be focused onto the target (for more details see <cit.>). For statistical reasons images from 700 subsequent laser shots were averaged.



 §.§ Image formation


Every lenslet focuses a part of the incident plane wave onto a spot on the detector. After passing the object under investigation, each beamlet is attenuated, displaced or broadened relative to the incident beam direction (see Fig. <ref>) depending on the local transmission, refraction or diffraction properties of the sample, respectively. By using the central result from Fourier optics that a slice of sample placed in the path of a converging spherical wave produces the optical Fourier transform of the part of the object intersected by the wavefront <cit.>, it is possible to intuitively relate the local object properties to the modification of intensity in the focal plane. The angular spectrum of a planar wavefront is a delta function, thus the spots of the unmodified incident beamlets are as sharp as possible, where its size only depends on the optical transfer function of the lenslet aperture.
We have to note that the wavefront in reality is not necessarily planar due to the finite distance and lateral extent of the X-ray source. In this case the smallest achievable focal spot is the demagnified image of the source. Furthermore, the focal distance changes according to optical equations. This effect, however, is  small in our case as the distance from the source (30 m) is two orders of magnitude larger than the focal distance.

In the case of local attenuation in the sample, the direction of the wavefront is not changed, only the magnitude of the local intensity is reduced. Refraction, on the other hand, introduces a local change of direction of the planar wavefront without changing the frequency content of the angular spectrum. The apparent width of the resulting intensity distribution in the focal plane is unchanged compared to the unmodified incident beamlets. Any other wavefront curvature may arise from local scattering on sub-aperture-sized discontinuities of the refractive index in the sample <cit.>. These have a high-frequency component in the angular spectrum. Since the intensity response is the convolution of the optical transfer function with the angular spectrum of the incident wavefront <cit.>, this introduces a broadening of the focal spot. 

Small-angle scattering of objects that fulfil the far field condition adds to the broadening. Here, diffraction contrast averages signals from all the objects that fall into the sensitivity size range <cit.>. By accessing several Fourier orders a crude size discrimination may be achieved <cit.>. In practice, the distinction between phase contrast and diffraction contrast originating from scattering is not strict and depends on the geometry of the experimental setup, as well as discretization effects. In the continuous case, when the lenslet separation tends to zero, the distinction vanishes <cit.>.

In summary, the array of CRLs acts as the X-ray equivalent of a Shack-Hartmann sensor konwn in the optical regime <cit.>, which allows for discrimination of absorption contrast, phase contrast and additionally a set of diffraction contrasts depending on the geometry and detector resolution. Thus, a second level of size information on the nanometer scale is overlaid on the morphology on the micrometer scale.



 §.§ Retrieval of the contrasts: amplitude, phase and diffraction



To retrieve the different contrasts of amplitude, differential phase and diffraction, each spot was fitted separately by a 2D-Gaussian:

    g = h · exp { - 1/2·[(x- μ_x/σ_x^2)^2 + (y-μ_y/σ_y^2)^2] } + o  ,

where h is the height, o the offset, σ_x,y the width and μ_x,y the peak position in horizontal and vertical direction, respectively. Any rotation of the Gaussian distribution was neglected to reduce calculation time and because no pronounced elliptical spot deformation was detected. The offset was only calculated for better fitting, but represents no separate contrast as long as it is just the residual background intensity. 

To determine the changes in spot properties due to the sample the Gaussian fit was performed for two images, one without the sample and one with the sample. This leads to six contrast images for the undisturbed (free) image (height: h_f, offset: o_f, width in horizontal and vertical direction: σ_x,f and σ_y,f, position in horizontal and vertical direction: μ_x,f and μ_y,f) and equivalent for the sample image (h_s, o_s, σ_x,s, σ_y,s, μ_x,s and μ_y,s). The changes where calculated by referring the undisturbed to the sample image for the transmission T, the differential phase dP and the diffraction D in horizontal and vertical direction by the following expressions:

    T    =   h_s/h_f
    
      dP_x = μ_x,s - μ_x,f   and    dP_y = μ_y,s - μ_y,f
    
      D_x = σ_x,s - σ_x,f   and     D_y = σ_y,s - σ_y,f

For calculating the actual beam deflection the differential phase was (in small-angle approximation) divided by the working distance. The latter being the distance between the SHARX or the detector when the sample is placed before the SHARX and the distance between the sample and the detector when the sample is placed between the SHARX and the detector.

Gaussian fits were used instead of Fourier-analysis, since this procedure prevents issues with phase wrapping <cit.>. Additionally, the Fourier analysis as introduced by Wen et al. <cit.> is known to introduce cross-talk between absorption and diffraction contrast <cit.>. However, depending on the numbers of pixels sampling each beamlet, the Fourier analysis is able to probe several Fourier orders, possibly yielding information on several length scales in diffraction contrast as mentioned above.

For wavefront sensing the angular resolution is of primary importance. It depends of the spatial resolution of the detector and the working distance. In setup (i) of section <ref> a resolution of 1.8 m / 20 cm = 9 rad is expected for a shift of a full pixel. Nevertheless, peak shifts smaller than one pixel can easily be detected, leading to a  higher angular resolution.

The reconstruction of the change in the wavefront was done by zonal wavefront estimation introduced by Southwell <cit.>. For better accuracy in phase estimation the modified Southwell algorithm, as proposed in <cit.> was used. Here, additionally to the horizontal and the vertical phase slope, also the diagonal ones are incorporated in the reconstruction.



§ RESULTS AND DISCUSSION




 §.§ Lens array characterization


	
	The focal length of biconcave X-ray lenses is calculated according to the basic optics formula <cit.>:

    f = R/2 · N ·δ ,

where R is the radius of curvature, N the number of stacked biconcave lenses and 1-δ the energy-dependent real part of the index of refraction of the material in the X-ray range. Here, each CRL consists of 20 crossed cylindrical lenses which leads to an effective total of 10 biconcave lenses for each direction. The radius R is 20 m and δ = 5.65 ·10^-6 at 9 keV. The expected focal length therefore is 177 mm. Scanning the distance between the detector and the array (see Fig. <ref> a)) yields the minimal width of the focal spots at (200 ± 15) mm which is in good agreement with the design value. Note that the real composition of the photoresist is unknown and therefore the exact value of δ may vary slightly. Additionally, lens errors may add to the shift. The slight difference in spot size in x- and y-direction originates from different sizes of the source in horizontal and vertical direction. The visibility V (V = (I_max-I_min)/(I_max+I_min) with I_max being the maximum and I_min the minimum intensity in the square zone of each spot) of the spots peaks at the same distance as the effective focal length.

The achieved visibility of 0.93 is considerably higher than typical values for grating interferometry <cit.>. Even in white-beam illumination (centered at 15 keV), we still obtained a high visibility of 0.5. Fig. <ref> b) shows three line plots in x-direction of different spot rows (solid lines, the upper two are shifted for clarity) in monochromatic mode. All of them show high gain factors compared to the free beam intensity (dashed lines). Two points have to be noted. First, not only intensity within the horizontal line is concentrated into the spots, but also in the perpendicular direction. Therefore intensity appears not to be conserved (area below dashed versus solid lines). Second, the minima in Fig. <ref> b) are not the global minima in each zone of the spot. The SHARX consists of crossed cylinders acting as 2D focusing lenses at the exact intersection point, while focusing works only in one direction at positions away from the cylinder cross points. The minima in Fig. <ref> b) reflect such positions leading to an increased background intensities compared to areas of no intensity increase or even intensity decrease. The true visibility, however, was derived from analysing a square area around each spot.  



 §.§ Reconstruction of the phase shift of a diamond lens



The SHARX can be used to analyze the changed wavefront of the X-ray beam in the exit plane of an X-ray optical element. For demonstration we analysed a CRL composed of 8 plano-concave parabolic lenses made from single-crystal diamond. Each unit lens has a design radius of curvature of 200 m at the vertex of the parabola and an entrance aperture of around 0.9 mm <cit.> leading to a focal length of some 2.5 m (at 8.5 keV). Fig. <ref> a) clearly illustrates the effect produced by the CRL inserted in the X-ray beam: the almost rectilinear grid of X-ray spots formed by the SHARX (left part) is bent toward the central axis of the CRL (right part). Shift in meridional and sagittal planes of each beamlet are derived by comparing reference and distorted image. The resulting local refraction angle in the CRL exit plane is shown as a vector field in Fig. <ref> b). As expected, spots at the peripheral part of the CRL experience larger shifts than peaks located near the lens central axis. From these shifts we could easily reconstruct the X-ray phase shift using, for instance, the modified Southwell algorithm <cit.> (see also section <ref>). The calculated phase shift (Fig. <ref> c)) is in reasonable agreement with the prediction based on the design values.

To verify the accuracy of the performed phase front metrology, the deflection of the X-rays introduced by the diamond CRL was derived independently from a microtomography (CT) with the CT rotation axis coinciding with the CRL optical axis<cit.>. The reconstructed distribution of material density was subjected to a segmentation procedure in order to outline the surfaces of the diamond. These in turn were used to calculate the thickness of diamond projected onto the CRL exit plane. Finally, this projected thickness was converted into X-ray deflection angles using the relation β(x,y)=Δ T(x,y)·δ_d, where Δ T(x,y) is the finite difference approximation to the first derivative of the thickness and 1-δ_d the real part of the index of refraction of diamond. The result is shown as dashed line in Fig. <ref> d) and compared to the SHARX deflection angles (red stars and black squares). Both angular shifts match well within a 4 rad confidence and relate to a CRL with effective radius of curvature at the apex of 205 m. A very good agreement between the measured beam deflection by the SHARX and the indirectly estimated deflection from the CT  demonstrates that a precise phasefront metrology can be performed with the SHARX device.

2D representations of the differential phase and diffraction contrast in x and y-direction are shown in Fig. <ref> e). Phase contrast is shown for both directions, reproducing the finding of phase front curvature. The diffraction patterns show no significant signal as expected from a smooth material. However, the edges display an elevated diffraction level. Here two effects appear: the abrupt beam displacement within the individual beamlets, leading to a high curvature of the wave front, and an enhanced roughness due to processing. The first effect is in line with the interdependence of diffraction and phase from objects sized within the spatial resolution <cit.>.



 §.§ In-situ imaging of laser-induced cavitation bubble



Pulsed laser ablation in liquid (PLAL) is employed  to produce ligand-free nanoparticles in (e.g. aqueous) suspension <cit.>. Due to its hierarchical processes spanning several length scales and the fast dynamics of structure formation the  mechanisms and control of the process are still under investigation <cit.>. Within the present nanosecond excitation it is known that ablated material is ejected from the target by phase explosion, while at the same time a plasma is ignited to modify particle formation and foster energy coupling into the water. The latter is easily visible as a millimeter-sized transient cavitation bubble <cit.>. This bubble develops into a hemispherical shape of some 1.5 mm radius within a typical 150 s lifetime. It is known that this bubble contains a major part of the ejected mass from the target as nanoparticles <cit.> and the mutual interaction governs particle ripening <cit.>. 

With the SHARX setup we could measure differential phase and diffraction in addition to the absorption and therefore gain information on the different hierarchical levels.
The bubble is imaged with the spatial resolution of the pitch of the SHARX (50 m). While the absorption and phase shift contain information about the bubble structure, the diffraction channel represents the scattering signal from structures in the nanometer up to the sub-micrometer scale as indicated by the sensitivity graph in fig. <ref>.
Due to the short exposure time an averaging over a number of (nominally reproducible) events is necessary for improving the signal-to-noise ratio.

Fig. <ref> f) shows an absorption contrast radiography without the SHARX for orientation <cit.> over a large field of view. Due to the limited size of the SHARX array we sampled separately two different positions relative to the cavitation bubble, which were merged afterwards to cover the region of interest (see Fig. <ref> f)). Fig. <ref> a)-d) shows the results at maximum bubble expansion for the different contrast channels. The differential phases in both directions clearly show the rim of the bubble, which coincides with the highest phase gradient. The directional sensitivity is unambiguous. A reconstruction of the absolute phase could be calculated as for the diamond lens (not shown here).

As known from investigations with Hartmann masks, the diffraction channel is referred to the scattering of spatially unresolvable structures<cit.>. As explained in theory by Lynch et al.<cit.> the sensitivity for different particle sizes is mainly dependent on the pitch of the spots, the X-ray energy and the distance between the sample and the detector. The sensitivity peak for the present setup is located at a structure size of 300 nm with a FWHM from 110 to 1500 nm <cit.> as shown in the sensitivity curve in Fig. <ref> e). From earlier studies it is known that the main size fraction of particles produced by PLAL is smaller than the present range<cit.>. However, ripening processes produce large particles that enter the visibility interval.

The small observed signals in the diffraction channel in fig <ref> d) (average of the two directional ones as no anisotropic scattering is expected) is of similar spatial distribution as the absorption signal. Therefore it is reasonable to consider a homogeneous filling of the bubble with scattering structures. On the other hand, some residual crosstalk from absorption or phase contrast <cit.> may contribute to the signal and is subject of ongoing investigations.
The high diffraction signal at the bubble boundary is attributed to the crosstalk from an unresolvable sharp phase change <cit.>.

A signal decomposition as in <cit.> could not shed more light on this issue. These signals are visible by Fourier analysis as well<cit.>. A further setup optimization and data analysis are required. It should be kept in mind that the integration time for the presented data only adds up to a total of 14 ms given the 33 s exposure of the individual frames.



§ SUMMARY AND CONCLUSION

Summarizing these results, the Shack-Hartmann approach with a 2D array of X-ray lenses, SHARX, is able to capture multi-contrast images of objects which discern absorption, phase and diffraction. This is quantitatively equivalent to the described Hartmann-mask approach <cit.> or even conceptually similar to other methods with structured X-ray beams <cit.>.  The only observable is the modification of a patterned X-ray field at the detector plane and its relation to the optical properties of the sample, but not to the way the patterning is produced. Sensitivity therefore relates to pure geometric design and pattern visibility.
The Shack-Hartmann sensor therefore increases usable X-ray flux compared to the Hartmann mask. An amplification factor of 6 has so far been achieved but with higher space filling fraction and at sources with smaller source size this can be improved further. Thus exposure times can be considerably lowered to study ultrafast processes. 
On the other hand, spatial resolution relative to the native detector pixel size is reduced. This can nevertheless be compensated for in some cases, where the sample resides closer to the detector (and is static). In that case the beamlets only illuminate a gridded part of the sample. A lateral rastering of the sample can then regain spatial resolution by interleaving several exposures. 

As a refraction-based element the SHARX is chromatic due to the fact that the refractive index of any lens material is energy-dependent. Consequently, the focal point will shift with X-ray energy. Nevertheless, this shift (estimated as 5.9 cm/keV at 9 keV) in our case is comparable to broadening due to the imaging of the source onto the detector as compared to the filtered spectrum with 50% bandwidth. With quasi-monochromatic sources (some 5% bandwidth) the broadening due to shifting the focal length would be about 0.2 m, which is small as compared to the source image in our case (7 m).

In this paper we have demonstrated the use of a 2D array of cylindrical holes inside a refractive polymer, called SHARX, as a compound refractive lens array for hard X-rays. We have produced a lens array of 50 m pitch with a focal length of 200 mm, which is sensitive to angular shifts of the wavefront down to 4 rad at a micrometer spatial resolution. At the same time structural length scales of some 300 nm are probed in diffraction contrast. The device requires low precision of alignment and no moving parts during the data acquisition. The performance is preserved at a relaxed bandwidth of the X-ray beam.

The sensitivity of phase modulation by an object (sample) has been shown to be very high and comparable to other phase-sensitive techniques such as the analyser-based wavefront sensing shown here. Adding the high flux, this approach is particularly useful for imaging irreversible or fast processes in materials or sensing wavefront (changes) of X-ray beams. The efficiency of the method can be further improved by increasing the filling factor for the lenses and focusing precision by shifting to individual parabolic lenses rather than lens cylinders.   

The availability of a diffraction contrast is able to fulfil demands imposed by the requirements of imaging hierarchical structure evolution. The simultaneously observable length scales combine the micrometer range in absorption contrast and the nanometer to sub-micrometer range in diffraction contrast at the spatial resolution of the absorption contrast. This results in the ability to assess contrast variations of interest on the microsecond time scale and to investigate how the structural dynamic influences each other.




§ ACKNOWLEDGMENTS

This research has been supported by the Helmholtz programme topic "from matter to materials and life - MML" and by the German Research Foundation DFG through grant PL 325/8-1. 
The Karlsruhe Nano Micro Facility KNMF is acknowledged for producing the lens array. We are grateful to TISNCM (Troitsk, Russian Federation) for providing  a single crystal diamond CRL for testing. The provision of beamtime at ANKA is acknowledged. We wish to thank T. Müller for support at the beamline and S. Barcikowski and A. Letzel for continuing discussions.  



§ REFERENCES


10
url<#>1#1urlprefixURL 

Wilkins1996
Wilkins S W, Gureyev T E, Gao D, Pogany A and Stevenson A W, 1996 Nature 384 335

Diemoz2012
Diemoz P, Bravin A and Coan P, 2012 Opt. Express 20 2789

Cloetens1996
Cloetens P, Barrett R, Baruchel J, Guigay J P and Schlenker M, 1996 J.
  Phys. D: Appl. Phys. 29 133.

  
Bonse1965
Bonse U and Hart M, 1965 Appl. Phys. Lett. 6 155

Momose2003
Momose A, Kawamoto S, Koyama I, Hamaishi Y, Takai K and Suzuki Y, 2003 Japan J. Appl. Phys. 42 L866

Weitkamp2005
Weitkamp T, Diaz A, David C, Pfeiffer F, Stampanoni M, Cloetens P and Ziegler
  E, 2005 Opt. Express 13 6296.

 
Pfeiffer2006
Pfeiffer F, Weitkamp T, Bunk O and David C, 2006 Nat. Phys. 2 258

Momose2009
Momose A, Yashiro W, Maikusa H and Takeda Y, 2009 Opt. Express 17
  12540.

 
Pfeiffer2008
Pfeiffer F, Bech M, Bunk O, Kraft P, Eikenberry E F, Brönnimann C,
  Grünzweig C and David C, 2008 Nat. Mater. 7 134

Kagias2017
Kagias M, Wang Z, Jefimovs K and Stampanoni M, 2017 Appl. Phys. Lett.
  110 014105

McDonald2009
McDonald S A, Marone F, Hintermüller C, Mikuljan G, David C, Pfeiffer F and
  Stampanoni M, 2009 J. Synchrotron Rad. 16 562.


Momose2011
Momose A, Yashiro W, Harasse S and Kuwabara H, 2011 Opt. Express 19 8423

Wen2010
Wen H H, Bennett E E, Kopace R, Stein A F and Pai V, 2010 Opt. Lett.
  35 1932

Morgan2011
Morgan K S, Paganin D M and Siu K K W, 2011 Opt. Express 19 19781

Rutishauser2012
Rutishauser S, Samoylova L, Krzywinski J, Bunk O, Grünert J, Sinn H,
  Cammarata M, Fritz D M and David C, 2012 Nat. Commun. 3 947

Rand2011
Rand D, Ortiz V, Liu Y, Derdak Z, Wands J R, Tatíček M and Rose-Petruck
  C, 2011 Nano Lett. 11 2678

Vittoria2015SR
Vittoria F A, Endrizzi M, Diemoz P C, Zamir A, Wagner U H, Rau C, Robinson I K
  and Olivo A, 2015 Sci. Rep. 5 16318

Wang2015PRL
Wang H, Kashyap Y and Sawhney K, 2015 Phys. Rev. Lett. 114 103901

Wang2016
Wang H, Kashyap Y and Sawhney K, 2016 Appl. Phys. Lett. 108
  124102

Kagias2016
Kagias M, Wang Z, Villanueva-Perez P, Jefimovs K and Stampanoni M, 2016 Phys. Rev. Lett. 116 093902.


Snigirev1996
Snigirev A, Kohn V, Snigireva I and Lengeler B, 1996 Nature 384
  49

Lengeler1998
Lengeler B, Tümmler J, Snigirev A, Snigireva I and Raven C, 1998 J.
  Appl. Phys. 84 5855

Lyubomirskiy2016
Lyubomirskiy M, Snigireva I, Kohn V, Kuznetsov S, Yunkin V, Vaughan G and
  Snigirev A, 2016 J. Synchrotron Rad. 23 1104.


Lynch2011
Lynch S K, Pai V, Auxier J, Stein A F, Bennett E E, Kemble C K, Xiao X, Lee
  W K, Morgan N Y and Wen H H, 2011 Appl. Opt. 50 4310

Strobl2014
Strobl M, 2014 Sci. Rep. 4 7243.


Ibrahimkutty2015
Ibrahimkutty S, Wagener P, dos Santos Rolo T, Karpov D, Menzel A, Baumbach T,
  Barcikowski S and Plech A, 2015 Sci. Rep. 5 16313

vonFreymann2010
Von Freymann G, Ledermann A, Thiel M, Staude I, Essig S, Busch K and Wegener
  M, 2010 Adv. Funct. Mater. 20 1038

Maruo1997
Maruo S, Nakamura O and Kawata S, 1997 Opt. Lett. 22 132

Rolo2014
dos Santos Rolo T, Ershov A, van de Kamp T and Baumbach T, 2014 Proc.
  Natl. Acad. Sci 111 3921

Reich2017JCIS
Reich S, Schönfeld P, Wagener P, Letzel A, Ibrahimkutty S, Gökce B,
  Barcikowski S, Menzel A, dos Santos Rolo T and Plech A, 2017 J.
  Colloid Interface Sci. 489 106

Reich2017CPC
Reich S, Schönfeld P, Letzel A, Kohsakowski S, Olbinado M, Gökce B,
  Barcikowski S and Plech A, 2017 ChemPhysChem 18 1084

Goodman2005
Goodman J W, Introduction to Fourier Optics (Roberts and Company, 2005)

Malecki2012
Malecki A, Potdevin G and Pfeiffer F, 2012 EPL 99 48001.


Wen2009
Wen H, Bennett E E, Hegedus M M and Rapacchi S, 2009 Radiology 251 910

Koenig2016
Koenig T, Zuber M, Trimborn B, Farago T, Meyer P, Kunka D, Albrecht F, Kreuer
  S, Volk T, Fiederle M and Baumbach T, 2016 Phys. Med. Biol. 61
  3427.

Artzner1992
Artzner G E, 1992 Opt. Eng. 31 1311

Vittoria2015APL
Vittoria F A, Kallon G K N, Basta D, Diemoz P C, Robinson I K, Olivo A and
  Endrizzi M, 2015 Appl. Phys. Lett. 106 224102

Southwell1980
Southwell W, 1980 J. Opt. Soc. Am. 70 998

Pathak2014JO
Pathak B and Boruah B R, 2014 J. Opt. 16 055403

David2002
David C, Nöhammer B, Solak H H and Ziegler E, 2002 Appl. Phys.
  Lett. 81 3287

Polyakov2017
Polyakov S N, Zholudev S I, Gasilov S V, Martyushov S Y, Denisov V N, Terentiev
  S A and Blank V D, 2017 Proc.SPIE 10243 102430X

Gasilov2017OE
Gasilov S, dos Santos Rolo T, Mittone A, Polyakov S, Terentyev S, Farago T,
  Blank V, Bravin A and Baumbach T, 2017 Opt. Express 25 25090

Gasilov2017JSR
Gasilov S, Mittone A, dos Santos Rolo T, Polyakov S, Zholudev S, Terentyev S,
  Blank V, Bravin A and Baumbach T, 2017 J. Synchrotron Rad. 24
  1137

Yashiro2015
Yashiro W and Momose A, 2015 Opt. Express 23 9233.


Zeng2012
Zeng H, Du X W, Singh S C, Kulinich S A, Yang S, He J and Cai W, 2012 Adv.
  Funct. Mater. 22 1333

Zhang2017CR
Zhang D, Gökce B and Barcikowski S, 2017 Chem. Rev. 117 3990

Reich2018APA
Reich S, Göttlicher J, Letzel A, Gökce B, Barcikowski S, dos
  Santos Rolo T, Baumbach T and Plech A, 2018 Appl. Phys. A 124
  71

Ibrahimkutty2012
Ibrahimkutty S, Wagener P, Menzel A, Plech A and Barcikowski S, 2012 Appl.
  Phys. Lett. 101 103104

Letzel2017
Letzel A, Gökce B, Wagener P, Ibrahimkutty S, Menzel A, Plech A and
  Barcikowski S, 2017 J. Phys. Chem. C 121 5356

Kaeppler2014
Kaeppler S, Bayer F, Weber T, Maier A, Anton G, Hornegger J, Beckmann M,
  Fasching P A, Hartmann A, Heindl F, Michel T, Oezguel G, Pelzer G, Rauh C,
  Rieger J, Schulz-Wendtland R, Uder M, Wachter D, Wenkel E and Riess C, Signal Decomposition for X-ray Dark-Field Imaging (Springer International
  Publishing, 2014), 170–177




