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.

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.

Vibrational modes ascribed to the stretching of X-H bonds from donor monomers (HXdonor) in complexes presenting hydrogen bonds (HF···HF, HCl···HCl, HCN···HCN, HNC···HNC, HCN···HF, HF···HCl and H2O···HF) exhibit large (4 to 7 times) infrared intensity increments during complexation according to CCSD/cc-pVQZ-mod calculations. These intensity increases are explained by the charge-charge flux-dipole flux (CCFDF) model based on multipoles from the Quantum Theory of Atoms in Molecules (QTAIM) as resulting from a reinforcing interaction between two contributions to the dipole moment derivatives with respect to the vibrational displacements: charge and charge flux. As such, variations that occur in their intensity cross terms in hydrogen bond formation correlate nicely with the intensity enhancements. These stretching modes of HXdonor bonds can be approximately modeled by sole displacement of the positively charged hydrogens towards the acceptor terminal atom with concomitant electronic charge transfers in the opposite direction that are larger than those occurring for the H atom displacements of their isolated donor molecules. This analysis indicates that the charge-charge flux interaction reinforcement on H-bond complexation is associated with variations of atomic charge fluxes in both parent molecules and small electronic charge transfers between them. The QTAIM/CCFDF model also indicates that atomic dipole flux contributions do not play a significant role in these intensity enhancements.

The Quantum Theory of Atoms In Molecules/Charge-Charge Flux-Dipole Flux (QTAIM/CCFDF) model has been used to investigate the electronic structure variations associated with intensity changes on dimerization for the vibrations of the water and hydrogen fluoride dimers as well as in the water-hydrogen fluoride complex. QCISD/cc-pVTZ wave functions applied in the QTAIM/CCFDF model accurately provide the fundamental band intensities of water and its dimer predicting symmetric and antisymmetric stretching intensity increases for the donor unit of 159 and 47 km mol(-1) on H-bond formation compared with the experimental values of 141 and 53 km mol(-1). The symmetric stretching of the proton donor water in the dimer has intensity contributions parallel and perpendicular to its C2v axis. The largest calculated increase of 107 km mol(-1) is perpendicular to this axis and owes to equilibrium atomic charge displacements on vibration. Charge flux decreases occurring parallel and perpendicular to this axis result in 42 and 40 km mol(-1) total intensity increases for the symmetric and antisymmetric stretches, respectively. These decreases in charge flux result in intensity enhancements because of the interaction contributions to the intensities between charge flux and the other quantities. Even though dipole flux contributions are much smaller than the charge and charge flux ones in both monomer and dimer water they are important for calculating the total intensity values for their stretching vibrations since the charge-charge flux interaction term cancels the charge and charge flux contributions. The QTAIM/CCFDF hydrogen-bonded stretching intensity strengthening of 321 km mol(-1) on HF dimerization and 592 km mol(-1) on HF:H2O complexation can essentially be explained by charge, charge flux and their interaction cross term. Atomic contributions to the intensities are also calculated. The bridge hydrogen atomic contributions alone explain 145, 237, and 574 km mol(-1) of the H-bond stretching intensity enhancements for the water and HF dimers and their heterodimer compared with total increments of 149, 321, and 592 km mol(-1), respectively.

No abstract

No abstract

No abstract

Some atomic multipoles (charges, dipoles and quadrupoles) from the Quantum Theory of Atoms in Molecules (QTAIM) and CHELPG charges are used to investigate interactions between a proton and a molecule (F2, Cl2, BF, AlF, BeO, MgO, LiH, H2CO, NH3, PH3, BF3, and CO2). Calculations were done at the B3LYP/6‐311G(3d,3p) level. The main aspect of this work is the investigation of polarization effects over electrostatic potentials and atomic multipoles along a medium to long range of interaction distances. Large electronic charge fluxes and polarization changes are induced by a proton mainly when this positive particle approaches the least electronegative atom of diatomic heteronuclear molecules. The search for simple equations to describe polarization on electrostatic potentials from QTAIM quantities resulted in linear relations with r − 4 (r is the interaction distance) for many cases. Moreover, the contribution from atomic dipoles to these potentials is usually the most affected contribution by polarization what reinforces the need for these dipoles to a minimal description of purely electrostatic interactions. Finally, CHELPG charges provide a description of polarization effects on electrostatic potentials that is in disagreement with physical arguments for certain of these molecules. © 2012 Wiley Periodicals, Inc.