Source author record

Timothy C. Berkelbach

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

24 published item(s)

preprint2022arXiv

Anharmonic Lattice Dynamics from Vibrational Dynamical Mean-Field Theory

We present a vibrational dynamical mean-field theory (VDMFT) of the dynamics of atoms in solids with anharmonic interactions. Like other flavors of DMFT, VDMFT maps the dynamics of a periodic anharmonic lattice of atoms onto those of a self-consistently defined impurity problem with local anharmonicity and coupling to a bath of harmonic oscillators. VDMFT is exact in the harmonic and molecular limits, nonperturbative, systematically improvable through its cluster extensions, usable with classical or quantum impurity solvers (depending on the importance of nuclear quantum effects), and can be combined with existing low-level diagrammatic theories of anharmonicity. When tested on models of anharmonic optical and acoustic phonons, we find that classical VDMFT gives good agreement with classical molecular dynamics, including the temperature dependence of phonon frequencies and lifetimes. Using a quantum impurity solver, signatures of nuclear quantum effects are observed at low temperatures. We test the description of nonlocal anharmonicity via cellular VDMFT and the combination with self-consistent phonon (SCPH) theory, yielding the powerful SCPH+VDMFT approach.

preprint2022arXiv

Correlation-consistent Gaussian basis sets for solids made simple

The rapidly growing interest in simulating condensed-phase materials using quantum chemistry methods calls for a library of high-quality Gaussian basis sets suitable for periodic calculations. Unfortunately, most standard Gaussian basis sets commonly used in molecular simulation show significant linear dependencies when used in close-packed solids, leading to severe numerical issues that hamper the convergence to the complete basis set (CBS) limit, especially in correlated calculations. In this work, we revisit Dunning's strategy for construction of correlation-consistent basis sets and examine the relationship between accuracy and numerical stability in periodic settings. Specifically, we find that limiting the number of primitive functions avoids the appearance of problematic small exponents while still providing smooth convergence to the CBS limit. As an example, we generate double-, triple-, and quadruple-zeta correlation-consistent Gaussian basis sets for periodic calculations with Goedecker-Teter-Hutter (GTH) pseudopotentials. Our basis sets cover the main-group elements from the first three rows of the periodic table. Especially for atoms on the left side of the periodic table, our basis sets are less diffuse than those used in molecular calculations. We verify the fast and reliable convergence to the CBS limit in both Hartree-Fock and post-Hartree-Fock (MP2) calculations, using a diverse test set of $19$ semiconductors and insulators.

preprint2022arXiv

Dark-exciton driven energy funneling into dielectric inhomogeneities in two-dimensional semiconductors

The optoelectronic and transport properties of two-dimensional transition metal dichalcogenide semiconductors (2D TMDs) are highly susceptible to external perturbation, enabling precise tailoring of material function through post-synthetic modifications. Here we show that nanoscale inhomogeneities known as nanobubbles can be used for both strain and, less invasively, dielectric tuning of exciton transport in bilayer tungsten disulfide (WSe2). We use ultrasensitive spatiotemporally resolved optical scattering microscopy to directly image exciton transport, revealing that dielectric nanobubbles are surprisingly efficient at funneling and trapping excitons at room temperature, even though the energies of the bright excitons are negligibly affected. Our observations suggest that exciton funneling in dielectric inhomogeneities is driven by momentum-indirect (dark) excitons whose energies are more sensitive to dielectric perturbations than bright excitons. These results reveal a new pathway to control exciton transport in 2D semiconductors with exceptional spatial and energetic precision using dielectric engineering of dark state energetic landscapes.

preprint2022arXiv

Ground-state properties of metallic solids from ab initio coupled-cluster theory

Metallic solids are a challenging target for wavefunction-based electronic structure theories and have not been studied in great detail by such methods. Here, we use coupled-cluster theory with single and double excitations (CCSD) to study the structure of solid lithium and aluminum using optimized Gaussian basis sets. We calculate the equilibrium lattice constant, bulk modulus, and cohesive energy and compare them to experimental values, finding accuracy comparable to common density functionals. Because the quantum chemical "gold standard" CCSD(T) (CCSD with perturbative triple excitations) is inapplicable to metals in the thermodynamic limit, we test two approximate improvements to CCSD, which are found to improve the predicted cohesive energies.

preprint2022arXiv

The normal state of attractive Fermi gases from coupled-cluster theory

We introduce coupled-cluster (CC) theory for the numerical study of the normal state of two-component, dilute Fermi gases with attractive, short-range interactions at zero temperature. We focus on CC theory with double excitations (CCD) and discuss its close relationship with -- and improvement upon -- the t-matrix approximation, i.e., the resummation of ladder diagrams via a random-phase approximation. We further discuss its relationship with Chevy's variational wavefunction ansatz for the Fermi polaron and argue that CCD is its natural extension to nonzero minority species concentrations. Studying normal state energetics for a range of interaction strengths below and above unitarity, we find that CCD yields good agreement with fixed-node diffusion Monte Carlo. We find that CCD does not converge for small polarizations and large interaction strengths, which we speculatively attribute to the nascent instability to a superfluid state.

preprint2021arXiv

Full-frequency dynamical Bethe-Salpeter equation without frequency and a study of double excitations

The Bethe-Salpeter equation (BSE) that results from the GW approximation to the self-energy is a frequency-dependent (nonlinear) eigenvalue problem due to the dynamically screened Coulomb interaction between electrons and holes. The computational time required for a numerically exact treatment of this frequency dependence is $O(N^6)$, where $N$ is the system size. To avoid the common static screening approximation, we show that the full-frequency dynamical BSE can be exactly reformulated as a frequency-independent eigenvalue problem in an expanded space of single and double excitations. When combined with an iterative eigensolver and the density fitting approximation to the electron repulsion integrals, this reformulation yields a dynamical BSE algorithm whose computational time is $O(N^5)$, which we verify numerically. Furthermore, the reformulation provides direct access to excited states with dominant double excitation character, which are completely absent in the spectrum of the statically screened BSE. We study the $2^1A_\mathrm{g}$ state of butadiene, hexatriene, and octatetraene and find that GW/BSE overestimates the excitation energy by about 1.5-2 eV and significantly underestimates the double excitation character.

preprint2020arXiv

Excitons in solids from periodic equation-of-motion coupled-cluster theory

We present an ab initio study of electronically excited states of three-dimensional solids using Gaussian-based periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). The explicit use of translational symmetry, as implemented via Brillouin zone sampling and momentum conservation, is responsible for a large reduction in cost. Our largest system studied, which samples the Brillouin zone using 64 k-points (a 4x4x4 mesh) corresponds to a canonical EOM-CCSD calculation of 768 electrons in 640 orbitals. We study eight simple semiconductors and insulators, with direct singlet excitation energies in the range of 3 to 15 eV. Our predicted excitation energies exhibit a mean absolute error of 0.27 eV when compared to experiment. We furthermore calculate the energy of excitons with nonzero momentum and compare the exciton dispersion of LiF with experimental data from inelastic X-ray scattering. By calculating excitation energies under strain, we extract hydrostatic deformation potentials in order to quantify the strength of interactions between excitons and acoustic phonons. Our results indicate that coupled-cluster theory is a promising method for the accurate study of a variety of exciton phenomena in solids.

preprint2020arXiv

Full-Frequency GW without Frequency

Efficient computer implementations of the GW approximation must approximate a numerically challenging frequency integral; the integral can be performed analytically, but doing so leads to an expensive implementation whose computational cost scales as $O(N^6)$ where $N$ is the size of the system. Here we introduce a new formulation of the full-frequency GW approximation by exactly recasting it as an eigenvalue problem in an expanded space. This new formulation (1) avoids the use of time or frequency grids, (2) naturally precludes the common "diagonal" approximation, (3) enables common iterative eigensolvers that reduce the canonical scaling to $O(N^5)$, and (4) enables a density-fitted implementation that reduces the scaling to $O(N^4)$. We numerically verify these scaling behaviors and test a variety of approximations that are motivated by this new formulation. In this new formulation, the relation of the GW approximation to configuration interaction, coupled-cluster theory, and the algebraic diagrammatic construction is made especially apparent, providing a new direction for improvements to the GW approximation.

preprint2019arXiv

A unification of the Holstein polaron and dynamic disorder pictures of charge transport in organic semiconductors

We present a unified and nonperturbative method for calculating spectral and transport properties of Hamiltonians with simultaneous Holstein (diagonal) and Peierls (off-diagonal) electron-phonon coupling. Our approach is motivated by the separation of energy scales in semiconducting organic molecular cystals, in which electrons couple to high-frequency intramolecular Holstein modes and to low-frequency intermolecular Peierls modes. We treat Peierls modes as quasi-classical dynamic disorder, while Holstein modes are included with a Lang-Firsov polaron transformation and no narrow-band approximation. Our method reduces to the popular polaron picture due to Holstein coupling and the dynamic disorder picture due to Peierls coupling. We derive an expression for efficient numerical evaluation of the frequency-resolved optical conductivity based on the Kubo formula and obtain the DC mobility from its zero-frequency component. We also use our method to calculate the electron-addition Green's function corresponding to the inverse photoemission spectrum. For realistic parameters, temperature-dependent DC mobility is largely determined by the Peierls-induced dynamic disorder with minor quantitative corrections due to polaronic band-narrowing, and an activated regime is not observed at relevant temperatures. In contrast, for frequency-resolved observables, a quantum mechanical treatment of the Holstein coupling is qualitatively important for capturing the phonon replica satellite structure.

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.

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.

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.

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.

preprint2015arXiv

Spectral Functions of the Uniform Electron Gas via Coupled-Cluster Theory and Comparison to the $GW$ and Related Approximations

We use, for the first time, ab initio coupled-cluster theory to compute the spectral function of the uniform electron gas at a Wigner-Seitz radius of $r_\mathrm{s}=4$. The coupled-cluster approximations we employ go significantly beyond the diagrammatic content of state-of-the-art $GW$ theory. We compare our calculations extensively to $GW$ and $GW$-plus-cumulant theory, illustrating the strengths and weaknesses of these methods in capturing the quasiparticle and satellite features of the electron gas. Our accurate calculations further allow us to address the long-standing debate over the occupied bandwidth of metallic sodium. Our findings indicate that the future application of coupled-cluster theory to condensed phase material spectra is highly promising.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.