Accuracy of relativistic energy-consistent pseudopotentials for superheavy elements 111–118: Molecular calibration calculations

Error of relativistic effective core potentials for closed-shell diatomic molecules of p-block heavy and superheavy elements in DFT and TDDFT calculations

Pseudopotentials (PP) are extensively used in electronic structure calculations, particularly for molecules containing heavy elements. Parameters in PPs are mainly determined from ab initio results, and errors of such PPs in density functional theory (DFT) calculations have been studied previously. However, PP errors on results with spin-orbit coupling and those in time-dependent DFT (TDDFT) calculations have not been reported previously. In this work, we investigate the error of the small-core energy-consistent Stuttgart/Koln pseudopotentials in DFT and TDDFT calculations with and without spin-orbit coupling. Ground state bond lengths, harmonic frequencies, dissociation energies, and vertical excitation energies for a series of closed-shell diatomic heavy and superheavy p-block molecules are calculated using several popular exchange-correlation functionals. PP errors are estimated by comparing with results using the all-electron Dirac-Coulomb (-Gaunt) Hamiltonian. Our results show that the difference between ground state properties and most excitation energies in scalar-relativistic calculations with the PP and those of all-electron calculations is quite small. This difference becomes somewhat larger when spin-orbit coupling (SOC) is present, especially for properties that are affected by SOC to some extent. In addition, the errors of the PPs are insensitive to the employed exchange-correlation functionals in most cases. Our results indicate that reasonable DFT and TDDFT results can be obtained using the small-core energy-consistent Stuttgart/Koln pseudopotentials for heavy and super-heavy p-block molecules.

CRYSTAL23: A Program for Computational Solid State Physics and Chemistry

The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.

Stable copernicium hexafluoride (CnF6) with an oxidation state of VI

As the heaviest group 12 element known currently, copernicium (Cn) often presents the oxidation states of I+, II+, and rarely IV+ as in its homologue mercury. In this work we systematically studied the stability of some oxides, fluorides, and oxyfluorides of Cn by two-component relativistic calculations and found that the CnF₆ molecule with an oxidation state of VI+ has an extraordinary stability. CnF₆ may decompose into CnF₄ by conquering an energy barrier of about 34 kcal mol^-1 without markedly releasing heat. Our results indicate that CnF₆ may exist under some special conditions.

The cohesive energy of superheavy element copernicium determined from accurate relativistic coupled-cluster theory

The cohesive energy of bulk copernicium is accurately determined using the incremental method within a relativistic coupled-cluster approach. For the lowest energy structure of hexagonal close-packed (hcp) symmetry, we obtain a cohesive energy of -36.3 kJ mol^-1 (inclusion of uncertainties leads to a lower bound of -39.6 kJ mol^-1), in excellent agreement with the experimentally estimated sublimation enthalpy of -38 kJ mol^-1 [R. Eichler et al., Angew. Chem. Int. Ed., 2008, 47, 3262]. At the coupled-cluster singles, doubles and perturbative triples level of theory, we find that the hcp structure is energetically quasi-degenerate with both face-centred and body-centred cubic structures. These results provide a basis for testing various density-functionals, of which the PBEsol functional yields a cohesive energy of -34.1 kJ mol^-1 in good agreement with our coupled-cluster value.

Electron and Nucleon Localization Functions of Oganesson: Approaching the Thomas-Fermi Limit

Fermion localization functions are used to discuss electronic and nucleonic shell structure effects in the superheavy element oganesson, the heaviest element discovered to date. Spin-orbit splitting in the 7p electronic shell becomes so large (∼10  eV) that Og is expected to show uniform-gas-like behavior in the valence region with a rather large dipole polarizability compared to the lighter rare gas elements. The nucleon localization in Og is also predicted to undergo a transition to the Thomas-Fermi gas behavior in the valence region. This effect, particularly strong for neutrons, is due to the high density of single-particle orbitals.