Encapsulation of charged halogens by the 512 water cage.

This study focuses on the encapsulation of the entire series of halides by the 512 cage of twenty water molecules and on the characterization of water to water and water to anion interactions. State-of-the-art computations are used to determine equilibrium geometries, energy related quantities, and thermal stability towards dissociation and to dissect the nature and strength of intermolecular interactions holding the clusters as stable units. Two types of structures are revealed: heavily deformed cages for F- indicating a preference for microsolvation, and slightly deformed cages for the remaining anions indicating a preference for encapsulation. The primary variable dictating the properties of the clusters is the charge density of the central halide, with the most severe effects observed for the F- case. For the remaining halides, the anion may be safely viewed as a sort of "big electron" with little local disruptive power, enough to affect the network of non-covalent hydrogen bonds in the cage, but not enough to break it. Gibbs energies for dissociation either into cavity and halide or into water molecules and halide suggest that, in a similar way as to methane clathrate, a more weakly bonded complex that has been detected in the gas phase, all halide containing clathrate-like structures should be amenable to experimental detection in the gas phase at moderate temperature and pressure conditions.

Interactions and reactivity in crystalline intermediates of mechanochemical cyclorhodation reactions.

There is experimental evidence that solid mixtures of the rhodium dimer [Cp*RhCl2]2 and benzo[h] quinoline (BHQ) produce two different polymorphic molecular cocrystals called 4α and 4ß under ball milling conditions. The addition of NaOAc to the mixture leads to the formation of the rhodacycle [Cp*Rh-(BHQ)Cl], where the central Rh atom retains its tetracoordinate character. Isolate 4ß reacts with NaOAc leading to the same rhodacycle while isolate 4α does not under the same conditions. We show that the puzzling difference in reactivity between the two cocrystals can be traced back to fundamental aspects of the intermolecular interactions between the BHQ and [Cp*RhCl2]2 fragments in the crystalline environment. To support this view, we report a number of descriptors of the nature and strength of chemical bonds and intermolecular interactions in the extended solids and in a cluster model. We calculate formal quantum mechanical descriptors based on electronic structure, electron density, and binding and interaction energies including an energy decomposition analysis. Without exception, all descriptors point to 4ß being a transient structure higher in energy than 4α with larger local and global electrophilic and nucleophilic powers, a more favorable spatial and energetic distribution of the frontier orbitals, and a more fragile crystal structure.

The Role of Crystalline Intermediates in Mechanochemical Cyclorhodation Reactions Elucidated by in-Situ X-ray Powder Diffraction and Computation.

The occurrence of crystalline intermediates in mechanochemical reactions might be more widespread than previously assumed. For example, a recent study involving the acetate-assisted C-H activation of N-Heterocycles with [Cp*RhCl2 ]2 by ball milling revealed the formation of transient cocrystals between the reagents prior to the C-H activation step. However, such crystalline intermediates were only observed through stepwise intervallic ex-situ analysis, and their exact role in the C-H activation process remained unclear. In this study, we monitored the formation of discrete, stoichiometric cocrystals between benzo[h]quinoline and [Cp*RhCl2 ]2 by ball milling using in-situ synchrotron X-ray powder diffraction. This continuous analysis revealed an initial cocrystal that transformed into a second crystalline form. Computational studies showed that differences in noncovalent interactions made the [Cp*RhCl2 ]2 unit in the later-appearing cocrystal more reactive towards NaOAc. This demonstrated the advantage of cocrystal formation before the acetate-assisted metalation-deprotonation step, and how the net cooperative action of weak interactions between the reagents in mechanochemical experiments can lead to stable supramolecular assemblies, which can enhance substrate activation under ball-milling conditions. This could explain the superiority of some mechanochemical reactions, such as acetate-assisted C-H activation, compared to their solution-based counterparts.

Encapsulation of Astatide by a water cage.

A 512 cage of (H2O)20 consisting of 30 hydrogen bonds encapsulates Astatide with little geometrical distortion. The cage is marginally destabilized but the non-covalent interactions are actually strengthened. Host⋯cage interactions in the [At@(H2O)20]- cluster are anti-electrostatic, placing both negatively charged atoms in direct contact as in Atδ-⋯δ-O-Hδ+. An orbital interaction analysis reveals that explicit host⋯cage contacts are "inverted" hydrogen bonds. That is, the same type of donor→acceptor charge transfer as in hydrogen bonding, with no proton bridging the two negative charges.

Microsolvation of electrons by a handful of ammonia molecules.

Microsolvation of electrons in ammonia is studied here via anionic NH3 n - clusters with n = 2-6. Intensive samplings of the corresponding configurational spaces using second-order perturbation theory with extended basis sets uncover rich and complex energy landscapes, heavily populated by many local minima in tight energy windows as calculated from highly correlated coupled cluster methods. There is a marked energetical preference for structures that place the excess electron external to the molecular frame, effectively coordinating it with the three protons from a single ammonia molecule. Overall, as the clusters grow in size, the lowest energy dimer serves as the basic motif over which additional ammonia molecules are attached via unusually strong charge-assisted hydrogen bonds. This is a priori quite unexpected because, on electrostatic grounds, the excess electron would be expected to be in contact with as many protons as possible. Accordingly, a full quantum mechanical treatment of the bonding interactions under the tools provided by the quantum theory of atoms in molecules is carried out in order to dissect and understand the nature of intermolecular contacts. Vertical detachment energies reveal bound electrons even for n = 2.

Microsolvation versus Encapsulation in Mono, Di, and Trivalent Cations.

The effects of the formal charge in the stability and bonding of water cavities when solvating a cation are studied here using [X(H2 O)20 ]q+ clusters starting with the well known 512 isomer of (water)20 , placing a single mono, di, or trivalent Xq+ cation at the interior, and then optimizing and characterizing the resulting clusters. Highly correlated interaction and deformation energies are calculated using the CCSD(T)-DLPNO formalism. Bonding interactions are characterized using the tools provided by the quantum theory of atoms in molecules, natural bond orbitals, and non-covalent surfaces. Our results indicate that water to water hydrogen bonds are sensibly strengthened resulting in strong cooperative effects, which amount to ≈ 2 ${ \approx 2}$  kcal/mol per hydrogen bond in the bare cavity and to larger values for the systems including the cations. Approximate encapsulation, that is, surrounding the cation by a network of hydrogen bonds akin to the well known methane clathrate seems to be preferred by cations with smaller charge densities while microsolvation, that is, cluster structures having explicit X⋯O contacts seem to be preferred by cations with larger charge densities which severely deform the cavity.

To Be or Not To Be? that is the Entropic, Enthalpic, and Molecular Interaction Dilemma in the Formation of (Water)20 Clusters and Methane Clathrate.

A detailed analysis under a comprehensive set of theoretical and computational tools of the thermodynamical factors and of the intermolecular interactions behind the stabilization of a well known set of (water)20 cavities and of the methane clathrate is offered in this work. Beyond the available reports of experimental characterization at extreme conditions of most of the systems studied here, all clusters should be amenable to experimental detection at 1 atm and moderate temperatures since 280 K marks the boundary at which, ignoring reaction paths, formation of all clusters is no longer spontaneous from the 20H2 O→(H2 O)20 and CH4 +20H2 O→CH4 @512 processes. As a function of temperature, a complex interplay leading to the free energy of formation occurs between the destabilizing entropic contributions, mostly due to cluster vibrations, and the stabilizing enthalpic contributions, due to intermolecular interactions and the PV term, is best illustrated by the highly symmetric 512 cage consistently showing signs of stronger intermolecular bonding despite having smaller binding energy than the other clusters. A fluxional wall of attractive non-covalent interactions, arising because of the cumulative effect of a large number of tiny individual charge transfers to the interstitial region, plays a pivotal role stabilizing the CH4 @512 clathrate.

Thermodynamics and Intermolecular Interactions during the Insertion of Anionic Naproxen into Model Cell Membranes.

The insertion process of Naproxen into model dimyristoylphosphatidylcholine (DMPC) membranes is studied by resorting to state-of-the-art classical and quantum mechanical atomistic computational approaches. Molecular dynamics simulations indicate that anionic Naproxen finds an equilibrium position right at the polar/nonpolar interphase when the process takes place in aqueous environments. With respect to the reference aqueous phase, the insertion process faces a small energy barrier of ≈5 kJ mol-1 and yields a net stabilization of also ≈5 kJ mol-1. Entropy changes along the insertion path, mainly due to a growing number of realizable microstates because of structural reorganization, are the main factors driving the insertion. An attractive fluxional wall of noncovalent interactions is characterized by all-quantum descriptors of chemical bonding (natural bond orbitals, quantum theory of atoms in molecules, noncovalent interaction, density differences, and natural charges). This attractive wall originates in the accumulation of tiny transfers of electron densities to the interstitial region between the fragments from a multitude of individual intermolecular contacts stabilizing the tertiary drug/water/membrane system.

Structure, energy, and bonding in anionic water tetramers obtained by exhaustive search.

An analysis of the structures, some energy related properties, and key aspects of the bonding nature of the microsolvated electron with four water molecules is presented. The study is based on an exhaustive potential energy surface scan of the ground state of (H2O)4 - at the UCCSD(T)/6-311(3+,4+)G(d,p)//UMP2/6-311(3+,4+)G(d,p) level. A total of 18 structures, most of them not reported before, spanning in an energy range of 8.8 kcal mol-1 were found. The energetic stability of the clusters is dictated by the effect of the excess electron on their structures, on their partial fragmentation, and on the hydrogen bonds' framework. Vertical detachment energies depend on the number of water molecules holding the excess electron in "direct contact" to their two protons at the same time and, to a lesser extent, also depend on the hydrogen bond sequence in the rest of the structure. In general, hydrogen bonds in (H2O)4 - are of closed shell character, and there are other less common interactions assisted by the excess electron.

Evolution of Bonding during the Insertion of Anionic Ibuprofen into Model Cell Membranes.

Descriptors of chemical bonding derived from five different analysis tools based on quantum mechanics (natural charges, electron density differences, atoms in molecules (AIM), natural bond orbitals (NBO), and non-covalent interactions (NCI) index) consistently afford a picture of a wall of weak, non-covalent intermolecular interactions separating anionic Ibuprofen from the environment. This wall, arising from the cumulative effect of a multitude of individual weak charge transfer interactions to the interstitial region between fragments, stabilizes the drug at all equilibrium positions in the free energy profile for its insertion into model cell membranes. The formal charge in anionic Ibuprofen strengthens all intermolecular interactions, having a particularly strong effect in the network of water to water hydrogen bonds in the solvent. Electron redistribution during the insertion process leads to a sensible reduction of electron delocalization in both the -CO2- group and the aromatic ring of Ibuprofen. Here, we conclusively show that, despite their purely classical origin, randomly chosen configurations from molecular dynamics simulations provide deep insight into the purely quantum nature of bonding interactions.

A Comprehensive Picture of the Structures, Energies, and Bonding in [SO4(H2O)n]2-, n = 1-6.

Two stochastic methods in conjunction with ab initio computations were used to explore the potential energy surfaces for the microsolvation of SO42- with up to six explicit water molecules. At least three water molecules are needed to stabilize the Coulomb repulsion that prevents the existence of isolated SO42-. The formal charge in SO42- is strong enough to induce water dissociation and subsequent microsolvation of the resulting HSO4-, OH- ionic pair. Hydrogen bonds characterized as having complex contributions from covalency and from ionicity are at play stabilizing [SO4(H2O)n]2- clusters. Ionicity and covalency act concomitantly rather than opposedly to strengthen both intermolecular interactions and the resulting O-H bond in HSO4- after proton abstraction.

Microsolvation of F.

A staggering structural diversity for the microsolvation of F- with up to six water molecules is uncovered in this work. Given the structural variety and the proximity in energy among several local minima, we show here that in order to match available experimental data, statistical averages over contributing structures are needed, rather than assigning experimental values to isolated structures. Our results suggest that the formal charge in F- is strong enough as to induce partial and total dissociation of water molecules and to alter the nature of the surrounding network of water to water hydrogen bonds. We provide an extensive analysis of bonding interactions under the NBO and QTAIM formalisms, our main results suggest a complex interplay between ionic and covalent characters for the FH interactions as a function of the separation between the atoms.

Adsorption of Nitrate and Bicarbonate on Fe-(Hydr)oxide.

In this work, we used density functional theory calculations to study the resulting complexes of adsorption and of inner- and outer-sphere adsorption-like of bicarbonate and nitrate over Fe-(hydr)oxide surfaces using acidic, neutral, and basic simulated pH conditions. High-spin states that follow the 5N + 1 (N is the number of Fe atoms, each having five unpaired electrons) rule are preferred. Monodentate mononuclear (MM1) surface complexes are shown to lead to the most favorable thermodynamic adsorption for both bicarbonate and nitrate with -63.91 and -28.25 kJ/mol, respectively, under neutral conditions. Our results suggest that four types of regular and charged-assisted hydrogen bonds are involved in the adsorption process; all of them can be classified as closed-shell (long-range or ionic). The formal charges induce unusually short and strong hydrogen bonds. The ability of high multiplicity states of Fe clusters to adsorb oxyanions in solvated environments arises from orbital interactions: the 4s virtual orbitals in Fe have a large affinity for the 2p-type electron pairs of oxygens.

Stochastic search of the quantum conformational space of small lithium and bimetallic lithium-sodium clusters.

In this paper we report the results obtained by an implementation and application of the simulated annealing optimization procedure to the exploration of the conformational space of small neutral and charged lithium clusters (Li(n)(q), n = 5, 6, 7; q = 0, +/-1) and of the bimetallic lithium/sodium clusters (Li5Na) in their lowest spin states. Our methodology eliminates the structure guessing procedure in the process of generating cluster configurations. We evaluate the quantum energy, typically with the Hartree-Fock Hamiltonian, of randomly generated points in the conformational space and use a modified Metropolis test in the annealing algorithm to generate candidate structures for atomic clusters. The structures are further optimized by analytical methods (gradient following) at the Møller-Plesset second order perturbation theory level (MP2), in conjunction with basis sets including polarization functions with and without diffuse functions. High accuracy ab initio energies at the coupled clusters level, with single, double, and triple substitutions from the Hartree-Fock determinant (CCSD(T)), on the MP2 geometries were calculated and used to establish the relative stability of the isomers within each potential energy surface. Various cluster properties were computed and compared to existing values in order to validate our methods. Our results show excellent agreement with previous experimental and theoretical reports. Even at these small sizes, evidence for 10 new structures never reported before for the lithium clusters and four new structures for the bimetallic clusters is presented.