chARpack: The Chemistry Augmented Reality Package.

Off-loading visualization and interaction into virtual reality (VR) using head-mounted displays (HMDs) has gained considerable popularity in simulation sciences, particularly in chemical modeling. Because of its unique way of soft immersion, augmented reality (AR) HMD technology has even more potential to be integrated into the everyday workflow of computational chemists. In this work, we present our environment to explore the prospects of AR in chemistry and general molecular sciences: The chemistry in Augmented Reality package (chARpack). Besides providing an extensible framework, our software focuses on a seamless transition between a 3D stereoscopic view with true 3D interactions and the traditional desktop PC setup to provide users with the best setup for all tasks in their workflow. Using feedback from domain experts, we discuss our design requirements for this kind of hybrid working environment (AR + PC), regarding input, features, degree of immersion, and collaboration.

Air-Stable Dinuclear Complexes of Four-Coordinate ZnII and NiII Ions with a Radical Bridge: A Detailed Look at Redox Activity and Antiferromagnetic Coupling.

Air-stable dinuclear complexes [(bmsab)NiII(tmsab)NiII(bmsab)]3- and [(bmsab)ZnII(tmsab)ZnII(bmsab)]3- (bmsab = bis(methanesulfoneamido)benzene, tmsab = tetra(methanesulfonamido)benzene) were prepared via a synthetic route based on heteroleptic precursor complexes. The new complexes combine a distorted tetrahedral coordination environment with an open-shell bridging ligand. The ZnII species was subjected to a detailed investigation of the (spectro-)electrochemical processes. The NiII species is a rare example of a complex that combines strong exchange coupling (J > 440 cm-1) with pronounced positive zero-field splitting (D = +72 cm-1). Combining SQUID magnetometry and (HF)EPR spectroscopy with ab initio calculations allowed for accurate quantification of the exchange interaction.

Performance Tests of the Second-Order Approximate Internally Contracted Multireference Coupled-Cluster Singles and Doubles Method icMRCC2.

Benchmark results are presented for the second-order approximation of the internally contracted multireference coupled-cluster method with single and double excitations, icMRCC2 [Köhn, Bargholz, J. Chem. Phys. 2019, 151, 041106], which was designed as a multireference analogue of the single-reference second-order approximate coupled-cluster method CC2 [Christiansen, Koch, Jørgensen, Chem. Phys. Lett. 1995, 243, 409-418]. Vertical excitation energies of various small to medium-sized organic molecules are investigated based on established test sets from the literature. Additionally, the spectroscopic constants of ground and excited states of diatomics and the geometric parameters of excited triatomic molecules were determined and compared to the experimental data. The results show that the method clearly extends the applicability of single-reference CC2, including doubly excited states, and also artifacts of CC2 like too low Rydberg excitations and too weak multiple bonds are eliminated. The method is computationally more demanding than standard multireference second-order perturbation theories but improves significantly in accuracy, as shown by the benchmark results. In addition, it is demonstrated that small active spaces are often sufficient to obtain accurate energies with icMRCC2. Example applications like the automerization of cyclobutadiene, the deactivation pathway of ethylene, and the excited states of an iron complex with a noninnocent nitrosyl ligand demonstrate the potential of icMRCC2 in cases with strong multireference character.

Toward an efficient implementation of internally contracted coupled-cluster methods.

A new implementation of the internally contracted multireference coupled-cluster with singles and doubles (icMRCCSD) method is presented. The new code employs an efficient tensor contraction kernel and can also avoid full four-external integral transformations, which significantly extends the scope of the applicability of icMRCCSD. The new implementation is currently restricted to the simple case of two active electrons in two orbitals and also supports the computation of spin-adapted doublet and triplet coupled-cluster wavefunctions. This contribution describes the basic approach for the automated derivation of working equations and benchmarks the current code against efficient implementations of standard methods, such as single-reference coupled-cluster singles and doubles (CCSD) and internally contracted multireference configuration interaction (icMRCI). Run times for linearized variants of icMRCCSD are only twice as long as comparable CCSD runs and similar to those of the icMRCI implementation, while non-linear terms of more complete variants of icMRCCSD lead to an order of magnitude longer computation times. Nevertheless, the new code allows for computations at larger scales than it was possible previously, with less demands on memory and disk-space resources. This is exemplified by numerical structure optimizations and harmonic force field determinations of NC2H5 isomers and the singlet and triplet states of m-benzyne. In addition, the exchange coupling of a dinuclear copper complex is determined. This work also defines a new commutator approximation for icMRCCSD, which includes all terms that are also present in the single-reference CCSD method, thus yielding a consistent pair of single-reference and multireference coupled-cluster methods.

Capabilities and limits of the unitary coupled-cluster approach with generalized two-body cluster operators.

Unitary cluster expansions of the electronic wavefunction have recently gained much interest because of their use in conjunction with quantum algorithms. In this contribution, we investigate some aspects of an ansatz, using generalized two-body excitation operators, which have been considered in some recent studies on quantum algorithms for quantum chemistry. Our numerical results show that, in particular, two-body operators with effective particle-hole excitation level of one in connection with the usual particle-hole double excitation operators lead to a very accurate, yet compact representation of the wavefunction. Generalized two-body operators with effective excitation rank zero have a considerably less pronounced effect. We compare with standard and unitary coupled-cluster expansions and show that the above mentioned approach matches or even surpasses the accuracy of expansions with three-body particle-hole excitations, in particular at the onset of strong correlation. A downside of the approach is that it is rather difficult to rigorously converge it to its variational minimum.

Thermally Averaged Magnetic Anisotropy Tensors via Machine Learning Based on Gaussian Moments.

We propose a machine learning method to model molecular tensorial quantities, namely, the magnetic anisotropy tensor, based on the Gaussian moment neural network approach. We demonstrate that the proposed methodology can achieve an accuracy of 0.3-0.4 cm-1 and has excellent generalization capability for out-of-sample configurations. Moreover, in combination with machine-learned interatomic potential energies based on Gaussian moments, our approach can be applied to study the dynamic behavior of magnetic anisotropy tensors and provide a unique insight into spin-phonon relaxation.

Limitations of perturbative coupled-cluster approximations for highly accurate investigations of Rb2.

We reveal limitations of several standard coupled-cluster (CC) methods with perturbation-theory based noniterative or approximate iterative treatments of triple excitations when applied to the determination of highly accurate potential energy curves (PECs) of ionic dimers, such as the XΣg+2 electronic ground state of Rb2 +. Such computations are of current interest for the understanding of ion-atom interactions in the ultracold regime. We demonstrate that these CC methods lead to an unphysical long-range barrier for the Rb2 + system. The barrier is small but spoils the long-range behavior of the PEC. The effect is also found for other X2 + systems, such as X = Li, Na, and K. Calculations using a flexible framework for obtaining leading perturbative triples corrections derived using an analytic CC singles and doubles energy derivative formulation demonstrate that the origin of this problem lies in the use of T̂3 amplitudes obtained from approximate CC singles, doubles, and triples amplitude equations. It is shown that the unphysical barrier is related to a symmetry instability of the underlying Hartree-Fock mean-field solution, leading to orbitals representing two +0.5-fold charged ions in the limit of separated fragments. This, in turn, leads to a wrong 1/R asymptote of the interaction potential computed by perturbation-based CC approximations. Physically meaningful perturbative corrections in the long-range tail of the PEC may instead be obtained using symmetry-broken reference determinants.

On the Accuracy of Mean-Field Spin-Orbit Operators for 3d Transition-Metal Systems.

We present an extensive study of the performance of mean-field approximations to the spin-orbit operators on realistic molecular systems, as widely used in applications like single-molecule magnets, molecular quantum bits, and molecular spintronic devices. The test systems feature a 3d transition-metal center ion (V, Cr, Mn, Fe, Co, and Ni) in various coordinations and a multitude of energetically close-lying open-shell configurations that can couple via the spin-orbit operator. We performed complete active space spin-orbit configuration interaction calculations and compared the full two-electron Breit-Pauli spin-orbit operator to different approximations: the one-center approximation, the spin-orbit mean-field approach with electron densities from different state-averaging procedures, and the atomic mean-field integral approximation. We show that the mean-field approaches can lead to significant errors in the spin-orbital coupling matrix elements, which becomes particularly visible for the computed zero-field splittings. The one-center approximation, keeping all relevant two-electron terms, seems to be a significantly more accurate choice for the examples from our test set.

The Molpro quantum chemistry package.

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

The Quasi-Binary Acetonitriletriide Sr3 [C2 N]2.

The first quasi-binary acetonitriletriide Sr3 [C2 N]2 has been synthesised and characterised. The nearly colourless crystals were obtained from the reaction of Sr metal, graphite, and elemental N2 , generated by decomposition of Sr(N3 )2 , in a sealed Ni ampoule with the aid of an alkali metal flux. The structure of this compound was analysed via single-crystal X-ray diffraction and the identity of the [C2 N]3- anion was confirmed by Raman spectroscopy and further investigated by quantum-chemical methods. Computed interatomic distances within the [C2 N]3- anion strikingly match the obtained experimental data.

Flavylium Salts: A Blooming Core for Bioinspired Ionic Liquid Crystals.

Thermotropic ionic liquid crystals based on the flavylium scaffold have been synthesized and studied for their structure-properties relationship for the first time. The mesogens were probed by differential scanning calorimetry (DSC), polarizing optical microscopy (POM), and X-ray diffraction (XRD). Low numbers of alkoxy side chains resulted in smectic (SmA) and lamello-columnar (LamCol ) phases, whereas higher substituted flavylium salts showed Colro as well as ordered and disordered columnar (Colho , Colhd ) mesophases. Mesophase width ranged from 13 K to 220 K, giving access to room temperature liquid crystals. The optical properties of the synthesized compounds were probed towards absorption and emission properties. Strong absorption with maxima between 444 and 507 nm was observed, and some chromophores were highly emissive with quantum yields up to 99 %. Ultimately, mesogenic and dye properties were examined by temperature-dependent emissive experiments in the solid state.

The second-order approximate internally contracted multireference coupled-cluster singles and doubles method icMRCC2.

The second-order approximate internally contracted multireference coupled-cluster singles and doubles method icMRCC2 is defined and tested. The method is designed to bridge the gap between multireference perturbation theory and single-reference second-order approximate coupled-cluster theory (CC2). By including semi-internal double excitations into the zeroth-order expansion, the new method is able to reliably describe the coupling between excitations within the active space and the entire single-excitation spectrum. This helps, for instance, to provide a balanced treatment of valence and ionic states in polyenes (as explicitly demonstrated for cyclopentadiene) and to arrive at a more complete coverage of the excitation spectrum without the need to include diffuse orbitals into the active space. Good performance is also seen for notoriously difficult molecules such as C2 and CN. Furthermore, the multireference extension removes the main failures of single-reference CC2 theory, such as in the case of ozone.

The extended explicitly-correlated second-order approximate coupled-cluster singles and doubles ansatz suitable for response theory.

We report the extended explicitly correlated approximate coupled-cluster singles and doubles CC2(F12*)-XSP method suitable for response properties. Equations are derived using an automated approach and have subsequently been hand-coded into the computer program KOALA, in which for all two-electron integrals, density fitting is employed. Numerical results are presented for the lowest two vertical singlet excitation energies of a set of selected molecules. The results show that the CC2(F12*)-XSP method provides the correct basis-set limit with no bias to the ground state, and an excellent agreement with reference CC2 values using large basis sets is found. Using Dunning's aug-cc-pVTZ basis, the CC2(F12*)-XSP method yields excitation energies which are converged within 1 mEh to the basis-set limit for valence excitations.

Linear and quadratic internally contracted multireference coupled-cluster approximations.

Linear and quadratic approximations to the internally contracted multireference coupled-cluster (icMRCC) method are implemented and analyzed by using the linked and unlinked coupled-cluster formalisms. This includes methods based on perturbation theory as well as the coupled-electron pair approximation, CEPA(0). The similarities and differences between all the approximations serve to highlight and provoke discussion about methodological peculiarities of the icMRCC ansatz. When calculating potential energy curves (PECs), discontinuities are observed for the linear icMRCC energies. Using a diagrammatic representation, the terms that cause but also reduce these discontinuities are identified. For benchmarking test cases such as calculating PECs, singlet-triplet splittings, and barrier heights, the multireference CEPA(0) approximation performs well; however, it suffers from a lack of size consistency and so cannot represent a step forward to the goal of developing a computationally cheap and accurate icMRCC method.

On the distinguishable cluster approximation for triple excitations.

The distinguishable cluster approximation applied to coupled cluster doubles equations greatly improves absolute and relative energies. We apply the same approximation to the triples equations and demonstrate that it can also improve the results of the coupled cluster method with singles, doubles, and triples. The resulting method has a nominal computational scaling of O(N7) in the real-space representation, and is orbital invariant, size extensive, and exact for three electrons.

Perturbation Expansion of Internally Contracted Coupled-Cluster Theory up to Third Order.

The internally contracted multireference coupled-cluster (icMRCC) method is analyzed through third order in perturbation theory. Up to second order, the icMRCC perturbation expansion is equivalent to that of the standard Rayleigh-Schrödinger perturbation theory, which is based on a linear ansatz for the wave function, and the resulting theory is, depending on the employed zeroth-order Hamiltonian, equivalent to either second-order complete active space perturbation theory (CASPT2), N-electron valence perturbation theory (NEVPT2), or Fink's retention of the excitation degree perturbation theory (REPT2). At third order, the icMRCC perturbation expansion features additional terms in comparison to the Rayleigh-Schrödinger perturbation theory, but these are shown to be nearly negligibly small by both analytic arguments and numerical examples. Considering these systematic cancellations, however, may be important in future work on approximations to icMRCC theory. In addition, we provide an extensive set of tests of the second and third-order perturbation theories based on three different zeroth-order Hamiltonians, namely, the projected effective Fock operator as used for CASPT, the Dyall Hamiltonian as used for NEVPT, and the Fink Hamiltonian used for REPT. While the third-order variant of REPT often gives absolute energies that are rather close to values from higher level calculations, the results for relative energies and spectroscopic constants such as harmonic frequencies, give a less clear picture and a general conclusion about any best zeroth-order Hamiltonian does not emerge from our data. For small active spaces, REPT is rather prone to intruder state problems.

Do CCSD and approximate CCSD-F12 variants converge to the same basis set limits? The case of atomization energies.

While the title question is a clear "yes" from purely theoretical arguments, the case is less clear for practical calculations with finite (one-particle) basis sets. To shed further light on this issue, the convergence to the basis set limit of CCSD (coupled cluster theory with all single and double excitations) and of different approximate implementations of CCSD-F12 (explicitly correlated CCSD) has been investigated in detail for the W4-17 thermochemical benchmark. Near the CBS ([1-particle] complete basis set) limit, CCSD and CCSD(F12*) agree to within their respective uncertainties (about ±0.04 kcal/mol) due to residual basis set incompleteness error, but a nontrivial difference remains between CCSD-F12b and CCSD(F12*), which is roughly proportional to the degree of static correlation. The observed basis set convergence behavior results from the superposition of a rapidly converging, attractive, CCSD[F12]-CCSD-F12b difference (consisting mostly of third-order terms) and a more slowly converging, repulsive, fourth-order difference between CCSD(F12*) and CCSD[F12]. For accurate thermochemistry, we recommend CCSD(F12*) over CCSD-F12b if at all possible. There are some indications that the nZaPa family of basis sets exhibits somewhat smoother convergence than the correlation consistent family.

A model hamiltonian tuned toward high level ab initio calculations to describe the character of excitonic states in perylenebisimide aggregates.

On the example of an aggregate of two perylenebisimide (PBI) molecules the character of the lowest excited electronic states in terms of charge transfer (CT) and Frenkel exciton (FE) configurations is investigated as a function of the intermolecular arrangement. A minimal model Hamiltonian based on two FE and two CT configurations at the frontier-orbitals CIS (FOCIS) level is shown to represent a simple and comprehensible approach providing insight into the physical significance of the model Hamiltonian matrix elements. The recently introduced analysis and diabatization procedure (Liu et al., J. Chem. Phys. 2015, 143, 084106 ) method is used to extract the energies of the configurations and their interactions (the model Hamiltonian parameters) also from the accurate CC2 approach. An analysis in terms of diabatic energy profiles and their interactions shows that the FOCIS parameters give a qualitatively correct description of the adiabatic excited state energy profiles. Comparison with CC2 reveals, however, the presence of avoided crossings at FOCIS level, associated with a large character change (CT/FE) of the excited states as a function of the aggregate structure, which represents the major drawback of FOCIS results. We show that proper amendment of the FOCIS-derived parameters allows to accurately represent the potential energy surfaces and crossings of the excited dimer states as a function of the aggregate structure. © 2018 Wiley Periodicals, Inc.

First-order properties from internally contracted multireference coupled-cluster theory with particular focus on hyperfine coupling tensors.

Internally contracted multireference coupled-cluster (icMRCC) theory is extended to the computation of first-order properties (expectation values). We use the previously defined Lagrange formulation of the energy functional to derive the required equations for the Lagrange multipliers and arrive at an expression for first-order properties according to the generalized Hellmann-Feynman theorem, analogous to single-reference coupled-cluster theory. The present formulation does not include orbital relaxation, but in line with previous experience in coupled-cluster theory, the single-excitation cluster operator can recover a significant portion of orbital relaxation. Further aspects of the theory that arise from the internal contraction approach are discussed. Using automated derivation techniques, we have implemented a pilot code for icMRCCSD and icMRCCSDT for testing the method numerically. We find good agreement with full configuration interaction for several properties of boron monohydride and dipole moment curves of hydrogen fluoride and chromium hydride. A particular focus is given to spin-dependent properties: The hyperfine coupling tensors of Σ and Π radicals have been computed and compared to experiment and previous computations. We discuss the problem of describing spin polarization with properly spin-adapted wavefunctions, which requires either including pseudo-triple excitations or employing sufficiently flexible reference functions.

Embedded Multireference Coupled Cluster Theory.

Internally contracted multireference coupled cluster (icMRCC) theory is embedded within multireference perturbation theory (MRPT) to calculate energy differences in large strongly correlated systems. The embedding scheme is based on partitioning the orbital spaces of a complete active space self-consistent field (CASSCF) wave function, with a truncated virtual space constructed by transforming selected projected atomic orbitals (PAOs). MRPT is applied to the environment using a subtractive embedding approach that also allows for multilayer embedding. Benchmark calculations are presented for biradical bond dissociation, spin splitting in a heterocyclic carbene and hydrated Fe(II), and for the super-exchange coupling constant in solid nickel oxide. The method is further applied to two large transition metal complexes with a triple-ζ basis set: an iron complex with 175 atoms and 2939 basis functions, and a nickel complex with 231 atoms, and 4175 basis functions.