BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

David R. Reichman

David R. Reichman contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phcond-mat.str-elcond-mat.dis-nncond-mat.mtrl-scicond-mat.stat-mechphysics.comp-phquant-phcond-mat.mes-hallcond-mat.softMachine Learning

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 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