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

Eran Rabani

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

ResearcherUnclaimed source record
cond-mat.mes-hallphysics.chem-phcond-mat.mtrl-scicond-mat.str-elphysics.comp-phcond-mat.stat-mechcond-mat.otherquant-phcond-mat.quant-gascond-mat.dis-nncond-mat.softOther Computer Science

Catalog footprint

What is connected

50works
12topics
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

50 published item(s)

preprint2025arXiv

Deep-learning atomistic semi-empirical pseudopotential model for nanomaterials

The semi-empirical pseudopotential method (SEPM) has been widely applied to provide computational insights into the electronic structure, photophysics, and charge carrier dynamics of nanoscale materials. We present "DeepPseudopot", a machine-learned atomistic pseudopotential model that extends the SEPM framework by combining a flexible neural network representation of the local pseudopotential with parameterized non-local and spin-orbit coupling terms. Trained on bulk quasiparticle band structures and deformation potentials from GW calculations, the model captures many-body and relativistic effects with very high accuracy across diverse semiconducting materials, as illustrated for silicon and group III-V semiconductors. DeepPseudopot's accuracy, efficiency, and transferability make it well-suited for data-driven in silico design and discovery of novel optoelectronic nanomaterials.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Simple Diagonal State Designs with Reconfigurable Real-Time Circuits

Unitary designs are widely used in quantum computation, but in many practical settings it suffices to construct a diagonal state design generated with unitary gates diagonal in the computational basis. In this work, we introduce a simple and efficient diagonal state 3-design based on real-time evolutions under 2-local Hamiltonians. Our construction is inspired by the classical Girard-Hutchinson trace estimator in that it involves the stochastic preparation of many random-phase states. Though the exact Girard-Hutchinson states are not tractably implementable on a quantum computer, we can construct states that match the statistical moments of the Girard-Hutchinson states with real-time evolution. Importantly, our random states are all generated using the same Hamiltonians for real-time evolution, with the randomness arising solely from stochastic variations in the durations of the evolutions. In this sense, the circuit is fully reconfigurable and thus suited for near-term realizations on both digital and analog platforms. Moreover, we show how to extend our construction to achieve diagonal state designs of arbitrarily high order.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Circumventing the Phonon Bottleneck by Multiphonon-Mediated Hot Exciton Cooling at the Nanoscale

In semiconductor materials, hot exciton cooling is the process by which highly excited carriers nonradiatively relax to form a band edge exciton. While cooling plays an important role in determining the thermal losses and quantum yield of a system, the timescales and mechanism of cooling are not well understood in confined semiconductor nanocrystals (NCs). A mismatch between electronic energy gaps and phonon frequencies in NCs has led to the hypothesis of a phonon bottleneck, in which cooling would be extremely slow, while enhanced electron-hole interactions in NCs have been used to explain cooling that would occur on ultrafast timescales. Here, we develop an atomistic approach for describing phonon-mediated exciton dynamics, and we use it to simulate hot exciton cooling in NCs of experimentally relevant sizes. Our framework includes electron-hole correlations as well as multiphonon processes, both of which are necessary to accurately describe the cooling process. We find that cooling occurs on timescales of tens of femtoseconds in CdSe cores, in agreement with experimental measurements, through a cascade of relaxation events that are mediated by efficient multiphonon emission. Cooling timescales increase with increasing NC size due to decreased exciton-phonon coupling (EXPC), and they are an order of magnitude larger in CdSe-CdS core-shell NCs because of reduced EXPC to low- and mid-frequency acoustic modes.

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Long-time memory effects in a localizable central spin problem

We study the properties of the Nakajima-Zwanzig memory kernel for a qubit immersed in a many-body localized (i.e., disordered and interacting) bath. We argue that the memory kernel decays as a power law in both the localized and ergodic regimes, and show how this can be leveraged to extract $t\to\infty$ populations for the qubit from finite time ($J t \leq 10^2$) data in the thermalizing phase. This allows us to quantify how the long-time values of the populations approach the expected thermalized state as the bath approaches the thermodynamic limit. This approach should provide a good complement to state-of-the-art numerical methods, for which the long-time dynamics with large baths are impossible to simulate in this phase. Additionally, our numerics on finite baths reveal the possibility for unbounded exponential growth in the memory kernel, a phenomenon rooted in the appearance of exceptional points in the projected Liouvillian governing the reduced dynamics. In small systems amenable to exact numerics, we find that these pathologies may have some correlation with delocalization.

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

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

Auger Recombination Lifetime Scaling for Type-I and Quasi-Type-II Core/Shell Quantum Dots

Having already achieved near-unity quantum yields, with promising properties for light-emitting diode, lasing, and charge separation applications, colloidal core/shell quantum dots have great technological potential. The shell thickness and band alignment of the shell and core materials are known to influence the efficiency of these devices. In many such applications, a key to improving the efficiency requires a deep understanding of multiexcitonic states. Herein, we elucidate the shell thickness and band alignment dependencies of the biexciton Auger recombination lifetime for quasi-type-II CdSe/CdS and type-I CdSe/ZnS core/shell quantum dots. We find that the biexciton Auger recombination lifetime increases with the total nanocrystal volume for quasi-type-II CdSe/CdS core/shell quantum dots and is independent of the shell thickness for type-I CdSe/ZnS core/shell quantum dots. In order to perform these calculations and compute Auger recombination lifetimes, we developed a low-scaling approach based on the stochastic resolution of identity. The numerical approach provided a framework to study the scaling of the biexciton Auger recombination lifetimes in terms of the shell thickness dependencies of the exciton radii, Coulomb couplings, and density of final states in quasi-type-II CdSe/CdS and type-I CdSe/ZnS core/shell quantum dots.

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
preprint2020arXiv

Structural relaxation in quantum supercooled liquids: A mode-coupling approach

We study supercooled dynamics in quantum hard-sphere liquid using quantum mode-coupling formulation. In the moderate quantum regime, classical cage effects lead to slower dynamics compared to strongly quantum regime, where tunneling overcomes classical caging, leading to faster relaxation. As a result, the glass transition critical density can become significantly higher than for the classical liquids. Perturbative approach is used to solve time dependent quantum mode-coupling equations to study in detail the dynamics of the supercooled liquid in moderate quantum regime. Similar to the classical case, relaxation time shows power-law increase with increasing density in the supercooled regime. However, the power-law exponent is found to be dependent on the quantumness; it increases linearly as the quantumness is increased in the moderate quantum regime.

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
preprint2015arXiv

Exact Calculation of the Time Convolutionless Master Equation Generator: Application to the Nonequilibrium Resonant Level Model

The generalized quantum master equation provides a powerful tool to describe the dynamics in quantum impurity models driven away from equilibrium. Two complementary approaches, one based on Nakajima--Zwanzig--Mori time-convolution (TC) and the other on the Tokuyama--Mori time-convolutionless (TCL) formulations provide a starting point to describe the time-evolution of the reduced density matrix. A key in both approaches is to obtain the so called "memory kernel" or "generator", going beyond second or fourth order perturbation techniques. While numerically converged techniques are available for the TC memory kernel, the canonical approach to obtain the TCL generator is based on inverting a super-operator in the \emph{full} Hilbert space, which is difficult to perform and thus, all applications of the TCL approach rely on a perturbative scheme of some sort. Here, the TCL generator is expressed using a reduced system propagator which can be obtained from system observables alone and requires the calculation of super-operators and their inverse in the \emph{reduced }Hilbert space rather than the full one. This makes the formulation amenable to quantum impurity solvers or to diagrammatic techniques, such as the nonequilibrium Green's function. We implement the TCL approach for the resonant level model driven away from equilibrium and compare the time scales for the decay of the generator with that of the memory kernel in the TC approach. Furthermore, the effects of temperature, source-drain bias, and gate potential on the TCL/TC generators are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Landauer Current and Mutual Information

We study quantum evolution of the entanglement of a quantum dot connected to left and right leads initially maintained at chemical potentials $μ_{L}$ and $μ_{R}$ respectively, within the non-interacting resonant-level model. The full nonequilirbium mixed state density matrix of the whole system is written down exactly, and entanglement is computed by recourse to the notion of mutual information. A strong and direct correlation is found between the Landauer current, and the entanglement at all times, the steady-state values in particular displaying a quadratic relationship at high temperatures. Strikingly, it is found that one can obtain a maximally entangled quantum dot by simply applying a sufficiently large `source-drain' voltage $V_{SD}$ even at high temperatures.

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

Sub-Ohmic to super-Ohmic crossover behavior in nonequilibrium quantum systems with electron-phonon interactions

The transition from weakly damped coherent motion to localization in the context of the spin-boson model has been the subject of numerous studies with distinct behavior depending on the form of the phonon-bath spectral density, $J\left(ω\right)\proptoω^{s}$. Sub-Ohmic ($s<1$) and Ohmic ($s=1$) spectral densities show a clear localization transition at zero temperature and zero bias, while for super-Ohmic ($s>1$) spectral densities this transition disappears. In this work, we consider the influence of the phonon-bath spectral density on the \emph{nonequilibrium} dynamics of a quantum dot with electron-phonon interactions described by the extended Holstein model. Using the reduced density matrix formalism combined with the multi-layer multiconfiguration time-dependent Hartree approach, we investigate the dynamic response, the time scales for relaxation, as well as the existence of multiple long-lived solutions as the system-bath coupling changes from the sub- to the super-Ohmic cases. Bistability is shown to diminish for increasing powers of $s$ similar to the spin-boson case. However, the physical mechanism and the dependence on the model parameters such as the typical bath frequency $ω_{c}$ and the polaron shift $λ$ are rather distinct.

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

A Quasi-Classical Mapping Approach to Vibrationally Coupled Electron Transport in Molecular Junctions

We develop a classical mapping approach suitable to describe vibrationally coupled charge transport in molecular junctions based on the Cartesian mapping for many-electron systems [J. Chem. Phys. 137, 154107 (2012)]. To properly describe vibrational quantum effects in the transport characteristics, we introduce a simple transformation rewriting the Hamiltonian in terms of occupation numbers and use a binning function to facilitate quantization. The approach provides accurate results for the nonequilibrium Holstein model for a range of bias voltages, vibrational frequencies and temperatures. It also captures the hallmarks of vibrational quantum effects apparent in step-like structure in the current-voltage characteristics at low temperatures as well as the phenomenon of Franck-Condon blockade.

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

Nonequilibrium quantum systems with electron-phonon interactions: Transient dynamics and approach to steady state

The nonequilibrium dynamics of a quantum dot with electron-phonon interactions described by a generalized Holstein model is presented. A combination of methodologies including the reduced density matrix formalism, the multilayer multiconfiguration time-dependent Hartree method, and a time-dependent nonequilibrium Green function approach, is used to explore the transient behavior on multiple timescales as the system approaches steady-state. The dot population dynamics on short to intermediate times is governed by the dot-lead hybridization parameter ($Γ$) and by the typical phonon frequency ($ω_{c}$) and depends on the location of the energy level of the dot relative to the bias window. At longer times, the dynamics show a distinct behavior depending on whether the system is in the adiabatic or non-adiabatic regime, with a quantum dot occupation that may depend on the initial preparation of the phonons degrees of freedom. A "phase" diagram of this localization effect as a function of the polaron shift ($λ$) for various phonon frequencies is derived, suggesting the existence of bistability on experimentally observable timescales.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phonon dynamics in correlated quantum systems driven away from equilibrium

A general form of a many-body Hamiltonian is considered, which includes an interacting fermionic sub-system coupled to non-interacting extended fermionic and bosonic systems. We show that the exact dynamics of the extended bosonic system can be derived from the reduced density matrix of the sub-system alone, despite the fact that the latter contains information about the sub-system only. The advantage of the formalism is immediately clear: While the reduced density matrix of the sub-system is readily available, the formalism offers access to observables contained in the full density matrix, which is often difficult to obtain. As an example, we consider an extended Holstein model and study the nonequilibrium dynamics of the, so called, "reaction mode" for different model parameters. The effects of the phonon frequency, the strength of the electron-phonon couplings, and the source-drain bias voltage on the phonon dynamics across the bistability are discussed.

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

Generalized projected dynamics for non-system observables of non-equilibrium quantum impurity models

The reduced dynamics formalism has recently emerged as a powerful tool to study the dynamics of non-equilibrium quantum impurity models in strongly correlated regimes. Examples include the non-equilibrium Anderson impurity model near the Kondo crossover temperature and the non-equilibrium Holstein model, for which the formalism provides an accurate description of the reduced density matrix of the system for a wide range of timescales. In this work, we generalize the formalism to allow for non-system observables such as the current between the impurity and leads. We show that the equation of motion for the reduced observable of interest can be closed with the equation of motion for the reduced density matrix and demonstrate the new formalism for a generic resonant level model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multiple steady-states in nonequilibrium quantum systems with electron-phonon interactions

The existence of more than one steady-state in a many-body quantum system driven out-of-equilibrium has been a matter of debate, both in the context of simple impurity models and in the case of inelastic tunneling channels. In this Letter, we combine a reduced density matrix formalism with the multilayer multiconfiguration time-dependent Hartree method to address this problem. This allows to obtain a converged numerical solution of the nonequilibrium dynamics. Considering a generic model for quantum transport through a quantum dot with electron-phonon interaction, we prove that a unique steady-state exists regardless of the initial electronic preparation of the quantum dot consistent with the converged numerical results. However, a bistability can be observed for different initial phonon preparations. The effects of the phonon frequency and strength of the electron-phonon couplings on the relaxation to steady-state and on the emergence of bistability is discussed.

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
preprint2013arXiv

The Electronic Structure of CdSe/CdS Core/Shell Seeded Nanorods: Type-I or Quasi-Type-II?

The electronic structure of CdSe/CdS core/shell seeded nanorods of experimentally relevant size is studied using a combination of molecular dynamics and semiempirical pseudopotential techniques, with the aim to address the transition from type-I to a quasi-type-II band alignment. The hole is found to be localized in the core region regardless of its size. The overlap of the electron density with the core region depends markedly on the size of the CdSe core: For small cores, we observe little overlap, consistent with type-II behavior. For large cores, significant core-overlap of a number of excitonic states can lead to type-I behavior. When electron-hole interactions are taken into account, the core-overlap is further increased. Our calculations indicate that the observed transition from type-II to type-I is largely due to simple volume effects, and not to band alignment.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analytical Continuation Approaches to Electronic Transport: The Resonant Level Model

The analytical continuation average spectrum method (ASM) and maximum entropy (MaxEnt) method are applied to the dynamic response of a noninteracting resonant level model within the framework of the Kubo formula for electric conductivity. The frequency dependent conductivity is inferred from the imaginary time current-current correlation function for a wide range of temperatures, gate voltages and spectral densities representing the leads, and compared with exact results. We find that the MaxEnt provides more accurate results compared to the ASM over the full spectral range.

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
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerically Exact Long Time Magnetization Dynamics at the Nonequilibrium Kondo Crossover of the Anderson Impurity Model

We investigate the dynamical and steady-state spin response of the nonequilibrium Anderson model to magnetic fields, bias voltage, and temperature using a numerically exact method combining a bold-line quantum Monte Carlo technique with the memory function formalism. We obtain converged results in a range of previously inaccessible regimes, in particular the crossover to the Kondo domain. We provide detailed predictions for novel nonequilibrium phenomena, including non-monotonic temperature dependence of observables at high bias voltage and oscillatory quench dynamics at high magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Periodic negative differential conductance in a single metallic nano-cage

We report a bi-polar multiple periodic negative differential conductance (NDC) effect on a single cage-shaped Ru nanoparticle measured using scanning tunneling spectroscopy. This phenomenon is assigned to the unique multiply-connected cage architecture providing two (or more) defined routes for charge flow through the cage. This, in turn, promotes a self- gating effect, where electron charging of one route affects charge transport along a neighboring channel, yielding a series of periodic NDC peaks. This picture is established and analyzed here by a theoretical model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulating Lattice Spin Models on Graphics Processing Units

Lattice spin models are useful for studying critical phenomena and allow the extraction of equilibrium and dynamical properties. Simulations of such systems are usually based on Monte Carlo (MC) techniques, and the main difficulty is often the large computational effort needed when approaching critical points. In this work, it is shown how such simulations can be accelerated with the use of NVIDIA graphics processing units (GPUs) using the CUDA programming architecture. We have developed two different algorithms for lattice spin models, the first useful for equilibrium properties near a second-order phase transition point and the second for dynamical slowing down near a glass transition. The algorithms are based on parallel MC techniques, and speedups from 70- to 150-fold over conventional single-threaded computer codes are obtained using consumer-grade hardware.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Steady state conductance in a double quantum dot array: The nonequilibrium equation-of-motion Green function approach

We study steady state transport through a double quantum dot array using the equation-of-motion approach to the nonequilibrium Green functions formalism. This popular technique relies on uncontrolled approximations to obtain a closure for a hierarchy of equations, however its accuracy is questioned. We focus on 4 different closures, 2 of which were previously proposed in the context of the single quantum dot system (Anderson impurity model) and were extended to the double quantum dot array, and develop 2 new closures. Results for the differential conductance are compared to those attained by a master equation approach known to be accurate for weak system-leads couplings and high temperatures. While all 4 closures provide an accurate description of the Coulomb blockade and other transport properties in the single quantum dot case, they differ in the case of the double quantum dot array, where only one of the developed closures provides satisfactory results.This is rationalized by comparing the poles of the Green functions to the exact many-particle energy differences for the isolate system. Our analysis provides means to extend the equation-of-motion technique to more elaborate models of large bridge systems with strong electronic interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Symmetry breaking and restoration using the equation-of-motion technique for nonequilibrium quantum impurity models

The description of the dynamics of correlated electrons in quantum impurity models is typically described within the nonequilibrium Green function formalism combined with a suitable approximation. One common approach is based on the equation-of-motion technique often used to describe different regimes of the dynamic response. Here, we show that this approach may violate certain symmetry relations that must be fulfilled by the definition of the Green functions. These broken symmetries can lead to unphysical behavior. To circumvent this pathological shortcoming of the equation-of-motion approach we provide a scheme to restore basic symmetry relations. Illustrations are given for the Anderson and double Anderson impurity models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transferable Pair Potentials for CdS and ZnS Crystals

A set of interatomic pair potentials is developed for CdS and ZnS crystals. We show that a simple energy function, which has been used to describe the properties of CdSe [J. Chem. Phys. 116, 258 (2002)], can be parametrized to accurately describe the lattice and elastic constants, and phonon dispersion relations of bulk CdS and ZnS in the wurtzite and rocksalt crystal structures. The predicted coexistence pressure of the wurtzite and rocksalt structures, as well as the equation of state are in good agreement with experimental observations. These new pair potentials enable the study of a wide range of processes in bulk and nanocrystalline II-VI semiconductor materials.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Application of a semiclassical model for the second-quantized many-electron Hamiltonian to nonequilibrium quantum transport: The resonant level model

A semiclassical (SC) approach is developed for nonequilibrium quantum transport in molecular junctions. Following the early work of Miller and White [J. Chem. Phys. 84, 5059 (1986)], the many-electron Hamiltonian in second quantization is mapped onto a classical model that preserves the fermionic character of electrons. The resulting classical electronic Hamiltonian allows for real-time molecular dynamics simulations of the many-body problem from an uncorrelated initial state to the steady state. Comparisons with exact results generated for the resonant level model reveal that a semiclassical treatment of transport provides a quantitative description of the dynamics at all relevant timescales for a wide range of bias and gate potentials, and for different temperatures. The approach opens a door to treating nontrivial quantum transport problems that remain far from the reach of fully quantum methodologies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Memory Effects In Nonequilibrium Quantum Impurity Models

Memory effects play a key role in the dynamics of strongly correlated systems driven out of equilibrium. In the present study, we explore the nature of memory in the nonequilibrium Anderson impurity model. The Nakajima--Zwanzig--Mori formalism is used to derive an exact generalized quantum master equation for the reduced density matrix of the interacting quantum dot, which includes a non-Markovian memory kernel. A real-time path integral formulation is developed, in which all diagrams are stochastically sampled in order to numerically evaluate the memory kernel. We explore the effects of temperature down to the Kondo regime, as well as the role of source--drain bias voltage and band width on the memory. Typically, the memory decays on timescales significantly shorter than the dynamics of the reduced density matrix itself, yet under certain conditions it develops a smaller long tail. In addition we address the conditions required for the existence, uniqueness and stability of a steady-state.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Theory and simulations of quantum glass forming liquids

A comprehensive microscopic dynamical theory is presented for the description of quantum fluids as they transform into glasses. The theory is based on a quantum extension of mode-coupling theory. Novel effects are predicted, such as reentrant behavior of dynamical relaxation times. These predictions are supported by path integral ring polymer molecular dynamics simulations. The simulations provide detailed insight into the factors that govern slow dynamics in glassy quantum fluids. Connection to other recent work on both quantum glasses as well as quantum optimization problems is presented.

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum fluctuations can promote or inhibit glass formation

The very nature of glass is somewhat mysterious: while relaxation times in glasses are of sufficient magnitude that large-scale motion on the atomic level is essentially as slow as it is in the crystalline state, the structure of glass appears barely different than that of the liquid that produced it. Quantum mechanical systems ranging from electron liquids to superfluid helium appear to form glasses, but as yet no unifying framework exists connecting classical and quantum regimes of vitrification. Here we develop new insights from theory and simulation into the quantum glass transition that surprisingly reveal distinct regions where quantum fluctuations can either promote or inhibit glass formation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Analytic continuation average spectrum method for quantum liquids

We revisit the problem of determining the real-frequency density response in quantum fluids via analytical continuation of imaginary-time quantum Monte Carlo data. We demonstrate that the average spectrum method (ASM) is capable of revealing resolved modes in the dynamic structure factor of both ortho-deuterium and liquid para-hydrogen, in agreement with experiments and quantum mode-coupling theories, while the maximum entropy approach yields only a smooth unimodal spectrum. Outstanding issues are discussed. Our work provides the first application of the ASM method in realistic off-lattice systems.

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

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

0 citations0 reviewsNo code linkedOpen access