A comparative ab initio study of Br2•− and Br2 water clusters

Energetics and spectroscopic studies of CNO (-) (H 2 O) n clusters and the temperature dependencies of the isomers: An approach based on a combined recipe of parallel tempering and quantum chemical methods

A system associated with several number of weak interactions supports numerous number of stable structures within a narrow range of energy. Often, a deterministic search method fails to locate the global minimum geometry as well as important local minimum isomers for such systems. Therefore, in this work, the stochastic search technique, namely parallel tempering, has been executed on the quantum chemical surface of theCNO (-) (H 2 O) n system for n = 1 -8 to generate global minimum as well as several number of local minimum isomers. IR spectrum can act as the fingerprint property for such system to be identified. Thus, IR spectroscopic features have also been included in this work. Vertical detachment energy has also been calculated to obtain clear information about number of water molecules in several spheres around the central anion. In addition, in a real experimental scenario, not only the global but also the local minimum isomers play an important role in determining the average value of a particular physically observable property. Therefore, the relative conformational populations have been determined for all the evaluated structures for the temperature range between 20K and 400K. Further to understand the phase change behavior, the configurational heat capacities have also been calculated for different sizes.

The anomalous halogen bonding interactions between chlorine and bromine with water in clathrate hydrates

Clathrate hydrate phases of Cl₂ and Br₂ guest molecules have been known for about 200 years. The crystal structure of these phases was recently re-determined with high accuracy by single crystal X-ray diffraction. In these structures, the water oxygen–halogen atom distances are determined to be shorter than the sum of the van der Waals radii, which indicates the action of some type of non-covalent interaction between the dihalogens and water molecules. Given that in the hydrate phases both lone pairs of each water oxygen atom are engaged in hydrogen bonding with other water molecules of the lattice, the nature of the oxygen–halogen interactions may not be the standard halogen bonds characterized recently in the solid state materials and enzyme–substrate compounds. The nature of the halogen–water interactions for the Cl₂ and Br₂ molecules in two isolated clathrate hydrate cages has recently been studied with ab initio calculations and Natural Bond Order analysis (Ochoa-Resendiz et al. J. Chem. Phys. 2016, 145, 161104). Here we present the results of ab initio calculations and natural localized molecular orbital analysis for Cl₂ and Br₂ guests in all cage types observed in the cubic structure I and tetragonal structure I clathrate hydrates to characterize the orbital interactions between the dihalogen guests and water. Calculations with isolated cages and cages with one shell of coordinating molecules are considered. The computational analysis is used to understand the nature of the halogen bonding in these materials and to interpret the guest positions in the hydrate cages obtained from the X-ray crystal structures.

Nucleophilicity and electrophilicity of the C(sp3)–H bond: methane and ethane binary complexes with iodine

Cooperativity is key to the ability of saturated hydrocarbons to act both as nucleophilic and electrophilic centres.

An adaptive mutation simulated annealing based investigation of Coulombic explosion and identification of dissociation patterns in (CO2)n2+ clusters

In this communication, we would like to discuss the advantages of adaptive mutation simulated annealing (AMSA) over standard simulated annealing (SA) in studying the Coulombic explosion of (CO₂)_n²⁺ clusters for n = 20–68, where ‘n’ is the size of the cluster.

An investigation on the structure, spectroscopy and thermodynamic aspects of Br2((-))(H2O)n clusters using a conjunction of stochastic and quantum chemical methods

In this work we obtained global as well as local structures of Br2((-))(H2O)n clusters for n = 2 to 6 followed by the study of IR-spectral features and thermochemistry for the structures. The way adopted by us to obtain structures is not the conventional one used in most cases. Here we at first generated excellent quality pre-optimized structures by exploring the suitable empirical potential energy surface using stochastic optimizer simulated annealing. These structures are then further refined using quantum chemical calculations to obtain the final structures, and spectral and thermodynamic features. We clearly showed that our approach results in very quick and better convergence which reduces the computational cost and obviously using the strategy we are able to get one [i.e. global] or more than one [i.e. global and local(s)] energetically lower structures than those which are already reported for a given cluster size. Moreover, IR-spectral results and the evolutionary trends in interaction energy, solvation energy and vertical detachment energy for global structures of each size have also been presented to establish the utility of the procedure employed.

Molecular dynamics simulation of halogen bonding in Cl2, BrCl, and mixed Cl2/Br2 clathrate hydrates

Clathrate hydrate phases of dihalogen molecules have properties that differ from those of other guest molecules of similar size. The water oxygen–chlorine distances in the structure I (sI) Cl2 hydrate are smaller than the sum of the van der Waals radii of oxygen and chlorine. Bromine hydrate forms a unique clathrate hydrate structure that is not seen in other guest substances. In mixed Cl2/Br2 structure I hydrate, the water oxygen–bromine distances are also smaller than the sum of the oxygen and bromine van der Waals radii. We previously studied the structure of three dihalogen clathrate hydrates using single crystal X-ray diffraction and described these structural features in terms of halogen bonding between the dihalogen and water molecules. In this work, we perform molecular dynamics simulations of cubic sI Cl2, mixed Cl2/Br2, and BrCl clathrate hydrate phases. We perform quantum chemical computations on the dihalogen molecules to determine the nature of σ-hole near the halogen atoms. We fit the electrostatic potential of the molecules to point charge models including dummy atoms that represent σ-holes adjacent to the halogen molecules. Molecular dynamics simulations are used to determine the lattice constants, radial distribution functions, and guest dynamics in these phases. We determine the effect of guest size and difference in halogen bonding on the properties of the clathrate hydrate phase. Simulations for the Cl2, BrCl, and mixed Cl2/Br2 hydrates are performed with small cages of the sI clathrate hydrate phases completely full or filled with experimental occupancies with Cl2 guests.

Microhydration effects on the electronic properties of protonated phenol: a theoretical study

The CC2 (second-order approximate coupled cluster method) has been employed to investigate microhydration effect on electronic properties of protonated phenol (PhH(+)) According to the CC2 calculation results on electronic excited states of microhydrated PhH(+), for the S1 and S2 electronic states, which are of (1)ππ* nature and belong to the A' representation of molecular Cs point group, a significant blue shift effect on the S1 and S2 electronic states, which are of 1ππ* nature and belong to the A' representation of molecular Cs point group, in comparison to corresponding transitions on bare cation (PhH(+)), has been predicted. Nevertheless, for the S3-S0 (1A'', 1σπ*) transition, a large red shift effect has been predicted. Furthermore, it has been found that the lowest (1)σπ* state plays a prominent role in the photochemistry of these systems. In the bare protonated phenol, the (1)σπ* state is a bound state with a broad potential curve along the OH stretching coordinate, while it is dissociative in microhydrated species. This indicates to a predissociation of the S1((1)ππ*) state by a low-lying (1)σπ* state, which leads the excited system to a concerted proton-transfer reaction from protonated chromophore to the solvent. The dissociative (1)σπ* state in monohydrated PhH(+) has small barrier, while increasing the solvent molecules up to three removes the barrier and consequently expedites the proton-transfer reaction dynamics.

Exploring water binding motifs to an excess electron via X2(-)(H2O) [X = O, F]

X(2)(-)(H(2)O) [X = O, F] is utilized to explore water binding motifs to an excess electron via ab initio calculations at the MP4(SDQ)/aug-cc-pVDZ + diffs(2s2p,2s2p) level of theory. X(2)(-)(H(2)O) can be regarded as a water molecule that binds to an excess electron, the distribution of which is gauged by X(2). By varying the interatomic distance of X(2), r(X1-X2), the distribution of the excess electron is altered, and the water binding motifs to the excess electron is then examined. Depending on r(X1-X2), both binding motifs of C(s) and C(2v) forms are found with a critical distance of ∼1.37 Å and ∼1.71 Å for O(2)(-)(H(2)O) and F(2)(-)(H(2)O), respectively. The energetic and geometrical features of O(2)(-)(H(2)O) and F(2)(-)(H(2)O) are compared. In addition, various electronic properties of X(2)(-)(H(2)O) are examined. For both O(2)(-)(H(2)O) and F(2)(-)(H(2)O), the C(s) binding motif appears to prevail at a compact distribution of the excess electron. However, when the electron is diffuse, characterized by the radius of gyration in the direction of the X(2) bond axis with a threshold of ∼0.84 Å, the C(2v) binding motif is formed.

Micro-solvation of the Zn2+ ion-a case study

The structures and thermodynamic parameters of hydrated zinc ion clusters incorporating a single zinc ion and up to eighteen water molecules have been determined with a quantum mechanical hybrid density functional, namely B3LYP using cc-PVDZ basis functions for H and O and a split valence 6-31G (d, p) basis function for Zn. The geometries for all the zinc ion water clusters are optimized with several initial guess structures and without imposing any initial symmetry restriction. Zinc metal ion is found to be preferably four coordinated for smaller sizes of hydrated cluster but attains an octahedral coordination for larger sizes (n >or= 11) of hydrated cluster. Structures with seven or more than seven water molecules attached directly to the central zinc ion are not found. The calculated gas phase coordination number in the first solvation sphere of a large hydrated zinc metal ion is found to be six and the same is also confirmed in the force field based classical molecular dynamics simulation study for an aqueous zinc ion and thus confirms the experimental findings. The equilibrium zinc-oxygen distance of 2.09328 A at the present B3LYP level of study is in excellent agreement with the X-ray diffraction result of 2.093 +/- 0.002 A for a hexahydrated zinc cluster.

Theoretical studies on photoelectron and IR spectral properties of Br2∙−(H2O)n clusters

We report vertical detachment energy (VDE) and IR spectra of Br2.-.(H2O)n clusters (n=1-8) based on first principles electronic structure calculations. Cluster structures and IR spectra are calculated at Becke's half-and-half hybrid exchange-correlation functional (BHHLYP) with a triple split valence basis function, 6-311++G(d,p). VDE for the hydrated clusters is calculated based on second order Moller-Plesset perturbation (MP2) theory with the same set of basis function. On full geometry optimization, it is observed that conformers having interwater hydrogen bonding among solvent water molecules are more stable than the structures having double or single hydrogen bonded structures between the anionic solute, Br2.-, and solvent water molecules. Moreover, a conformer having cyclic interwater hydrogen bonded network is predicted to be more stable for each size hydrated cluster. It is also noticed that up to four solvent H2O units can reside around the solute in a cyclic interwater hydrogen bonded network. The excess electron in these hydrated clusters is localized over the solute atoms. Weighted average VDE is calculated for each size (n) cluster based on statistical population of the conformers at 150 K. A linear relationship is obtained for VDE versus (n+3)(-1/3) and bulk VDE of Br2.- aqueous solution is calculated as 10.01 eV at MP2 level of theory. BHHLYP density functional is seen to make a systematic overestimation in VDE values by approximately 0.5 eV compared to MP2 data in all the hydrated clusters. It is observed that hydration increases VDE of bromine dimer anion system by approximately 6.4 eV. Calculated IR spectra show that the formation of Br2.--water clusters induces large shifts from the normal O-H stretching bands of isolated water keeping bending modes rather insensitive. Hydrated clusters, Br2.-.(H2O)n, show characteristic sharp features of O-H stretching bands of water in the small size clusters.

Microhydration of Cs+ ion: A density functional theory study on Cs+–(H2O)n clusters (n=1–10)

Structure, energy enthalpy, and IR frequency of hydrated cesium ion clusters, Cs+-(H2O)n (n=1-10), are reported based on all electron calculations. Calculations have been carried out with a hybrid density functional, namely, Becke's three-parameter nonlocal hybrid exchange-correlation functional B3LYP applying cc-PVDZ correlated basis function for H and O atoms and a split valence 3-21G basis function for Cs atom. Geometry optimizations for all the cesium ion-water clusters have been carried out with several possible initial guess structures following Newton-Raphson procedure leading to many conformers close in energy. The calculated values of binding enthalpy obtained from present density functional based all electron calculations are in good agreement with the available measured data. Binding enthalpy profile of the hydrated clusters shows a saturation behavior indicating geometrical shell closing in hydrated structure. Significant shifts of O-H stretching bands with respect to free water molecule in IR spectra of hydrated clusters are observed in all the hydrated clusters.

Microhydration of cytosine and its radical anion: Cytosine∙(H2O)n (n=1–5)

Microhydration effects on cytosine and its radical anion have been investigated theoretically, by explicitly considering various structures of cytosine complexes with up to five water molecules. Each successive water molecule (through n=5) is bound by 7-10 kcal mol(-1) to the relevant cytosine complex. The hydration energies are uniformly higher for the analogous anion systems. While the predicted vertical detachment energy (VDE) of the isolated cytosine is only 0.48 eV, it is predicted to increase to 1.27 eV for the lowest-lying pentahydrate of cytosine. The adiabatic electron affinity (AEA) of cytosine was also found to increase from 0.03 to 0.61 eV for the pentahydrate, implying that the cytosine anion, while questionable in the gas phase, is bound in aqueous solution. Both the VDE and AEA values for cytosine are smaller than those of uracil and thymine for a given hydration number. These results are in qualitative agreement with available experimental results from photodetachment-photoelectron spectroscopy studies of Schiedt et al. [Chem. Phys. 239, 511 (1998)].

Structure, energy, and IR spectra of I2∙−.nH2O clusters (n=1–8): A theoretical study

The authors report theoretical results on structure, bonding, energy, and infrared spectra of iodine dimer radical anion hydrated clusters, I(2) (-).nH(2)O (n=1-8), based on a systematic study following density functional theory. Several initial guess structures are considered for each size cluster to locate minimum energy conformers with a Gaussian 6-311++G(d,p) split valence basis function (triple split valence 6-311 basis set is applied for iodine). It is observed that three different types of hydrogen bonded structures, namely, symmetrical double hydrogen bonding, single hydrogen bonding, and interwater hydrogen bonding structures, are possible in these hydrated clusters. But conformers having interwater hydrogen bonding arrangements are more stable compared to those of double or single hydrogen bonded structures. It is also noticed that up to four solvent H(2)O units can reside around the solute in interwater hydrogen bonding network. At the maximum six H(2)O units are independently linked to the dimer anion having four double hydrogen bonding and two single hydrogen bonding, suggesting the hydration number of I(2) (-) to be 6. However, conformers having H(2)O units independently linked to the iodine dimer anion are not the most stable structures. In all these hydrated clusters, the odd electron is found to be localized over two I atoms and the two atoms are bound by a three-electron hemi bond. The solvation, interaction, and vertical detachment energies are calculated for all I(2) (-).nH(2)O clusters. Energy of interaction and vertical detachment energy profiles show stepwise saturation, indicating geometrical shell closing in the hydrated clusters, but solvation energy profile fails to show such behavior. A linear correlation is observed between the calculated energy of interaction and vertical detachment energy. It is observed that formation of I(2) (-)-water cluster induces significant shifts from the normal O-H stretching modes of isolated H(2)O. However, bending mode of H(2)O remains insensitive to the successive addition of solvent H(2)O units. Weighted average energy profiles and IR spectra are reported for all the hydrated clusters based on the statistical population of individual conformers at room temperature.

Microhydration shell structure in Cl2∙−∙nH2O clusters: A theoretical study

We present the results of a detailed study on structure and electronic properties of hydrated cluster Cl2*-.nH2O (n = 1-7) based on a nonlocal density functional, namely, Becke's [J. Chem. Phys. 98, 1372 (1993)] half and half hybrid exchange-correlation functional with a split valence 6-311++G(d,p) basis function. Geometry optimizations for all the clusters are carried out with various possible initial guess structures without any symmetry restriction. Several minimum energy structures (conformers) are predicted with a small difference in total energy. There is a competition between the binding of solvent H2O units with Cl2*- dimer radical anion directly through ion-molecule interaction and forming interwater hydrogen-bonding network in Cl2*-.nH2O (n > or = 2) hydrated cluster. Structure having interwater H-bonded network is more stable over the structure where H2O units are connected to the solute dimer radical anion Cl2*- rather independently either by single or double H bonding in a particular size (n) of hydrated cluster Cl2*-.nH2O. At the maximum four solvent H2O units reside in interwater H-bonding network present in these hydrated clusters. It is observed that up to six H2O units are independently linked to the anion having four double H bondings and two single H bondings suggesting the primary hydration number of Cl2*- to be 6. In all these clusters, the odd electron is found to be mostly localized over the two Cl atoms and these two atoms are bound by a three-electron hemibond. Calculated interaction (between solute and different water clusters) and vertical detachment energy profiles show saturation at n = 6 in the hydrated cluster Cl2*-.nH2O (n = 1-7). However, calculated solvation energy increases with the increase in number of solvent H2O molecules in the cluster. Interaction energy varies linearly with vertical detachment energy for the hydrated clusters Cl2*-.nH2O (n < or = 6). Calculation of the vibration frequencies show that the formation of Cl2*(-)-water clusters induces significant shifts from the normal stretching modes of isolated water. A clear difference in the pattern of IR spectra is observed in the O-H stretching region of water from hexa- to heptahydrated cluster.