Electron superhalogens as positronium superhalogens.

Positronium (Ps) exhibits the ability to form energetically stable complexes with atoms and molecules before annihilation occurs. In particular, F, a halogen, shows the highest reported positronium binding energy (2.95 eV) in the periodic table. Superhalogens are defined as molecules with electron affinities exceeding that of Cl (3.61 eV), the atom with the highest electron affinity. Building upon the concept of superhalogens, we can define Ps-superhalogens as molecules with Ps binding energies surpassing that of F. This study explores structural and energetic aspects of positronium and positron binding to neutral and anionic superhalogen molecules of the MXk+1 family (M = Li, Na, Be, Mg, B, Al, Si, P; X = F, Cl, Br), respectively and where k represents the highest formal valence of M. We perform multicomponent MP2 calculations for positron systems, which reveal how positron affinities vary with the type and number of halogen atoms present. The analysis of the results emphasizes the predominant role of electrostatic interactions in determining the positron affinity, with negligible effects of electronic and geometric relaxation upon positron attachment. We predict the energetic stability of 22 of the 24 PsMXk+1 complexes with respect to the chemically relevant dissociation channels: e+ emission, Ps emission and M-X bond breaking. Our findings reveal six MFk+1 systems that qualify as Ps-superhalogens, showing a positronium binding energy exceeding 2.95 eV. Of these, AlF4 stands out by setting a new record for the highest positronium binding energy among neutral molecules, reaching 4.36 eV.

Four-Dimensional Scaling of Dipole Polarizability in Quantum Systems.

Polarizability is a key response property of physical and chemical systems, which has an impact on intermolecular interactions, spectroscopic observables, and vacuum polarization. The calculation of polarizability for quantum systems involves an infinite sum over all excited (bound and continuum) states, concealing the physical interpretation of polarization mechanisms and complicating the derivation of efficient response models. Approximate expressions for the dipole polarizability, α, rely on different scaling laws α∝R^{3}, R^{4}, or R^{7}, for various definitions of the system radius R. Here, we consider a range of single-particle quantum systems of varying spatial dimensionality and having qualitatively different spectra, demonstrating that their polarizability follows a universal four-dimensional scaling law α=C(4µq^{2}/ℏ^{2})L^{4}, where µ and q are the (effective) particle mass and charge, C is a dimensionless excitation-energy ratio, and the characteristic length L is defined via the L^{2} norm of the position operator. This unified formula is also applicable to many-particle systems, as shown by accurately predicting the dipole polarizability of 36 atoms, 1641 small organic molecules, and Bloch electrons in periodic systems.

Binding Matter with Antimatter: The Covalent Positron Bond.

We report sufficient theoretical evidence of the energy stability of the e+ ⋅H22- molecule, formed by two H- anions and one positron. Analysis of the electronic and positronic densities of the latter compound undoubtedly points out the formation of a positronic covalent bond between the otherwise repelling hydride anions. The lower limit for the bonding energy of the e+ ⋅H22- molecule is 74 kJ mol-1 (0.77 eV), accounting for the zero-point vibrational correction. The formation of a non electronic covalent bond is fundamentally distinct from positron attachment to stable molecules, as the latter process is characterized by a positron affinity, analogous to the electron affinity.

Fast and accurate prediction of proton affinities: revisiting the extended Koopmans' theorem for protons.

In this work we propose schemes based on the extended Koopmans' theorem for quantum nuclei (eKT), in the framework of the any particle molecular orbital approach (APMO/KT), for the quantitative prediction of gas phase proton affinities (PAs). The performance of these schemes has been tested on a set of 300 organic molecules containing diverse functional groups. The APMO/KT scheme scaled by functional group (APMO/KT-SC-FG) displays an overall mean absolute error of 1.1 kcal mol-1 with respect to experimental data. Its performance in PA calculations is similar to that of post-Hartree-Fock composite methods or that of the APMO second order proton propagator (APMO/PP2) approach. The APMO/KT-SC-FG scheme is also employed to predict PAs of polyfunctional molecules such as the Nerve Agent VX and the 20 common α-amino acids, finding excellent agreement with available theoretical and/or experimental data. The accuracy of the predictions demonstrates that the APMO/KT-SC-FG scheme is a low-cost alternative to adiabatic methods for the calculation of accurate PAs. One of the most appealing features of the APMO/KT-SC-FG scheme, is that PAs can be derived from one single-point APMO Hartree-Fock calculation.

On the physical interpretation of the nuclear molecular orbital energy.

Recently, several groups have extended and implemented molecular orbital (MO) schemes to simultaneously obtain wave functions for electrons and selected nuclei. Many of these schemes employ an extended Hartree-Fock approach as a first step to find approximate electron-nuclear wave functions and energies. Numerous studies conducted with these extended MO methodologies have explored various effects of quantum nuclei on physical and chemical properties. However, to the best of our knowledge no physical interpretation has been assigned to the nuclear molecular orbital energy (NMOE) resulting after solving extended Hartree-Fock equations. This study confirms that the NMOE is directly related to the molecular electrostatic potential at the position of the nucleus.

Calculation of positron binding energies using the generalized any particle propagator theory.

We recently extended the electron propagator theory to any type of quantum species based in the framework of the Any-Particle Molecular Orbital (APMO) approach [J. Romero, E. Posada, R. Flores-Moreno, and A. Reyes, J. Chem. Phys. 137, 074105 (2012)]. The generalized any particle molecular orbital propagator theory (APMO/PT) was implemented in its quasiparticle second order version in the LOWDIN code and was applied to calculate nuclear quantum effects in electron binding energies and proton binding energies in molecular systems [M. Díaz-Tinoco, J. Romero, J. V. Ortiz, A. Reyes, and R. Flores-Moreno, J. Chem. Phys. 138, 194108 (2013)]. In this work, we present the derivation of third order quasiparticle APMO/PT methods and we apply them to calculate positron binding energies (PBEs) of atoms and molecules. We calculated the PBEs of anions and some diatomic molecules using the second order, third order, and renormalized third order quasiparticle APMO/PT approaches and compared our results with those previously calculated employing configuration interaction (CI), explicitly correlated and quantum Montecarlo methodologies. We found that renormalized APMO/PT methods can achieve accuracies of ~0.35 eV for anionic systems, compared to Full-CI results, and provide a quantitative description of positron binding to anionic and highly polar species. Third order APMO/PT approaches display considerable potential to study positron binding to large molecules because of the fifth power scaling with respect to the number of basis sets. In this regard, we present additional PBE calculations of some small polar organic molecules, amino acids and DNA nucleobases. We complement our numerical assessment with formal and numerical analyses of the treatment of electron-positron correlation within the quasiparticle propagator approach.

The three-center two-positron bond.

Computational studies have shown that one or more positrons can stabilize two repelling atomic anions through the formation of two-center positronic bonds. In the present work, we study the energetic stability of a system containing two positrons and three hydride anions, namely 2e+[H3 3-]. To this aim, we performed a preliminary scan of the potential energy surface of the system with both electrons and positrons in a spin singlet state, with a multi-component MP2 method, that was further refined with variational and diffusion Monte Carlo calculations, and confirmed an equilibrium geometry with D 3h symmetry. The local stability of 2e+[H3 3-] is demonstrated by analyzing the vertical detachment and adiabatic energy dissociation channels. Bonding properties of the positronic compound, such as the equilibrium interatomic distances, force constants, dissociation energies, and bonding densities are compared with those of the purely electronic H3 + and Li3 + systems. Through this analysis, we find compelling similarities between the 2e+[H3 3-] compound and the trilithium cation. Our results strongly point out the formation of a non-electronic three-center two-positron bond, analogous to the well-known three-center two-electron counterparts, which is fundamentally distinct from the two-center two-positron bond [D. Bressanini, J. Chem. Phys., 2021, 155, 054306], thus extending the concept of positron bonded molecules.

Covalent bonds in positron dihalides.

We report a computational study on homo- and heteronuclear e+[X-Y-] compounds formed by two halide anions (X-, Y- = F-, Cl-, Br-) and one positron. Our results indicate the formation of energetically stable positronic molecules in all cases. Analysis of the electron and positron densities points out that the formation of positron covalent bonds underlies the stabilization of the otherwise repelling dihalides, revealing that positronic bonding can reach far beyond the previously addressed e+[H-H-] molecule [J. Charry, M. T. do N. Varella and A. Reyes, Angew. Chem. Int. Ed., 2018, 57, 8859-8864.]. To a significant extent, the properties of the positron dihalides are similar to those of the purely electronic analogs, e-[A+B+], molecular cations with isoelectronic atomic cores (A+, B+ = Na+, K+, Rb+) bound by one electron. The positron bonds in the e+[X-Y-] complexes are however stronger than those in the isoelectronic e-[A+B+] counterparts, as the former have shorter bond lengths and higher bond energies. While an energy decomposition analysis points out that both electronic and positronic bonds essentially arise from electrostatic interactions, the more stable positron bonds are partly due to the higher polarizabilities of the dihalide anions, and partly to more significant contributions from correlation and relaxation effects.