New correlation consistent basis sets optimized using the all-electron third-order Douglas-Kroll-Hess (DKH3) scalar relativistic Hamiltonian are reported for the actinide elements Ac and Np through Lr. These complete the series of sets reported previously for Th-U [K. A. Peterson, J. Chem. Phys. 142, 074105 (2015); M. Vasiliu et al., J. Phys. Chem. A 119, 11422 (2015)]. The new sets range in size from double- to quadruple-zeta and encompass both those optimized for valence (6s6p5f7s6d) and outer-core electron correlations (valence + 5s5p5d). The final sets have been contracted for both the DKH3 and eXact 2-component (X2C) Hamiltonians, yielding cc-pVnZ-DK3/cc-pVnZ-X2C sets for valence correlation and cc-pwCVnZ-DK3/cc-pwCVnZ-X2C sets for outer-core correlation (n = D, T, Q in each case). In order to test the effectiveness of the new basis sets, both atomic and molecular benchmark calculations have been carried out. In the first case, the first three atomic ionization potentials (IPs) of all the actinide elements Ac-Lr have been calculated using the Feller-Peterson-Dixon (FPD) composite approach, primarily with the multireference configuration interaction (MRCI) method. Excellent convergence towards the respective complete basis set (CBS) limits is achieved with the new sets, leading to good agreement with experiment, where these exist, after accurately accounting for spin-orbit effects using the 4-component Dirac-Hartree-Fock method. For a molecular test, the IP and atomization energy (AE) of PuO have been calculated also using the FPD method but using a coupled cluster approach with spin-orbit coupling accounted for using the 4-component MRCI. The present calculations yield an IP for PuO of 159.8 kcal/mol, which is in excellent agreement with the experimental electron transfer bracketing value of 162 ± 3 kcal/mol. Likewise, the calculated 0 K AE of 305.6 kcal/mol is in very good agreement with the currently accepted experimental value of 303.1 ± 5 kcal/mol. The ground state of PuO is predicted to be the Σ0g+5 state.

The developments of the open-source chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes can address, while showing that is an attractive platform for state-of-the-art atomistic computer simulations.

Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab initio wavefunction calculations at the multireference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin–orbit coupling was treated via RAS state interaction (SO–RASSI). The method was tested by calculations of electron–nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO–RASSI implementation for hyperfine coupling that relied on nonrelativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538–549] and therefore had limited applicability. Results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.

A comprehensive ab initio study of periodic actinide−ligand bonding trends for trivalent actinides is performed. Relativistic density functional theory (DFT) and complete activespace (CAS) self-consistent field wavefunction calculations are used to dissect the chemical bonding in the [AnCl 6 ] 3− , [An-(CN) 6 ] 3− , [An(NCS) 6 ] 3− , [An(S 2 PMe 2 ) 3 ], [An(DPA) 3 ] 3− , and [An(HOPO)] − series of actinide (An = U−Es) complexes. Except for some differences for the early actinide complexes with DPA, bond orders and excess 5f-shell populations from donation bonding show qualitatively similar trends in 5f n active-space CAS vs DFT calculations. The influence of spin−orbit coupling on donation bonding is small for the tested systems. Along the actinide series, chemically soft vs chemically harder ligands exhibit clear differences in bonding trends. There are pronounced changes in the 5f populations when moving from Pu to Am or Cm, which correlate with previously noted "breaks" in chemical trends. Bonding involving 5f becomes very weak beyond Cm/Bk. We propose that Cm(III) is a borderline case among the trivalent actinides that can be meaningfully considered to be involved in ground-state 5f covalent bonding.