A Re-evaluation of Factors Controlling the Nature of Complementary Hydrogen-Bonded Networks

The energy profiles of hydrogen-bonded heterocyclic aromatics have been decomposed into atomistic energy contributions using the Interacting Quantum Atoms (IQA) method. The resulting energy contributions have been sequenced by the Relative Energy Gradient (REG) approach to determine their influence upon the shape of these energy profiles. The results show inadequacies in Jorgensen's secondary interaction hypothesis (SIH). A novel method of finding a condensed analogy for the interaction between the molecules is presented. The findings of this work further doubt the validity of the SIH, and reinforce previous warnings against its misguided use.     

Hydrogen‐Bond Cooperative Effects in Small Cyclic Water Clusters as Revealed by the Interacting Quantum Atoms Approach

The cooperative effects of hydrogen bonding in small water clusters (H₂O)_n (n=3–6) have been studied by using the partition of the electronic energy in accordance with the interacting quantum atoms (IQA) approach. The IQA energy splitting is complemented by a topological analysis of the electron density (ρ( r )) compliant with the quantum theory of atoms‐in‐molecules (QTAIM) and the calculation of electrostatic interactions by using one‐ and two‐electron integrals, thereby avoiding convergence issues inherent to a multipolar expansion. The results show that the cooperative effects of hydrogen bonding in small water clusters arise from a compromise between: 1) the deformation energy (i.e., the energy necessary to modify the electron density and the configuration of the nuclei of the isolated water molecules to those within the water clusters), and 2) the interaction energy (E_int) of these contorted molecules in (H₂O)_n. Whereas the magnitude of both deformation and interaction energies is enhanced as water molecules are added to the system, the augmentation of the latter becomes dominant when the size of the cluster is increased. In addition, the electrostatic, classic, and exchange components of E_int for a pair of water molecules in the cluster (H₂O)_n?1 become more attractive when a new H₂O unit is incorporated to generate the system (H₂O)_n with the last‐mentioned contribution being consistently the most important part of E_int throughout the hydrogen bonds under consideration. This is opposed to the traditional view, which regards hydrogen bonding in water as an electrostatically driven interaction. Overall, the trends of the delocalization indices, δ(Ω,Ω′), the QTAIM atomic charges, the topology of ρ( r ), and the IQA results altogether show how polarization, charge transfer, electrostatics, and covalency contribute to the cooperative effects of hydrogen bonding in small water clusters. It is our hope that the analysis presented in this paper could offer insight into the different intra‐ and intermolecular interactions present in hydrogen‐bonded systems.     

On the Relationship between Hydrogen Bond Strength and the Formation Energy in Resonance-Assisted Hydrogen Bonds

Resonance-assisted hydrogen bonds (RAHB) are intramolecular contacts that are characterised by being particularly energetic. This fact is often attributed to the delocalisation of π electrons in the system. In the present article, we assess this thesis via the examination of the effect of electron-withdrawing and electron-donating groups, namely −F, −Cl, −Br, −CF₃, −N(CH₃)₂, −OCH₃, −NHCOCH₃ on the strength of the RAHB in malondialdehyde by using the Quantum Theory of Atoms in Molecules (QTAIM) and the Interacting Quantum Atoms (IQA) analyses. We show that the influence of the investigated substituents on the strength of the investigated RAHBs depends largely on its position within the π skeleton. We also examine the relationship between the formation energy of the RAHB and the hydrogen bond interaction energy as defined by the IQA method of wave function analysis. We demonstrate that these substituents can have different effects on the formation and interaction energies, casting doubts regarding the use of different parameters as indicators of the RAHB formation energies. Finally, we also demonstrate how the energy density can offer an estimation of the IQA interaction energy, and therefore of the HB strength, at a reduced computational cost for these important interactions. We expected that the results reported herein will provide a valuable understanding in the assessment of the energetics of RAHB and other intramolecular interactions.     

Molecular Photochemistry: Recent Developments in Theory

Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light‐induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non‐adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.     

The Activation Strain Model in the Light of Real Space Energy Partitions

The distortion/interaction or activation strain model (ASM) of chemical reactivity is examined in real space through the interacting quantum atoms (IQA) approach. Attention is paid to the role that the geometrically constrained ASM structures of the fragments play in the chemical interpretation of the driving forces that lead to a given reaction channel. These fictitious intermediate states are necessary in the ASM, but IQA may or may not use them at will. Similarities and differences are highlighted by studying the endo/exo preference rules of simple [4+2] Diels–Alder cycloadditions. Although overall the agreement is reasonable, we warn about a blind use of the plain ASM if no further energy decomposition analyses of its interaction energy are done.     

Cooperativity and Anticooperativity in Ion-Water Interactions: Implications for the Aqueous Solvation of Ions

Non-additive effects in hydrogen bonds (HB) take place as a consequence of electronic charge transfers. Therefore, it is natural to expect cooperativity and anticooperativity in ion-water interactions. Nevertheless, investigations on this matter are scarce. This paper addresses the interactions of (i) the cations Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺, and Ca²⁺ together with (ii) the anions F⁻, Cl⁻, Br⁻, NO₃⁻ and SO₄²⁻ with water clusters (H₂O)_n, n=1–8, and the effects of these ions on the HBs within the complete molecular adducts. We used quantum chemical topology tools, specifically the quantum theory of atoms in molecules and the interacting quantum atoms energy partition to investigate non-additive effects among the interactions studied herein. Our results show a decrease on the interaction energy between ions and the first neighbouring water molecules with an increment of the coordination number. We also found strong cooperative effects in the interplay between HBs and ion-dipole interactions within the studied systems. Such cooperativity affects considerably the interactions among ions with their first and second solvation shells in aqueous environments. Overall, we believe this article provides valuable information about how ion-dipole contacts interact with each other and how they relate to other interactions, such as HBs, in the framework of non-additive effects in aqueous media.     

On the strength of hydrogen bonding within water clusters on the coordination limit

Hydrogen bonds (HB) are arguably the most important noncovalent interactions in chemistry. We study herein how differences in connectivity alter the strength of HBs within water clusters of different sizes. We used for this purpose the interacting quantum atoms energy partition, which allows for the quantification of HB formation energies within a molecular cluster. We could expand our previously reported hierarchy of HB strength in these systems (Phys. Chem. Chem. Phys., 2016, 18 , 19557) to include tetracoordinated monomers. Surprisingly, the HBs between tetracoordinated water molecules are not the strongest HBs despite the widespread occurrence of these motifs (e.g., in ice I_h). The strongest HBs within H₂O clusters involve tricoordinated monomers. Nonetheless, HB tetracoordination is preferred in large water clusters because (a) it reduces HB anticooperativity associated with double HB donors and acceptors and (b) it results in a larger number of favorable interactions in the system. Finally, we also discuss (a) the importance of exchange-correlation to discriminate among the different examined types of HBs within H₂O clusters, (b) the use of the above-mentioned scale to quickly assess the relative stability of different isomers of a given water cluster, and (c) how the findings of this research can be exploited to indagate about the formation of polymorphs in crystallography. Overall, we expect that this investigation will provide valuable insights into the subtle interplay of tri- and tetracoordination in HB donors and acceptors as well as the ensuing interaction energies within H₂O clusters.     

When hydroquinone meets methoxy radical: Hydrogen abstraction reaction from the viewpoint of interacting quantum atoms

Interacting Quantum Atoms methodology is used for a detailed analysis of hydrogen abstraction reaction from hydroquinone by methoxy radical. Two pathways are analyzed, which differ in the orientation of the reactants at the corresponding transition states. Although the discrepancy between the two barriers amounts to only 2 kJ/mol, which implies that the two pathways are of comparable probability, the extent of intra‐atomic and inter‐atomic energy changes differs considerably. We thus demonstrated that Interacting Quantum Atoms procedure can be applied to unravel distinct energy transfer routes in seemingly similar mechanisms. Identification of energy components with the greatest contribution to the variation of the overall energy (intra‐atomic and inter‐atomic terms that involve hydroquinone's oxygen and the carbon atom covalently bound to it, the transferring hydrogen and methoxy radical's oxygen), is performed using the Relative energy gradient method. Additionally, the Interacting Quantum Fragments approach shed light on the nature of dominant interactions among selected fragments: both Coulomb and exchange‐correlation contributions are of comparable importance when considering interactions of the transferring hydrogen atom with all other atoms, whereas the exchange‐correlation term dominates interaction between methoxy radical's methyl group and hydroquinone's aromatic ring. This study represents one of the first applications of Interacting Quantum Fragments approach on first order saddle points. © 2018 Wiley Periodicals, Inc.     

The N‐Arylamino Conjugation Effect in the Photochemistry of Fluorescent Protein Chromophores and Aminostilbenes

To understand the nonradiative decay mechanism of fluorescent protein chromophores in solutions, a systematic comparison of a series of (Z)‐4‐(N‐arylamino)benzylidene‐2,3‐imidazolinones (ABDIs: 2P , 2PP , 2OM , and 2OMB ) and the corresponding trans‐4‐(N‐arylamino)‐4′‐cyanostilbenes (ACSs: 1P , 1PP , 1OM , and 1OMB ) was performed. We have previously shown that the parameter Φ_f+2 Φ_tc, in which Φ_f and Φ_tc are the quantum yields of fluorescence and trans→cis photoisomerization, respectively, is an effective probe for evaluating the contribution of twisted intramolecular charge transfer (TICT) states in the excited decays of trans‐aminostilbenes, including the push–pull ACSs. One of the criteria for postulating the presence of a TICT state is Φ_f+2 Φ_tc?1.0, because its formation is decoupled with the C?C bond (τ) torsion pathway and its decay is generally nonradiative. Our results show that the same concept also applies to ABDIs 2 with the parameter Φ_f+2 Φ_ZE in which Φ_ZE is the quantum yield of Z→E photoisomerization. We conclude that the τ torsion rather than the C? C bond (φ) torsion is responsible for the nonradiative decays of ABDIs 2 in aprotic solvents (hexane, THF, acetonitrile). The phenyl‐arylamino C? N bond (ω) torsion that leads to a nonradiative TICT state is important only for 2OM in THF and acetonitrile. If the solvent is protic (methanol and 10–20 % H₂O in THF), a new nonradiative decay channel is present for ABDIs 2 , but not for ACSs 1 . It is attributed to internal conversion (IC) induced by solvent (donor)–solute (acceptor) hydrogen‐bonding (HB) interactions. The possible HB modes and the concept of τ torsion‐coupled proton transfer are also discussed.     

Charged versus Neutral Hydrogen‐Bonded Complexes: Is There a Difference in the Nature of the Hydrogen Bonds?

A theoretical study on some carboxylic acid dimers formed by positively or negatively charged molecules has been carried out by using DFT methods. The resulting dimers possess either a charge of +2 or ?2. In addition, the corresponding neutral complexes have also been considered. The electron density distribution described by the atoms in molecules and the natural bond orbital methods, as well as the electric field maps of the systems, have been analyzed and compared without finding significant differences between the neutral and ionic complexes. The interaction energy along the dissociation path of the charged dimers shows both a local minimum and a local maximum, defining a stability region between them. When this energetic profile is recalculated by removing the repulsion between the charged groups, it resembles to those of the neutral molecules. Hence, the characteristics of the charged dimers are similar to those of the neutral ones: the addition of a repulsion term for the charged groups permits to retrieve the energetic profiles dependence with the distance in the charged system. The interacting quantum atom (IQA) method has been used to calculate the interaction energy terms, including the classic Coulombic term between the whole molecules and the corresponding of the carboxylic acid groups. The IQA results show repulsive electrostatic interactions when the whole molecules are considered in the ionic complexes, but attractive ones between the carboxylic groups in both neutral and ionic complexes.     

Statistical Equivalence of Quantum Chemical Methods for Energy Distribution in Water Clusters

Methods of quantum chemistry are instrumental in understanding molecular structures and properties. However, the results demonstrate significant variability, which is difficult to predict and rationalize. The fundamental question is whether some molecular systems exhibit properties invariant with respect to the computational method. The idea explored here is that collective properties of statistical ensembles should be more robust than characteristics of individual molecules and their arbitrary sets. This effect is demonstrated for the complete set of hydrogen-bond topologies of the dodecahedral water cluster (H₂O)₂₀. Non-Gaussian energy distributions produced by various methods have the same functional form despite strong differences in mean values and standard deviations. The conclusion is tested on methods of different complexity and origin employing a number of criteria. A linear mapping between the energies produced by different methods is discussed. The significance of the results is in establishing a collective equivalence property of quantum chemical methods.     

A Quantitative View of Charge Transfer in the Hydrogen Bond: The Water Dimer Case

The hydrogen bond represents a fundamental intermolecular interaction that binds molecules in vapor and liquid water. A crucial and debated aspect of its electronic structure and chemistry is the charge transfer (CT) accompanying it. Much effort has been devoted, in particular, to the study of the smallest prototype system, the water dimer, but even here results and interpretations differ widely. In this paper, we reassess CT in the water dimer by using charge‐displacement analysis. Besides a reliable estimate of the amount of CT (14.6 me) that characterizes the system, our study provides an unambiguous context, and very useful bounds, within which CT effects may be evaluated, crucially including the associated energy stabilization.     

On the role of substituent in noncovalent functionalization of graphene and organophosphor recognition: IQA and SAPT perspective

Despite importance of integrating organic molecules with graphene to fabricate graphene‐based electronic devices, the role of substituents and interface stabilizing forces are poorly understood. In this work, the interactions of 7,7,8,8‐tetracyanoquinodimethane (TCNQ), 2,3,5,6‐tetrafluoro‐7,7,8,8‐tetracyanoquinodimethane (F4TCNQ), hydroquinone (Q), and tetrafluorohydroquinone (TFQ) with graphene have been investigated by means of interacting quantum atoms and SAPT(DFT). In addition, in context of potential design of a graphene‐based sensor for detection of the nerve agent sarin, we studied the interaction of graphene and the organic molecules with the dimethyl methylphosphonate (DMMP)—the molecule that mimics sarin. The results show that the organic molecules attach to graphene via C(sp²)?C(sp²), C(sp²)?C(sp) and H?π bonds. In addition, they trap DMMP via various linkages such as hydrogen, lonepair?π and H?π . The quantum effects play a significant role. The Pauli repulsion is responsible for p‐doping of graphene. The substituents are stabilized on graphene by the exchange‐correlation energy. The fluorination of the benzenoid ring raises the electron‐sharing . The through space and through bond effects of the fluorine atoms (‐F) increase the classical attraction of the cyano groups and benzenoid ring with graphene, respectively. When comparing performance of the ab initio and DFT methods, MP2 predicts too much attraction due to well‐known overestimation of the dispersion energy by the uncoupled dispersion component for benzene rings, while ω B97xD functional and SAPT(DFT) provide weaker interaction energies, in good agreement with each other.     

Directionality of Halogen Bonds: An Interacting Quantum Atoms (IQA) and Relative Energy Gradient (REG) Study

Interacting Quantum Atoms (IQA) and Interacting Quantum Fragments (IQF) analyses are used to study (X=Cl and Br) model complexes in order to determine the origin of halogen bond directionality. IQA allows for the calculation of intra- and interatomic classical and exchange-correlation energies, which can be used to determine the energetic nature of the changes that occur when deviating from the preferred halogen bond approach. The Relative Energy Gradient (REG) method is also applied to rank the IQA energies and reveal which energy contributions best describe the total behavior of the system. Indeed, all the pairwise interactions and atomic self-energies are angularly dependent; some terms favor the linear structure and some tend toward nonlinear arrangements. For instance, when the C−X−N angle is altered, the halogen-nitrogen interaction energy behaves like the total energy of the system while the carbon-nitrogen interaction works against the total energy profile. Furthermore, the REG values reveal that the contribution of the halogen-nitrogen interaction to the total behavior of the system is small. Instead, the secondary interactions (e. g., fluorine-nitrogen and carbon-hydrogen interactions) and atomic self-energies are mainly responsible for the angular preference of these halogen bonds. Finally, IQF calculations followed by REG analysis reveal the importance of the self-energy of the fragments.     

Fluorescence Quenching of Benzaldehyde in Water by Hydrogen Atom Abstraction

We computed the mechanism of fluorescence quenching of benzaldehyde in water through relaxed potential energy surface scans. Time‐dependent density functional theory calculations along the protonation coordinate from water to benzaldehyde reveal that photoexcitation to the bright ππ* (S₃) state is immediately followed by ultrafast decay to the nπ* (S₁) state. Evolving along this state, benzaldehyde (BA) abstracts a hydrogen atom, resulting in a BAH^. and OH^. radical pair. Benzaldehyde does not act as photobase in water, but abstracts a hydrogen atom from a nearby solvent molecule. The system finally decays back to the ground state by non‐radiative decay and an electron transfers back to the OH^. radical. Proton transfer from BAH⁺ to OH^? restores the initial situation, BA in water.     

Stereoselectivity of the Hydrogen‐Atom Transfer in Benzophenone–Tyrosine Dyads: An Intramolecular Kinetic Solvent Effect

To be or not to be solvated is the decisive parameter that controls the photoinduced hydrogen‐atom transfer in diastereomeric ketone/phenol dyads. A kinetic solvent effect that refers to hydrogen bonding between the phenol and the solvent is suggested to be the main source of the stereoselective discrimination in the hydrogen transfer (see figure).