- Greta Lupi (2023)
Superconductivity in Multilayer Molybdenum Disulfide: the Role of Long-Range Coulomb Interaction (pdf, 1,3 MB)
Superconductivity in n-doped molybdenum disulfide (𝑀𝑜𝑆2) is heavily dependent to the number of layers, as established in previous studies. A striking feature, observed experimentally, is the sudden enhancement of critical temperature when transitioning from a monolayer to a bilayer configuration, a phenomenon also common in various other transition metal dichalcogenides, e.g., 𝑊𝑆2 and 𝑁𝑏𝑆𝑒2. The influence of specific Fermi surface pockets, termed 𝑄 pockets, in the Brillouin zone, is conjectured to further enhance superconductivity.
This study aims to investigate the effects of Coulomb interaction on superconducting properties when an additional layer is introduced to a single layer of 𝑀𝑜𝑆2. The central objective is to calculate the Tolmachev-Morel-Anderson pseudopotential 𝜇∗, in an endeavor to derive a multi-valley parameter and assess the contribution of the Q pockets. Furthermore, EELS spectra are calculated and analyzed to get an insight on plasmonic properties, which might affect superconductivity. To accomplish this, the long-range Coulomb interaction is calculated from first principles for monolayer and bilayer 𝑀𝑜𝑆2 at varied doping levels. Specifically, Density Functional Theory calculations are conducted in combination with constrained Random Phase Approximation calculations and extended to estimate the Coulomb interaction contribution to superconductivity, via a consistent treatment of the screening channels.
Gerco Ganzevoort (July 2023)
Towards plasmarons: modelling the non-local Coulomb interactions in (bilayer) graphene (pdf, 11 MB)
When electrons interact via the Coulomb interaction, a composite particle called the plasmon emerges. This particle can subsequently couple to bare electrons, forming a dressed electron which is known as a plasmaron. In this thesis, the energy-momentum relation of this dressed electron is calculated for electrons in monolayer and bilayer graphene using G0W0 theory. In order to do so, the Coulomb interaction between various pz-orbitals in the graphene layers is derived first, taking into account their shape. Non-local creening stemming from the presence of a substrate is taken into account as well to make the model applicable for any insulating substrate. It is found that, at the K-point, the two bands are visible whereas a noninteracting system shows only a single band. This finding is in agreement with other literature and is the cause of the apparent spectral weight transfer (also known as mass enhancement). Further research can focus on the role of these plasmons on the superconducting state of the graphene layer(s).
- Mika van Zijl (March 2023)
Graphene in a Molecular Liquid & under Uniaxial Strain (pdf, 8,8 MB)
Atomistic simulations have shown that the statistical mechanical properties of graphene, a prototype of a 2D crystal, are well described by the phenomenological theories of flexible membranes. In many realistic applications, graphene is immersed in a liquid. Atomistic simulations of graphene, however, have only been performed for graphene in vacuum and the effect that a solvent might have on the structural an thermal properties is not yet fully understood. In the present work we present NPT Molecular Dynamics (MD) simulations of graphene immersed in a benzene liquid based on LCHBOP. A negative Poisson ratio (PR) is one of the few predictions from the phenomenological membrane theory that has not yet been confirmed for graphene from atomistic simulations. Recently proposed detailed theoretical predictions regarding the PR of a crystalline membrane require verification from accurate numerical simulations. We find from NVT MD simulations that both the absolute and differential PR ratio are strictly positive for graphene in the limit of small stress.
- Wouter Bos (February 2023)
An Ab Initio Study of Superconductivity in Mono- and bilayer Molybdenum disulfide (pdf, 6,5 MB)
The two-dimensional electron doped form of molybdenumdisulfide (MoS2) exhibits superconducting properties around liquid-helium temperatures. Ion-liquid gating experiments have shown a layer dependent increase in critical temperature from 1-2 Kelvin (K) to around 7 K in the transition from monolayer to bilayer MoS2 converging to 10 K in bulk MoS2 [1, 2, 3]. In this thesis an attempt has been made to increase the understanding of conventional superconductivity in doped monolayer and bilayer MoS2 by using state of the art ab initio techniques such as, DFT, DFPT, Wannier interpolation and Eliashberg theory. These techniques were ultimately used to predict trends in critical temperatures differentiating by layer thickness, spin polarization and charge distribution. The valley and band resolved electron-phonon coupling strengths were calculated using a self-written Python script, which showed electron occupation of the Q-region to be crucial to superconductivity in bilayer MoS2. The role of spin polarized electronic states originating from asymmetric spin orbit coupling in high doping monolayer MoS2 was also examined. All scattering channels requiring spin flipping were found to diminish due to the prohibited nature of this process in electron-phonon interactions. Additionally, spin orbit coupling increases coupling to an optical mode at the M high symmetry point. A heterogenous charge distribution by applying an out-of-plane electric field induces spin-polarization in bilayer MoS2 as the global inversion symmetry is broken. The DFT calculations showed Q valley occupation might be reached at a lower charge carrier concentration, which could explain the earlier onset of superconductivity compared to homogeneously doped samples observed in experiments .
- Matthijs Vernooij (August 2021)
Skyrmions are specific patterns in the magnetisation direction of a ferromagnet. They are the subject of a lot of active research and a lot of ideas for their applications exist. Many of these applications require an efficient mechanism to measure the number of skyrmions. The topological Hall effect, which means that the Hall resistance is proportional to the number of skyrmions, is a candidate for this mechanism. In this thesis we theoretically calculate the resistivity for a 2D Rashba ferromagnet. We show that the Hall resistance is indeed proportional to the number of skyrmions, proving that the topological Hall effect is present in this model. This improves existing descriptions of the topological Hall effect by incorporating spin-orbit coupling via the Rashba effect, which had not been studied before.
- Yann in 't Veld (October 2020)
Plasmons in two-dimensional materials have a gapless square-root-like dispersion which can be eﬀiciently tuned by external dielectric screening or charge doping. At the same time, plasmons can couple strongly to electrons, which renders them particularly interesting for superconducting properties in layered materials. In this thesis we use an extended version of density functional theory for superconductors (SC-DFT) to account for both the full dynamic Coulomb interaction and phonon contributions in two dimensions. We find that interplaying phonons strongly enhance the superconducting critical temperature, especially at small interaction strengths. The interplay between phonons and plasmons shows both suppression and enhancement depending on the doping level. These results are promising in the view of exciton mediated superconductivity in sandwich structures and in various layered materials like transition metal dichalcogenides and twisted bilayer graphene. In all of these materials plasmonic contributions have so far been neglected.
- Tom Westerhout (August 2020)
Recent advances in experiments with cold atoms confined to wire geometries using strong optical lattices allow for a realization of tunable strongly interacting bosonic and fermionic models. In such experiments probably the most distinguishing feature of 1D fermionic systems can be observed, namely, the Luttinger liquid behavior. t-V-V' model is a very popular model hosting a Luttinger liquid phase, and even though it has been studied extensively, there remains some controversy. Around 2000, it was suggested that the Luttinger liquid phase breaks for large V = 2V', but more recent studies did not confirm this statement. Here, using VUMPS, a modern tensor-network algorithm which works directly in the thermodynamic limit, we thoroughly analyze the t-V-V' model and resolve the contradiction.
- Matteo Vandelli (October 2018)
Three-Particle Response Functions in Condensed Matter Physics: Application to Optical Second Harmonic Generation and Chiral Three-Spin Interaction (pdf, 2,9 MB)Starting from the expression for the partition function of a quantum electronic system, three-particle correlation functions are derived and analyzed the in the framework of a quantum field theory in condensed matter physics. Corresponding threepoint diagrams are calculated in the Matsubara scheme using the formalism of Green’s functions at finite temperatures. Afterwards, obtained results are applied to the following physical problems of interest. The first problem is related to optical properties of two-dimensional Dirac-like systems and consists of the calculation of the optical second harmonic response from these materials. The latter is presented in systems with the broken inversion symmetry and can be revealed, for example, in the Graphene model with the nonzero mass term. This situation is relevant, for instance, when the Graphene is disposed on top of an insulating substrate with the honeycomb lattice structure. Here, all components of the second harmonic response function have been calculated for the model of Graphene on the hexagonal boron nitride substrate, and their relation to the crystal symmetry of the obtained system has been investigated. The second problem concerns the study of the three-spin interaction that appears in the system of localized spins coupled to itinerant electrons. Here, the general expression for the threeparticle correction to the energy of the system is explicitly calculated in the framework of the s-d model. As the result, the relevant three-spin interaction term that enters an effective spin model has been derived. It has been found that the chiral term corresponding to the triple product of spins emerges in the presence of a magnetic flux. Other terms that arise in the presence of strong electronspin interaction or spin-orbit coupling have also been investigated. Furthermore, a continuous limit of the spin model in the presence of the chiral three-spin interaction (micromagnetic model) has been developed, starting a preliminary study of the stability of topological spin textures that can appear in the system due to the novel twisted-exchange interaction derived in this work.
- Claudius Müller (June 2018)
In this project the TSMs under investigation are the Dirac nodal line semimetals ZrSiS and HfSiS. The goal is to thoroughly study both materials experimentally and theoretically by means of magnetic quantum oscillations and DFT calculations. The magnetic quantum oscillation studies were performed by measuring the Shubnikov - de Haas (SdH) effect and the de Haas - van Alphen (dHvA) effect by means of magnetotransport and capacitive magnetotometry experiments up to 35T at the High Field Magnet Laboratory (HFML). For both materials a complex oscillatory behaviour of the SdH and dHvA oscillations is observed and compared to the bandsturcture obtained by DFT calculations. By Fourier analyzing the experimental data it becomes evident that the complex oscillatory behaviour is due to magnetic breakdown (MB) occuring in these two materials. Also, for ZrSiS a discrepancy between the frequency spectrum of the SdH and dHvA studies is observed.
- Giacomo Sesti (October 2017)
Time-Evolution of Magnetic Vector Fields through Disorder-Averaged Green Functions in Layered Antiferromagnets (pdf, 1,3 MB)
- Davide Tisi (September 2017)
In my master project I developed a code to compute, numerically, the phonons of graphene and a coefficient called bending rigidity which is a measure of the energy needed to bend the graphene layer. It has been demonstrated that, near the Γ point, the two in-plane acoustic modes of graphene exhibit a linear dispersion, while the out-of-plane acoustic mode has a quadratic dispersion. From this quadratic dispersion we can compute directly the bending rigidity. Since it has been shown that the usual quasi-harmonic approximation does not work for graphene due to the large anharmonicity, the phonons will be obtained from the peaks of the power spectral density. This will be computed as the temporal Fourier Transform of the k-space velocity autocorrelation function (kVACS). To calculate kVACS, we need the positions and the velocities of all particles from a molecular dynamics simulation, using an existing code. This code uses an empirical potential, called long-range Carbon Bond-Order Potential, to estimate the interactions between the carbon atoms. Previous numerical studies suggested that the temperature dependence of the bending rigidity changes, using a different torsion term in the potential of the molecular dynamics code. Therefore, I have implemented in the molecular dynamics code a new torsion term, and calculated the phonon spectrum with both models of the torsion.
- Matteo Stifano (May 2017)
- Rob Ouwersloot (December 2016)
The currently most promising methods to produce graphene in large amounts tend to create graphene with relatively small grain sizes. Small grains mean many grains and many grain boundaries. While pristine graphene has been studied and characterized extensively both theoretically and experimentally, little is known about the properties of graphene with grain boundaries. To improve the understanding of the effect of grain boundaries on the elastic properties of graphene, we have done simulations on several model systems containing specific grain boundaries. Earlier work found the lowest energy structures of -tilt grain boundary structures with varying misorientation angles. Here we take three of those structures as a starting point and go on to characterize several of their elastic properties, both for graphene with grain boundaries and for carbon nanotubes with grain boundaries.
Han van der Pluijm (August 2016)
Since the birth of quantum mechanics the theory has proven to generate reliable experimental predictions. But the status of the wave function is still at question: is it a real property of nature or is it the knowledge one has of nature? If we consider this question in an ontological model then the wave function would be called Ψ-ontic and Ψ-epistemic respectively. Here an ontological model is a model that assumes there exist ontic states with fixed properties. In this thesis the Pussy-Barrett-Rudolph no-go theorem is introduced which is an argument in favour of the Ψ-ontic view of quantum states. Furthermore, Spekkens' toy theory, an example of a Ψ-epistemic model is discussed. With this toy model it is possible to reproduce part of the phenomenology that is typical for quantum theory and therefore it motivates the Ψ-epistemic view of the wave function. Finally, a restricted description of the Mach-Zehnder interferometer is given in terms of Spekkens' toy theory.
- Robert Sokolewicz (May 2016)
 C. W. J. Beenakker, Phys. Rev. Lett. 112, 70604 (2014)
- Remco Hens (August 2015)
In this thesis we study two topics related to graphene. The first topic is about the dimer in an asymmetric periodic potential. We introduce it as a simple model to study the friction of graphene flakes on a graphite substrate. The second topic is about the mechanical deformation properties of graphene. We examine the differences between 1D (longstrips) and 2D graphene.
- Nicole Orval (July 2015)
In this thesis we study a spin lattice where a domain wall may behave as a soliton. The domain wall then consists of two straight parts with a kink in the middle. We look at the dynamics of domain wall kinks in the presence of an external magnetic field. We are interested in the nontrivial relation between the speed of the kinks and the momentum which is applied by an external magnetic field.
- Annelot Schuring (August 2014)
We present a theoretical/computational study of moir´e patterns in graphene on hexagonal substrates, in particular graphene on graphene and graphene on hexagonal boron nitride (hBN). We first consider rotated graphene on graphene systems. By varying the angle of rotation between the two layers, one can create beautiful super periodicity: moir´e patterns. One period of this hexagonal pattern is called a moiron. Our ambition is to answer this main question: How does minimization of the energy affect the moir´e patterns?
- Erik van Loon (October 2013)
Collective Excitations in Strongly Correlated Systems -- Charge conservation in the dual boson approach (pdf, 3,5 MB)This thesis contains the first step towards implentation and application of the dual boson approach. Quantum field theory methods are used to study the susceptibility of the (extended) Hubbard model. This model describes the electrons in a strongly correlated material. Depending on the choice of parameters, such a system can be metallic, insulating or even undergo a charge-ordering transition. This thesis shows how the charge susceptibility changes in these phases.The further developments in this project can be found here.
- Maria Patelkou (August 2013)
Pushing and Sliding allotropes of Carbon (two files due to pdf size)first pages (1 to 77)
last pages (70 to 107) (pdf, 13 MB)The focus of this master thesis lies on sliding friction and wear. What is the effect of pushing solid surfaces together and moving them with respect to each other? How can the effect be analyzed? We suggest a molecular dynamics simulation of pushing and sliding different allotropes of carbon bulk structures using the long-range order carbon bond order potential (LCBOP) and analyze the formation and the structure of the interface.
- Lennert van Tilburg (May 2013)
- Linde van Heeringen (December 2012)
- Guus Slotman (December 2012)
Part of this thesis has been publised in GJ Slotman, A. Fasolino, Structure, stability and defects of single layer h-BN in comparison to graphene. J Phys Condens Matter. 2013 Jan 30;25(4):045009
- Inka Locht (August 2012)
- Jonas Sweep (April 2012)
- Joost van den Ende (October 2011)
Friction at the atomic scale, sliding velocity studied with the Frenkel-Kontorova model (pdf (pdf, 9,7 MB)Friction is an important force, which affects both daily life and industrial processes. The origins of friction at the micro or nanoscale are not yet understood. At the nanoscale, the structure of the interface of the contacting surfaces (commensurate or incommensurate) is important. The Frenkel-Kontorova model, which describes the motion of a chain onto a periodic substrate, is suitable to understand the role of commensurability on friction. The periodic force due to the substrate can excite resonantly the phonons of the sliding chain and thereby dissipate mechanical energy into heat. We have studied this dissipation as a function of sliding velocity for different temperatures. We have found that resonant excitations remain important up to high temperatures for incommensurate interfaces and that the relation between friction and temperature is opposite for (in)commensurate interfaces. These effects could be studied with real systems, like solid xenon sliding on a silver substrate.
Part of this master thesis has been published in J. A. van den Ende, A. S. de Wijn, A. Fasolino, Effect of temperature and velocity on superlubricit. J. Phys.: Condens. Matter 24 (2012) 445009. A synopsis is given in Superlubricity in slow motion.
- Merel van Wijk (August 2011)
- Koen Reijnders (August 2011)
Semiclassical theory of chiral tunneling in grapheneThis thesis presents a systematic study of tunneling effects in graphene in graphene n-p and n-p-n junctions. Using a variety of semiclassical methods, we derive very precise analytical formulae for the transmission and reflection coefficients for the case of angular scattering off a one-dimensional potential barrier. These formulae turn out to be in good agreement with numerical calculations.
This master thesis has been awarded the``Universitaire studieprijs 2012'' of the Radboud university Nijmegen (see movie).
- Leendertjan Karssemeijer (August 2010)
Part of this thesis has been published in L.J. Karssemeijer, A Fasolino, Phonons of graphene and graphitic materials derived from the empirical potential LCBOPII. Surface Science 605 (17), 1611-1615
- Jaap Kroes (May 2010)
Part of this thesis has been published in J.M. H. Kroes, M.A. Akhukov, J.H. Los, N. Pineau, and A. Fasolino, Mechanism and free-energy barrier of the type-57 reconstruction of the zigzag edge of graphene. Phys. Rev. B 83, 165411 (2011)
- Titus van Erp (1999)
A. Fasolino, O. Radulescu, and T. Janssen Pinning and phonon localization in Frenkel-Kontorova models on quasiperiodic substrates. Phys. Rev. B 60, 6522