Source author record

Leeor Kronik

Leeor Kronik 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

26works
9topics
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

26 published item(s)

preprint2026arXiv

Microscopic theory of phonon polaritons and long wavelength dielectric response

We present a first-principles approach for calculating phonon-polariton dispersion relations. In this approach, phonon-photon interaction is described by quantization of a Hamiltonian that describes harmonic lattice vibrations coupled with the electromagnetic field inside the material. All Hamiltonian parameters are obtained from first-principles calculations, with diagonalization leading to non-interacting polariton quasiparticles. This method naturally includes retardation effects and resolves non-analytical behavior and ambiguities in phonon frequencies at the Brillouin zone center, especially in non-cubic and optically anisotropic materials. Furthermore, by incorporating higher-order terms in the Hamiltonian, we also account for quasiparticle interactions and spectral broadening. Specifically, we show how anharmonic effects in phonon polaritons lead to a dielectric response that challenges traditional models. The accuracy and consequences of the approach are demonstrated on GaP and GaN as harmonic test systems and PbTe and $β$-Ga$_2$O$_3$ as anharmonic test systems.

preprint2022arXiv

An Optimally-Tuned Starting Point for Single-Shot $GW$ Calculations of Solids

The dependence of ab initio many-body perturbation theory within the $GW$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $GW$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely PBE, PBE0, and HSE, as well as to self-consistent $GW$ schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $GW$ (QS$GW$) and eigenvalue self-consistent $GW$ (ev$GW$). We also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

preprint2019arXiv

Mechanical and Tribological Properties of Layered Materials Under High Pressure: Assessing the Importance of Many-Body Dispersion Effects

The importance of many-body dispersion effects in layered materials subjected to high external loads is evaluated. State-of-the-art many-body dispersion density functional theory calculations performed for graphite, hexagonal boron nitride, and their hetero-structures were used to fit the parameters of a classical registry-dependent interlayer potential. Using the latter, we performed extensive equilibrium molecular dynamics simulations and studied the mechanical response of homogeneous and heterogeneous bulk models under hydrostatic pressures up to 30 GPa. Comparison with experimental data demonstrates that the reliability of the many-body dispersion model extends deep into the sub-equilibrium regime. Friction simulations demonstrate the importance of many-body dispersion effects for the accurate description of the tribological properties of layered materials interfaces under high pressure.

preprint2016arXiv

Interplay of bias-driven charging and the vibrational Stark effect in molecular junctions

We observe large, reversible, bias driven changes in the vibrational energies of PCBM, based on simultaneous transport and surface-enhanced Raman spectroscopy (SERS) measurements on PCBM-gold junctions. A combination of linear and quadratic shifts in vibrational energies with voltage is analyzed and compared with similar measurements involving C60-gold junctions. A theoretical model based on density functional theory (DFT) calculations suggests that both a vibrational Stark effect and bias-induced charging of the junction contribute to the shifts in vibrational energies. In the PCBM case, a linear vibrational Stark effect is observed due to the permanent electric dipole moment of PCBM. The vibrational Stark shifts shown here for PCBM junctions are comparable to or larger than the charging effects that dominate in C60 junctions.

preprint2016arXiv

Optical Phonons in Methylammonium Lead Halide Perovskites and Implications for Charge Transport

Lead-halide perovskites are promising materials for opto-electronic applications. Recent reports indicated that their mechanical and electronic properties are strongly affected by the lattice vibrations. Herein we report far-infrared spectroscopy measurements of CH$_{3}$NH$_{3}$Pb(I/Br/Cl)$_{3}$ thin films and single crystals at room temperature and a detailed quantitative analysis of the spectra. We find strong broadening and anharmonicity of the lattice vibrations for all three halide perovskites, which indicates dynamic disorder of the lead-halide cage at room temperature. We determine the frequencies of the transversal and longitudinal optical phonons, and use them to calculate, via appropriate models, the static dielectric constants, polaron masses, electron-phonon coupling constants, and upper limits for the phonon-scattering limited charge carrier mobilities. Within the limitations of the model used, we can place an upper limit of 200$\,$cm$^{2}$V$^{-1}$s$^{-1}$ for the room temperature charge carrier mobility in MAPbI$_{3}$ single crystals. Our findings are important for the basic understanding of charge transport processes and mechanical properties in metal halide perovskites.

preprint2016arXiv

Structural and excited-state properties of oligoacene crystals from first principles

Molecular crystals are a prototypical class of van der Waals (vdW) bound organic materials with excited state properties relevant for optoelectronics applications. Predicting the structure and excited state properties of molecular crystals presents a challenge for electronic structure theory, as standard approximations to density functional theory (DFT) do not capture long range vdW dispersion interactions and do not yield excited state properties. In this work, we use a combination of DFT including vdW forces) using both non local correlation functionals and pair wise correction methods (together with many body perturbation theory (MBPT) to study the geometry and excited states, respectively, of the entire series of oligoacene crystals, from benzene to hexacene. We find that vdW methods can predict lattice constants within 1 percent of the experimental measurements, on par with the previously reported accuracy of pairwise approximations for the same systems. We further find that excitation energies are sensitive to geometry, but if optimized geometries are used MBPT can yield excited state properties within a few tenths of an eV from experiment. We elucidate trends in MBPT computed charged and neutral excitation energies across the acene series and discuss the role of common approximations used in MBPT.

preprint2015arXiv

Atomically wired molecular junctions: Connecting a single organic molecule by chains of metal atoms

Using a break junction technique, we find a clear signature for the formation of conducting hybrid junctions composed of a single organic molecule (benzene, naphthalene or anthracene) connected to chains of platinum atoms. The hybrid junctions exhibit metallic-like conductance (~0.1-1G0), which is rather insensitive to further elongation by additional atoms. At low bias voltage the hybrid junctions can be elongated significantly beyond the length of the bare atomic chains. Ab initio calculations reveal that benzene based hybrid junctions have a significant binding energy and high structural flexibility that may contribute to the survival of the hybrid junction during the elongation process. The fabrication of hybrid junctions opens the way for combining the different properties of atomic chains and organic molecules to realize a new class of atomic scale interfaces.

preprint2015arXiv

Effect of ensemble generalization on the highest-occupied Kohn-Sham eigenvalue

There are several approximations to the exchange-correlation functional in density-functional theory that accurately predict total energy-related properties of many-electron systems, such as binding energies, bond lengths, and crystal structures. Other approximations are designed to describe potential-related processes, such as charge transfer and photoemission. However, the development of a functional which can serve the two purposes simultaneously is a long-standing challenge. Trying to address it, we employ in the current work the ensemble generalization procedure proposed in Phys. Rev. Lett. 110, 126403 (2013). Focusing on the prediction of the ionization potential via the highest occupied Kohn-Sham eigenvalue, we examine a variety of exchange-correlation approximations: the local spin-density approximation, semi-local generalized gradient approximations, and global and local hybrid functionals. Results for a test set of 26 diatomic molecules and single atoms are presented. We find that the aforementioned ensemble generalization systematically improves the prediction of the ionization potential, for various systems and exchange-correlation functionals, without compromising the accuracy of total energy-related properties. We specifically examine hybrid functionals. These depend on a parameter controlling the ratio of semi-local to non-local functional components. The ionization potential obtained with ensemble-generalized functionals is found to depend only weakly on the parameter value, contrary to common experience with non-generalized hybrids, thus eliminating one aspect of the so-called `parameter dilemma' of hybrid functionals.

preprint2015arXiv

Elimination of the asymptotic fractional dissociation problem in Kohn-Sham density functional theory using the ensemble-generalization approach

Many approximations within density-functional theory spuriously predict that a many-electron system can dissociate into fractionally charged fragments. Here, we revisit the case of dissociated diatomic molecules, known to exhibit this problem when studied within standard approaches, including the local spin-density approximation (LSDA). By employing our recently proposed [E. Kraisler and L. Kronik, Phys. Rev. Lett. 110, 126403 (2013)] ensemble-generalization we find that asymptotic fractional dissociation is eliminated in all systems examined, even if the underlying exchange-correlation (xc) is still the LSDA. Furthermore, as a result of the ensemble generalization procedure, the Kohn-Sham potential develops a spatial step between the dissociated atoms, reflecting the emergence of the derivative discontinuity in the xc energy functional. This step, predicted in the past for the exact Kohn-Sham potential and observed in some of its more advanced approximate forms, is a desired feature that prevents any fractional charge transfer between the system's fragments. It is usually believed that simple xc approximations such as the LSDA cannot develop this step. Our findings show, however, that ensemble generalization to fractional electron densities automatically introduces the desired step even to the most simple approximate xc functionals and correctly predicts asymptotic integer dissociation.

preprint2015arXiv

Experimental and theoretical electronic structure of quinacridone

The energy positions of frontier orbitals in organic electronic materials are often studied experimentally by (inverse) photoemission spectroscopy and theoretically within density functional theory. However, standard exchange-correlation functionals often result in too small fundamental gaps, may lead to wrong orbital energy ordering, and do not capture polarization-induced gap renormalization. Here, we examine these issues and a strategy for overcoming them by studying the gas phase and bulk electronic structure of the organic molecule quinacridone (5Q), a promising material with many interesting properties for organic devices. Experimentally, we perform angle-resolved photoemission spectroscopy (ARUPS) on thin films of the crystalline $β$-phase of 5Q. Theoretically, we employ an optimally-tuned range-separated hybrid functional (OT-RSH) within density functional theory. For the gas phase molecule, our OT-RSH result for the ionization potential (IP) represents a substantial improvement over the semi-local PBE and the PBE0 hybrid functional results, producing an IP in quantitative agreement with experiment. For the bulk crystal, we take into account the correct screening in the bulk, using the recently developed OT-SRSH approach, while retaining the optimally-tuned parameters for the range-separation and the short-range Fock exchange. This leads to a band gap narrowing due to polarization effects and results in a valence band spectrum in excellent agreement with experimental ARUPS data, with respect to both peak positions and heights. Finally, full-frequency $G_0W_0$ results based on a hybrid functional starting point are shown to agree with the OT-SRSH approach, improving substantially on the PBE-starting point.

preprint2015arXiv

Fundamental gaps with approximate density functionals: the derivative discontinuity revealed from ensemble considerations

The fundamental gap is a central quantity in the electronic structure of matter. Unfortunately, the fundamental gap is not generally equal to the Kohn-Sham gap of density functional theory (DFT), even in principle. The two gaps differ precisely by the derivative discontinuity, namely, an abrupt change in slope of the exchange-correlation (xc) energy as a function of electron number, expected across an integer-electron point. Popular approximate functionals are thought to be devoid of a derivative discontinuity, strongly compromising their performance for prediction of spectroscopic properties. Here we show that, in fact, all exchange-correlation functionals possess a derivative discontinuity, which arises naturally from the application of ensemble considerations within DFT, without any empiricism. This derivative discontinuity can be expressed in closed form using only quantities obtained in the course of a standard DFT calculation of the neutral system. For small, finite systems, addition of this derivative discontinuity indeed results in a greatly improved prediction for the fundamental gap, even when based on the most simple approximate exchange-correlation density functional - the local density approximation (LDA). For solids, the same scheme is exact in principle, but when applied to LDA it results in a vanishing derivative discontinuity correction. This failure is shown to be directly related to the failure of LDA in predicting fundamental gaps from total energy differences in extended systems.

preprint2015arXiv

One-electron self-interaction and the asymptotics of the Kohn-Sham potential: an impaired relation

One-electron self-interaction and an incorrect asymptotic behavior of the Kohn-Sham exchange-correlation potential are among the most prominent limitations of many present-day density functionals. However, a one-electron self-interaction-free energy does not necessarily lead to the correct long-range potential. This is here shown explicitly for local hybrid functionals. Furthermore, carefully studying the ratio of the von Weizsäcker kinetic energy density to the (positive) Kohn-Sham kinetic energy density, $τ_\mathrm{W}/τ$, reveals that this ratio, which frequently serves as an iso-orbital indicator and is used to eliminate one-electron self-interaction effects in meta-generalized-gradient approximations and local hybrid functionals, can fail to approach its expected value in the vicinity of orbital nodal planes. This perspective article suggests that the nature and consequences of one-electron self-interaction and some of the strategies for its correction need to be reconsidered.

preprint2015arXiv

Probing the Orbital Origin of Conductance Oscillations in Atomic Chains

We investigate periodical oscillations in the conductance of suspended Au and Pt atomic chains during elongation under mechanical stress. Analysis of conductance and shot noise measurements reveals that the oscillations are mainly related to variations in a specific conduction channel as the chain undergoes transitions between zigzag and linear atomic configurations. The calculated local electronic structure shows that the oscillations originate from varying degrees of hybridization between the atomic orbitals along the chain as a function of the zigzag angle. These variations are highly dependent on the directionality and symmetry of the relevant orbitals, in agreement with the order-of-magnitude difference between the Pt and Au oscillation amplitudes observed in experiment. Our results demonstrate that the sensitivity of conductance to structural variations can be controlled by designing atomic-scale conductors in view of the directional interactions between atomic orbitals.

preprint2015arXiv

Solid-state optical absorption from optimally-tuned time-dependent range-separated hybrid density functional theory

We present a framework for obtaining reliable solid-state charge and optical excitations and spectra from optimally-tuned range-separated hybrid density functional theory. The approach, which is fully couched within the formal framework of generalized Kohn-Sham theory, allows for accurate prediction of exciton binding energies. We demonstrate our approach through first principles calculations of one- and two-particle excitations in pentacene, a molecular semiconducting crystal, where our work is in excellent agreement with experiments and prior computations. We further show that with one adjustable parameter, set to produce an accurate bandgap, this method accurately predicts band structures and optical spectra of silicon and lithium flouride, prototypical covalent and ionic solids. Our findings indicate that for a broad range of extended bulk systems, this method may provide a computationally inexpensive alternative to many-body perturbation theory, opening the door to studies of materials of increasing size and complexity.

preprint2014arXiv

A multi-scale approach to the electronic structure of doped semiconductor surfaces

The inclusion of the global effects of semiconductor doping poses a unique challenge for first-principles simulations, because the typically low concentration of dopants renders an explicit treatment intractable. In systems which do not display long-range band bending, a satisfactory remedy is offered by the use of "pseudoatoms", with a fractional nuclear charge matching the bulk doping concentration. However, this alone is not always sufficient for charged surfaces, where the width of the space-charge region (SCR) often exceeds realistic supercell dimensions. One generalization of the pseudoatom approach which overcomes this difficulty relies on the introduction of an artificially high doping level within a slab calculation, in conjunction with a multi-scale electrostatic energy correction. Here, we present an alternative technique that naturally extends the pseudoatom approach while bypassing the need for calculations with an unrealistically high doping level. It is based on the introduction of excess charge, mimicking free charge carriers from the SCR, along with a fixed sheet of counter-charge mimicking the SCR-related field. Self-consistency is obtained by imposing charge conservation and Fermi level equilibration between the bulk, treated semi-classically, and the electronic states of the slab/surface, which are treated quantum-mechanically. The method, which we call CREST - the Charge-Reservoir Electrostatic Sheet Technique - can be used with standard electronic structure codes. We validate the approach using a simple tight-binding model, which allows for comparison of its results with calculations encompassing the full SCR explicitly. We then employ it with density functional theory, where it is used to obtain insights into the electronic structures of the "clean-cleaved" Si(111) surface and its buckled (2x1) reconstruction, at various doping densities.

preprint2014arXiv

A real-space method for highly parallelizable electronic transport calculations

We present a real-space method for first-principles nano-scale electronic transport calculations. We use the non-equilibrium Green's function method with density functional theory and implement absorbing boundary conditions (ABCs, also known as complex absorbing potentials, or CAPs) to represent the effects of the semi-infinite leads. In real space, the Kohn-Sham Hamiltonian matrix is highly sparse. As a result, the transport problem parallelizes naturally and can scale favorably with system size, enabling the computation of conductance in relatively large molecular junction models. Our use of ABCs circumvents the demanding task of explicitly calculating the leads' self-energies from surface Green's functions, and is expected to be more accurate than the use of the jellium approximation. In addition, we take advantage of the sparsity in real space to solve efficiently for the Green's function over the entire energy range relevant to low-bias transport. We illustrate the advantages of our method with calculations on several challenging test systems and find good agreement with reference calculation results.

preprint2014arXiv

A self-interaction-free local hybrid functional: Accurate binding energies vis-à-vis accurate ionization potentials from Kohn-Sham eigenvalues

We present and test a new approximation for the exchange-correlation (xc) energy of Kohn-Sham density functional theory. It combines exact exchange with a compatible non-local correlation functional. The functional is by construction free of one-electron self-interaction, respects constraints derived from uniform coordinate scaling, and has the correct asymptotic behavior of the xc energy density. It contains one parameter that is not determined ab initio. We investigate whether it is possible to construct a functional that yields accurate binding energies and affords other advantages, specifically Kohn-Sham eigenvalues that reliably reflect ionization potentials. Tests for a set of atoms and small molecules show that within our local-hybrid form accurate binding energies can be achieved by proper optimization of the free parameter in our functional, along with an improvement in dissociation energy curves and in Kohn-Sham eigenvalues. However, the correspondence of the latter to experimental ionization potentials is not yet satisfactory, and if we choose to optimize their prediction, a rather different value of the functional's parameter is obtained. We put this finding in a larger context by discussing similar observations for other functionals and possible directions for further functional development that our findings suggest.

preprint2014arXiv

A State Representation Approach for Atomistic Time-Dependent Transport Calculations in Molecular Junctions

We propose a new method for simulating electron dynamics in open quantum systems out of equilibrium, using a finite atomistic model. The proposed method is motivated by the intuitive and practical nature of the driven Liouville von-Neumann equation approach of Sánchez et al. [J. Chem. Phys., 124, 214708 (2006)]. A key ingredient of our approach is a transformation of the Hamiltonian matrix from an atomistic to a state representation of the molecular junction. This allows us to uniquely define the bias voltage across the system while maintaining a proper thermal electronic distribution within the finite lead models. Furthermore, it allows us to investigate complex molecular junctions, including non-linear setups and multi-lead configurations. A heuristic derivation of our working equation leads to explicit expressions for the damping and driving terms, which serve as appropriate electron sources and sinks that effectively "open" the finite model system. Although the method does not forbid it, in practice we find neither violation of Pauli's exclusion principles nor deviation from density matrix positivity throughout our numerical simulations of various tight-binding model systems. We believe that the new approach offers a practical and physically sound route for performing atomistic time-dependent transport calculations in realistic molecular junctions.

preprint2014arXiv

Deviations from piecewise linearity in the solid-state limit with approximate density functionals

In exact density functional theory (DFT) the total ground-state energy is a series of linear segments between integer electron points, a condition known as "piecewise linearity". Deviation from this condition is indicative of poor predictive capabilities for electronic structure, in particular of ionization energies, fundamental gaps, and charge transfer. In this article, we take a new look at the deviation from linearity (i.e., curvature) in the solid-state limit by considering two different ways of approaching it: a large finite system of increasing size and a crystal represented by an increasingly large reference cell with periodic boundary conditions. We show that the curvature approaches vanishing values in both limits, even for functionals which yield poor predictions of electronic structure, and therefore can not be used as a diagnostic or constructive tool in solids. We find that the approach towards zero curvature is different in each of the two limits, owing to the presence of a compensating background charge in the periodic case. Based on these findings, we present a new criterion for functional construction and evaluation, derived from the size-dependence of the curvature, along with a practical method for evaluating this criterion. For large finite systems we further show that the curvature is dominated by the self-interaction of the highest occupied eigenstate. These findings are illustrated by computational studies of various solids, semiconductor nanocrystals, and long alkane chains.

preprint2013arXiv

Inter-layer Potential for Hexagonal Boron Nitride

A new interlayer force-field for layered hexagonal boron nitride (h-BN) based structures is presented. The force-field contains three terms representing the interlayer attraction due to dispersive interactions, repulsion due to anisotropic overlaps of electron clouds, and monopolar electrostatic interactions. With appropriate parameterization, the potential is able to simultaneously capture well the binding and lateral sliding energies of planar h-BN based dimer systems as well as the interlayer telescoping and rotation of double walled boron-nitride nanotubes of different crystallographic orientations. The new potential thus enables the accurate and efficient modeling and simulation of large-scale h-BN based layered structures.

preprint2012arXiv

Curvature and frontier orbital energies in density functional theory

Perdew et al. [Phys. Rev. Lett 49, 1691 (1982)] discovered and proved two different properties of exact Kohn-Sham density functional theory (DFT): (i) The exact total energy versus particle number is a series of linear segments between integer electron points; (ii) Across an integer number of electrons, the exchange-correlation potential may "jump" by a constant, known as the derivative discontinuity (DD). Here, we show analytically that in both the original and the generalized Kohn-Sham formulation of DFT, the two are in fact two sides of the same coin. Absence of a derivative discontinuity necessitates deviation from piecewise linearity, and the latter can be used to correct for the former, thereby restoring the physical meaning of the orbital energies. Using selected small molecules, we show that this results in a simple correction scheme for any underlying functional, including semi-local and hybrid functionals as well as Hartree-Fock theory, suggesting a practical correction for the infamous gap problem of DFT. Moreover, we show that optimally-tuned range-separated hybrid functionals can inherently minimize both DD and curvature, thus requiring no correction, and show that this can be used as a sound theoretical basis for novel tuning strategies.

preprint2012arXiv

Quasiparticle and Optical Spectroscopy of Organic Semiconductors Pentacene and PTCDA from First Principles

The broad use of organic semiconductors for optoelectronic applications relies on quantitative understanding and control of their spectroscopic properties. Of paramount importance are the transport gap - the difference between ionization potential and electron affinity - and the exciton binding energy - inferred from the difference between the transport and optical absorption gaps. Transport gaps are commonly established via photoemission and inverse photoemission spectroscopy (PES/IPES). However, PES and IPES are surface-sensitive, average over a dynamic lattice, and are subject to extrinsic effects, leading to significant uncertainty in gaps. Here, we use density functional theory and many-body perturbation theory to calculate the spectroscopic properties of two prototypical organic semiconductors, pentacene and 3,4,9,10-perylene tetracarboxylic dianhydride (PTCDA), quantitatively comparing with measured PES, IPES, and optical absorption spectra. For bulk pentacene and PTCDA, the computed transport gaps are 2.4 and 3.0 eV, and optical gaps are 1.7 and 2.1 eV, respectively. Computed bulk quasiparticle spectra are in excellent agreement with surface-sensitive photoemission measurements over several eV only if the measured gap is reduced by 0.6 eV for pentacene, and 0.6-0.9 eV for PTCDA. We attribute this redshift to several physical effects, including incomplete charge screening at the surface, static and dynamical disorder, and experimental resolution. Optical gaps are in excellent agreement with experiment, with solid-state exciton binding energies of ~0.5 eV for both systems; for pentacene, the exciton is delocalized over several molecules and exhibits significant charge transfer character. Our parameter-free calculations provide new interpretation of spectroscopic properties of organic semiconductors critical to optoelectronics.

preprint2012arXiv

Quasiparticle spectra from a non-empirical optimally-tuned range-separated hybrid density functional

We present a method for obtaining outer valence quasiparticle excitation energies from a DFT-based calculation, with accuracy that is comparable to that of many-body perturbation theory within the GW approximation. The approach uses a range-separated hybrid density functional, with asymptotically exact and short-range fractional Fock exchange. The functional contains two parameters - the range separation and the short-range Fock fraction. Both are determined non-empirically, per system, based on satisfaction of exact physical constraints for the ionization potential and many-electron self-interaction, respectively. The accuracy of the method is demonstrated on four important benchmark organic molecules: perylene, pentacene, 3,4,9,10-perylene-tetracarboxylic-dianydride (PTCDA) and 1,4,5,8-naphthalene-tetracarboxylic dianhydride (NTCDA). We envision that for finite systems the approach could provide an inexpensive alternative to GW, opening the door to the study of presently out of reach large-scale systems.

preprint2011arXiv

Charge transport across metal/molecular (alkyl) monolayer-Si junctions is dominated by the LUMO level

We compare the charge transport characteristics of heavy doped p- and n-Si-alkyl chain/Hg junctions. Photoelectron spectroscopy (UPS, IPES and XPS) results for the molecule-Si band alignment at equilibrium show the Fermi level to LUMO energy difference to be much smaller than the corresponding Fermi level to HOMO one. This result supports the conclusion we reach, based on negative differential resistance in an analogous semiconductor-inorganic insulator/metal junction, that for both p- and n-type junctions the energy difference between the Fermi level and LUMO, i.e., electron tunneling, controls charge transport. The Fermi level-LUMO energy difference, experimentally determined by IPES, agrees with the non-resonant tunneling barrier height deduced from the exponential length-attenuation of the current.

preprint2010arXiv

Fundamental gaps of finite systems from the eigenvalues of a generalized Kohn-Sham method

We present a broadly-applicable, physically-motivated first-principles approach to determining the fundamental gap of finite systems. The approach is based on using a range-separated hybrid functional within the generalized Kohn-Sham approach to density functional theory. Its key element is the choice of a range-separation parameter such that Koopmans' theorem for both the neutral and anionic is obeyed as closely as possible. We demonstrate the validity, accuracy, and advantages of this approach on first, second, and third row atoms, the oligoacene family of molecules, and a set of hydrogen-passivated silicon nanocrystals. This extends the quantitative usage of density functional theory to an area long believed to be outside its reach.

preprint2010arXiv

Stacking and Registry Effects in Layered Materials: The Case of Hexagonal Boron Nitride

The interlayer sliding energy landscape of hexagonal boron nitride (h-BN) is investigated via a van der Waals corrected density functional theory approach. It is found that the main role of the van der Waals forces is to "anchor" the layers at a fixed distance, whereas the electrostatic forces dictate the optimal stacking mode and the interlayer sliding energy. A nearly free-sliding path is identified, along which bandgap modulations of ~0.6 eV are obtained. We propose a simple geometrical model that quantifies the registry matching between the layers and captures the essence of the corrugated h-BN interlayer energy landscape. The simplicity of this phenomenological model opens the way to the modeling of complex layered structures, such as carbon and boron nitride nanotubes.