Skip to content

Commit

Permalink
final
Browse files Browse the repository at this point in the history
  • Loading branch information
fzimmermann89 committed Nov 5, 2021
1 parent 7a9eec2 commit 325cd37
Show file tree
Hide file tree
Showing 7 changed files with 41 additions and 44 deletions.
2 changes: 1 addition & 1 deletion Tex/conclusion.tex
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
\chapter{Conclusion and Outlook}
There is still an ongoing discussion about the feasibility of \textit{Incoherent Diffractive Imaging} and about in which areas it can provide advantages about already established methods. It has not yet been experimentally proven that the assumptions made about the statistics of inner shell fluorescence are valid at the atomic level. Under the assumption of chaotic light, classical wave theory can illustrate the basic working principle. In the quantum mechanical description of the underlying two-photon interference, indistinguishability between the quantum paths leading from the emission of two photons to simultaneous detection in two detectors creates the spatial bunching in the fluorescence patterns. Both explanations predict that a successful IDI experiment has to be designed with few distinguishable paths or modes.
There is still an ongoing discussion about the feasibility of \textit{Incoherent Diffractive Imaging} and its potential advantages over already established methods. It has not yet been experimentally proven that the assumptions made about the statistics of inner shell fluorescence are valid at the atomic level. Under the assumption of chaotic light, classical wave theory can illustrate the basic working principle. In the quantum mechanical description of the underlying two-photon interference, indistinguishability between the quantum paths leading from the emission of two photons to simultaneous detection in two detectors creates the spatial bunching in the fluorescence patterns. Both explanations predict that a successful IDI experiment has to be designed with few distinguishable paths or modes.

The present work contributes simulations to validate theoretical considerations about the Signal-to-Noise characteristics and to optimize the experimental setup. Furthermore, previously published results indicating that IDI can extract spatial information out of X-ray fluorescence intensity patterns were validated.

Expand Down
43 changes: 21 additions & 22 deletions Tex/experiment.tex

Large diffs are not rendered by default.

Binary file modified Tex/images/sim_foil5umCu_shared.pdf
Binary file not shown.
7 changes: 4 additions & 3 deletions Tex/main.tex
Original file line number Diff line number Diff line change
Expand Up @@ -251,12 +251,13 @@

\begin{minipage}[b][12cm]{\textwidth}
\begin{center}
\begin{tabular}{m{12cm}}
\begin{tabular}{m{14.5cm}}
\large{Hiermit erkläre ich, dass ich die vorliegende Arbeit selbstständig und eigenhändig sowie ohne unerlaubte fremde Hilfe und ausschließlich unter Verwendung der aufgeführten Quellen und Hilfsmittel angefertigt habe. } \\
\\
\\
\\
\midrule
\textbf{Datum / Unterschrift} \\
\textbf{Unterschrift} \\
\end{tabular}
\end{center}
\end{minipage}
Expand Down Expand Up @@ -324,7 +325,7 @@ \chapter*{List of Abbreviations}


%%%%%%%%%%%%%%%% KAPITEL %%%%%%%%%%%%%%%%%
\setlength{\parskip}{1mm plus 0.25mm minus 0.25mm}
\setlength{\parskip}{1mm plus 0.25mm minus 0.2mm}
\pagenumbering{arabic}
\pagestyle{headings}
\fancyhf{}
Expand Down
23 changes: 11 additions & 12 deletions Tex/simulation.tex
Original file line number Diff line number Diff line change
Expand Up @@ -76,7 +76,7 @@ \section{Detector Effects}
where $\Delta x$, $\Delta y$ denote the distance of the pixel to the sampled photon center. Afterwards, a Gaussian readout noise is added. The effect of this degradation on the spectrum is illustrated in \fref{fig:degrad}.

\section{Reconstruction of the Structure Factor}
In the approximation introduced in \fref{sec:idi}, the normalized intensity correlation $g^2(\Delta q)-1$ has to be calculated to recover the structure factor from the simulated fluorescence speckle patterns. The task can be split up into two logical units, the calculation of the autocorrelation and the application of a sufficient normalization.
In the approximation introduced in \fref{sec:idi}, the normalized intensity correlation $g^2(\Delta\vec{q})-1$ has to be calculated to recover the structure factor from the simulated fluorescence speckle patterns. The task can be split up into two logical units, the calculation of the autocorrelation and the application of a sufficient normalization.

\subsection{Calculation of the Autocorrelation}
\label{sec:corr}
Expand Down Expand Up @@ -205,8 +205,7 @@ \subsection{Photon Counting}

\subsection{Accessible Reciprocal Space}

As IDI is based on $g^2(\Delta \vec{q})$, for a set of accessible $\vec{q}$ determined by the experimental geometry ($k$, detector size and distance), IDI can give information about $S(\vec{q})$ at higher $\left|\vec{q}\right|$ than a traditional scattering setup (see also \fref{fig:scatteringvectors}), increasing the numerical aperture, in a small angle regime up to a factor of two. With a flat detector, compared IDI can give access to a three-dimensional volume in reciprocal space, as shown in \fref{fig:accesibleq} and, with greater $q_z$ coverage the greater the curvature of the Ewald sphere is. This, for example, gives access to multiple Bragg peaks in a single crystal experiment.
as shown in the simulation in \fref{fig:accesiblebraggq}, but leads to nontrivial shape of the accessible volume.
As IDI is based on $g^2(\Delta \vec{q})$, for a set of accessible $\vec{q}$ determined by the experimental geometry ($k$, detector size and distance), IDI can give information about $S(\vec{q})$ at higher $\left|\vec{q}\right|$ than a traditional scattering setup (see also \fref{fig:scatteringvectors}), increasing the numerical aperture, in a small angle regime up to a factor of two. With a flat detector IDI can give access to a three-dimensional volume in reciprocal space, as shown in \fref{fig:accesibleq}, with greater $q_z$ coverage the greater the curvature of the Ewald sphere is. This, for example, gives access to multiple Bragg peaks in a single crystal experiment as shown in the simulation in \fref{fig:accesiblebraggq}, but leads to nontrivial shape of the accessible volume.

\begin{figure}
\centering
Expand Down Expand Up @@ -252,14 +251,14 @@ \subsection{Comparison of Normalization Techniques}
\subsection{Parameters Influencing the Signal and Noise Characteristics}
\paragraph{Detector Size}

To assess the influence of the number of pixels of a detector on the SNR, a simulation for a 1\,\textmu m thick Copper foil in an 100\,nm FWHM focus was performed. The detector size was varied from 64x64 to 3072x3072 pixels and always placed at the same distance of 1\,m, keeping the mean photon count per pixel constant. The number of correlation averaged over in the reconstruction increases linearly with the number of pixels. For the SNR calculations, the signal is defined as peak intensity, the noise as the standard deviation over independent simulations. As shown in \fref{fig:SNRdetsize}, under these conditions, the SNR is proportional to the square root of the number of pixels of the detector.
Outside of the center area of the reconstruction, the number of different pixel pairs resulting in the same $\Delta{q}$ varies. In a small angle setup, the number of pairs resulting from a finite detector decrease linear in the distance from the center of the reconstruction. In a wide-angle setup, i.e., for a crystal sample, the varying solid angle covered by each pixel results in a nontrivial dependence of the number of pairs and hence of the SNR on $\Delta{q}$ as shown in \fref{fig:pairnoise}.
To assess the influence of the number of pixels of a detector on the SNR, a simulation for a 1\,\textmu m thick Copper foil in an 100\,nm FWHM focus was performed. The detector size was varied from 64x64 to 3072x3072 pixels and always placed at the same distance of 1\,m, keeping the mean photon count per pixel constant. The number of correlation pairs averaged over in the reconstruction of a particular $\vec{q}$ increases linearly with the number of pixels. For the SNR calculations, the signal is defined as peak intensity, the noise as the standard deviation over independent simulations. As shown in \fref{fig:SNRdetsize}, under these conditions, the SNR is proportional to the square root of the number of pixels of the detector.
Outside of the center area of the reconstruction, the number of different pixel pairs resulting in the same $\Delta{q}$ varies. In a small angle setup, the number of pairs resulting from a finite detector decrease linearly in the distance from the center of the reconstruction. In a wide-angle setup, i.e., for a crystal sample, the varying solid angle covered by each pixel results in a nontrivial dependence of the number of pairs and hence of the SNR on $\Delta{q}$ as shown in \fref{fig:pairnoise}.




\paragraph{Number of Images}
As shown in \fref{fig:SNRNimages} if only considering shot and phase noise, uncorrelated between different images, the SNR scales with the square root of the number of images. This can be used to estimate the number of shots necessary to achieve an SNR $>5$ (Rose Criterion) to be able to verify IDI as an imaging method experimentally. Correlated noise by systematic noise sources does not follow this relation, and should be minimized, for example, by masking out affected areas of the detector.
As shown in \fref{fig:SNRNimages} if only considering shot and phase noise, uncorrelated between different images, the SNR scales with the square root of the number of images. This can be used to estimate the number of shots necessary to achieve an SNR\,$>5$ (Rose Criterion) to be able to verify IDI as an imaging method experimentally \cite{rose}. Correlated noise by systematic noise sources does not follow this relation, and should be minimized, for example, by masking out affected areas of the detector.

\paragraph{Number of Modes}
To Illustrate the effect of the number of modes on the signal and the SNR, simulations were performed by averaging the intensity over $M$ realizations of the random phases. In the results of the simulation of a crystalline sample, (\fref{fig:modes}) the predicted $1/M$ scaling of the signal can be seen. Furthermore, as expected, in \fref{fig:SNRNimages} the SNR shows the same scaling.
Expand Down Expand Up @@ -287,9 +286,9 @@ \subsection{Parameters Influencing the Signal and Noise Characteristics}

A simulation of a cubic single crystal with a simple cubic lattice of varying size from 20$^3$ to 200$^3$\,atoms was performed. The lattice constant is chosen as 5.7\,\AA, the fluorescence energy as 8\,keV. The simulated detector has 1024x1024 50\,µm sized pixels and is placed 8\,cm from the sample. The simulation of the fluorescence patterns was performed with 4x4 oversampling and rebinning to the detector size. Only a single coherence mode was simulated. The number of photons emitted by the sample in 4$\pi$ was chosen as half the number of atoms in the sample. Hence, the mean number of photons per pixel is especially for small crystal sizes very small, still higher than the expected photon yield in an experiment.

The peak signal-to-noise ratio was calculated by simulating 2000 independent images and taking the mean intensity at the visible Bragg peaks positions as signal and the standard deviations at those positions over the independent simulations as noise, resulting in an estimated of the peak SNR of a single image.
Due to the low photon numbers, the Poisson noise is dominating the noise characteristic, and with an increase in atoms in the focus, the SNR increases linear (as shown in \fref{fig:SNRNatoms}). This holds up to the point where the Bragg peaks are no longer fully sampled and the signal decreases linear with a further increase in the number of atoms in the crystal, resulting in a nearly constant peak SNR.
If instead of considering the peak value of the reconstruction as the signal, the three-dimensional integral over the Bragg peak would be considered as the signal, the decreasing width of the Bragg peak with increasing size of the crystal would result in a nearly constant SNR under these low photon count conditions.
The peak signal-to-noise ratio was calculated by simulating 2000 independent images and taking the mean intensity at the visible Bragg peaks positions as signal and the standard deviations at those positions over the independent simulations as noise, resulting in an estimation of the peak SNR of a single image.
Due to the low photon numbers, the Poisson noise is dominating the noise characteristic, and with an increase in atoms in the focus, the SNR increases linearly (as shown in \fref{fig:SNRNatoms}). This holds up to the point where the Bragg peaks are no longer fully sampled and the signal decreases linearly with a further increase in the number of atoms in the crystal, resulting in a nearly constant peak SNR.
If instead of considering the peak value of the reconstruction as the signal, the three-dimensional integral over the Bragg peak was considered the signal, the decreasing width of the Bragg peak with increasing size of the crystal would result in a nearly constant SNR under these low photon count conditions.



Expand All @@ -306,7 +305,7 @@ \subsection{Parameters Influencing the Signal and Noise Characteristics}
\caption{SNR dependence on detector size}
\label{fig:SNRdetsize}
\end{subfigure}
\caption[SNR dependence on crystal size and detector size]{SNR dependence on the crystal size (a): As the crystal size increases, the SNR increases linearly (as shown by the linear regression) with the number of atoms in the sample as long as the Bragg peaks is sufficiently sampled (blue points), in undersampling conditions the signal strength decreases linear with the number of atoms resulting in a near constant SNR (right side of the dashed line). The error bars show the standard deviation of the SNR calculated for the 12 visible Bragg peaks. For details on the simulation parameters, see main text. The SNR dependence on the detector size for a simulated foil sample (b) shows a linear relationship between the square root of the number of pixels and the SNR.}
\caption[SNR dependence on crystal size and detector size]{SNR dependence on the crystal size (a): As the crystal size increases, the SNR increases linearly (as shown by the linear regression) with the number of atoms in the sample as long as the Bragg peaks is sufficiently sampled (blue points), in undersampling conditions the signal strength decreases linearly with the number of atoms resulting in a near constant SNR (right side of the dashed line). The error bars show the standard deviation of the SNR calculated for the 12 visible Bragg peaks. For details on the simulation parameters, see main text. The SNR dependence on the detector size for a simulated foil sample (b) shows a linear relationship between the square root of the number of pixels and the SNR.}

\end{figure}

Expand All @@ -324,9 +323,9 @@ \subsection{Parameters Influencing the Signal and Noise Characteristics}

A simulation with a varying number of spheres placed randomly inside the focal volume while ensuring a minimum distance between neighboring spheres with no other interaction was performed to decide if measuring the fluorescence of more than one nanoparticle in one shot is advantageous.
The minimum distance was chosen according to the typical size of nanoparticle capping agents; the number of spheres was increased up to a simulated Poisson sphere distribution (see appendix \fref{algo:bridson} for details) as densest random placement of particles, leading to 25\% of the volume filled by possible emitters.
Three factors determine the structure factor of these samples: The structure factor of the focus, the structure factor of points following the poison sphere distribution, and the structure factor of a single sphere (see \fref{fig:multisphere1} and \fref{fig:multisphere3}).
Three factors determine the structure factor of these samples: The structure factor of the focus, the structure factor of points following a Poisson sphere distribution, and the structure factor of a single sphere (see \fref{fig:multisphere1} and \fref{fig:multisphere3}).
For spheres with 20\, nm radius, a spacing layer of 5\,nm around each sphere, with 50000 excited atoms per sphere on average, a focus of 200\,nm FWHM, the fluorescence on a 1024x1024 pixel (pixel size 100\,\textmu m) detector placed 30\,cm was simulated. In this geometry, assuming constant distance to the sample for each detector pixel, approximately 1\% of the emitted photons reach the detector. For each number of spheres, 5000 images were used for a radial reconstruction using the direct method (see \fref{fig:multisphere2}).
As the number of photons emitted by each sphere was kept constant, opposing effects occur: With an increasing number of particles, more are photons are emitted, and the Poisson noise is reduced, but the signal strength in low scattering angles decreases as the structure factor changes from a single sphere to a hard-sphere model and the influence of the focal volume and the distribution of the spheres on the structure factor increases,
As the number of photons emitted by each sphere was kept constant, opposing effects occur: With an increasing number of particles, more photons are emitted, and the Poisson noise is reduced. On the other hand, the signal strength in low scattering angles decreases as the structure factor changes from a single sphere to a hard-sphere model and the influence of the focal volume and the distribution of the spheres on the structure factor increases.

In the case of 16 randomly positioned spheres, the mean distance to the nearest neighbor in the simulation is approximately three times the diameter. Compared to a single sample, the noise as measured by the standard error in the reconstruction is reduced by a factor of 10. At the same time, the signal strength at the first maximum decreases by a factor of only 3, leading to an overall gain in SNR. The dense random packing, with a mean nearest simulated neighbor distance of 1.4 times the diameter, reduces the noise further by an additional factor of 5, but reduces the signal by a factor of 10. Hence, the medium-dense sample shows the best SNR in this simulation.

Expand Down
2 changes: 1 addition & 1 deletion Tex/thanks.tex
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
\chapter*{Acknowledgements}
\onehalfspacing
I would like to thank my supervisors, Tais Gorkhover and Thomas Möller, for their continued support; the Gorkhover group at SLAC, i.e., Stephan Kuschel, Matt Ware, and Siqi Li, for the warm welcome to California, for allowing me to learn from them (not only about the life of Selachimorpha), for joined discussions, planning, preparation, and execution of the experiments. \textcolor{white}{Baby shark, do do, do do do do}
I would like to thank my supervisors, Tais Gorkhover and Thomas Möller, for their continued support; the Gorkhover group at SLAC, i.e., Stephan Kuschel, Matt Ware, and Siqi Li, for the warm welcome to California, for allowing me to learn from them (not only about the life of Selachimorpha\footnote{\textcolor{white}{Baby shark, do do, do do do do}}), for joined discussions, planning, preparation, and execution of the experiments.

Additionally, I would like to thank An Marshal for the help in getting the TEM images, Ben Reeves for his preparation of the GaAs samples, Carina Hedrich for the AAO sample preparation, Amy Cordones-Hahn for guidance in chemistry, and Cindy Patty for general awesomeness.

Expand Down
Loading

0 comments on commit 325cd37

Please sign in to comment.