Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userIterative Configuration Interaction with Selection
66
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Even when starting with a very poor initial guess, the iterative configuration interaction (iCI) approach can converge from above to full CI very quickly by constructing and diagonalizing a small Hamiltonian matrix at each macro/micro-iteration. However, iCI scales exponentially with respect to the numbers of electrons and orbitals. The problem can be mitigated by observing that a vast number of configurations have little weights in the wave function and do not contribute to the correlation energy. The real questions are then (a) how to identify important configurations in the early stage of the calculation and (b) how to account for the residual contributions of those unimportant configurations. It is generally true that if a high-quality yet compact variational space can be determined for describing the static correlation, a low-order treatment of the residual dynamic correlation would be sufficient. While this is common to all selected CI schemes, the `iCI with selection scheme presented here has the following distinctive features: (1) Full spin symmetry is maintained. (2) Although the selection is performed on individual CSFs, it is orbital configurations (oCFG) that are used as the organizing units. (3) Given a coefficient pruning-threshold, the selection of important oCFGs/CSFs is performed iteratively until convergence. (4) At each iteration for the growth of wave function, the first-order interacting space is decomposed into disjoint subspaces to reduce memory requirement and facilitate parallelization. (5) Upper bounds for the interactions between oCFG pairs are used to screen each first-order interacting subspace before the first-order coefficients of individual CSFs are evaluated. (6) Upon termination of the selection, dynamic correlation is estimated by using the ENPT2 theory. The efficacy of the iCIPT2 scheme is demonstrated numerically by taking several examples.
rate research
Read More
We extend the recently proposed heat-bath configuration interaction (HCI) method [Holmes, Tubman, Umrigar, J. Chem. Theory Comput. 12, 3674 (2016)], by introducing a semistochastic algorithm for performing multireference Epstein-Nesbet perturbation theory, in order to completely eliminate the severe memory bottleneck of the original method. The proposed algorithm has several attractive features. First, there is no sign problem that plagues several quantum Monte Carlo methods. Second, instead of using Metropolis-Hastings sampling, we use the Alias method to directly sample determinants from the reference wavefunction, thus avoiding correlations between consecutive samples. Third, in addition to removing the memory bottleneck, semistochastic HCI (SHCI) is faster than the deterministic variant for many systems if a stochastic error of 0.1 mHa is acceptable. Fourth, within the SHCI algorithm one can trade memory for a modest increase in computer time. Fifth, the perturbative calculation is embarrassingly parallel. The SHCI algorithm extends the range of applicability of the original algorithm, allowing us to calculate the correlation energy of very large active spaces. We demonstrate this by performing calculations on several first row dimers including F2 with an active space of (14e, 108o), Mn-Salen cluster with an active space of (28e, 22o), and Cr2 dimer with up to a quadruple-zeta basis set with an active space of (12e, 190o). For these systems we were able to obtain better than 1 mHa accuracy with a wall time of merely 55 seconds, 37 seconds, and 56 minutes on 1, 1, and 4 nodes, respectively.
An iterative configuration interaction (iCI)-based multiconfigurational self-consistent field (SCF) theory, iCISCF, is proposed to handle systems that require large complete active spaces (CAS). The success of iCISCF stems from three ingredients: (1) efficient selection of individual configuration state functions spanning the CAS, meanwhile maintaining full spin symmetry; (2) the use of Jacobi rotation for the optimization of active orbitals, in conjunction with a quasi-Newton algorithm for the core/active-virtual and core-active orbital rotations; (3) a second-order perturbative treatment of the residual space left over by the selection procedure (i.e., iCISCF(2)). Just like selected iCI being a very accurate approximation to CASCI, iCISCF(2) is a very accurate approximation to CASSCF. Several examples that go beyond the capability of CASSCF are taken as showcases to reveal the performances of iCISCF and iCISCF(2).
We introduce vibrational heat-bath configuration interaction (VHCI) as an accurate and efficient method for calculating vibrational eigenstates of anharmonic systems. Inspired by its origin in electronic structure theory, VHCI is a selected CI approach that uses a simple criterion to identify important basis states with a pre-sorted list of anharmonic force constants. Screened second-order perturbation theory and simple extrapolation techniques provide significant improvements to variational energy estimates. We benchmark VHCI on four molecules with 12 to 48 degrees of freedom and use anharmonic potential energy surfaces truncated at fourth and sixth order. For all molecules studied, VHCI produces vibrational spectra of tens or hundreds of states with sub-wavenumber accuracy at low computational cost.
Full configuration interaction (FCI) solvers are limited to small basis sets due to their expensive computational costs. An optimal orbital selection for FCI (OptOrbFCI) is proposed to boost the power of existing FCI solvers to pursue the basis set limit under a computational budget. The optimization problem coincides with that of the complete active space SCF method (CASSCF), while OptOrbFCI is algorithmically quite different. OptOrbFCI effectively finds an optimal rotation matrix via solving a constrained optimization problem directly to compress the orbitals of large basis sets to one with a manageable size, conducts FCI calculations only on rotated orbital sets, and produces a variational ground-state energy and its wave function. Coupled with coordinate descent full configuration interaction (CDFCI), we demonstrate the efficiency and accuracy of the method on the carbon dimer and nitrogen dimer under basis sets up to cc-pV5Z. We also benchmark the binding curve of the nitrogen dimer under the cc-pVQZ basis set with 28 selected orbitals, which provide consistently lower ground-state energies than the FCI results under the cc-pVDZ basis set. The dissociation energy in this case is found to be of higher accuracy.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
