Systematic Study of Different Types of Interactions in α-, ß- and γ-Cyclodextrin: Quantum Chemical Investigation.

In this work, comprehensive ab initio quantum chemical calculations using the DFT level of theory were performed to characterize the stabilization interactions (H-bonding and hyperconjugation effects) of two stable symmetrical conformations of α-, ß-, and γ-cyclodextrins (CDs). For this purpose, we analyzed the electron density using "Atom in molecules" (AIM), "Natural Bond Orbital" (NBO), and energy decomposition method (CECA) in 3D and in Hilbert space. We also calculated the H-bond lengths and OH vibrational frequencies. In every investigated CD, the quantum chemical descriptors characterizing the strength of the interactions between the H-bonds of the primary OH (or hydroxymethyl) and secondary OH groups are examined by comparing the same quantity calculated for ethylene glycol, α-d-glucose (α-d-Glcp) and a water cluster as reference systems. By using these external standards, we can characterize more quantitatively the properties of these bonds (e.g., strength). We have demonstrated that bond critical points (BCP) of intra-unit H-bonds are absent in cyclodextrins, similar to α-d-Glcp and ethylene glycol. In contrast, the CECA analysis showed the existence of an exchange (bond-like) interaction between the interacting O…H atoms. Consequently, the exchange interaction refers to a chemical bond, namely the H-bond between two atoms, unlike BCP, which is not suitable for its detection.

Chloride ions as integral parts of hydrogen bonded networks in aqueous salt solutions: the appearance of solvent separated anion pairs.

Hydrogen bonding to chloride ions has been frequently discussed over the past 5 decades. Still, the possible role of such secondary intermolecular bonding interactions in hydrogen bonded networks has not been investigated in any detail. Here we consider computer models of concentrated aqueous LiCl solutions and compute the usual hydrogen bond network characteristics, such as distributions of cluster sizes and of cyclic entities, both for models that take and do not take chloride ions into account. During the analysis of hydrogen bonded rings, a significant amount of 'solvent separated anion pairs' have been detected at high LiCl concentrations. It is demonstrated that taking halide anions into account as organic constituents of the hydrogen bonded network does make the interpretation of structural details significantly more meaningful than when considering water molecules only. Finally, we compare simulated structures generated by 'good' and 'bad' potential sets on the basis of the tools developed here, and show that this novel concept is, indeed, also helpful for distinguishing between reasonable and meaningless structural models.

On Dipole Moments and Hydrogen Bond Identification in Water Clusters.

It is demonstrated that the localized orbitals calculated for a water cluster have small delocalization tails along the hydrogen bonds, that are crucial in determining the resulting dipole moments of the system. (By cutting them, one gets much smaller dipole moments for the individual monomers-close to the values one obtains by using a Bader-type analysis.) This means that the individual water monomers can be delimited only in a quite fuzzy manner, and the electronic charge density in a given point cannot be assigned completely to that or another molecule. Thus, one arrives to the brink of breaking the concept of a water cluster consisting of individual molecules. The analysis of the tails of the localized orbitals can also be used to identify the pairs of water molecules actually forming hydrogen bonds.

Hierarchy of the Collective Effects in Water Clusters.

The results of dipole moment as well as of intra- and intermolecular bond order calculations indicate the big importance of collective electrostatic effects caused by the nonimmediate environment in liquid water models. It is also discussed how these collective effects are built up as consequences of the electrostatic and quantum chemical interactions in water clusters.

Large polarization but small electron transfer for water around Al(3+) in a highly hydrated crystal.

Precise molecular-level information on the water molecule is precious, since it affects our interpretation of the role of water in a range of important applications of aqueous media. Here we propose that electronic structure calculations for highly hydrated crystals yield such information. Properties of nine structurally different water molecules (19 independent OO hydrogen bonds) in the Al(NO3)3·9H2O crystal have been calculated from DFT calculations. We combine the advantage of studying different water environments using one and the same compound and method (instead of comparing a set of independent experiments, each with its own set of errors) with the advantage of knowing the exact atomic positions, and the advantage of calculating properties that are difficult to extract from experiment. We find very large Wannier dipole moments for H2O molecules surrounding the cations: 4.0-4.3 D (compared to our calculated value of 1.83 D in the gas phase). These are induced by the ions and the H-bonds, while other water interactions and the relaxation of the internal water geometry in fact decrease the dipole moments. We find a good correlation between the water dipole moment and the OO distances, and an even better (non-linear) correlation with the average electric field over the molecule. Literature simulation data for ionic aqueous solutions fit quite well with our crystalline 'dipole moment vs. OO distance' curve. The progression of the water and cation charges from 'small clusters ⇒ large clusters ⇒ the crystal' helps explain why the net charges on all the water molecules are so small in the crystal.

Outstanding Properties of the Hydration Shell around ß-d-Glucose: A Computational Study.

Ab initio molecular dynamics (AIMD) simulations have been performed on aqueous solutions of four simple sugars, α-d-glucose, ß-d-glucose, α-d-mannose, and α-d-galactose. Hydrogen-bonding (HB) properties, such as the number of donor- and acceptor-type HB-s, and the lengths and strengths of hydrogen bonds between sugar and water molecules, have been determined. Related electronic properties, such as the dipole moments of water molecules and partial charges of the sugar O atoms, have also been calculated. The hydrophilic and hydrophobic shells were characterized by means of spatial distribution functions. ß-d-Glucose was found to form the highest number of hydrophilic and the smallest number of hydrophobic connections to neighboring water molecules. The average sugar-water H-bond length was the shortest for ß-d-glucose, which suggests that these are the strongest such H-bonds. Furthermore, ß-d-glucose appears to stand out in terms of the symmetry properties of both its hydrophilic and hydrophobic hydration shells. In summary, in all aspects considered here, there seems to be a correlation between the distinct characteristics of ß-d-glucose reported here and its outstanding solubility in water. Admittedly, our findings represent only some of the important factors that influence the solubility.

Reactivity models of hydrogen activation by frustrated Lewis pairs: synergistic electron transfers or polarization by electric field?

Two alternative qualitative reactivity models have recently been proposed to interpret the facile heterolytic cleavage of H2 by frustrated Lewis pairs (FLPs). Both models assume that the reaction takes place via reactive intermediates with preorganized acid/base partners; however, they differ in the mode of action of the active centers. In the electron transfer (ET) model, the hydrogen activation is associated with synergistic electron donation processes with the simultaneous involvement of active centers and the bridging hydrogen, showing similarity to transition-metal-based and other H2-activating systems. In contrast, the electric field (EF) model suggests that the heterolytic bond cleavage occurs as a result of polarization by the strong EF present in the cavity of the reactive intermediates. To assess the applicability of the two conceptually different mechanistic views, we examined the structural and electronic rearrangements as well as the EFs along the H2 splitting pathways for a representative set of reactions. The analysis reveals that electron donations developing already in the initial phase are general characteristics of all studied reactions, and the related ET model provides qualitative interpretation for the main features of the reaction pathways. On the other hand, several arguments have emerged that cast doubt on the relevance of EF effects as a conceptual basis in FLP-mediated hydrogen activation.

Hydrogen bond network topology in liquid water and methanol: a graph theory approach.

Networks are increasingly recognized as important building blocks of various systems in nature and society. Water is known to possess an extended hydrogen bond network, in which the individual bonds are broken in the sub-picosecond range and still the network structure remains intact. We investigated and compared the topological properties of liquid water and methanol at various temperatures using concepts derived within the framework of graph and network theory (neighbour number and cycle size distribution, the distribution of local cyclic and local bonding coefficients, Laplacian spectra of the network, inverse participation ratio distribution of the eigenvalues and average localization distribution of a node) and compared them to small world and Erdos-Rényi random networks. Various characteristic properties (e.g. the local cyclic and bonding coefficients) of the network in liquid water could be reproduced by small world and/or Erdos-Rényi networks, but the ring size distribution of water is unique and none of the studied graph models could describe it. Using the inverse participation ratio of the Laplacian eigenvectors we characterized the network inhomogeneities found in water and showed that similar phenomena can be observed in Erdos-Rényi and small world graphs. We demonstrated that the topological properties of the hydrogen bond network found in liquid water systematically change with the temperature and that increasing temperature leads to a broader ring size distribution. We applied the studied topological indices to the network of water molecules with four hydrogen bonds, and showed that at low temperature (250 K) these molecules form a percolated or nearly-percolated network, while at ambient or high temperatures only small clusters of four-hydrogen bonded water molecules exist.

Connecting Diffraction Experiments and Network Analysis Tools for the Study of Hydrogen-Bonded Networks.

A self-consistent scheme is presented that is applicable for revealing details of the microscopic structure of hydrogen-bonded liquids, including the description of the hydrogen-bonded network. The scheme starts with diffraction measurements, followed by molecular dynamics simulations. Computational results are compared with the experimentally accessible information on the structure, which is most frequently the total scattering structure factor. In the case of an at least semiquantitative agreement between experiment and simulation, sets of particle coordinates from the latter may be exploited for revealing nonmeasurable structural details. Calculations of some properties concerning the hydrogen-bonded network are also described, in the order of increasing complexity: starting with the definition of a hydrogen bond, first and second neighborhoods are described via spatial correlations functions. Attention is then turned to cyclic and noncyclic hydrogen-bonded clusters, before cluster size distributions and percolation are discussed. We would like to point out that, as a result of applying the novel protocol, these latter, rather abstract, quantities become consistent with diffraction data: it may thus be argued that the approach reviewed here is the first one that establishes a direct link between measurements and elements of network theories. Applications for liquid water, simple alcohols, and alcohol-water liquid mixtures demonstrate the usefulness of the aforementioned characteristics. The procedure can readily be applied to more complicated hydrogen-bonded networks, like mixtures of polyols (diols, triols, sugars, etc.) and water, and complex aqueous solutions of even larger molecules (even of proteins).

A gold cyano complex in nitromethane: MD simulation and X-ray diffraction.

The solvation structure around the dicyanoaurate(I) anion (Au(CN)2-) in a dilute nitromethane (CH3NO2) solution is presented from X-ray diffraction measurements and molecular dynamics simulation (NVT ensemble, 460 nitromethane molecules at room temperature). The simulations are based on a new solute-solvent force-field fitted to a training set of quantum-chemically derived interaction energies. Radial distribution functions from experiment and simulation are in good agreement. The solvation structure has been further elucidated from MD data. Several shells can be identified. We obtain a solvation number of 13-17 nitromethane molecules with a strong preference to be oriented with their methyl groups towards the solute.

Two Faces of the Two-Phase Thermodynamic Model.

The quantum harmonic model and the two-phase thermodynamic method (2PT) are widely used to obtain quantum-corrected properties such as isobaric heat capacities or molar entropies. 2PT heat capacities were calculated inconsistently in the literature. For water, the classical heat capacity was also considered, but for organic liquids, it was omitted. We reanalyzed the performance of different quantum corrections on the heat capacities of common organic solvents against experimental data. We have pointed out serious flaws in previous 2PT studies. The vibrational density of states was calculated incorrectly causing a 39% relative error in diffusion coefficients and 45% error in the 2PT heat capacities. The wrong conversion of isobaric and isochoric heat capacities also caused about 40% error but in the other direction. We have introduced the concept of anharmonic correction (AC), which is simply the deviation of the classical heat capacity from that of the harmonic oscillator model. This anharmonic contribution is around +30 to 40 J/(mol K) for water depending on the water model and -8 to -10 J/(mol K) for hydrocarbons and halocarbons. AC is unrealistically large, +40 J/(K mol) for alcohols and amines, indicating some deficiency of the OPLS force field. The accuracy of the computations was also assessed with the determination of the self-diffusion coefficients.

Properties of Hydrogen-Bonded Networks in Ethanol-Water Liquid Mixtures as a Function of Temperature: Diffraction Experiments and Computer Simulations.

New X-ray and neutron diffraction experiments have been performed on ethanol-water mixtures as a function of decreasing temperature, so that such diffraction data are now available over the entire composition range. Extensive molecular dynamics simulations show that the all-atom interatomic potentials applied are adequate for gaining insight into the hydrogen-bonded network structure, as well as into its changes on cooling. Various tools have been exploited for revealing details concerning hydrogen bonding, as a function of decreasing temperature and ethanol concentration, like determining the H-bond acceptor and donor sites, calculating the cluster-size distributions and cluster topologies, and computing the Laplace spectra and fractal dimensions of the networks. It is found that 5-membered hydrogen-bonded cycles are dominant up to an ethanol mole fraction xeth = 0.7 at room temperature, above which the concentrated ring structures nearly disappear. Percolation has been given special attention, so that it could be shown that at low temperatures, close to the freezing point, even the mixture with 90% ethanol (xeth = 0.9) possesses a three-dimensional (3D) percolating network. Moreover, the water subnetwork also percolates even at room temperature, with a percolation transition occurring around xeth = 0.5.

Interaction of carbon monoxide with Au(111) modified by ion bombardment: a surface spectroscopy study under elevated pressure.

Gold based model systems exhibiting the structural versatility of nanoparticle ensembles and being accessible for surface spectroscopic investigations are expected to provide new information about the adsorption of carbon monoxide, a key process influencing the CO oxidation activity of this noble metal in nanoparticulate form. Accordingly, in the present work the interaction of CO is studied with an ion bombardment modified Au(111) surface by means of a combination of photoelectron spectroscopy (XPS and UPS), sum frequency generation vibrational spectroscopy (SFG), and scanning tunneling microscopy (STM). While no adsorption was found on intact Au(111), data collected on the ion bombarded surface at cryogenic temperatures indicated the presence of stable CO adsorbates below 190 K. A quantitative evaluation of the C 1s XPS spectra and the surface morphology explored by STM revealed that the step edge sites created by ion bombardment are responsible for CO adsorption. The identification of the CO binding sites was confirmed by density functional theory (DFT) calculations. Annealing experiments up to room temperature showed that at temperatures above 190 K unstable adsorbates are formed on the surface under dynamic exposure conditions that disappeared immediately when gaseous CO was removed from the system. Spectroscopic data as well as STM records revealed that prolonged CO exposure at higher pressures of up to 1 mbar around room temperature facilitates massive atomic movements on the roughened surface, leading to its strong reordering toward the structure of the intact Au(111) surface, accompanied by the loss of the CO binding capacity.

Hydrogen bonded network properties in liquid formamide.

Molecular dynamics simulations have been performed for liquid formamide using two different types of potential model (OPLS, Cordeiro). The structural results obtained from simulation were compared to experimental (x-ray and neutron diffraction measurements) outcomes. A generally good agreement for both models examined has been found, but in the hydrogen bonded region (2.9 A) the Cordeiro model shows a slightly better fit. Besides the evaluation of partial radial distribution functions, orientational correlation functions and energy distribution functions, describing the hydrogen bonded structure, have been calculated based on the statistical analysis of configurations, resulting into a new insight in the clustering properties and topology of hydrogen bonded network. It has been shown that in liquid formamide exists a continuous hydrogen bonded network and from the analysis of the distribution of small rings revealed the ring size distribution in liquid formamide. Our study resulted that the ring size distribution of the hydrogen bonded liquid formamide shows a broad distribution with a maximum around 11. It has been found that the topology in formamide is significantly different than in water.

Nuclear Quantum Effects from the Analysis of Smoothed Trajectories: Pilot Study for Water.

Nuclear quantum effects have significant contributions to thermodynamic quantities and structural properties; furthermore, very expensive methods are necessary for their accurate computation. In most calculations, these effects, for instance, zero-point energies, are simply neglected or only taken into account within the quantum harmonic oscillator approximation. Herein, we present a new method, Generalized Smoothed Trajectory Analysis, to determine nuclear quantum effects from molecular dynamics simulations. The broad applicability is demonstrated with the examples of a harmonic oscillator and different states of water. Ab initio molecular dynamics simulations have been performed for ideal gas up to the temperature of 5000 K. Classical molecular dynamics have been carried out for hexagonal ice, liquid water, and vapor at atmospheric pressure. With respect to the experimental heat capacity, our method outperforms previous calculations in the literature in a wide temperature range at lower computational cost than other alternatives. Dynamic and structural nuclear quantum effects of water are also discussed.

Modeling tetraalkylammonium halide salts in water: how hydrophobic and electrostatic interactions shape the thermodynamic properties.

The explicit water molecular dynamics simulation was used to study tetramethylammonium and tetraethylammonium chloride and bromide solutions in water at 298 K. The outcome of the simulations in the form of various distribution functions was used to construct the solvent-averaged potentials between interacting molecules. In the next step, which involved the Ornstein-Zernike integral equation theory in the hypernetted chain approximation, these potentials were used to calculate the osmotic coefficients. We showed that this approach is able to explain the experimental results for the osmotic pressure of these salts.

Roles of the ether oxygen in hydration of tetrahydrofuran studied by IR, NMR, and DFT calculation methods.

We studied the concentration dependence of nu(C-H)'s in IR and (1)J(C,H) in NMR for binary water-tetrahydrofuran (THF) mixtures and found different trends for the two types of CH(2) groups in the five-membered ring. The changes of the nu(C-O) spectra showed that complexes of THF associated with water are formed, in which the number of water molecules increases with the water concentration. We suggested that hydration proceeds through the formation of 1:1, and 1:2 complexes of [THF:water] up to X(H(2)O) approximately 0.9, where X(H)((2))(O) is the mole fraction of the water in the mixtures. We carried out ab initio MO and DFT calculations to optimize the geometries of a THF dimer as a model of THF molecules in pure liquid, and 1:1 and 1:2 complexes of [THF:water] to simulate observed concentration dependence of nu(C-H)'s in IR and (1)J(C,H) in NMR. The changes of the calculated nu(C-H) spectra and (1)J(C,H) values for the optimized complexes are in agreement with those observed with varying X(H)((2))(O), supporting our proposal. From the vibrational and NBO analyses of the optimized complexes, the observed blue shift of nu(C-H)'s and the increase of (1)J(C,H) for the CH(2) groups neighboring to the ether oxygen were explained in terms of the changes in the stereoelectronic effect, resulting from HO-H...O< hydrogen bonding. The optimized 1:2-complex contains two weak C-H...OH(2) hydrogen bonds, and blue shift of nu(C-H)'s and increase of (1)J(C,H) were demonstrated from the same analyses of the complexes. This result of simulation also supports that the blue shift of nu(C-H)'s and increase of (1)J(C,H) observed for both the type of CH(2) groups at 0.6 X(H)((2))(O) < 0.9 are attributed to these interactions. On the basis of all these results, we propose that the formation of the 1:2-complex involving weak C-H...OH(2) hydrogen bonds is responsible dominantly for the hydrophobic hydration of THF.

Solution structure of NaNO3 in water: diffraction and molecular dynamics simulation study.

The structure of a series of aqueous sodium nitrate solutions (1.9-7.6 M) was studied using a combination of experimental and theoretical methods. The results obtained from diffraction (X-ray, neutron) and molecular dynamics simulation have been compared and the capabilities and limitations of the methods in describing solution structure are discussed. For the solutions studied, diffraction methods were found to perform very well in description of hydration spheres of the sodium ion but do not yield detailed structural information on the anion's hydration structure. Molecular dynamics simulations proved to be a suitable tool in the detailed interpretation of the hydration sphere of ions, ion pair formation, and bulk structure of solutions.

Benchmarking polarizable molecular dynamics simulations of aqueous sodium hydroxide by diffraction measurements.

Results from molecular dynamics simulations of aqueous hydroxide of varying concentrations have been compared with experimental structural data. First, the polarizable POL3 model was verified against neutron scattering using a reverse Monte Carlo fitting procedure. It was found to be competitive with other simple water models and well suited for combining with hydroxide ions. Second, a set of four polarizable models of OH- were developed by fitting against accurate ab initio calculations for small hydroxide-water clusters. All of these models were found to provide similar results that robustly agree with structural data from X-ray scattering. The present force field thus represents a significant improvement over previously tested nonpolarizable potentials. Although it cannot in principle capture proton hopping and can only approximately describe the charge delocalization within the immediate solvent shell around OH-, it provides structural data that are almost entirely consistent with data obtained from scattering experiments.

Molecular Dynamics Simulation Studies of the Temperature-Dependent Structure and Dynamics of Isopropanol-Water Liquid Mixtures at Low Alcohol Content.

Series of molecular dynamics simulations for 2-propanol-water mixtures, as a function of temperature (between freezing and room temperature) and composition (xip = 0, 0.5, 0.1, and 0.2), have been performed for temperatures reported in the only available experimental structure study. It is shown that when the all-atom optimized potentials for liquid simulations interatomic potentials for the alcohol are combined with the TIP4P/2005 water model, then the near-quantitative agreement with measured X-ray data, in the reciprocal space, can be achieved. Such an agreement justifies detailed investigations of structural, energetic, and dynamic properties on the basis of simulation trajectories. Here, we focus on characteristics related to hydrogen bonds (HB): cluster-, and in particular, ring formation, energy distributions, and lifetimes of HB-s have been scrutinized for the entire system, as well as for the water and isopropanol subsystems. It is demonstrated that similar to ethanol-water mixtures, the occurrence of 5-membered-hydrogen-bonded rings are significant, particularly at higher alcohol concentrations. Concerning HB energetics, an intriguing double maximum appears on the alcohol-alcohol HB energy distribution function. HB lifetimes have been found significantly longer in the mixtures than they are in pure liquids.