Quantitative Neutron Dark-Field Imaging of Milk: A Feasibility Study

: Scattering studies of milk and milk products, which are highly relevant food products on the global market, are often utilized and reported in literature to investigate and understand the subtle microscopic structural differences between dairy samples. These structural features determine the physical properties and ultimately the texture of milk products and, thus, also inﬂuence the consumer’s experience. Small-angle neutron scattering is a prominent example, which enables observations of length scales, which convey proteins and fat globules in food-grade milk. In addition, deuteration enables contrast variations between the constituents of dairy products. In this study, we investigate the potential of probing small-angle neutron scattering from milk samples through quantitative neutron dark-ﬁeld imaging using grating interferometry, to establish the feasibility of studying, in particular, fat globules and milk gel structures with this spatially resolved scattering technique.


Introduction
Dairy milk has a very long history of use and forms a basic part of global nutrition, which is reflected in global production in 2019 of 881 hundred million tons [1]. Thus, the structure of milk during processing is of considerable commercial and consumer interest in the sense of creating a demand for new dairy products; optimizing processing; human nutrition; and oral drug delivery. The most basic and main constituents of milk, apart from water, salts, and whey proteins, are casein micelles (CM) and fat in the form of fat globules (FG) with sizes ranging between 10-10 2 nm for CM and approximately 10 3 nm for FG, respectively. Dairy products contain a hierarchy of length scales in the arrangements of these structures. Milk coagulation in the formation of dairy gels, by, for example, microbiological action, pH variation, or enzymatic action, results in the formation of continuous networks. In milk, the CM is a colloidal and electrostatically stabilized sol, where glycomacropeptide strands, κ-casein, decorate the CM surface in the form of brushes, providing electrostatic or steric repulsion. To enable gelation, the repulsion can either be neutralized or, in the case of rennet gels, the brush strands can be removed by enzymes.

Samples
The milk powders used were products 21,447 and 22,058 manufactured by LACTALIS with 0% (skim milk) and 26% fat (whole milk), respectively. The milk powders were dissolved in deuterated water (1 g of powder for 10 mL D 2 O). The two solutions were mixed in the ratios 5%, 10%, 25%, and 50% in volume of fat-milk solution (we will refer to the samples with respect to these ratios). This corresponds to 1.3 g/L, 2.6 g/L, 6.5 g/L, and 13 g/L of fat, respectively. D 2 O enables a higher neutron scattering length density (nSLD) contrast with casein and even higher with fat.
For the experiments, the deuterated dairy samples were filled into cuvettes as shown in Figure 1. These were arranged on a sample holder in groups that were probed simultaneously, profiting from the spatial resolution of the method. The deuterated milk samples were contained and sealed in these standard quartz cuvettes of 5 mm thickness (t).

Samples
The milk powders used were products 21,447 and 22,058 manufactured LACTALIS with 0% (skim milk) and 26% fat (whole milk), respectively. The milk pow were dissolved in deuterated water (1 g of powder for 10 mL D2O). The two solut were mixed in the ratios 5%, 10%, 25%, and 50% in volume of fat-milk solution (we refer to the samples with respect to these ratios). This corresponds to 1.3 g/L, 2.6 g/L g/L, and 13 g/L of fat, respectively. D2O enables a higher neutron scattering length den (nSLD) contrast with casein and even higher with fat.
For the experiments, the deuterated dairy samples were filled into cuvettes as sh in Figure 1. These were arranged on a sample holder in groups that were pro simultaneously, profiting from the spatial resolution of the method. The deuterated m samples were contained and sealed in these standard quartz cuvettes of 5 mm thick (t). Gelation of the milk samples was achieved using the same initial diluted m powders, by applying two different mechanisms: acidification on the one hand rennet-induction on the other hand. The rennet-induced gel was prepared for fat-(skim) and 25%-fat milk and the acidified gel for fat-free and 10%-fat milk (due to difficulty associated with gelling the 25% milk sample). The acidification was gener by adding Glucono delta-lactone to the solution, which progressively lowers the similar to lactic bacillae. The surface pendent chains of CM, hair-like κ-caseins, bear C H+ charged groups, which induce electrostatic inter-casein repulsion. The low pH chan the COO-H+ charged groups into COOH uncharged groups and thus neutralizes repulsion of the casein MCs and facilitates coagulation. The rennet gel uses re enzymes, which shave off the hair-like κ-casein, suppressing repulsion as well. GP (P "Présure") and GA (A for "Acide") in Figure 1 correspondingly denote rennet and gels, respectively. Figure 2 shows the principal set-up of a Talbot-Lau grating interferometer (TLI). TLI consists of a source grating, G0, a phase grating, G1, and an analyzer grating, G2. G an absorption grating that is positioned downstream of the source to generate mul collimated beams, providing a partially coherent beam. G1 is a phase grating generates an interference pattern, called a Talbot pattern, further downstream. G Gelation of the milk samples was achieved using the same initial diluted milk powders, by applying two different mechanisms: acidification on the one hand and rennet-induction on the other hand. The rennet-induced gel was prepared for fat-free (skim) and 25%-fat milk and the acidified gel for fat-free and 10%-fat milk (due to the difficulty associated with gelling the 25% milk sample). The acidification was generated by adding Glucono deltalactone to the solution, which progressively lowers the pH, similar to lactic bacillae. The surface pendent chains of CM, hair-like κ-caseins, bear COO-H+ charged groups, which induce electrostatic inter-casein repulsion. The low pH changes the COO-H+ charged groups into COOH uncharged groups and thus neutralizes the repulsion of the casein MCs and facilitates coagulation. The rennet gel uses rennet enzymes, which shave off the hair-like κ-casein, suppressing repulsion as well. GP (P for "Présure") and GA (A for "Acide") in Figure 1 correspondingly denote rennet and acid gels, respectively. Figure 2 shows the principal set-up of a Talbot-Lau grating interferometer (TLI). The TLI consists of a source grating, G 0 , a phase grating, G 1 , and an analyzer grating, G 2 . G 0 is an absorption grating that is positioned downstream of the source to generate multiple collimated beams, providing a partially coherent beam. G 1 is a phase grating that generates an interference pattern, called a Talbot pattern, further downstream. G 2 is another absorption grating that is positioned at a fractional Talbot distance and matches the period of the interference pattern. The geometry of G 0 and G 1 is such that the interference patterns of the individual coherent beams from G 0 add constructively at G 2 . The imaging detector system is positioned directly after G 2 and resolves the local interference pattern pixel-wise through the phase stepping of one of the gratings, typically G 0 .

Dark-Field Contrast Imaging
another absorption grating that is positioned at a fractional Talbot distance and matc the period of the interference pattern. The geometry of G0 and G1 is such that interference patterns of the individual coherent beams from G0 add constructively at The imaging detector system is positioned directly after G2 and resolves the l interference pattern pixel-wise through the phase stepping of one of the gratings, typic G0. Small-angle neutron scattering from a sample convolutes with the inten modulation induced by the gratings. This results in the damping of the interfere fringes at G2, i.e., a loss of visibility, . The relative loss of visibility, ⁄ , is referre as the DFI contrast. The scattering distribution at G2 varies with the sample distance and the wavelength , and it can be shown that these parameters, together with the be modulation period , define the real-space length probed in the scattering structure [ This length is referred to as instrumental or probed correlation length, , in full anal to the spin-echo length in SESANS. It can be written as: where coincides with the period of G2 in the TLI setup. Correspondingly, the DFI contrast can be described quantitatively with regard the real-space correlation function ( ) of a scattering structure: where Σ is the total small-angle neutron scattering cross section, containing the volu fraction , the nSLD contrast ∆ , the characteristic structure size (often referred t the specific correlation length of a structure), and the utilized wavelength , while the sample thickness [19]. The real-space correlation function, ( ), which is a o dimensional projection of the three-dimensional scattering-length density correla function of the scattering structure, is related to the small-angle scattering intensity of the structure through the Hankel transform [26]: where is the zeroth-order Bessel function of the first kind, and = (4π sin )⁄ is magnitude of the scattering vector, with being the scattering angle. The character structure size can be expressed as: Small-angle neutron scattering from a sample convolutes with the intensity modulation induced by the gratings. This results in the damping of the interference fringes at G 2 , i.e., a loss of visibility, v s . The relative loss of visibility, v s /v o , is referred to as the DFI contrast. The scattering distribution at G 2 varies with the sample distance L s and the wavelength λ, and it can be shown that these parameters, together with the beam modulation period p, define the real-space length probed in the scattering structure [19]. This length is referred to as instrumental or probed correlation length, ξ, in full analogy to the spin-echo length in SESANS. It can be written as:

Data Acquisition
where p coincides with the period of G 2 in the TLI setup. Correspondingly, the DFI contrast can be described quantitatively with regards to the real-space correlation function G(ξ) of a scattering structure: where Σ s is the total small-angle neutron scattering cross section, containing the volume fraction ϕ, the nSLD contrast ∆ρ, the characteristic structure size ζ (often referred to as the specific correlation length of a structure), and the utilized wavelength λ, while t is the sample thickness [19]. The real-space correlation function, G(ξ), which is a one-dimensional projection of the three-dimensional scattering-length density correlation function of the scattering structure, is related to the small-angle scattering intensity I(q) of the structure through the Hankel transform [26]: where J 0 is the zeroth-order Bessel function of the first kind, and q = (4π sin θ)/λ is the magnitude of the scattering vector, with θ being the scattering angle. The characteristic structure size ζ can be expressed as: Appl. Sci. 2022, 12, 833 5 of 13

Data Acquisition
The TLI was applied at the cold neutron beamline, BOA, at the Swiss neutron source, SINQ, of the Paul Scherrer Institute [27]. The TLI was a symmetric setup, using the first fractional Talbot distance. The TLI worked at the effective neutron wavelength of 4.1 Å and utilized the cold spectrum of the beamline, which was cut off at short wavelengths by a cooled beryllium filter installed upstream. The detailed parameters of TLI are presented at Table 1. Inter-grating distance (mm) Period of gratings (µm) p g (p 0 = p 1 = p 2 ) 50 Height of gratings (µm) Duty cycle of gratings G 0 has the one-dimensional line structure of gadolinium oxysulfide (Gadox, Gd 2 O 2 S), G 1 has the two-dimensional checker-board structure of silicon, and G 2 has the twodimensional mesh-grid structure of Gadox. Each grating has the fabricated-structure area of 100 mm × 100 mm on a 15 cm diameter silicon wafer. G 0 and G 2 were fabricated by a Gadox-particle-filling method [28], and G 1 was fabricated by silicon deep wet etching.
The detector used was an Andor iKon-M CCD camera operated with a standard PSI Midi-box detector design, with a mirror and commercial optical lens system (100 mm Nikkor). The scintillator was a 200 µm thick LiF:ZnS screen. The image area was 1024 × 1024 pixels corresponding to an approximately 67 × 67 mm 2 field of view.
The phase-stepping method was adopted by a G 1 scan performed in eight phase steps for each DFI image. The exposure time was 30 s for each raw image, and three images were recorded per phase step, which were combined by a median filter operation during data reduction. The DFI contrast was extracted from the resulting images pixel by pixel using Fourier analysis. In addition, the sample distance, L s , was scanned to probe the local correlation function, G(ξ), for the correlation length ξ (Equation (1)) from 130 nm to 3.13 µm.

Data Analyses
The deuterated milk samples comprised of CM and FG scattering structures, which are both spherical entities of colloidal dimensions, in aqueous solvent. The radius of CM is typically less than 70 nm within a total range of 20-200 nm, and FG are expected to fall into a size range of 100-15,000 nm. Both are highly polydisperse in the solvent [6,9]. The CM range at the lower limit of the resolved range and provided a lower nSLD contrast than the FG in the used solvent. Thus, they were not expected to be resolved, at least when they were not aggregated by gelation.
In particular, the FG provided a high nSLD contrast, as expected. This also manifested in a significant DFI contrast increase for samples of increasing fat content. The scattering of the FG can be modeled assuming hard spheres, a model applied regularly in quantitative DFI analyses [19,20]. However, it is known from the literature that the FG do not appear as being monodisperse, but rather as being polydisperse, a situation that has, in contrast, not been modelled in the context of DFI so far. The polydispersity of spheres can be considered using a model of spheres weighted by a size distribution P(r), which is expressed by [26]: According to the literature for the FG, a log-normal size distribution can be assumed according to [8] with: where µ is the mean, and σ is the standard deviation of sphere radius r. For the log-normal distribution, the median radius of the sphere is r 0 = e µ , and the mean radius of the sphere is r = r 0 e σ 2 /2 . In addition, particularly for high volume fractions of fat, a structure factor of interparticle correlations must be accounted for in modelling G(ξ) of the system. For this, we first consider the Fourier space description of small-angle scattering, where the scattered intensity distribution, I(q), of a dense system of hard spheres with radius r, sphere volume V(r), volume fraction ϕ, and nSLD contrast ∆ρ is described by: where F is the form factor, and S is the structure factor. The form factor of a sphere of is: The structure factor of such system, describing the scattering contribution arising due to the arrangement of the spheres with respect to each other, can be written as the so-called Percus-Yevick expression [29,30]: where the function Y(x) is given by with Considering, on the other hand, the gelled samples, it is widely accepted that, during gelation, CMs aggregate and form a self-affine network with characteristic fractal dimensions in a manner analogous to colloidal gels [17], as shown in Figure 3.

= .
Considering, on the other hand, the gelled samples, it is widely accepted that gelation, CMs aggregate and form a self-affine network with characteristic dimensions in a manner analogous to colloidal gels [17], as shown in Figure 3. As outlined above the electrostatic repulsion between CMs is diminished, mechanistically, through acidification, or through shaving of surface κ-caseins, leading to a decrease in colloidal stability and aggregation and the formation of a network structure with a number of physical factors determining the scaling of the network [31]. The structural characteristics of a fractal network such as milk gel can be expressed by coarseness and space-filling capacity related to the fractal dimension D f . These are parametrized by a, the characteristic length of inhomogeneity, and the Hurst exponent, H, which is related to fractal dimension through [8]: To model the fractal network of the CM of the gels, we employed a self-affine-structure model referred to as a random two-phase media model [26], reading: where Γ is the gamma function, and K is the modified Bessel function of the second kind. However, to better reflect the structure of the network of CM, the model was modified by substituting the exponential function (e −( ξ 2a ) ) for the modified Bessel function of the second kind (K H+1/2 ξ a ), accounting in particular for the long length-scale correlations of the fractal network of CM according to Ref. [8].
Thus, we used: for modelling the data of the skim milk gels. Figure 4 shows conventional neutron attenuation contrast images (AI) alongside darkfield contrast images (DFI) recorded at two different instrumental settings and probed lengths (130 nm ad 3.13 µm) of the deuterated milk and milk gel samples. The conventional AI images, which, in terms of contrast, do not depend on the changing experimental settings, display only minor contrast variations with respect to fat content and gelation. This is clearly different for the DFI images, where the contrast increases obviously with increasing fat content, as well as with an increasing probed correlation length, ξ. This is the first indication that the contribution of the scattering by FG dominates the scattering signal and that this scattering contrast is picked up for all samples, respectively.  In order to analyze the detected local scattering, i.e., the dark-field contrast over probed correlation length, the DFI curves were extracted for different regions of inte (ROI). These curves can be modelled and fitted according to Section 2.4. In Figure 5a, plotted DFI contrast values are median values of ROIs consisting of 60 × 180 pixel indicated exemplarily for two samples in Figure 4. Figure 5 displays the correspond DFI curves of the deuterated milk samples. Initial visual inspection suggested dominating structure sizes of up to approximately 1-2 μm by observing where the cur became flat, as ( ) tended to zero (compare Equation (2)). The increasing contrast, the deeper levels of flat tails, indicates increased total scattering cross sections, Σ , du increasing volume fractions of the dominating scatterers, i.e., FG. The dashed l represent initial fits with a simple isolated monodisperse hard-sphere model [19,20 and Table 2 presents the fit results in terms of fitted parameters. It is observed that s fits, as expected, do not fully account for the measured slopes, which is especially the c where these turn from a steep decay to flat tails. The more extended bending region in measured curves shows the influence of polydispersity. The solid lines, theref represent subsequent nonlinear least square fits with the polydisperse hard-sphere mo where the projected real-space correlation function ( ) is calculated by a numer Hankel transform, according to Ref. [32] (Equation (3)), of the analytic expression of scattered intensity according to Equations (5)- (12). The model thus accommodates form factor (FG) of the spheres, i.e., the hard-sphere structure factor, for their h concentration and thus correlation, as well as their polydispersity. These fits result in g agreement with the measured data of the DFI contrast for different probed correla lengths in the ROIs (Figure 5a) and result in the fitted structural parameters for the as shown in Table 3. The goodness of fits of FGs in the deuterated milk for the monopolydisperse hard-sphere models were considered using R 2 and the adjusted evaluation. The polydisperse hard-sphere model improves the quality of fit, in term both (R 2 ) and adjusted R 2 , compared to the monodisperse hard-sphere model, as sho in Table 3. Both reach values > 0.95. We also found that the use of a structure factor d not further improve these values, which might be expected due to the low actual volume fraction. In order to analyze the detected local scattering, i.e., the dark-field contrast over the probed correlation length, the DFI curves were extracted for different regions of interest (ROI). These curves can be modelled and fitted according to Section 2.4. In Figure 5a, the plotted DFI contrast values are median values of ROIs consisting of 60 × 180 pixels as indicated exemplarily for two samples in Figure 4. Figure 5 displays the corresponding DFI curves of the deuterated milk samples. Initial visual inspection suggested the dominating structure sizes of up to approximately 1-2 µm by observing where the curves became flat, as G(ξ) tended to zero (compare Equation (2)). The increasing contrast, i.e., the deeper levels of flat tails, indicates increased total scattering cross sections, Σ s , due to increasing volume fractions of the dominating scatterers, i.e., FG. The dashed lines represent initial fits with a simple isolated monodisperse hard-sphere model [19,20,26], and Table 2 presents the fit results in terms of fitted parameters. It is observed that such fits, as expected, do not fully account for the measured slopes, which is especially the case where these turn from a steep decay to flat tails. The more extended bending region in the measured curves shows the influence of polydispersity. The solid lines, therefore, represent subsequent nonlinear least square fits with the polydisperse hard-sphere model, where the projected real-space correlation function G(ξ) is calculated by a numerical Hankel transform, according to Ref. [32] (Equation (3)), of the analytic expression of the scattered intensity according to Equations (5)- (12). The model thus accommodates the form factor (FG) of the spheres, i.e., the hard-sphere structure factor, for their high concentration and thus correlation, as well as their polydispersity. These fits result in good agreement with the measured data of the DFI contrast for different probed correlation lengths in the ROIs (Figure 5a) and result in the fitted structural parameters for the FGs as shown in Table 3. The goodness of fits of FGs in the deuterated milk for the mono-and polydisperse hard-sphere models were considered using R 2 and the adjusted R 2 evaluation. The polydisperse hard-sphere model improves the quality of fit, in terms of both (R 2 ) and adjusted R 2 , compared to the monodisperse hard-sphere model, as shown in Table 3. Both reach values > 0.95. We also found that the use of a structure factor does not further improve these values, which might be expected due to the low actual fat volume fraction.   Table 3. Structural parameters and goodness of fits for fat globules (FG) in the deuterat samples with different fat contents derived using a model for polydisperse interacting (Equations (7)- (13)). The median radius, , and the mean radius, 〈 〉, of FG are calculated by ext the structural parameters from the fitted log-normal polydispersity, chosen based literature about corresponding systems [33]. The radii are calculated from the m and the standard deviation of the log-normal distribution, respectively. We fou the median radius and the mean radius of FG both increased with an increase in v fraction. The mean size of FG was found to increase from approximately 0.36 approximately 1.36 μm in mean diameter (2〈 〉) when considering fat concentratio to 50% fat milk.

Milk Sample
When considering the homogeneity of the samples by choosing several smalle from top to bottom for each sample, as this imaging technique enables, a subtle in of total scattering from bottom to top was found. Figure 6 shows an evaluation of th neutron scattering with an increment of approximately 2.5 mm from the bottom to of the samples. The corresponding ROIs are indicated in the exemplary image on th hand side of Figure 6. The total neutron scattering increases clearly, in particular   Table 3. Structural parameters and goodness of fits for fat globules (FG) in the deuterated milk samples with different fat contents derived using a model for polydisperse interacting spheres (Equations (7)-(13)). The median radius, r 0 , and the mean radius, r, of FG are calculated by extracting the structural parameters from the fitted log-normal polydispersity, chosen based on the literature about corresponding systems [33]. The radii are calculated from the mean µ and the standard deviation σ of the log-normal distribution, respectively. We found that the median radius and the mean radius of FG both increased with an increase in volume fraction. The mean size of FG was found to increase from approximately 0.36 µm to approximately 1.36 µm in mean diameter (2r) when considering fat concentrations of 5 to 50% fat milk.

Milk Sample
When considering the homogeneity of the samples by choosing several smaller ROIs from top to bottom for each sample, as this imaging technique enables, a subtle increase of total scattering from bottom to top was found. Figure 6 shows an evaluation of the total neutron scattering with an increment of approximately 2.5 mm from the bottom to the top of the samples. The corresponding ROIs are indicated in the exemplary image on the right-hand side of Figure 6. The total neutron scattering increases clearly, in particular for the samples with higher fat concentrations. This indicates the presence of larger volume fractions and/or larger sizes of FG towards the top, which is intuitive, considering the lower density of fat. The fits do not, however, provide a clear conclusion about increasing radii towards the top of the samples. It is, therefore, possible that FG concentrates at the top without coalescing, due to a persistent slight repulsion.
Appl. Sci. 2022, 12, x FOR PEER REVIEW Figure 6. Total neutron scattering dependence on sample height. Right-hand-side illu illustrates the ROIs utilized for height dependent data points.
The results of the deuterated milk gel samples are presented in Figure 7. D analyses strategies and models were required to be chosen for the skim milk gels fat milk gels. This was due to the fact that once again the FG signal dominated in of the fat milk gels. The fractal network of CM in the skim milk gels was modelle an exponential function substituting the Bessel function in Equation (13), as desc Section 2.4. Good fits were achieved this way for both the rennet gel (GA) and the a gel (GP). The extracted structural parameters of the fractal CM networks are summ in Table 4 and appear in reasonable agreement with the literature [34,35].

Gel Sample
(μm) Adj The results of the deuterated milk gel samples are presented in Figure 7. Different analyses strategies and models were required to be chosen for the skim milk gels and the fat milk gels. This was due to the fact that once again the FG signal dominated in the case of the fat milk gels. The fractal network of CM in the skim milk gels was modelled using an exponential function substituting the Bessel function in Equation (13), as described in Section 2.4. Good fits were achieved this way for both the rennet gel (GA) and the acidified gel (GP). The extracted structural parameters of the fractal CM networks are summarized in Table 4 and appear in reasonable agreement with the literature [34,35].
Appl. Sci. 2022, 12, x FOR PEER REVIEW 10 Figure 6. Total neutron scattering dependence on sample height. Right-hand-side illustra illustrates the ROIs utilized for height dependent data points.
The results of the deuterated milk gel samples are presented in Figure 7. Differ analyses strategies and models were required to be chosen for the skim milk gels and fat milk gels. This was due to the fact that once again the FG signal dominated in the c of the fat milk gels. The fractal network of CM in the skim milk gels was modelled us an exponential function substituting the Bessel function in Equation (13), as described Section 2.4. Good fits were achieved this way for both the rennet gel (GA) and the acidi gel (GP). The extracted structural parameters of the fractal CM networks are summari in Table 4 and appear in reasonable agreement with the literature [34,35].    For the acidified gel (GA Skim) and the rennet gel (GP Skim), one can notice at first sight that the curves are flatter than those representing the results for milk, as expected from the spatial correlations of a fractal network. The typical size of inhomogeneities was found to be approximately 0.2170 µm and 0.3545 µm for the two gels, respectively. The Hurst exponent, representing the space-filling capacity, i.e., the dimensionality of the fractal structure, assists in quantifying the fractal dimension, D f , as 2.7598 for the acid gel (GA Skim) and 2.5427 for the rennet gel (GP Skim).
For the fat milk gels, GA Fat and GP Fat, respectively, the scattering of which was dominated by the 10-and 25%-fat milk content, only the assumption of an unaltered fractal network (as compared to the skim gels) enabled us to use an approach to characterize the FG structure. This is enabled, in a crude force approach, by dividing the corresponding curves by those of the skim milk gels and subsequently fitting a polydisperse hard-sphere model with a hard-sphere structure factor, in analogy to the fat milk samples. The results are shown in Figure 7 and are summarized in Table 5. The results suggest FG sizes approximately double those of the milk samples with a corresponding fat content, namely approximately 1.5 µm diameter.

Conclusions
Dark-field imaging was applied to study the structure of deuterated milk and milk gel samples using a symmetric Talbot-Lau neutron-grating interferometer. In order to model the data, a numerical Hankel transform of the analytic expression of the smallangle scattering intensity in the case of polydisperse spheres was applied. This enabled us to fit the projected real-space correlation function resulting from the quantitative DFI scans, taking into account a polydisperse size distribution for the form factor. Moreover, a structure factor was considered, but showed little impact with regard to the relatively low concentration of fat globules. As a result, the mean size and log-normal size distribution of the fat globules, as well as their dependence on concentration, could be quantified. In addition, the mapping of the total scattering cross-section as a function of sample height revealed a subtle gradient in the concentration of fat globules towards the top surface of the samples, increasing with fat concentration. Moreover, the structure of the milk gel could be modelled as a fractal network with long-range correlations for skim milk gels related to the casein micelle aggregation enabled by acid and rennet enzymes, respectively. Slightly different structure parameters and fractal dimensions were modelled. For fat milk gels, the fat globules in the fat milk gels were analyzed by approximating their scattering contribution through normalization with the respective skim milk gel data. This approximation implied somewhat larger fat globules in the gels than in the milk samples with the same fat concentrations.
In conclusion, this paper has demonstrated that dark-field imaging data enables the extraction of quantitative information for complex systems involving polydispersity and structure-factor contributions. In particular, quantitative dark-field contrast neutron imag-ing was found to be a well-suited quantitative tool for investigations of dairy system structures. The proven ability to analyze local variations through the imaging approach allows us to add not only the identification but also the assessment of macroscopic variations of microstructures resulting from external stimuli or processing and mixing of constituents to the body of research in this field. This approach holds the promise to be extended to other soft-matter and material science problems.