Mixed finite-element methods for elliptic convection-dominated problems arising in semiconductor physics

In this work we develop and analyze an adaptive numerical scheme for simulating a class of macroscopic semiconductor models. At first the numerical modelling of semiconductors is reviewed in order to classify the Energy-Transport models for semiconductors that are later simulated in 2D. In this class of models the flow of charged particles, that are negatively charged electrons and so-called holes, which are quasi-particles of positive charge, as well as their energy distributions are described by a coupled system of nonlinear partial differential equations. A considerable difficulty in simulating these convection-dominated equations is posed by the nonlinear coupling as well as due to the fact that the local phenomena such as "hot electron effects" are only partially assessable through the given data. The primary variables that are used in the simulations are the particle density and the particle energy density. The user of these simulations is mostly interested in the current flow through parts of the domain boundary - the contacts. The numerical method considered here utilizes mixed finite-elements as trial functions for the discrete solution. The continuous discretization of the normal fluxes is the most important property of this discretization from the users perspective. It will be proven that under certain assumptions on the triangulation the particle density remains positive in the iterative solution algorithm. Connected to this result an a priori error estimate for the discrete solution of linear convection-diffusion equations is derived. The local charge transport phenomena will be resolved by an adaptive algorithm, which is based on a posteriori error estimators. At that stage a comparison of different estimations is performed. Additionally a method to effectively estimate the error in local quantities derived from the solution, so-called "functional outputs", is developed by transferring the dual weighted residual method to mixed finite elements. For a model problem we present how this method can deliver promising results even when standard error estimator fail completely to reduce the error in an iterative mesh refinement process.

Synthesis and surface modification of semiconductor nanocrystals

The last decade has witnessed an exponential growth of activities in the field of nanoscience and nanotechnology worldwide, driven both by the excitement of understanding new science and by the potential hope for applications and economic impacts. The largest activity in this field up to date has been in the synthesis and characterization of new materials consisting of particles with dimensions in the order of a few nanometers, so-called nanocrystalline materials. [1-8] Semiconductor nanomaterials such as III/V or II/VI compound semiconductors exhibit strong quantum confinement behavior in the size range from 1 to 10 nm. Therefore, preparation of high quality semiconductor nanocrystals has been a challenge for synthetic chemists, leading to the recent rapid progress in delivering a wide variety of semiconducting nanomaterials. Semiconductor nanocrystals, also called quantum dots, possess physical properties distinctly different from those of the bulk material. Typically, in the size range from 1 to 10 nm, when the particle size is changed, the band gap between the valence and the conduction band will change, too. In a simple approximation a particle in a box model has been used to describe the phenomenon[9]: at nanoscale dimensions the degenerate energy states of a semiconductor separate into discrete states and the system behaves like one big molecule. The size-dependent transformation of the energy levels of the particles is called “quantum size-effect”. Quantum confinement of both the electron and hole in all three dimensions leads to an increase in the effective bandgap of the material with decreasing crystallite size. Consequently, both the optical absorption and emission of semiconductor nanaocrystals shift to the blue (higher energies) as the size of the particles gets smaller. This color tuning is well documented for CdSe nanocrystals whose absorption and emission covers almost the whole visible spectral range. As particle sizes become smaller the ratio of surface atoms to those in the interior increases, which has a strong impact on particle properties, too. Prominent examples are the low melting point [8] and size/shape dependent pressure resistance [10] of semiconductor nanocrystals. Given the size dependence of particle properties, chemists and material scientists now have the unique opportunity to change the electronic and chemical properties of a material by simply controlling the particle size. In particular, CdSe nanocrystals have been widely investigated. Mainly due to their size-dependent optoelectronic properties [11, 12] and flexible chemical processibility [13], they have played a distinguished role for a number of seminal studies [11, 12, 14, 15]. Potential technical applications have been discussed, too. [8, 16-27] Improvement of the optoelectronic properties of semiconductor nanocrystals is still a prominent research topic. One of the most important approaches is fabricating composite type-I core-shell structures which exhibit improved properties, making them attractive from both a fundamental and a practical point of view. Overcoating of nanocrystallites with higher band gap inorganic materials has been shown to increase the photoluminescence quantum yields by eliminating surface nonradiative recombination sites. [28] Particles passivated with inorganic shells are more robust than nanocrystals covered by organic ligands only and have greater tolerance to processing conditions necessary for incorporation into solid state structures or for other applications. Some examples of core-shell nanocrystals reported earlier include CdS on CdSe [29], CdSe on CdS, [30], ZnS on CdS, [31] ZnS on CdSe[28, 32], ZnSe on CdSe [33] and CdS/HgS/CdS [34]. The characterization and preparation of a new core-shell structure, CdSe nanocrystals overcoated by different shells (CdS, ZnS), is presented in chapter 4. Type-I core-shell structures as mentioned above greatly improve the photoluminescence quantum yield and chemical and photochemical stability of nanocrystals. The emission wavelengths of type-I core/shell nanocrystals typically only shows a small red-shift when compared to the plain core nanocrystals. [30, 31, 35] In contrast to type-I core-shell nanocrystals, only few studies have been conducted on colloidal type-II core/shell structures [36-38] which are characterized by a staggered alignment of conduction and valence bands giving rise to a broad tunability of absorption and emission wavelengths, as was shown for CdTe/CdSe core-shell nanocrystals. [36] The emission of type-II core/shell nanocrystals mainly originates from the radiative recombination of electron-hole pairs across the core-shell interface leading to a long photoluminescence lifetime. Type-II core/shell nanocrystals are promising with respect to photoconduction or photovoltaic applications as has been discussed in the literature.[39] Novel type-II core-shell structures with ZnTe cores are reported in chapter 5. The recent progress in the shape control of semiconductor nanocrystals opens new fields of applications. For instance, rod shaped CdSe nanocrystals can enhance the photo-electro conversion efficiency of photovoltaic cells, [40, 41] and also allow for polarized emission in light emitting diodes. [42, 43] Shape control of anisotropic nanocrystals can be achieved by the use of surfactants, [44, 45] regular or inverse micelles as regulating agents, [46, 47] electrochemical processes, [48] template-assisted [49, 50] and solution-liquid-solution (SLS) growth mechnism. [51-53] Recently, formation of various CdSe nanocrystal shapes has been reported by the groups of Alivisatos [54] and Peng, [55] respectively. Furthermore, it has been reported by the group of Prasad [56] that noble metal nanoparticles can induce anisotropic growth of CdSe nanocrystals at lower temperatures than typically used in other methods for preparing anisotropic CdSe structures. Although several approaches for anisotropic crystal growth have been reported by now, developing new synthetic methods for the shape control of colloidal semiconductor nanocrystals remains an important goal. Accordingly, we have attempted to utilize a crystal phase control approach for the controllable synthesis of colloidal ZnE/CdSe (E = S, Se, Te) heterostructures in a variety of morphologies. The complex heterostructures obtained are presented in chapter 6. The unique optical properties of nanocrystals make them appealing as in vivo and in vitro fluorophores in a variety of biological and chemical investigations, in which traditional fluorescence labels based on organic molecules fall short of providing long-term stability and simultaneous detection of multiple emission colours [References]. The ability to prepare water soluble nanocrystals with high stability and quantum yield has led to promising applications in cellular labeling, [57, 58] deep-tissue imaging, [59, 60] and assay labeling [61, 62]. Furthermore, appropriately solubilized nanocrystals have been used as donors in fluorescence resonance energy transfer (FRET) couples. [63-65] Despite recent progress, much work still needs to be done to achieve reproducible and robust surface functionalization and develop flexible (bio-) conjugation techniques. Based on multi-shell CdSe nanocrystals, several new solubilization and ligand exchange protocols have been developed which are presented in chapter 7. The organization of this thesis is as follows: A short overview describing synthesis and properties of CdSe nanocrystals is given in chapter 2. Chapter 3 is the experimental part providing some background information about the optical and analytical methods used in this thesis. The following chapters report the results of this work: synthesis and characterization of type-I multi-shell and type-II core/shell nanocrystals are described in chapter 4 and chapter 5, respectively. In chapter 6, a high–yield synthesis of various CdSe architectures by crystal phase control is reported. Experiments about surface modification of nanocrystals are described in chapter 7. At last, a short summary of the results is given in chapter 8.

Numerical coupling of thermal-electric network models and energy-transport equations including optoelectronic semiconductor devices

In this work the numerical coupling of thermal and electric network models with model equations for optoelectronic semiconductor devices is presented. Modified nodal analysis (MNA) is applied to model electric networks. Thermal effects are modeled by an accompanying thermal network. Semiconductor devices are modeled by the energy-transport model, that allows for thermal effects. The energy-transport model is expandend to a model for optoelectronic semiconductor devices. The temperature of the crystal lattice of the semiconductor devices is modeled by the heat flow eqaution. The corresponding heat source term is derived under thermodynamical and phenomenological considerations of energy fluxes. The energy-transport model is coupled directly into the network equations and the heat flow equation for the lattice temperature is coupled directly into the accompanying thermal network. The coupled thermal-electric network-device model results in a system of partial differential-algebraic equations (PDAE). Numerical examples are presented for the coupling of network- and one-dimensional semiconductor equations. Hybridized mixed finite elements are applied for the space discretization of the semiconductor equations. Backward difference formluas are applied for time discretization. Thus, positivity of charge carrier densities and continuity of the current density is guaranteed even for the coupled model.

Detecting interfaces in a parabolic-elliptic problem from surface measurements

Assuming that the heat capacity of a body is negligible outside certain inclusions the heat equation degenerates to a parabolic-elliptic interface problem. In this work we aim to detect these interfaces from thermal measurements on the surface of the body. We deduce an equivalent variational formulation for the parabolic-elliptic problem and give a new proof of the unique solvability based on Lions’s projection lemma. For the case that the heat conductivity is higher inside the inclusions, we develop an adaptation of the factorization method to this time-dependent problem. In particular this shows that the locations of the interfaces are uniquely determined by boundary measurements. The method also yields to a numerical algorithm to recover the inclusions and thus the interfaces. We demonstrate how measurement data can be simulated numerically by a coupling of a finite element method with a boundary element method, and finally we present some numerical results for the inverse problem.

Fluorescence close to dielectric interfaces

The aim of the work was to study the correlation between the orientation and excited-state lifetimes of organic dyes close to dielectric interfaces. For this purpose, an experimental setup was designed and built, guiding the light through a prism in total internal reflection geometry. Fluorescence intensities and lifetimes for an ensemble of dye molecules were analyzed as a function of the excitation and detection polarizations. Working close to the total internal reflection angle, the differences between polarization combinations were enhanced. A classical electromagnetic model that assumes a chromophore as a couple of point-like electrical dipoles was developed. A numerical method to calculate the excitation and emission of dye molecules embedded in a multilayer system was implemented, by which full simulation of the time resolved fluorescence experiments was achieved. Free organic dyes and organic dyes covalently bound to polyelectrolyte chains were used. The polymer functionalization process avoided aggregation and provided control over the dyes position, within a few nanometers to the interface. Moreover, by varying the pH, the polymer chains could be deposited on different substrates with different conformations and the resulting fluorescence characteristics analyzed. Initially the fluorescence of organic dyes embedded in a polymer matrix was studied as a function of the distance between the fluorophores and the polymer-air interface. The non-radiative decay rate, vacuum decay rate and the relative angle between the excitation and emission dipoles of the chromophores could be determined. Different free organic dyes were deposited onto different dielectric spacers, as close as possible to the air-dielectric interface. Surprisingly, the fluorescence characteristics of dyes deposited onto polyelectrolyte layer were in good agreement with theoretical predictions of dyes in a polymer matrix, even when the layer was only 2 nm thick. When functionalized chains were deposited at low pH, on top of a polyelectrolyte spacer, the fluorescence had the characteristics of emitters embedded in a polymer matrix as well. Surface deposition at high pH showed an intermediate behaviour between emitters embedded in polymer and on top of the surface, in air. In general, for low pH values, the chains are deposited on a substrate in a train-like conformation. For high pH values, the chains are deposited in a loop-like conformation. As a consequence at low pH the functionalized polymer strongly interdigitates with the polyelectrolyte chains of the spacer, bringing most of the dyes inside the polymer. Thus, the fluorophores may experience the polymer as surrounding environment. On the other hand, for high pH values the dye-loaded chains adsorbed have a conformational arrangement of dense loops that extend away from the surface. Therefore many fluorophores experience the air as surrounding environment. Changing the spacer from polyelectrolyte to negatively charged silane produced contradictory results for lifetimes and intensities. The fluorescence intensities indicated the behaviour of emitters embedded in a polymer matrix, regardless of the pH value. On the other hand, for low pH values, the excited-state lifetimes showed that the emitters behaved as in air. For higher pH values, an intermediate behaviour between fluorophores located within and above of a dielectric film was observed. The poor agreement between theoretical and experimental data may be due to the simplified model utilized, by which the dipoles are assumed either in one side or in the other with respect to a geometrical air-dielectric interface. In the case when the dielectric film is constituted by the functionalized polymer chains themselves, reality is more complex and a different model may apply. Nevertheless, possible applications of the technique arise from a qualitative analysis.

Thermal Casimir effect at fluid interfaces: Fluctuation-induced forces between anisotropic colloids

We study the effective interaction between two ellipsoidal particles at the interface of two fluid phases which are mediated by thermal fluctuations of the interface. Within a coarse-grained picture, the properties of fluid interfaces are very well described by an effective capillary wave Hamiltonian which governs both the equilibrium interface configuration and the thermal fluctuations (capillary waves) around this equilibrium (or mean-field) position. As postulated by the Goldstone theorem the capillary waves are long-range correlated. The interface breaks the continuous translational symmetry of the system, and in the limit of vanishing external fields - like gravity - it has to be accompanied by easily excitable long wavelength (Goldstone) modes – precisely the capillary waves. In this system the restriction of the long-ranged interface fluctuations by particles gives rise to fluctuation-induced forces which are equivalent to interactions of Casimir type and which are anisotropic in the interface plane. Since the position and the orientation of the colloids with respect to the interface normal may also fluctuate, this system is an example for the Casimir effect with fluctuating boundary conditions. In the approach taken here, the Casimir interaction is rewritten as the interaction between fluctuating multipole moments of an auxiliary charge density-like field defined on the area enclosed by the contact lines. These fluctuations are coupled to fluctuations of multipole moments of the contact line position (due to the possible position and orientational fluctuations of the colloids). We obtain explicit expressions for the behavior of the Casimir interaction at large distances for arbitrary ellipsoid aspect ratios. If colloid fluctuations are suppressed, the Casimir interaction at large distances is isotropic, attractive and long ranged (double-logarithmic in the distance). If, however, colloid fluctuations are included, the Casimir interaction at large distances changes to a power law in the inverse distance and becomes anisotropic. The leading power is 4 if only vertical fluctuations of the colloid center are allowed, and it becomes 8 if also orientational fluctuations are included.

Adsorption of macromolecules on interfaces studied by isothermal titration calorimetry

The mixing of nanoparticles with polymers to form composite materials has been applied for decades. They combine the advantages of polymers (e.g., elasticity, transparency, or dielectric properties) and inorganic nanoparticles (e.g., specific absorption of light, magneto resistance effects, chemical activity, and catalysis etc.). Nanocomposites exhibit several new characters that single-phase materials do not have. Filling the polymeric matrix with an inorganic material requires its homogeneous distribution in order to achieve the highest possible synergetic effect. To fulfill this requirement, the incompatibility between the filler and the matrix, originating from their opposite polarity, has to be resolved. A very important parameter here is the strength and irreversibility of the adsorption of the surface active compound on the inorganic material. In this work the Isothermal titration calorimetry (ITC) was applied as a method to quantify and investigate the adsorption process and binding efficiencies in organic-inorganic–hybrid-systems by determining the thermodynamic parameters (ΔH, ΔS, ΔG, KB as well as the stoichiometry n). These values provide quantification and detailed understanding of the adsorption process of surface active molecules onto inorganic particles. In this way, a direct correlation between the adsorption strength and structure of the surface active compounds can be achieved. Above all, knowledge of the adsorption mechanism in combination with the structure should facilitate a more rational design into the mainly empirically based production and optimization of nanocomposites.

Tracer diffusion at water-oil interfaces studied by fluorescence correlation spectroscopy

The adsorption of particles and surfactants at water-oil interfaces has attracted continuous attention because of its emulsion stabilizing effect and the possibility to form two-dimensional materials. Herein, I studied the interfacial diffusion of single molecules and nanoparticles at water-oil interfaces using fluorescence correlation spectroscopy. rnrnFluorescence correlation spectroscopy (FCS) is a promising technique to study diffusion of fluorescent tracers in diverse conditions. This technique monitors and analyzes the fluorescence fluctuation caused by single fluorescent tracers coming in and out of a diffraction-limited observation volume “one at a time”. Thus, this technique allows a combination of high precision, high spatial resolution and low tracer concentration. rnrnIn chapter 1, I discussed some controversial questions regarding the properties of water-hydrophobic interfaces and also introduced the current progress on the stability and dynamic of single nanoparticles at water-oil interfaces. The materials and setups I used in this thesis were summarized in chapter 2. rnrnIn chapter 3, I presented a new strategy to study the properties of water-oil interfaces. The two-dimensional diffusion of isolated molecular tracers at water/n-alkane interfaces was measured using fluorescence correlation spectroscopy. The diffusion coefficients of larger tracers with a hydrodynamic radius of 4.0 nm agreed well with the values calculated from the macroscopic viscosities of the two bulk phases. However, for small molecule tracers with hydrodynamic radii of only 1.0 and 0.6 nm, notable deviations were observed, indicating the existence of an interfacial region with a reduced effective viscosity. rnrnIn chapter 4, the interfacial diffusion of nanoparticles at water-oil interfaces was investigated using FCS. In stark contrast to the interfacial diffusion of molecular tracers, that of nanoparticles at any conditions is slower than the values calculated in accordance to the surrounding viscosity. The diffusion of nanoparticles at water-oil interfaces depended on the interfacial tension of liquid-liquid interfaces, the surface properties of nanoparticles, the particle sizes and the viscosities of surrounding liquid phases. In addition, the interfacial diffusion of nanoparticles with Janus motif is even slower than that of their symmetric counterparts. Based on the experimental results I obtained, I drew some possibilities to describe the origin of nanoparticle slowdown at water-oil interfaces.

Präparation von Bismutsulfid-Absorberschichten für Dünnschichtsolarzellen mittels physikalischer Gasphasenabscheidung

Angesichts der sich abzeichnenden Erschöpfung fossiler Ressourcen ist die Erforschung alternativer Energiequellen derzeit eines der meistbeachteten Forschungsgebiete. Durch ihr enormes Potential ist die Photovoltaik besonders im Fokus der Wissenschaft. Um großflächige Beschichtungsverfahren nutzen zu können, wird seit einigen Jahren auf dem Gebiet der Dünnschichtphotovoltaik intensiv geforscht. Jedoch sind die gegenwärtigen Solarzellenkonzepte allesamt durch die Verwendung giftiger (Cd, As) oder seltener Elemente (In, Ga) oder durch eine komplexe Phasenbildung in ihrem Potential beschränkt. Die Entwicklung alternativer Konzepte erscheint daher naheliegend.rnAufgrund dessen wurde in einem BMBF-geförderten Verbundprojekt die Abscheidung von Dünnschichten des binären Halbleiters Bi2S3 mittels physikalischer Gasphasenabscheidung mit dem Ziel der Etablierung als quasi-intrinsischer Absorber in Solarzellenstrukturen mit p-i-n-Schichtfolge hin untersucht.rnDurch sein von einem hochgradig anisotropen Bindungscharakter geprägtes Kristallwachstum war die Abscheidung glatter, einphasiger und für die Integration in eine Multischichtstruktur geeigneter Schichten mit Schichtdicken von einigen 100 nm eine der wichtigsten Herausforderungen. Die Auswirkungen der beiden Parameter Abscheidungstemperatur und Stöchiometrie wurden hinsichtlich ihrer Auswirkungen auf die relevanten Kenngrößen (wie Morphologie, Dotierungsdichte und Photolumineszenz) untersucht. Es gelang, erfolgreich polykristalline Schichten mit geeigneter Rauigkeit und einer Dotierungsdichte von n ≈ 2 1015cm-3 auf anwendungsrelevanten Substraten abzuscheiden, wobei eine besonders starke Abhängigkeit von der Gasphasenzusammensetzung ermittelt werden. Es konnten weiterhin die ersten Messungen der elektronischen Zustandsdichte unter Verwendung von Hochenergie-Photoemissionsspektroskopie durchgeführt werden, die insbesondere den Einfluss variabler Materialzusammensetzungen offenbarten.rnZum Nachweis der Eignung des Materials als Absorberschicht standen innerhalb des Projektes mit SnS, Cu2O und PbS prinzipiell geeignete p-Kontaktmaterialien zur Verfügung. Es konnten trotz der Verwendung besonders sauberer Abscheidungsmethoden im Vakuum keine funktionstüchtigen Solarzellen mit Bi2S3 deponiert werden. Jedoch war es unter Verwendung von Photoemissionspektroskopie möglich, die relevanten Grenzflächen zu spektroskopieren und die Ursachen für die Beobachtungen zu identifizieren. Zudem konnte erfolgreich die Notwendigkeit von Puffermaterialien bei der Bi2S3-Abscheidung nachgewiesen werden, um Oberflächenreaktionen zu unterbinden und die Transporteigenschaften an der Grenzfläche zu verbessern.rn

Structure and dynamics of ionic liquid,electrode interfaces

In the early 20th century, Gouy, Chapman, and Stern developed a theory to describe the capacitance and the spatial ion distribution of diluted electrolytes near an electrode. After a century of research, considerable progress has been made in the understanding of the electrolyte/electrode interface. However, its molecular-scale structure and its variation with an applied potential is still under debate. In particular for room-temperature ionic liquids, a new class of solventless electrolytes, the classical theories for the electrical double layer are not applicable. Recently, molecular dynamics simulations and phenomenological theories have attempted to explain the capacitance of the ionic liquid/electrode interface with the molecular-scale structure and dynamics of the ionic liquid near the electrode. rnHowever, experimental evidence is very limited. rnrnIn the presented study, the ion distribution of an ionic liquid near an electrode and its response to applied potentials was examined with sub-molecular resolution. For this purpose, a new sample chamber was constructed, allowing in situ high energy X-ray reflectivity experiments under potential control, as well as impedance spectroscopy measurements. The combination of structural information and electrochmical data provided a comprehensive picture of the electric double layer in ionic liquids. Oscillatory charge density profiles were found, consisting of alternating anion- and cation-enriched layers at both, cathodic and anodic, potentials. This structure was shown to arise from the same ion-ion correlations dominating the liquid bulk structure that were observed as a distinct X-ray diffraction peak. Therefore, existing physically motivated models were refined and verified by comparison with independent measurements. rnrnThe relaxation dynamics of the interfacial structure upon potential variation were studied by time resolved X-ray reflectivity experiments with sub-millisecond resolution. The observed relaxation times during charging/discharging are consistent with the impedance spectroscopy data revealing three processes of vastly different characteristic time-scales. Initially, the ion transport normal to the interface happens on a millisecond-scale. Another 100-millisecond-scale process is associated with molecular reorientation of electrode-adsorbed cations. Further, a minute-scale relaxation was observed, which is tentatively assigned to lateral ordering within the first layer.