Multi-scale multireference configuration interaction calculations for large systems using localized orbitals: Partition in zones

Model protein excited states: MRCI calculations with large active spaces vs CC2 method

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.

Generalization of Intrinsic Orbitals to Kramers-Paired Quaternion Spinors, Molecular Fragments, and Valence Virtual Spinors

Localization of molecular orbitals finds its importance in the representation of chemical bonding (and antibonding) and in the local correlation treatments beyond mean-field approximation. In this paper, we generalize the intrinsic atomic and bonding orbitals [G. Knizia, J. Chem. Theory Comput. 2013, 9, 11, 4834-4843] to relativistic applications using complex and quaternion spinors, as well as to molecular fragments instead of atomic fragments only. By performing a singular value decomposition, we show how localized valence virtual orbitals can be expressed on this intrinsic minimal basis. We demonstrate our method on systems of increasing complexity, starting from simple cases such as benzene, acrylic acid, and ferrocene molecules, and then demonstrate the use of molecular fragments and inclusion of relativistic effects for complexes containing heavy elements such as tellurium, iridium, and astatine. The aforementioned scheme is implemented into a standalone program interfaced with several different quantum chemistry packages.

Electronic Structure and Magnetic Interactions in the Radical Salt [BEDT-TTF]2[CuCl4]

The magnetic behavior and electric properties of the hybrid radical salt [BEDT-TTF]₂[CuCl₄] have been revisited through extended experimental analyses and DDCI and periodic DFT plane waves calculations. Single crystal X-ray diffraction data have been collected at different temperatures, discovering a phase transition occurring in the 250-300 K range. The calculations indicate the presence of intradimer, interdimer, and organic-inorganic π-d interactions in the crystal, a magnetic pattern much more complex than the Bleaney-Bowers model initially assigned to this material. Although this simple model was good enough to reproduce the magnetic susceptibility data, our calculations demonstrate that the actual magnetic structure is significantly more intricate, with alternating antiferromagnetic 1D chains of the organic BEDT-TTF⁺ radical, connected through weak antiferromagnetic interactions with the CuCl₄^2- ions. Combination of experiment and theory allowed us to unambiguously determine and quantify the leading magnetic interactions in the system. The density-of-states curves confirm the semiconductor nature of the system and the dominant organic contribution of the valence and conduction band edges. This general and combined approach appears to be fundamental in order to properly understand the magnetic structure of these complex materials, where experimental data can actually be fitted from a variety of models and parameters.

Switching magnetic interactions in the NiFe Prussian Blue Analogue: an ab initio inspection

Wave function embedded cluster calculations are performed to stress a cation-coupled electron transfer process in the NiFe Prussian Blue Analogue.