Source author record

Claudia Draxl

Claudia Draxl 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

44works
14topics
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

44 published item(s)

preprint2026arXiv

From Knowledge to Action: Outcomes of the 2025 Large Language Model (LLM) Hackathon for Applications in Materials Science and Chemistry

Large language models (LLMs) are rapidly changing how researchers in materials science and chemistry discover, organize, and act on scientific knowledge. This paper analyzes a broad set of community-developed LLM applications in an effort to identify emerging patterns in how these systems can be used across the scientific research lifecycle. We organize the projects into two complementary categories: Knowledge Infrastructure, systems that structure, retrieve, synthesize, and validate scientific information; and Action Systems, systems that execute, coordinate, or automate scientific work across computational and experimental environments. The submissions reveal a shift from single-purpose LLM tools toward integrated, multi-agent workflows that combine retrieval, reasoning, tool use, and domain-specific validation. Prominent themes include retrieval-augmented generation as grounding infrastructure, persistent structured knowledge representations, multimodal and multilingual scientific inputs, and early progress toward laboratory-integrated closed-loop systems. Together, these results suggest that LLMs are evolving from general-purpose assistants into composable infrastructure for scientific reasoning and action. This work provides a community snapshot of that transition and a practical taxonomy for understanding emerging LLM-enabled workflows in materials science and chemistry.

preprint2022arXiv

A consistent picture of excitations in cubic BaSnO$_{3}$ revealed by combining theory and experiment

Among the transparent conducting oxides, the perovskite barium stannate is most promising for various electronic applications due to its outstanding carrier mobility achieved at room temperature. However, most of its important characteristics, such as band gaps, effective masses, and absorption edge, remain controversial. Here, we provide a fully consistent picture by combining state-of-the-art {\it ab initio} methodology with forefront electron energy-loss spectroscopy and optical absorption measurements. Valence electron energy-loss spectra, featuring signals originating from band gap transitions, are acquired on defect-free sample regions of a BaSnO$_{3}$ single crystal. These high-energy-resolution measurements are able to capture also very weak excitations below the optical gap, attributed to indirect transitions. By temperature-dependent optical absorption measurements, we assess band-gap renormalization effects induced by electron-phonon coupling. Overall, we find for the effective electronic mass, the direct and the indirect gap, the optical gap, as well as the absorption onsets and spectra, excellent agreement between both experimental techniques and the theoretical many-body results, supporting also the picture of a phonon-mediated mechanism where indirect transitions are activated by phonon-induced symmetry lowering. This work demonstrates a fruitful connection between different high-level theoretical and experimental methods for exploring the characteristics of advanced materials.

preprint2022arXiv

All-electron many-body approach to resonant inelastic x-ray scattering

We present a formalism for the resonant inelastic x-ray scattering (RIXS) cross section. The resulting compact expression in terms of polarizability matrix elements, particularly lends itself to the implementation in an all-electron many-body perturbation theory (MBPT) framework, which is realized in the full-potential package exciting. With the carbon K edge RIXS of diamond and the oxygen K edge RIXS of $β-\mathrm{Ga_2 O_3}$, respectively, we demonstrate the importance of electron-hole correlation and atomic coherence in the RIXS spectra.

preprint2022arXiv

Density-of-states similarity descriptor for unsupervised learning from materials data

We develop a materials descriptor based on the electronic density of states and investigate the similarity of materials based on it. As an application example, we study the Computational 2D Materials Database that hosts thousands of two-dimensional materials with their properties calculated by density-functional theory. Combining our descriptor with a clustering algorithm, we identify groups of materials with similar electronic structure. We characterize these clusters in terms of their crystal structure, their atomic composition, and the respective electronic configurations to rationalize the found (dis)similarities.

preprint2022arXiv

Ehrenfest dynamics implemented in the all-electron package exciting

Ehrenfest Dynamics combined with real-time time-dependent density functional theory has proven to be a reliable tool to study non-adiabatic molecular dynamics with a reasonable computational cost. Among other possibilities, it allows for assessing in real time electronic excitations generated by ultra-fast laser pulses, as e.g., in pump-probe spectroscopy, and their coupling to the nuclear vibrations even beyond the linear regime. In this work, we present its implementation in the all-electron full-potential package exciting. Three cases are presented as examples: diamond and cubic BN relaxed after an initial lattice distortion, and cubic BN exposed to a laser pulse. Comparison with the Octopus code exhibits good agreement.

preprint2022arXiv

FAIR data enabling new horizons for materials research

The prosperity and lifestyle of our society are very much governed by achievements in condensed matter physics, chemistry and materials science, because new products for sectors such as energy, the environment, health, mobility and information technology (IT) rely largely on improved or even new materials. Examples include solid-state lighting, touchscreens, batteries, implants, drug delivery and many more. The enormous amount of research data produced every day in these fields represents a gold mine of the twenty-first century. This gold mine is, however, of little value if these data are not comprehensively characterized and made available. How can we refine this feedstock; that is, turn data into knowledge and value? For this, a FAIR (findable, accessible, interoperable and reusable) data infrastructure is a must. Only then can data be readily shared and explored using data analytics and artificial intelligence (AI) methods. Making data 'findable and AI ready' (a forward-looking interpretation of the acronym) will change the way in which science is carried out today. In this Perspective, we discuss how we can prepare to make this happen for the field of materials science.

preprint2022arXiv

How ferroelectric BaTiO$_3$ can tune a two-dimensional electron gas at the interface of LaInO3 and BaSnO$_3$: a first-principles study

The emerging interest in two-dimensional electron gases (2DEGs), formed at interfaces between two insulating oxide perovskites poses crucial fundamental question in view of future electronic devices. In the framework of density-functional theory, we investigate the possibility to control the characteristics of the 2DEG formed at the LaInO$_3$/BaSnO$_3$ interface by including a ferroelectric BaTiO$_3$ layer. To do so, we examine how the orientation of the ferroelectric polarization impacts density and confinement of the 2DEG. We find that aligning the ferroelectric polarization toward (outward) the LaInO$_3$/BaSnO$_3$ interface leads to an accumulation (depletion) of the interfacial 2DEG. Varying its magnitude, we find a linear effect on the 2DEG charge density that is confined within the BaSnO$_3$ side. Analysis of the optimized geometries revels that inclusion of the BaTiO$_3$ block makes structural distortions at the LaInO$_3$/BaSnO$_3$ less pronounced, which, in turn, enhances the 2DEG density. Thicker ferroelectric layers allow for reaching higher polarization magnitudes. We discuss the mechanisms behind all these findings and rationalize how the characteristics of both 2DEG and 2D hole gases can be controlled in the considered heterostructures. Overall, our results can be generalized to other combinations of ferroelectric, polar, and nonpolar materials.

preprint2022arXiv

Influence of spin-orbit coupling on chemical bonding

The influence of spin-orbit interaction on chemical bonds in elemental solids and homonuclear dimers is analyzed by means of density-functional-theory calculations. Employing highly precise all-electron full-potential methodology, our results represent benchmark quality. Comparison of the scalar- and fully-relativistic approaches for elemental solids shows that the spin-orbit interaction may contract or expand the volume of the considered material. The largest variation of the volume is obtained for Au, Tl, I, Bi, Po and Hg, exhibiting changes between 1.0--7.6\%. Using the tight-binding model, we show for diatomic molecules that the nature of this effect lies in the angular rearrangement of bonding and antibonding orbitals introduced by spin-orbit coupling. Such an angular rearrangement appears in partially filled $p$- or $d$-orbitals in heavy elements. Finally, we discuss the impact of the relativistic effects on the chemical bonding in single-layer iodides and transition metal dichalcogenides

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.

preprint2022arXiv

Robust excitons across the phase transition of two-dimensional hybrid perovskites

Two-dimensional halide perovskites are among intensely studied materials platforms profiting from solution based growth and chemical flexibility. They feature exceptionally strong interactions among electronic, optical as well as vibrational excitations and hold a great potential for future optoelectronic applications. A key feature for these materials is the occurrence of structural phase transitions that can impact their functional properties, including the electronic band gap and optical response dominated by excitons. However, to what extent the phase-transitions in two-dimensional perovskites alter the fundamental exciton properties remains barely explored so far. Here, we study the influence of the phase transition on both exciton binding energy and exciton diffusion, demonstrating their robust nature across the phase transition. These findings are unexpected in view of the associated substantial changes of the free carrier masses, strongly contrast broadly considered effective mass and drift-diffusion transport mechanisms, highlighting the unusual nature of excitons in two-dimensional perovskites.

preprint2022arXiv

Tuning two-dimensional electron (hole) gases at LaInO$_{3}$/BaSnO$_{3}$ interfaces: Impact of polar distortions, termination, and thickness

Two-dimensional election gases (2DEG), arising due to quantum confinement at interfaces between transparent conducting oxides, have received tremendous attention in view of electronic applications. The challenge is to find a material system that exhibits both a high charge-carrier density and mobility, at and even above room temperature. Here, we explore the potential of interfaces formed by two lattice-matched wide-gap oxides of emerging interest, $\textit{i.e.}$, the polar, orthorhombic perovskite LaInO$_{3}$ and the non-polar, cubic perovskite BaSnO$_{3}$, employing density-functional theory and many-body theory. We demonstrate that this material combination exhibits all the key features for reaching the goal. For periodic heterostructures, we find that the polar discontinuity at the interface is mainly compensated by electronic relaxation through charge transfer from the LaInO$_{3}$ to the BaSnO$_{3}$ side. This leads to the formation of a 2DEG hosted by the highly-dispersive Sn-$s$-derived conduction band and a 2D hole gas of O-$p$ character, strongly localized inside LaInO$_{3}$. Remarkably, structural distortions through octahedra tilts induce a depolarization field counteracting the polar discontinuity, and thus increasing the $critical$ (minimal) LaInO$_{3}$ thickness, $t_c$, required for the formation of a 2DEG. These polar distortions decrease with increasing LaInO$_{3}$ thickness, enhancing the polar discontinuity and leading to a 2DEG density of 0.5 electron per unit-cell surface. Interestingly, in non-periodic heterostructures, these distortions lead to a decrease of $t_c$, thereby enhancing and delocalizing the 2DEG. We rationalize how polar distortions, termination, and thickness can be exploited in view of tailoring the 2DEG characteristics, and why this material is superior to the most studied prototype LaAlO$_{3}$/SrTiO$_{3}$.

preprint2022arXiv

Two-dimensional electron gas at LaInO$_3$/BaSnO$_3$ interfaces controlled by a ferroelectric layer

With the example of LaInO$_{3}$/BaSnO$_3$, we demonstrate how both density and distribution of a two-dimensional electron gas (2DEG) formed at the interface between these perovskite oxides, can be efficiently controlled by a ferroelectric functional material. A polarization induced in a BaTiO$_3$ layer pointing toward the interface enhances the polar discontinuity which, in turn, significantly increases the 2DEG density and confinement, while, the opposite polarization depletes the 2DEG. Our predictions and analysis, based on first-principles calculations, can serve as a guide for designing such material combinations to be used in electronic devices.

preprint2021arXiv

Elastic stability of Ga$_2$O$_3$: Addressing the $β$ to $α$ phase transition from first principles

Elastic and structural properties of $β$-Ga$_2$O$_3$ and $α$-Ga$_2$O$_3$ are investigated from first principles. The full elastic tensors and elastic moduli of both phases at $0$ K are computed in the framework of semi-local density-functional theory. We determine mechanical instabilities of $β$-Ga$_2$O$_3$ by evaluating the full stiffness tensor under load for a range of hydrostatic pressure values. While a phase transition from the $β$ to $α$ phase is found to be energetically favored at $2.6$ GPa, we show that the $β$ phase is only mechanically unstable for much higher pressures ($>30$ GPa), which agrees well with experimental results. Our employed approach is based on the Born stability criterion, is independent of crystal symmetry, and thus can be readily applied to different materials.

preprint2021arXiv

Fingerprints of optical absorption in the perovskite LaInO$_{3}$: Insight from many-body theory and experiment

We provide a combined theoretical and experimental study of the electronic structure and the optical absorption edge of the orthorhombic perovskite LaInO$_{3}$. Employing density-functional theory and many-body perturbation theory, we predict a direct electronic quasiparticle band gap of about 5 eV and an effective electron (hole) mass of 0.31 (0.48) m$_{0}$. We find the lowest-energy excitation at 0.2 eV below the fundamental gap, reflecting a sizeable electron-hole attraction. Since the transition from the valence band maximum (VBM, $Γ$ point) is, however, dipole forbidden the onset is characterized by weak excitations from transitions around it. The first intense excitation appears about 0.32 eV above. Interestingly, this value coincides with an experimental value obtained by ellipsometry (4.80 eV) which is higher than the onset from optical absorption spectroscopy (4.35 eV). The latter discrepancy is attributed to the fact that the weak transitions that define the optical gap are not resolved by the ellipsometry measurement. The absorption edge shows a strong dependency on the light polarization, reflecting the character of the involved valence states. Temperature-dependent measurements show a redshift of the optical gap by about 120 meV by increasing the temperature from 5 to 300 K. Renormalization due to zero-point vibrations is extrapolated from the latter measurement to amount to 150 meV. By adding the excitonic binding energy of 0.2 eV obtained theoretically to the experimental optical absorption onset, we determine the fundamental band gap at room temperature to be 4.55 eV.

preprint2021arXiv

Rashba and Dresselhaus effects in 2D Pb-I-based perovskites

Bulk hybride halide perovskites are governed by significant Rashba and Dresselhaus splitting. This indicates that such effects will not only affect their optoelectronic properties but also those of their two dimensional layered relatives. This work aims at understanding how different ways of symmetry breaking influence these effects in those materials. For this purpose, model structures are adopted where the organic compounds are replaced by Cs atoms. Disregarding possible distortions in the inorganic layers, results in structures with composition Cs$_{n+1}$Pb$_n$I$_{3n+1}$. Using the all-electron full-potential density-functional-theory code \texttt{exciting}, the impact of atomic displacement on the band structure is systematically studied for $n=1$, 2, 3 and $\infty$. The displacement patterns that yield Rashba or Dresselhaus splitting are identified, and the amount of the splitting is determined as a function of displacement. Furthermore, the spin textures in the electronic states around the band gap are analyzed to differentiate between Rashba and Dresselhaus effects. This study reveals in-plane Pb displacements as the origin of the strongest effects.

preprint2020arXiv

Assessment of approaches for dispersive forces employing graphone as a case study

We have studied two interchange layer systems, (i) free standing partly hydrogenated graphene (graphone), and (ii) graphone on the Nickel (111) surface, to assess various density functional theory based computational schemes incorporating van der Waals forces. The various van der Waals methods have been employed ranging from the semiempirical force-field-like correction of Grimme, through non-local van der Waals density functionals, up to the functionals involving exact exchange and the random phase approximation for the correlation. Generally, all computational schemes lead to a similar qualitative picture of hydrogen layer physisorption and chemisorption to graphene. The largest discrepancies between the approaches emerge for the energetics of the investigated systems. Our studies shed light on the physical mechanisms of graphene hydrogenation both in vacuum and in the proximity of metallic surface. In particular, it is revealed that the adsorption of hydrogen atoms affects the nature of the bonding between graphene and the Ni(111) surface, from the weak to strong semi-covalent bonding. On the other hand, it turns out that the adsorption of hydrogen layer to graphene is stronger in the presence of the metallic surface.

preprint2020arXiv

Band gap renormalization, carrier mobilities, and the electron-phonon self-energy in crystalline naphthalene

Organic molecular crystals are expected to feature appreciable electron-phonon interactions that influence their electronic properties at zero and finite temperature. In this work, we report first-principles calculations and an analysis of the electron-phonon self-energy in naphthalene crystals. We compute the zero-point renormalization and temperature dependence of the fundamental band gap, and the resulting scattering lifetimes of electronic states near the valence- and conduction-band edges employing density functional theory. Further, our calculated phonon renormalization of the $GW$-corrected quasiparticle band structure predicts a fundamental band gap of 5 eV for naphthalene at room temperature, in good agreement with experiments. From our calculated phonon-induced electron lifetimes, we obtain the temperature-dependent mobilities of electrons and holes in good agreement with experimental measurements at room temperatures. Finally, we show that an approximate energy self-consistent computational scheme for the electron-phonon self-energy leads to the prediction of strong satellite bands in the electronic band structure. We find that a single calculation of the self-energy can reproduce the self-consistent results of the band gap renormalization and electrical mobilities for naphthalene, provided that the on-the-mass-shell approximation is used, i.e., if the self-energy is evaluated at the bare eigenvalues.

preprint2020arXiv

High Performance Solution of Skew-symmetric Eigenvalue Problems with Applications in Solving the Bethe-Salpeter Eigenvalue Problem

We present a high-performance solver for dense skew-symmetric matrix eigenvalue problems. Our work is motivated by applications in computational quantum physics, where one solution approach to solve the so-called Bethe-Salpeter equation involves the solution of a large, dense, skew-symmetric eigenvalue problem. The computed eigenpairs can be used to compute the optical absorption spectrum of molecules and crystalline systems. One state-of-the art high-performance solver package for symmetric matrices is the ELPA (Eigenvalue SoLvers for Petascale Applications) library. We extend the methods available in ELPA to skew-symmetric matrices. This way, the presented solution method can benefit from the optimizations available in ELPA that make it a well-established, efficient and scalable library, such as GPU support. We compare performance and scalability of our method to the only available high-performance approach for skew-symmetric matrices, an indirect route involving complex arithmetic. In total, we achieve a performance that is up to 3.67 higher than the reference method using Intel's ScaLAPACK implementation. The runtime to solve the Bethe-Salpeter-Eigenvalue problem can be improved by a factor of 10. Our method is freely available in the current release of the ELPA library.

preprint2020arXiv

MoTe2 as a natural hyperbolic material across the visible and the ultraviolet region

Hyperbolic materials are of particular interest for the next generation of photonic and optoelectronic devices. Since artificial metamaterials are intrinsically limited by the size of their nanostructured components, there has been a hunt for natural hyperbolic materials in the last few years. In a first-principles work based on density-functional theory and many-body perturbation theory, we investigate the fundamental dielectric response of MoTe2 in monolayer, bilayer, and bulk form, and find that it is a natural type-II hyperbolic material with low losses between 3 and 6 eV. Going from the monolayer to the bulk, the energy window of hyperbolic dispersion is blue-shifted by a few tenths of an eV. We show that excitonic effects and optical anisotropy play a major role in the hyperbolic behavior of MoTe2. Our results confirm the potential of layered materials as hyperbolic media for opto-electronics, photonics, and nano-imaging applications.

preprint2020arXiv

Robust mixing in self-consistent linearized augmented planewave calculations

We devise a mixing algorithm for full-potential (FP) all-electron calculations in the linearized augmented planewave (LAPW) method. Pulay's direct inversion in the iterative subspace is complemented with the Kerker preconditioner and further improvements to achieve smooth convergence, avoiding charge sloshing and noise in the exchange-correlation potential. As the Kerker preconditioner was originally designed for the planewave basis, we have adapted it to the FP-LAPW method and implemented in the exciting code. Applications to the $2\times 2$ Au(111) surface with a vacancy and to the Pd(111) surface demonstrate that this approach and our implementation work reliably with both density and potential mixing.

preprint2020arXiv

Structural, electronic, and optical properties of periodic graphene/h-BN van der Waals heterostructures

The emerging interest in van der Waals heterostructures as new materials for opto-electronics and photonics poses questions about their stability and structure-property relations. In the framework of density-functional and many-body perturbation theory, we investigate the structural, electronic, and optical properties of periodic heterostructures formed by graphene and hexagonal boron nitride (h-BN). To understand how the constituents affect each other depending on the layer stacking, we examine 12 commensurate arrangements. We find that interaction with graphene improves the stability of bulk h-BN also in those configurations that are predicted to be energetically metastable. In return, the interaction with h-BN can open a band gap of a few hundred meV in graphene. Its actual size can be tuned by the arrangement of the layers. In the semiconducting configurations, the character and spatial distribution of optical excitations are affected by the specific stacking, that determines the electronic states involved in the transitions. Remarkably, six out of the 12 explored heterostructures remain semi-metallic.

preprint2020arXiv

The role of the interface in controlling the epitaxial relationship between orthorhombic $\text{LaInO}_\text{3}$ and cubic $\text{BaSnO}_\text{3}$

Epitaxial perovskite oxide interfaces with different symmetry of the epitaxial layers have attracted considerable attention due to the emergence of novel behaviors and phenomena. In this paper, we show by aberration corrected transmission electron microscopy that orthorhombic $\text{LaInO}_\text{3}$ films grow in form of three different types of domains on the cubic $\text{BaSnO}_\text{3}$ pseudosubstrate. Quantitative evaluation of our TEM data shows that $c_{pc}$-oriented and $a_{pc}/b_{pc}$-oriented domains are present with similar probability. While continuum elasticity theory suggests that $c_{pc}$-oriented domains should exhibit a significantly higher strain energy density than $a_{pc}/b_{pc}$-oriented domains, density functional calculations confirm that $c_{pc}$- and $a_{pc}$-oriented domains on $\text{BaSnO}_\text{3}$ have similar energies.

preprint2020arXiv

Work-function modification of PEG(thiol) adsorbed on the Au(111) surface: A first-principles study

The possibility of modifying the work function of electrodes is important for optimizing the energy barriers for charge-injection (extraction) at the interface to an organic material. In this study, we perform density-functional-theory calculations to investigate the impact of dithiol-terminated polyethylene glycol (PEG(thiol)) based self-assembled monolayers (SAMs) with different numbers of PEG repeat units on the work function of the Au(111) surface. We find that a monolayer of PEG(thiol) decreases the work function of the Au(111) surface, where the magnitude of this reduction strongly depends on the length of the PEG backbone. The main contribution arises from the dipole due to the adsorption-induced charge rearrangement at the interface. Our work reveals a pronounced odd-even effect, which can be traced back to the dipole moment of the PEG(thiol) layer.

preprint2019arXiv

Maximally localized Wannier functions within the (L)APW+LO method

We present a robust algorithm that computes (maximally localized) Wannier functions (WFs) without the need of providing an initial guess. Instead, a suitable starting point is constructed automatically from so-called local orbitals which are fundamental building blocks of the basis set within (linearized) augmented planewave methods. Our approach is applied to a vast variety of materials such as metals, bulk and low-dimensional semiconductors, and complex inorganic-organic hybrid interfaces. For the interpolation of electronic single-particle energies, an accuracy in the meV range can be easily achieved. We exemplify the capabilities of our method by the calculation of the joint density of states in aluminum, (generalized) Kohn-Sham and quasi-particle band structures in various semiconductors, and the electronic structure of $β$-Ga$_2$O$_3$, including electron and hole effective masses.

preprint2019arXiv

Photoemission Signatures of Non-Equilibrium Carrier Dynamics from First Principles

Time- and angle-resolved photoemission spectroscopy (tr-ARPES) constitutes a powerful tool to inspect the dynamics and thermalization of hot carriers. The identification of the processes that drive the dynamics, however, is challenging even for the simplest systems owing to the coexistence of several relaxation mechanisms. Here, we devise a Green's function formalism for predicting the tr-ARPES spectral function and establish the origin of carrier thermalization entirely from first principles. The predictive power of this approach is demonstrated by an excellent agreement with experiments for graphene over time scales ranging from a few tens of femtoseconds up to several picoseconds. Our work provides compelling evidence of a non-equilibrium dynamics dominated by the establishment of a hot-phonon regime.

preprint2016arXiv

Ab-initio approach for gap plasmonics

Gap plasmonics deals with the properties of surface plasmons in the narrow region between two metallic nanoparticles forming the gap. For sub-nanometer gap distances electrons can tunnel between the nanoparticles leading to the emergence of novel charge transfer plasmons. These are conveniently described within the quantum corrected model by introducing an artificial material with a tunnel conductivity inside the gap region. Here we develop a methodology for computing such tunnel conductivities within the first-principles framework of density functional theory, and apply our approach to a jellium model representative for sodium. We show that the frequency dependence of the tunnel conductivity at infrared and optical frequencies can be significantly more complicated than previously thought.

preprint2016arXiv

Accurate all-electron $G_0W_0$ quasiparticle energies employing the full-potential augmented planewave method

The $GW$ approach of many-body perturbation theory (MBPT) has become a common tool for calculating the electronic structure of materials. However, with increasing number of published results, discrepancies between the values obtained by different methods and codes become more and more apparent. For a test set of small- and wide-gap semiconductors, we demonstrate how to reach the numerically \emph{best} electronic structure within the framework of the full-potential linearized augmented planewave (FLAPW) method. We first evaluate the impact of local orbitals in the Kohn-Sham eigenvalue spectrum of the underlying starting point. The role of the basis-set quality is then further analyzed when calculating the $G_0W_0$ quasiparticle energies. Our results, computed with the \exciting{} code, are compared to those obtained using the projector-augmented planewave (PAW) formalism, finding overall, good agreement between both methods. We also provide data produced with a typical FLAPW basis set as a benchmark for other $G_0W_0$ implementations.

preprint2016arXiv

Atomic signatures of local environment from core-level spectroscopy in $β$-Ga$_2$O$_3$

We present a joint theoretical and experimental study on core-level excitations from the oxygen $K$ edge of $β$-Ga$_2$O$_3$. A detailed analysis of the electronic structure reveals the importance of O-Ga hybridization effects in the conduction region. The spectrum from O 1$s$ core electrons is dominated by excitonic effects, which overall redshift the absorption onset by 0.5 eV, and significantly redistribute the intensity to lower energies. Analysis of the spectra obtained within many-body perturbation theory reveals atomic fingerprints of the inequivalent O atoms. From the comparison of energy-loss near-edge fine-structure (ELNES) spectra computed with respect to different crystal planes, with measurements recorded under the corresponding diffraction conditions, we show how the spectral contributions of specific O atoms can be enhanced while quenching others. These results suggest ELNES, combined with ab initio many-body theory, as a very powerful technique to characterize complex systems, with sensitivity to individual atomic species and to their local environment.

preprint2016arXiv

Mapping atomic orbitals with the transmission electron microscope: Images of defective graphene predicted from first-principles theory

Transmission electron microscopy has been a promising candidate for mapping atomic orbitals for a long time. Here, we explore its capabilities by a first principles approach. For the example of defected graphene, exhibiting either an isolated vacancy or a substitutional nitrogen atom, we show that three different kinds of images are to be expected, depending on the orbital character. To judge the feasibility of visualizing orbitals in a real microscope, the effect of the optics aberrations is simulated. We demonstrate that, by making use of energy-filtering, it should indeed be possible to map atomic orbitals in a state-of-the-art transmission electron microscope.

preprint2016arXiv

Optical properties of azobenzene-functionalized self-assembled monolayers: Intermolecular coupling and many-body interactions

In a joint theoretical and experimental work the optical properties of azobenzene-functionalized self-assembled monolayers (SAMs) are studied at different molecular packing densities. Our results, based on density-functional and many-body perturbation theory, as well as on differential reflectance (DR) spectroscopy, shed light on the microscopic mechanisms ruling photo-absorption in these systems. While the optical excitations are intrinsically excitonic in nature, regardless of the molecular concentration, in densely-packed SAMs intermolecular coupling and local-field effects are responsible for a sizable weakening of the exciton binding strength. Through a detailed analysis of the character of the electron-hole pairs, we show that distinct excitations involved in the photo-isomerization at low molecular concentrations are dramatically broadened by intermolecular interactions. Spectral shifts in the calculated DR spectra are in good agreement with the experimental results. Our findings represent an important step forward to rationalize the excited-state properties of these complex materials.

preprint2016arXiv

Polymorphism in $α$-sexithiophene crystals: Relative stability and transition path

We present a joint theoretical and experimental study to investigate polymorphism in $α$-sexithiophene (6T) crystals. By means of density-functional theory calculations, we clarify that the low-temperature phase is favorable over the high-temperature one, with higher relative stability by about 50 meV/molecule. This result is in agreement with our thermal desorption measurements. We also propose a transition path between the high- and low-temperature 6T polymorphs, estimating an upper bound for the energy barrier of about 1 eV/molecule. The analysis of the electronic properties of the investigated 6T crystal structures complements our study.

preprint2016arXiv

Probing the LDA-1/2 method as a starting point for $G_0W_0$ calculations

Employing the $G_0W_0$ approximation of Hedin's $GW$ approach one can obtain quasi-particle energies of extended systems and molecules with good accuracy. However, for many materials, semi-local exchange-correlation functionals are unsatisfactory starting points for $G_0W_0$ calculations. Hybrid functionals often improve upon them, but at a substantially higher computational cost. As an alternative, we suggest the LDA-1/2 method, which provides reasonable band gaps, without being computationally involved. In this work, we systematically compare 3 starting points for $G_0W_0$: LDA, PBE0, and LDA-1/2. A selection of solids is chosen for this benchmark: C, Si, SiC, AlP, LiF, MgO, Ne, Ar, GaN, GaAs, CdS, ZnS, and ZnO. We demonstrate that LDA-1/2 is a good starting point in most cases, reducing the mean absolute error of band gaps by 50% when compared to the other 2 functionals.

preprint2015arXiv

Big Data of Materials Science - Critical Role of the Descriptor

Statistical learning of materials properties or functions so far starts with a largely silent, non-challenged step: the choice of the set of descriptive parameters (termed descriptor). However, when the scientific connection between the descriptor and the actuating mechanisms is unclear, causality of the learned descriptor-property relation is uncertain. Thus, trustful prediction of new promising materials, identification of anomalies, and scientific advancement are doubtful. We analyse this issue and define requirements for a suited descriptor. For a classical example, the energy difference of zincblende/wurtzite and rocksalt semiconductors, we demonstrate how a meaningful descriptor can be found systematically.

preprint2015arXiv

Bound excitons and many-body effects in x-ray absorption spectra of azobenzene-functionalized self-assembled monolayers

We study x-ray absorption spectra of azobenzene-functionalized self-assembled monolayers (SAMs), investigating excitations from the nitrogen K edge. Azobenzene with H-termination and functionalized with CF3 groups is considered. The Bethe-Salpeter equation is employed to compute the spectra, including excitonic effects, and to determine the character of the near-edge resonances. Our results indicate that core-edge excitations are intense and strongly bound: Their binding energies range from about 6 to 4 eV, going from isolated molecules to densely-packed SAMs. Electron-hole correlation rules these excitations, while the exchange interaction plays a negligible role.

preprint2015arXiv

Computing Equilibrium Shapes of Wurtzite Crystals: The Example of GaN

Crystal morphologies are important for the design and functionality of devices based on low-dimensional nanomaterials. The equilibrium crystal shape (ECS) is a key quantity in this context. It is determined by surface energies, which are hard to access experimentally but can generally be well predicted by first-principles methods. Unfortunately, this is not necessarily so for polar and semipolar surfaces of wurtzite crystals. By extending the concept of Wulff construction, we demonstrate that the ECSs can nevertheless be obtained for this class of materials. For the example of GaN, we identify different crystal shapes depending on the chemical potential, shedding light on experimentally observed GaN nanostructures.

preprint2015arXiv

Estimating Excitonic Effects in the Absorption Spectra of Solids: Problems and Insight from a Guided Iteration Scheme

A major obstacle for computing optical spectra of solids is the lack of reliable approximations for capturing excitonic effects within time-dependent density-functional theory. We show that the trustful prediction of strongly bound electron-hole pairs within this framework using simple approximations is still a challenge and that available promising results have to be revisited. Deriving a set of analytical formula we analyze and explain the difficulties. We deduce an alternative approximation from an iterative scheme guided by previously available knowledge, significantly improving the description of exciton binding energies. Finally, we show how one can "read" exciton binding energies from spectra determined in the random phase approximation, without any further calculation.

preprint2015arXiv

Fingerprint of fractional charge transfer at metal/organic interface

Although physisorption is a widely occurring mechanism of bonding at the organic/metal interface, contradictory interpretations of this phenomenon are often reported. Photoemission and X-ray absorption spectroscopy investigations of nanorods of a substituted pentacene, 2,3,9,10-tetrafluoropentacene, deposited on gold single crystals reveal to be fundamental to identify the bonding mechanisms. We find fingerprints of a fractional charge transfer from the clean metal substrate to the physisorbed molecules. This phenomenon is unambiguously recognizable by a non-rigid shift of the core-level main lines while the occupied states at the interface stay mostly unperturbed, and the unoccupied states experience pronounced changes. The experimental results are corroborated by first-principles calculations.

preprint2015arXiv

LayerOptics: Microscopic modeling of optical coefficients in layered materials

Theoretical spectroscopy is a powerful tool to describe and predict optical properties of materials. While nowadays routinely performed, first-principles calculations only provide bulk dielectric tensors in Cartesian coordinates. These outputs are hardly comparable with experimental data, which are typically given by macroscopic quantities, crucially depending on the laboratory setup. Even more serious discrepancies can arise for anisotropic materials, e.g., organic crystals, where off-diagonal elements of the dielectric tensor can significantly contribute to the spectral features. Here, we present LayerOptics, a versatile and user-friendly implementation, based on the solution of the Maxwell's equations for anisotropic materials, to compute optical coefficients in anisotropic layered materials. We apply this tool for post-processing full dielectric tensors of molecular materials, including excitonic effects, as computed from many-body perturbation theory using the exciting code. For prototypical examples, ranging from optical to X-ray frequencies, we show the importance of combining accurate ab initio methods to obtain dielectric tensors, with the solution of the Maxwell's equations to compute optical coefficients accounting for optical anisotropy of layered systems. Good agreement with experimental data supports the potential of our approach, in view of achieving microscopic understanding of spectroscopic properties in complex materials.

preprint2015arXiv

Optical spectra from molecules to crystals: Insight from many-body perturbation theory

Time-dependent density-functional theory (TDDFT) often successfully reproduces excitation energies of finite systems, already in the adiabatic local-density approximation (ALDA). Here we show for prototypical molecular materials, i.e., oligothiophenes, that ALDA largely fails and explain why this is so. By comparing TDDFT with an in-depth analysis based on many-body perturbation theory, we demonstrate that correlation effects crucially impact energies and character of the optical excitations not only for molecules of increasing length and in crystalline environment, but even for isolated small molecules. We argue that only high-level methodologies, which explicitly include correlation effects, can reproduce optical spectra of molecular materials with equal accuracy from gas phase to crystal structures.

preprint2013arXiv

Dynamic structure factors of Cu, Ag, and Au: A comparative study from first principles

We present a comparative theoretical and experimental study of dynamic structure factors (momentum-dependent loss functions) of the noble metals Cu, Ag, and Au in the energy range of 0$-$60 eV. The emphasis is on theoretical results that are compared with new as well as available experimental data. Dynamic structure factors are calculated within the linear-response formalism of time-dependent density-functional theory, using the full-potential linearized augmented plane-wave (FP-LAPW) method. For the studied energy range, local-field effects are found to be very important for Ag and Au and only marginally relevant for Cu. We present an explanation for this surprising behavior. Loss functions of all three metals possess a complex multi-peak structure. We classify the features in the loss function as being related to collective excitations, interband transitions, or mixed modes. The impact of short-range correlations on the dynamic response functions are evaluated by comparing the results of the random-phase approximation to those of the time-dependent local-density approximation. Exchange correlation effects are found to be weak for small momentum transfers, but increasingly important for larger momenta. The calculated structure factors agree well with experiments.

preprint2013arXiv

Excited State Properties of Organic Semiconductors: Breakdown of the Tamm-Dancoff Approximation

The solution of the Bethe-Salpeter equation within the framework of many-body perturbation theory has turned out to be a benchmark for ab-initio calculations of optical properties of semiconductors and insulators. Commonly, however, the coupling between the resonant and anti-resonant excitations is neglected which is referred to as the Tamm-Dancoff approximation (TDA). This is well justified in cases where the exciton binding energy is much smaller than the band gap. Here, we report on the optical properties of a representative series of organic semiconductors where we find the TDA to no longer hold. We observe an increase of the exciton binding energy of up to 44% thereby improving the agreement with experiment.

preprint2013arXiv

Magnetic control of spin-orbit fields: a first-principles study of Fe/GaAs junctions

The microscopic structure of spin-orbit fields for the technologically important Fe/GaAs interface is uncovered from first principles. A symmetry based method allows to obtain the spin-orbit fields---both their magnitude and orientation---for a generic Bloch state, from the electronic band structure for any in-plane magnetization orientation. It is demonstrated that the spin-orbit fields depend not only on the electric field across the interface, but also surprisingly strongly on the Fe magnetization orientation, opening prospects for their magnetic control. These results give important clues in searching for spin-orbit transport and optical phenomena in ferromagnetic/nonmagnetic systems.

preprint2013arXiv

Predicting the electronic structure of weakly interacting hybrid systems: The example of nanosized pea-pod structures

We provide a simple scheme for predicting the electronic structure of van-der-Waals bound systems, based on the mere knowledge of the electronic structure of the subunits. We demonstrate this with the example of nano-peapods, consisting of polythiophene encapsulated in single-wall carbon nanotubes. Using density functional theory we disentangle the contributions to the level alignment. The main contribution is shown to be given by the ionization potential of the polymer inside the host, which, in turn is determined by the curvature of the tube. Only a small correction arises from charge redistributions within the domains of the constituents. Polarization effects turn out to be minor due to the cylindrical geometry of the peapods and their dielectric characteristics. Our findings open a perspective towards designing opto-electronic properties of such complex materials.

preprint2012arXiv

Band renormalization of a polymer physisorbed on graphene investigated by many-body perturbation theory

Many-body perturbation theory at the $G_0W_0$ level is employed to study the electronic properties of poly(\emph{para}-phenylene) (PPP) on graphene. Analysis of the charge density and the electrostatic potential shows that the polymer-surface interaction gives rise to the formation of only weak surface dipoles with no charge transfer between the polymer and the surface. In the local-density approximation (LDA) of density-functional theory, the band structure of the combined system appears as a superposition of the eigenstates of its constituents. Consequently, the LDA band gap of PPP remains unchanged upon adsorption onto graphene. $G_0W_0$ calculations, however, renormalize the electronic levels of the weakly physisorbed polymer. Thereby, its band gap is considerably reduced compared to that of the isolated PPP chain. This effect can be understood in terms of image charges induced in the graphene layer, which allows us to explain the quasi-particle gap of PPP versus polymer-graphene distance by applying a classical image-potential model. For distances below 4.5 Å, however, deviations from this simple classical model arise which we qualitatively explain by taking into account the polarizablity of the adsorbate. For a quantitative description with predictive power, however, we emphasize the need for an accurate ab-initio description of the electronic structure for weakly coupled systems at equilibrium bonding distances.