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

James J. Shepherd

James J. Shepherd contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.comp-phphysics.chem-phcond-mat.mtrl-scicond-mat.quant-gas

Trust snapshot

Quick read

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

6 published item(s)

preprint2022arXiv

Machine learning for a finite size correction in periodic coupled cluster theory calculations

We introduce a straightforward Gaussian process regression (GPR) model for the transition structure factor of metal periodic coupled cluster singles and doubles (CCSD) calculations. This is inspired by the method introduced by Liao and Grüneis for interpolating over the transition structure factor to obtain a finite size correction for CCSD [J. Chem. Phys. 145, 141102 (2016)], and by our own prior work using the transition structure factor to efficiently converge CCSD for metals to the thermodynamic limit [Nat. Comput. Sci. 1, 801 (2021)]. In our CCSD-FS-GPR method to correct for finite size errors, we fit the structure factor to a 1D function in the momentum transfer, $G$. We then integrate over this function by projecting it onto a k-point mesh to obtain comparisons with extrapolated results. Results are shown for lithium, sodium, and the uniform electron gas.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Piecewise Interaction Picture Density Matrix Quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact $N$-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC method (IP-DMQMC) which overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Effective Hamiltonians for the study of real metals using quantum chemical theories

Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are at the heart of computational materials science. In that respect the coupled cluster hierarchy of methods plays a central role in molecular quantum chemistry because of its systematic improvability and computational efficiency. In this hierarchy, coupled cluster singles and doubles (CCSD) is one of the most important steps in moving towards chemical accuracy and, in recent years, its scope has successfully been expanded to the study of insulating surfaces and solids. Here, we show that CCSD theory can also be applied to real metals. In so doing, we overcome the limitation of needing extremely large supercells to capture long range electronic correlation effects. An effective Hamiltonian can be found using the transition structure factor--a map of electronic excitations from the Hartree--Fock wavefunction--which has fewer finite size effects than conventional periodic boundary conditions. This not only paves the way of applying coupled cluster methods to real metals but also reduces the computational cost by two orders of magnitude compared to previous methods. Our applications to phases of lithium and silicon show a resounding success in reaching the thermodynamic limit, taking the first step towards a truly universal quantum chemical treatment of solids.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Power laws used to extrapolate the coupled cluster correlation energy to the thermodynamic limit

Recent calculations using coupled cluster on solids have raised discussion of using a $N^{-1/3}$ power law to fit the correlation energy when extrapolating to the thermodynamic limit, an approach which differs from the more commonly used $N^{-1}$ power law which is (for example) often used by quantum Monte Carlo methods. In this paper, we present one way to reconcile these viewpoints. Coupled cluster doubles calculations were performed on uniform electron gases reaching system sizes of $922$ electrons for an extremely wide range of densities ($0.1<r_s<100.0$) to study how the correlation energy approaches the thermodynamic limit. The data were corrected for basis set incompleteness error and use a selected twist angle approach to mitigate finite size error from shell filling effects. Analyzing these data, we initially find that a power law of $N^{-1/3}$ appears to fit the data better than a $N^{-1}$ power law in the large system size limit. However, we provide an analysis of the transition structure factor showing that $N^{-1}$ still applies to large system sizes and that the apparent $N^{-1/3}$ power law occurs only at low $N$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accelerating convergence to the thermodynamic limit with twist angle selection applied to methods beyond many-body perturbation theory

We recently developed a scheme to use low-cost calculations to find a single twist angle where the couple cluster doubles energy of a single calculation matches the twist-averaged coupled cluster doubles energy in a finite unit cell. We used initiator full configuration interaction quantum Monte Carlo ($i$-FCIQMC) as an example of an exact method beyond coupled cluster doubles theory to show that this selected twist angle approach had comparable accuracy in methods beyond coupled cluster. Further, at least for small system sizes, we show that the same twist angle can also be found by comparing the energy directly (at the level of second-order Moller-Plesset theory) suggesting a route toward twist angle selection which requires minimal modification to existing codes which can perform twist averaging.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

NECI: N-Electron Configuration Interaction with emphasis on state-of-the-art stochastic methods

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24000 CPU cores. In this paper, we describe the core functionalities of NECI and recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green&#39;s functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs and it is licensed under GPL-3.0.

0 citations0 reviewsNo code linkedOpen access