Source author record

Karsten W. Jacobsen

Karsten W. Jacobsen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Computational exfoliation of atomically thin 1D materials with application to Majorana bound states

We introduce a computational database with calculated structural, thermodynamic, electronic, magnetic, and optical properties of 820 one-dimensional materials. The materials are systematically selected and exfoliated from experimental databases of crystal structures based on a dimensionality scoring parameter. The database is furthermore expanded by chemical element substitution in the materials. The materials are investigated in both their bulk form and as isolated one-dimensional components. We discuss the methodology behind the database, give an overview of some of the calculated properties, and look at patterns and correlations in the data. The database is furthermore applied in computational screening to identify materials, which could exhibit Majorana bound states.

preprint2022arXiv

Numerical Quality Control for DFT-based Materials Databases

Electronic-structure theory is a strong pillar of materials science. Many different computer codes that employ different approaches are used by the community to solve various scientific problems. Still, the precision of different packages has only recently been scrutinized thoroughly, focusing on a specific task, namely selecting a popular density functional, and using unusually high, extremely precise numerical settings for investigating 71 monoatomic crystals. Little is known, however, about method- and code-specific uncertainties that arise under numerical settings that are commonly used in practice. We shed light on this issue by investigating the deviations in total and relative energies as a function of computational parameters. Using typical settings for basis sets and k-grids, we compare results for 71 elemental and 63 binary solids obtained by three different electronic-structure codes that employ fundamentally different strategies. On the basis of the observed trends, we propose a simple, analytical model for the estimation of the errors associated with the basis-set incompleteness. We cross-validate this model using ternary systems obtained from the NOMAD Repository and discuss how our approach enables the comparison of the heterogeneous data present in computational materials databases.

preprint2021arXiv

Room temperature structure and energetics of water-hydroxyl layers on Pt(111)

The interactions between water and hydroxyl species on Pt(111) surfaces have been intensely investigated due to their importance to fuel cell electrocatalysis. Here we present a room temperature molecular dynamics study of their structure and energetics using an ensemble of neural network potentials, which allow us to obtain unprecedented statistical sampling. We first study the energetics of hydroxyl formation, where we find a near-linear adsorption energy profile, which exhibits a soft and gradual increase in the differential adsorption energy at high hydroxyl coverages. This is strikingly different from the predictions of the conventional bilayer model, which displays a kink at 1/3ML OH coverage indicating a sizeable jump in differential adsorption energy, but within the statistical uncertainty of previously reported ab initio molecular dynamics studies. We then analyze the structure of the interface, where we provide evidence for the water-OH/Pt(111) interface being hydrophobic at high hydroxyl coverages. We furthermore explain the observed adsorption energetics by analyzing the hydrogen bonding in the water-hydroxyl adlayers, where we argue that the increase in differential adsorption energy at high OH coverage can be explained by a reduction in the number of hydrogen bonds from the adsorbed water molecules to the hydroxyls.

preprint2019arXiv

Minimum-Strain Symmetrization of Bravais Lattices

Bravais lattices are the most fundamental building blocks of crystallography. They are classified into groups according to their translational, rotational, and inversion symmetries. In computational analysis of Bravais lattices, fulfilment of symmetry conditions is usually determined by analysis of the metric tensor, using either a numerical tolerance to produce a binary (i.e. yes or no) classification, or a distance function which quantifies the deviation from an ideal lattice type. The metric tensor, though, is not scale-invariant, which complicates the choice of threshold and the interpretation of the distance function. Here, we quantify the distance of a lattice from a target Bravais class using strain. For an arbitrary lattice, we find the minimum-strain transformation needed to fulfil the symmetry conditions of a desired Bravais lattice type; the norm of the strain tensor is used to quantify the degree of symmetry breaking. The resulting distance is invariant to scale and rotation, and is a physically intuitive quantity. By symmetrizing to all Bravais classes, each lattice can be placed in a 14 dimensional space, which we use to create a map of the space of Bravais lattices and the transformation paths between them. A software implementation is available online under a permissive license.

preprint2016arXiv

Defect Tolerant Monolayer Transition Metal Dichalcogenides

Localized electronic states formed inside the band gap of a semiconductor due to crystal defects can be detrimental to the material's optoelectronic properties. Semiconductors with lower tendency to form defect induced deep gap states are termed defect tolerant. Here we provide a systematic first principles investigation of defect tolerance in 29 monolayer transition metal dichalcogenides (TMDs) of interest for nanoscale optoelectronics. We find that the TMDs based on group VI and X metals form deep gap states upon creation of a chalcogen (S, Se, Te) vacancy while the TMDs based on group IV metals form only shallow defect levels and are thus predicted to be defect tolerant. Interestingly, all the defect sensitive TMDs have valence and conduction bands with very similar orbital composition. This indicates a bonding/anti-bonding nature of the gap which in turn suggests that dangling bonds will fall inside the gap. These ideas are made quantitative by introducing a descriptor that measures the degree of similarity of the conduction and valence band manifolds. Finally, the study is generalized to non-polar nanoribbons of the TMDs where we find that only the defect sensitive materials form edge states within the band gap.

preprint2015arXiv

Anharmonic stabilization and band gap renormalization in the perovskite CsSnI$_3$

Amongst the X(Sn,Pb)Y$_3$ perovskites currently under scrutiny for their photovoltaic applications, the cubic B-$α$ phase of CsSnI$_3$ is arguably the best characterized experimentally. Yet, according to the standard harmonic theory of phonons, this deceptively simple phase should not exist at all due to rotational instabilities of the SnI$_6$ octahedra. Here, employing self-consistent phonon theory we show that these soft modes are stabilized at experimental conditions through anharmonic phonon-phonon interactions between the Cs ions and their iodine cages. We further calculate the renormalization of the electronic energies due to vibrations and find an unusual opening of the band gap, estimated as 0.24 and 0.11 eV at 500 and 300 K, which we attribute to the stretching of Sn-I bonds. Our work demonstrates the important role of temperature in accurately describing these materials.

preprint2015arXiv

Plasmons on the edge of MoS2 nanostructures

Using ab initio calculations we predict the existence of one-dimensional (1D), atomically confined plasmons at the edges of a zigzag MoS2 nanoribbon. The strongest plasmon originates from a metallic edge state localized on the sulfur dimers decorating the Mo edge of the ribbon. A detailed analysis of the dielectric function reveals that the observed deviations from the ideal 1D plasmon behavior result from single-particle transitions between the metallic edge state and the valence and conduction bands of the MoS2 sheet. The Mo and S edges of the ribbon are clearly distinguishable in calculated spatially resolved electron energy loss spectrum owing to the different plasmonic properties of the two edges. The edge plasmons could potentially be utilized for tuning the photocatalytic activity of MoS2 nanoparticles.

preprint2013arXiv

Spatially resolved quantum plasmon modes in metallic nano-films from first principles

Electron energy loss spectroscopy (EELS) can be used to probe plasmon excitations in nanostructured materials with atomic-scale spatial resolution. For structures smaller than a few nanometers quantum effects are expected to be important, limiting the validity of widely used semi-classical response models. Here we present a method to identify and compute spatially resolved plasmon modes from first principles based on a spectral analysis of the dynamical dielectric function. As an example we calculate the plasmon modes of 0.5-4 nm thick Na films and find that they can be classified as (conventional) surface modes, sub-surface modes, and a discrete set of bulk modes resembling standing waves across the film. We find clear effects of both quantum confinement and non-local response. The quantum plasmon modes provide an intuitive picture of collective excitations of confined electron systems and offer a clear interpretation of spatially resolved EELS spectra.

preprint2012arXiv

Ab-initio non-equilibrium quantum transport and forces with the real space projector augmented wave method

We present an efficient implemention of a non-equilibrium Green function (NEGF) method for self-consistent calculations of electron transport and forces in nanostructured materials. The electronic structure is described at the level of density functional theory (DFT) using the projector augmented wave method (PAW) to describe the ionic cores and an atomic orbital basis set for the valence electrons. External bias and gate voltages are treated in a self-consistent manner and the Poisson equation with appropriate boundary conditions is solved in real space. Contour integration of the Green function and parallelization over k-points and real space makes the code highly efficient and applicable to systems containing several hundreds of atoms. The method is applied to a number of different systems demonstrating the effects of bias and gate voltages, multiterminal setups, non-equilibrium forces, and spin transport.

preprint2012arXiv

Conventional and Acoustic Surface Plasmons on Noble Metal Surfaces: A Time-dependent Density Functional Theory Study

First-principles calculations of the conventional and acoustic surface plasmons (CSPs and ASPs) on the (111) surfaces of Cu, Ag, and Au are presented. The effect of $s-d$ interband transitions on both types of plasmons is investigated by comparing results from the local density approximation and an orbital dependent exchange-correlation (xc) potential that improves the position and width of the $d$ bands. The plasmon dispersions calculated with the latter xc-potential agree well with electron energy loss spectroscopy (EELS) experiments. For both the CSP and ASP, the same trend of Cu$<$Au$<$Ag is found for the plasmon energies and is attributed to the reduced screening by interband transitions from Cu, to Au and Ag. This trend for the ASP, however, contradicts a previous model prediction. While the ASP is seen as a weak feature in the EELS, it can be clearly identified in the static and dynamic dielectric band structure.

preprint2012arXiv

First-principles study of the phonon-limited mobility in n-type single-layer MoS2

In the present work we calculate the phonon-limited mobility in intrinsic n-type single-layer MoS2 as a function of carrier density and temperature for T > 100 K. Using a first-principles approach for the calculation of the electron-phonon interaction, the deformation potentials and Fröhlich interaction in the isolated MoS2 layer are determined. We find that the calculated room-temperature mobility of ~410 cm^2 V^-1 s^-1 is dominated by optical phonon scattering via deformation potential couplings and the Fröhlich interaction with the deformation potentials to the intravalley homopolar and intervalley longitudinal optical phonons given by 4.1 x 10^8 eV/cm and 2.6 x 10^8 eV/cm, respectively. The mobility is weakly dependent on the carrier density and follows a μ~ T^-1 temperature dependence with γ= 1.69 at room temperature. It is shown that a quenching of the characteristic homopolar mode which is likely to occur in top-gated samples, boosts the mobility with 70 cm^2 V^-1 s^-1 and can be observed as a decrease in the exponent to γ= 1.52. Our findings indicate that the intrinsic phonon-limited mobility is approached in samples where a high-kappa dielectric that effectively screens charge impurities is used as gate oxide.

preprint2012arXiv

Optical properties of bulk semiconductors and graphene/boron-nitride: The Bethe-Salpeter equation with derivative discontinuity-corrected DFT energies

We present an efficient implementation of the Bethe-Salpeter equation (BSE) for optical properties of materials in the projector augmented wave method GPAW. Single-particle energies and wave functions are obtained from the GLLBSC functional which explicitly includes the derivative discontinuity, is computationally inexpensive, and yields excellent fundamental gaps. Electron-hole interactions are included through the BSE using the statically screened interaction evaluated in the random phase approximation. For a representative set of semiconductors and insulators we find excellent agreement with experiments for the dielectric functions, onset of absorption, and lowest excitonic features. For the two-dimensional systems of graphene and hexagonal boron-nitride (h-BN) we find good agreement with previous many-body calculations. For the graphene/h-BN interface, we find that the fundamental and optical gaps of the h-BN layer are reduced by 2.0 eV and 0.7 eV, respectively, compared to freestanding h-BN. This reduction is due to image charge screening which shows up in the GLLBSC calculation as a reduction (vanishing) of the derivative discontinuity.

preprint2012arXiv

Unraveling the acoustic electron-phonon interaction in graphene

Using a first-principles approach we calculate the acoustic electron-phonon couplings in graphene for the transverse (TA) and longitudinal (LA) acoustic phonons. Analytic forms of the coupling matrix elements valid in the long-wavelength limit are found to give an almost quantitative description of the first-principles based matrix elements even at shorter wavelengths. Using the analytic forms of the coupling matrix elements, we study the acoustic phonon-limited carrier mobility for temperatures 0-200 K and high carrier densities of 10^{12}-10^{13} cm^{-2}. We find that the intrinsic effective acoustic deformation potential of graphene is Ξ_eff = 6.8 eV and that the temperature dependence of the mobility μ~ T^{-α} increases beyond an α= 4 dependence even in the absence of screening when the full coupling matrix elements are considered. The large disagreement between our calculated deformation potential and those extracted from experimental measurements (18-29 eV) indicates that additional or modified acoustic phonon-scattering mechanisms are at play in experimental situations.

preprint2011arXiv

Ab initio van der Waals interactions in simulations of water alter structure from mainly tetrahedral to high-density-like

The structure of liquid water at ambient conditions is studied in ab initio molecular dynamics simulations using van der Waals (vdW) density-functional theory, i.e. using the new exchange-correlation functionals optPBE-vdW and vdW-DF2. Inclusion of the more isotropic vdW interactions counteracts highly directional hydrogen-bonds, which are enhanced by standard functionals. This brings about a softening of the microscopic structure of water, as seen from the broadening of angular distribution functions and, in particular, from the much lower and broader first peak in the oxygen-oxygen pair-correlation function (PCF), indicating loss of structure in the outer solvation shells. In combination with softer non-local correlation terms, as in the new parameterization of vdW-DF, inclusion of vdW interactions is shown to shift the balance of resulting structures from open tetrahedral to more close-packed. The resulting O-O PCF shows some resemblance with experiment for high-density water (A. K. Soper and M. A. Ricci, Phys. Rev. Lett., 84:2881, 2000), but not directly with experiment for ambient water. However, an O-O PCF consisting of a linear combination of 70% from vdW-DF2 and 30% from experiment on low-density liquid water reproduces near-quantitatively the experimental O-O PCF for ambient water, indicating consistency with a two-liquid model with fluctuations between high- and low-density regions.

preprint2011arXiv

Computational design of chemical nanosensors: Transition metal doped single-walled carbon nanotubes

We present a general approach to the computational design of nanostructured chemical sensors. The scheme is based on identification and calculation of microscopic descriptors (design parameters) which are used as input to a thermodynamic model to obtain the relevant macroscopic properties. In particular, we consider the functionalization of a (6,6) metallic armchair single-walled carbon nanotube (SWNT) by nine different 3d transition metal (TM) atoms occupying three types of vacancies. For six gas molecules (N_{2}, O_{2}, H_{2}O, CO, NH_{3}, H_{2}S) we calculate the binding energy and change in conductance due to adsorption on each of the 27 TM sites. For a given type of TM functionalization, this allows us to obtain the equilibrium coverage and change in conductance as a function of the partial pressure of the "target" molecule in a background of atmospheric air. Specifically, we show how Ni and Cu doped metallic (6,6) SWNTs may work as effective multifunctional sensors for both CO and NH_{3}. In this way, the scheme presented allows one to obtain macroscopic device characteristics and performance data for nanoscale (in this case SWNT) based devices.

preprint2011arXiv

Linear density response function in the projector-augmented wave method: Applications to solids, surfaces, and interfaces

We present an implementation of the linear density response function within the projector-augmented wave (PAW) method with applications to the linear optical and dielectric properties of both solids, surfaces, and interfaces. The response function is represented in plane waves while the single-particle eigenstates can be expanded on a real space grid or in atomic orbital basis for increased efficiency. The exchange-correlation kernel is treated at the level of the adiabatic local density approximation (ALDA) and crystal local field effects are included. The calculated static and dynamical dielectric functions of Si, C, SiC, AlP and GaAs compare well with previous calculations. While optical properties of semiconductors, in particular excitonic effects, are generally not well described by ALDA, we obtain excellent agreement with experiments for the surface loss function of the Mg(0001) surface with plasmon energies deviating by less than 0.2 eV. Finally, we apply the method to study the influence of substrates on the plasmon excitations in graphene. On SiC(0001), the long wavelength $π$ plasmons are significantly damped although their energies remain almost unaltered. On Al(111) the $π$ plasmon is completely quenched due to the coupling to the metal surface plasmon.

preprint1992arXiv

Melting a copper cluster: Critical droplet theory

We simulate the melting of a 71 A diameter cluster of Cu. At low temperatures the crystal exhibits facets. With increasing temperatures the open facets pre-melt, the melted regions coalesce into a liquid envelope containing a crystalline nucleus, and the nucleus finally goes unstable to the supercooled liquid. Using critical droplet theory and experimental data for Cu, we explain the thermodynamics of the coexistence region. The width of the transition scales as (Number of particles) to the power (-1/4).