Researcher profile

Timothy C. Berkelbach

Timothy C. Berkelbach contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
6topics
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

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

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.