Opportunities for B.Sc., M.Sc. or Ph.D. theses

If you are interested in doing a Bachelor, Master or PhD thesis with us, please do not hesitate to inquire either via email to Prof. Karsten Reuter or by stopping by at the secretary of the Chair of Theoretical Chemistry. Topics for Bachelor and Master theses are in general not individually announced, but are continuously available in one of our research focus areas. Please consult our research pages for more information.

At present we specifically invite applications for the following projects (other projects can be assigned on personal interest and availability).



Ph.D. projects

I) PhD position: Multi-scale modelling approach to the growth of 2D materials

The goal of the project is to explore a novel route to the fabrication of two-dimensional (2D) materials such as graphene, hexagonal boron nitride, silicene, and germanane involving gas-phase precursors reacting on a molten catalyst surface. The successful candidate will apply a multi-scale modelling approach. The elementary processes at the atomic scale will be studied using density functional theory (DFT) or its tight-binding version (DFTB). The obtained mechanistic insight will then serve as input for kinetic Monte Carlo simulations of the growth process. The project is part of a larger European collaboration involving experimental researchers from the Netherlands, France and Greece.

A background in chemistry, physics, materials science or similar is required. Prior knowledge and experience in electronic structure theory, UNIX based operating systems and a scripting language (e.g. Python) would be beneficial, but is not necessary. Evaluation of applications will begin immediately. The starting date should preferably be in the fall of 2017 and no later than Jan 1st 2018.

Enquiries and applications including CV, motivation letter, recent transcripts and recommendation letters (if available) should be directed to Dr. Mie Andersen (mie.andersen .at. ch.tum.de).


M.Sc. projects

I) Implicit solvation studies using the Multipole Expansion model

One of the major goals in present energy science is to relieve society from its strong dependancy on limited and unrenewable energy sources, such as coal or oil. In this context, sunlight is very promising since abundantly available and easily accessible. However, for most applications photoenergy needs to be converted into chemical energy that can be stored, e.g. in batteries, by chemical reduction of carbon dioxide, or by splitting water into hydrogen and oxygen. In our project we want to gain a mechanistical understanding of the latter reaction on a titanium dioxide surface which has been observed experimentally.

Previous first principles studies commonly neglected effects of the surrounding water molecules to avoid costly ab-initio molecular dynamics simulations. Faced with possibly charged intermediate states that are readily stabilized by a polar solvent like water, this approximation seems questionable. The interested student will investigate solvation effects in suitable, simplified model reactions using an implicit solvation scheme based on the Multipole Expansion (MPE) method in the framework of Density Functional Theory (DFT). Existing knowledge on UNIX based operating systems and programming is desirable, but not required.


II) Tackling complexity in multiscale kinetic simulations: Kinetic Monte Carlo software performance analysis. 

The kinetic Monte Carlo (kMC) method is ideal for tackling problems that require atomic scale detail but whose extent is beyond capabilities of state of the art quantum chemistry methods. In particular, it is routinely used for problems in crystal growth, heterogeneous catalysis, and solid diffusion, among others. Currently, one the most pressing challenge for the advancement of the method resides in the need for a more efficient treatment of problems of high complexity. In our group we are currently working towards solving these problems in two main areas of research: reactions in surfaces (heterogeneous catalysis) and lithium-ion diffusion in LTO (battery materials). In both cases, complexity becomes an issue when the microscopic rates (of diffusion or other elementary reaction events) are dependent on the instantaneous state of the system (lateral interaction effects). For heterogeneous catalysis, this problem is (combinatorially) exacerbated when the number of reactants considered is high, like in technological reactions. In the case of battery materials, typical crystal structures allow for a wide variety of local configurations leading to strong effects on microscopic diffusion rates.

The main tool to be used for this project is the `kmos' kMC framework; an extensible, modular software package, actively developed and used in our group. Recently, an extension has been implemented to better deal with the class of problems mentioned above. The interested student will be expected to familiarize him or herself with the concept of kMC simulations and with the `kmos' package and to implement and run a collection of kMC models of varying levels of complexity. The models used will include some developed in the Reuter group and some obtained from the scientific literature. The performance data obtained will be used to identify bottlenecks and guide the further development of `kmos'. Basic knowledge of UNIX based operating systems and some experience with programming (and/or scripting) languages (preferably Python) is desirable.


III) Surface morphology in complex reaction environments: 3D ab initio thermodynamics phase diagrams. 

Ab initio thermodynamics represent a solid theoretical foundation through which “T = 0K” results from electronic structure calculations can be used to determine the equilibrium phase diagram of surfaces in a wide range of pressures and temperatures. In the context of heterogeneous catalysis, this provides an extensive amount of valuable information for the subsequent analysis of reactivity in realistic (finite pressure and temperature) conditions. This project will focus on the study of the Pd(100) surface. Palladium is widely used for automotive catalytic converters and represents a cheaper alternative to more widespread platinum. The work here proposed will build upon previous studies, in which phase diagrams for this surface under two simplified gas mixtures have been constructed: either under mixtures of CO and O2 or mixtures of NO and O2.

For this project, the interested student will construct an ab initio thermodynamics phase diagram for the Pd(100) surface under exposure to a mixture of O2, CO and NO gases and evaluate the extent of thermodynamic stability of coadsorption phases containing the three adsorbates, as well as the conditions for oxide formation under such complex gas mixtures. This work will complement our recent theoretical prediction suggesting that it is not possible to infer kinetic behavior in complex gas mixture from the separate analysis of simplified gas mixtures. In this project, the student will learn how to perform total energy calculation using the modern density functional theory (DFT) code CASTEP. Existing knowledge of UNIX based operating systems and basic scripting is desirable, but not required. 


IV) Higher alcohol synthesis from syngas on metal catalysts

Higher alcohols are attractive fuel additives since they can be blended directly into gasoline. At present, the sustainable production of ethanol from biomass occurs primarily via fermentation. The synthesis of higher alcohols from biomass­derived syngas (mainly CO and H2) represents an interesting alternative, since it allows for the use of a wider range of biomass sources. Still, the lack of efficient catalysts limits the industrial use of this approach.

The focus of the present project is on the catalytic properties of extended metal surfaces for ethanol synthesis. Experimentally, Rh catalysts show some selectivity and activity towards higher alcohol synthesis, though often in modified forms. With the use of computational modeling we aim for an improved understanding of the reaction mechanism on this type of catalysts. In the proposed M.Sc. project, the student will use density functional theory (DFT) to investigate key reaction steps such as C­C coupling steps on different metal surfaces and at different active sites. Also the effect of the local environment in terms of adsorbate­adsorbate interactions will be investigated. The obtained DFT calculations will be used to establish energetic trends in the form of scaling relations and in the end provide the input for microkinetic models such as mean­field and kinetic Monte Carlo simulations. Prior experience with UNIX based operating systems and a scripting language (e.g. Python) is advantageous, but not a prerequisite.

V) Germanium Intercalation of Graphene on Silicon Carbide: Can We Understand the Interface?

Germanium intercalated quasi-freestanding monolayer graphene (Ge-QFMLG) grown on SiC has recently been proven to be an ideal material to form ballistic graphene pn-junctions because the graphene layer doping level depends on the Ge layer thickness intercalated between the SiC substrate and the graphene sheet. The p-type doped interface contains double as many intercalated Ge atoms than the n-doped structure. So far all attempts to identify the interface structure were based on density-functional theory (DFT) using smaller approximated unit cells (3 SiC cells). However, for a reliable structure prediction of the two different interfaces the large (6sqrt3x6sqrt3)R30 unit cell (108 SiC and 169 graphene unit cells) is necessary. To correctly address the van der Waals (vdW) bonded graphene layer state-of-the-art dispersion corrections have to be included. A structure search on the bases of DFT for structure sizes with hundreds of atoms is currently computationally too demanding. A promising compromise between computational cost and accuracy is the Density Functional based Tight Binding (DFTB) method incorporating a state-of-the-art vdW correction scheme.

The proposed MSc thesis is part of a collaborative research project with the group of Prof. Dr. F. S. Tautz (Forschungszentrum Juelich). We will perform a multilevel structure search to address the challenges of predicting an interface structure for very large systems. First, under consideration of the substrate symmetry, we will generate trial interface structures. Then, for these structures, we will construct a surface phase diagram using the semi-empiric DFTB method including dispersion corrections to identify the most promising candiates. The search will then be refined by recalculating the lowest energy structure candidates in DFT. In the process, the MSc student will deepen his/her theoretical background in the interdisciplinary field of materials science and become familiar with state-of-the-art computational methods and modeling techniques. Existing knowledge on UNIX based operating systems and basic programming skills (preferably in Python) is desirable.


B.Sc. projects

I) Photo-catalytic CO2 reduction

Among other things humanity is facing two problems in the not too far future, global warming due to an overabundance of greenhouse gasses such as CO2 in the atmosphere and the need to transition away from finite fossil fuels to renewable energies. A way to tackle both these problems presents itself in the form of photochemical carbon dioxide reduction, where sunlight is harnessed to convert CO2 to higher energy products such as methanol. Unfortunately, currently available CO2 photo-reduction catalysts all suffer from very low turnover rates and are therefore not suitable for large scale industrial application. In the proposed Bachelor's project, candidates will employ state of the art computational---such as embedded density functional theory/classical mechanics simulations---and theoretical methods developed in the Reuter group to identify efficiency limiting steps and search for more effective photo-catalysts. Existing knowledge of UNIX based operating systems and programming is desirable, but not a necessity.


II) Battery materials: LTO for LIBs

Rechargeable lithium-ion batteries (LIBs) are key components of today's technology that power a range of devices from mobile phones to electric vehicles. During the discharge/charge cycles the electrode materials of a LIB take up or release Li ions and electrons, thereby undergoing changes in chemical composition, structure and electronic structure. In case of lithium-transition metal-oxides, which represent an important class of LIB electrode materials, reversible (de)intercalation of lithium ions in their structures can take place, and the redox-reaction involves a change in the oxidation state of the transition metal atoms. The focus of our project is on lithium-titanium-oxide materials, especially on Li4Ti5O12 (LTO), which can be used as an anode material for LIBs. The computational studies aim to contribute to a better atomic scale understanding of materials properties, such as Li ion mobility and electronic conductivity. Key issues to be dealt with include structural disorder and the occurrence of different oxidation states of the titanium atoms as the Li content is varied. In the BSc project you will use density functional theory techniques to study relations between the chemical composition, structural features and the electronic structure of Li-Ti-O materials. Prior experience with Linux and a scripting language is advantageous but not a prerequisite. 


III) Bistability in NO oxidation reactions in Pd(100)

Palladium is a very interesting material for automotive catalytic converters and their properties have been extensively studied. However, the actual composition of the surface during operation, whether it is the pristine metal termination, a monolayer surface oxide or a bulk oxide layer, is still unknown. Recently, theoretical and experimental analysis on conditions of CO oxidation have suggested that actually two distinct terminations (pristine metal and surface oxide) might be present simultaneously, i.e. the surface is bistable. Recently, we have developed first-principles kinetic Monte Carlo models to analyze the NO oxidation properties of Pd(100), a reaction of interest for catalytic converters for lean-burn and diesel engines. In this project, the student will use these models to analyze the (relative) kinetic stability of the monolayer surface oxide and the pristine metal surface for different conditions of NO and O_{2} partial pressures, with the aim of determining whether a bistability region analogous to the one observed for CO oxidation exists also for NO oxidation.


IV) A chemical puzzle: How the structure of a molecule influences the polarizability

Despite the rapid advancement in computing power and methods, the calculation of large chemical systems remains a challenge for modern ab initio electronic structure theory. To investigate and understand processes in proteins, polymers or covalent organic frameworks, empirical force field methods have been used with some success. A main shortcoming of existing force fields is the ability to describe the polarization of the electron density with respect to changes in the environment. To overcome this problem, a new method to describe the polarizability is developed in our group (titled 'enhanced ACKS2'). One of the remaining obstacles is a property called 'response matrix', which depends in an unknown way on the geometry of the molecule, making it necessary to perform expensive computations that countermand the efficiency of the force field approach. In the proposed Bachelor's project the candidate will employ state of the art computational methods to perform a systematic investigation of the response matrix for a wide range of chemically diverse molecules and molecular fragments. This data is then analyzed to shed light on the structure-property relationship and to predict response matrices for different molecules, eventually employing machine learning methods. Existing knowledge of UNIX based operating systems and programming (Python) is highly desirable, but not a necessity.