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

David R. Reichman

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

ResearcherUnclaimed source record
cond-mat.stat-mechphysics.chem-phcond-mat.mes-hallcond-mat.str-elcond-mat.mtrl-scicond-mat.dis-nncond-mat.softquant-phcond-mat.quant-gasphysics.comp-phcond-mat.otherMachine Learningphysics.soc-phSocial and Information Networks

Catalog footprint

What is connected

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

68 published item(s)

preprint2023arXiv

Microscopic observation of two-level systems in a metallic glass model

The low-temperature quasi-universal behavior of amorphous solids has been attributed to the existence of spatially-localized tunneling defects found in the low-energy regions of the potential energy landscape. Computational models of glasses can be studied to elucidate the microscopic nature of these defects. Recent simulation work has demonstrated the means of generating stable glassy configurations for models that mimic metallic glasses using the swap Monte Carlo algorithm. Building on these studies, we present an extensive exploration of the glassy metabasins of the potential energy landscape of a variant of the most widely used model of metallic glasses. We carefully identify tunneling defects and reveal their depletion with increased glass stability. The density of tunneling defects near the experimental glass transition temperature appears to be in good agreement with experimental measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Localized-Orbital Energy Evaluation for Auxiliary-Field Quantum Monte Carlo

Phaseless Auxiliary-Field Quantum Monte Carlo (ph-AFQMC) has recently emerged as a promising method for the production of benchmark-level simulations of medium to large-sized molecules, due to its accuracy and favorable polynomial scaling with system size. Unfortunately the memory footprint of standard energy evaluation algorithms are non-trivial, which can significantly impact timings on graphical processing units (GPUs) where memory is limited. Previous attempts to reduce scaling by taking advantage of the low rank structure of the Coulombic integrals have been successful, but are significantly limited by high prefactors, rendering the utility limited to very large systems. Here, we present a complementary, cubic scaling route to reduce memory and computational scaling based on the low rank of the Coulombic interactions between localized orbitals, focusing on the application to phaseless AFQMC. We show that the error due to this approximation, which we term Localized Orbital AFQMC (LO-AFQMC), is systematic and controllable via a single variable, and is computationally favorable even for small systems. We present results demonstrating a robust retention of accuracy versus both experiment and full ph-AFQMC for a variety of test cases chosen for their potential difficulty for localized orbital based methods, including the singlet-triplet gaps of polyacenes benzene through pentacene, the heats of formation for a set of platonic hydrocarbon cages, and the total energy of ferrocene (Fe(Cp)$_2$). Finally, we reproduce our previous result of the gas phase ionization energy of Ni(Cp)$_2$, agreeing with full ph-AFQMC to within statistical error while using less than a fifteenth of the computer time.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cumulant methods for electron-phonon problems. I. Perturbative expansions

In this work we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CE for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as well as the approach to the thermodynamic limit, drawing distinctions and connections between the behavior of systems in and away from the thermodynamic limit that enable a greater understanding of the ability of the CE to capture real-frequency information across the full range of wave vectors. We find that for zero electronic momentum, the spectral function is well described by the second-order CE at low and high temperatures. However, for non-zero electronic momenta, the CE is only accurate at high temperature. We analyze the fourth-order cumulant, and find that while it improves the description of the short-time dynamics encoded in the one-particle Green's function, it can introduce divergences in the time domain as well as unphysical negative spectral weight in the spectral function. When well-behaved, the fourth-order CE does provide notable accurate corrections to the second-order CE. Finally, we use our results to comment on the use of the CE as a tool for calculating transport behavior in the realistic ab initio modeling of materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cumulant methods for electron-phonon problems. II. The self-consistent cumulant expansion

In this work we present a self-consistent cumulant expansion (SC-CE) and investigate its accuracy for the one-dimensional Holstein model with and without phonon dispersion. We show that for finite lattices sizes, the numerical integration of the SC-CE equations becomes unstable at long times. This defect is partially ameliorated when studying systems in the thermodynamic limit, enabling the demonstration that the SC-CE corrects many deficits of the standard perturbative CE in the (non-dispersive) Holstein model. The natural phonon damping that arises in the more realistic dispersed Holstein model renders the SC-CE stable, allowing for a complete assessment of the method. Here we find that self-consistency dramatically corrects many of the failures found in the perturbative CE, but also introduces some unphysical features. Finally, we comment on the potential use of SC-CE as a tool for calculating Green's functions in generic many-body problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Response to "Exponential challenges in unbiasing quantum Monte Carlo algorithms with quantum computers"

A recent preprint by Mazzola and Carleo numerically investigates exponential challenges that can arise for the QC-QMC algorithm introduced in our work, "Unbiasing fermionic quantum Monte Carlo with a quantum computer." As discussed in our original paper, we agree with this general concern. However, here we provide further details and numerics to emphasize that the prospects for practical quantum advantage in QC-QMC remain open. The exponential challenges in QC-QMC are dependent on (1) the choice of QMC methods, (2) the underlying system, and (3) the form of trial and walker wavefunctions. While one can find difficult examples with a specific method, a specific system, and a specific walker/trial form, for some combinations of these choices, the approach is potentially more scalable than other near-term quantum algorithms. Future research should aim to identify examples for which QC-QMC enables practical quantum advantage.

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

Strongly correlated ladders in K-doped $p$-terphenyl crystals

Potassium-doped terphenyl has recently attracted attention as a potential host for high-transition-temperature superconductivity. Here, we elucidate the many-body electronic structure of recently synthesized potassium-doped terphenyl crystals. We show that this system may be understood as a set of weakly coupled one-dimensional ladders. Depending on the strength of the inter-ladder coupling the system may exhibit spin-gapped valence-bond solid or antiferromagnetic phases, both of which upon hole doping may give rise to superconductivity. This terphenyl-based ladder material serves as a new platform for investigating the fate of ladder phases in presence of three-dimensional coupling as well as for novel superconductivity.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Unbiasing Fermionic Quantum Monte Carlo with a Quantum Computer

Many-electron problems pose some of the greatest challenges in computational science, with important applications across many fields of modern science. Fermionic quantum Monte Carlo (QMC) methods are among the most powerful approaches to these problems. However, they can be severely biased when controlling the fermionic sign problem using constraints, as is necessary for scalability. Here we propose an approach that combines constrained QMC with quantum computing tools to reduce such biases. We experimentally implement our scheme using up to 16 qubits in order to unbias constrained QMC calculations performed on chemical systems with as many as 120 orbitals. These experiments represent the largest chemistry simulations performed on quantum computers (more than doubling the size of prior electron correlation calculations), while obtaining accuracy competitive with state-of-the-art classical methods. Our results demonstrate a new paradigm of hybrid quantum-classical algorithm, surpassing the popular variational quantum eigensolver in terms of potential towards the first practical quantum advantage in ground state many-electron calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Attractive versus truncated repulsive supercooled liquids: The dynamics is encoded in the pair correlation function

We compare glassy dynamics in two liquids that differ in the form of their interaction potentials. Both systems have the same repulsive interactions but one has also an attractive part in the potential. These two systems exhibit very different dynamics despite having nearly identical pair correlation functions. We demonstrate that a properly weighted integral of the pair correlation function, which amplifies the subtle differences between the two systems, correctly captures their dynamical differences. The weights are obtained from a standard machine learning algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Depletion of two-level systems in ultrastable computer-generated glasses

Amorphous solids exhibit quasi-universal low-temperature anomalies whose origin has been ascribed to localized tunneling defects. Using an advanced Monte Carlo procedure, we create {\it in silico} glasses spanning from hyperquenched to ultrastable glasses. Using a multidimensional path-finding protocol, we locate tunneling defects with energy splittings smaller than $k_{B}T_Q$, with $T_Q$ the temperature below which quantum effects are relevant ($T_Q \approx 1 \,$K in most experiments). We find that as the stability of a glass increases, its energy landscape as well as the manner in which it is probed tend to deplete the density of tunneling defects, as observed in recent experiments. We explore the real-space nature of tunneling defects, and find that they are mostly localized to a few atoms, but are occasionally dramatically delocalized.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Microscopic model of the doping dependence of line widths in monolayer transition metal dichalcogenides

A fully microscopic model of the doping-dependent exciton and trion line widths in the absorption spectra of monolayer transition metal dichalcogenides in the low temperature and low doping regime is explored. The approach is based on perturbation theory and avoids the use of phenomenological parameters. In the low-doping regime, we find that the trion line width is relatively insensitive to doping levels while the exciton line width increases monotonically with doping. On the other hand, we argue that the trion line width shows a somewhat stronger temperature dependence. The magnitudes of the line widths are likely to be masked by phonon scattering for $T \geq 20$ K in encapsulated samples in the low doping regime. We discuss the breakdown of perturbation theory, which should occur at relatively low doping levels and low temperatures. Our work also paves the way towards understanding a variety of related scattering processes, including impact ionization and Auger scattering in clean 2D samples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo

Transition metal complexes are ubiquitous in biology and chemical catalysis, yet they remain difficult to accurately describe with ab initio methods due to the presence of a large degree of dynamic electron correlation, and, in some cases, strong static correlation which results from a manifold of low-lying states. Progress has been hindered by a scarcity of high quality gas-phase experimental data, while exact ab initio predictions are usually computationally unaffordable due to the large size of the systems. In this work, we present a data set of 34 3d metal-containing complexes with gas-phase ligand-dissociation energies that have reported uncertainties of $\leq$ 2 kcal/mol. We perform all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multi-determinant trial wavefunctions selected by a blackbox procedure. We compare the results with those from DFT with various functionals, and DLPNO-CCSD(T). We find MAE of 1.09 $\pm$ 0.28 kcal/mol for our best ph-AFQMC method, vs 2.89 kcal/mol for DLPNO-CCSD(T) and 1.57 - 3.87 kcal/mol for DFT. We find maximum errors of 2.96 $\pm$ 1.71 kcal/mol for our best ph-AFQMC method, vs 9.15 kcal/mol for DLPNO-CCSD(T) and 5.98 - 13.69 kcal/mol for DFT. The reasonable performance of several functionals is in stark contrast to the much poorer accuracy previously demonstrated for diatomics, suggesting a moderation in electron correlation due to ligand coordination. However, the unpredictably large errors for a small subset of cases with both DFT and DLPNO-CCSD(T) leave cause for concern, especially due to the unreliability of common multi-reference indicators. In contrast, the robust and, in principle, systematically improvable results of ph-AFQMC for these realistic complexes establish it as a useful tool for elucidating the electronic structure of transition metal-containing complexes and predicting their gas-phase properties.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stochastic Resolution-of-the-Identity Auxiliary-Field Quantum Monte Carlo: Scaling Reduction without Overhead

We explore the use of the stochastic resolution-of-the-identity (sRI) with the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) method. sRI is combined with four existing local energy evaluation strategies in ph-AFQMC, namely (1) the half-rotated electron repulsion integral tensor (HR), (2) Cholesky decomposition (CD), (3) tensor hypercontraction (THC), or (4) low-rank factorization (LR). We demonstrate that HR-sRI achieves no scaling reduction, CD-sRI scales as $\mathcal O(N^3)$, and THC-sRI and LR-sRI scale as $\mathcal O(N^2)$, albeit with a potentially large prefactor. Furthermore, the walker-specific extra memory requirement in CD is reduced from $\mathcal O(N^3)$ to $\mathcal O(N^2)$ with sRI, while sRI-based THC and LR algorithms lead to a reduction from $\mathcal O(N^2)$ extra memory to $\mathcal O(N)$. Based on numerical results for one-dimensional hydrogen chains and water clusters, we demonstrated that, along with the use of a variance reduction technique, CD-sRI achieves cubic-scaling {\it without overhead}. In particular, we find for the systems studied the observed scaling of standard CD is $\mathcal O(N^{3-4})$ while for CD-sRI it is reduced to $\mathcal O(N^{2-3})$. Once a memory bottleneck is reached, we expect THC-sRI and LR-sRI to be preferred methods due to their quadratic-scaling memory requirements and their quadratic-scaling of the local energy evaluation (with a potentially large prefactor). The theoretical framework developed here should facilitate large-scale ph-AFQMC applications that were previously difficult or impossible to carry out with standard computational resources.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The performance of phaseless auxiliary-field quantum Monte Carlo on the ground state electronic energy of benzene

The ground state electronic energy of benzene is the focus of a recent blind test by Eriksen and co-workers [arXiv:2008.02678]. In this note, we report the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) results on the identical problem. We examine trial wavefunctions of the spin-restricted Hartree-Fock (RHF), and complete active space self-consistent field (CASSCF) with an active space of 6-electron and 6-orbital (i.e., CAS(6,6)) form for the cc-pVDZ basis set. ph-AFQMC+RHF deviates from the value where many methods agreed on (i.e., -863 m$E_h$) by -3.1(3) m$E_h$ whereas ph-AFQMC+CAS(6,6) deviates from the same value by -1.3(4) m$E_h$. In addition to this, we report the frozen core correlation energy of ph-AFQMC+RHF in the cc-pVTZ and cc-pVQZ basis sets as well as their complete basis set limit. Our findings highlight the accuracy, flexibility, and scalability of ph-AFQMC with simple trial wavefunctions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Description of quasiparticle and satellite properties via cumulant expansions of the retarded one-particle Green's function

A new cumulant-based $GW$ approximation for the retarded one-particle Green's function is proposed, motivated by an exact relation between the improper Dyson self-energy and the cumulant generating function. Qualitative aspects of this method are explored within a simple one-electron independent phonon model, where it is seen that the method preserves the energy moment of the spectral weight while also reproducing the exact Green's function in the weak coupling limit. For the three-dimensional electron gas, this method predicts multiple satellites at the bottom of the band, albeit with inaccurate peak spacing. However, its quasiparticle properties and correlation energies are more accurate than both previous cumulant methods and standard $G_0W_0$. Our results point to new features that may be exploited within the framework of cumulant-based methods and suggest promising directions for future exploration and improvement of cumulant-based $GW$ approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Many-Body Localization in System with a Completely Delocalized Single-Particle Spectrum

Many-body localization (MBL) in a one-dimensional Fermi Hubbard model with random on-site interactions is studied. While for this model all single-particle states are trivially delocalized, it is shown that for sufficiently strong disordered interactions the model is many-body localized. It is therefore argued that MBL does not necessary rely on localization of the single-particle spectrum. This model provides a convenient platform to study pure MBL phenomenology, since Anderson localization in this model does not exist. By examining various forms of the interaction term a dramatic effect of symmetries on charge transport is demonstrated. A possible realization in a cold atom experiments is proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the accuracy of surface hopping dynamics in condensed phase non-adiabatic problems

We perform extensive benchmark comparisons of surface hopping dynamics with numerically exact calculations for the spin-boson model over a wide range of energetic and coupling parameters as well as temperature. We find that deviations from golden-rule scaling in the Marcus regime are generally small and depend sensitively on the energetic bias between electronic states. Fewest switched surface hopping (FSSH) is found to be surprisingly accurate over a large swath of parameter space. The inclusion of decoherence corrections via the augmented FSSH (A-FSSH) algorithm improves the accuracy of dynamical behavior compared to exact simulations, but the effects are generally not dramatic, at least for the case of an environment modeled with the commonly used Debye spectral density.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the accuracy of the Pade-resummed master equation approach to dissipative quantum dynamics

Well--defined criteria are proposed for assessing the accuracy of quantum master equations whose memory functions are approximated by Padé resummation of the first two moments in the electronic coupling. These criteria partition the parameter space into distinct levels of expected accuracy, ranging from quantitatively accurate regimes to regions of parameter space where the approach is not expected to be applicable. Extensive comparison of Padé--resummed master equations with numerically exact results in the context of the spin--boson model demonstrate that the proposed criteria correctly demarcate the regions of parameter space where the Padé approximation is reliable. The applicability analysis we present is not confined to the specifics of the Hamiltonian under consideration and should provide guidelines for other classes of resummation techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Slow Dynamics in a Two-Dimensional Anderson-Hubbard Model

We study the real-time dynamics of a two-dimensional Anderson--Hubbard model using nonequilibrium self-consistent perturbation theory within the second-Born approximation. When compared with exact diagonalization performed on small clusters, we demonstrate that for strong disorder this technique approaches the exact result on all available timescales, while for intermediate disorder, in the vicinity of the many-body localization transition, it produces quantitatively accurate results up to nontrivial times. Our method allows for the treatment of system sizes inaccessible by any numerically exact method and for the complete elimination of finite size effects for the times considered. We show that for a sufficiently strong disorder the system becomes nonergodic, while for intermediate disorder strengths and for all accessible time scales transport in the system is strictly subdiffusive. We argue that these results are incompatible with a simple percolation picture, but are consistent with the heuristic random resistor network model where subdiffusion may be observed for long times until a crossover to diffusion occurs. The prediction of slow finite-time dynamics in a two-dimensional interacting and disordered system can be directly verified in future cold atoms experiments

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Anderson--Holstein Model in Two Flavors of the Non--Crossing Approximation

The dynamical interplay between electron-electron interactions and electron-phonon coupling is investigated within the Anderson-Holstein model, a minimal model for open quantum systems that embody these effects. The influence of phonons on spectral and transport properties is explored in equilibrium, for non-equilibrium steady state and for transient dynamics after a quench. Both the particle-hole symmetric and the more generic particle-hole asymmetric cases are studied. The treatment is based on two complementary non-crossing approximations, the first of which is constructed around the weak-coupling limit and the second around the polaron limit. In general, the two methods disagree in nontrivial ways, indicating that more reliable approaches to the problem are needed. The frameworks used here can form the starting point for numerically exact methods based on bold-line continuous-time quantum Monte Carlo algorithms capable of treating open systems simultaneously coupled to multiple fermionic and bosonic baths.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Absence of diffusion in an interacting system of spinless fermions on a one-dimensional disordered lattice

We study the infinite temperature dynamics of a prototypical one-dimensional system expected to exhibit many-body localization. Using numerically exact methods, we establish the dynamical phase diagram of this system based on the statistics of its eigenvalues and its dynamical behavior. We show that the nonergodic phase is reentrant as a function of the interaction strength, illustrating that localization can be reinforced by sufficiently strong interactions even at infinite temperature. Surprisingly, within the accessible time range, the ergodic phase shows subdiffusive behavior, suggesting that the diffusion coefficient vanishes throughout much of the phase diagram in the thermodynamic limit. Our findings strongly suggest that Wigner-Dyson statistics of eigenvalue spacings may appear in a class of ergodic but subdiffusive systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Binding energies and spatial structures of small carrier complexes in monolayer transition metal dichalcogenides via diffusion Monte Carlo

Ground state diffusion Monte Carlo is used to investigate the binding energies and carrier probability distributions of excitons, trions, and biexcitons in a variety of two-dimensional transition metal dichalcogenide materials. We compare these results to approximate variational calculations, as well as to analogous Monte Carlo calculations performed with simplified carrier interaction potentials. Our results highlight the successes and failures of approximate approaches as well as the physical features that determine the stability of small carrier complexes in monolayer transition metal dichalcogenide materials. Lastly, we discuss points of agreement and disagreement with recent experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bright and dark singlet excitons via linear and two-photon spectroscopy in monolayer transition-metal dichalcogenides

We discuss the linear and two-photon spectroscopic selection rules for spin-singlet excitons in monolayer transition-metal dichalcogenides. Our microscopic formalism combines a fully $k$-dependent few-orbital band structure with a many-body Bethe-Salpeter equation treatment of the electron-hole interaction, using a model dielectric function. We show analytically and numerically that the single-particle, valley-dependent selection rules are preserved in the presence of excitonic effects. Furthermore, we definitively demonstrate that the bright (one-photon allowed) excitons have $s$-type azimuthal symmetry and that dark $p$-type excitons can be probed via two-photon spectroscopy. The screened Coulomb interaction in these materials substantially deviates from the $1/\varepsilon_0 r$ form; this breaks the "accidental" angular momentum degeneracy in the exciton spectrum, such that the 2$p$ exciton has a lower energy than the 2$s$ exciton by at least 50 meV. We compare our calculated two-photon absorption spectra to recent experimental measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Extending the applicability of Redfield theories into highly non-Markovian regimes

We present a new, computationally inexpensive method for the calculation of reduced density matrix dynamics for systems with a potentially large number of subsystem degrees of freedom coupled to a generic bath. The approach consists of propagation of weak-coupling Redfield-like equations for the high frequency bath degrees of freedom only, while the low frequency bath modes are dynamically arrested but statistically sampled. We examine the improvements afforded by this approximation by comparing with exact results for the spin-boson model over a wide range of parameter space. The results from the method are found to dramatically improve Redfield dynamics in highly non--Markovian regimes, at a similar computational cost. Relaxation of the mode-freezing approximation via classical (Ehrenfest) evolution of the low frequency modes results in a dynamical hybrid method. We find that this Redfield-based dynamical hybrid approach, which is computationally more expensive than bare Redfield dynamics, yields only a marginal improvement over the simpler approximation of complete mode arrest.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generalized mode-coupling theory of the glass transition: schematic results at finite and infinite order

We present an extensive treatment of the generalized mode-coupling theory (GMCT) of the glass transition, which seeks to describe the dynamics of glass-forming liquids using only static structural information as input. This theory amounts to an infinite hierarchy of coupled equations for multi-point density correlations, the lowest-order closure of which is equivalent to standard mode-coupling theory. Here we focus on simplified schematic GMCT hierarchies, which lack any explicit wavevector-dependence and therefore allow for greater analytical and numerical tractability. For one particular schematic model, we derive the unique analytic solution of the infinite hierarchy, and demonstrate that closing the hierarchy at finite order leads to uniform convergence as the closure level increases. We also show numerically that a similarly robust convergence pattern emerges for more generic schematic GMCT models, suggesting that the GMCT framework is generally convergent, even though no small parameter exists in the theory. Finally, we discuss how different effective weights on the high-order contributions ultimately control whether the transition is continuous, discontinuous, or strictly avoided, providing new means to relate structure to dynamics in glass-forming systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Microscopic dynamics of supercooled liquids from first principles

Glasses are solid materials whose constituent atoms are arranged in a disordered manner. The transition from a liquid to a glass remains one of the most poorly understood phenomena in condensed matter physics, and still no fully microscopic theory exists that can describe the dynamics of supercooled liquids in a quantitative manner over all relevant time scales. Here we present such a theoretical framework that yields near-quantitative accuracy for the time-dependent correlation functions of a supercooled system over a broad density range. Our approach requires only simple static structural information as input and is based entirely based on first principles. Owing to this first-principles nature, the framework offers a unique platform to study the relation between structure and dynamics in glass-forming matter, and paves the way towards a systematically correctable and ultimately fully quantitative theory of microscopic glassy dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Algorithm independent bounds on community detection problems and associated transitions in stochastic block model graphs

We derive rigorous bounds for well-defined community structure in complex networks for a stochastic block model (SBM) benchmark. In particular, we analyze the effect of inter-community "noise" (inter-community edges) on any "community detection" algorithm's ability to correctly group nodes assigned to a planted partition, a problem which has been proven to be NP complete in a standard rendition. Our result does not rely on the use of any one particular algorithm nor on the analysis of the limitations of inference. Rather, we turn the problem on its head and work backwards to examine when, in the first place, well defined structure may exist in SBMs.The method that we introduce here could potentially be applied to other computational problems. The objective of community detection algorithms is to partition a given network into optimally disjoint subgraphs (or communities). Similar to k-SAT and other combinatorial optimization problems, "community detection" exhibits different phases. Networks that lie in the "unsolvable phase" lack well-defined structure and thus have no partition that is meaningful. Solvable systems splinter into two disparate phases: those in the "hard" phase and those in the "easy" phase. As befits its name, within the easy phase, a partition is easy to achieve by known algorithms. When a network lies in the hard phase, it still has an underlying structure yet finding a meaningful partition which can be checked in polynomial time requires an exhaustive computational effort that rapidly increases with the size of the graph. When taken together, (i) the rigorous results that we report here on when graphs have an underlying structure and (ii) recent results concerning the limits of rather general algorithms, suggest bounds on the hard phase.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Correlation of Local Order with Particle Mobility in Supercooled Liquids is Highly System Dependent

We investigate the connection between local structure and dynamical heterogeneity in supercooled liquids. Through the study of four different models we show that the correlation between a particle's mobility and the degree of local order in nearby regions is highly system dependent. Our results suggest that the correlation between local structure and dynamics is weak or absent in systems that conform well to the mean-field picture of glassy dynamics and strong in those that deviate from this paradigm. Finally, we investigate the role of order-agnostic point-to-set correlations and reveal that they provide similar information content to local structure measures, at least in the system where local order is most pronounced.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Critical dynamical heterogeneities close to continuous second-order glass transitions

We analyse, using Inhomogenous Mode-Coupling Theory, the critical scaling behaviour of the dynamical susceptibility at a distance epsilon from continuous second-order glass transitions. We find that the dynamical correlation length xi behaves generically as epsilon^{-1/3} and that the upper critical dimension is equal to six. More surprisingly, we find activated dynamic scaling, where xi grows with time as [ln(t)]^2 exactly at criticality. All these results suggest a deep analogy between the glassy behaviour of attractive colloids or randomly pinned supercooled liquids and that of the Random Field Ising Model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Crossovers in the dynamics of supercooled liquids probed by an amorphous wall

We study the relaxation dynamics of a binary Lennard-Jones liquid in the presence of an amorphous wall generated from equilibrium particle configurations. In qualitative agreement with the results presented in Nature Phys. {\bf 8}, 164 (2012) for a liquid of harmonic spheres, we find that our binary mixture shows a saturation of the dynamical length scale close to the mode-coupling temperature $T_c$. Furthermore we show that, due to the broken symmetry imposed by the wall, signatures of an additional change in dynamics become apparent at a temperature well above $T_c$. We provide evidence that this modification in the relaxation dynamics occurs at a recently proposed dynamical crossover temperature $T_s > T_c$, which is related to the breakdown of the Stokes-Einstein relation. We find that this dynamical crossover at $T_s$ is also observed for a system of harmonic spheres as well as a WCA liquid, showing that it may be a general feature of glass-forming systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamics of Many-Body Localization

Following the field theoretic approach of Basko et al., Ann. Phys. 321, 1126 (2006), we study in detail the real-time dynamics of a system expected to exhibit many-body localization. In particular, for time scales inaccessible to exact methods, we demonstrate that within the second-Born approximation that the temporal decay of the density-density correlation function is non-exponential and is consistent with a finite value for $t\to\infty$, as expected in a non-ergodic state. This behavior persists over a wide range of disorder and interaction strengths. We discuss the implications of our findings with respect to dynamical phase boundaries based both on exact diagonalization studies and as well as those established by the methods of Ref. 1.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equilibrium ultrastable glasses produced by random pinning

Ultrastable glasses have risen to prominence due to their potentially useful material properties and the tantalizing possibility of a general method of preparation via vapor deposition. Despite the importance of this novel class of amorphous materials, numerical studies have been scarce because achieving ultrastability in atomistic simulations is an enormous challenge. Here we bypass this difficulty and establish that randomly pinning the position of a small fraction of particles inside an equilibrated supercooled liquid generates ultrastable configurations at essentially no numerical cost, while avoiding undesired structural changes due to the preparation protocol. Building on the analogy with vapor-deposited ultrastable glasses, we study the melting kinetics of these configurations following a sudden temperature jump into the liquid phase. In homogeneous geometries, we find that enhanced kinetic stability is accompanied by large scale dynamic heterogeneity, while a competition between homogeneous and heterogeneous melting is observed when a liquid boundary invades the glass at constant velocity. Our work demonstrates the feasibility of large-scale, atomistically resolved, and experimentally relevant simulations of the kinetics of ultrastable glasses.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exciton Binding Energy and Nonhydrogenic Rydberg Series in Monolayer WS2

We have experimentally determined the energies of the ground and first four excited excitonic states of the fundamental optical transition in monolayer WS2, a model system for the growing class of atomically thin two-dimensional semiconductor crystals. From the spectra, we establish a large exciton binding energy of 0.32 eV and a pronounced deviation from the usual hydrogenic Rydberg series of energy levels of the excitonic states. We explain both of these results using a microscopic theory in which the non-local nature of the effective dielectric screening modifies the functional form of the Coulomb interaction. These strong but unconventional electron-hole interactions are expected to be ubiquitous in atomically thin materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Green's functions from real-time bold-line Monte Carlo

We present two methods for computing two-time correlation functions or Green's functions from real time bold-line continuous time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real time methods to analytically continued spectral functions obtained from imaginary time Monte Carlo, thus probing the limits of analytic continuation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Green's functions from real-time bold-line Monte Carlo: spectral properties of the nonequilibrium Anderson impurity model

The nonequilibrium spectral properties of the Anderson impurity model with a chemical potential bias are investigated within a numerically exact real time quantum Monte Carlo formalism. The two-time correlation function is computed in a form suitable for nonequilibrium dynamical mean field calculations. Additionally, the evolution of the model's spectral properties are simulated in an alternative representation, defined by a hypothetical but experimentally realizable weakly coupled auxiliary lead. The voltage splitting of the Kondo peak is confirmed and the dynamics of its formation after a coupling or gate quench are studied. This representation is shown to contain additional information about the dot's population dynamics. Further, we show that the voltage-dependent differential conductance gives a reasonable qualitative estimate of the equilibrium spectral function, but significant qualitative differences are found including incorrect trends and spurious temperature dependent effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relaxation Patterns in Supercooled Liquids from Generalized Mode-Coupling Theory

The mode-coupling theory of the glass transition treats the dynamics of supercooled liquids in terms of two-point density correlation functions. Here we consider a generalized, hierarchical formulation of schematic mode-coupling equations in which the full basis of multipoint density correlations is taken into account. By varying the parameters that control the effective contributions of higher-order correlations, we show that infinite hierarchies can give rise to both sharp and avoided glass transitions. Moreover, small changes in the form of the coefficients result in different scaling behaviors of the structural relaxation time, providing a means to tune the fragility in glass-forming materials. This demonstrates that the infinite-order construct of generalized mode-coupling theory constitutes a powerful and unifying framework for kinetic theories of the glass transition.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A small subset of normal modes mimics the properties of dynamical heterogeneity in a model supercooled liquid

In this work, we study the nature of transitions between inherent structures of a two-dimensional model supercooled liquid. We demonstrate that these transitions occur predominately along a small number of directions on the energy landscape. Moreover, we show that the number of such directions decreases as the temperature of the liquid is decreased in the mildly supercooled regime, in concert with earlier studies on an athermal jamming system. We show that this decrease happens in parallel with a change in character of the transitions as dynamics in the system become more heterogeneous and localized. We investigate the origin of these trends, which suggests interesting connections between jamming and thermal glassy phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coherent quantum dynamics in donor-bridge-acceptor systems: Beyond the hopping and super-exchange mechanisms

The population transfer dynamics of model donor-bridge-acceptor systems is studied by comparing a recently developed polaron-transformed quantum master equation (PQME) with the well-known Redfield and Forster theories of quantum transport. We show that the PQME approach reduces to these two theories in their respective limits of validity and naturally interpolates between them as a function of the system-bath coupling strength. By exploring the parameter space of our model problem, we identify novel regimes of transport dynamics in bridged systems like those encountered in biological and organic energy transfer problems. Furthermore, we demonstrate that three-level systems like the ones studied herein represent ideal minimal models for the identification of quantum coherent transport as embodied in super-exchange phenomena that cannot be captured by Forster-like hopping approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamic Length Scales in Glass-Forming Liquids: A Inhomogeneous Molecular Dynamics Simulation Approach

In this work we numerically investigate a new method for the characterization of growing length scales associated with spatially heterogeneous dynamics of glass-forming liquids. This approach, motivated by the formulation of the inhomogeneous mode-coupling theory (IMCT) [Biroli G. et al. Phys. Rev. Lett. 2006 97, 195701], utilizes inhomogeneous molecular dynamics simulations in which the system is perturbed by a spatially modulated external potential. We show that the response of the two-point correlation function to the external field allows one to probe dynamic correlations. We examine the critical properties shown by this function, in particular the associated dynamic correlation length, that is found to be comparable to the one extracted from standardly-employed four-point correlation functions. Our numerical results are in qualitative agreement with IMCT predictions but suggest that one has to take into account fluctuations not included in this mean-field approach in order to reach quantitative agreement. Advantages of our approach over the more conventional one based on four-point correlation functions are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Grains and grain boundaries in highly crystalline monolayer molybdenum disulfide

Recent progress in large-area synthesis of monolayer molybdenum disulfide, a new two-dimensional direct-bandgap semiconductor, is paving the way for applications in atomically thin electronics. Little is known, however, about the microstructure of this material. Here we have refined chemical vapor deposition synthesis to grow highly crystalline islands of monolayer molybdenum disulfide up to 120 um in size with optical and electrical properties comparable or superior to exfoliated samples. Using transmission electron microscopy, we correlate lattice orientation, edge morphology, and crystallinity with island shape to demonstrate that triangular islands are single crystals. The crystals merge to form faceted tilt and mirror boundaries that are stitched together by lines of 8- and 4- membered rings. Density functional theory reveals localized mid-gap states arising from these 8-4 defects. We find that mirror boundaries cause strong photoluminescence quenching while tilt boundaries cause strong enhancement. In contrast, the boundaries only slightly increase the measured in-plane electrical conductivity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Microscopic theory to quantify the competing kinetic processes in photoexcited surface-coupled quantum dots

We present a self-contained theoretical and computational framework for dynamics following photoexcitation in quantum dots near planar interfaces. A microscopic Hamiltonian parameterized by first principles calculations is merged with a reduced density matrix formalism that allows for the prediction of time-dependent charge and energy transfer processes between the quantum dot and the electrode. While treating charge and energy transfer processes on an equal footing, the non-perturbative effects of sudden charge transitions on the Fermi sea of the electrode are included. We illustrate the formalism with calculations of an InAs quantum dot coupled to the Shockley state on an Au[111] surface, and use it to concretely discuss the wide range of kinetics possible in these systems and their implications for photovoltaic systems and tunnel junction devices. We discuss the utility of this framework for the analysis of recent experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Theory of Neutral and Charged Excitons in Monolayer Transition Metal Dichalcogenides

We present a microscopic theory of neutral excitons and charged excitons (trions) in monolayers of transition metal dichalcogenides, including molybdenum disulfide. Our theory is based on an effective mass model of excitons and trions, parametrized by ab initio calculations and incorporating a proper treatment of screening in two dimensions. The calculated exciton binding energies are in good agreement with high-level many-body computations based on the Bethe-Salpeter equation. Furthermore, our calculations for the more complex trion species compare very favorably with recent experimental measurements, and provide atomistic insight into the microscopic features which determine the trion binding energy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effects of the interfacial polarization on tunneling in surface coupled quantum dots

Polarization effects are included exactly in a model for a quantum dot in close proximity to a planar interface. Efficient incorporation of this potential into the Schrödinger equation is utilized to map out the influence of the image potential effects on carrier tunneling in such heterostructures. In particular, the interplay between carrier mass and the dielectric constants of a quantum dot, its surrounding matrix, and the electrode is studied. We find that the polarizability of the planar electrode structure can significantly increase the tunneling rates for heavier carriers, potentially resulting in a qualitative change in the dependence of tunneling rate on mass. Our method for treating polarization can be generalized to the screening of two particle interactions, and can thus be applied to calculations such as exciton dissociation and the Coulomb blockade. In contrast to tunneling via intermediate surface localized states of the quantum dot, our work identifies the parameter space over which volume states undergo significant modification in their tunneling characteristics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite-Temperature Auxiliary-Field Quantum Monte Carlo for Bose-Fermi Mixtures

We present a quantum Monte Carlo (QMC) technique for calculating the exact finite-temperature properties of Bose-Fermi mixtures. The Bose-Fermi Auxiliary-Field Quantum Monte Carlo (BF-AFQMC) algorithm combines two methods, a finite-temperature AFQMC algorithm for bosons and a variant of the standard AFQMC algorithm for fermions, into one algorithm for mixtures. We demonstrate the accuracy of our method by comparing its results for the Bose-Hubbard and Bose-Fermi-Hubbard models against those produced using exact diagonalization for small systems. Comparisons are also made with mean-field theory and the worm algorithm for larger systems. As is the case with most fermion Hamiltonians, a sign or phase problem is present in BF-AFQMC. We discuss the nature of these problems in this framework and describe how they can be controlled with well-studied approximations to expand BF-AFQMC's reach. The new algorithm can serve as an essential tool for answering many unresolved questions about many-body physics in mixed Bose-Fermi systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Growing point-to-set length scale correlates with growing relaxation times in model supercooled liquids

It has been demonstrated recently that supercooled liquids sharing simple structural features (e.g. pair distribution functions) may exhibit strikingly distinct dynamical behavior. Here we show that a more subtle structural feature correlates with relaxation times in three simulated systems that have nearly identical radial distribution functions but starkly different dynamical behavior. In particular, for the first time we determine the thermodynamic "point-to-set" length scale in several canonical model systems and demonstrate the quantitative connection between this length scale and the growth of relaxation times. Our results provide clues necessary for distinguishing competing theories of the glass transition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Microscopic theory of singlet exciton fission. I. General formulation

Singlet fission, a spin-allowed energy transfer process generating two triplet excitons from one singlet exciton, has the potential to dramatically increase the efficiency of organic solar cells. However, the dynamical mechanism of this phenomenon is not fully understood and a complete, microscopic theory of singlet fission is lacking. In this work, we assemble the components of a comprehensive microscopic theory of singlet fission that connects excited state quantum chemistry calculations with finite-temperature quantum relaxation theory. We elaborate on the distinction between localized diabatic and delocalized adiabatic bases for the interpretation of singlet fission experiments in both the time and frequency domains. We discuss various approximations to the exact density matrix dynamics and propose Redfield theory as an ideal compromise between speed and accuracy for the detailed investigation of singlet fission in dimers, clusters, and crystals. Investigations of small model systems based on parameters typical of singlet fission demonstrate the numerical accuracy and practical utility of this approach.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Microscopic theory of singlet exciton fission. II. Application to pentacene dimers and the role of superexchange

In the preceding paper, we assembled the theoretical components necessary for a unified framework of singlet fission, a type of multiexciton generation producing two triplet excitons from one singlet exciton. In this paper, we apply our methodology to molecular dimers of pentacene, a widely studied material that exhibits singlet fission. We address a longstanding theoretical issue, namely whether singlet fission proceeds via two sequential electron transfer steps mediated by a charge-transfer state or via a direct two-electron transfer process. We find evidence for a superexchange mediated mechanism, whereby the fission process proceeds through virtual charge-transfer states which may be very high in energy. In particular, this mechanism predicts efficient singlet fission on the sub-picosecond timescale, in reasonable agreement with experiment. We investigate the role played by molecular vibrations in mediating relaxation and decoherence, finding that different physically reasonable forms for the bath relaxation function give similar results. We also examine the competing direct coupling mechanism and find it to yield fission rates slower in comparison with the superexchange mechanism for the dimer. We discuss implications for crystalline pentacene, including the limitations of the dimer model.

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
preprint2011arXiv

Dynamical Mean-Field Theory for Quantum Chemistry

The dynamical mean-field concept of approximating an unsolvable many-body problem in terms of the solution of an auxiliary quantum impurity problem, introduced to study bulk materials with a continuous energy spectrum, is here extended to molecules, i.e., finite systems with a discrete energy spectrum. The application to small clusters of hydrogen atoms yields ground state energies which are competitive with leading quantum chemical approaches at intermediate and large interatomic distances as well as good approximations to the excitation spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerically Exact Long Time Behavior of Nonequilibrium Quantum Impurity Models

A Monte Carlo sampling of diagrammatic corrections to the non-crossing approximation is shown to provide numerically exact estimates of the long-time dynamics and steady state properties of nonequilibrium quantum impurity models. This `bold' expansion converges uniformly in time and significantly ameliorates the sign problem that has heretofore limited the power of real-time Monte Carlo approaches to strongly interacting real-time quantum problems. The new approach enables the study of previously intractable problems ranging from generic long time nonequilibrium transport characteristics in systems with large onsite repulsion to the direct description of spectral functions on the real frequency axis in Dynamical Mean Field Theory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum quench spectroscopy of a Luttinger liquid: Ultrarelativistic density wave dynamics due to fractionalization in an XXZ chain

We compute the dynamics of localized excitations produced by a quantum quench in the spin 1/2 XXZ chain. Using numerics combining the density matrix renormalization group and exact time evolution, as well as analytical arguments, we show that fractionalization due to interactions in the pre-quench state gives rise to "ultrarelativistic" density waves that travel at the maximum band velocity. The system is initially prepared in the ground state of the chain within the gapless XY phase, which admits a Luttinger liquid (LL) description at low energies and long wavelengths. The Hamiltonian is then suddenly quenched to a band insulator, after which the chain evolves unitarily. Through the gapped dispersion of the insulator spectrum, the post-quench dynamics serve as a "velocity microscope," revealing initial state particle correlations via space time density propagation. We show that the ultrarelativistic wave production is tied to the particular way in which fractionalization evades Pauli-blocking in the zero-temperature initial LL state.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Reduced density matrix hybrid approach: An efficient and accurate method for adiabatic and non-adiabatic quantum dynamics

We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into "core" and "reservoir" modes, with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Reduced density matrix hybrid approach: Application to electronic energy transfer

Electronic energy transfer in the condensed phase, such as that occurring in photosynthetic complexes, frequently occurs in regimes where the energy scales of the system and environment are similar. This situation provides a challenge to theoretical investigation since most approaches are accurate only when a certain energetic parameter is small compared to others in the problem. Here we show that in these difficult regimes, the Ehrenfest approach provides a good starting point for a dynamical description of the energy transfer process due to its ability to accurately treat coupling to slow environmental modes. To further improve on the accuracy of the Ehrenfest approach, we use our reduced density matrix hybrid framework to treat the faster environmental modes quantum mechanically, at the level of a perturbative master equation. This combined approach is shown to provide an efficient and quantitative description of electronic energy transfer in a model dimer and the Fenna-Matthews-Olson complex and is used to investigate the effect of environmental preparation on the resulting dynamics.

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
preprint2011arXiv

Visualizing Individual Nitrogen Dopants in Monolayer Graphene

In monolayer graphene, substitutional doping during growth can be used to alter its electronic properties. We used scanning tunneling microscopy (STM), Raman spectroscopy, x-ray spectroscopy, and first principles calculations to characterize individual nitrogen dopants in monolayer graphene grown on a copper substrate. Individual nitrogen atoms were incorporated as graphitic dopants, and a fraction of the extra electron on each nitrogen atom was delocalized into the graphene lattice. The electronic structure of nitrogen-doped graphene was strongly modified only within a few lattice spacings of the site of the nitrogen dopant. These findings show that chemical doping is a promising route to achieving high-quality graphene films with a large carrier concentration.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Bold Line Diagrammatic Monte Carlo Method: General formulation and application to expansion around the Non-Crossing Approximation

We present a general framework for performing "bold-line" diagrammatic Monte Carlo calculations based on using an analytical partial resummation as a starting point for stochastically summing all diagrams. As an example we solve the one orbital Anderson model using the non-crossing approximation as a starting point. We assess the accuracy of the method by using it to solve the equations of single-site dynamical mean field theory. We establish the validity of the starting approximations and show that the bold method provides a very accurate treatment of the Mott-insulating phase.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Conductivity of disordered quantum lattice models at infinite temperature: Many-body localization

We reinvestigate the behavior of the conductivity of several disordered quantum lattice models at infinite temperature using exact diagonalization. Contrary to the conclusion drawn in a recent investigation of similar quantities in identical systems, we find evidence of a localized regime for strong random fields. We estimate the location of the critical field for the many-body localization transition for the random-field XXZ spin chain, and compare our findings with recent investigations in related systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effect of a Coulombic dot-lead coupling on the dynamics of a quantum dot

The effect of a Coulombic coupling on the dynamics of a quantum dot hybridized to leads is determined. The calculation treats the interaction between charge fluctuations on the dot and the dynamically generated image charge in the leads. A formally exact solution is presented for a dot coupled to a Luttinger liquid and an approximate solution, equivalent to treating the lead dynamics within a random phase approximation, is given for a dot coupled to a two- or three-dimensional metallic lead. The leading divergences arising from the long-ranged Coulomb interaction are found to cancel, so that in the two- and three-dimensional cases the quantum-dot dynamics is equivalent to that obtained by neglecting both the dot-lead Coulomb coupling and the Coulomb renormalization of the lead electrons, while in the one-dimensional case the dot-lead mixing is enhanced relative to the non-interacting case. Explicit results are given for the short-time dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerically exact path integral simulation of nonequilibrium quantum transport and dissipation

We develop an iterative, numerically exact approach for the treatment of nonequilibrium quantum transport and dissipation problems that avoids the real-time sign problem associated with standard Monte Carlo techniques. The method requires a well-defined decorrelation time of the non-local influence functional for proper convergence to the exact limit. Since finite decorrelation times may arise either from temperature or from a voltage drop at zero temperature, the approach is well suited for the description of the real-time dynamics of single-molecule devices and quantum dots driven to a steady-state via interaction with two or more electron leads. We numerically investigate two non-trivial models: the evolution of the nonequilibrium population of a two-level system coupled to two electronic reservoirs, and quantum transport in the nonequilibrium Anderson model. For the latter case, two distinct formulations are described. Results are compared to those obtained by other techniques.

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
preprint2010arXiv

Transport through a quantum dot with excitonic dot-lead coupling

We study the effect of a dot-lead interaction on transport through a quantum dot hybridized to two semi-infinite Luttinger-liquid leads. A bosonization approach is applied to treat the interaction between charge fluctuations on the dot and the dynamically generated image charge in the leads. The nonequilibrium distribution function of the dot and the tunneling current are computed within a master-equation approach. The presence of the excitonic dot-lead coupling is found to enhance transport in the vicinity of the Coulomb-blockade threshold. This behavior is in contrast to the usual power-law suppression of electronic tunneling which is found if this interaction is ignored.

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
preprint2009arXiv

Dynamical Heterogeneity in Lattice Glass Models

In this paper we consider in detail the properties of dynamical heterogeneity in lattice glass models (LGMs). LGMs are lattice models whose dynamical rules are based on thermodynamic, as opposed to purely kinetic, considerations. We devise a LGM that is not prone to crystallization and displays properties of a fragile glass-forming liquid. Particle motion in this model tends to be locally anisotropic on intermediate time scales even though the rules governing the model are isotropic. The model demonstrates violations of the Stokes-Einstein relation and the growth of various length scales associated with dynamical heterogeneity. We discuss future avenues of research comparing the predictions of lattice glass models and kinetically constrained models to atomistic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Localised soft modes and the supercooled liquid's irreversible passage through its configuration space

Using computer simulations, we show that the localized low frequency normal modes of a configuration in a supercooled liquid are strongly correlated with the irreversible structural reorganization of the particles within that configuration. Establishing this correlation constitutes the identification of the aspect of a configuration that determines the heterogeneity of the subsequent motion. We demonstrate that the spatial distribution of the summation over the soft local modes can persist in spite of particle reorganization that produces significant changes in individual modes. Along with spatial localization, the persistent influence of soft modes in particle relaxation results in anisotropy in the displacements of mobile particles over the timescale referred to as beta relaxation.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Non-Equilibrium Quantum Dissipation

Dissipative processes in non-equilibrium many-body systems are fundamentally different than their equilibrium counterparts. Such processes are of great importance for the understanding of relaxation in single molecule devices. As a detailed case study, we investigate here a generic spin-fermion model, where a two-level system couples to two metallic leads with different chemical potentials. We present results for the spin relaxation rate in the nonadiabatic limit for an arbitrary coupling to the leads, using both analytical and exact numerical methods. The non-equilibrium dynamics is reflected by an exponential relaxation at long times and via complex phase shifts, leading in some cases to an "anti-orthogonality" effect. In the limit of strong system-lead coupling at zero temperature we demonstrate the onset of a Marcus-like Gaussian decay with {\it voltage difference} activation. This is analogous to the equilibrium spin-boson model, where at strong coupling and high temperatures the spin excitation rate manifests temperature activated Gaussian behavior. We find that there is no simple linear relationship between the role of the temperature in the bosonic system and a voltage drop in a non-equilibrium electronic case. The two models also differ by the orthogonality-catastrophe factor existing in a fermionic system, which modifies the resulting lineshapes. Implications for current characteristics are discussed. We demonstrate the violation of pair-wise Coulomb gas behavior for strong coupling to the leads. The results presented in this paper form the basis of an exact, non-perturbative description of steady-state quantum dissipative systems.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Zeno and anti-Zeno dynamics in spin-bath models

We investigate the quantum Zeno and anti-Zeno effects in spin bath models: the spin-boson model and a spin-fermion model. We show that the Zeno-anti-Zeno transition is critically controlled by the system-bath coupling parameter, the same parameter that determines spin decoherence rate. We also discuss the crossover in a biased system, at high temperatures, and for a nonequilibrium spin-fermion system, manifesting the counteracting roles of electrical bias, temperature, and magnetic field on the spin decoherence rate.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Mode-Coupling Theory (MCT) Lecture Notes

In this set of lecture notes we review the mode-coupling theory of the glass transition from several perspectives. First, we derive mode-coupling equations for the description of density fluctuations from microscopic considerations with the use the Mori-Zwanzig projection operator technique. We also derive schematic mode-coupling equations of a similar form from a field-theoretic perspective. We review the successes and failures of mode-coupling theory, and discuss recent advances in the applications of the theory.

0 citations0 reviewsNo code linkedOpen access