Browsing by study line "Computational Material Physics"

Halide perovskites are a promising materials class for solar energy production. The photovoltaic efficiency of halide perovskites is remarkable but their toxicity and instability have prevented commercialization. These problems could be addressed through compositional engineering in the halide perovskite materials space but the number of different materials that would need to be considered is too large for conventional experimental and computational methods. Machine learning can be used to accelerate computations to the level that is required for this task. In this thesis I present a machine learning approach for compositional exploration and apply it to the composite halide perovskite CsPb(Cl, Br)3 . I used data from density functional theory (DFT) calculations to train a machine learning model based on kernel ridge regression with the many-body tensor representation for the atomic structure. The trained model was then applied to predict the decomposition energies of CsPb(Cl, Br)3 materials from their atomic structure. The main part of my work was to derive and implement gradients for the machine learning model to facilitate efficient structure optimization. I tested the machine learning model by comparing its decomposition energy predictions to DFT calculations. The prediction accuracy was under 0.12 meV per atom and the prediction time was five orders of magnitude faster than DFT. I also used the model to optimize CsPb(Cl, Br)3 structures. Reasonable structures were obtained, but the accuracy was qualitative. Analysis on the results of the structural optimizations exposed shortcomings in the approach, providing important insight for future improvements. Overall, this project makes a successful step towards the discovery of novel perovskite materials with designer properties for future solar cell applications.

Ion beams have been the subject of significant industry interest since the 1950s. They have gained usage in many fields for their ability to modify material properties in a controlled manner. Most important has been the application to semiconductor devices such as diodes and transistors, where the necessary doping is commonly achieved by irradiation with appropriate ions, allowing the development of the technology that we see in everyday use. With the ongoing transition to ever smaller semiconductor devices, the precision required of the manufacturing process correspondingly increases. A strong suite of modeling tools is therefore needed to advance the understanding and application of ion beam methods. The binary collision approximation (BCA) as a simulation tool was first introduced in the 1950s. It allows the prediction of many radiation-related phenomena for single collision cascades, and has been adopted in many experimental laboratories and industries due to its efficiency. However, it fails to describe chemical and thermodynamic effects, limiting its usefulness where ballistic effects are not a sufficient description. Parallel to BCA, the molecular dynamics (MD) simulation algorithm was developed. It allows a more accurate and precise description of interatomic forces and therefore chemical effects. It is, however, orders of magnitude slower than the BCA method. In this work, a new variant of the MD algorithm is developed to combine the advantages of both the MD and the BCA methods. The activation and deactivation of atoms involved in atomic cascades is introduced as a way to save computational effort, concentrating the performed computations in the region of interest around the cascade and ignoring surrounding equilibrium regions. By combining this algorithm with a speedup scheme limiting the number of necessary relaxation simulations, a speedup of one order of magnitude is reached for covalent materials such as Si and Ge, for which the algorithm was validated. The developed algorithm is used to explain the behavior of Ge nanowires under Xe ion irradiation. The nanowires were shown experimentally to bend towards or away from the ion beam, and computational simulations might help with the understanding of the underlying physical processes. In this thesis, the high-fluence irradiation of a Ge nanowire is simulated and the resulting defect structure analyzed to study the bending, doubling as a second test of the developed algorithm.

The rapidly increasing global energy demand has led to the necessity of finding sustainable alternatives for energy production. Fusion power is seen as a promising candidate for efficient and environmentally friendly energy production. One of the main challenges in the development of fusion power plants is finding suitable materials for the plasma-facing components in the fusion reactor. The plasma-facing components must endure extreme environments with high heat fluxes and exposure to highly energetic ions and neutral particles. So far the most promising materials for the plasma-facing components are tungsten (W) and tungsten-based alloys. A promising class of materials for the plasma-facing components is high-entropy alloys. Many high-entropy alloys have been shown to exhibit high resistance to radiation and other wanted properties for many industrial and high-energy applications. In materials research, both experimental and computational methods can be used to study the materials’ properties and characteristics. Computational methods can be either quantum mechanical calculations, that produce accurate results while being computationally extremely heavy, or more efficient atomistic simulations such as classical molecular dynamics simulations. In molecular dynamics simulations, interatomic potentials are used to describe the interactions between particles and are often analytical functions that can be fitted to the properties of the material. Instead of fixed functional forms, interatomic potentials based on machine learning methods have also been developed. One such framework is the Gaussian approximation potential, which uses Gaussian process regression to estimate the energies of the simulation system. In this thesis, the current state of fusion reactor development and the research of high-entropy alloys is presented and an overview of the interatomic potentials is given. Gaussian approximation potentials for WMoTa concentrated alloys are developed using different number of sparse training points. A detailed description of the training database is given and the potentials are validated. The developed potentials are shown to give physically reasonable results in terms of certain bulk and surface properties and could be used in atomistic simulations.

The plasma-facing materials of future fusion reactors are exposed to high doses of radiation. The characterization of the radiation damage is an essential part in the study of fusion relevant materi- als. Electron microscopy is one of the most important tools used for characterization of radiation damage, as it provides direct observations of the microstructure of materials. However, the char- acterization of defects from electron microscope images remains difficult. Simulated images can be used to bridge the gap between experimental results and models. In this thesis, scanning transmission electron microscope (STEM) images of radiation damage were simulated. Molecular dynamics simulations were employed in order to create defects in tungsten. STEM images were simulated based on the created systems using the multislice method. A data- base of images of h001i dislocation loops and defects produced from collision cascade simulations was generated. The simulated images provide insight into the observed contrast of the defect structures. Differences in the image contrast between vacancy and interstitial h001i dislocation loops were reported. In addition to this, the results were compared against experimental images and used in identification of a dislocation loop. The simulated images demonstrate that it is feasible to simulate STEM images of radiation damage produced from collision cascade simulations.