Source author record

Qiming Sun

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

ResearcherUnclaimed source record

physics.chem-ph cond-mat.str-el physics.comp-ph quant-ph cond-mat.mtrl-sci

Catalog footprint

What is connected

8works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Relativistic Self-Consistent $GW$: Exact Two-Component Formalism with One-Electron Approximation for Solids

We present a formulation of relativistic self-consistent $GW$ for solids based on the exact two-component formalism with one-electron approximation (X2C1e) and non-relativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and the interplay of relativistic effects with electron correlation without adjustable parameters. Our all-electron implementation is fully $ab$ $initio$ and does not require a pseudopotential constructed from atomic calculations. We examine the effect of the X2C1e approximation by comparison to the established four-component formalism and reach excellent agreement. The simplicity of X2C1e enables the construction of higher order theories, such as embedding theories, on top of perturbative calculations.

preprint2020arXiv

Accurate many-body electronic structure near the basis set limit: application to the chromium dimer

We describe a method for computing near-exact energies for correlated systems with large Hilbert spaces. The method efficiently identifies the most important basis states (Slater determinants) and performs a variational calculation in the subspace spanned by these determinants. A semistochastic approach is then used to add a perturbative correction to the variational energy to compute the total energy. The size of the variational space is progressively increased until the total energy converges to within the desired tolerance. We demonstrate the power of the method by computing a near-exact potential energy curve (PEC) for a very challenging molecule -- the chromium dimer.

preprint2019arXiv

Electronic landscape of the P-cluster of nitrogenase as revealed through many-electron quantum wavefunctions

The electronic structure of the nitrogenase metal cofactors is central to nitrogen fixation. However, the P-cluster and iron molybdenum cofactor, each containing eight irons, have resisted detailed characterization of their electronic properties. Through exhaustive many-electron wavefunction simulations enabled by new theoretical methods, we report on the low-energy electronic states of the P-cluster in three oxidation states. The energy scales of orbital and spin excitations overlap, yielding a dense spectrum with features we trace to the underlying atomic states and recouplings. The clusters exist in superpositions of spin configurations with non-classical spin correlations, complicating interpretation of magnetic spectroscopies, while the charges are mostly localized from reorganization of the cluster and its surroundings. Upon oxidation, the opening of the P-cluster significantly increases the density of states, which is intriguing given its proposed role in electron transfer. These results demonstrate that many-electron simulations stand to provide new insights into the electronic structure of the nitrogenase cofactors.

preprint2019arXiv

Electronic structure of bulk manganese oxide and nickel oxide from coupled cluster theory

We describe the ground- and excited-state electronic structure of bulk MnO and NiO, two prototypical correlated electron materials, using coupled cluster theory with single and double excitations (CCSD). As a corollary, this work also reports the first implementation of unrestricted periodic ab initio equation-of motion CCSD. Starting from a Hartree-Fock reference, we find fundamental gaps of 3.46 eV and 4.83 eV for MnO and NiO respectively for the 16 unit supercell, slightly overestimated compared to experiment, although finite-size scaling suggests that the gap is more severely overestimated in the thermodynamic limit. From the character of the correlated electronic bands we find both MnO and NiO to lie in the intermediate Mott/charge-transfer insulator regime, although NiO appears as a charge transfer insulator when only the fundamental gap is considered. While the lowest quasiparticle excitations are of metal 3d and O 2p character in most of the Brillouin zone, near the Γ point, the lowest conduction band quasiparticles are of s character. Our study supports the potential of coupled cluster theory to provide high level many-body insights into correlated solids.

preprint2016arXiv

A practical guide to density matrix embedding theory in quantum chemistry

Density matrix embedding theory (DMET) provides a theoretical framework to treat finite fragments in the presence of a surrounding molecular or bulk environment, even when there is significant correlation or entanglement between the two. In this work, we give a practically oriented and explicit description of the numerical and theoretical formulation of DMET. We also describe in detail how to perform self-consistent DMET optimizations. We explore different embedding strategies with and without a self-consistency condition in hydrogen rings, beryllium rings, and a sample S$_{\text{N}}$2 reaction. The source code for the calculations in this work can be obtained from \url{https://github.com/sebwouters/qc-dmet}.

preprint2016arXiv

N-electron valence state perturbation theory based on a density matrix renormalization group reference function, with applications to the chromium dimer and poly-p-phenylene vinylene oligomer

The strongly-contracted variant of second order N -electron valence state perturbation theory (NEVPT2) is an efficient perturbative method to treat dynamic correlation without the problems of intruder states or level shifts, while the density matrix renormalization group (DMRG) provides the capability to tackle static correlation in large active spaces. We present a combination of the DMRG and strongly-contracted NEVPT2 (DMRG-SC-NEVPT2) that uses an efficient algorithm to compute high order reduced density matrices from DMRG wave functions. The capabilities of DMRG-SC-NEVPT2 are demonstrated on calculations of the chromium dimer potential energy curve at the basis set limit, and the excitation energies of poly-p-phenylene vinylene trimer (PPV(n=3)).

preprint2015arXiv

Stochastic multi-configurational self-consistent field theory

The multi-configurational self-consistent field theory is considered the standard starting point for almost all multireference approaches required for strongly-correlated molecular problems. The limitation of the approach is generally given by the number of strongly-correlated orbitals in the molecule, as its cost will grow exponentially with this number. We present a new multi-configurational self-consistent field approach, wherein linear determinant coefficients of a multi-configurational wavefunction are optimized via the stochastic full configuration interaction quantum Monte Carlo technique at greatly reduced computational cost, with non-linear orbital rotation parameters updated variationally based on this sampled wavefunction. This extends this approach to strongly-correlated systems with far larger active spaces than it is possible to treat by conventional means. By comparison with this traditional approach, we demonstrate that the introduction of stochastic noise in both the determinant amplitudes and the gradient and Hessian of the orbital rotations does not preclude robust and reliable convergence of the orbital optimization. It can even improve the ability to avoid convergence to local minima in the orbital space, and therefore aid in finding variationally lower-energy solutions. We consider the effect on the convergence of the orbitals as the number of walkers and the sampling time within the active space increases, as well as the effect on the final energy and error. The scope of the new protocol is demonstrated with a study of the increasingly strongly correlated electronic structure in a series of polycyclic aromatic hydrocarbons, up to the large coronene molecule in a complete active space of 24 pi electrons in 24 orbitals, requiring only modest computational resources.

preprint2014arXiv

Libcint: An efficient general integral library for Gaussian basis functions

An efficient integral library Libcint was designed to automatically implement general integrals for Gaussian-type scalar and spinor basis functions. The library can handle arbitrary integral expressions on top of $\mathbf{p}$, $\mathbf{r}$ and $σ$ operators with one-electron overlap and nuclear attraction, two-electron Coulomb and Gaunt operators. Using a symbolic algebra tool, new integrals are derived and translated to C code programmatically. The generated integrals can be used in various types of molecular properties. In the present work, we computed the analytical gradients and NMR shielding constants at both non-relativistic and four-component relativistic Hartree-Fock level to demonstrate the capability of the integral library. Due to the use of kinetically balanced basis and gauge including atomic orbitals, the relativistic analytical gradients and shielding constants requires the integral library to handle the fifth-order electron repulsion integral derivatives. The generality of the integral library is achieved without losing efficiency. On the modern multi-CPU platform, Libcint can easily reach the overall throughput being many times of the I/O bandwidth. On a 20-core node, we are able to achieve an average output 7.9 GB/s for C$_{60}$ molecule with cc-pVTZ basis.