Source author record

Roi Baer

Roi Baer 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

33works
8topics
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

33 published item(s)

preprint2022arXiv

Forces from stochastic density functional theory under nonorthogonal atom-centered basis sets

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis-set representation (Fabian et al. WIREs Comput Mol Sci. 2019;e1412. https://doi.org/10.1002/wcms.1412) and apply it to Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of $0.065\,eV/Å$ ($\sim1.2\times10^{-3}E_{h}/a_{0}$) when 120 stochastic orbitals are used, independently of systems size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

preprint2022arXiv

High frequency limit of spectroscopy

We consider an arbitrary quantum mechanical system, initially in its ground-state, exposed to a time-dependent electromagnetic pulse with a carrier frequency $ω_0$ and a slowly varying envelope of finite duration. By working out a solution to the time-dependent Schrödinger equation in the high-$ω_0$ limit, we find that, to the leading order in $ω_0^{-1}$, a perfect self-cancellation of the system's linear response occurs as the pulse switches off. Surprisingly, the system's observables are, nonetheless, describable in terms of a combination of its linear density response function and nonlinear functions of the electric field. An analysis of jellium slab and jellium sphere models reveals a very high surface sensitivity of the considered setup, producing a richer excitation spectrum than accessible within the conventional linear response regime. On this basis, we propose a new spectroscopic technique, which we provisionally name the Nonlinear High-Frequency Pulsed Spectroscopy (NLHFPS). Combining the advantages of the extraordinary surface sensitivity, the absence of constraints by the traditional dipole selection rules, and the clarity of theoretical interpretation utilizing the linear response time-dependent density functional theory, NLHFPS has a potential to evolve into a powerful characterization method for nanoscience and nanotechnology.

preprint2022arXiv

Linear scalability of density functional theory calculations without imposing electron localization

Linear scaling density functional theory approaches to electronic structure are often based on the tendency of electrons to localize even in large atomic and molecular systems. However, in many cases of actual interest, for example in semiconductor nanocrystals, system sizes can reach very large extension before significant electron localization sets in and the scaling of the numerical methods may deviate strongly from linear. Here, we address this class of systems, by developing a massively parallel density functional theory (DFT) approach which doesn't rely on electron localizationa and is formally quadratic scaling, yet enables highly efficient linear wall-time complexity in the weak scalability regime. The approach extends from the stochastic DFT method described in Fabian et. al. WIRES: Comp. Mol. Science, e1412 2019 but is fully deterministic. It uses standard quantum chemical atom-centered Gaussian basis sets for representing the electronic wave functions combined with Cartesian real space grids for some of the operators and for enabling a fast solver for the Poisson equation. Our main conclusion is, that when a processor-abundant high performance computing (HPC) infrastructure is available, this type of approach has the potential to allow the study of large systems in regimes where quantum confinement or electron delocalization prevents linear-scaling.

preprint2022arXiv

Time dependent second order Green's function theory for neutral excitations

We develop a time dependent second order Green's function theory (GF2) for calculating neutral excited states in molecules. The equation of motion for the lesser Green's function (GF) is derived within the adiabatic approximation to the Kadanoff-Baym (KB) equation using the second order Born approximation for the self-energy. In the linear response regime, we recast the time dependent KB equation into a Bethe-Salpeter-like equation (GF2-BSE), with a kernel approximated by the second order Coulomb self-energy. We then apply our GF2-BSE to a set of molecules and atoms and find that GF2-BSE is superior to configuration interaction with singles (CIS) and/or time dependent Hartree-Fock (TDHF), particularly for charge transfer excitations, and is comparable to CIS with perturbative doubles (CIS(D)) in most cases.

preprint2021arXiv

Tuning the range separation parameter in periodic systems

Kohn-Sham DFT with optimally tuned range-separated hybrid (RSH) functionals provides accurate and nonempirical fundamental gaps for a wide variety of finite-size systems. The standard tuning procedure relies on calculation of total energies of charged systems and thus cannot be applied to periodic solids. Here, we develop a framework for tuning the range separation parameter that can be used for periodic and open boundary conditions. The basic idea is to choose the range parameter that results in a stationary point where the fundamental gap obtained by RSH matches the gap obtained from a G0W0 over RSH calculation. The proposed framework is therefore analogous to eigenvalue self-consistent GW (scGW). We assess the method for various solids and obtain very good agreement with scGW results.

preprint2020arXiv

Dopant levels in large nanocrystals using stochastic optimally tuned range-separated hybrid density functional theory

We apply a stochastic version of an optimally tuned range-separated hybrid functional to provide insight on the electronic properties of P- and B- doped Si nanocrystals of experimentally relevant sizes. We show that we can use the range-separation parameter for undoped systems to calculate accurate results for dopant activation energies. We apply this strategy for tuning functionals to study doped nanocrystals up to 2.5 nm in diameter at the hybrid functional level. In this confinement regime, the P- and B- dopants have large activation energies and have strongly localized states that lie deep within the energy gaps. Structural relaxation plays a greater role for B-substituted dopants and contributes to the increase in activation energy when the B dopant is near the nanocrystal surface.

preprint2020arXiv

Efficient Langevin dynamics for "noisy" forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Stochastic approaches provide a gentler system-size dependency at the cost of introducing "noisy" forces, which serve to limit the efficiency of the sampling. In the first-order Langevin dynamics (FOLD), efficient sampling is achievable by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator (Mazzola et al., Phys. Rev. Lett., 118, 015703 (2017)). However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure which severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

preprint2020arXiv

Range-Separated Stochastic Resolution of Identity: Formulation and Application to Second Order Green's Function Theory

We develop a range-separated stochastic resolution of identity approach for the $4$-index electron repulsion integrals, where the larger terms (above a predefined threshold) are treated using a deterministic resolution of identity and the remaining terms are treated using a stochastic resolution of identity. The approach is implemented within a second-order Greens function formalism with an improved $O(N^3)$ scaling with the size of the basis set, $N$. Moreover, the range-separated approach greatly reduces the statistical error compared to the full stochastic version ({\it J. Chem. Phys.} {\bf 151}, 044144 (2019)), resulting in computational speedups of ground and excited state energies of nearly two orders of magnitude, as demonstrated for hydrogen dimer chains.

preprint2020arXiv

Stochastically Realized Observables for Excitonic Molecular Aggregates

We show that a stochastic approach enables calculations of the optical properties of large 2-dimensional and nanotubular excitonic molecular aggregates. Previous studies of such systems relied on numerically diagonalizing the dense and disordered Frenkel Hamiltonian, which scales approximately as $\mathcal{O}(N^3)$ for $N$ dye molecules. Our approach scales much more efficiently as $\mathcal{O}(N\log(N))$, enabling quick study of systems with a million of coupled molecules on the micron size scale. We calculate several important experimental observable including the optical absorption spectrum and density of states, and develop a stochastic formalism for the participation ratio. Quantitative agreement with traditional matrix diagonalization methods is demonstrated for both small- and intermediate-size systems. The stochastic methodology enables the study of the effects of spatial-correlation in site energies on the optical signatures of large 2D aggregates. Our results demonstrate that stochastic methods present a path forward for screening structural parameters and validating experiments and theoretical predictions in large excitonic aggregates.

preprint2019arXiv

Stochastic embedding DFT: theory and application to p-nitroaniline

Over this past decade, we combined the idea of stochastic resolution of identity with a variety of electronic structure methods. In our stochastic Kohn-Sham DFT method, the density is an average over multiple stochastic samples, with stochastic errors that decrease as the inverse square root of the number of sampling orbitals. Here we develop a stochastic embedding density functional theory method (se-DFT) that selectively reduces the stochastic error (specifically on the forces) for a selected sub-system(s). The motivation, similar to that of other quantum embedding methods, is that for many systems of practical interest the properties are often determined by only a small sub-system. In stochastic embedding DFT two sets of orbitals are used: a deterministic one associated with the embedded subspace, and the rest which is described by a stochastic set. The method is exact in the limit of large number of stochastic samples. We apply se-DFT to study a p-nitroaniline molecule in water, where the statistical errors in the forces on the system (the p-nitroaniline molecule) are reduced by an order of magnitude compared with non-embedding stochastic DFT.

preprint2018arXiv

Stochastic density functional theory

Linear-scaling implementations of density functional theory (DFT) reach their intended efficiency regime only when applied to systems having a physical size larger than the range of their Kohn-Sham density matrix (DM). This causes a problem since many types of large systems of interest have a rather broad DM range and are therefore not amenable to analysis using DFT methods. For this reason, the recently proposed stochastic DFT (sDFT), avoiding exhaustive DM evaluations, is emerging as an attractive alternative linear-scaling approach. This review develops a general formulation of sDFT in terms of a (non)orthogonal basis representation and offers an analysis of the statistical errors (SEs) involved in the calculation. Using a new Gaussian-type basis-set implementation of sDFT, applied to water clusters and silicon nanocrystals, it demonstrates and explains how the standard deviation and the bias depend on the sampling rate and the system size in various types of calculations. We also develop basis-set embedded-fragments theory, demonstrating its utility for reducing the SEs for energy, density of states and nuclear force calculations. Finally, we discuss the algorithmic complexity of sDFT, showing it has CPU wall-time linear-scaling. The method parallelizes well over distributed processors with good scalability and therefore may find use in the upcoming exascale computing architectures.

preprint2016arXiv

Spontaneous charge carrier localization in extended one-dimensional systems

Charge carrier localization in extended atomic systems has been described previously as being driven by disorder, point defects or distortions of the ionic lattice. Here we show for the first time by means of first-principles computations that charge carriers can spontaneously localize due to a purely electronic effect in otherwise perfectly ordered structures. Optimally-tuned range-separated density functional theory and many-body perturbation calculations within the GW approximation reveal that in trans-polyacetylene and polythiophene the hole density localizes on a length scale of several nanometers. This is due to exchange-induced translational symmetry breaking of the charge density. Ionization potentials, optical absorption peaks, excitonic binding energies and the optimally-tuned range parameter itself all become independent of polymer length as it exceeds the critical localization scale. Moreover, lattice disorder and the formation of a polaron result from the charge localization in contrast to the traditional view that lattice distortions precede charge localization. Our results can explain experimental findings that polarons in conjugated polymers form instantaneously after exposure to ultrafast light pulses.

preprint2016arXiv

Stochastic self-consistent Green's function second-order perturbation theory (sGF2)

The second-order Green's function method (GF2) was shown recently to be an accurate self-consistent approach for electronic structure of correlated systems since the self-energy accounts for both the weak and some of the strong correlation. The numerical scaling of GF2 is quite steep however, $O({N^5})$ (where the pre-factor is often hundreds), effectively preventing its application to large systems. Here, we develop a stochastic approach to GF2 (sGF2) where the self-energy is evaluated by a random-vector decomposition of Green's functions so that the dominant part of the calculation scales quasi linearly with system size. A study of hydrogen chains shows that the resulting approach is numerically efficient and accurate, as the stochastic errors are very small, 0.05% of the correlation energy for large systems with only a moderate computational effort. The method also yields automatically efficient MP2 energies and is automatically temperature dependent.

preprint2015arXiv

Smoothing and extrapolating shifted-contour auxiliary-field Monte Carlo signals using discrete Laguerre functions

We develop a new smoothing or extrapolating method, based on discrete Laguerre functions, for systematically analyzing the stochastic signal of shifted-contour auxiliary-field Monte Carlo. We study the statistical errors and extrapolation errors using full configuration-interaction energies for the doubly stretched water molecule. The only free parameter is the order N of the fit. We show that low N emphasizes stability while higher N enable improved extrapolation, at the cost of increased statistical errors. Typically, one should use low order for signals based on a small number of iterations while higher order is efficacious for signals based on large number of iterations. We provide a heuristic algorithm for determining the order to be used and show its utility.

preprint2015arXiv

Stochastic Optimally-Tuned Ranged-Separated Hybrid Density Functional Theory

We develop a stochastic formulation of the optimally-tuned range-separated hybrid density functional theory which enables significant reduction of the computational effort and scaling of the non-local exchange operator at the price of introducing a controllable statistical error. Our method is based on stochastic representations of the Coulomb convolution integral and of the generalized Kohn-Sham density matrix. The computational cost of the approach is similar to that of usual Kohn-Sham density functional theory, yet it provides much more accurate description of the quasiparticle energies for the frontier orbitals. This is illustrated for a series of silicon nanocrystals up to sizes exceeding 3000 electrons. Comparison with the stochastic GW many-body perturbation technique indicates excellent agreement for the fundamental band gap energies, good agreement for the band-edge quasiparticle excitations, and very low statistical errors in the total energy for large systems. The present approach has a major advantage over one-shot GW by providing a self-consistent Hamiltonian which is central for additional post-processing, for example in the stochastic Bethe-Salpeter approach.

preprint2015arXiv

Time-dependent Stochastic Bethe-Salpeter Approach

A time-dependent formulation for electron-hole excitations in extended finite systems, based on the Bethe-Salpeter equation (BSE), is developed using a stochastic wave function approach. The time-dependent formulation builds on the connection between time-dependent Hartree-Fock (TDHF) theory and configuration-interaction with single substitution (CIS) method. This results in a time dependent Schrödinger-like equation for the quasiparticle orbital dynamics based on an effective Hamiltonian containing direct Hartree and screened exchange terms, where screening is described within the Random Phase Approximation (RPA). To solve for the optical absorption spectrum, we develop a stochastic formulation in which the quasiparticle orbitals are replaced by stochastic orbitals to evaluate the direct and exchange terms in the Hamiltonian as well as the RPA screening. This leads to an overall quadratic scaling, a significant improvement over the equivalent symplectic eigenvalue representation of the BSE. Application of the time-dependent stochastic BSE (TDsBSE) approach to silicon and CdSe nanocrystals up to size of ~3000 electrons is presented and discussed.

preprint2014arXiv

Breaking the theoretical scaling limit for predicting quasi-particle energies: The stochastic GW approach

We develop a formalism to calculate the quasi-particle energy within the GW many-body perturbation correction to the density functional theory (DFT). The occupied and virtual orbitals of the Kohn-Sham (KS) Hamiltonian are replaced by stochastic orbitals used to evaluate the Green function, the polarization potential, and thereby the GW self-energy. The stochastic GW (sGW) relies on novel theoretical concepts such as stochastic time-dependent Hartree propagation, stochastic matrix compression and spatial/temporal stochastic decoupling techniques. Beyond the theoretical interest, the formalism enables linear scaling GW calculations breaking the theoretical scaling limit for GW as well as circumventing the need for energy cutoff approximations. We illustrate the method for silicon nanocrystals of varying sizes with over 3000 electrons.

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.

preprint2014arXiv

Embedded fragment stochastic density functional theory

We develop a method in which the electronic densities of small fragments determined by Kohn-Sham density functional theory (DFT) are embedded using stochastic DFT to form the exact density of the full system. The new method preserves the scaling and the simplicity of the stochastic DFT but cures the slow convergence that occurs when weakly coupled subsystems are treated. It overcomes the spurious charge fluctuations that impair the applications of the original stochastic DFT approach. We demonstrate the new approach on a fullerene dimer and on clusters of water molecules and show that the density of states and the total energy can be accurately described with a relatively small number of stochastic orbitals.

preprint2014arXiv

Sublinear scaling for time-dependent stochastic density functional theory

A stochastic approach to time-dependent density functional theory (TDDFT) is developed for computing the absorption cross section and the random phase approximation (RPA) correlation energy. The core idea of the approach involves time-propagation of a small set of stochastic orbitals which are first projected on the occupied space and then propagated in time according to the time-dependent Kohn-Sham equations. The evolving electron density is exactly represented when the number of random orbitals is infinite, but even a small number (? 16) of such orbitals is enough to obtain meaningful results for absorption spectrum and the RPA correlation energy per electron. We implement the approach for silicon nanocrystals (NCs) using real-space grids and find that the overall scaling of the algorithm is sublinear with computational time and memory.

preprint2013arXiv

A guided stochastic energy-domain formulation of the second order Møller-Plesset perturbation theory

We develop an alternative formulation in the energy-domain to calculate the second order Møller-Plesset (MP2) perturbation energies. The approach is based on repeatedly choosing four random energies using a non-separable guiding function, filtering four random orbitals at these energies, and averaging the resulting Coulomb matrix elements to obtain a statistical estimate of the MP2 correlation energy. In contrast to our time-domain formulation, the present approach is useful for both quantum chemistry and real-space/plane wave basis sets. The scaling of the MP2 calculation is roughly linear with system size, providing a useful tool to study dispersion energies in large systems. This is demonstrated on a structure of 64 fullerenes within the SZ basis as well as on silicon nanocrystals using real-space grids.

preprint2013arXiv

Self-averaging stochastic Kohn-Sham density functional theory

We formulate the Kohn-Sham density functional theory (KS-DFT) as a statistical theory in which the electron density is deter-mined from an average of correlated stochastic densities in a trace formula. The key idea is that it is sufficient to converge the total energy per electron to within a predefined statistical error in order to obtain reliable estimates of the electronic band structure, the forces on nuclei, the density and its moments, etc. The fluctuations in the total energy per electron are guaranteed to decay to zero as the system size increases. This facilitates "self-averaging" which leads to the first ever report of sublinear scaling KS-DFT electronic structure. The approach sidesteps calculation of the density matrix and thus is insensitive to its evasive sparseness, as demonstrated here for silicon nanocrystals. The formalism is not only appealing in terms of its promise to far push the limits of application of KS-DFT, but also represents a cognitive change in the way we think of electronic structure calculations as this stochastic theory seamlessly converges to the thermodynamic limit.

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

Expeditious stochastic approach for MP2 energies in large electronic systems

A fast stochastic method for calculating the 2nd order Møller-Plesset (MP2) correction to the correlation energy of large systems of electrons is presented. The approach is based on reducing the exact summation over occupied and unoccupied states to a time-dependent trace formula amenable to stochastic sampling. We demonstrate the abilities of the method to treat systems of thousands electrons using hydrogen passivated silicon spherical nanocrystals represented on a real space grids, much beyond capabilities of present day MP2 implementations.

preprint2012arXiv

Expeditious Stochastic Calculation of Random-Phase Approximation Energies for Thousands of Electrons in 3 Dimensions

A fast method is developed for calculating the Random-Phase-Approximation (RPA) correlation energy for density functional theory. The correlation energy is given by a trace over a projected RPA response matrix and the trace is taken by a stochastic approach using random perturbation vectors. The method scales, at most, quadratically with the system size but in practice, due to self-averaging, requires less statistical sampling as the system grows and the performance is close to linear scaling. We demonstrate the method by calculating the RPA correlation energy for cadmium selenide and silicon nanocrystals with over 1500 electrons. In contrast to 2nd order Møller-Plesset correlation energies, we find that the RPA correlation energies per electron are largely independent on the nanocrystal size.

preprint2012arXiv

Multiexciton generation in IV-VI nanocrystals: The role of carrier effective mass, band mixing, and phonon emission

We study the role of the effective mass, band mixing and phonon emission on multiexciton generation in IV-VI nanocrystals. A 4-band k dot p effective mass model, which allows for an independent variation of these parameters, is adopted to describe the electronic structure of the nanocrystals. Multiexciton generation efficiencies are calculated using a Green's function formalism, providing results that are numerically similar to impact excitation. We find that multiexciton generation efficiencies are maximized when the effective mass of the electron and hole are small and similar. Contact with recent experimental results for multiexciton generation in PbS and PbSe is made.

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.

preprint2010arXiv

Can impact excitation explain efficient carrier multiplication in carbon nanotube photodiodes?

We address recent experiments (Science 325, 1367 (2009)) reporting on highly efficient multiplication of electron-hole pairs in carbon nanotube photodiodes at photon energies near the carrier multiplication threshold (twice the quasi-particle band gap). This result is surprising in light of recent experimental and theoretical work on multiexciton generation in other confined materials, such as semiconducting nanocrystals. We propose a detailed mechanism based on carrier dynamics and impact excitation resulting in highly efficient multiplication of electron-hole pairs. We discuss the important time and energy scales of the problem and provide analysis of the role of temperature and the length of the diode.

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

Ground-state degeneracies leave recognizable topological scars in the one-particle density

In Kohn-Sham density functional theory (KS-DFT) a fictitious system of non-interacting particles is constructed having the same ground-state (GS) density as the physical system of interest. A fundamental open question in DFT concerns the ability of an exact KS calculation to spot and characterize the GS degeneracies in the physical system. In this article we provide theoretical evidence suggesting that the GS density, as a function of position on a 2D manifold of parameters affecting the external potential, is "topologically scarred" in a distinct way by degeneracies. These scars are sufficiently detailed to enable determination of the positions of degeneracies and even the associated Berry phases. We conclude that an exact KS calculation can spot and characterize the degeneracies of the physical system.

preprint2008arXiv

A density functional theory for symmetric radical cations from bonding to dissociation

It is known for quite some time that approximate density functional (ADF) theories fail disastrously when describing the dis-sociative symmetric radical cations R2+. Considering this dissociation limit, previous work has shown that Hartree-Fock (HF) theory favors the R+1--R0 charge distribution while DF approximations favor the R+0.5 -- R+0.5. Yet, general quantum mechanical principles indicate that both these (as well as all intermediate) average charge distributions are asymptotically energy degenerate. Thus HF and ADF theories mistakenly break the symmetry but in a contradicting way. In this letter we show how to construct system-dependent long-range corrected (LC) density functionals that can successfully treat this class of molecules, avoiding the spurious symmetry breaking. Examples and comparisons to experimental data is given for R=H, He and Ne and it is shown that the new LC theory improves considerably the theoretical description of the R2+ bond properties, the long range form of the asymptotic potential curve as well as the atomic polarizability. The broader impact of this finding is discussed as well and it is argued that the widespread semi-empirical approach which advocates treating the LC parameter as a system-independent parameter is in fact inappropriate under general circumstances.

preprint2008arXiv

Distribution of carrier multiplication rates in CdSe and InAs nanocrystals

The distribution of rates of carrier multiplication (CM) following photon absorption is calculated for semiconductor nanocrystals (NCs). The NC electronic structure is described using a screened pseudopotential method known to give reliable description of NC excitons. The rates of biexciton generation are calculated using the Fermi golden rule with all relevant Coulomb matrix elements, taking into account proper selection rules. In CdSe and InAs NCs we find a broad distribution biexciton generation rates depending strongly on the exciton energy and size of the NC. The process becomes inefficient for NC exceeding 3 nm in diameter in the photon energy range of 2-3 times the band gap.

preprint2006arXiv

Inelastic effects in Aharonov-Bohm molecular interferometer

Inelastic effects arising from electron-phonon coupling in molecular Aharonov-Bohm (AB) interferometers are studied using the nonequilibrium Green's function method. Results for the magnetoconductance are compared for different values of the electron-phonon coupling strength. At low bias voltages, the coupling to the phonons does not change the lifetime and leads mainly to scattering phase shifts of the conducting electrons. Surprisingly, opposite to the behavior of an electrical gate, the magnetoconductance of the molecular AB interferometer becomes more sensitive to the threading magnetic flux as the electron-phonon coupling is increased.