Dynamical configuration interaction: Quantum embedding that combines wave functions and Green's functions

We present the concept, derivation, and implementation of dynamical configuration interaction, a quantum embedding theory that combines Green's function methodology with the many-body wave function. In a strongly-correlated active space, we use full configuration interaction (CI) to describe static correlation exactly. We add energy dependent corrections to the CI Hamiltonian which, in principle, include all remaining correlation derived from the bath space surrounding the active space. Next, we replace the exact Hamiltonian in the bath with one of excitations defined over a correlated ground state. This transformation is naturally suited to the methodology of many-body Green's functions. In this space, we use a modified $GW$/Bethe-Salpeter equation procedure to calculate excitation energies. Combined with an estimate of the ground state energy in the bath, we can efficiently compute the energy dependent corrections, which correlate the full set of orbitals, for very low computational cost. We present dimer dissociation curves for H$_2$ and N$_2$ in good agreement with exact results. Additionally, excited states of N$_2$ and C$_2$ are in excellent agreement with benchmark theory and experiment. By combining the strengths of two disciplines, we achieve a balanced description of static and dynamic correlation in a fully ab-initio, systematically improvable framework

simecol: An Object-Oriented Framework for Ecological Modeling in R

The simecol package provides an open structure to implement, simulate and share ecological models. A generalized object-oriented architecture improves readability and potential code re-use of models and makes simecol-models freely extendable and simple to use. The simecol package was implemented in the S4 class system of the programming language R. Reference applications, e.g. predator-prey models or grid models are provided which can be used as a starting point for own developments. Compact example applications and the complete code of an individual-based model of the water flea Daphnia document the efficient usage of simecol for various purposes in ecological modeling, e.g. scenario analysis, stochastic simulations and individual based population dynamics. Ecologists are encouraged to exploit the abilities of simecol to structure their work and to use R and object-oriented programming as a suitable medium for the distribution and share of ecological modeling code.

A quantum embedding theory in the screened Coulomb interaction: Combining configuration interaction with GW/BSE

We present a new quantum embedding theory called dynamical configuration interaction (DCI) that combines wave function and Green's function theories. DCI captures static correlation in a correlated subspace with configuration interaction and couples to high-energy, dynamic correlation outside the subspace with many-body perturbation theory based on Green's functions. In the correlated subspace, we use a wave function description to avoid embedding the two-particle vertex, which greatly simplifies the frequency structure of the embedding. DCI takes the strengths of both theories to balance static and dynamic correlation in a single, fully ab-initio embedding concept. We show that treating high-energy correlation up to the $GW$ and Bethe-Salpeter equation level is sufficient even for challenging multi-reference problems. Our theory treats ground and excited states on equal footing, and we compute the dissociation curve of N$_2$, vertical excitation energies of N$_2$ and C$_2$, and the ionization spectrum of benzene in excellent agreement with high level quantum chemistry methods and experiment

Thermodynamic equilibrium conditions of graphene films on SiC

First-principles surface phase diagrams reveal that epitaxial monolayer graphene films on the Si side of 3C-SiC(111) can exist as thermodynamically stable phases in a narrow range of experimentally controllable conditions, defining a path to the highest-quality graphene films. Our calculations are based on a van der Waals corrected density functional. The full, experimentally observed (6 sqrt(3)x 6 sqrt(3))-R30 supercells for zero- to trilayer graphene are essential to describe the correct interface geometries and the relative stability of surface phases and possible defects

Piecewise linearity in the GWGW approximation for accurate quasiparticle energy predictions

We identify the deviation from the straight line error (DSLE) -- i.e., the spurious non-linearity of the total energy as a function of fractional particle number -- as the main source for the discrepancy between experimental vertical ionization energies and theoretical quasiparticle energies, as obtained from the $GW$ and $GW$+SOSEX approximations to many-body perturbation theory (MBPT). For self-consistent calculations, we show that $GW$ suffers from a small DSLE. Conversely, for perturbative $G_0W_0$ and $G_0W_0$+SOSEX calculations the DSLE depends on the starting point. We exploit this starting-point dependence to reduce (or completely eliminate) the DSLE. We find that the agreement with experiment increases as the DSLE reduces. DSLE-minimized schemes, thus, emerge as promising avenues for future developments in MBPT

Exact-exchange based quasiparticle energy calculations for the band gap, effective masses and deformation potentials of ScN

The band gaps, longitudinal and transverse effective masses, and deformation potentials of ScN in the rock-salt structure have been calculated employing G0W0-quasiparticle calculations using exact-exchange Kohn-Sham density functional theory one-particle wavefunctions and energies as input. Our quasiparticle gaps support recent experimental observations that ScN has a much lower indirect band gap than previously thought. The results are analyzed in terms of the influence of different approximations for exchange and correlation taken in the computational approach on the electronic structure of ScN.Comment: Comments: 9 pages including 5 figures; related publications can be found at http://www.fhi-berlin.mpg.de/th/th.htm

Self-Consistent Green Function Embedding for Advanced Electronic Structure Methods Based on a Dynamical Mean-Field Concept

We present an embedding scheme for periodic systems that facilitates the treatment of the physically important part (here the unit cell) with advanced electronic-structure methods, that are computationally too expensive for periodic systems. The rest of the periodic system is treated with computationally less demanding approaches, e.g., Kohn-Sham density-functional theory, in a self- consistent manner. Our scheme is based on the concept of dynamical mean-field theory (DMFT) formulated in terms of Green functions. In contrast to the original DMFT formulation for correlated model Hamiltonians, we here consider the unit cell as local embedded cluster in a first-principles way, that includes all electronic degrees of freedom. Our real-space dynamical mean-field embedding (RDMFE) scheme features two nested Dyson equations, one for the embedded cluster and another for the periodic surrounding. The total energy is computed from the resulting Green functions. The performance of our scheme is demonstrated by treating the embedded region with hybrid functionals and many-body perturbation theory in the GW approach for simple bulk systems. The total energy and the density of states converge rapidly with respect to the computational parameters and approach their bulk limit with increasing cluster (i.e., unit cell) size

Supersymmetric Boundary Conditions for the N=2 Sigma Model

We clarify the discussion of N=2 supersymmetric boundary conditions for the classical d=2, N=(2,2) Non-Linear Sigma Model on an infinite strip. Our conclusions about the supersymmetric cycles match the results found in the literature. However, we find a constraint on the boundary action that is not satisfied by many boundary actions used in the literature.Comment: 23 pages, no figures. minor changes in text, a reference adde

Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schr\"odinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.Comment: 16 pages, 11 figure