Are "GAPT Charges" Really Just Charges?

Symmetry-Constrained Properties Behave Differently for 2D or 3D Structures under the Same Point Group

In chemistry and physics, two molecules belonging to the same point group are expected to behave similarly regarding various properties, following their character tables. Here, we show that the derivative of the dipole moment with respect to the normal coordinate of vibration might show different symmetry constraints if the molecule is planar, even if these molecules belong to the same point group. Examples of pairs of molecules featuring these conditions are presented. These findings open a new path toward a much deeper understanding of how 2D materials behave so differently compared to 3D materials featuring the very same atoms and arrangements (graphene and graphite, for example); chemists and physicists dealing with 2D materials could benefit from looking more deeply into pure mathematical relations that might be governing 2D systems in a different way when compared to 3D systems. The aid from mathematicians is welcomed.

Molecular Pseudorotation in Phthalocyanines as a Tool for Magnetic Field Control at the Nanoscale

Metal phthalocyanines, a highly versatile class of aromatic, planar, macrocyclic molecules with a chelated central metal ion, are topical objects of ongoing research and particularly interesting due to their magnetic properties. However, while the current focus lies almost exclusively on spin-Zeeman-related effects, the high symmetry of the molecule and its circular shape suggests the exploitation of light-induced excitation of 2-fold degenerate vibrational states in order to generate, switch, and manipulate magnetic fields at the nanoscale. The underlying mechanism is a molecular pseudorotation that can be triggered by infrared pulses and gives rise to a quantized, small, but controllable magnetic dipole moment. We investigate the optical stimulation of vibrationally induced molecular magnetism and estimate changes in the magnetic shielding constants for confirmation by future experiments.

On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment

Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges q on 1-particle basis set and (for WFT methods) n -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence.

Assessment of an Electrostatic Energy-Based Charge Model for Modeling the Electrostatic Interactions in Water Solvent

The protein force field based on the restrained electrostatic potential (RESP) charges has limitations in accurately describing hydrogen bonding interactions in proteins. To address this issue, we propose an alternative approach called the electrostatic energy-based charges (EEC) model, which shows improved performance in describing electrostatic interactions (EIs) of hydrogen bonds in proteins. In this study, we further investigate the performance of the EEC model in modeling EIs in water solvent. Our findings demonstrate that the fixed EEC model can effectively reproduce the quantum mechanics/molecular mechanics (QM/MM)-calculated EIs between a water molecule and various water solvent environments. However, to achieve the same level of computational accuracy, the electrostatic potential (ESP) charge model needs to fluctuate according to the electrostatic environment. Our analysis indicates that the requirement for charge adjustments depends on the specific mathematical and physical representation of EIs as a function of the environment for deriving charges. By comparing with widely used empirical water models calibrated to reproduce experimental properties, we confirm that the performance of the EEC model in reproducing QM/MM EIs is similar to that of general purpose TIP4P-like water models such as TIP4P-Ew and TIP4P/2005. When comparing the computed 10,000 distinct EI values within the range of -40 to 0 kcal/mol with the QM/MM results calculated at the MP2/aug-cc-pVQZ/TIP3P level, we noticed that the mean unsigned error (MUE) for the EEC model is merely 0.487 kcal/mol, which is remarkably similar to the MUE values of the TIP4P-Ew (0.63 kcal/mol) and TIP4P/2005 (0.579 kcal/mol) models. However, both the RESP method and the TIP3P model exhibit a tendency to overestimate the EIs, as evidenced by their higher MUE values of 1.761 and 1.293 kcal/mol, respectively. EEC-based molecular dynamics simulations have demonstrated that, when combined with appropriate van der Waals parameters, the EEC model can closely reproduce oxygen-oxygen radial distribution function and density of water, showing a remarkable similarity to the well-established TIP4P-like empirical water models. Our results demonstrate that the EEC model has the potential to build force fields with comparable accuracy to more sophisticated empirical TIP4P-like water models.

On the quest for a molecular vibration whose absolute intensity is described solely by fluctuations of atomic dipoles: Can we find it?

The conditions for a molecular vibration to be active in the infrared spectrum while its absolute intensity is described solely by fluctuations of atomic dipoles are presented. A quick recipe of features of such a system is presented, to guide the seek for it. If such a system can be found, it will then be proved that population analyses based solely on a point-charge approximation are unrealistic as they cannot properly describe this intensity, a real, measurable, unambiguous property.

Berry Population Analysis: Atomic Charges from the Berry Curvature in a Magnetic Field

The Berry curvature is essential in Born-Oppenheimer molecular dynamics, describing the screening of the nuclei by the electrons in a magnetic field. Parts of the Berry curvature can be understood as the external magnetic field multiplied by an effective charge so that the resulting Berry force behaves like a Lorentz force during the simulations. Here, we investigate whether these effective charges can provide insight into the electronic structure of a given molecule or, in other words, whether we can perform a population analysis based on the Berry curvature. To develop our approach, we first rewrite the Berry curvature in terms of charges that partially capture the effective charges and their dependencies on the nuclear velocities. With these Berry charges and charge fluctuations, we then construct our population analysis yielding atomic charges and overlap populations. Calculations at the Hartree-Fock level reveal that the atomic charges are similar to those obtained from atomic polar tensors. However, since we additionally obtain an estimate for the fluctuations of the charges and a partitioning of the atomic charges into contributions from all atoms, we conclude that the Berry population analysis is a useful alternative tool to analyze the electronic structures of molecules.

How deeply should we analyze non-covalent interactions?

CONTEXT

Just how much effort and detail should we invest in analyzing interactions of the order of 5 kcal mol^-1? This comment attempts to provide a conciliatory overview of what is often a contentious field and to pose some questions that I hope will eventually lead at least to some consensus.

METHODS

This is an opinion article without calculations or data.

Energy Electronegativity and Chemical Bonding

Historical development of the concept of electronegativity (EN) and its significance and prospects for physical and structural chemistry are discussed. The current cutting-edge results are reviewed: new methods of determining the ENs of atoms in solid metals and of bond polarities and effective atomic charges in molecules and crystals. The ENs of nanosized elements are calculated for the first time, enabling us to understand their unusual reactivity, particularly the fixation of N₂ by nanodiamond. Bond polarities in fluorides are also determined for the first time, taking into account the peculiarities of the fluorine atom's electronic structure and its electron affinity.

Infrared intensities of [Formula: see text]: a true challenge for DFT methods

Absolute infrared intensities of [Formula: see text] were evaluated with a great variety of DFT and ab initio methods and basis sets. It is shown that the intensities calculated by different levels of theory may not agree with each other even in the qualitative (weak/strong) sense. Geometrical parameters, as well as net atomic charges evaluated from multiple partition schemes, did not vary as much as the intensities and thus cannot explain the tremendous differences found for the latter. As there are no experimental estimates for the intensities to guide the theoretical evaluation, deciding on the best level of theory is reduced to comparisons between the different DFT methods and QCISD or CCSD, believed to be the best theoretical estimates in the set. The differences found among the various DFT methods suggest the development of new methods, instead of converging to a focal point, is rather diverging.

Unavoidable failure of point charge descriptions of electronic density changes for out-of-plane distortions

Population analyses based on point charge approximations accurately estimating the equilibrium dipole moment will systematically fail when predicting infrared intensities of out-of-plane vibrations of planar molecules, whereas models based on both charges and dipoles will always succeed. It is not a matter of how the model is devised but rather how many degrees of freedom are available for the calculation. Population analyses based on point charges are very limited in terms of the amount of meaningful chemical information they provide, whereas models employing both atomic charges and atomic dipoles should be preferred for molecular distortions. A good model should be able to correctly describe not only static, equilibrium structures but also distorted geometries in order to correctly assess information from vibrating molecules. The limitations of point charge models also hold for distortions much larger than those encountered vibrationally.

Electrostatics Explains the Reverse Lewis Acidity of BH3 and Boron Trihalides: Infrared Intensities and a Relative Energy Gradient (REG) Analysis of IQA Energies

The reaction path for the formation of BX₃–NH₃ (X = H, F, Cl, Br) complexes was divided into two processes: (i) rehybridization of the acid while adopting a pyramidal geometry, and (ii) the complex formation from the pyramidal geometries of the acid and base. The interacting quantum atom (IQA) method was used to investigate the Lewis acidity trend of these compounds. This topological analysis suggests that the boron–halogen bond exhibits a considerable degree of ionicity. A relative energy gradient (REG) analysis on IQA energies indicates that the acid–base complex formation is highly dependent on electrostatic energy. With increasing halogen electronegativity, a higher degree of ionicity of the B–X is observed, causing an increase in the absolute value of X and B charges. This increases not only the attractive electrostatic energy between the acid and base but also enhances the repulsive energy. The latter is the main factor behind the acidity trend exhibited by trihalides. Changes in geometry are relevant only for complexes where BH₃ acts as an acid, where lower steric hindrance facilitates the adoption of the pyramidal geometry observed in the complex. The CCTDP analysis shows that infrared intensities of BX₃–NH₃ are determined mostly by the atomic charges and not by the charge transfer or polarization. The opposite is observed in covalent analogues.