No Arabic abstract
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 Greens 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.
We expand upon the recent semi-stochastic adaptation to full configuration interaction quantum Monte Carlo (FCIQMC). We present an alternate method for generating the deterministic space without a priori knowledge of the wave function and present stochastic efficiencies for a variety of both molecular and lattice systems. The algorithmic details of an efficient semi-stochastic implementation are presented, with particular consideration given to the effect that the adaptation has on parallel performance in FCIQMC. We further demonstrate the benefit for calculation of reduced density matrices in FCIQMC through replica sampling, where the semi-stochastic adaptation seems to have even larger efficiency gains. We then combine these ideas to produce explicitly correlated corrected FCIQMC energies for the Beryllium dimer, for which stochastic errors on the order of wavenumber accuracy are achievable.
We present a new approach to calculate excited states with the full configuration interaction quantum Monte Carlo (FCIQMC) method. The approach uses a Gram-Schmidt procedure, instantaneously applied to the stochastically evolving distributions of walkers, to orthogonalize higher energy states against lower energy ones. It can thus be used to study several of the lowest-energy states of a system within the same symmetry. This additional step is particularly simple and computationally inexpensive, requiring only a small change to the underlying FCIQMC algorithm. No trial wave functions or partitioning of the space is needed. The approach should allow excited states to be studied for systems similar to those accessible to the ground-state method, due to a comparable computational cost. As a first application we consider the carbon dimer in basis sets up to quadruple-zeta quality, and compare to existing results where available.
We study several approaches to orbital optimization in selected configuration interaction plus perturbation theory (SCI+PT) methods, and test them on the ground and excited states of three molecules using the semistochastic heatbath configuration interaction (SHCI) method. We discuss the ways in which the orbital optimization problem in SCI resembles and differs from that in complete active space self-consistent field (CASSCF). Starting from natural orbitals, these approaches divide into three classes of optimization methods according to how they treat coupling between configuration interaction (CI) coefficients and orbital parameters, namely uncoupled, fully coupled, and quasi-fully coupled methods. We demonstrate that taking the coupling into account is crucial for fast convergence and recommend two quasi-fully coupled methods for such applications: accelerated diagonal Newton and Broyden-Fletcher-Goldfarb-Shanno (BFGS).
Computational modeling is an important aspect of the research on nuclear waste materials. In particular, atomistic simulations, when used complementary to experimental efforts, contribute to the scientific basis of safety case for nuclear waste repositories. Here we discuss the state-of-the-art and perspectives of atomistic modeling for nuclear waste management on a few cases of successful synergy of atomistic simulations and experiments. In particular, we discuss here: (1) the potential of atomistic simulations to investigate the uranium oxidation state in mixed valence uranium oxides and (2) the ability of cementitious barrier materials to retain radionuclides such as 226Ra and 90Sr, and of studtite/metastudtite secondary peroxide phases to incorporate actinides such as Np and Am. The new contribution we make here is the computation of the incorporation of Sr by C-S-H (calcium silicate hydrate) phases.
We present a Flexible Ansatz for N-body Configuration Interaction (FANCI) that includes any multideterminant wavefunction. This ansatz is a generalization of the Configuration Interaction (CI) wavefunction, where the coefficients are replaced by a specified function of certain parameters. By making an appropriate choice for this function, we can reproduce popular wavefunction structures like CI, Coupled-Cluster, Tensor Network States, and geminal-product wavefunctions. The universality of this framework suggests a programming structure that allows for the easy construction and optimization of arbitrary wavefunctions. Here, we will discuss the structures of the FANCI framework and its implications for wavefunction properties, particularly accuracy, cost, and size-consistency. We demonstrate the flexibility of this framework by reconstructing popular wavefunction ans{a}tze and modifying them to construct novel wavefunction forms. FANCI provides a powerful framework for exploring, developing, and testing new wavefunction forms.