Ab initio study of the complexes of halogen-containing molecules RX (X=Cl, Br, and I) and NH3: towards understanding the nature of halogen bonding and the electron-accepting propensities of covalently bonded halogen atoms

Ab initio calculations have been performed on a series of complexes formed between halogen-containing molecules and ammonia to gain a deeper insight into the nature of halogen bonding. It appears that the dihalogen molecules form the strongest halogen-bonded complexes with ammonia, followed by HOX; the charge-transfer-type contribution has been demonstrated to dominate the halogen bonding in these complexes. For the complexes involving carbon-bound halogen molecules, our calculations clearly indicate that electrostatic interactions are mainly responsible for their binding energies. Whereas the halogen-bond strength is significantly enhanced by progressive fluorine substitution, the substitution of a hydrogen atom by a methyl group in the CH(3)X...NH(3) complex weakened the halogen bonding. Moreover, remote substituent effects have also been noted in the complexes of halobenzenes with different para substituents. The influence of the hybridization state of the carbon atom bonded to the halogen atom has also been examined and the results reveal that halogen-bond strengths decrease in the order HC triple bond CX > H(2)C=CHX approximately O=CHX approximately C(6)H(5)X > CH(3)X. In addition, several excellent linear correlations have been established between the interaction energies and both the amount of charge transfer and the electrostatic potentials corresponding to an electron density of 0.002 au along the R-X axis; these correlations provide good models with which to evaluate the electron-accepting abilities of the covalently bonded halogen atoms. Finally, some positively charged halogen-bonded systems have been investigated and the effect of the charge has been discussed.     

Analysis of competing bonding parameters. Part 2. The structure of halosilanes and halogermanes (MH4−nXn,n=1–4; M=Si,Ge; X=F,Cl, Br)

A qualitative rationalization of bonding patterns in halosilanes and halogermanes (MH_4−nX_n, n=1–4; M=Si, Ge; X=F, Cl, Br) is presented. Geometrical and bonding properties in these molecules are discussed on the basis of ab initio molecular orbital calculations employing the natural bond orbital population analysis. The results have been compared with data derived previously for halomethanes. Differences in the n-dependence of the M–Cl and M–Br bond lengths for M=C, Si, Ge are explained by a significant reduction in the closed-shell repulsion between the halide atoms. As M gets larger, a continuous decrease in the X–M–X bond angle is observed. Small bond angles (for n=2, 3) are favoured by the p-rich M orbitals in the M–X bonds. They are opposed, however, by the X⋯X repulsion. As M gets larger, the X⋯X separation for a given bond angle increases. A reduction in the X–M–X bond angle is therefore accomplished without overcompensation due to the X⋯X repulsion energy. The variation in the charge density at M as a function of n has been rationalized by differences in the electronegativity of the terminal atoms H and X. Dipole moments have been computed for the molecules in the series. As in the fluoromethanes, a maximum in the dipole moments at n=2 is explained by a combination of geometric and electronic properties unique to the fluoro-compounds. These are, an n-independent charge density at the F sites and a significant decrease in the M–F bond distance as n increases.     

Electronic structure and metal-metal interactions in trinuclear face-shared [M3X12]3- (M = Mo, W; X = F, Cl, Br, I) systems

The molecular and electronic structures of trinuclear face-shared [M3X12]3-species of Mo (X = F, Cl, Br, I) and W (X = Cl), containing linear chains of metal atoms, have been investigated using density functional theory. The possibility of variations in structure and bonding has been explored by considering both symmetric (D3d) and unsymmetric (C3v) forms, the latter having one long and one short metal-metal distance. Analysis of the bonding in the structurally characterized [Mo3I12]3- trimer reveals that the metal-metal interaction qualitatively corresponds to a two-electron three-center sigma bond between the Mo atoms and, consequently, a formal Mo-Mo bond order of 0.5. However, the calculated spin densities suggest that the electrons in the metal-metal sigma bond are not fully decoupled and therefore participate in the antiferromagnetic interactions of the metal cluster. Although the same observation applies to [Mo3X12]3- (X = Br, Cl, F) and [W3Cl12]3-, both the spin densities and shorter distances between the metal atoms indicate that the metal-metal interaction is stronger in these systems. The broken-symmetry approach combined with spin projection has been used to determine the energy of the low-lying spin multiplets arising from the magnetic coupling between the metal centers. Either the symmetric and unsymmetric S = 3/2 state is predicted to be the ground state for all five systems. For [Mo3X12]3- (X = Cl, Br, I), the symmetric form is more stable but the unsymmetric structure, where two metal centers are involved in a metal-metal triple bond while the third center is decoupled, lies close in energy and is thermally accessible. Consequently, at room temperature, interconversion between the two energetically equivalent configurations of the unsymmetric form should result in an averaged structure that is symmetric. This prediction is consistent with the reported structure of [Mo3I12]3-, which, although symmetric, indicates significant movement of the central Mo atom toward the terminal Mo atoms on either side. In contrast, unsymmetric structures with a triple bond between two metal centers are predicted for [Mo3F2]3- and [W3C12]3-, as the symmetric structure lies too high in energy to be thermally accessible.     

Molecular orbital investigation of the protonated H2X2AlNHn(CH3)3-n+ (X = F, Cl, and Br; n = 0-3) complexes

Structures of protonated alane-Lewis base donor-acceptor complexes H2X2AlNHn(CH3)(3-n)+ (X = F, Cl, and Br; n = 0-3) as well as their neutral parents were investigated. All the monocations H2X2AlNHn(CH3)(3-n)+ are Al-H protonated involving hypercoordinated alane with a three-center two-electron bond and adopt the C(s) symmetry arrangement. The energetic results show that the protonated alane-Lewis complexes are more stable than the neutral ones. They also show that this stability decreases on descending in the corresponding periodic table column from fluorine to bromine atoms. The calculated protonation energies of HX2AlNHn(CH3)(3-n) to form H2X2AlNHn(CH3)(3-n)+ were found to be highly exothermic. The possible dissociation of the cations H2X2AlNHn(CH3)(3-n)+ into X2AlNHn(CH3)(3-n)+ and molecular H2 is calculated to be endothermic.     

Bonding analysis using localized relativistic orbitals: water, the ultrarelativistic case and the heavy homologues H2X (X = Te, Po, eka-Po)

We report the implementation of Pipek-Mezey [J. Chem. Phys. 90, 4916 (1989)] localization of molecular orbitals in the framework of a four-component relativistic molecular electronic structure theory. We have used an exponential parametrization of orbital rotations which allows the use of unconstrained optimization techniques. We demonstrate the strong basis set dependence of the Pipek-Mezey localization criterion and how it can be eliminated. We have employed localization in conjunction with projection analysis to study the bonding in the water molecule and its heavy homologues. We demonstrate that in localized orbitals the repulsion between hydrogens in the water molecule is dominated by electrostatic rather than exchange interactions and that freezing the oxygen 2s orbital blocks polarization of this orbital rather than hybridization. We also point out that the bond angle of the water molecule cannot be rationalized from the potential energy alone due to the force term of the molecular virial theorem that comes into play at nonequilibrium geometries and which turns out to be crucial in order to correctly reproduce the minimum of the total energy surface. In order to rapidly assess the possible relativistic effects we have carried out the geometry optimizations of the water molecule at various reduced speed of light with and without spin-orbit interaction. At intermediate speeds, the bond angle is reduced to around 90 degrees , as is known experimentally for H(2)S and heavier homologues, although our model of ultrarelativistic water by construction does not allow any contribution from d orbitals to bonding. At low speeds of light the water molecule becomes linear which is in apparent agreement with the valence shell electron pair repulsion (VSEPR) model since the oxygen 2s12 and 2p12 orbitals both become chemically inert. However, we show that linearity is brought about by the relativistic stabilization of the (n + 1)s orbital, the same mechanism that leads to an electron affinity for eka-radon. Actual calculations on the series H2X (X = Te, Po, eka-Po) show the spin-orbit effects for the heavier species that can be rationalized by the interplay between SO-induced bond lengthening and charge transfer. Finally, we demonstrate that although both the VSEPR and the more recent ligand close packing model are presented as orbital-free models, they are sensitive to orbital input. For the series H2X (X = O, S, Se, Te) the ligand radius of the hydrogen can be obtained from the covalent radius of the central atom by the simple relation r(lig)(H) = 0.67r(cov)(X) + 27 (in picometers).     

Structures and stabilities of group 13 adducts [(NHC)(EX3)] and [(NHC)2(E2X(n))] (E=B to In; X=H, Cl; n=4, 2, 0; NHC=N-heterocyclic carbene) and the search for hydrogen storage systems: a theoretical study

Quantum chemical calculations using density functional theory at the BP86/TZVPP level and ab initio calculations at the SCS-MP2/TZVPP level have been carried out for the group 13 complexes [(NHC)(EX(3))] and [(NHC)(2)(E(2)X(n))] (E=B to In; X=H, Cl; n=4, 2, 0; NHC=N-heterocyclic carbene). The monodentate Lewis acids EX(3) and the bidentate Lewis acids E(2) X(n) bind N-heterocyclic carbenes rather strongly in donor-acceptor complexes [(NHC)(EX(3))] and [(NHC)(2)(E(2)X(n))]. The equilibrium structures of the bidentate complexes depend on the electronic reference state of E(2)X(n), which may vary for different atoms E and X. All complexes [(NHC)(2)(E(2)X(4))] possess C(s) symmetry in which the NHC ligands bind in a trans conformation to the group 13 atoms E. The complexes [(NHC)(2)(E(2)H(2))] with E=B, Al, Ga have also C(s) symmetry with a trans arrangement of the NHC ligands and a planar CE(H)E(H)C moiety that has a E=E π bond. In contrast, the indium complex [(NHC)(2)(In(2) H(2))] has C(i) symmetry with pyramidal-coordinated In atoms in which the hydrogen atoms are twisted above and below the CInInC plane. The latter C(i) form is calculated for all chloride systems [(NHC)(2)(E(2)Cl(2))], but the boron complex [(NHC)(2)(B(2)Cl(2))] deviates only slightly from C(s) symmetry. The B(2) fragment in the linear coordinated complex [(NHC)(2)(B(2))] has a highly excited (3)(1)Σ(g)(-) reference state, which gives an effective B≡B triple bond with a very short interatomic distance. The heavier homologues [(NHC)(2)(E(2))] (E=Al to In) exhibit a anti-periplanar arrangement of the NHC ligands in which the E(2) fragments have a (1)(1) Δ(g) reference state and an E=E double bond. The calculated energies suggest that the dihydrogen release from the complexes [(NHC)(EH(3))] and [(NHC)(2)(E(2)H(n))] becomes energetically more favourable when atom E becomes heavier. The indium complexes should therefore be the best candidates of the investigated series for hydrogen-storage systems that could potentially deliver dihydrogen at close to ambient temperature. The hydrogenation reaction of the dimeric magnesium(I) compound [LMgMgL] (L=β-diketiminate) with [(NHC)(EH(3))] becomes increasingly exothermic with the trend B相似文献   

Computational study on the characteristics of the interaction in naphthalene...(H2X)n=1,2 (X = O,S) clusters

The characteristics of the interaction between the pi cloud of naphthalene and up to two H2O or H2S molecules were studied. Calculations show that clusters formed by naphthalene and one H2O or H2S molecule have similar geometric features, and also present similar interaction energies. Our best estimates for the interaction energy amount to -2.95 and -2.92 kcal/mol for H2O and H2S, respectively, as obtained with the CCSD(T) method. Calculations at the MP2 level employing large basis sets should be avoided because they produce highly overestimated interaction energies, especially for hydrogen sulfide complexes. The MPWB1K functional, however, provides values pretty similar to those obtained with the CCSD(T) method. Although the magnitude of the interaction is similar with both H2X molecules, its nature is somewhat different: the H2O complex presents electrostatic and dispersion contributions of similar magnitude, whereas for H2S the interaction is dominated by dispersion. In clusters containing two H2X molecules several minima were characterized. In water clusters, the total interaction energy is dominated by the presence of a O-H...O hydrogen bond and, as a consequence, structures where this contact is present are the most stable. However, clusters containing H2S show structures with no interaction between H2S moieties which are as stable as the hydrogen bonded ones, because they allow an optimal H2S...naphthalene interaction, which is stronger than the S-H...S contact. Therefore it is possible that in larger polycycles hydrogen sulfide molecules will be spread onto the surface maximizing S-H...pi interactions rather than aggregated, forming H2S clusters.     

Theoretical study of (CH...C)- hydrogen bonds in CH(4-n)X(n) (X = F, Cl; n = 0, 1, 2) systems complexed with their homoconjugate and heteroconjugate carbanions

This work deals with a theoretical study of the (CH...C)- hydrogen bonds in CH4, CH3X, and CH2X2 (X = F, Cl) complexed with their homoconjugate and heteroconjugate carbanions. The properties of the complexes are calculated with the B3LYP method using the 6-311++G(d,p) or 6-311++G(2df,2p) basis sets. The deprotonation enthalpies (DPE) of the CH bond or the proton affinities of the carbanions (PA(C-) are calculated as well. All the systems with the exception of the CH4...CHCl2(-) one are characterized by a double minimum potential. In some of the complexes, the (CH(b)...C)- hydrogen bond is linear. In other systems, such as CH3F...CH2F- and CH3F...CHF2(-), there is a large departure from linearity, the systems being stabilized by electrostatic interactions between the nonbonded H of the neutral molecule and the F atom of the carbanion. In the transition state, the (CH(b)...C)- bond is linear, and there is a large contraction of the intermolecular C...C distance. The binding energies vary within a large range, from -1.4 to -11.1 kcal mol(-1) for the stable complexes and -8.6 to -44.1 kcal mol(-1) for the metastable complexes. The energy barriers to proton transfer are between 5 and 20 kcal mol(-1) for the heteroconjugate systems and between 3.8 and 8.3 kcal mol(-1) for the homoconjugate systems. The binding energies of the linear complexes depend exponentially on 1.5DPE - PA(C-), showing that the proton donor is more important than the proton acceptor in determining hydrogen bond strength. The NBO analysis indicates an important electronic reorganization in the two partners. The elongations of the CH bond resulting from the interaction with the carbanion depend on the occupation of the sigma*(CH(b)) antibonding orbitals and on the hybridization of the C bonded to H(b). The frequency shifts of the nu(CH)(A1) stretching vibration range between 15 and 1150 cm(-1). They are linearly correlated to the elongation of the CH(b) bond.     

Use of oxidation-state differences and molecular orbitals to interpret bonding in the series ONXYZ (X, Y, Z=H, F, Cl), HNNX3, HNNX2Y, and HNNXY2 (X, Y=H, F) and OCX3-, OCX2Y-, and OCXY2- (X, Y=H, F)

Bonding in the series ONXYZ (X, Y, Z=H, F, Cl), HNNX3, HNNX2Y, HNNXY2 (X, Y=H, F), and OCX3-, OCX2Y-, OCXY2- (X, Y=H, F) shows evidence of a significant ionic contribution modifying the underlying covalent bonding. Increased ionic character can be correlated with oxidation-state differences between the bound atoms and is expressed in terms of shorter bond lengths. All members of the series, with the exception of ONH3, HNNH3, and OCH3-, possess a multiple O-N, N-N, or C-O bond modified by the ionic character of that bond. The O-N, N-N, and O-C single bonds in ONH3, HNNH3, and OCH3-, respectively, show some variation in length relative to "typical" single bonds of these types due to differences in ionic character. The two highest-occupied molecular orbitals in the ONXYZ or OCXYZ- (X, Y, Z=H, F) series which are piNO or piCO (when X=Y=Z=H) exhibit a distinct shift in their nodal plane as hydrogen is replaced by fluorine. The nodal plane moves from a location between the oxygen and the nitrogen or carbon to between the nitrogen or carbon and the fluorines impacting on the nature and length of the bonds joining these atoms. The pattern of N-F and C-F bond lengths in the series, ONH3-ONF3 and OCH3--OCF3-, respectively, lends support to the idea of resonance structures of the form ONXY+ F- or OCXY F- (where X, Y=H, F).     

Accounting for the differences in the structures and relative energies of the highly homoatomic np pi-np pi (n > or = 3)-bonded S2I4 2+, the Se-I pi-bonded Se2I4 2+, and their higher-energy isomers by AIM, MO, NBO, and VB methodologies

The bonding in the highly homoatomic np pi-np pi (n > or = 3)-bonded S2I42+ (three sigma + two pi bonds), the Se-I pi-bonded Se2I42+ (four sigma + one pi bonds), and their higher-energy isomers have been studied using modern DFT and ab initio calculations and theoretical analysis methods: atoms in molecules (AIM), molecular orbital (MO), natural bond orbital (NBO), and valence bond (VB) analyses, giving their relative energies, theoretical bond orders, and atomic charges. The aim of this work was to seek theory-based answers to four main questions: (1) Are the previously proposed simple pi*-pi* bonding models valid for S2I42+ and Se2I42+? (2) What accounts for the difference in the structures of S2I42+ and Se2I42+? (3) Why are the classically bonded isolobal P2I4 and As2I4 structures not adopted? (4) Is the high experimentally observed S-S bond order supported by theoretical bond orders, and how does it relate to high bond orders between other heavier main group elements? The AIM analysis confirmed the high bond orders and established that the weak bonds observed in S2I42+ and Se2I42+ are real and the bonding in these cations is covalent in nature. The full MO analysis confirmed that S2I42+ contains three sigma and two pi bonds, that the positive charge is essentially equally distributed over all atoms, that the bonding between S2 and two I2+ units in S2I42+ is best described by two mutually perpendicular 4c2e pi*-pi* bonds, and that in Se2I42+, two SeI2+ moieties are joined by a 6c2e pi*-pi* bond, both in agreement with previously suggested models. The VB treatment provided a complementary approach to MO analysis and provided insight how the formation of the weak bonds affects the other bonds. The NBO analysis and the calculated AIM charges showed that the minimization of the electrostatic repulsion between EI2+ units (E = S, Se) and the delocalization of the positive charge are the main factors that explain why the nonclassical structures are favored for S2I42+ and Se2I42+. The difference in the structures of S2I42+ and Se2I42+ is related to the high strength of the S-S pi bond compared to the weak S-I sigma bond and the additional stabilization from increased delocalization of positive charge in the structure of S2I42+ compared to the structure of Se2I42+. The investigation of the E2X42+ series (E = S, Se, Te; X = Cl, Br, I) revealed that only S2I42+ adopts the highly np pi-np pi (n > or = 3)-bonded structure, while all other dications favor the pi-bonded Se2I42+ structure. Theoretical bond order calculations for S2I42+ confirm the previously presented experimentally based bond orders for S-S (2.1-2.3) and I-I (1.3-1.5) bonds. The S-S bond is determined to have the highest reported S-S bond order in an isolated compound and has a bond order that is either similar to or slightly less than the Si-Si bond order in the proposed triply bonded [(Me3Si)2CH]2(iPr)SiSi triple bond SiSi(iPr)[CH(SiMe3)2]2 depending on the definition of bond orders used.     

Theoretical Density Functional Theory insights into the nature of chalcogen bonding between CX2 (X = S,Se, Te) and diazine from monomer to supramolecular complexes

Chalcogen bonding is a noncovalent interaction, highly similar to halogen and hydrogen bonding, occurring between a chalcogen atom and a nucleophilic region. Two density functional theory (DFT) approaches B3LY-D3 and B97-D3 were performed on a series of complexes formed between CX₂ (X = S, Se, Te) and diazine (pyridazine, pyrimidine and pyrazine). Chalcogen atoms prefer interacting with the lone pair of a nitrogen atom rather than with the π-cloud of an aromatic ring. CTe₂ and CSe₂ form a stronger chalcogen bond than CS₂. The electrostatic potential of CX₂ (X = S, Se and Te) reveals the presence of two equivalent σ-holes, one on each chalcogen atom. These CX₂ molecules interact with diazine giving rise to supramolecular interactions. Wiberg bond index and second-order perturbation theory analysis in NBO were performed to better understand the nature of the chalcogen bond interaction.     

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.     

Cationic P-S-X cages (X=Br, I)

The first condensed-phase preparation of ternary P-Ch-X cations (Ch=O-Te, X=F-I) is reported: [P5S3X2]+, [P5S2X2]+, and [P4S4X]+ (X=Br, I). [P5S3X2]+ is formed from the reaction of the Ag+/PX3 reagent with P4S3. The [P5S3X2]+ ions have a structure that is related to P4S5 by replacing P=S by P+--X and S in the four-membered ring by P(X). We provide evidence that the active ingredient of the Ag+/PX3 reagent is the (H2CCl2)Ag-X-PX2+ cation. The latter likely reacts with the HOMO of P4S3 in a concerted HOMO-LUMO addition to give the P5S3X2+ ion as the first species visible in situ in the low-temperature 31P NMR spectrum. The [P5S3X2]+ ions are metastable at -78 degrees C and disproportionate at slightly higher temperatures to give [P5S2X2]+ and [P4S4X]+, probably with the extrusion of 1/n (PX)n (X=Br, I). All six new cage compounds have been characterized by multinuclear NMR spectroscopy and, in part, by IR or Raman spectroscopy. The [P5S2X2]+ salts have a nortricyclane skeleton and were also characterized by X-ray crystallography. The structure of the [P4S4X]+ ion is related to that of P4S5 in that the exo-cage P=S bond is replaced by an isoelectronic P+--X moiety.     

Cooperative effects and strengths of hydrogen bonds in open-chain cis-triaziridine clusters (n = 2-8): a DFT investigation

We employ DFT/B3LYP method to investigate linear open-chain clusters (n = 2-8) of the cis-triaziridine molecule that is a candidate molecule for high energy density materials (HEDM). Our calculations indicate that the pervasive phenomena of cooperative effects are observed in the clusters of n = 3-8, which are reflected in changes in lengths of N...H hydrogen bonds, stretching frequencies, and intensities of N-H bonds, dipole moments, and charge transfers as cluster size increases. The n(N) --> sigma*(N-H) interactions, i.e., the charge transfers from lone pairs (n(N)) of the N atoms into antibonds (sigma*) of the N-H bonds acting as H-donors, can be used to explain the observed cooperative phenomena. The approaches based upon natural bond orbital (NBO) method and theory of atoms in molecule (AIM) to evaluating N...H strengths are found to be equivalent. In the process of N...H bonding, cooperative nature of n(N) --> sigma*(N-H) interactions promotes formation of stronger N...H bonds as reflected in increases in the capacities of cis-triaziridine clusters to concentrate electrons at the bond critical points of N...H bonds. The calculated nonadditive energies also show that the cooperative effects due to n(N) --> sigma*(N-H) interactions indeed provide additional stabilities for the clusters.     

Halogen-hydride interaction between Z-X (Z = CN, NC; X = F, Cl, Br) and H-Mg-Y (Y = H, F, Cl, Br, CH3)

Halogen-hydride interactions between Z-X (Z = CN, NC and X = F, Cl, Br) as halogen donor and H-Mg-Y (Y = H, F, Cl, Br, CH(3)) as electron donor have been investigated through the use of Becke three-parameter hybrid exchange with Lee-Yang-Parr correlation (B3LYP), second-order M?ller-Plesset perturbation theory (MP2), and coupled-cluster single and double excitation (with triple excitations) [CCSD(T)] approaches. Geometry changes during the halogen-hydride interaction are accompanied by a mutual polarization of both partners with some charge transfer occurring from the electron donor subunit. Interaction energies computed at MP2 level vary from -1.23 to -2.99 kJ/mol for Z-F···H-Mg-Y complexes, indicating that the fluorine interactions are relatively very weak but not negligible. Instead, for chlorine- and bromine-containing complexes the interaction energies span from -5.78 to a maximum of -26.42 kJ/mol, which intimate that the interactions are comparable to conventional hydrogen bonding. Moreover, the calculated interaction energy was found to increase in magnitude with increasing positive electrostatic potential on the extension of Z-X bond. Analysis of geometric, vibrational frequency shift and the interaction energies indicates that, depending on the halogen, CN-X···H interactions are about 1.3-2.0 times stronger than NC-X···H interactions in which the halogen bonds to carbon. We also identified a clear dependence of the halogen-hydride bond strength on the electron-donating or -withdrawing effect of the substituent in the H-Mg-Y subunits. Furthermore, the electronic and structural properties of the resulting complexes have been unveiled by means of the atoms in molecules (AIM) and natural bond orbital (NBO) analyses. Finally, several correlative relationships between interaction energies and various properties such as binding distance, frequency shift, molecular electrostatic potential, and intermolecular density at bond critical point have been checked for all studied systems.     

Study of the interaction in clusters formed by phenol and CH3X (X=CN,F,Cl) molecules

The characteristics of the interaction between phenol and acetonitrile, methyl fluoride and methyl chloride were studied. The most stable structures for clusters containing one or two CH3X molecules and one phenol moiety were located by means of ab initio and density functional theory calculations. Phenol-acetonitrile dimer presents two almost equally stable structures; one of them is a typical linearly hydrogen bonded minimum, whereas in the other one, a C-H...pi contact is established accompanied by a distorted O-H...N hydrogen bond. Although the latter minimum presents the larger interaction energy, deformation effects favor the formation of the linear hydrogen bonded one. In complexes with methyl fluoride and methyl chloride, this arrangement is the most stable structure and no linear hydrogen bonded structures were located. Our best estimates for the interaction energies amount to -27.8, -21.6, and -19.7 kJ/mol for clusters of phenol with acetonitrile, methyl fluoride, and methyl chloride, respectively. The main contribution to the stabilization of these clusters is of electrostatic nature, although in structures where a C-H...pi contact is present, the dispersion contribution is also significant. In clusters formed by phenol and two CH3X units, the most stable arrangement corresponds to a head to tail disposal with O-H...X, C-H...X, and C-H...pi contacts forming a cycle. Only for this type of arrangement, three body effects are non-negligible even though they constitute a minor effect. The results also indicate that interactions with methyl fluoride and methyl chloride are of similar intensity, although weaker than with acetonitrile. Significant frequency shifts are predicted for the O-H stretching, which increase when increasing the number of CH3X molecules.     

Redox energetics of hypercloso boron hydrides B(n)H(n) (n = 6-13) and B12X12 (X = F, Cl, OH, and CH3)

The reduction potentials (E°(Red) versus SHE) of hypercloso boron hydrides B(n)H(n) (n = 6-13) and B(12)X(12) (X = F, Cl, OH, and CH(3)) in water have been computed using the Conductor-like Polarizable Continuum Model (CPCM) and the Solvation Model Density (SMD) method for solvation modeling. The B3LYP/aug-cc-pvtz and M06-2X/aug-cc-pvtz as well as G4 level of theory were applied to determine the free energies of the first and second electron attachment (ΔG(E.A.)) to boron clusters. The solvation free energies (ΔG(solv)) greatly depend on the choice of the cavity set (UAKS, Pauling, or SMD) while the dependence on the choice of exchange/correlation functional is modest. The SMD cavity set gives the largest ΔΔG(solv) for B(n)H(n)(0/-) and B(n)H(n)(-/2-) while the UAKS cavity set gives the smallest ΔΔG(solv) value. The E°(Red) of B(n)H(n)(-/2-) (n = 6-12) with the G4/M06-2X(Pauling) (energy/solvation(cavity)) combination agrees within 0.2 V of experimental values. The experimental oxidative stability (E(1/2)) of B(n)X(n)(2-) (X = F, Cl, OH, and CH(3)) is usually located between the values predicted using the B3LYP and M06-2X functionals. The disproportionation free energies (ΔG(dpro)) of 2B(n)H(n)(-) → B(n)H(n) + B(n)H(n)(2-) reveal that the stabilities of B(n)H(n)(-) (n = 6-13) to disproportionation decrease in the order B(8)H(8)(-) > B(9)H(9)(-) > B(11)H(11)(-) > B(10)H(10)(-). The spin densities in B(12)X(12)(-) (X = F, Cl, OH, and CH(3)) tend to delocalize on the boron atoms rather than on the exterior functional groups. The partitioning of ΔG(solv)(B(n)H(n)(2-)) over spheres allows a rationalization of the nonlinear correlation between ΔG(E.A.) and E°(Red) for B(6)H(6)(-/2-), B(11)H(11)(-/2-), and B(13)H(13)(-/2-).     

Localized molecular orbital studies on B_4C1_4,1,5-C_2B_3H_5 and B_nH_n~(2-)(n = 6-10, 12)

The wave functions, level energies and Mulliken population analysis of localized molecular orbitals (LMO's) for B4Cl4, 1,5-C2B3H5 and the closo-BnHn2- (n = 6-10, 12) are calculated by using the Edmiston-Reudenberg energy localization scheme under the CNDO/2 approximation in order to reveal the nature of quasi-aromaticity of the closo-BnHn2- (n > 5). It has been found that all the B-H or B-Cl LMO's are highly localized between the B and H (or Cl) atoms, corresponding to B-H or B-Cl o-bond, while the Bn framework bonding is formed mainly by the three-centered two-electron B-B-B bonds on the polyhedral faces. In the cases of B4Cl4 and 1,5-C2B3H5, these three-centered B-B-B bonds just fill their polyhedral faces; however, for the framework bonding of the closo-BnHn2- (n > 5), the valence electron deficiency leads to the delocalization of their three-centered B-B-B bonds, and as delocalizability of this three-centered B-B-B bond increases, some three-centered B-B-B bonds are further delocalized to become a f     

Geometric, electronic, and bonding properties of AuNM (N = 1-7, M = Ni, Pd, Pt) clusters

Employing first-principles methods, based on density functional theory, we report the ground state geometric and electronic structures of gold clusters doped with platinum group atoms, Au(N)M (N = 1-7, M = Ni, Pd, Pt). The stability and electronic properties of Ni-doped gold clusters are similar to that of pure gold clusters with an enhancement of bond strength. Due to the strong d-d or s-d interplay between impurities and gold atoms originating in the relativistic effects and unique properties of dopant delocalized s-electrons in Pd- and Pt-doped gold clusters, the dopant atoms markedly change the geometric and electronic properties of gold clusters, and stronger bond energies are found in Pt-doped clusters. The Mulliken populations analysis of impurities and detailed decompositions of bond energies as well as a variety of density of states of the most stable dopant gold clusters are given to understand the different effects of individual dopant atom on bonding and electronic properties of dopant gold clusters. From the electronic properties of dopant gold clusters, the different chemical reactivity toward O(2), CO, or NO molecule is predicted in transition metal-doped gold clusters compared to pure gold clusters.