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

Kernel-Based Minimal Distributed Charges: A Conformationally Dependent ESP-Model for Molecular Simulations

A kernel-based method (kernelized minimal distributed charge model (kMDCM)) to represent the molecular electrostatic potential (ESP) in terms of off-center point charges is introduced. The positions of the charges adapt to the molecular geometry and allow the description of intramolecular charge flow. Using Gaussian kernels and atom-atom distances as the features, the ESPs for water and methanol are shown to improve by at least a factor of 2 compared with point charge models fit to an ensemble of structures. The conformationally fluctuating molecular dipole moment of water is reproduced almost twice as accurately using kMDCM compared with static PCs, despite not fitting to the dipole directly. The role of hyperparameters in the kernelization is investigated and their implication on model performance and simulation stability is discussed. Combining kMDCM for the electrostatics and reproducing kernels for the bonded terms allows energy-conserving simulations of 2000 water molecules with periodic boundary conditions on the nanosecond time scale. These MD simulations sample geometries outside the training set but remain stable, which demonstrates the robustness of the model and its implementation.

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.

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.

Electronic Distribution of SN2 IRC and TS Structures: Infrared Intensities of Imaginary Frequencies

A novel IRC-TS-CCTDP method to investigate transition states (TS) is proposed in which changes in the molecular geometry follow atomic displacements corresponding to the imaginary frequency normal coordinate. Electronic charge structure changes can be analyzed using the charge-charge-transfer-dipolar polarization (CCTDP) model. An application is presented for the gas-phase S_N2 reaction transition state structures for nine NuCX₃LG^- systems, with Nu and LG = H, F, Cl and X = H, F. Using quantum theory of atoms in molecules (QTAIM) at the QCISD/aug-cc-pVTZ level, atomic charges and atomic dipoles were obtained and applied to calculate the CCTDP contributions to their imaginary normal mode intensities. The results show that the imaginary bands are exceptionally strong, ranging from 1217 to 16 086 km·mol^-1, much higher than the stretching intensities found in the methyl halides (that are all less than 100 km·mol^-1). For all systems, the CT contributions are responsible for 63% of the total dipole moment derivatives. The charge contributions are slightly higher for transition states where X = F. Dipolar polarization contributions are always small and only reflect the molecular orientation change when the nucleophile displaces the leaving group and, therefore, can be neglected. The same occurs for contributions from the X atoms. Only atoms aligned with the reaction axis Nu^--C-LG contribute to the total intensity. Almost all of the infrared intensities are determined by electron transfers from the nucleophile to carbon and subsequently from carbon to the leaving group. The mechanism of charge transfer revealed by the CCTDP model is consistent with the well-accepted reaction mechanism. Open-access codes for performing the IRC-TS-CCTDP analysis are described and provided for potential users in the Supporting Information.