2
1. Introduction
There is a continued need for improved 3D imaging methods to analyse complex multi-length-scale
structures of a plethora of technologically or biologically important materials. In this connection, neutron
imaging is of interest due to the good penetration capabilities, the non-continuous dependence of
interaction cross sections with atomic number and the possibilities for adjusting material contrast by
controlling isotope compositions. The primary principle of image formation with neutron imaging is
attenuation contrast from absorption. This is not always the best way to provide contrast, especially for
samples consisting of materials with very similar absorption cross sections or very low absorption. For
this case, it is useful to apply the principle of phase contrast that harnesses the real part of the refractive
index. Phase contrast imaging can be performed in several ways, using various kinds of interferometry
(Pushin et al., 2017; Strobl et al., 2019) or propagation-based phase contrast (Fiori et al., 2006; Paganin et
al., 2019). The benefit of the latter is that no gratings or other specialised elements are needed to measure
the phase shift. Therefore, propagation-based phase contrast imaging has evolved to a very powerful tool
in X-ray imaging (Alloo et al., 2022; Bidola et al., 2017; Wieland et al., 2021; Yu et al., 2021) but has
only been investigated for neutrons in a very few cases (Allman et al., 2000; Jacobson et al., 2004;
Lehmann et al., 2005; McMahon et al., 2003).
When neutrons pass through a sample, they are refracted due to local variations in the refractive index.
Features in the sample with different refractive indices will act as neutron lenses that either focus or
diverge the neutrons from that point in the sample. These changes of the neutron directions are small and
thus difficult to observe at short distances between sample and detector but they do become progressively
easier to detect when increasing the distance between sample and detector. The build-up of phase contrast
requires a coherent beam, meaning that the beam divergence must be smaller than the changes caused by
refractive features. The propagation-based phase contrast imaging only gives relative information about
the sample, which means that a phase-retrieval technique is required in order to obtain the projected
density of the sample (Paganin et al., 2019).
Here we explore propagation-based phase contrast imaging first in a model sample with very low
absorption and secondly in bone, which is a hierarchically structured material.
Bone is replete with blood vessels and cells situated in lacunae interconnected by canaliculi only a few
hundred nm in diameter (Wittig et al., 2022). Together, these form a vast fluid-containing network.
Transport of fluid in bone is very important since the fluid contains nutrients, signal molecules, and ions
essential for the entire body (Cowin & Cardoso, 2015). In addition, liquid transport is proposed to be the
main mechanism of stress sensing in bone, suggesting that the osteocytes sense changes in shear liquid
flow through the canaliculi (van Tol et al., 2020; Robling & Bonewald, 2020; Burger & Klein-Nulend,