API tools faq
paste
OreganoHauch

Bachelor thesis as of Jan 4, 2021

OreganoHauch
Jan 3rd, 2021
100
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 79.63 KB | None | 0 0
raw download clone embed print report
  1. \documentclass[10pt]{article}
  2. \usepackage{geometry}
  3. \usepackage[utf8]{inputenc}
  4. \geometry{a4paper, textheight = 720pt}
  5. %\usepackage[parfill]{parskip} % Activate to begin paragraphs with an empty line rather than an indent
  6. \usepackage{graphicx}
  7. \usepackage{subcaption} % This (and the package below) allows for having Fig a), Fig b) as subfigures
  8. \usepackage[labelformat=parens,labelsep=quad,skip=3pt]{caption}
  9. \renewcommand{\thesubfigure}{\thefigure\alph{subfigure})}
  10. \captionsetup[sub]{labelformat=simple}
  11. \usepackage{amssymb}
  12. \usepackage{amsmath}
  13. \usepackage{siunitx} % e.g. for the \mu
  14. \usepackage{comment} % for multiple lines comments
  15. \usepackage{url} % for referencing URLs (e.g. Zurich Instruments lock-in detection)
  16. \usepackage{hyperref} % for red borders around references
  17. \usepackage{tikz}
  18. \usepackage{epstopdf}
  19. \usepackage{titlesec}
  20. \newcommand{\sectionbreak}{\clearpage} % Create a page break for each section
  21. \usepackage[hang,flushmargin]{footmisc}
  22. \usepackage[page]{appendix}
  23. \renewcommand{\appendixpagename}{Appendix}
  24. \setcounter{secnumdepth}{4}
  25. \titleformat{\paragraph}
  26. {\normalfont\normalsize\bfseries}{\theparagraph}{1em}{}
  27. \titlespacing*{\paragraph}
  28. {0pt}{3.25ex plus 1ex minus .2ex}{1.5ex plus .2ex}
  29. \usepackage[backend=biber,
  30. sorting=none]{biblatex}
  31. \addbibresource{lit.bib}
  32. \DeclareGraphicsRule{.tif}{png}{.png}{`convert #1 `dirname #1`/`basename #1 .tif`.png}
  33. \numberwithin{equation}{section}
  34. \usepackage{setspace}
  35. \setstretch{1.3} % for custom spacing
  36.  
  37. \begin{document}
  38. \pagenumbering{arabic}
  39.  
  40. \begin{minipage}[t]{0.25\linewidth}
  41. \includegraphics[width=1.3\linewidth]{Figures/Logos/up-uhh-logo-u-2010-u-png_resized.png}
  42. \end{minipage}
  43. \hspace*{9cm}
  44. \begin{minipage}[t]{0.25\linewidth}
  45. \includegraphics[width=0.4\linewidth]{Figures/Logos/DESY_logo_3C_web_ger.png}
  46. \end{minipage}
  47.  
  48. \begin{center}
  49. \large
  50.  
  51. \textsc{University of Hamburg}\\
  52. \vspace*{2mm}
  53. \textsc{Faculty of Mathematics, Informatics and Natural Sciences}\\
  54. \vspace*{2mm}
  55. \textsc{Department of Physics}\\
  56.  
  57. \end{center}
  58.  
  59. \vspace*{1cm}
  60.  
  61. \begin{center}
  62. \Large\textsc{Bachelor Thesis}
  63. \end{center}
  64.  
  65.  
  66.  
  67. \hrulefill
  68. \vspace*{-.187cm}
  69.  
  70. \begin{center}
  71. \LARGE\textbf{Analysis of spatially resolved current-voltage dependence in CIGS solar cells}
  72. \end{center}
  73.  
  74. \vspace*{-.187cm}
  75. \hrulefill
  76.  
  77. \vspace*{2cm}
  78.  
  79. \begin{center}
  80. \end{center}
  81.  
  82. \begin{center}
  83. \Large submitted by:\\
  84. \vspace*{1cm}
  85. Jan Hense\\
  86. Matr. No: 6956738
  87. \end{center}
  88.  
  89. \vspace*{2cm}
  90.  
  91.  
  92. \begin{center}
  93. First Supervisor: Dr. Michael E. Stuckelberger\\
  94. Second Supervisor: Prof. Dr. Christian Schroer
  95. \end{center}
  96.  
  97. \vfill
  98.  
  99. \begin{center}
  100. Submitted at: January 8, 2021
  101. \end{center}
  102.  
  103. %\pagestyle{empty}
  104.  
  105. \restoregeometry
  106. %\title{Analysis of spatially resolved current-voltage dependence in CIGS solar cells\\ \small{Supervisor: Dr. Michael Stückelberger}}
  107. %\author{Jan Hense}
  108. %\date{\today} % Activate to display a given date or no date
  109.  
  110.  
  111. %\begin{document}
  112.  
  113. \pagenumbering{gobble}
  114. %\maketitle
  115.  
  116. \clearpage
  117. \begin{abstract}
  118. The IV curve of a solar cell carries information about its performance and about limiting performance factors. However, not all regions of a solar cell manifest the same performance due to localized electronic defects or inhomogeneities in the elemental distribution and thickness. Hitherto, IV curve analysis has been done in several ways to extract a solar cell's parameters under full or partially shaded illumination. Here, however, the solar cell is illuminated by a focused X-ray beam, which allows for spatially-resolved IV curve analysis.
  119. The examined solar cell has a copper indium gallium diselenide (CIGS) thin-film absorber material which over the years has become a competitive alternative to established materials such as crystalline silicon. Lock-in amplification filters out the voltage-induced current from a modulated signal, leaving only the X-ray beam induced current (XBIC). This way, the charge carrier behaviour of single regions can be examined under different bias voltages.
  120.  
  121. In this analysis, the performance of the solar cell is both resolved in its spatial and voltage dimension, yielding an IV curve for each scanned pixel. A first distinction into pixels with high and low current is made. It is found that pixels with high current for the most part retain a high current over different bias voltages and the same is found for the pixels with low current. A second distinction into pixels with good and bad diodes is made based on their behaviour at forward bias. A low correlation between these two distinctions leads to the hypothesis that the performance is mainly based on the sample's thickness, whereas the quality of the diodes is based on the charge collection efficiency.
  122.  
  123. \end{abstract}
  124. \tableofcontents
  125. \clearpage
  126. \pagenumbering{arabic}
  127.  
  128. \section{Introduction}
  129. \subsection{Motivation}
  130. Within the last decades, global warming and anthropogenic CO\textsubscript{2} emissions caught the attention of scientists, policy-makers and the media \cite{rich-2018}. Since then, the Intergovernmental Panel on Climate Change (IPCC) has published its Fifth Assessment Report in which it asserts with high confidence the anthropogenic CO\textsubscript{2} gas as the main cause for global warming.\cite{pachauri-2014-IPCC} In 2016, the world's energy usage took up 73\% of the total greenhouse gas emissions \cite{climatewatchdata-2016}, most of which stem from coal and gas.\cite{iea-2020} Heat and electricity mainly determine the energy usage, but renewable energy sources still play a comparatively small role. They have been on the rise in the global electricity mix, where their contribution has increased by a rate more than twice as their contribution to the global energy mix \cite{our-world-in-data-electricity}\cite{our-world-in-data-energy}, but this can be rather attributed to renewable energies substituting nuclear energy; fossil fuel consumption has barely been reduced over the last 20 years.
  131. Still, even slight improvements of solar cell conversion efficiencies from sunlight into electric energy can have an important impact on both economics and the environment. The type of material examined is copper indium gallium diselenide (CIGS) which has several interesting properties. Being a thin-film material and not being dependent on a glass substrate, it can be rolled during manufacturing and application. Of all thin-film materials, CIGS has the highest efficiency\cite[p.195]{smets-2016}, with a current record efficiency of over 23\%.\cite{green-2020-pip-january}
  132. Due to its thin material, it can be produced cheaply. The scarcity of indium however limits the advantages of CIGS.
  133.  
  134. \subsection{Outline}
  135. In this thesis, the spatially resolved X-ray beam induced current of a CIGS solar cell is examined under different bias voltages. This allows to differentiate between voltage-dependent and voltage-independent effects, determining factors for regions with high and low currents.
  136. In \textbf{Chapter 2}: \textit{Theory of solar cells}, the basic structure of a solar cell, as well as the specific structure of CIGS, and its behaviour under bias voltage and illumination are explained.
  137. In \textbf{Chapter 3}: \textit{Experiment}, the experimental setup environment at DESY's beamline and the design of the experiment are sketched. Also, lock-in amplification is shortly summarized, since it plays a big role in extracting the XBIC signal from the directly measured current.
  138. In \textbf{Chapter 4}: \textit{Data preparation}, it is explained how image registration and Gaussian filtering is applied to the data and which part of the data had to be discarded for proper IV curve analysis.
  139. In \textbf{Chapter 5}: \textit{Results and Discussion}, the IV curves are grouped in different ways which are then compared and discussed.
  140. In \textbf{Chapter 6}: \textit{Conclusions}, conclusions are drawn with respect to the type of analysis employed here and with regards to the results.
  141.  
  142. \section{Theory of solar cells}
  143. \subsection{The \textit{pn}-junction in equilibrium}
  144. The \textit{pn}-junction is the heart of a solar cell. It is crucial for generating electric current under illumination. A \textit{pn}-junction is composed of two different layers stacked together. Those layers differ in their \textit{doping} and can also differ in their materials being used.
  145. Doping means that additional \textit{donor} or \textit{acceptor} atoms are being added to the layer. In \textit{n}-type doping, donor atoms are added in a way that the layer's free electrons amount to \textit{n}, whereas in \textit{p}-type doping, acceptor atoms are added in a way that the layer's free holes amount to \textit{p}. These additional atoms are also referred to as \textit{impurity atoms}. A semiconductor with insignificantly many impurity atoms still has 'more' than the intrinsic number \textit{n\textsubscript{i}} of free holes/electrons. Electrons are called the \textit{majority carriers} in the \textit{n}-type layers, since they constitute the majority of free charge carriers there. In the \textit{p}-type layer, they are the \textit{minority carriers}, and the same goes for holes vice versa.
  146.  
  147. In the language of band theory, electrons are situated in the conduction band (CB), whereas holes are situated in the valence band (VB). It follows from the law of mass action and band theory\footnote{The derivation is shortly summarized in \cite[p.56]{smets-2016}.} that
  148. \begin{equation}
  149. \label{eqn:np}
  150. np = n_{\text{i}}^{2} = N_{\text{C}}N_{\text{V}}\cdot exp\left(\frac{E_{\text{V}}-E_{\text{C}}}{kT}\right),
  151. \end{equation}
  152. where N\textsubscript{C} and N\textsubscript{V} is the number of states in the CB and VB respectively. E\textsubscript{C}$-$E\textsubscript{V} is the band gap energy and $kT$ the thermal energy. Interesting to note about this equation is that the product $np$ is independent of the Fermi energy. Since doping shifts the Fermi energy level, the equation will still be valid for semiconductors that are doped. For example, if we dope 1 cm\textsuperscript{3} of Si with 10\textsuperscript{17} P atoms, we can calculate from Equation \ref{eqn:np} that there will be about $2.22\cdot 10^3$ holes. This huge difference underlines the distinction made between majority and minority carriers.
  153. Both, free holes and free electrons, can move inside the layer and also enter the adjacent layer. The movement, however, is restricted by the \textit{diffusion length} and \textit{lifetime}. The diffusion length is defined for minority carriers only.\footnote{The majority charge carriers are much less likely to recombine, so their diffusion length would exceed by far the dimension of the \textit{pn}-junction. Therefore, minority charge carriers determine the current.} If the layer materials only differ by the elements injected during doping, then the junction is called a \textit{homojunction}. A typical example would be a crystalline silicon solar cell (c-Si).\footnote{Amorphous silicon solar cells (a-Si) have a p-i-n junction, where there is an intermediary intrinsic layer which is not doped, because the absence of an ordered crystalline lattice structure introduces dangling bonds that limit the diffusion length of charge carriers.\cite[p.181]{smets-2016}} If the layers differ by the materials being used, or more precisely their band gaps differ, then the junction is called a \textit{heterojunction}. As we will see later, the layers in CIGS solar cells are composed of different materials.
  154. \\
  155. \\
  156. Now, let us consider an \textit{n}-type and \textit{p}-type layer forming a \textit{pn}-junction. For clarity of reading, the processes are explained here and in the following subsection for the flow of electrons only, noting that the same process\footnote{By principle only, since the level of \textit{n}-type doping is often different from the level of \textit{p}-type doping.} is happening with the holes, only vice versa.
  157.  
  158. A free electron leaves a positively ionized donor atom. Both charges are equally opposite to each other and therefore the bulk of both the n- and \textit{p}-type region remains electrically neutral during doping. However, the huge difference in concentration of electrons allows for a significant diffusion across the \textit{pn}-junction, where the free electrons recombine with the holes that are abundant in the \textit{p}-type region. Because of this recombination, this flow of electrons is called the \textit{recombination current} $I\textsubscript{rec}$. The electrons close to the depletion region will probably be the first to recombine right at the opposite side of the junction. The recombination will create a negatively ionized \textit{space charge} in the \textit{p}-type region. "Space charge" refers to the fact that, in contrast to the majority of ionized acceptor atoms in the \textit{p}-type region, this charge will remain fixed at its position and remain electrically charged, for there aren't enough holes in its vicinity that the electron could again recombine with.
  159. This process will repeat until the space charges have built up an electric field that is exactly opposite to the recombination current. Now, that there is a region depleted of free charge carriers, we refer to it as the \textit{depletion region} and to the electric field as the \textit{internal electric field}. Due to thermally generated free electrons in the \textit{n}-type region\footnote{Also in the depletion region, but that is rather irrelevant, because the size of the depletion region is very small compared to the bulk regions.}, that will randomly move into the depletion region where they experience the force of the internal electric field, there must be a separate flow of electrons be accounted for. It is called the \textit{drift current} or \textit{thermal-generation current} $I\textsubscript{gen}$.\footnote{The electrons of the recombination current will experience the force of the internal electric field of course as well, but this decrease in current is already included in the $I\textsubscript{rec}$ variable.} The \textit{pn}-junction is now in equilibrium, when there is no net current $I$ flowing through the depletion region:
  160. \begin{equation}
  161. \label{eqn:currents-in-equilibrium}
  162. I = I_{\text{rec}} - I_{\text{gen}} = 0
  163. \end{equation}
  164.  
  165. \begin{figure}[!ht]
  166. \centering
  167. \includegraphics[width=0.6\textwidth]{Figures/Fig8.2_p-n_junction_in_equilibrium.png}
  168. \caption{Schematic of the \textit{pn}-junction in equilibrium with an internal electric field. Two import remarks must be made about the illustration: First, the minus in the \textit{p}-type region symbolizes local charges, but it is not fixed in space; unlike the minus at the space charge region that symbolize a fixed charge. Second, the space charge region is not exactly sharp, but has an exponential fall-off which can be approximated to a sharp edge, neglecting however certain effects mentioned in subsection \ref{ssec:important-parameters}. Image source: \cite[p.85]{smets-2016}}
  169. \label{fig:pn-junction-equilibrium}
  170. \end{figure}
  171. \subsection{The \textit{pn}-junction under bias voltage and illumination}
  172. \label{ssec:pn-junction-bias-voltage}
  173. A solar cell with a \textit{pn}-junction as displayed in Figure \ref{fig:pn-junction-equilibrium} is, of course, not used for producing power, but the solar cell will be part of an array of solar cells in both series and parallel circuits and experience external potentials from the adjacent solar cells and the external circuit. More on that at the end of the subsection.
  174.  
  175. If an external \textit{bias voltage} is applied to the front and back contact of the solar cell, the internal electric field will become smaller or larger, depending on whether the potential is positive (\textit{forward bias}) or negative (\textit{reverse bias}) on the p-doped side with respect to the n-doped side. The generation current will not be affected by the voltage, as long as it is in a moderate range, whereas the recombination current increases with the Boltzmann factor:
  176. \begin{subequations}
  177. \begin{align}
  178. \label{eqn:dark-IV-I_gen}
  179. I_{\text{gen}}(V) &= I_{\text{gen}}(V=0) = I_{\text{0}}\\
  180. \label{eqn:dark-IV-I_rec}
  181. I_{\text{rec}}(V) &= I_{\text{rec}}(V=0)\cdot exp\left(\frac{eV}{kT}\right)
  182. \end{align}
  183. \end{subequations}
  184. $I\textsubscript{0}$ is called the \textit{saturation current} and since it is a measure of the amount of recombination in a \textit{pn}-junction, a good diode should have it as low as possible.
  185.  
  186. Under forward bias condition, a negative electrode will be attached to the contact of the \textit{n}-type region. This means that electrons flow through the cathode into the \textit{n}-type layer. Since there are now more electrons than before, the diffusion across the depletion region will become larger and the recombination current $I\textsubscript{rec}$ increases. The external electric field adds to the internal electric field, such that the electric field across the depletion region is decreased. Due to the increased diffusion current, many more electrons will be present in the \textit{p}-type region; the amount of ionized dopants will become larger than the number of majority charge carriers. Electrons are now the majority charge carriers in the \textit{p}-type region and the negative space charges are excited to a level that breaks the bond of ionized dopants. Therefore, the depletion region will shrink. Since the generation current only depends on charge carriers reaching the depletion region within the minority charge carrier diffusion length and shrinkage of the depletion region is minor compared to the dimension of the \textit{pn}-junction, $I\textsubscript{gen}$ remains constant. A current flows despite the circumstance that there are so many ionized dopants which could act as a countering force, because the anode will constantly supply holes for the electrons to recombine with.
  187. The forward bias will move the Fermi energy for the \textit{n}-type region $E_{\text{FC}}$ closer to the CB edge.
  188.  
  189. Under reverse bias condition, electrons are injected into the \textit{p}-type region, where holes are the majority carriers, and recombine with the holes. Electrons experience an electric force away from the cathode and a from the negatively charged edge of the depletion region which poses as a potential barrier for electrons. This way, the electrons will be more likely to recombine at the edge of the depletion region and thereby widen it. The larger the reverse bias voltage, the less diffusion will take place. Since diffusion leads to restoration of equilibrium by recombination, it is detrimental for the overall flow of current. A reverse bias voltage will lead to less recombination by diffusion and the recombination current $I\textsubscript{rec}$ will be reduced. A larger electric field will not affect the generation current. For that reason, for high reverse bias voltages, there is a small remaining saturation current.
  190. The reverse bias will move the Fermi energy for the \textit{p}-type region band $E_{\text{FV}}$ closer to the VB edge.
  191.  
  192. Under illumination, additional energy is given to the thermal excitation energy and more electron-hole pairs are created. This significantly increases the number of minority charge carriers and adds to the generation current $I\textsubscript{gen}$. Putting all of these contributions together finally yields the \textit{diode} or \textit{Shockley equation}:
  193. \begin{equation}
  194. \label{eqn:shockley}
  195. I = I_{\text{ph}} + I_{\text{gen}} - I_{\text{rec}} = I_{\text{ph}} + I_{\text{0}}\cdot \left[1 - exp\left(\frac{eV}{kT}\right)\right]
  196. \end{equation}
  197. The equation is visualized in Figure \ref{fig:I-V-characteristics}. There it can be seen that the photo-generated current corresponds to an upward shift of the dark IV curve, because in this model it is voltage-independent and has the opposite direction to the voltage-induced current.\\
  198. Later, when we analyze the current-voltage dependence in CIGS solar cells, it will be seen that the model does not hold and voltage-dependent effects of the (X-ray beam induced) photo-current become visible. Note that the direction of the photo-current is determined by the internal electric field only, since there the electron-hole pairs are separated. This will later help to understand Figure \subref{fig:violin_currents-direct-zeroed_numbered} (subsection \ref{ssec:current-distribution}).\\
  199. It depends on the energy of the photon, how much an electron will be excited. It can be excited beyond the CB edge, but then it will dissipate its energy through \textit{thermalization}. A currently researched future technology are hot-carrier solar cells (HCSC) which seek to slow down thermalization and make use of that additional thermalization energy. The maximum energy of such an electron could then exceeds $e \cdot V_{\text{OC}}$, but in normal (non-hot carrier) solar cells this is a negligible rare process.
  200. \\
  201. How does all of this pertain to the real-world application of a solar cell? As briefly mentioned above, solar cells are merged into a series and parallel circuit array. Where the solar cells are connected in series, the voltages add up. Variation in the load will change the voltage of the individual cell, between 0 (for zero resistance)\footnote{Diffusion would no longer happen across the depletion region, but through the conducting wire and recombination would not take place in the other doped region, but in the wire. If this was not the case, the recombination in the other doped region would leave even more ionized dopants and thereby increase the voltage drastically.} and the voltage that is created by the illumination under open-circuit condition. But this is only a partial contribution to the overall voltage that an individual solar cell will experience inside the array: The adjacent solar cells will also produce an additional voltage under illumination which is contrary to the internal electric field voltage, because the photo-generated electron-hole pairs will elicit a counteracting recombination current that, just like during forward-biasing, will lower the potential difference across the depletion region.
  202. The series circuit will connect the p-doped region with the n-doped region. This comes akin to reverse-biasing the solar cell. So, under ideal working circumstances, the array's solar cells will be reverse-biased. However, under non-ideal circumstances, there can be two detrimental effects of reverse-bias induced by adjacent solar cells. One is \textit{partial shading}: If one of the cell does not receive any light, it will not create any power, but because its p-side is attached to an n-side and its n-side to an p-side, it is still forced to have a voltage. This will inevitably lead to heat dissipation\cite[p.256]{smets-2016}. Similarly, reverse currents occur due to ohmic shunt resistances which is the most frequent source of detrimental reverse currents\cite{breitenstein-2011}.
  203.  
  204. \subsection{Important parameters for characterization of solar cells}
  205. \label{ssec:important-parameters}
  206. When the p- and \textit{n}-type region are short-circuited, the voltage will be zero and the current will be determined by the photo-generated electron-hole pairs. As $I$ approaches $I\textsubscript{SC}$, it is the maximum current of a solar cell that still produces power.
  207. In contrast, when the p- and \textit{n}-type region are not connected or the resistance is infinitely high, then the voltage will be the \textit{open-circuit voltage} $V\textsubscript{OC}$. The $V\textsubscript{OC}$ could also be defined as the voltage established by the internal electric field or the bias voltage that will be necessary to make $I$ zero. As $V$ approaches $V\textsubscript{OC}$, it is the maximum voltage of a solar cell that still produces power. It can be calculated from rearranging Equation \ref{eqn:shockley}:
  208. \begin{subequations}
  209. \begin{align}
  210. \label{eqn:ISC}
  211. I_{\text{SC}} &= I_{\text{ph}}\\
  212. \label{eqn:VOC}
  213. V_{\text{OC}} &= \frac{kT}{e}\cdot ln\left(\frac{I_{\text{ph}}}{I_{\text{0}}}+1 \right)
  214. \end{align}
  215. \end{subequations}
  216. As can be seen from Equation \ref{eqn:VOC}, the $V\textsubscript{OC}$ mainly depends on the saturation current, because it varies much more over different diodes (by orders of magnitude) than the photo-generated current changes over different weather conditions. $I\textsubscript{0}$, in turn, depends on the recombination and therefore $V\textsubscript{OC}$ is highly sensitive to the amount of recombination in the solar cell. It is one of the most important parameters of a solar cell, since it carries information of how much power eventually can be produced and how much will be lost to recombination.
  217. Only when $V\in (0, V\textsubscript{OC})$ will the solar cell create power. Otherwise, it will consume power. The point at which the solar cell produces the maximum amount of power is called the \textit{Maximum Power Point} $(I\textsubscript{MP}, V\textsubscript{MP})$.
  218. The variables introduced so far are contained in the \textit{Fill Factor}:
  219. \begin{equation}
  220. \label{eqn:FF}
  221. FF = \frac{I_{\text{MP}}\cdot V_{\text{MP}}}{I_{\text{SC}}\cdot V_{\text{OC}}}
  222. \end{equation}
  223. The fill factor is a measure of how efficiently the solar cell produces power, compared to the power $P\textsubscript{in}$ of the irradiance it receives. The \textit{conversion efficiency} $\eta$ of the solar cell is then measured under standard (ideal) test conditions, yielding $\eta\textsubscript{STC}$ and under normal operating (realistic) cell temperature conditions, yielding $\eta\textsubscript{NOCT}$\footnote{Under NOCT conditions, the cell is tested as if there was no load, so the voltage will be the $V\textsubscript{OC}$ and since energy is not dissipated in the load, this will cause more heating and therefore a higher temperature than under standard test conditions.}:
  224. \begin{equation}
  225. \label{eqn:conversion-efficiency}
  226. \eta = \frac{I_{\text{MP}}\cdot V_{\text{MP}}}{P_{\text{in}}} = \frac{I_{\text{SC}}\cdot V_{\text{OC}}\cdot FF}{P_{\text{in}}}
  227. \end{equation}
  228. The V\textsubscript{OC} could be thought of as proportional to $E_{\text{g}}$, since this is the maximum electric energy an electron can have in the semiconductor, but because of the Urbach energy\cite{solomon-2015}, $e\cdot V_{\text{OC}}$ will be lowered by a subtraction term. This has to do with the fact that the space charge region is not exactly sharp which causes certain effects under strong illumination.\cite{crandall-1996} I\textsubscript{SC} decreases with the band gap energy.\footnote{Why the V\textsubscript{OC} increases with E\textsubscript{g}, whereas I\textsubscript{SC} decreases, can also be explained this way: With higher E\textsubscript{g} more energy is needed to excite an electron-hole pair, so this will limit the current, but the V\textsubscript{OC} is measured under forward bias condition, where the external electric field diminishes the potential barrier across the junction more than the higher E\textsubscript{g} increases it.} Thus, in CIGSS solar cells, V\textsubscript{OC} and E\textsubscript{g} have an affine linear relationship. In CIGS solar cells, for reasons that are not understood yet\footnote{Models suggest that it has to do with an unfavourable band alignment at the n-CIGS and p-CdS interface. See figure \subref{fig:CIGS-layer-structure}. Models cited in: Weinhardt et. al (Ibid.)}, this cannot be said as easily\cite{weinhardt-2017}.
  229.  
  230. \begin{figure}[!ht]
  231. \centering
  232. \includegraphics[width=0.6\textwidth]{Figures/Fig8.9_I-V characteristics_mirrored.jpg}
  233. \caption{IV characteristic. The lower black line shows the current due to thermal-generation and diffusion. It approaches $I\textsubscript{0}$ under reverse bias voltage, as can only be seen barely, and otherwise shows an exponential increase. The upper black line is the illuminated IV curve, i.e. the dark IV curve shifted by the photo-generation current (blue line). The product of I and V yields the power output, which has a maximum at the Maximum Power Point ($J\textsubscript{MP}$,$V\textsubscript{MP}$). Edited from this image source: \cite[p.99]{smets-2016}}
  234. \label{fig:I-V-characteristics}
  235. \end{figure}
  236.  
  237. If one takes into account the non-ideality of the solar cell, new variables, such as for shunt and ohmic (parallel) resistance are introduced, and there are different models as well as different methods to account for this and to determine these variables\cite{cotfas-2013}.
  238.  
  239. \subsection{CIGS solar cells}
  240. \label{ssec:cigs-solar-cells}
  241. Copper indium gallium diselenide (CIGS) belongs to the second generation of solar cells. In contrast to older materials, it is cheaper to produce and can be deposited on flexible materials which allows to build thinner solar cells.\footnote{For c-Si, all of the light is absorbed within a few micrometers, but the absorber is about \SI{180}{\micro\meter} thick for mechanical support\cite{ramanujam-2017}}. Thereby, the energy payback time ($\sim$ 1 year) is less than that of c-Si ($\sim$ 2 years) now. CIGS was an improvement of the \textit{copper indium diselenide} (CIS: CuInSe\textsubscript{2}) material which was first synthesized in 1953, by Harry Hahn et. al. It received major attention in the 80s, when efficiencies of over 9\% were reported. Gallium and Sulfide (for CIGSS) were added then\cite{ramanujam-2017}.
  242. \\
  243. Chemically, CIGS is a solid solution of CIS and \textit{copper gallium diselenide} (CGS: CuGaSe\textsubscript{2}) with the chemical formula CuIn\textsubscript{x}Ga\textsubscript{1-x}Se\textsubscript{2}, where $x\in [0,1]$ is a measure of the indium/gallium ratio. Varying x from 0 to 1 allowed Xin Cui et. al to \textit{tune} the band gap from 1.08 to 1.53 eV\footnote{If sulfur is present as well, the Se:S ratio y can be tuned as well, allowing for a band gap range of 1.0 to 1.7 eV\cite[p.192]{smets-2016}.}.\cite{cui-2015} Tuning the band gap with gallium allows to better match the solar spectrum ($\lambda_{\text{red}}\approx$ 1.77 eV)\cite{ramanujam-2017} and a higher gallium to indium ratio coincides with a lower charge collection efficiency, except for cells that might have a larger amount of defects and inhomogeneities\cite{pyrlik-2020-bachelor}. CIGS has a chalcopyrite structure which means that its tetrahedrally bonded. It has intrinsic copper defects which is part of the reason why it is intrinsically p-doped.\footnote{The physical description of this and other defects in CIGS solar cells is still pending more ongoing research.} CIGS is a direct band gap material and therefore has a higher absorption coefficient which means that an excited charge carrier will recombine earlier. This enables a thinner material and allows to classify CIGS as a \textit{thin film technology}.
  244. \begin{figure}
  245. \centering
  246. \begin{subfigure}{0.4\linewidth}
  247. \centering\includegraphics[width=0.7\linewidth]{Figures/Fig13.22a_layer_structure.png}
  248. \caption{Layer structure of CIGS solar cell. The \textit{pn}-junction is also referred to as a \textit{buried junction}, because the separation of electron-hole pairs happens one layer above. Image source: \cite[p.193]{smets-2016}}
  249. \label{fig:CIGS-layer-structure}
  250. \end{subfigure}
  251. \hspace{1cm}
  252. \centering
  253. \begin{subfigure}{0.4\linewidth}
  254. \centering\includegraphics[width=0.7\linewidth]{Figures/Fig13.22b_band_diagram.png}
  255. \caption{Band diagram of CIGS solar cell. Image source: \cite[p.193]{smets-2016}}
  256. \label{fig:CIGS-band-diagram}
  257. \end{subfigure}
  258. \end{figure}
  259. Figure \subref{fig:CIGS-layer-structure} shows the different layers that a CIGS solar cell is composed of. The substrate is glass in this case, and then a Molybdenum layer of about \SI{500}{nm} thickness serves as an electric back contact. Molybdenum is chosen over other electrically conductive materials, not only due to its high conductivity, but also because it is chemically and mechanically stable during the CIGS growth (\textit{selenization})\cite{ong-2018}. Next, a \textit{p}-type absorber layer of about \SI{2}{\micro\meter} and a thin \textit{n}-type layer (for example with a high Indium/Gallium ratio) are deposited. The \textit{n}-type CIGS prevents the absorption of holes. Next, another thin layer, consisting of Cadmium Sulfide, with around \SI{50}{nm} thickness, is deposited. The CdS serves as a buffer layer to shield the absorber surface during the further layer deposition process and also to smoothen the band offset across the heterojunction (see Figure \subref{fig:CIGS-band-diagram}).\footnote{This CdS layer has been much thicker in the past, before it was replaced with the subsequent i-ZnO and Al-ZnO layers\cite{ramanujam-2017}}\cite{rana-2017}. This aims at better utilizing the solar spectrum, because the CdS acts as a \textit{window layer} that transmits low-energy photons and absorbs high-energy photons. A non-optimal CdS can cause severe reductions to the V\textsubscript{OC}\cite{weinhardt-2017}. \textit{Shockley Read Hall} (SRH) recombination becomes dominant at the n-CIGS/CdS interface, whereas for the \textit{n}-type CIGS it is dominant in the bulk region. An intrinsic zinc oxide (i-ZnO) layer is then followed by an aluminium-doped layer (ZnO:Al) which is effectively n-doped. The i-ZnO is highly resistive to prevent aluminium from the ZnO:Al layer to diffuse into the CIGS absorber. layer\cite{ramanujam-2017}. The ZnO:Al is transparent to light and acts as the \textit{Transparent Conductive Oxide} (TCO). As seen in Figure \subref{fig:CIGS-band-diagram}, with its band gap energy of 3.3 eV, no visible photon will be absorbed there.
  260. The photo-generated electron-hole pairs are separated at the CdS-CIGS interface, to be collected at the Al-ZnO.
  261. \section{Experiment}
  262. \subsection{Experimental environment}
  263. This work analyzes data which have been acquired amongst several other measurements that were led by Michael Stuckelberger in November 2018 at Deutsches Elektronen-Synchrotron (DESY) in Hamburg, Germany. There, special interest is taken in the nanoscale structure of solar cell samples. The microprobe station at the P06 beamline at DESY's synchrotron PETRA III offers such a unique opportunity for examining X-ray beam induced current (XBIC)\footnote{XBIC was first applied to solar cells in 2000 by Hieslmaier et. al.\cite{hieslmair-2000-csscmp}(Cited in: \cite{stuckelberger-2017-jmr})}. It is positioned in the \textit{microhutch} where the experiments are conducted. It receives a hard X-ray beam from the optics hutch, about 40 to 60 meters from the microhutch, with an energy that can be continuously tuned between 6 and 18 keV when using the channel-cut crystal monochromator (CCM). The prefocusing system of the optics hutch is based on beryllium compound refractive lenses (Be CRL) and is used to further adapt the photon flux and coherence properties for the micro- and nanohutch\cite{schropp-2020-jac}. By using corrective phase plates, the FWHM lies around \SI{100}{nm} and diffraction-limited X-rays are achieved\cite{seiboth-2017-ncomm}. This setup is often used when the microhutch is operated in conjunction with the nanohutch (for Ptychography\footnote{In the nanohutch, resolving samples below \SI{10}{nm} is the current goal, and to focus a \SI{50}{nm} FWHM X-ray beam necessitates advanced X-ray optics which cumulated in the ptychographic nano-analytical microscope (PtyNAMi).}).
  264. The microhutch's microscope focuses the X-ray beam with CRLs onto the sample that is fixed at the stage. The stage then is able to move steps with a size of \SI{50}{nm}. For short dwell times, the stage can also move continuously in a so-called \textit{fly-scan}, as was done for this experiment.
  265. Furthermore, the P06 employs many different techniques at the same time to get a much better imaging representation of samples with complex composition such as CIGS solar cells, that have multiple layers (compare subsection \ref{ssec:cigs-solar-cells}). This \textit{multi-modal} approach is a unique possibility that the P06 can provide and has successfully combined so far XBIC with X-ray fluorescence (XRF)\cite{west-2017-jsr}\cite{west-2017-nanoen}, X-ray excited optical luminescence (XEOL), X-ray diffraction (XRD), X-ray beam induced voltage (XBIV)\cite{stuckelberger-2017-pvsc-xbiv} and ptychography \cite{stuckelberger-2018-mimi}. This vast combination of different techniques under hard X-ray illumination allows for point-by-point correlations, however still limited by the necessity of thorough 3D simulations to complement the measurements, increased perturbations due to X-ray beam induced temperature, bias voltage, bias light and locally generated electric fields and, most importantly, the sheer amount of data that needs to processed and still has pending reconstruction algorithms, possibly being alleviated by machine learning in the future\cite{stuckelberger-2017-jmr}.
  266. The hard X-ray beam allows for deeper penetration depth, revealing more of the nanoscale ingredients.
  267.  
  268. \subsection{Experimental setup (Nov 2018)}
  269. \label{ssec:experimental-setup}
  270. For this experiment, the X-ray beam was focused to \SI{105}{nm} $\times$ \SI{108}{nm} (vertical $\times$ horizontal) at FWHM with an energy per photon of 15.25 keV. Its direction defines the \textit{downstream} X-axis, and the axis parallel/perpendicular to the ground defines then the Y-/Z- axis respectively. The sample was mounted on a stage that was moved, amongst others, by a PI piezo motor. The movement is a so-called \textit{snake scan}, where the stage moves 200 steps, each \SI{50}{nm}, in the Y direction, then simultaneously 1 step in the Y and Z direction, and then reverses the Y movement. One step has a \textit{dwell time} of 9.7 ms on average.\footnote{The scans were not optimized for simultaneous XRF measurements which require about ten times as many photons for good signal/noise ratio. As the chopper halves the amount of photons, it is moved out of the beam for XRF-optimized scans.} A scan during one step is here referred to as a \textit{pixel}. Thereby, a 201 $\times$ 200 grid can be extracted.
  271. \begin{figure}[!ht]
  272. \centering
  273. \includegraphics[width=\textwidth]{Figures/Cablingoverview29012019hellblau.png}
  274. \caption{Setup scheme of the experiment with focus on the signal paths. The current is converted into a voltage and pre-amplified. Then it is split up into two signal paths. One path is lock-in amplified, the other is not. Both paths then lead into the voltage-to-frequency converters (V2F) and are separated into positive and negative signals where the data acquirement unit (DAQ) picks them up. Image source: \cite[p.9]{ossig-2019-jove}}
  275. \label{fig:setup-of-experiment}
  276. \end{figure}
  277. Figure \ref{fig:setup-of-experiment} shows the setup of the experiment. The beam is first modulated by a chopper with 8.015 kHz and then focused with KB mirrors onto the device under test (DUT), i.e. CIGS sample. The sample was provided by Empa. The front contact (\textit{n}-type, compare subsection \ref{ssec:cigs-solar-cells}) was grounded, because the X-ray photons have enough energy to eject charge carriers, close to surface, from the side facing upstream. These are of interest for total electron yield (TEY) measurements, but not for XBIC, since a replacement current will be caused, mainly by the majority charge carriers there that then find minority charge carriers to recombine with\cite{stuckelberger-2017-jmr}. The voltage to frequency converters (V2f) do not accept negative voltages, therefore an inverter is used to also capture negative signals. The frequency can be calculated back to the solar cell current, as described by Ossig et. al (compare section \ref{sec:data-preparation})\cite{ossig-2019-jove}.
  278. \subsubsection{Pre-amplification and Lock-in amplification}
  279. \label{sssec:pa-and-lia}
  280. Usually, XBIC is measured under short-circuit condition. However, a voltage series was conducted here; alternating from negative to positive bias voltage to avoid hysteresis. The bias voltage was supplied by the Stanford Research Systems SR570 current-sensitive pre-amplifier which converts the current into a voltage, with an amplification factor of 10\textsuperscript{6} V/A. Its signal was both dispensed to the lock-in amplifier as well as directly to the DAQ. In the latter case, the signal is referred to as the \textit{direct signal}.
  281. \\
  282. Lock-in amplification is a technique that especially plays to its strengths when an original AC signal needs to measured which is smaller than the background in the input signal. It improves the signal/noise ratio by orders of magnitude, decreases the sensitivity to measurement artifacts, enables applying bias voltage and measuring in environments with light noise or intentionally applied bias light\cite{stuckelberger-2018-mimi}.
  283. For this high-precision experiment, a UHFLI lock-in amplifier (LIA) by Zurich Instruments that is certified to only have 4 $nV/\sqrt{Hz}$ noise was used with an amplification factor 1 V/V. Figure \ref{fig:lia-scheme} sketches its working principle. Detailed explanations are found in \cite{zurich-instruments-2016} and \cite{meade-1983}. The steps described in the following can also be achieved by analog means, however as this adds noise, they are performed digitally in high-precision instruments such as the UHFLI LIA. The sample needed to be small, because increasing the size would increase the voltage-induced (DC) current which is regarded as noise by the LIA (compare \cite{stuckelberger-2017-pvsc-xbiv}). First, the input (AC) signal $V_{\text{s}}(t)$ is split. The reference signal $V_{\text{r}}(t)$ is a square wave modulation supplied by the optical chopper with frequency $\omega_{\text{r}}$. Since there can be a delay between the phase of the carrier and the reference signal, V\textsubscript{r} is also shifted by $\frac{\pi}{2}$ which later reveals the phase difference $\phi$. For demodulation, both signals are then multiplied by the input signal and averaged over a certain time span that can be specified by the number of modulation cycles/time. Here, the modulation cycle is represented b The calculation is simplified when representing V\textsubscript{s} and V\textsubscript{r} as complex numbers\footnote{If one does not prefer to calculate with complex numbers, the in-phase and quadrature components can also be calculated by separately multiplying with $sin(\omega_{\text{r}}t)$ and $cos(\omega_{\text{r}}t)$. Terms will then cancel out during averaging due to the orthogonality of sine and cosine. $R$ is the root-mean-square of the amplitude, hence the prefactor $\sqrt{2}R$.}:
  284. \begin{align}
  285. \label{eqn:lock-in-multiplication}
  286. \left[V_{\text{s}}\right]\cdot \left[V_{\text{r}}\right] &= \left[\sqrt{2}R\cdot cos(\omega_{\text{s}}t + \phi)\right]\cdot \left[\sqrt{2}cos(\omega_{\text{r}}t) - i\sqrt{2}sin(\omega_{\text{r}}t)\right] \nonumber \\
  287. &= \left[\frac{R}{\sqrt{2}}e^{+i(\omega_{\text{s}}t+\phi)} + \frac{R}{\sqrt{2}}e^{-i(\omega_{\text{s}}t+\phi)}\right]\cdot \left[\sqrt{2}e^{-i\omega_{\text{r}}t}\right] \nonumber \\
  288. &= R\left[e^{i[(\omega_{\text{s}}-\omega_{\text{r}})t+\phi]}+e^{-i[(\omega_{\text{s}}+\omega_{\text{r}})t+\phi]}\right]
  289. \end{align}
  290. Since the signal frequency $\omega_{\text{s}}$ is modulated by the chopper frequency $\omega_{\text{r}}$, $\omega_{\text{r}} \approx \omega_{\text{s}}$. Therefore, the first term in equation \ref{eqn:lock-in-multiplication} will become time-independent and since the second term is harmonic, it can be time-averaged to zero. The product of V\textsubscript{s} and V\textsubscript{r} then simplifies to
  291. \begin{equation}
  292. \label{eqn:lock-in-final-product}
  293. V_{\text{s}}\cdot V_{\text{r}} = R\cdot e^{i\phi}.
  294. \end{equation}
  295. The real part of equation \ref{eqn:lock-in-final-product} is called the \textit{in-phase component} and often denoted by $X$, whereas the imaginary part is called the \textit{quadrature component Y}. R then represents the root-mean-square amplitude of the original signal for the modulation cycle (i.e. time interval over which has been averaged). The time-averaging electronically corresponds to connecting an RC low-pass filter of n-th order that specifies the roll-off behaviour. Averaging over longer periods will yield a trade-off between less noise and a slower response of R. The modulation frequency should be set below the frequency of the noise, as the low-pass filter cuts off frequencies that are higher than the reference signal. Lastly, it should be remarked that LIA itself is prone to noise. For example, all electrical instruments in the experiments should be grounded to the same point to avoid ground loops.
  296.  
  297. \begin{figure}[!ht]
  298. \centering
  299. \includegraphics[width=0.6\textwidth]{Figures/lockin_primer_fig_02_dual_phase_corrected.png}
  300. \caption{The modulated input signal V\textsubscript{s} is (digitally) split and multiplied with the reference signal V\textsubscript{r} and a $\frac{\pi}{2}$ phase-shifted reference signal respectively. The result is fed to a low-pass filter with a specified cut-off frequency at which the signal power ($\sim$ amplitude\textsuperscript{2}) is diminished by half. Edited from this image source: \cite[p.2]{zurich-instruments-2016}}
  301. \label{fig:lia-scheme}
  302. \end{figure}
  303.  
  304. \section{Data preparation}
  305. \label{sec:data-preparation}
  306. \subsection{Raw data}
  307. The data gathered from the DAQ has the HDF5 file format. It contains several dimensions of data, of which this analysis uses:
  308. \begin{itemize}
  309. \item Y position: Fast scan, encoded.
  310. \item Z position: Slow scan, encoded.
  311. \item Frequency counts from the V2f converter / amplitude R: $n_{\text{XBIC}}$.
  312. \item Frequency counts from the V2f converter / Negative direct signal.
  313. \item Frequency counts from the V2f converter / Positive direct signal.
  314. \item QBPM counts / normalized photon flux: $n_{\text{QBPM}}$.
  315. \item Dwell time / time of exposure during which the detectors integrate the signal (according to a continuous fly-scan, as described in subsection \ref{ssec:experimental-setup}): $T$.
  316. \end{itemize}
  317.  
  318. The DAQ data can be converted to a current by following the formula described by Christina Ossig \cite[p.4]{ossig-2019-jove}:
  319. \begin{subequations}
  320. \begin{align}
  321. \label{eqn:XBIC-conversion}
  322. I_{\text{XBIC}} &= \alpha_{\text{norm}}\cdot H_{\text{conv}}\cdot n_{\text{XBIC}}\\[10pt]
  323. \label{eqn:XBIC-conversion_normalization}
  324. \alpha_{\text{norm}} &= \frac{1}{T}\cdot \frac{\overline{n_{\text{QBPM}}}}{(n_{\text{QBPM}}/T)}\\[10pt] \label{eqn:XBIC-conversion_linear-factor}
  325. H_{\text{conv}} &= \frac{2\cdot W_{\text{ff}}}{C_{\text{V2f}}\cdot A_{\text{PA}}\cdot A_{\text{LIA}}}
  326. \end{align}
  327. \end{subequations}
  328.  
  329. where $\alpha_{\text{norm}}$ is the normalization factor for each pixel to account for a varying photon flux and $H_{\text{conv}}$ is a linear factor to take into account the amplification and conversion settings of the pre-amplifier, lock-in amplifier and V2f converter. The photon flux is normalized to the mean QBPM count rate, $\overline{n_{\text{QBPM}}}$. Therefore, $I_{\text{XBIC}}$ is a measure of the current relative to the mean photon flux, and not with regards to any total number of photons. The linear conversion factor contains the factor for the sinusoidal chopper modulation $W_{\text{ff}} = \sqrt{2}$, the amplification factors of the pre-amplifier $A_{\text{PA}}$ and lock-in amplifier $A_{\text{LIA}}$, as well as the ratio of the V2f's frequency range to its voltage acceptance range $C_{\text{V2f}}$.
  330.  
  331. The QBPM (\textit{Quadrupole Beam Position Monitor}) does not measure the absolute number of photons, but instead a current proportional to it. So, again, all XBIC measurements have to be seen relative to the sample and measurement condition. It is not possible to compute from this data a quantity such as the EQE (\textit{External Quantum Efficiency})\cite[p.11]{stuckelberger-2017-jmr}, nor is it possible to interpolate how much current or voltage the tested solar cell would produce if a larger area was illuminated with the same overall photon energy, because the physical processes are much different when a small part of the cell is illuminated by an X-Ray, compared to operation under normal sunlight illumination.
  332. Next, the normalized current is interpolated from a 201 $\times$ 200 px to a so-called \textit{supergrid} with 1000 $\times$ 1000 px. This allows for subpixel translations during image registration and more precise application of the Gaussian filter (see subsection \ref{sssec:image-registration}). The interpolation uses the nearest-neighbour method to account for deviations, caused by the fly-scan, of the dialed positions from the encoded positions (compare \cite{pyrlik-2020-pvsc}).
  333.  
  334. \begin{figure}[!ht]
  335. \centering
  336. \includegraphics[width=\textwidth]{Figures/AllXBICinput.png}
  337. \caption{All channels from the DAQ: a) lock-in X: -10 to 0V, b) lock-in X: 0 to +10V, c) lock-in R, d) direct signal -10 to 0V, e) direct signal 0 to +10V. Each column represents a bias in the voltage series -1V, -500mV, -200mV, -100mV, 0mV (no bias), +100mV, +200mV, +500mV and +600mV.}
  338. \label{fig:AllXBICinput}
  339. \end{figure}
  340.  
  341. If this is done, one can plot the input from the DAQ as shown in Figure \ref{fig:AllXBICinput}. Each column represents one scan. From left to right, the voltage series -1 V, -500 mV, -200 mV, -100 mV, 0 mV (no bias), +100 mV, +200 mV, +500 mV and +600 mV is shown. As mentioned in section \ref{sssec:pa-and-lia}, this order does not correspond to the order in time during the experiment, to be sensitive to hysteresis artifacts. The actual order of measurement was 0 (no bias), 100, -100, -200 and 200 mV. Since the in-phase X component of the lock-in signal can have both negative and positive values, it is separated into two channels. The first row shows the channel which records X within the interval from -10 to 0 V and the second row shows the channel which records X within the interval from 0 to +10 V. As can be seen from the overview, mainly one of the two channels contains a signal. This is because the phase shift remains relatively constant. If this was not the case, the Y channel would have to be taken into account. The third row contains the R amplitude of the lock-in signal and is the most valuable data for this analysis. The last two rows contain the direct signals from the pre-amplifier (compare Figure \ref{fig:setup-of-experiment}), where the fourth contains the data from -10 to 0V and the fifth contains the data from 0V to +10V. Here, only data is recorded in either one of those and the other contains noise, because the channels split the positive from the negative currents.
  342. Furthermore, one can see in the figure that sensible data is recorded only in the range from -200 to +200 mV. This is because larger voltages lead to higher direct currents such that a recovery of the AC signal by the LIA was not possible. Therefore, this analysis relies only on the scans from -200 to +200 mV.
  343.  
  344. \subsection{Gaussian Filter and Image Registration}
  345. \label{sssec:image-registration}
  346. The profile of the X-Ray beam is not quadratic like the pixels, but has a Gaussian profile which leads to overlapping for each scanned pixel. To account for this oversampling, a Gaussian filter needs to be applied. The height and width of the Gaussian filter kernel size must be odd, but can have different sizes, to properly define a middle point of the filter matrix. The standard deviation $\sigma$ of the Gaussian filter can be calculated in openCV from the kernel size K:
  347. \begin{equation}
  348. \label{eqn:sigma}
  349. \sigma = 0.3\cdot \left(\frac{K-1}{2}-1\right) + 0.8
  350. \end{equation}
  351.  
  352. When $\sigma = 1$, the filter will consider the next pixel only by half the amount of the middle pixel. Since the beam's FWHM is \SI{105}{nm} in one direction and one pixel in the supergrid corresponds to \SI{10}{nm}, $\sigma = 5.3$ or K = 33 is a good approximation. Thus, a Gaussian filter with a kernel size of 33 pixels is applied.
  353.  
  354. During such high-precision measurements, it is inevitable that the scans do not overlap precisely, due to unaccounted beam shift, sample or stage movement. From this the need for \textit{image registration} arises. The openCV for Python library offers the solution with its \texttt{findtransformEEC} function which calculates six affine transformation coefficients. The Gaussian filter helps to separate characteristic features from noise. The border values of the current grids where applying the filter would have to take into account pixels exceeding the grid are set to NaN.\footnote{One must be careful here that the Gaussian filter does not reflect the border values which would intensify the random distribution of border values. To be safe, a border of 1 px width can be manually set to NaN. Thereby, the filter will only be applied where the filter size does not exceed the grid.} The six coefficients are applied to each pixel (Y,Z) by the openCV \texttt{warpAffine} function as following:
  355.  
  356. \begin{equation}
  357. \label{eqn:warp-matrix}
  358. \begin{pmatrix}
  359. Y' \\
  360. Z' \\
  361. \end{pmatrix} =
  362. \begin{bmatrix}
  363. a_{11} & a_{12} \\
  364. a_{21} & a_{22} \\
  365. \end{bmatrix} \cdot
  366. \begin{pmatrix}
  367. Y \\
  368. Z \\
  369. \end{pmatrix} +
  370. \begin{bmatrix}
  371. b_{1} \\
  372. b_{2} \\
  373. \end{bmatrix}
  374. \end{equation}
  375.  
  376. a\textsubscript{11} to a\textsubscript{22} are the coefficients for rotation and shearing. b\textsubscript{1} and b\textsubscript{2} are the coefficients for translation. For this data, no shearing or rotation, but only translation needs to be corrected, thus a\textsubscript{11} = a\textsubscript{22} = 1 and a\textsubscript{12} = a\textsubscript{21} = 0, whereas b\textsubscript{1} and b\textsubscript{2} contain the pixel values for translation in Y and Z direction.
  377. \texttt{warpAffine} applies the coefficients to the supergrid by shifting it b\textsubscript{1} pixels in Y and b\textsubscript{2} pixels in Z position and inserts NaN values at the vacant pixels.
  378. For more details about the effect of different Gaussian filter sizes on the correlation coefficients and about affine transformations see \cite{pyrlik-2020-pvsc}. At last, all pixel positions (Y,Z), where any of the scans is NaN, all other scans at (Y,Z) will be set to NaN as well. Effectively, 894,915 or 89.5\% of the supergrid pixels are used for analysis.
  379.  
  380. \section{Results and Discussion}
  381. \label{sec:results-and-discussion}
  382.  
  383. \subsection{Voltage dependence of spatially-resolved XBIC measurements}
  384. \label{ssec:voltage-dependence}
  385.  
  386. \begin{figure}
  387. \centering
  388.  
  389. \begin{minipage}[t]{0.4\textwidth}
  390. \centering \hspace{-0.8cm} \textbf{a)} \\
  391. \includegraphics[trim={0 1cm 0 0},clip, width=\textwidth]{Figures/singleplots/Scan185.png}
  392. \end{minipage}
  393. \hspace{1cm}
  394. \begin{minipage}[t]{0.4\textwidth}
  395. \centering \hspace{-1cm} \textbf{b)} \\
  396. \includegraphics[trim={0 1cm 0 0},clip, width=\textwidth]{Figures/singleplots/Scan185_direct.png}
  397. \end{minipage}
  398.  
  399. \begin{minipage}[t]{0.4\textwidth}
  400. \includegraphics[trim={0 1cm 0 0},clip, width=\textwidth]{Figures/singleplots/Scan183.png}
  401. \end{minipage}
  402. \hspace{1cm}
  403. \begin{minipage}[t]{0.4\textwidth}
  404. \includegraphics[trim={0 1cm 0 0},clip, width=\textwidth]{Figures/singleplots/Scan183_direct.png}
  405. \end{minipage}
  406.  
  407. \begin{minipage}[t]{0.4\textwidth}
  408. \includegraphics[trim={0 1cm 0 0},clip, width=\textwidth]{Figures/singleplots/Scan181.png}
  409. \end{minipage}
  410. \hspace{1cm}
  411. \begin{minipage}[t]{0.4\textwidth}
  412. \includegraphics[trim={0 1cm 0 0},clip, width=\textwidth]{Figures/singleplots/Scan181_direct.png}
  413. \end{minipage}
  414.  
  415. \begin{minipage}[t]{0.4\textwidth}
  416. \includegraphics[trim={0 1cm 0 0},clip, width=\textwidth]{Figures/singleplots/Scan182.png}
  417. \end{minipage}
  418. \hspace{1cm}
  419. \begin{minipage}[t]{0.4\textwidth}
  420. \includegraphics[trim={0 1cm 0 0},clip, width=\textwidth]{Figures/singleplots/Scan182_direct.png}
  421. \end{minipage}
  422.  
  423. \begin{minipage}[t]{0.4\textwidth}
  424. \includegraphics[width=\textwidth]{Figures/singleplots/Scan186.png}
  425. \end{minipage}
  426. \hspace{1cm}
  427. \begin{minipage}[t]{0.4\textwidth}
  428. \includegraphics[width=\textwidth]{Figures/singleplots/Scan186_direct.png}
  429. \end{minipage}
  430.  
  431. \caption{The data used for this analysis: a) XBIC from the lock-in amplifier, b) direct signal after pre-amplification. The direct signal carries partly the voltage-induced and partly the X-ray beam induced current. For the unbiased scan \#181, the same NaN border is added.}
  432. \label{fig:lia-and-pa-XBICplot}
  433.  
  434. \end{figure}
  435. After having performed the preparatory steps described in section \ref{sec:data-preparation}, the current signal from the lock-in and pre-amplifier channels can be plotted. Figure \ref{fig:lia-and-pa-XBICplot} shows the X-ray beam induced and the direct signal of the solar cell. The direct signal does not only carry the voltage-induced current, but also the current induced by the X-ray, therefore the same patterns are visible there. One notices that the lock-in signal exhibits a relatively small change of currents, whereas the direct signal changes from large positive to large negative values. Characteristic features are maintained in both sequences.
  436. The unbiased scan was taken as a reference in image registration for calculating the transformation coefficients that have been applied in the same way for the direct signal plots. The direct signal has a time component which was filtered out during the lock-in amplification. At least part of it can well be observed in the direct signal at bias; there is a regular pattern of dotted lines, possibly due to a varying irradiation of the X-ray beam. A possible explanation as to why the dotted lines exhibit relatively higher absolute currents both at reverse and forward bias will be given in subsection \ref{ssec:current-distribution}. Another time-dependent effect that might be fully or partially ascribed to external sources can be noted: Comparing the direct signal plots for +100 and +200 mV, it can be observed that with the latter there seems to be a distinction into the top \SI{3}{\micro\metre} and lower \SI{7}{\micro\metre} position of the solar cell sample, where both parts have a lower contrast than at +100 mV, but with lower absolute currents at the top and higher absolute currents at the bottom of the sample. Again, this could be caused by external sources in the microhutch. In the next subsection, this feature will be well visible in the current distribution, illustrated by Figure \subref{fig:violin_currents-direct-zeroed_numbered}. Whereas here the artifacts are visible as spatially resolved, in the following they are instead visible as quantitatively resolved.
  437.  
  438. \subsection{Current distribution}
  439. \label{ssec:current-distribution}
  440.  
  441. \begin{figure}
  442. \centering
  443. \begin{subfigure}{\linewidth}
  444. \includegraphics[width=0.8\textwidth]{Figures/Violins/violin_currents.png}
  445. \caption{Current distribution for each bias, -200 to +200mV from left to right. Violins are scaled by the number of counts and counts are taken from Seaborn's kernel density estimator. In the middle of each violin, a box plot with the 25\% and 75\% percentile as well as the median is shown.}
  446. \label{fig:violin_currents}
  447. \end{subfigure}
  448. \hspace{0.46cm}
  449. \centering
  450. \begin{subfigure}{\linewidth}
  451. \includegraphics[width=0.8\textwidth]{Figures/Violins/violin_currents_direct_zeroed_numbered.png}
  452. \caption{Same display style, but for the direct signal. Since there is a strong offset, the currents are subtracted by the offset of each scan from zero respectively, shown above each violin in (nA).}
  453. \label{fig:violin_currents-direct-zeroed_numbered}
  454. \end{subfigure}
  455. \end{figure}
  456.  
  457. In the previous subsection, it was observed that the characteristic features are maintained over different voltages. In the direct signal, both artifacts and peculiar features were observed.
  458. To confirm the maintenance of features, the distribution of currents is visualized in a histogram (as a violin plot) in Figure \subref{fig:violin_currents}. The distribution of discrete current values is interpolated to a continuous distribution in Seaborn with a kernel density estimator. We observe that most of the values are situated close to the median. With respect to the minimum and maximum, they are slightly asymmetrically distributed with a trend to lower currents. Outliers are found in both low and high currents, but they are still close to the main distribution. The shape of the distribution is apparently only slightly dependent on the bias voltage. The violin at +200 mV displays the largest deviation from the others, both in its shape and its median. The latter observation that the median shows a larger deviation at the highest forward bias of this measurement is important to remember. It will be further discussed in subsection \ref{sssec:good-and-bad-diodes}.
  459. Figure \subref{fig:violin_currents-direct-zeroed_numbered} illustrates the distribution of the direct signal and allows for a better examination of the peculiar features than Figure \ref{fig:lia-and-pa-XBICplot} from the previous subsection. The violins are plotted without an offset from zero, to allow for direct comparison.\footnote{Remember that in the forward bias regime the top/bottom of the violin is inverted in comparison to Figure \ref{fig:lia-and-pa-XBICplot}b), such that the bottom part of the last two violins here corresponds to high currents there.}
  460. We can observe that the (min,max) range of direct current values is larger than for the lock-in amplified currents. Above it was noted that at +200 mV the pixels appear to concentrate around lower currents at the top and around higher currents at the bottom position of the solar cell sample. Here, the increased number of pixels with low currents can be seen as the top of the +200 mV violin. The shape of the distributions is essentially retained until 0 mV, with a tendency towards less outliers, and in the forward bias regime inverted, with a tendency again towards more outliers. To explain the peculiar feature that the dotted lines stand out as lower under reverse and as higher under forward bias, an important observation must be first noted: The shapes of distributions from Figure \subref{fig:violin_currents} are not the same as in Figure \subref{fig:violin_currents-direct-zeroed_numbered}! The voltage-induced current reduces all violins by a certain amount, but apparently not all pixels are affected in the same way by the voltage. As mentioned in subsection \ref{ssec:pn-junction-bias-voltage}, the direction of the photo-current remains the same during reverse and forward bias, but the voltage-induced (dark) current changes its direction from positive to negative, because the diffusion current dominates the thermal-generation current. Therefore, the violin shape could be expected to be exactly retained. Thus, there seem to be three factors at play:
  461. \begin{itemize}
  462. \item external sources such as varying beam intensity or electromagnetic field perturbations in the microhutch during the measurement causing arbitrary artifacts
  463. \item voltage-dependent factor causing the shape to be not retained
  464. \item voltage-independent causing the shape to be \textit{still vaguely} similar
  465. \end{itemize}
  466. An investigation of external sources that might cause or have caused artifacts in the measurement is beyond the scope of this study. It can be mentioned again though that the dotted lines are caused by the varying X-ray beam intensity, as cross-checked with the QBPM count data.
  467. Voltage-dependent factors are a hint that some pixels display \textit{good} and some \textit{bad diodes}. The distinction will be further elaborated in subsection \ref{sssec:good-and-bad-diodes}. For now it suffices to say that the this behaviour can possibly be ascribed to the fact that the X-ray beam strongly increases the conductivity of the solar cell, perhaps introducing shunts (compare \cite{stuckelberger-2017-pvsc-xbiv}). However, the variations in X-ray beam do not seem to have an effect on the conductivity, since cross-checking Figure \ref{fig:lia-and-pa-XBICplot}b) with the QBPM counts shows that the dots represent a \textit{reduced} X-ray beam intensity.
  468. Main features are still vaguely visible in Figure \subref{fig:violin_currents-direct-zeroed_numbered} and even better by comparing Figure \ref{fig:lia-and-pa-XBICplot}a) and \ref{fig:lia-and-pa-XBICplot}b). Therefore, conclusions can be drawn with regards to the localization of voltage-independent regions on the solar cell sample, but not yet about regions of voltage-dependent regions. This will be attempted in the next section.
  469. The main cause for voltage-independent behaviour in solar cells is the variation in thickness of the absorber layer. For thin-film materials such as CIGS and sufficiently low photon energy, the factor can be assumed to be linear in first approximation\cite{stuckelberger-2017-pvsc-xbiv}.
  470.  
  471. Furthermore, the unchanged shape from the lock-in signal in Figure \subref{fig:violin_currents} leads to the question whether the pixels of the solar cell with high currents retain high currents and if the same is valid for the pixels with low currents. To answer this with certainty, the data needs to be viewed \textit{both} as spatially resolved \textit{and} quantitatively resolved.
  472.  
  473. \subsubsection{Pixels with high and low current}
  474. \label{sssec:pixels-with-high-and-with-low-current}
  475.  
  476. \begin{figure}[!ht]
  477. \centering
  478. \includegraphics[width=\textwidth]{Figures/MedianPlots/MedianGroups.png}
  479. \caption{The values above the median are grouped in group 1 and the values below the median in group 2 for each bias voltage, from a) -200mV to e) +200mV. The Pearson correlation coefficient lies between 0.91 and 0.94.}
  480. \label{fig:MedianGroups}
  481. \end{figure}
  482.  
  483. While subsection \ref{ssec:voltage-dependence} gave a first impression on the spatial resolution of currents, it didn't allow for direct quantitative comparison. This was complemented with subsection \ref{ssec:current-distribution}, but the spatial resolution was lost in the one-dimensional distribution. Since one might arrive at the same result by looking at \ref{fig:lia-and-pa-XBICplot}a), the aim of this subsection is to exactly confirm this conjecture.
  484. Here, each pixel is defined as having a low or high current with respect to the median. All pixels above the median are said to have a high current and all pixels below the median are said to have a low current. As illustrated in Figure \ref{fig:MedianGroups}, the pixels with high currents are in group 1 and the pixels with low currents in group 2.\footnote{The groups are called \textit{cluster groups} in accordance with the clustering method in subsection \ref{ssec:local-iv-curves}, where there are also distinctions into two groups where group 1 largely coincides with the regions with high currents.} When the plots are compared successively, Pearson correlation coefficients between 0.91 and 0.94 are yielded. In other words, 81.1\% of the examined pixels display the behaviour that they don't change their allocation to either group 1 or 2 at all. The other 18.9\% are situated at the edges of group 1 and 2, as shown in Figure \ref{fig:MedianGroups_changing}.
  485. \begin{figure}[!ht]
  486. \centering
  487. \includegraphics[width=0.25\textwidth]{Figures/MedianPlots/MedianGroups_changing.png}
  488. \caption{All values that change at least once, from Figure \ref{fig:MedianGroups}a) to e), their allocation to group 1 or 2 are grouped in group 2, the other 81.1\% in group 1.}
  489. \label{fig:MedianGroups_changing}
  490. \end{figure}
  491. The latter can probably be attributed to the fact that those values of group 2 in Figure \ref{fig:MedianGroups_changing} are situated much closer around the median than those from group 2 in Figure \ref{fig:MedianGroups}a) to e).
  492.  
  493. \subsection{Local \textit{IV} curves}
  494. \label{ssec:local-iv-curves}
  495.  
  496. \begin{figure}
  497. \centering
  498. \begin{subfigure}{0.46\linewidth}
  499. \includegraphics[trim={0 0 34.5cm 0},clip, width=\linewidth]{Figures/plotIV.png}
  500. \caption{In the background, a sample of 64 IV curves, uniformly spread across the current grid, and in bold, the mean IV curve for \textit{all} pixels are displayed.}
  501. \label{fig:plotIV}
  502. \end{subfigure}
  503. \hspace{0.6cm}
  504. \centering
  505. \begin{subfigure}{0.46\linewidth}
  506. \includegraphics[trim={0 0 34.5cm 0},clip, width=\linewidth]{Figures/plotIV_direct.png}
  507. \caption{The same display style as for Figure \subref{fig:plotIV}, but for the direct signal.}
  508. \label{fig:plotIV_direct}
  509. \end{subfigure}
  510. \end{figure}
  511.  
  512. Hitherto, the scans have been plotted as spatially resolved for different voltages into the dimension 5 $\times$ 1000 $\times$ 1000 pixels in subsection \ref{ssec:voltage-dependence}, as quantitatively resolved into the dimension 5 $\times$ 1,000,000 current values in subsection \ref{ssec:current-distribution} and as a combination of both, but only with regards to a single criterion, in subsection \ref{sssec:pixels-with-high-and-with-low-current}.\footnote{Barring the 105,085 NaN values for each of the five scans (subsection \ref{sssec:image-registration}).}
  513. Now, the currents are examined not as single values anymore, but as 1000 $\times$ 1000 $\times$ 5 curves, retaining the (Y,Z) coordinate for each IV curve.
  514. Each of the curves has a different offset from zero and a different relative shape. Both the absolute and relative values are considered here. To get a first impression, the mean IV curves for all pixels, including exemplary single IV curves, are plotted in Figures \subref{fig:plotIV} and \subref{fig:plotIV_direct}. Both the XBIC and In Figure \subref{fig:plotIV_direct}, which shows the direct signal, it can be seen that the curve differs from the exponential curve of the Shockley equation (Equation \ref{eqn:shockley}.
  515. (this is due to the higher conductivity)
  516. However, Figure \subref{fig:plotIV} neither represents the Shockley equation. Fitting the curve to Equation \ref{eqn:shockley} was not possible.
  517.  
  518. A rather open research question was asked at the beginning of the data analysis, how to classify all IV curves. To answer this, the unsupervised machine learning k-means clustering algorithm has been used. As an input, the algorithm takes data points with the same number of (any) dimension and randomly placed k points, so-called \textit{centroids}. Each data point is assigned to the centroid where the Euclidian distance is minimized. Then each centroid is replaced by a new one that corresponds to the mean of all of its data points. This algorithm proceeds until the centroids remain the same for all groups.
  519. \begin{figure}
  520. \centering
  521. \begin{subfigure}{0.46\linewidth}
  522. \includegraphics[width=1.2\linewidth]{Figures/ClusterPlots/ClusterCurves_2-k_1-kmeans-input_False-normalize-samples.png}
  523. \caption{In the background, a sample of 26/38 IV curves for group 1/2, uniformly spread across the current grid, and in bold, the mean IV curve for \textit{all} pixels of the respective group are displayed. K-means input: absolute values.}
  524. \label{fig:plotClusterCurves1}
  525. \end{subfigure}
  526. \hspace{0.6cm}
  527. \centering
  528. \begin{subfigure}{0.46\linewidth}
  529. \includegraphics[width=\linewidth]{Figures/ClusterPlots/ClusterPosition_2-k_1-kmeans-input_False-normalize-samples.png}
  530. \caption{Group 1 and 2 for the absolute values as k-means input.}
  531. \label{fig:plotClusterPosition1}
  532. \end{subfigure}
  533. \end{figure}
  534.  
  535. \subsubsection{Good and bad diodes}
  536. \label{sssec:good-and-bad-diodes}
  537.  
  538. --> reference to \ref{ssec:current-distribution}: high conductivity (maximal: shunts) is basically the same as good and bad diode, but the behaviour here in the IV curves is still different
  539.  
  540. Different derivations of the original XBIC data can be used as an input for kmeans clustering. If the IV curves with their five absolute values are taken as an input for k-means clustering, the result in Figure \subref{fig:plotClusterPosition1} resembles the original XBIC map from Figure \ref{fig:lia-and-pa-XBICplot} with a Pearson correlation coefficient of -0.82. This is clearly the case, because the relative changes of the I(V) values are roughly about ten times smaller compared to the absolute current values. However, using the slope as a relative measure did not result in any meaningful clustering: The mean IV curves did not show any significant difference. The Pearson correlation coefficient between the plot of good and bad diodes in Figure \ref{fig:MedianGroups} and Figure \subref{fig:plotClusterPosition1} is 0.87.
  541. \begin{figure}
  542. \centering
  543. \begin{subfigure}{0.46\linewidth}
  544. \includegraphics[width=1.2\linewidth]{Figures/ClusterPlots/ClusterCurves_2-k_1-kmeans-input_False-normalize-samples_True-normalize-curves.png}
  545. \caption{Same display style as Figure \subref{fig:plotClusterCurves1}, but the mean IV curves are normalized to their value at I(V=0mV).}
  546. \label{fig:plotClusterCurves1_normalized}
  547. \end{subfigure}
  548. \hspace{0.6cm}
  549. \centering
  550. \begin{subfigure}{0.46\linewidth}
  551. \includegraphics[width=1.2\linewidth]{Figures/ClusterPlots/ClusterCurves_2-k_7-kmeans-input_False-normalize-samples_True-normalize-curves.png}
  552. \caption{Same display style as Figure \subref{fig:plotClusterCurves1_normalized}, but here the slope from +100 to +200 mV serves as an input for k-means clustering.}
  553. \label{fig:plotClusterCurves7_normalized}
  554. \end{subfigure}
  555. \end{figure}
  556. Therefore, Figure \subref{fig:plotClusterCurves1_normalized} can also serve as an example of the mean IV curves of good and bad diodes. In there, it can be seen that the current difference between +100 and +200 mV is smaller than in Figure \subref{fig:plotClusterCurves7_normalized}. If it is assumed that the classification into good and bad diodes should be made on the basis of the slope within this region, then Figure \subref{fig:plotClusterPosition7} serves as a better classification than the grouping with the absolute current values (Figure \subref{fig:plotClusterPosition1}) or the grouping with respect to the median (Figure \subref{fig:MedianGroups}). However, the plots clearly differ from each other (with a Pearson correlation coefficient of 0.25 to 0.27). This leads to the hypothesis that only for the pixels in Figure \subref{fig:plotClusterPosition7}
  557.  
  558.  
  559. Generally, combining different slopes as an input for the clustering didn't yield much different results as when using a one-dimensional input, such as the total slope from the current at -200 mV to +200 mV. For this, the slope and intercept were fitted by minimizing the sum of squared distances to the fitted line of all five points.
  560.  
  561. \begin{figure}
  562. \centering
  563. \begin{subfigure}{0.46\linewidth}
  564. \includegraphics[width=1.2\linewidth]{Figures/ClusterPlots/ClusterCurves_2-k_7-kmeans-input_False-normalize-samples.png}
  565. \caption{In the background, a sample of 24/40 IV curves for group 1/2, uniformly spread across the current grid, and in bold, the mean IV curve for \textit{all} pixels of the respective group are displayed. K-means input: current difference between +100 and +200mV.}
  566. \label{fig:plotClusterCurves7}
  567. \end{subfigure}
  568. \hspace{0.6cm}
  569. \centering
  570. \begin{subfigure}{0.46\linewidth}
  571. \includegraphics[width=\linewidth]{Figures/ClusterPlots/ClusterPosition_2-k_7-kmeans-input_False-normalize-samples.png}
  572. \caption{Group 1 and 2 for the current difference between +100 and +200mV as k-means input.}
  573. \label{fig:plotClusterPosition7}
  574. \end{subfigure}
  575. \end{figure}
  576.  
  577. In Figure \subref{fig:plotClusterCurves8}, it can be seen that those curves which show a smaller slope (group 1) also have a higher current. Figure \subref{fig:plotClusterPosition8} therefore resembles the current map. Taking a single slope from one (I,V) point to its next, does not yield any significant difference in the mean IV curves for each group, except for the slope from +100 to +200 mV. This is not surprising, because the performance of a good diode differs by its behaviour in the forward bias region.
  578.  
  579. \begin{figure}[!ht]
  580. \centering
  581. \includegraphics[width=0.7\textwidth]{Figures/ShockleyCurves/ShockleyCurve_normalized.png}
  582. \caption{Comparison of two curves from the Shockley equation. The "Bad diode" has half of the saturation current $I_{\text{0}}$ as the "Good diode". Two random examples of IV curves are included, from group 1 and 2 in Figure \ref{fig:MedianGroups}.}
  583. \label{fig:ShockleyCurve}
  584. \end{figure}
  585.  
  586. Figure \ref{fig:ShockleyCurve} illustrates the difference of a high- and low-performing diode. Two different saturation currents $I_{\text{0}}$ (compare subsection \ref{ssec:pn-junction-bias-voltage} and \ref{ssec:important-parameters}) were chosen. As can be seen in the figure, the curves only significantly differ within the forward bias region from +100 to +200 mV.
  587.  
  588.  
  589. \begin{figure}
  590. \centering
  591. \begin{subfigure}{0.46\linewidth}
  592. \includegraphics[trim={0 0 0 1.5cm},clip, width=1.2\linewidth]{Figures/ClusterPlots/ClusterCurves_2-k_7-kmeans-input_False-normalize-samples.png}
  593. \caption{In the background, a sample of 24/40 IV curves for group 1/2, uniformly spread across the current grid, and in bold, the mean IV curve for \textit{all} pixels of the respective group are displayed. K-means input: slope from I(+100mV) to I(+200mV).}
  594. \label{fig:plotClusterCurves8}
  595. \end{subfigure}
  596. \hspace{0.6cm}
  597. \centering
  598. \begin{subfigure}{0.46\linewidth}
  599. \includegraphics[trim={0 0 0 1.5cm},clip, width=\linewidth]{Figures/ClusterPlots/ClusterPosition_2-k_7-kmeans-input_False-normalize-samples.png}
  600. \caption{Group 1 and 2 for the slope from I(+100mV) to I(+200mV) as k-means input.}
  601. \label{fig:plotClusterPosition8}
  602. \end{subfigure}
  603. \end{figure}
  604.  
  605. Interestingly, the curves are clearly shifted, although the clustering only received the relative slope as an input. To reconfirm that, a violin plot can be made that shows
  606.  
  607. \begin{figure}[!ht]
  608. \centering
  609. \includegraphics[width=0.8\textwidth]{Figures/Violins/violin_groups_7.png}
  610. \caption{The current distribution for each scan, separated into group 1 and 2 for clustering with the }
  611. \label{fig:violin_groups_7}
  612. \end{figure}
  613.  
  614. \section{Conclusions}
  615. --- Das ist eher ein Ausblick:
  616. This vast combination of different techniques under hard X-ray illumination allows for point-by-point correlation in 3D space, however still limited by the necessity of thorough 3D simulations to complement the measurements, increased perturbations due to X-ray beam induced temperature, bias voltage, bias light and locally generated electric fields and, most importantly, the sheer amount of data that needs to processed and still have pending reconstruction algorithms, possibly being alleviated by machine learning in the future\cite{stuckelberger-2017-jmr}.
  617. The hard X-ray beam allows for deeper penetration depth, revealing more of the nanoscale ingredients.
  618.  
  619.  
  620. \section*{References}
  621. \addcontentsline{toc}{section}{References}
  622. \printbibliography[heading=none]
  623.  
  624.  
  625. \section*{Acknowledgements}
  626. Above all, I would like to express my gratitude to Michael E. Stueckelberger who welcomed me as a student of special education to the research group, guided me through the different steps of the data analysis, always having an open ear for my questions, and who I will especially remember for creating such a positive work environment.
  627. \\~\\
  628. I want to show my appreciation to Christian G. Schroer for the opportunity to conduct research in the FS-PETRA group.
  629. \\~\\
  630. Furthermore, I feel obliged to pay my special regards to Christina Ossig with whom I delved into long discussions of solar cell physics and lock-in amplification. I am indebted to Niklas Pyrlik, my fellow bachelor thesis student, who helped whenever I encountered seemingly insoluble problems in Python programming and whose research work served as an inspiration for this thesis.
  631. I also want to thank the rest of the research group for the discussions and inputs at our Zoom meetings which includes Giovanni Fevola, Catharina Ziska, Leutrim Pacarizi, Anna Barinskaya, Jackson Barp, Gero Falkenberg and Tizian Kaiser.
  632. \\~\\
  633. For taking the time to explain me the working principle of the quadrupole beam position monitor (QBPM) at the P06, I thank Frank Seiboth.
  634. \\~\\
  635. I acknowledge Empa for providing the solar cell.
  636. \\~\\
  637. Last, the research group acknowledges DESY (Hamburg, Germany), a member of the Helmholtz Association HGF, for the provision of experimental facilities. This research was carried out at PETRA III.
  638.  
  639. \newpage
  640. \section*{Erklärung}
  641. \addcontentsline{toc}{section}{Erklärung} % add section* to the ToC
  642. \thispagestyle{empty}
  643. Hiermit versichere ich an Eides statt, dass ich die Arbeit eigenständig verfasst habe. Ich habe keine anderen als die angegebenen Quellen und Hilfsmittel benutzt sowie wörtliche und sinngemäße Zitate kenntlich gemacht. Die Arbeit hat in gleicher oder ähnlicher Form noch keiner Prüfungsbehörde vorgelegen und die eingereichte schriftliche Fassung entspricht der auf dem elektronischen Speichermedium.
  644.  
  645. Ich bin damit \emph{nicht} einverstanden, dass die Bachelorarbeit ver\"offentlicht wird.
  646. \vspace{1cm}
  647. \begin{center}
  648. \begin{tabular}[h]{lp{2cm}p{5.5cm}}
  649. Hamburg, \today & & \\
  650. \cline{1-1}\cline{3-3}
  651. Ort, Datum& & Jan Hense\\
  652. \end{tabular}
  653. \end{center}
  654.  
  655.  
  656.  
  657. %\newpage
  658. %\begin{appendices}
  659. %\section{Appendix A}
  660. %\label{app:pyth}
  661.  
  662. %\end{appendices}
  663. %\pagebreak
  664.  
  665. %\printbibliography
  666. \end{document}
Add Comment
Please, Sign In to add comment
Public Pastes
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!