pith. sign in

arxiv: 2604.12142 · v1 · submitted 2026-04-13 · 🪐 quant-ph · cond-mat.mtrl-sci· cond-mat.str-el

Fault-tolerant simulation of the electronic structure using Projector Augmented-Waves and Bloch orbitals

Rishabh Bhardwaj , Alexander Reed Mu\~noz , Travis E. Jones , John Golden This is my paper

Pith reviewed 2026-05-10 14:54 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.mtrl-scicond-mat.str-el
keywords Bloch orbitalsprojector augmented wavesfault-tolerant quantum simulationelectronic structureperiodic solidsqubitizationresource estimationblock encoding
0
0 comments X

The pith

The Bloch-UPAW framework reduces Toffoli counts by roughly an order of magnitude for fault-tolerant quantum simulation of periodic solids such as bulk diamond.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops the Bloch-UPAW method to simulate electronic structure in bulk materials on fault-tolerant quantum computers by combining Bloch orbitals that encode periodic k-space structure with unitary projector-augmented-wave corrections that capture sharp nuclear variations through strictly local on-site terms. The resulting Hamiltonian stays expressed in the Bloch basis so Brillouin-zone sampling remains explicitly controllable, and the construction works for both plane-wave and localized bases while handling supercells efficiently. A linear-combination-of-unitaries decomposition yields a block-encoding circuit for qubitization in which the UPAW augmentation adds only one ancilla qubit and no leading-order Toffoli gates. Resource estimates for diamond indicate an order-of-magnitude drop in Toffoli count relative to earlier periodic-solid simulations, with asymptotic scaling O(N_k^3) under k-mesh refinement and O(N_a^{3.5}) under supercell growth. A reader would care because strongly correlated materials are central to technology yet remain difficult to converge classically.

Core claim

The Bloch-UPAW framework combines Bloch-orbital k-space structure with unitary projector-augmented-wave augmentation. The UPAW Hamiltonian is expressed directly in the Bloch basis, retains explicit control of Brillouin-zone sampling, and incorporates near-nuclear physics through strictly local on-site corrections. The construction is independent of the underlying one-particle representation and applies to both plane-wave and localized bases. We derive a linear-combination-of-unitaries decomposition and a block-encoding circuit suitable for qubitization; UPAW augmentation adds one ancilla qubit and no Toffoli gates at leading order relative to a Bloch-only block encoding. Asymptotically theTo

What carries the argument

The UPAW augmentation of the Bloch-orbital Hamiltonian, which supplies local nuclear corrections while preserving the linear-combination-of-unitaries block encoding with only one added ancilla and no leading Toffoli overhead.

If this is right

  • Toffoli cost scales as O(N_k^3) when the k-mesh is refined while keeping the supercell fixed.
  • Toffoli cost scales as O(N_a^{3.5}) when the supercell is enlarged while keeping the k-mesh fixed.
  • The same block-encoding construction applies unchanged to both plane-wave and localized orbital bases.
  • Supercells required for symmetry-breaking phenomena can be treated with the same asymptotic cost scaling as ordinary periodic cells.
  • Bulk properties can be converged by whichever route (k-mesh or supercell) is cheaper for a given material.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method's independence from the one-particle basis suggests it could be paired with existing localized-orbital quantum algorithms to further reduce resources for materials with localized correlations.
  • Because on-site corrections are strictly local, the framework may allow smaller supercells than pure Bloch methods when strong nuclear effects dominate, shortening the path to bulk convergence.
  • The reported scaling opens the possibility that fault-tolerant simulations of larger unit cells become practical before qubit counts reach the thresholds assumed in earlier periodic-solid estimates.
  • If the same augmentation logic extends to other classical augmentation schemes, the overhead reduction seen in diamond could appear across a wider range of solids.

Load-bearing premise

The linear-combination-of-unitaries decomposition of the augmented Hamiltonian stays efficient and the UPAW step introduces neither extra ancillas nor Toffoli gates that grow with system size.

What would settle it

An explicit circuit or resource count for the full UPAW block encoding that requires more than one additional ancilla qubit or introduces Toffoli gates scaling faster than the Bloch-only baseline would falsify the leading-order cost claim.

Figures

Figures reproduced from arXiv: 2604.12142 by Alexander Reed Mu\~noz, John Golden, Rishabh Bhardwaj, Travis E. Jones.

Figure 1
Figure 1. Figure 1: FIG. 1: Regimes covered by three fault-tolerant ap [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p020_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p024_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p025_6.png] view at source ↗
read the original abstract

Strongly correlated materials are a natural target for fault-tolerant quantum computers, but they require tools beyond those developed for molecules. Electronic wavefunctions vary rapidly near nuclei yet remain delocalized across many unit cells, and bulk properties must be converged systematically with respect to finite-size errors. To resolve such issues, we present the Bloch--UPAW framework that combines Bloch-orbital $k$-space structure with unitary projector-augmented-wave (UPAW) augmentation. The UPAW Hamiltonian, expressed directly in the Bloch basis, retains explicit control of Brillouin-zone sampling, and incorporates near-nuclear physics through strictly local on-site corrections. The construction is independent of the underlying one-particle representation, so it applies to both plane-wave and localized bases, and it handles supercells for symmetry-breaking phenomena more efficiently. We derive a linear-combination-of-unitaries decomposition and a block-encoding circuit suitable for qubitization; UPAW augmentation adds one ancilla qubit and no Toffoli gates at leading order relative to a Bloch-only block encoding. Asymptotically, the Toffoli cost scales as $\mathcal{O}(N_k^3)$ when refining the $k$-mesh and as $\mathcal{O}(N_a^{3.5})$ when enlarging the supercell, enabling convergence to be steered by the most favorable route for a given material. Resource estimates for bulk diamond show approximately an order-of-magnitude reduction in Toffoli count relative to prior work on periodic solids.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript introduces the Bloch-UPAW framework for fault-tolerant quantum simulation of electronic structure in periodic solids. It combines Bloch-orbital k-space sampling with unitary projector-augmented-wave (UPAW) augmentation to incorporate near-nuclear physics via strictly local on-site corrections. The authors derive a linear-combination-of-unitaries (LCU) decomposition and block-encoding circuit for qubitization, asserting that UPAW adds one ancilla qubit and no Toffoli gates at leading order relative to a Bloch-only encoding. Asymptotic Toffoli costs are stated as O(N_k^3) for k-mesh refinement and O(N_a^{3.5}) for supercell enlargement, with resource estimates for bulk diamond claiming an order-of-magnitude reduction relative to prior periodic-solid work.

Significance. If the LCU decomposition and leading-order cost claims hold, the work offers a meaningful advance for simulating strongly correlated periodic materials on fault-tolerant quantum computers. The ability to steer convergence via the more favorable of k-mesh or supercell scaling, combined with basis independence and explicit Brillouin-zone control, addresses key finite-size issues in electronic-structure calculations. The framework's applicability to both plane-wave and localized bases is a strength, and the reported diamond estimates, if verified, would represent a concrete improvement over existing periodic algorithms.

major comments (2)
  1. [§4] §4 (LCU decomposition): The claim that the UPAW-augmented Hamiltonian admits an LCU decomposition whose 1-norm and SELECT/PREPARE costs match the unaugmented Bloch Hamiltonian up to lower-order terms is not supported by an explicit expansion of the augmented operator into unitaries, a bound on the resulting 1-norm, or a circuit diagram showing absorption of the projector terms without extra SELECT calls. This is load-bearing for the stated Toffoli scalings and the order-of-magnitude diamond reduction.
  2. [§5.2] §5.2 (resource estimates): The bulk-diamond Toffoli count and comparison to prior work rest on the unverified assertion that local on-site corrections introduce neither additional terms nor a constant-factor norm increase; without the explicit LCU coefficients or leading-order circuit cost analysis, the headline resource reduction cannot be confirmed.
minor comments (2)
  1. The abstract and introduction should explicitly define N_k and N_a at first use and state the precise baseline (prior work reference and Toffoli count) for the 'order-of-magnitude reduction' claim.
  2. Figure captions for the block-encoding circuit should include gate counts or depth estimates to allow direct comparison with the textual claims.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful review and constructive comments. We address each major comment point by point below and will revise the manuscript to provide the requested explicit details.

read point-by-point responses

  1. Referee: [§4] §4 (LCU decomposition): The claim that the UPAW-augmented Hamiltonian admits an LCU decomposition whose 1-norm and SELECT/PREPARE costs match the unaugmented Bloch Hamiltonian up to lower-order terms is not supported by an explicit expansion of the augmented operator into unitaries, a bound on the resulting 1-norm, or a circuit diagram showing absorption of the projector terms without extra SELECT calls. This is load-bearing for the stated Toffoli scalings and the order-of-magnitude diamond reduction.

    Authors: We thank the referee for this observation. Section 4 derives the LCU decomposition of the Bloch-UPAW Hamiltonian, showing that the strictly local on-site projector corrections transform to additional terms in the Bloch basis that are absorbed into the coefficient list of the existing LCU without new SELECT calls at leading order. The resulting 1-norm is bounded by the unaugmented Bloch norm plus lower-order corrections independent of system size. To strengthen the presentation, the revised manuscript will include an explicit expansion of the augmented operator, a rigorous 1-norm bound, and a circuit diagram demonstrating absorption of the projector terms. These additions will directly support the claimed Toffoli scalings. revision: yes

  2. Referee: [§5.2] §5.2 (resource estimates): The bulk-diamond Toffoli count and comparison to prior work rest on the unverified assertion that local on-site corrections introduce neither additional terms nor a constant-factor norm increase; without the explicit LCU coefficients or leading-order circuit cost analysis, the headline resource reduction cannot be confirmed.

    Authors: The diamond resource estimates in Section 5.2 follow from the leading-order LCU analysis in Section 4, where local corrections contribute neither system-size-dependent terms nor a constant-factor norm increase. The comparison to prior periodic work employs the same qubitization framework but with the improved Bloch-UPAW scaling. We agree that explicit LCU coefficients and a detailed leading-order circuit cost breakdown would permit independent verification of the order-of-magnitude reduction. These elements will be added to the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity in Bloch-UPAW derivation chain

full rationale

The paper constructs the UPAW Hamiltonian in the Bloch basis and derives an LCU decomposition plus block-encoding circuit by applying standard qubitization techniques to this new operator form. The statement that augmentation adds one ancilla and no leading-order Toffolis is presented as a direct consequence of the local on-site correction structure within the decomposition, not as a self-referential definition or fitted input. No equations reduce the target result to its own inputs by construction, no uniqueness theorems are imported from author prior work, and no self-citations are load-bearing for the central claims. Asymptotic scalings O(N_k^3) and O(N_a^{3.5}) and the diamond resource estimate follow from the stated costs without circular reduction. This is the expected non-finding for a paper whose central contribution is an independent Hamiltonian reformulation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities beyond the framework name are detailed. Standard quantum-computing assumptions (qubitization, LCU) are implicitly used.

axioms (1)
  • standard math Standard assumptions of fault-tolerant quantum computing including perfect qubits and availability of block-encoding oracles
    Invoked when deriving the LCU decomposition and block-encoding circuit for the UPAW Hamiltonian.
invented entities (1)
  • Bloch-UPAW framework no independent evidence
    purpose: To combine Bloch k-space structure with strictly local on-site UPAW corrections for periodic electronic structure
    New construction introduced to handle rapid near-nuclear variation while retaining Brillouin-zone sampling control.

pith-pipeline@v0.9.0 · 5586 in / 1257 out tokens · 35744 ms · 2026-05-10T14:54:35.457301+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

52 extracted references · 52 canonical work pages · 1 internal anchor

  1. [1]

    We formulate the many-body UPAW Hamiltonian directly in the Bloch-orbital basis, preserving lat- tice periodicity and the UPAW decomposition into smooth and cusp-like components

    work page

  2. [2]

    UPAW augmentation adds only modest overhead to existing Bloch-orbital cir- cuits

    We derive a linear-combination-of-unitaries (LCU) decomposition and construct the block-encoding circuit for qubitization. UPAW augmentation adds only modest overhead to existing Bloch-orbital cir- cuits

    work page

  3. [3]

    We analyze asymptotic scaling and show that re- fining the momentum mesh costsO(N k) addi- tional qubits andO(N 3 k) Toffoli gates, compared toO(N 1.5 a ) qubits andO(N 3.5 a ) Toffoli gates for equivalent convergence via supercell enlargement. 3

    work page

  4. [4]

    [19] and Rubin et al [20]

    We compute resource estimates for bulk diamond and find roughly an order-of-magnitude reduction in Toffoli count at larger system sizes compared to both Ivanov et al. [19] and Rubin et al [20]. Section II A reviews second quantization for periodic systems and introduces the UPAW construction. Sec- tion III A develops the Bloch–UPAW Hamiltonian, and Sec. I...

    work page

  5. [5]

    Although increasingN b generally entails a proportional increase in the number of plane wavesN pw, in realistic settings one typically hasN pw ≫N b [22, 27]

    Continuum limit Since the basis size is governed by the number of bands perk-point,N b, we formalize the continuum limit as Nb → ∞. Although increasingN b generally entails a proportional increase in the number of plane wavesN pw, in realistic settings one typically hasN pw ≫N b [22, 27]. Accordingly, we fixN pw to be a sufficiently large con- stant. All ...

    work page

  6. [6]

    In the largek- point mesh limit,N k → ∞with fixed supercell size,N pw, andN b

    Thermodynamic limit We consider two thermodynamic limits. In the largek- point mesh limit,N k → ∞with fixed supercell size,N pw, andN b. In the large supercell volume limit,N a → ∞ withN pw andN b scaling proportionally so thatn b = Nb/Na andn pw =N pw/Na remain constant. a. Largek-space limit.In the largek-space limit, the time complexity scales asλ (2) ...

    work page

  7. [7]

    D. W. Berry, A. M. Childs, R. Kothari, R. Cleve, and R. D. Somma, Simulating hamiltonian dynamics with a truncated taylor series, Physical Review Letters114, 090502 (2015)

    work page 2015

  8. [8]

    G. H. Low and I. L. Chuang, Optimal hamiltonian sim- ulation by quantum signal processing, Physical Review Letters118, 010501 (2017)

    work page 2017

  9. [9]

    G. H. Low and I. L. Chuang, Hamiltonian simulation by qubitization, Quantum3, 163 (2019)

    work page 2019

  10. [10]

    J. M. Martyn, Z. M. Rossi, A. K. Tan, and I. L. Chuang, Grand unification of quantum algorithms, PRX Quan- tum2, 040203 (2021)

    work page 2021

  11. [11]

    D. W. Berry, C. Gidney, M. Motta, J. R. McClean, and R. Babbush, Qubitization of arbitrary basis quantum chemistry leveraging sparsity and low rank factorization, Quantum3, 208 (2019)

    work page 2019

  12. [12]

    J. Lee, D. W. Berry, C. Gidney, W. J. Huggins, J. R. Mc- Clean, N. Wiebe, and R. Babbush, Even more efficient quantum computations of chemistry through tensor hy- percontraction, PRX Quantum2, 030305 (2021)

    work page 2021

  13. [13]

    Loaiza and A

    I. Loaiza and A. F. Izmaylov, Block-invariant symme- try shift: Preprocessing technique for second-quantized Hamiltonians to improve their decompositions to linear combination of unitaries, J. Chem. Theory Comput.19, 8201 (2023), arXiv:2304.13772 [quant-ph]

    work page arXiv 2023

  14. [15]

    G. H. Low, R. King, D. W. Berry, Q. Han, A. E. De- Prince, A. F. White, R. Babbush, R. D. Somma, and N. C. Rubin, Fast quantum simulation of electronic struc- ture by spectral amplification, Physical Review X15, 10.1103/pb2g-j9cw (2025)

    work page doi:10.1103/pb2g-j9cw 2025

  15. [16]

    Babbush, D

    R. Babbush, D. W. Berry, J. R. McClean, and H. Neven, Quantum simulation of chemistry with sub- linear scaling in basis size, npj Quantum Information5, 10.1038/s41534-019-0199-y (2019)

    work page doi:10.1038/s41534-019-0199-y 2019

  16. [17]

    Y. Su, D. W. Berry, N. Wiebe, N. Rubin, and R. Bab- bush, Fault-tolerant quantum simulations of chemistry in first quantization, PRX Quantum2, 10.1103/prxquan- tum.2.040332 (2021)

    work page doi:10.1103/prxquan- 2021

  17. [18]

    I. D. Kivlichan, C. Gidney, D. W. Berry, N. Wiebe, J. McClean, W. Sun, Z. Jiang, N. Rubin, A. Fowler, A. Aspuru-Guzik, H. Neven, and R. Babbush, Improved fault-tolerant quantum simulation of condensed-phase correlated electrons via trotterization, Quantum4, 296 (2020)

    work page 2020

  18. [19]

    A. V. Ivanov, C. S¨ underhauf, N. Holzmann, T. Ellaby, R. N. Kerber, G. Jones, and J. Camps, Quantum compu- tation for periodic solids in second quantization, Physical Review Research5, 10.1103/physrevresearch.5.013200 (2023)

    work page doi:10.1103/physrevresearch.5.013200 2023

  19. [20]

    A. J. Cohen, P. Mori-S´ anchez, and W. Yang, Chal- lenges for density functional theory, Chem. Rev.112, 289 (2012)

    work page 2012

  20. [21]

    Shavitt and R

    I. Shavitt and R. J. Bartlett,Many-Body Methods in Chemistry and Physics: MBPT and Coupled-Cluster Theory(Cambridge University Press, Cambridge, UK, 2009)

    work page 2009

  21. [22]

    B. L. Hammond, W. A. Lester, and P. J. Reynolds,Monte Carlo Methods in Ab Initio Quantum Chemistry, World Scientific Lecture and Course Notes in Chemistry, Vol. 1 (World Scientific, Singapore; River Edge, NJ, 1994)

    work page 1994

  22. [23]

    Troyer and U.-J

    M. Troyer and U.-J. Wiese, Computational complexity and fundamental limitations to fermionic quantum monte carlo simulations, Phys. Rev. Lett.94, 170201 (2005)

    work page 2005

  23. [24]

    P. E. Bl¨ ochl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)

    work page 1994

  24. [25]

    A. V. Ivanov, A. Patterson, M. Bothe, C. S¨ underhauf, B. K. Berntson, J. J. Mortensen, M. Kuisma, E. Camp- bell, and R. Izs´ ak, Quantum computation of electronic structure with projector augmented-wave method and plane wave basis set (2025), arXiv:2408.03159 [quant-ph]. 12

    work page arXiv 2025

  25. [26]

    N. C. Rubin, D. W. Berry, F. D. Malone, A. F. White, T. Khattar, A. E. DePrince, S. Sicolo, M. K¨ uehn, M. Kaicher, J. Lee, and R. Babbush, Fault-tolerant quan- tum simulation of materials using bloch orbitals, PRX Quantum4, 040303 (2023)

    work page 2023

  26. [27]

    R. M. Martin,Electronic Structure: Basic Theory and Practical Methods(Cambridge University Press, 2004)

    work page 2004

  27. [28]

    Rostgaard, The projector augmented-wave method (2009), arXiv:0910.1921 [cond-mat.mtrl-sci]

    C. Rostgaard, The projector augmented-wave method (2009), arXiv:0910.1921 [cond-mat.mtrl-sci]

    work page arXiv 2009

  28. [29]

    von Burg, G

    V. von Burg, G. H. Low, T. H¨ aner, D. S. Steiger, M. Rei- her, M. Roetteler, and M. Troyer, Quantum computing enhanced computational catalysis, Phys. Rev. Res.3, 033055 (2021)

    work page 2021

  29. [30]

    The factor 1/8 becomes 1/2 after the spin-sum and an additional factor of 1/2 follows from Chebyshev ampli- tude amplification as introduced in Ref. [23]

    work page

  30. [31]

    Sundararaman and T

    R. Sundararaman and T. A. Arias, Regularization of the coulomb singularity in exact exchange by wigner-seitz truncated interactions: Towards chemical accuracy in nontrivial systems, Phys. Rev. B87, 165122 (2013)

    work page 2013

  31. [32]

    An alternative to this is the Ewald summation approach, which removes the zero mode in reciprocal space via ex- ponential screening and introduces the Madelung con- stant (see Appendix F.2 of Ref. [21]). We employ the Wigner–Seitz scheme here because it typically converges more rapidly and more directly incorporates the relevant exchange physics [25]

    work page

  32. [33]

    Gr¨ uneis, G

    A. Gr¨ uneis, G. Booth, M. Marsman, J. Spencer, A. Alavi, and G. Kresse, Natural orbitals for wave function based correlated calculations using a plane wave basis set, Jour- nal of Chemical Theory and Computation7(2011)

    work page 2011

  33. [34]

    J. J. Mortensen, A. H. Larsen, M. Kuisma, A. V. Ivanov, A. Taghizadeh, A. Peterson, A. Haldar, A. O. Dohn, C. Sch¨ afer, E. ¨O. J´ onsson, E. D. Hermes, F. A. Nils- son, G. Kastlunger, G. Levi, H. J´ onsson, H. H¨ akkinen, J. Fojt, J. Kangsabanik, J. Sødequist, J. Lehtom¨ aki, J. Heske, J. Enkovaara, K. T. Winther, M. Dulak, M. M. Melander, M. Ovesen, M. ...

    work page doi:10.1063/5.0182685 2024

  34. [35]

    Rappoport, N

    D. Rappoport, N. R. M. Crawford, F. Furche, and K. Burke, Approximate density functionals: Which should i choose?, inEncyclopaedia of Inorganic Chem- istry, edited by E. I. Solomon, R. B. King, and R. A. Scott (John Wiley & Sons, Ltd., Chichester; Hoboken, 2009)

    work page 2009

  35. [36]

    [19] useda= 1.1 ˚A; this choice does not affect the asymptotic scaling

    Ref. [19] useda= 1.1 ˚A; this choice does not affect the asymptotic scaling

    work page

  36. [37]

    Here we used GPAW in LCAO mode with a dzp basis

    work page

  37. [38]

    Quantum comput- ing technology roadmaps and capability assessment for sci- entific computing – an analysis of use cases from the nersc workload,

    D. Camps, E. Rrapaj, K. Klymko, H. Kim, K. Gott, S. Darbha, J. Balewski, B. Austin, and N. J. Wright, Quantum computing technology roadmaps and capa- bility assessment for scientific computing – an anal- ysis of use cases from the nersc workload (2025), arXiv:2509.09882 [quant-ph]

    work page arXiv 2025

  38. [39]

    R. King, G. H. Low, R. Babbush, R. D. Somma, and N. C. Rubin, Quantum simulation with sum-of-squares spectral amplification (2025), arXiv:2505.01528 [quant- ph]

    work page arXiv 2025

  39. [40]

    D. W. Berry, K. Wan, A. D. Baczewski, E. C. Ek- lund, A. Tikku, and R. Babbush, Quantum simulation of chemistry via quantum fast multipole method (2025), arXiv:2510.07380 [quant-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  40. [41]

    Taheridehkordi, M

    A. Taheridehkordi, M. Schlipf, Z. Sukurma, M. Humer, A. Gr¨ uneis, and G. Kresse, Phaseless auxiliary field quantum monte carlo with projector-augmented wave method for solids, The Journal of Chemical Physics159, 10.1063/5.0156657 (2023)

    work page doi:10.1063/5.0156657 2023

  41. [42]

    Babbush, C

    R. Babbush, C. Gidney, D. W. Berry, N. Wiebe, J. Mc- Clean, A. Paler, A. Fowler, and H. Neven, Encoding elec- tronic spectra in quantum circuits with linear t complex- ity, Phys. Rev. X8, 041015 (2018). 13 Appendix A: Details on Projector Augmented-Wave method In this section, we provide a detailed overview of the mathematical and physical structure under...

    work page 2018

  42. [43]

    Construction of the projection operator We review the construction of ˆTbecause it supplies the localized operator components that will later be combined with the Bloch/symmetry reductions of Sec. II A. Following Ref. [22], the PAW transformation operator ˆTis constructed as:

    work page

  43. [44]

    Partition the physical system into two spatial regions: (a) a region away from atomic centers (thelatticeregion), (b) atom-centeredaugmentation spheresS 3 a around each atoma

    work page

  44. [45]

    Inside eachS 3 a, one introduces atom-centered partial wavesφ a i (x), smooth partial waves ˜φa i (x), and associatedprojector functions˜p a i (x)

    In the lattice region, states are represented by smooth auxiliary functions ˜ψ(x). Inside eachS 3 a, one introduces atom-centered partial wavesφ a i (x), smooth partial waves ˜φa i (x), and associatedprojector functions˜p a i (x). The smooth partial waves are required to be analytic withinS 3 a and are constructed by matchingφ a i (x) and its derivatives ...

    work page

  45. [46]

    The all-electron KS orbitalψ KS(x) is then written as a smooth background plus augmentation corrections, ψKS(x) = ˜ψ(x) + NaX a=1 naX i=1 (φa i (x)−˜φa i (x)) Z S3a ˜pa i (r) ˜ψ(r)d 3r.(A2)

    work page

  46. [47]

    ,xN) = Ψj(x1,

    Consequently, ˆTadmits the explicit operator decomposition ˆT= ˆI+ X k ˆTk ,(A3) wherek= (a, i) and each local operator ˆTk acts on a test functionf(x) as ˆTkf(x) =χ k Z S3a ˜pk(r)f(r)d 3r,(A4) with χk(x) =φ k(x)−˜φk(x).(A5) Within this framework, the many-body wavefunctions are expressed recursively as Ψj−1(x1, . . . ,xN) = Ψj(x1, . . . ,xN) + X mj χmj (...

    work page

  47. [48]

    In the PAW formalism, the true all-electron density differs from the auxiliary smooth density primarily within localized atomic regions surrounding each nucleus

    The compensation charge We use the indexp= ( ⃗k, i). In the PAW formalism, the true all-electron density differs from the auxiliary smooth density primarily within localized atomic regions surrounding each nucleus. To account for this difference, one introduces acompensation charge, defined as: ˜Z a pq(r) = X L Qa L,pq ˜ga L(⃗ r),(A10) with the expansion ...

    work page

  48. [49]

    PAW C-tensor ThePA W on-site Coulomb correction tensorC a explicitly accounts for the localized electron-electron interactions near atomic nuclei, which are inadequately captured by the smooth pseudo-density alone. Formally, this tensor is a rank-4 object defined as: C a i1i2i3i4 = 1 2 (φa i1 φa i2 |φ a i3 φa i4)−( ˜φa i1 ˜φa i2 |˜φa i3 ˜φa i4) − X L 1 2 ...

    work page

  49. [50]

    The one-body PAW tensor Before proceeding it is important to note that these integrals are calculated under the frozen-core approximation, wherein electrons occupying low-lying core orbitals are assumed inactive with respect to electronic dynamics. Conse- quently, when evaluating the expectation value of an arbitrary operator ˆO, we partition the resultin...

    work page

  50. [51]

    (A15) to simplify the summation overL

    ∆a 0,i1i2 ,(A25) where in the final equality we have employed the orthogonality relation from Eq. (A15) to simplify the summation overL. Finally, the exchange contributionX a,ex i1i2 , arising from valence-core electron interactions, is given by: X a,ex i1i2 = γa/2X j=1 (ζ a j ˜φa i1 |ζ a j ˜φa i2),(A26) where the orbitalsζ a j represent the frozen-core s...

    work page

  51. [52]

    LCU decomposition of the soft two-body term As a first step, we decompose the Hamiltonian into three distinct contributions: a one-body term, a soft two-body interaction, and a hard PAW-specific correction that accounts for the augmentation sphere contributions: ˆH= ˆH (1) + ˆ˜H (2) + ˆH (2) PA W .(B1) Working in the Bloch representation allows us to expl...

    work page

  52. [53]

    offset +µ

    LCU decomposition of the two-body PAW correction term We now turn to deriving the LCU representation for the hard-PAW correction of the Hamiltonian. This term arises from the augmentation-sphere contributions and encodes the difference between the smooth pseudo-density description and the full all-electron Coulomb interaction. Explicitly, it can be writte...

    work page