BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Daniel Neuhauser

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

ResearcherUnclaimed source record
physics.chem-phcond-mat.mtrl-sciphysics.comp-phcond-mat.mes-hallcond-mat.othercond-mat.str-elphysics.optics

Catalog footprint

What is connected

20works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

20 published item(s)

preprint2022arXiv

Bethe Salpeter Equation Spectra for Very Large Systems

We present a highly efficient method for the extraction of optical properties of very large molecules via the Bethe-Salpeter equation. The crutch of this approach is the calculation of the action of the effective Coulombic interaction, $W$, through a stochastic TD Hartree propagation, which uses only 10 stochastic orbitals rather than propagating the full sea of occupied states. This leads to a scaling that is at most cubic in system size, with trivial MPI parallelization. We apply this new method to calculate the spectra and electronic density of the dominant excitons of a carbon-nanohoop bound fullerene system with 520 electrons, using less than 4000 core hours.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum-Spillover-Enhanced Surface-Plasmonic Absorption at the Interface of Silver and High-Index Dielectrics

We demonstrate an unexpectedly strong surface-plasmonic absorption at the interface of silver and high-index dielectrics based on electron and photon spectroscopy. The measured bandwidth and intensity of absorption deviate significantly from the classical theory. Our density-functional calculation well predicts the occurrence of this phenomenon. It reveals that due to the low metal-to-dielectric work function at such interfaces, conduction electrons can display a drastic quantum spillover, causing the interfacial electron-hole pair production to dominate the decay of surface plasmons. This finding can be of fundamental importance in understanding and designing quantum nano-plasmonic devices that utilize noble metals and high-index dielectrics.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access