The present thesis is dedicated to the efficient computation of electrostatic interactions in particle systems, which is of great importance in the field of molecular dynamics simulations. In order to compute the therefor required physical quantities with only O(N log N) arithmetic operations, so called particle-mesh methods make use of the well-known Ewald summation approach and the fast Fourier transform (FFT). Typically, such methods are able to handle systems of point charges subject to periodic boundary conditions in all spatial directions. However, periodicity is not always desired in all three dimensions and, moreover, also interactions to dipoles play an important role in many applications. Within the scope of the present work, we consider the particle-particle NFFT method (P²NFFT), a particle-mesh approach based on the fast Fourier transform for nonequispaced data (NFFT). An extension of this method for mixed periodic as well as open boundary conditions is presented. Furthermore, the method is appropriately modified in order to treat particle systems containing both charges and dipoles. Consequently, an efficient algorithm for mixed charge-dipole systems, that additionally allows a unified handling of various types of periodic boundary conditions, is presented for the first time. Appropriate error estimates as well as parameter tuning strategies are developed and verified by numerical examples. / Die vorliegende Arbeit widmet sich der Berechnung elektrostatischer Wechselwirkungen in Partikelsystemen, was beispielsweise im Bereich der molekulardynamischen Simulationen eine zentrale Rolle spielt. Um die dafür benötigten physikalischen Größen mit lediglich O(N log N) arithmetischen Operationen zu berechnen, nutzen sogenannte Teilchen-Gitter-Methoden die Ewald-Summation sowie die schnelle Fourier-Transformation (FFT). Typischerweise können derartige Verfahren Systeme von Punktladungen unter periodischen Randbedingungen in allen Raumrichtungen handhaben. Periodizität ist jedoch nicht immer bezüglich aller drei Dimensionen erwünscht. Des Weiteren spielen auch Wechselwirkungen zu Dipolen in vielen Anwendungen eine wichtige Rolle. Zentraler Gegenstand dieser Arbeit ist die Partikel-Partikel-NFFT Methode (P²NFFT), ein Teilchen-Gitter-Verfahren, welches auf der schnellen Fouriertransformation für nichtäquidistante Daten (NFFT) basiert. Eine Erweiterung dieses Verfahrens auf gemischt periodische sowie offene Randbedingungen wird vorgestellt. Außerdem wird die Methode für die Behandlung von Partikelsystemen, in denen sowohl Ladungen als auch Dipole vorliegen, angepasst. Somit wird erstmalig ein effizienter Algorithmus für gemischte Ladungs-Dipol-Systeme präsentiert, der zusätzlich die Behandlung sämtlicher Arten von Randbedingungen mit einem einheitlichen Zugang erlaubt. Entsprechende Fehlerabschätzungen sowie Strategien für die Parameterwahl werden entwickelt und anhand numerischer Beispiele verifiziert.

The present thesis provides a modularized view on the structure of fast numerical methods for computing Coulomb interactions between charged particles in three-dimensional space. Thereby, the common structure is given in terms of three self-contained algorithmic frameworks that are built on top of each other, namely fast Fourier transform (FFT), nonequispaced fast Fourier transform (NFFT) and NFFT based particle-mesh methods (P²NFFT). For each of these frameworks algorithmic enhancement and parallel implementations are presented with special emphasis on scalability up to hundreds of thousands of parallel processes. In the context of FFT massively parallel algorithms are composed from hardware adaptive low level modules provided by the FFTW software library. The new algorithmic NFFT concepts include pruned NFFT, interlacing, analytic differentiation, and optimized deconvolution in Fourier space with respect to a mean square aliasing error. Enabled by these generalized concepts it is shown that NFFT provides a unified access to particle-mesh methods. Especially, mixed-periodic boundary conditions are handled in a consistent way and interlacing can be incorporated more efficiently. Heuristic approaches for parameter tuning are presented on the basis of thorough error estimates. / Die vorliegende Dissertation beschreibt einen modularisierten Blick auf die Struktur schneller numerischer Methoden für die Berechnung der Coulomb-Wechselwirkungen zwischen Ladungen im dreidimensionalen Raum. Die gemeinsame Struktur ist geprägt durch drei selbstständige und auf einander aufbauenden Algorithmen, nämlich der schnellen Fourier-Transformation (FFT), der nicht äquidistanten schnellen Fourier-Transformation (NFFT) und der NFFT-basierten Teilchen-Gitter-Methode (P²NFFT). Für jeden dieser Algorithmen werden Verbesserungen und parallele Implementierungen vorgestellt mit besonderem Augenmerk auf massiv paralleler Skalierbarkeit. Im Kontext der FFT werden parallele Algorithmen aus den Hardware adaptiven Modulen der FFTW Softwarebibliothek zusammengesetzt. Die neuen NFFT-Konzepte beinhalten abgeschnittene NFFT, Versatz, analytische Differentiation und optimierte Entfaltung im Fourier-Raum bezüglich des mittleren quadratischen Aliasfehlers. Mit Hilfe dieser Verallgemeinerungen bietet die NFFT einen vereinheitlichten Zugang zu Teilchen-Gitter-Methoden. Insbesondere gemischt periodische Randbedingungen werden einheitlich behandelt und Versatz wird effizienter umgesetzt. Heuristiken für die Parameterwahl werden auf Basis sorgfältiger Fehlerabschätzungen angegeben.

This work contributes to the structure investigation on the basis of small-angle neutron scattering (SANS). A new analytical scattering function for polydispers precipitates with diffusion zones is presented and used in SANS experiments. For diluted and dense packed systems structure describing parameter values were obtained. These results lead to a deeper understanding of the process of nanocristallization of amorphous alloys. The investigation of SANS on Fe73.5Si15.5B7Cu1Nb3 shows that the Fe3Si type nanocrystals created in the amorphous matrix during annealing are covered by Nb-atoms. The accumulation of Nb-atoms or Nb-B-aggregates acting as inhibitors at the surface of the nanocrystals is assumed to be the basic mechanism controlling the evolution of the precipitates. For the first time this inhibitor-model is shown to be correct without doubts. In the Zr32Ti7.5Al10Cu20Ni8 amorphous alloy the formation of ultrafine nanocystals of about 2-3 nm in diameter was observed. The nanocrystallization starts after ordered clusters achieved particular sizes and a certain packing fraction. This leads to a new model for the microscopic formation procedure of ultrafine nanocrystals in this amorphous alloy. Theoretical models of fractal systems are applied to complicated polydisperse materials. Both the theory for an exact surface fractal of Hermann (1994)and the model for coupled volume and surface fractals in the formulation of Wong (1992) are shown to be applicable. The latter approach is applied to experimental data here for the first time. With computer simulations conditions for scattering experiments were optained therewith predictions about the quality and grade of fractality in real specimens become possible. / Die vorliegende Arbeit ist ein Beitrag zur Strukturaufklärung mittels Neutronen-Kleinwinkel-Streuung (SANS). Es wird eine neu entwickelte analytische Streufunktion für polydisperse Ausscheidungen mit Diffusionszonen genutzt, um SANS Experimente auszuwerten. Sowohl für verdünnte, als auch für dicht gepackte Systeme werden auf diese Weise quantitative Strukturparameter gewonnen. Diese liefern einen Beitrag zum Verständnis des Nanokristallisationsverhaltens amorpher metallischer Gläser. Die Auswertung der Experimente an on Fe73.5Si15.5B7Cu1Nb3 zeigt, dass Fe3Si-artige Nanokristalle, die während der Temperaturbehandlung in der amorphen Matrix entstehen, von Nb-Atomen bedeckt werden. Diese Ansammlung von Nb-Atomen oder von entsprechenden Nb-B-Aggregaten auf der Oberfläche dieser Ausscheidungen hemmt das Größenwachstum der entstehenden Nanokristalle. Dieses Inhibitor-Modell wurde hier erstmals zweifelsfrei bestätigt. In Proben des amorphen metallischen Glases Zr32Ti7.5Al10Cu20Ni8 werden ultrafeine Ausscheidungen mit Durchmessern von 2-3 nm beobachtet. Diese entstehen verzögert nach der Ausprägung dicht gepackter Gebiete mit erhöhter Nahordnungsstruktur. Es wird ein Modell vorgeschlagen, das diesen Prozess erklären kann. Theoretisch diskutierte Modelle für fraktale Systeme werden auf komplizierte polydisperse Materialien angewendet. Sowohl die Formulierung von Hermann (1994) für ein exaktes Oberflächenfraktal, als auch der erstmals auf experimentelle Daten angewendete Ansatz von Wong (1992) für ein gekoppeltes Volumen- und Oberflächenfraktal erweisen sich als praktisch nutzbar. Mittels Computersimulationen wurden Bedingungen abgeleitet, die an Streuexperimente zu stellen sind, damit Aussagen über Qualität und Grad von Fraktalität in realen Proben getroffen werden können.

Die Entdeckung eines geladenen Higgs-Bosons, H+, wäre ein unbestreitbarer Nachweis von Physik jenseits des Standardmodells. In der vorliegenden Arbeit wird die Suche nach dem H+ mit Hilfe von Proton-Proton-Kollisionen, welche im Jahr 2011 mit dem ATLAS Experiment am Large Hadron Collider, LHC, des CERN aufgenommen wurden, beschrieben. Im Rahmen dieser Arbeit wurde eine überarbeitete Analyse der Suche nach geladenen Higgs-Bosonen, die eine Verhältnismethode anwendet und damit die Sensitivität des traditionell direkten Suchansatzes stark verbessert, durchgeführt. Leichte geladene Higgs-Bosonen, welche eine Masse geringer als die des Top-Quarks aufweisen, können aus einem Top-Quark-Zerfall hervorgehen. Im Gegensatz zu den schweren geladenen Higgs-Bosonen sind die leichten aufgrund des hohen Produktionswirkungsquerschnitts von Top-Quark-Paaren am LHC potenziell mit den ersten Daten des Experiments beobachtbar. In den meisten Theorien und Szenarien sowie dem größten Bereich ihres Phasenraumes zerfallen leichte geladene Higgs-Bosonen meist im H± → τ±ν Kanal. Demzufolge spielen sowohl die τ-Identifikation als auch die τ-Fehlidentifikation eine besondere Rolle für die Suche nach geladenen Higgs-Bosonen. Eigens für die Ermittlung der Fehlidentifikationswahrscheinlichkeiten von Elektronen als hadronisch zerfallende τ-Leptonen wurde eine “tag-and-probe”-Methode, basierend auf Z → ee Ereignissen, entwickelt. Diese Messungen sind mit den allerersten Daten durchgeführt worden. Dabei haben diese einerseits für alle Analysen, welche die Elektronenveto-Algorithmen der τ-Identifikation nutzen, essenzielle Skalenfaktoren hervorgebracht. Andererseits wurde, beruhend auf diesen Ergebnissen, eine datenbasierte Abschätzungsmethode entwickelt und für die Untergründe der geladenen Higgs-Boson-Suche, die von der Fehlidentifikation von Elektronen als hadronisch zerfallende τ-Leptonen stammen, erfolgreich implementiert. Im Rahmen dieser Arbeit wurden Triggerstudien, mit dem Ziel höchstmögliche Signaleffizienzen zu gewährleisten, durchgeführt. Neuartige Triggerobjekte, basierend auf einer Kombination aus τ-Trigger und fehlender transversaler Energie-Trigger, wurden entworfen, überprüft und in das Triggermenü für die Datennahme im Jahr 2012 aufgenommen. Eine direkte Suche nach dem geladenen Higgs-Boson wurde in drei Kanälen mit einem τ-Lepton im Endzustand unter Berücksichtigung des gesamten Datensatzes des Jahres 2011 durchgeführt. Da kein signifikanter Überschuss, der von den Vorhersagen des Standardmodells abweicht, in den Daten beobachtet wurde, sind obere Ausschlussgrenzen auf B(t → bH+) gesetzt worden. Letztlich ist die Analyse des Kanals mit einem hadronisch zerfallenden τ-Lepton und einem Myon oder Elektron im Endzustand des tt ̄-Zerfalls, unter Anwendung der sogenannten Verhältnismethode, wiederholt worden. Diese Methode misst Verhältnisse von Ereignisausbeuten, anstatt die Verteilungen diskriminierender Variablen zu evaluieren. Folglich kürzen sich die meisten dominant beitragenden systematischen Unsicherheiten intrinsisch heraus. Die Daten stimmen mit den Vorhersagen des Standardmodells überein. Durch Zuhilfenahme der Verhältnismethode wurden die oberen Ausschlussgrenzen, im Vergleich zur direkten Suche, signifikant verbessert. Die Resultate der Verhältnismethode sind mit denen der direkten Suche, welche ein hadronisch zerfallendes τ-Lepton und zwei Jets im Endzustand des tt ̄-Zerfalls nutzt, kombiniert worden. Auf diese Art und Weise konnten obere Ausschlussgrenzen auf B(t → bH+) in einem Bereich von 0,8 %–3,4 % für geladene Higgs-Bosonen in einem Massenbereich für m_H+ zwischen 90 GeV und 160 GeV gesetzt werden. Sollte das Minimal Supersymmetrische Standardmodell (MSSM) in der Natur realisiert sein, so haben die hier ermittelten oberen Ausschlussgrenzen auf B(t → bH+) direkte Konsequenzen für die Identität des Higgs-Boson-ähnlichen Teilchens, welches im Jahr 2012 am LHC entdeckt wurde.:Kurzfassung v Abstract vii Contents ix 1 Introduction 1 2 Theoretical Framework 3 2.1 The Standard Model of Particle Physics 3 2.1.1 Particles, Fields and Interactions 3 2.1.2 Gauge Theory in a Nutshell 6 2.1.3 Brout-Englert-Higgs Mechanism 9 2.2 Supersymmetry 11 2.2.1 Sources of Supersymmetry Breaking 20 2.2.2 Two Higgs Doublet Model 21 2.2.3 Charged Higgs Boson Production and Decay 22 2.3 Current Status of charged Higgs Boson Searches 24 3 Monte Carlo Simulation 33 3.1 Methodology of Monte Carlo Simulation 33 3.2 Monte Carlo Simulation for Electron to t Mis-identification Analysis 35 3.3 Monte Carlo Simulation for H± Analysis with Data taken 2010 35 3.4 Monte Carlo Simulation for H± Analyses with Data taken 2011 37 4 LHC and the ATLAS Detector 41 4.1 The Large Hadron Collider 41 4.2 The ATLAS Detector 42 4.2.1 Magnet Systems 46 4.2.2 Inner Detector 47 4.2.3 Calorimeter Systems 50 4.2.4 Muon Spectrometer 60 4.2.5 Forward Detectors 61 4.2.6 Trigger and Data Acquisition 63 4.3 Data Taking 67 5 Event Selection and Data-Driven Background Estimation Techniques 73 5.1 Event Cleaning 74 5.2 Trigger for the Charged Higgs Boson Analyses 74 5.2.1 Trigger for the thad+Lepton and tlep+Jets Channels 75 5.2.2 Trigger for the thad+Jets Channel 77 5.3 Physics Object Reconstruction 77 5.3.1 Muons 77 5.3.2 Electrons 77 5.3.3 Jets 78 5.3.4 b-Tagging 79 5.3.5 Hadronically decaying t Leptons 79 5.3.6 Missing Transverse Momentum 79 5.3.7 Removal of Geometric Overlaps between Objects 80 5.4 Selection and Cut Optimisation 80 5.4.1 thad+Lepton Analysis Selection 80 5.4.2 tlep+Jets Analysis Selection 82 5.4.3 thad+Jets Analysis Selection 83 5.5 Background Estimations 84 5.5.1 Measurement of the t Lepton Mis-identification Probability from Electrons 84 5.5.2 Backgrounds with Electrons and Jets Mis-identified as t Leptons 88 5.5.3 Embedding Method 92 5.5.4 Multi-Jet Background 92 5.5.5 Backgrounds with Mis-identified Leptons 93 6 Direct Searches for the Charged Higgs Boson 95 6.1 Analysis of the thad+Lepton Channel 95 6.2 Analysis of the tlep+Jets Channel 95 6.3 Analysis of the thad+Jets Channel 98 6.4 Systematic Uncertainties 101 6.4.1 Systematic Uncertainties of Detector Simulation 101 6.4.2 Systematic Uncertainties of Generation of tt Events 103 6.4.3 Systematic Uncertainties of Data-Driven Background Estimates 103 6.5 Results 107 7 Indirect Search for the Charged Higgs Boson – The Ratio-Method 111 7.1 Ratio-Method: A Re-Analysis of the thad+Lepton Channel 111 7.2 Event Selection 112 7.3 Measured Event Yield Ratios 115 7.4 Systematic Uncertainties 115 7.4.1 Systematic Uncertainties of Detector Simulation 115 7.4.2 Systematic Uncertainties of Generation of tt Events 115 7.4.3 Systematic Uncertainties of Data-Driven Background Estimates 118 7.5 Results 119 7.5.1 Upper Limits obtained from Results of the Ratio-Method 119 7.5.2 Combination of Upper Limits obtained from Direct Searches for Charged Higgs Bosons in the thad+jets final state and the Ratio-Method Results 124 8 Comparison and Discussion of the Results 131 9 Summary and Outlook 137 A Monte Carlo Simulation Samples 141 A.1 Monte Carlo Simulation Samples for 2010 Analyses 141 A.1.1 Monte Carlo Simulation Samples for the Electron to t Lepton Mis-identification Analysis 141 A.1.2 Monte Carlo Simulation Samples for the H± Analysis 143 A.2 Monte Carlo Simulation Samples for 2011 Analyses 147 A.2.1 Monte Carlo Samples for H± Ratio-Method Search Analysis 147 List of Figures 149 List of Tables 153 Bibliography 155 Acknowledgements 179

Mining will remain indispensable for the foreseeable future. For millennia, our society has been exploring and exploiting mineral deposits. Consequently, most of the easily exploitable high-grade deposits, which were of primary interest given their obvious technical and economic advantages, have already been depleted. For the future, the mining sector will have to efficiently produce metals and minerals from low-grade orebodies with complex mineralogical and microstructural properties -- these are generally referred to as complex orebodies. The exploitation of such complex orebodies carries significant technical risks. However, these risks may be reduced by applying modelling tools that are reliable and robust. In a broad sense, modelling techniques are already applied to estimate the resources and reserves contained in a deposit, and to evaluate the potential recovery (i.e., behaviour in comminution and separation processes) of these materials. This thesis focusses on the modelling of recovery processes, more specifically mineral separation processes, suited to complex ores. Despite recent developments in the fields of process mineralogy and geometallurgy, current mineral separation modelling methods do not fully incorporate the available information on ore complexity. While it is well known that the mineralogical and microstructural properties of individual particles control their process behaviour, currently widely applied modelling methods consider only distributions of bulk particle properties, which oftentimes require much simplification of the particle data available. Moreover, many of the methods used in industrial plant design and process modelling are based on the chemical composition of the samples, which is only a proxy for the mineralogical composition of the ores. A modelling method for mineral separation processes suited to complex ores should be particle-based, taking into consideration all quantifiable particle properties, and capable of estimating uncertainties. Moreover, to achieve a method generalizable to diverse mineral separation units (e.g., magnetic separation or flotation) with minimal human bias, strategies to independently weight the importance of different particle properties for the process(es) under investigation should be incorporated. This dissertation introduces a novel particle-based separation modelling method which fulfills these requirements. The core of the method consists of a least absolute shrinkage and selection operator-regularized (multinomial) logistic regression model trained with a balanced particle dataset. The required particle data are collected with scanning electron microscopy-based automated mineralogy systems. Ultimately, the method can quantify the recovery probability of individual particles, with minimal human input, considering the joint influence of particle shape, size, and modal and surface compositions, for any separation process. Three different case studies were modelled successfully using this new method, without the need for case-specific modifications: 1) the industrial recovery of pyrochlore from a carbonatite deposit with three froth flotation and one magnetic separation units, 2) the laboratory-scale magnetic separation of a complex skarn ore, and 3) the laboratory-scale separation of apatite from a sedimentary ore rich in carbonate minerals by flotation. Moreover, the generalization potential of the method was tested by predicting the process outcome of samples which had not been used in the model training phase, but came from the same geometallurgical domain of a specific ore deposit. In each of these cases, the method obtained high predictive accuracy. In addition to its predictive power, the new particle-based separation modelling method provides detailed insights into the influence of specific particle properties on processing behaviour. To name a couple, the influence of size on the recovery of different carbonate minerals by flotation in an industrial operation; and a comparison to traditional methodologies demonstrated the limitation of only considering particle liberation in process mineralogy studies -- the associated minerals should be evaluated, too. Finally, the potential application of the method to minimize the volume of test work required in metallurgical tests was showcased with a complex ore. The approach developed here provides a foundation for future developments, which can be used to optimize mineral separation processes based on particle properties. The opportunity exists to develop a similar approach to model the comminution of single particles and ultimately allow for the full prediction of the recovery potential of complex ores.:1 Introduction 1 1.1 Background . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 1.2 Objectives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 1.3 State-of-the-art in particle-based separation models . . . . . . . . . . . 11 1.4 Moving forward . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 1.4.1 Particle data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 1.4.2 Mathematical tools required for the particle-based separation model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 1.4.3 Workflow . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 1.5 Structure of the thesis . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 2 The method and its application to industrial operations 23 2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24 2.1.1 Background . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24 2.2 Methodology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 2.2.1 Assumptions and limitations . . . . . . . . . . . . . . . . . . . . 26 2.2.2 Data structure and required pre-treatment . . . . . . . . . . . . 27 2.2.3 Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 2.3 Demonstration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31 2.3.1 Artificial test cases . . . . . . . . . . . . . . . . . . . . . . . . . 31 2.3.2 Real case study . . . . . . . . . . . . . . . . . . . . . . . . . . . 35 2.4 Discussion and final considerations . . . . . . . . . . . . . . . . . . . . 39 3 The robustness of the method towards compositional variations of new feed samples 45 3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46 3.2 Generalization potential of current Particle-based Separation Model (PSM) methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 3.3 Case study . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 3.3.1 Samples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 3.3.2 Dry magnetic separation tests . . . . . . . . . . . . . . . . . . . 53 3.3.3 Sample characterization . . . . . . . . . . . . . . . . . . . . . . 53 3.3.4 Particle-based separation models . . . . . . . . . . . . . . . . . 54 3.4 Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56 3.5 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62 3.6 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63 4 Flotation kinetics of individual particles 67 4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68 4.2 Method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70 4.2.1 Data collection . . . . . . . . . . . . . . . . . . . . . . . . . . . 71 4.2.2 Cumulative recovery probability . . . . . . . . . . . . . . . . . . 72 4.2.3 Particle-based kinetic flotation model . . . . . . . . . . . . . . . 74 4.3 Demonstration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75 4.3.1 Materials and methods . . . . . . . . . . . . . . . . . . . . . . . 75 4.3.2 Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76 4.4 Discussion and final thoughts . . . . . . . . . . . . . . . . . . . . . . . 80 5 Conclusions and outlook 85 5.1 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85 5.2 Outlook . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86 Bibliography 89