Search

Examines the excited states of self- and size-consistent configuration interactions observed in a complete active space. Overview of the related studies on configuration interactions; Description of formalism and interpretation; Analysis of CH[sub 2] and BH excitation energies.

Reports on the finding that the mean value total dressing (MVTD) method had shown that MVTD energies for closed shell systems are generally better than coupled cluster approach. Method's basis on the single reference intermediate Hamiltonian theory; Separate contributions of linked and unlinked triples and quadruples.

C 4n cyclacenes exhibit strong bond-alternation in their equilibrium geometry. In the two equivalent geometries, the system keeps an essentially closed-shell character. The two energy minima are separated by a transition state suppressing the bond-alternation, where the wave function is strongly diradical. This paper discusses the physical factors involved in this energy difference and possible evaluations of the barrier height. The barrier given as the energy difference between the restricted density functional theory (DFT)/B3LYP for the equilibrium and the broken symmetry DFT/B3LYP of the transition state is either negative or small, in contradiction with the most reliable Wave Function Theory calculations. The minimal (two electrons in two molecular orbitals) Complete Active Space self-consistent field (CASSCF) overestimates the barrier, and the subsequent second-order perturbation cancels it. Due to the collective character of the spin-polarization effect, it is necessary to perform a full π CASSCF + second-order perturbation to reach a reasonable value of the barrier, but this type of treatment cannot be applied to large molecules. DFT procedures treating on an equal foot the closed-shell and open-shell geometries have been explored, such as Mixed-Reference Spin-Flip Time-dependent-DFT and a new spin-decontamination proposal, namely, DFT-dressed configuration interaction, but the results still depend on the density functional. M06-2X without or with spin-decontamination gives the best agreement with the accurate wave function results.

Benchmarking calculations on excited states of models of phenylalanine protein chains are presented to assess the ability of alternative methods to the standard and most commonly used multiconfigurational wave function-based method, the complete active space self-consistent field (CASSCF), in recovering the non-dynamical correlation for systems that become not affordable by the CASSCF. The exploration of larger active spaces beyond the CASSCF limit is benchmarked through three strategies based on the reduction in the number of determinants: the restricted active space self-consistent field, the generalized active space self-consistent field (GASSCF), and the occupation-restricted multiple active space (ORMAS) schemes. The remaining dynamic correlation effects are then added by the complete active space second-order perturbation theory and by the multireference difference dedicated configuration interaction methods. In parallel, the approximate second-order coupled cluster (CC2), already proven to be successful for small building blocks of model proteins in one of our previous works [Ben Amor et al., J. Chem. Phys. 148, 184105 (2018)], is investigated to assess its performances for larger systems. Among the different alternative strategies to CASSCF, our results highlight the greatest efficiency of the GASSCF and ORMAS schemes in the systematic reduction of the configuration interaction expansion without loss of accuracy in both nature and excitation energies of both singlet ππ * and nπ * CO excited states with respect to the equivalent CASSCF calculations. Guidelines for an optimum applicability of this scheme to systems requiring active spaces beyond the complete active space limit are then proposed. Finally, the extension of the CC2 method to such large systems without loss of accuracy is demonstrated, highlighting the great potential of this method to treat accurately excited states, mainly single reference, of very large systems.

Taking as an example the simple CH 3 radical, this work demonstrates the cooperative character of the spin-polarization phenomenon of the closed-shell core in free radicals. Spin polarization of CH σ bonds is not additive here, as spin polarization of one bond enhances that of the next bond. This cooperativity is demonstrated by a series of configuration interaction calculations converging to the full valence limit and is rationalized by analytic developments. The same phenomenon is shown to take place in those diradicals where spin polarization plays a major role, as illustrated in square planar carbo-cyclobutadiene C 12 H 4 . The treatment of cooperativity represents a challenge for usual post-Hatree-Fock methods.

A benchmark set of relevant geometries of a model protein, the N-acetylphenylalanylamide, is presented to assess the validity of the approximate second-order coupled cluster (CC2) method in studying low-lying excited states of such bio-relevant systems. The studies comprise investigations of basis-set dependence as well as comparison with two multireference methods, the multistate complete active space 2nd order perturbation theory (MS-CASPT2) and the multireference difference dedicated configuration interaction (DDCI) methods. First of all, the applicability and the accuracy of the quasi-linear multireference difference dedicated configuration interaction method have been demonstrated on bio-relevant systems by comparison with the results obtained by the standard MS-CASPT2. Second, both the nature and excitation energy of the first low-lying excited state obtained at the CC2 level are very close to the Davidson corrected CAS+DDCI ones, the mean absolute deviation on the excitation energy being equal to 0.1 eV with a maximum of less than 0.2 eV. Finally, for the following low-lying excited states, if the nature is always well reproduced at the CC2 level, the differences on excitation energies become more important and can depend on the geometry.

A new multireference configuration interaction method using localised orbitals is proposed, in which a molecular system is divided into regions of unequal importance. The advantage of dealing with local orbitals, i.e., the possibility to neglect long range interaction is enhanced. Indeed, while in the zone of the molecule where the important phenomena occur, the interaction cut off may be as small as necessary to get relevant results, in the most part of the system it can be taken rather large, so that results of good quality may be obtained at a lower cost. The method is tested on several systems. In one of them, the definition of the various regions is not based on topological considerations, but on the nature, σ or π, of the localised orbitals, which puts in evidence the generality of the approach.

A selected multireference configuration interaction (CI) method and the corresponding code are presented. It is based on a procedure of localization that permits to obtain well localized occupied and virtual orbitals. Due to the local character of the electron correlation, using local orbitals allows one to neglect long range interactions. In a first step, three topological matrices are constructed, which determine whether two orbitals must be considered as interacting or not. Two of them concern the truncation of the determinant basis, one for occupied/virtual, the second one for dispersive interactions. The third one concerns the truncation of the list of two electron integrals. This approach permits a fine analysis of each kind of approximation and induces a huge reduction of the CI size and of the computational time. The procedure is tested on linear polyene aldehyde chains, dissociation potential energy curve, and reaction energy of a pesticide-Ca(2+) complex and finally on transition energies of a large iron system presenting a light-induced excited spin-state trapping effect.

We report fixed-node diffusion Monte Carlo (FN-DMC) calculations of the singlet n-->pi( *) (CO) vertical transition of acrolein. The impact of the fixed-node approximation on the excitation energy is investigated. To do that, trial wave functions corresponding to various nodal patterns are used. They are constructed by using either a minimal complete-active-space self-consistent field (CASSCF) calculation involving an oxygen lone pair n and the pi( *) (CO) molecular orbitals or a more complete set involving all the molecular orbitals expected to play a significant role in the excitation process. Calculations of both states have been performed with molecular orbitals optimized separately for each state via standard "state specific" CASSCF calculations or by using a common set of optimized orbitals ["state averaged" CASSCF calculations] whose effect is to introduce some important correlation between the nodal patterns of the two electronic states. To investigate the role of the basis set three different basis of increasing size have been employed. The comparative study based on the use of all possible combinations of basis sets, active spaces, and type of optimized molecular orbitals shows that the nodal error on the difference of energies is small when chemically relevant active space and state-averaged-type CASSCF wave functions are used, although the fixed-node error on the individual total energies involved can vary substantially. This remarkable result obtained for the acrolein suggests that FN-DMC calculations based on a simple strategy (use of standard ab initio wave functions and no Monte Carlo optimization of molecular orbital parameters) could be a working computational tool for computing electronic transition energies for more general systems.

The present paper presents a revised version of a size-consistency correction to the multireference configuration interaction techniques previously proposed by Szalay et al. [J. Phys. Chem. 100, 6288 (1996)]. The method assumes a complete active space reference space and separates the nonreference determinants in several classes according to their number of inactive holes and particles. The correction is formulated as a dressing of the diagonal energies of these determinants, which depends on their class, as originally proposed by Ruttink et al. [J. Chem. Phys. 94, 7212 (1991)]. The exclusion principle violating corrections are evaluated through a simple counting of the various excitation processes which remain possible on each class. The efficiency of the method has been tested on a series of multireference problems for which full configuration interaction results are available (OH(2) bond breaking, Be insertion in H(2), excited states of CH(2)). The dressing of a given state not only provides excellent results for this state but also provides accurate excited roots. The efficiency of state-specific dressings is dramatic. The adaptation of this proposal to difference-dedicated configuration interactions can be extremely fruitful, as illustrated in the calculation of the 1 (1)A(g)-1(1)B(u)(pi->pi(*)) transition energy of the trans-butadiene molecule.