Toward Large-Scale AFQMC Calculations: Large Time Step Auxiliary-Field Quantum Monte Carlo.

We report modifications of the ph-AFQMC algorithm that allow the use of large time steps and reliable time step extrapolation. Our modified algorithm eliminates size-consistency errors present in the standard algorithm when large time steps are employed. We investigate various methods to approximate the exponential of the one-body operator within the AFQMC framework, distinctly demonstrating the superiority of Krylov methods over the conventional Taylor expansion. We assess various propagators within AFQMC and demonstrate that the Split-2 propagator is the optimal method, exhibiting the smallest time-step errors. For the HEAT set molecules, the time-step extrapolated energies deviate on average by only 0.19 kcal/mol from the accurate small time-step energies. For small water clusters, we obtain accurate complete basis-set binding energies using time-step extrapolation with a mean absolute error of 0.07 kcal/mol compared to CCSD(T). Using large time-step ph-AFQMC for the N2 dimer, we show that accurate bond lengths can be obtained while reducing CPU time by an order of magnitude.

Benchmark Phaseless Auxiliary-Field Quantum Monte Carlo Method for Small Molecules.

We report a scalable Fortran implementation of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) and demonstrate its excellent performance and beneficial scaling with respect to system size. Furthermore, we investigate modifications of the phaseless approximation that can help to reduce the overcorrelation problems common to the ph-AFQMC. We apply the method to the 26 molecules in the HEAT set, the benzene molecule, and water clusters. We observe a mean absolute deviation of the total energy of 1.15 kcal/mol for the molecules in the HEAT set, close to chemical accuracy. For the benzene molecule, the modified algorithm despite using a single-Slater-determinant trial wavefunction yields the same accuracy as the original phaseless scheme with 400 Slater determinants. Despite these improvements, we find systematic errors for the CN, CO2, and O2 molecules that need to be addressed with more accurate trial wavefunctions. For water clusters, we find that the ph-AFQMC yields excellent binding energies that differ from CCSD(T) by typically less than 0.5 kcal/mol.

Phaseless auxiliary field quantum Monte Carlo with projector-augmented wave method for solids.

We implement the phaseless auxiliary field quantum Monte Carlo method using the plane-wave based projector augmented wave method and explore the accuracy and the feasibility of applying our implementation to solids. We use a singular value decomposition to compress the two-body Hamiltonian and, thus, reduce the computational cost. Consistent correlation energies from the primitive-cell sampling and the corresponding supercell calculations numerically verify our implementation. We calculate the equation of state for diamond and the correlation energies for a range of prototypical solid materials. A down-sampling technique along with natural orbitals accelerates the convergence with respect to the number of orbitals and crystal momentum points. We illustrate the competitiveness of our implementation in accuracy and computational cost for dense crystal momentum point meshes compared to a well-established quantum-chemistry approach, the coupled-cluster ansatz including singles, doubles, and perturbative triple particle-hole excitation operators.

Approaching the basis-set limit of the dRPA correlation energy with explicitly correlated and projector augmented-wave methods.

The direct random-phase approximation (dRPA) is used to calculate and compare atomization energies for the HEAT set and ten selected molecules of the G2-1 set using both plane waves and Gaussian-type orbitals. We describe detailed procedures to obtain highly accurate and well converged results for the projector augmented-wave method as implemented in the Vienna Ab initio Simulation Package as well as the explicitly correlated dRPA-F12 method as implemented in the TURBOMOLE package. The two approaches agree within chemical accuracy (1 kcal/mol) for the atomization energies of all considered molecules, both for the exact exchange as well as for the RPA. The root mean-square deviation is 0.41 kcal/mol for the exact exchange (evaluated using density functional theory orbitals) and 0.33 kcal/mol for exact exchange plus correlation from the RPA.

Dynamic Rashba-Dresselhaus Effect.

The Rashba-Dresselhaus effect is the splitting of doubly degenerate band extrema in semiconductors, accompanied by the emergence of counterrotating spin textures and spin-momentum locking. Here we investigate how this effect is modified by lattice vibrations. We show that, in centrosymmetric nonmagnetic crystals, for which a bulk Rashba-Dresselhaus effect is symmetry-forbidden, electron-phonon interactions can induce a phonon-assisted, dynamic Rashba-Dresselhaus spin splitting in the presence of an out-of-equilibrium phonon population. In particular, we show how Rashba, Dresselhaus, or Weyl spin textures can selectively be established by driving coherent infrared-active phonons, and we perform ab initio calculations to quantify this effect for halide perovskites.

Carrier Lifetimes and Polaronic Mass Enhancement in the Hybrid Halide Perovskite CH_{3}NH_{3}PbI_{3} from Multiphonon Fröhlich Coupling.

We elucidate the nature of the electron-phonon interaction in the archetypal hybrid perovskite CH_{3}NH_{3}PbI_{3} using ab initio many-body calculations and an exactly solvable model. We demonstrate that electrons and holes near the band edges primarily interact with three distinct groups of longitudinal-optical vibrations, in order of importance: the stretching of the Pb-I bond, the bending of the Pb-I-Pb bonds, and the libration of the organic cations. These polar phonons induce ultrafast intraband carrier relaxation over timescales of 6-30 fs and yield polaron effective masses 28% heavier than the bare band masses. These findings allow us to rationalize previous experimental observations and provide a key to understanding carrier dynamics in halide perovskites.

Origin of Superconductivity and Latent Charge Density Wave in NbS_{2}.

We elucidate the origin of the phonon-mediated superconductivity in 2H-NbS_{2} using the ab initio anisotropic Migdal-Eliashberg theory including Coulomb interactions. We demonstrate that superconductivity is associated with Fermi surface hot spots exhibiting an unusually strong electron-phonon interaction. The electron-lattice coupling is dominated by low-energy anharmonic phonons, which place the system on the verge of a charge density wave instability. We also provide definitive evidence for two-gap superconductivity in 2H-NbS_{2}, and show that the low- and high-energy peaks observed in tunneling spectra correspond to the Γ- and K-centered Fermi surface pockets, respectively. The present findings call for further efforts to determine whether our proposed mechanism underpins superconductivity in the whole family of metallic transition metal dichalcogenides.

Melting of Pb Charge Glass and Simultaneous Pb-Cr Charge Transfer in PbCrO3 as the Origin of Volume Collapse.

A metal to insulator transition in integer or half integer charge systems can be regarded as crystallization of charges. The insulating state tends to have a glassy nature when randomness or geometrical frustration exists. We report that the charge glass state is realized in a perovskite compound PbCrO3, which has been known for almost 50 years, without any obvious inhomogeneity or triangular arrangement in the charge system. PbCrO3 has a valence state of Pb(2+)(0.5)Pb(4+)(0.5)Cr(3+)O3 with Pb(2+)-Pb(4+) correlation length of three lattice-spacings at ambient condition. A pressure induced melting of charge glass and simultaneous Pb-Cr charge transfer causes an insulator to metal transition and ∼10% volume collapse.

Spin-polarized surface state in EuO(100).

High-quality films of the ferromagnetic semiconductor EuO are grown on epitaxial graphene on Ir(111) and investigated in situ with scanning tunneling microscopy and spectroscopy. Electron scattering at defects leads to standing-wave patterns, manifesting the existence of a surface state in EuO. The surface state is analyzed at different temperatures and energies. We observe a pronounced energy shift of the surface state when cooling down below the Curie temperature TC, which indicates a spin polarization of this state at low temperatures. The experimental results are in agreement with corresponding density functional theory calculations.

Hybrid functionals and GW approximation in the FLAPW method.

We present recent advances in numerical implementations of hybrid functionals and the GW approximation within the full-potential linearized augmented-plane-wave (FLAPW) method. The former is an approximation for the exchange­correlation contribution to the total energy functional in density-functional theory, and the latter is an approximation for the electronic self-energy in the framework of many-body perturbation theory. All implementations employ the mixed product basis, which has evolved into a versatile basis for the products of wave functions, describing the incoming and outgoing states of an electron that is scattered by interacting with another electron. It can thus be used for representing the nonlocal potential in hybrid functionals as well as the screened interaction and related quantities in GW calculations. In particular, the six-dimensional space integrals of the Hamiltonian exchange matrix elements (and exchange self-energy) decompose into sums over vector­matrix­vector products, which can be evaluated easily. The correlation part of the GW self-energy, which contains a time or frequency dependence, is calculated on the imaginary frequency axis with a subsequent analytic continuation to the real axis or, alternatively, by a direct frequency convolution of the Green function G and the dynamically screened Coulomb interaction W along a contour integration path that avoids the poles of the Green function. Hybrid-functional and GW calculations are notoriously computationally expensive. We present a number of tricks that reduce the computational cost considerably, including the use of spatial and time-reversal symmetries, modifications of the mixed product basis with the aim to optimize it for the correlation self-energy and another modification that makes the Coulomb matrix sparse, analytic expansions of the interaction potentials around the point of divergence at k = 0, and a nested density and density-matrix convergence scheme for hybrid-functional calculations. We show CPU timings for prototype semiconductors and illustrative results for GdN and ZnO.