On the supramolecular interactions into a pH- and Metal-Actuated Molecular Shuttle: some insights from QTAIM modeling.

Supramolecular contacts responsible for chemical interaction of cucurbit[7]uril (CB[7]) macrocycle on a Tolyl-Viologen-Phenylene-Imidazole (T-VPI) molecular thread, at acid pH (T-VPI-H+) or after Ag+ cation addition (T-VPI-Ag+), are analytically addressed in a computational framework combining Quantum Theory of Atoms in Molecules (QTAIM) with Density Functional Theory (DFT). In this respect, the crystallographic structure (CCDC number 2217466) is taken as reference condition for addressing the nature of the chemical interactions driving the shuttling of the CB[7] between T and P stations recently observed in dilute water solutions. Beside the host(CB[7]) vs guest(T-VPI-H+ or T-VPI-Ag+) complexation, the coordination sphere of the Ag+ cation is also investigated by means of local electronic energy density - H(r) - descriptors. The derived non-covalent interaction patterns are found to support diagnostic 1H NMR signals used for detecting the mutual position of the CB[7] along the axle. This work highlights the potentialities of a QTAIM based approach in the characterization of supramolecular and metal-complexation effects in molecular aggregates such as not-interlocked synthetic molecular shuttles.

Quantum Information Patterns Between Atoms in a Molecule.

Quantum information theory provides a powerful toolbox of descriptors that characterize many-electron systems based on quantum information patterns between open quantum systems. Despite the wealth of insights gained in the con- densed matter community, the use of these descriptors to study interactions between atoms in a molecule remains limited. In this study, we develop a quantum information framework for molecules that characterizes the quantum in- formation patterns between quantum atoms as defined in the Quantum Theory of Atoms in Molecules. We show that quantum information analyses capture key properties of quantum atoms and how they interact with their molec- ular environment. Additionally, we show that the presence of bond critical points can remain invariant despite large changes in the quantum information patterns between the quantum atoms. Our findings indicate that quantum infor- mation theory can shed a new light on molecular electronic structure.

Effect of Terminal Hydroxy Groups on the Structural Properties of Acyl Thioureas.

The 1-acyl thiourea family [R1C(O)NHC(S)NR2R3] exhibits the flexibility to incorporate a wide variety of substituents into their structure. The structural attributes of these compounds are intricately tied to the type and extent of substitution. In the case of 3-mono-substituted thioureas (R2 = H), the conformational behavior is predominantly shaped by the presence of an intramolecular N-H···O=C hydrogen bond. This study delves into the structural consequences stemming from the inclusion of substituents possessing hydrogen-donor capabilities within four novel 1-acyl-3-mono-substituted thiourea derivatives. A comprehensive suite of analytical techniques, encompassing FTIR, Raman spectroscopy, multinuclear (1H and 13C) NMR spectroscopy, single-crystal X-ray diffraction, and supported by computational methods, notably NBO (Natural Bond Orbital) population analysis, Hirshfeld analysis, and QTAIM (Quantum Theory of Atoms in Molecules), was harnessed to scrutinize and characterize these compounds. In the crystalline state, these compounds exhibit an intricate interplay of intermolecular interactions, prominently featuring an expansive network of hydrogen bonds between the hydroxy (-OH) groups and the carbonyl and thiocarbonyl bonds within the 1-acyl thiourea fragment. Notably, the topological analysis underscores significant distinctions in the properties of the acyl thiourea fragment and the intramolecular >C=O···H-N bond when transitioning from the isolated molecule to the crystalline environment.

Exchange-correlation effects in interatomic energies for pure density functionals and their application to the molecular energy prediction.

In this proof-of-concept paper, we show how exchange-correlation effects can be simply recovered for interatomic energies within the interacting quantum atoms decomposition when local, gradient generalized, or meta-gradient generalized approximations are used in density functional theory (DFT) calculations. We also demonstrate how inhomogeneity and non-local effects can be introduced even from a pure local scheme, without resorting to any orbital information. Finally, we provide numerical evidence on a database of selected energetic molecules that this decomposition scheme can be efficiently used to build accurate models for the prediction of molecular energies from an initial "cheap" DFT calculation.

Does Aromaticity Play a Role in Electronic and Structural Properties of YBn (n=2-14) Clusters?

Nanoclusters exhibit electronic, optical, and magnetic properties that differ significantly from those of extended and molecular systems with comparable stoichiometries. In this work, we examined the structural, energetic, and electronic characteristics of yttrium-doped boron clusters (YBn,  n =2-14) with robust wavefunction analysis tools. Special emphasis is placed on the elucidation of the potential aromatic character exhibited by the resultant molecules and how it can affect their chemical bonding and stability. Our results revealed that the YBn stability is governed by the maximization of the ionic Y-B interactions. This is evidenced from the lowest-energy conformations, which manifest as half-sandwich structures wherein the majority of boron atoms are bonded to yttrium. The stabilization of such chemical contacts comes at the expense of a notorious depletion of the Y local electron density, crystallizing in a considerable ionic character, close to Y2+ + Bn2-. Such a  charge transfer is coupled to the enhancement of the electron delocalization within the YBn lattice, resulting in quite remarkable local and global aromatic characters. Altogether, this study shows how the toolkit of real space chemical bonding descriptors can offer valuable insights into the structural and electronic properties,  of YBn clusters, contributing to a better understanding of their behavior.

Chemical bonding in phase-change chalcogenides.

Almost all phase-change memory materials (PCM) contain chalcogen atoms, and their chemical bonds have been denoted both as 'electron-deficient' [sometimes referred to as 'metavalent'] and 'electron-rich' ['hypervalent', multicentre]. The latter involve lone-pair electrons. We have performed calculations that can discriminate unambiguously between these two classes of bond and have shown that PCM have electron-rich, 3c-4e ('hypervalent') bonds. Plots of charge transferred between (ET) and shared with (ES) neighbouring atoms cannot on their own distinguish between 'metavalent' and 'hypervalent' bonds, both of which involve single-electron bonds. PCM do not exhibit 'metavalent' bonding and are not electron-deficient; the bonding is electron-rich of the 'hypervalent' or multicentre type.

Determining the Factors Accounting for Reaction Selectivity: A Relative Energy Gradient - Interacting Quantum Atoms and Natural Bonding Orbitals Study.

Identifying the main physicochemical properties accounting for the course of a reaction is of utmost importance to rationalize chemical syntheses. To this aim, the relative energy gradient (REG) method is an appealing approach because it is an unbiased and automatic process to extract the most relevant pieces of energy information. Initially formulated within the interacting quantum atoms (IQA) framework for a single reaction, here we extend the REG method to natural bond orbitals (NBO) analysis and to the case of two competitive processes. This development enables the determination of the driving forces of any chemical selectivity. We illustrate the extended REG method on the case study of ring opening in cyclobutenes, which is an important instance of the so-called torquoselectivity.

Exploring the Efficiency of Algerian Kaolinite Clay in the Adsorption of Cr(III) from Aqueous Solutions: Experimental and Computational Insights.

The current study comprehensively investigates the adsorption behavior of chromium (Cr(III)) in wastewater using Algerian kaolinite clay. The structural and textural properties of the kaolinite clay are extensively characterized through a range of analytical methods, including XRD, FTIR, SEM-EDS, XPS, laser granulometry, N2 adsorption isotherm, and TGA-DTA. The point of zero charge and zeta potential are also assessed. Chromium adsorption reached equilibrium within five minutes, achieving a maximum removal rate of 99% at pH 5. Adsorption equilibrium is modeled using the Langmuir, Freundlich, Temkin, Elovich, and Dubinin-Radushkevitch equations, with the Langmuir isotherm accurately describing the adsorption process and yielding a maximum adsorption capacity of 8.422 mg/g for Cr(III). Thermodynamic parameters suggest the spontaneous and endothermic nature of Cr(III) sorption, with an activation energy of 26.665 kJ/mol, indicating the importance of diffusion in the sorption process. Furthermore, advanced DFT computations, including COSMO-RS, molecular orbitals, IGM, RDG, and QTAIM analyses, are conducted to elucidate the nature of adsorption, revealing strong binding interactions between Cr(III) ions and the kaolinite surface. The integration of theoretical and experimental data not only enhances the understanding of Cr(III) removal using kaolinite but also demonstrates the effectiveness of this clay adsorbent for wastewater treatment. Furthermore, this study highlights the synergistic application of empirical research and computational modeling in elucidating complex adsorption processes.

Theoretical Insights into Bifurcated Intramolecular Dihydrogen Bonds.

Two-ring intramolecular π-electron delocalization assisted dihydrogen bonds existing in (1Z,4Z)-1,4-dipentene-3-bora-1,5-diol and its symmetrically substituted derivatives have been analysed here since the MP2/6-311++G(d,p) calculations on these systems were performed. The influence of the coexistence of two intramolecular dihydrogen bonded rings in these molecular structures on properties of intramolecular dihydrogen bonds as well as on the π-electron delocalization within these rings was investigated. The comparison with corresponding structures of typical two-ring, so-called resonance-assisted, RAHB, systems was performed. The results of calculations show that such rings' coexistence leads to the weakening of dihydrogen bonds, similarly as for the typical two-ring RAHB systems. The Quantum Theory of ''Atoms in Molecules'' (QTAIM) was also applied here to get more details about the nature of dihydrogen bonds. Correlations between dihydrogen bond strength measures and other energetic, geometrical and topological parameters were also analysed. It was found that characteristics of bond critical points as well as of ring critical points are useful to estimate the strength of intramolecular dihydrogen bonds in two-ring dihydrogen bonded systems discussed here. The Natural Bond Orbital, NBO, approach parameters are also discussed as useful ones to describe properties of dihydrogen bonded systems.

Response of the mechanical and chiral character of ethane to ultra-fast laser pulses.

A pair of simulated left and right circularly polarized ultra-fast laser pulses of duration 20 femtoseconds that induce a mixture of excited states are applied to ethane. The response of the electron dynamics is investigated within the next generation quantum theory of atoms in molecules (NG-QTAIM) using third-generation eigenvector-trajectories which are introduced in this work. This enables an analysis of the mechanical and chiral properties of the electron dynamics of ethane without needing to subject the C-C bond to external torsions as was the case for second-generation eigenvector-trajectories. The mechanical properties, in particular, the bond-flexing and bond-torsion were found to increase depending on the plane of the applied laser pulses. The bond-flexing and bond-torsion, depending on the plane of polarization, increases or decreases after the laser pulses are switched off. This is explainable in terms of directionally-dependent effects of the long-lasting superpositions of excited states. The chiral properties correspond to the ethane molecule being classified as formally achiral consistent with previous NG-QTAIM investigations. Future planned investigations using ultra-fast circularly polarized lasers are briefly discussed.

Implementation of machine learning protocols to predict the hydrolysis reaction properties of organophosphorus substrates using descriptors of electron density topology.

Prediction of catalytic reaction efficiency is one of the most intriguing and challenging applications of machine learning (ML) algorithms in chemistry. In this study, we demonstrated a strategy for utilizing ML protocols applied to Quantum Theory of Atoms In Molecules (QTAIM) parameters to predict the ability of the A17 L47K catalytic antibody to covalently capture organophosphate pesticides. We found that the novel "composite" DFT functional B97-3c could be effectively employed for fast and accurate initial geometry optimization, aligning well with the input dataset creation. QTAIM descriptors proved to be well-established in describing the examined dataset using density-based and hierarchical clustering algorithms. The obtained clusters exhibited correlations with the chemical classes of the input compounds. The precise physical interpretation of the QTAIM properties simplifies the explanation of feature impact for both supervised and unsupervised ML protocols. It also enables acceleration in the search for entries with desired properties within large databases. Furthermore, our findings indicated that Ridge Regression with Laplacian kernel and CatBoost Regressor algorithms demonstrated suitable performance in handling small datasets with non-trivial dependencies. They were able to predict the actual reaction barrier values with a high level of accuracy. Additionally, the CatBoost Classifier proved reliable in discriminating between "active" and "inactive" compounds.

Aza-Michael Addition in Explicit Solvent: A Relative Energy Gradient-Interacting Quantum Atoms Study.

Aza-Michael additions are key reactions in organic synthesis. We investigate, from a theoretical and computational point of view, several examples ranging from weak to strong electrophiles in dimethylsulfoxide treated as explicit solvent. We use the REG-IQA method, which is a quantum topological energy decomposition (Interacting Quantum Atoms, IQA) coupled to a chemical-interpretation calculator (Relative Energy Gradient, REG). We focus on the rate-limiting addition step in order to unravel the different events taking place in this step, and understand the influence of solvent on the reaction, with an eye on predicting the Mayr electrophilicity. For the first time, a link is established between an REG-IQA analysis and experimental values.

Theoretical investigation of anion perfluorocubane.

CONTEXT: In this work, we did a theoretical exploration of C8F8 (Ib) and its anion radical analogue (IIb) in this work. By investigating the thermochemistry of electron capture, we find that the free energy associated with the conversion of C8H8 (Ia) into its anion radical analogue IIa is of the order of + 92.83 kcal.mol-1, while the conversion of Ib into IIb is - 6.42 kcal.mol-1. Therefore, species IIb is thermodynamically more stable than its neutral analogue. Natural bond orbitals (NBO) analyses revealed that compound Ib exhibits a relative electronic stability as a function of intramolecular delocalisations of the type [Formula: see text] of the order of 2.70 kcal.mol-1. Similar delocalizations for Ia are energetically lower (1.45 kcal.mol-1). Topological analyses of compounds Ib and IIb indicate that the addition of an electron to Ib enhances the covalency of the C-C bond, as can be seen by the reduction in the ellipticity of the C-C bond. The opposite is observed for Ia, whose addition of the electron (leading to IIa) reduces the covalency of the C-C bond. By comparing the free and packaged forms of the species, it is found that, in the crystalline form, the system will present greater relative stability due to the dispersive interactions involved, as evidenced by non-covalent interactions (NCI) analysis. Finally, it was possible to verify that the manifestation of the current density with a lower paratropic and less antiaromatic character in Ib and IIb point to C8F8 as a strong candidate for electron capture. METHODS: Geometry optimization calculations were carried out, for all monomer structures using the hybrid functional B3LYP-D3 and the 6-31+G(d,p) basis set. To determine the formation thermochemistry of the ions, electronic energy corrections was performed using the DLPNO-CCSD(T)/aug-cc-pVTZ/C method. Starting from the optimised forms, shielding, nuclear magnetic resonance (NMR) spectra employing gauge-independent atomic orbital (GIAO), and NBO calculations were performed for these monomers, using the PBE0 functional and the pCSseg-2 atomic basis set. The magnetochemical analysis of ring currents was performed using the GIMIC formalism. For the topological analysis, it was applied the combination DLPNO-CCSD(T)/aug-cc-pVTZ/C, previously used for correcting the electronic energy.

Ga···C Triel Bonds-Why They Are Not Strong Enough to Change Trigonal Configuration into Tetrahedral One: DFT Calculations on Dimers That Occur in Crystal Structures.

Structures characterized by the trigonal coordination of the gallium center that interacts with electron rich carbon sites are described. These interactions may be classified as Ga···C triel bonds. Their properties are analyzed in this study since these interactions may be important in numerous chemical processes including catalytical activities; additionally, geometrical parameters of corresponding species are described. The Ga···C triel bonds discussed here, categorized also as the π-hole bonds, do not change the trigonal configuration of the gallium center into the tetrahedral one despite total interactions in dimers being strong; however, the main contribution to the stabilization of corresponding structures comes from the electrostatic forces. The systems analyzed theoretically here come from crystal structures since the Cambridge Structural Database, CSD, search was performed to find structures where the gallium center linked to CC bonds of Lewis base units occurs. The majority structures found in CSD are characterized by parallel, stacking-like arrangements of species containing the Ga-centers. The theoretical results show that interactions within dimers are not classified as the three-centers links as in a case of typical hydrogen bonds and numerous other interactions. The total interactions in dimers analyzed here consist of several local intermolecular atom-atom interactions; these are mainly the Ga···C links. The DFT results are supported in this study by calculations with the use of the quantum theory of atoms in molecules, QTAIM, the natural bond orbital, NBO, and the energy decomposition analysis, EDA, approaches.

Electronic and dynamical properties of non-covalent diatomic aggregates formed by He with neutral and ionic Li and Be.

CONTEXT AND RESULTS: This work aims to study the influence of the absence and presence of permanent charges on the electronic and dynamical properties of the non-covalent bound diatomic systems involving He and Li, Be as neutral and ionic partners. The charge displacement results suggest that in the formation of HeLi[Formula: see text], HeBe[Formula: see text], and HeBe[Formula: see text], the neutral He atom undergoes, in the electric field of the ion, a pronounced electronic polarization, and the natural bond order theoretical approach indicates that in the formation of the molecular orbital He acts as a weak electron donor. The energy decomposition analysis provides the dispersion and induction components as the attractive leading terms controlling the stability of all systems, confirming that the formed bond substantially maintains a non-covalent nature which is also supported by the Quantum Theory of Atoms in Molecules (QTAIM) analysis. Finally, it was found that the HeLi and HeBe neutral systems are unstable under any condition, HeLi[Formula: see text] and HeBe[Formula: see text] ionic systems are stable below 317K and 138K, respectively, while the HeBe[Formula: see text] system becomes unstable only after 3045K. COMPUTATIONAL AND THEORETICAL TECHNIQUES: The potential energy curves and interactions in all systems were studied theoretically based on coupled-cluster singles and doubles method with perturbative inclusion of triples CCSD(T) method with an aug-cc-pV5Z basis set. More precisely, it was determined the potential energy curves describing the stability of the HeLi, HeLi[Formula: see text], HeBe, HeBe[Formula: see text], and HeBe[Formula: see text] systems, the charge displacement within the formed adducts, the decomposition of their total interaction energy, the topological analysis of their bonds, their rovibrational energies, their spectroscopic constants and lifetimes.

Dynamic and Static Nature of XH-∗-π and YX-∗-π (X = F, Cl, Br, and I; Y = X and F) in the Distorted π-System of Corannulene Elucidated with QTAIM Dual Functional Analysis.

The dynamic and static nature of the XH-∗-π and YX-∗-π (X = F, Cl, Br, and I; Y = X and F) interactions in the distorted π-system of corannulene (π(C20H10)) is elucidated with a QTAIM dual functional analysis (QTAIM-DFA), where asterisks emphasize the presence of bond critical points (BCPs) on the interactions. The static and dynamic nature originates from the data of the fully optimized and perturbed structures, respectively, in QTAIM-DFA. On the convex side, H in F-H-∗-π(C20H10) and each X in Y-X-∗-π(C20H10) join to C of the central five-membered ring in π(C20H10) through a bond path (BP), while each H in X-H-∗-π(C20H10) does so to the midpoint of C=C in the central five-membered ring for X = Cl, Br, or I. On the concave side, each X in F-X-∗-π(C20H10) also joins to C of the central five-membered ring with a BP for X = H, Cl, Br, and I; however, the interactions in other adducts are more complex than those on the convex side. Both H and X in X-H-∗-π(C20H10) (X = Cl and Br) and both Fs in F-F-∗-π(C20H10) connect to the three C atoms in each central five-membered ring (with three BPs). Two, three, and five BPs were detected for the Cl-Cl, I-H, Br-Br, and I-I adducts, where some BPs do not stay on the central five-membered ring in π(C20H10). The interactions are predicted to have a vdW to CT-MC nature. The interactions on the concave side seem weaker than those on the convex side for X-H-∗-π(C20H10), whereas the inverse trend is observed for Y-X-∗-π(C20H10) as a whole. The nature of the interactions in the π(C20H10) adducts of the convex and concave sides is examined in more detail, employing the adducts with X-H and F-X placed on their molecular axis together with the π(C24H12) and π(C6H6) adducts.

Origin of structure and stability of M@C18 (M = Cu, Ag, and Au) complexes with D9h point group.

Theoretical predictions and recent experimental studies lead to the discovery of an exciting new member of the carbon allotrope family polyynic cyclo[18]carbon (C18 ). Present investigation aims to probe the structure, stability, and properties of coinage metal (M)@C18 complexes using density functional theory (DFT) calculations. The DFT results unequivocally show that even Cu@C18 , Ag@C18 , and Au@C18 complexes substantially preserve the ground state polyynic structure of C18 . It is also worth to mention that only Au@C18 is a stable D9h structure, however the symmetry is distorted in the case of Cu@C18 and Ag@C18 . Due to computational limitations, in this investigation the M@C18 complexes were scrutinized using the C2v sub abelian group of D9h . The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of D9h conformers are a singlet a1 and two same value singlets a1 ⊕ b1 generated from doublet e, respectively. The non-covalent interaction index (NCI), quantum theory of atoms in molecule (QTAIM), and energy decomposition analysis (EDA) vividly explains the interaction between a coinage metal atom and C18 ring. It is found from the results that the stability of Cu@C18 Ag@C18 , and Au@C18 is governed by the attractive electrostatic, orbital and dispersion interaction.

An atoms-in-molecules characterization of the nature of the OO bond in peroxides and nitroxide dimers.

The quantum theory of atoms-in-molecules (QTAIM) method is used to examine the OO bond in peroxides (RO-OR) and nitroxide dimers (R2 NO-ONR2 ), including Fremy's salt. The electron density (ρ), electron kinetic energy density [K(ρ)], and Laplacian of the electron density (∇2 ρ) at bond critical points characterize the nature of the OO bond. The data distinguish OO bonding of two kinds. Large values of ρ and positive ∇2 ρ and K(ρ) suggest that simple peroxides have charge-shift bonds. Nitroxide dimers, with smaller ρ, positive ∇2 ρ, and near-zero K(ρ), show a lack of shared electron density, suggesting there is no conventional OO bonding in these molecules. QTAIM analysis at the B3LYP/6-311+G(d,p) level of theory gives results in agreement with valence-bond theory and X-ray diffraction characterizations of peroxide OO bonds as charge-shift bonds. In contrast, CCSD/cc-pVDZ calculations fail to agree with previous results because of an insufficient, single-determinant treatment of the charge-shift bond.

Inverse Versus Normal Behavior of Interactions, Elucidated Based on the Dynamic Nature with QTAIM Dual-Functional Analysis.

In QTAIM dual-functional analysis, Hb(rc) is plotted versus Hb(rc) - Vb(rc)/2 for the interactions, where Hb(rc) and Vb(rc) are the total electron energy densities and potential energy densities, respectively, at the bond critical points (BCPs) on the interactions in question. The plots are analyzed by the polar (R, θ) coordinate representation for the data from the fully optimized structures, while those from the perturbed structures around the fully optimized structures are analyzed by (θp, κp). θp corresponds to the tangent line of the plot, and κp is the curvature; θ and θp are measured from the y-axis and y-direction, respectively. The normal and inverse behavior of interactions is proposed for the cases of θp > θ and θp < θ, respectively. The origin and the mechanism for the behavior are elucidated. Interactions with θp < θ are typically found, although seldom for [F-I-∗-F]-, [MeS-∗-TeMe]2+, [HS-∗-TeH]2+ and CF3SO2N-∗-IMe, where the asterisks emphasize the existence of BCPs in the interactions and where [Cl-Cl-∗-Cl]- and CF3SO2N-∗-BrMe were employed as the reference of θp > θ. The inverse behavior of the interactions is demonstrated to arise when Hb(rc) - Vb(rc)/2 and when the corresponding Gb(rc), the kinetic energy densities at BCPs, does not show normal behavior.

Distributed functional-group polarizabilities in polypeptides and peptide clusters toward accurate prediction of electro-optical properties of biomacromolecules.

CONTEXT: Aiming at accurately predicting electro-optical properties of biomolecules, this work presents distributed atomic and functional-group polarizability tensors for a series of polypeptides and peptide clusters constructed from glycine and its residuals. By partitioning the electron density using the quantum theory of atoms in molecules, we demonstrated a very good transferability of the group polarizabilities. We were able to identify and extract the most efficient functional groups capable of generating the largest electrical susceptibility in condensed phases. Both the isotropic polarizability and its anisotropy were used to understand the way functional groups act as sources of linear optical responses, how they interact with each other reinforcing the macroscopic optical behavior within the material, and how covalent bonds and non-covalent interactions, such as hydrogen bonds, determine refractive indices and birefringence. Particular attention is devoted to the peptide bonds as they provide links to build biomacromolecules or polymers. An adequate quantum-mechanical treatment of at least the first interaction sphere of a given functional group is required to properly describe the effects of mutual polarization, but we identified optimum cluster size and shape to better estimate polarizabilities and dipole moments of larger molecules or molecular aggregates from the knowledge of the electron density of a central molecule or amino acid residual that is representative of the bulk. The strategy outlined here is a fast yet effective tool for estimating the optical properties of proteins but could eventually find application in the rational design of optical organic materials as well. METHODS: Electronic-structure calculations were performed on the Gaussin16 program at the DFT level using the CAMB3LYP functional and the double-ζ quality Dunning basis set aug-cc-pVDZ. Electron density partitioning followed the concepts of the Quantum Theory of Atoms and Molecules (QTAIM) and was performed using the AIMAll program. The locally developed Polaber routine was applied to calculate dipole moment vectors and polarizability tensors. It was amended to include the effects of the local field on a given central molecule by means of a modified Atom-Dipole Interaction Model (ADIM).