First, exploratory calculations were performed to investigate the presence of non-nuclear maxima (NNMs) in ground-state electron densities of homonuclear diatomic molecules from hydrogen up to calcium at their equilibrium geometries. In a second stage, only for the cases in which these features were previously detected, a rigorous analysis was carried out by several combinations of theoretical methods and basis sets in order to ensure that they are not only calculation artifacts. Our best results support that Li2, B2, C2, and P2 are molecules that possess true NNMs. A NNM was found in values obtained from the largest basis sets for Na2, but it disappeared at the experimental geometry because optimized bond lengths are significantly inaccurate for this case (deviations of 0.10 Å). Two of these maxima are also observed in Si2 with CCSD and large basis sets, but they are no longer detected as core-valence correlation or multiconfigurational wave functions are taken into account. Therefore, the NNMs in Si2 can be considered unphysical features due to an incomplete treatment of electron correlation. Finally, we show that a NNM is encountered in LiNa, representing the first discovery of such electron density maxima in a heteronuclear diatomic system at its equilibrium geometry, to our knowledge. Some results for LiNa, found in variations in internuclear distances, suggest that molecular electric moments, such as dipole and quadrupole, are sensitive to the presence of NNMs.

The nuclear electric quadrupole moment (NQM) of the (63)Cu nucleus was determined from an indirect approach by combining accurate experimental nuclear quadrupole coupling constants (NQCCs) with relativistic Dirac-Coulomb coupled cluster calculations of the electric field gradient (EFG). The data obtained at the highest level of calculation, DC-CCSD-T, from 14 linear molecules containing the copper atom give rise to an indicated NQM of -198(10) mbarn. Such result slightly deviates from the previously accepted standard value given by the muonic method, -220(15) mbarn, although the error bars are superimposed.

Nonnuclear attractors (NNAs) are observed in the electron density of a variety of systems, but the factors governing their appearance and their contribution to the system's properties remain a mystery. The NNA occurring in homo- and heteronuclear diatomics of main group elements with atomic numbers up to Z = 38 is investigated computationally (at the UCCSD/cc-pVQZ level of theory) by varying internuclear separations. This was done to determine the NNA occurrence window along with the evolution of the respective pseudoatomic basin properties. Two distinct categories of NNAs were detected in the data analyzed by means of catastrophe theory. Type "a" implies electronic charge transfer between atoms mediated by a pseudoatom. Type "b" shows an initial relocation of some electronic charge to a pseudoatom, which posteriorly returns to the same atom that donated this charge in the first place. A small difference of polarizability between the atoms that compose these heteronuclear diatomics seems to favor NNA formation. We also show that the NNA arising tends to result in some perceptible effects on molecular dipole and/or quadrupole moment curves against internuclear distance. Finally, successive cationic ionization results in the fast disappearance of the NNA in Li2 indicating that its formation is mainly governed by the field generated by the quantum mechanical electronic density and only depends parametrically on the bare nuclear field/potential at a given molecular geometry.

This study reports a new relativistic prolapse-free Gaussian basis set series of quadruple-ζ quality, RPF-4Z, and an augmented version that includes extra diffuse functions, aug-RPF-4Z, for all the s- and p-block elements. The relativistic adapted Gaussian basis sets (RAGBSs), which are free of variational prolapse, were used as the starting primitive sets. Exponents of correlating/polarization functions were taken from a polynomial version of the generator coordinate Dirac-Fock (p-GCDF) method, in which the previously optimized RAGBS parameters are applied. By using such procedure we aimed to reduce the computational demand of these sets in comparison with fully optimized ones. The effect of these basis set increments on the correlation energy was evaluated by atomic multireference configuration interaction calculations with single and double excitations out of the valence shell. Finally, atomic and molecular calculations of fundamental properties (bond lengths, vibrational frequencies, dipole moments and electron affinities) corroborate the quadruple-ζ quality of these new sets that are also about half-time-consuming than the correspondent Dyall's v4z sets. The read-to use format of these (aug-)RPF-4Fz sets are available as Supporting Information files and can also be found at http://basis-sets.iqsc.usp.br/ .

No abstract