Electron density and energy decomposition analysis in hydrogen-bonded complexes of azabenzenes with water, acetamide, and thioacetamide

Ab initio and density functional theoretical studies on hydrogen-bonded complexes of azabenzenes with water, acetamide, and thioacetamide have been carried out to explore the controversy involved in the relative order of their stability in a systematic way. The interaction energies of these complexes have been analyzed using the Morokuma energy decomposition method, and the nature of the various hydrogen bonds formed has been investigated through topological aspects using Bader's atom in a molecule (AIM) theory. Morokuma energy decomposition analysis reveals that the major contributions to the energetics are from the polarization (PL) and charge transfer (CT) energies. From the calculated topological results, excellent linear correlation is shown to exist between the hydrogen-bond length, electron density [rho(r)], and its Laplacian [nabla(2)rho(r)] at the bond critical points for all the complexes considered.     

Bond length and local energy density property connections for non-transition-metal oxide-bonded interactions

For a variety of molecules and earth materials, the theoretical local kinetic energy density, G(r(c)), increases and the local potential energy density, V(r(c)), decreases as the M-O bond lengths (M = first- and second-row metal atoms bonded to O) decrease and the electron density, rho(r(c)), accumulates at the bond critical points, r(c). Despite the claim that the local kinetic energy density per electronic charge, G(r(c))/rho(r(c)), classifies bonded interactions as shared interactions when less than unity and closed-shell when greater, the ratio was found to increase from 0.5 to 2.5 au as the local electronic energy density, H(r(c)) = G(r(c)) + V(r(c)), decreases and becomes progressively more negative. The ratio appears to be a measure of the character of a given M-O bonded interaction, the greater the ratio, the larger the value of rho(r(c)), the smaller the coordination number of the M atom and the more shared the bonded interaction. H(r(c))/rho(r(c)) versus G(r(c))/rho(r(c)) scatter diagrams categorize the M-O bonded interactions into domains with the local electronic energy density per electron charge, H(r(c))/rho(r(c)), tending to decrease as the electronegativity differences for the bonded pairs of atoms decrease. The values of G(r(c)) and V(r(c)), estimated with a gradient-corrected electron gas theory expression and the local virial theorem, are in good agreement with theoretical values, particularly for the bonded interactions involving second-row M atoms. The agreement is poorer for shared C-O and N-O bonded interactions.     

Experimental and theoretical charge density distribution of the colossal magnetoresistive transition metal sulfide FeCr2S4

The total charge density distribution rho(r) of the colossal magnetoresistive transition metal sulfide FeCr(2)S(4) was evaluated through a multipole formalism from a set of structure factors obtained both experimentally, by means of single crystal high-quality x-ray diffraction data collected at T=23 K, and theoretically, with an extended-basis unrestricted Hartree-Fock periodic calculation on the experimental geometry. A full topological analysis, followed by the calculation of local energy density values and net atomic charges, was performed using the quantum theory of atoms in molecules. The experimental and theoretical results were compared. Good agreement was found for the topological properties of the system, as well as for the atomic net charges and the nature of the chemical bonds. An analysis of the electron density rho(r), its Laplacian nabla(2)[rho(r)], and the total energy density H(r) at the bond critical points was employed to classify all the interactions that resulted as predominantly closed shell (ionic) in nature. The topological indicators of the bonded interactions for Fe are distinct from those for Cr. The Fe-S bond distances were found to be 0.145 A shorter than the ideal values computed on the basis of Shannon's crystal radii, much shorter than the Cr-S distances with respect to their ideal Shannon lengths. Concomitantly, rho(r) and |H(r)| at the bond critical points are greater for Fe-S interactions, indicating that the local concentration of charge density in the internuclear region is larger for the tetrahedrally coordinated iron than for the octahedrally coordinated chromium. The isosurface in the real space for nabla(2)[rho(r)]=0 was plotted for both iron and chromium, pointing out the local zones of valence shell charge concentration and relating them to the partial d-orbital occupancy of the two transition metal atoms.     

Understanding the electronic structure, reactivity, and hydrogen bonding for a 1,2-diphosphonium dication

The experimental charge density for hexamethyldiphosphonium ditriflate has been determined from low-temperature high-resolution X-ray diffraction data. These results have been compared with theoretically calculated values for the isolated gas-phase compound. Analysis of the topological and atomic basin properties has provided insight into the exact nature of the P-P bond in both the crystalline and the gas-phase structures. The rho(b)(r) and nabla2rho(b)(r) values highlight the covalent nature of the P-P bond, while the atomic charges indicate a localization of the positive charges on the two phosphorus atoms. This seems to indicate that a covalent bond is formed despite a strong electrostatic repulsion between these two heteroatoms. The topological properties and electrostatic potentials have also been shown to provide significant insight into the chemical reactivity of the title compound. A topological analysis of P2Me4, P2Me5(+), and P2Me6(+2) species has provided information about the progression of the P-P bond in the synthesis of the title compound. An investigation of the different hydrogen-bonding networks present in the crystalline and gas-phase structures, along with their affect on the electronic structure of the title compound has also been investigated. This has all led to significant new insight into the electronic structure, reactivity, and weak hydrogen bonding in prototypical 1,2-diphosphonium dications.     

On the perturbation of the H-bonding interaction in ethylene glycol clusters upon hydration

Ab initio and density functional methods have been employed to study the structure, stability, and spectral properties of various ethylene glycol (EG(m)) and ethylene glycol-water (EG(m)W(n)) (m = 1-3, n = 1-4) clusters. The effective fragment potential (EFP) approach was used to explore various possible EG(m)W(n) clusters. Calculated interaction energies of EG(m)W(n) clusters confirm that the hydrogen-bonding interaction between EG molecules is perturbed by the presence of water molecules and vice versa. Further, energy decomposition analysis shows that both electrostatic and polarization interactions predominantly contribute to the stability of these clusters. It was found from the same analysis that ethylene glycol-water interaction is predominant over the ethylene glycol-ethylene glycol and water-water interactions. Overall, the results clearly illustrate that the presence of water disrupts the ethylene glycol-ethylene glycol hydrogen bonds.     

Synthesis and characterization of stable hypervalent carbon compounds (10-C-5) bearing a 2,6-bis(p-substituted phenyloxymethyl)benzene ligand

X-ray analysis of bis(p-fluorophenyl)methyl cation bearing a 2,6-bis(p-tolyloxymethyl)benzene ligand showed a symmetrical structure (10-C-5) where the two C-O distances are identical, although the distance (2.690(4) A) is longer than those (2.43(1) and 2.45(1) A) of 1,8-dimethoxy-9-dimethoxymethylanthracene monocation, which was recently reported by us. However, X-ray analysis of the more stable aromatic xanthylium cation with the same benzene ligand showed the tetracoordinate carbon structure where only one of the two oxygen ligands is coordinated with the central carbon atom. These results clearly indicate that the carbocations (10-C-5) bearing the sterically flexible benzene ligand were quite sensitive to the electronic effect on the central carbon atom. The electron distribution analysis by accurate X-ray measurements and the density functional calculation on the initially mentioned bis(p-fluorophenyl)methyl cation clearly show that the central carbon atom and the two oxygen atoms are bonded even if the bond is weak and ionic based on the small value of the electron density (rho(r)) and the small positive Laplacian value (nabla(2)rho(r)) at the bond critical points.     

Spatial shape of electron delocalization: structure of the laplacian of the negative exchange-correlation density

The Laplacian of the negative exchange-correlation density (with respect to coordinate r(2)), nabla<(r)2>(2)[-Gamma(sigma1)(sigma2)(XC) (r(1),r(2))] = nabla(r)2(2)X(sigma1)(sigma2)(r(1),r(2)), is proposed as an instrument for the analysis of electron delocalization in real space. It determines local concentrations in the amount of electrons that are delocalized from a reference point r(1) over space. Integration of the reference coordinate r(1) over an atomic basin Omega(n) gives the function nabla(2)X(sigma1)(sigma2)(Omega(n);r), which contains detailed information about the spatial shape of the delocalization that originates from an atom in a molecule. Its isosurface representations are richly structured and resemble molecular orbitals in their complexity and partly also in their shape. The sum over all nabla(2)X(sigma1)(sigma2)(Omega(n);r) functions of a molecule equals the Laplacian of the electron density nabla(2)rho(r), for which it provides a meaningful partitioning into atomic contributions.     

Inhibited phenol ionization in reverse micelles: confinement effect at the nanometer scale

We found that the absorption spectra of 2-acetylphenol (2-HAP), 4-acetylphenol (4-HAP), and p-nitrophenol (p-NPh) in water/sodium 1,4-bis(2-ethylhexyl)sulfosuccinate (AOT)/n-heptane reverse micelles (RMs) at various W(0) (W(0) = [H(2)O]/[surfactant]) values studied changed with time if (-)OH ions were present in the RM water pool. There is an evolution of ionized phenol (phenolate) bands to nonionized phenol absorption bands with time and this process is faster at low W(0) values and with phenols with higher bulk water pK(a) values. That is, in bulk water and at the hydroxide anion concentration used, only phenolate species are observed, whereas in AOT RMs at this fixed hydroxide anion concentration, ionized phenols convert into nonionized phenol species over time. Furthermore, we demonstrate that, independent of the (-)OH concentration used to prepare the AOT RMs, the nonionized phenols are the more stable species in the RM media. We explain our results by considering that strong hydrogen-bonding interactions between phenols and the AOT polar head groups result in the existence of only nonionized phenols at the AOT RM interface. The situation is quite different when the phenols are dissolved in cationic benzyl-n-hexadecyldimethylammonium chloride RMs. Therein, only phenolates species are present at the (-)OH concentrations used. The results clearly demonstrate that the classical definition of pH does not apply in a confined environment, such as in the interior of RMs and challenge the general idea that pH can be determined inside RMs.     

Experimental electron density analysis of Mn2(CO)10: metal-metal and metal-ligand bond characterization

The experimental electron density rho(r) of Mn2(CO)10 was determined by a multipole analysis of accurate X-ray diffraction data at 120 K. The quantum theory of atoms in molecules (QTAM) was applied to rho(r) and its Laplacian [symbol: see text] 2 rho(r). The QTAM analysis of rho(r) showed the presence of a bond critical point (rc); its associated bond path connects the two Mn atoms, but no cross interaction line was found between one manganese and the equatorial carbonyls of the other. The distribution of [symbol: see text] 2 rho(r) indicated "closed-shell" interactions for the metallic Mn-Mn bond and the dative Mn-CO bonds. The values of the topological parameters of the density at rc, rho(rc), [symbol: see text] 2 rho(rc), G(rc) (kinetic energy density), and V(rc) (potential energy density), characterize the bonds and are intermediate to those corresponding to typical ionic and covalent bonds.     

Topological properties of the electrostatic potential in weak and moderate N...H hydrogen bonds

The topological analyses of the electrostatic potential phi(r) and the electron density distribution rho(r) have been performed for a set of 20 neutral complexes with weak and moderate N...H bonds. In all cases, a zero flux surface of the electrostatic potential containing a saddle point analogous to the bond critical point of the electron density distribution is observed. These surfaces define an equivalent of the atomic basin of rho(r) for the electrostatic potential, which exhibits zero net charge and can be regarded as an electrostatically isolated region if its volume is finite. The phi(r) and rho(r) zero flux surfaces divide the hydrogen-bonding region in three parts, being the central one related to the electrostatic interaction between donor and acceptor. This central region exhibits a relative size of approximately 13-14% of the N...H distance dNH, it belongs to the outermost shell of the nitrogen and is mainly associated with its lone pair. Topological properties of both rho(r) and phi(r), as well as the electron kinetic (G) and potential (V) energy densities, show similar dependences with dNH at both bond critical points (phi-BCP and rho-BCP). Phenomenological proportionalities between the rho(r) curvatures and G and V are also found at the electrostatic potential critical point. The curvatures of the electrostatic potential, which are interpreted in terms of the electrostatic forces in the bonding region, present the same exponential dependency as the electron density distribution, to which they are related by Poisson's equation.     

Theoretical study of dihydrogen bonds between (XH)2, X = Li, Na, BeH, and MgH, and weak hydrogen bond donors (HCN, HNC, and HCCH)

The dihydrogen-bonded (DHB) complexes formed by (XH)2, with X = Li, Na, BeH, and MgH, with one, two, and four protonic molecules (HCN, HNC, and HCCH) have been studied. These complexes have been compared to those of the XH monomers with the same hydrogen bond donor molecules. The energetic results have been rationalized based on the electrostatic potential of the isolated hydridic systems. The electron density properties have been analyzed within the AIM methodology, both at the bond critical points and the integrated values at the atomic basins. Exponential relationships between several properties calculated at the bond critical points (rho,nabla2rho, lambdai, G, and V) and variation of integrated properties (energy, charge, and volume) vs the DHB distance have been obtained.     

Cooperative effects, strengths of hydrogen bonds, and intermolecular interactions in circular cis, trans-cyclotriazane clusters (n = 3-8)

Based on Becke's three parameter functional [J. Chem. Phys. 98, 5648 (1993)] of density functional theory (DFT) with the correlation of Lee-Yang-Parr [Phys. Rev. B 37, 785 (1988)] (DFT/B3LYP), the natural bond orbital (NBO) analysis, the Bader's theory of atoms in molecule (AIM), our calculations indicate that as cluster size (n) increases, the n-dependent cooperative changes in the lengths of the N...H H bonds (HBs) and N-H bonds, the N-H stretching frequencies and intensities, and the n(N)-->sigma*(N-H) charge transfers are observed to be pervasive in the circular cis, trans-cyclotriazane clusters (n = 3-8), which is very different from the linear cis, trans-cyclotriazane clusters reported in previous work. According to the NBO and AIM theories, the cooperativity of the intermolecular n(N)-->sigma*(N-H) interaction leads to the n-dependent N...H contractions. In this way, the stronger N...H bond is formed, as reflected in the increase in their rho(r(cp)) values. This increased electron density is translated into the improved capacity to concentrate electrons at the HB bond critical point (BCP), i.e., a higher potential energy V(r(cp)). On the other hand, stronger repulsion is also activated to counteract the contraction, which is reflected in the increased G(r(cp)) value that gives the tendency of the system to dilute electrons at the HB BCP. In terms of the three-body symmetry-adapted perturbation theory (three-body SAPT), the induction nonadditivity accounts for up to 97% of the nonadditive energy in the circular trimer. It can believed that the marked cooperativity of the n(N)-->sigma*(N-H) interactions is of nonadditive induction in nature. The N...H formation and nature of cooperativity in the circular clusters differ from those in the linear clusters that have been reported previously. According to the SAPT(DFT) method which is a combination of SAPT with the asymptotically corrected DFT, the cis, trans-cyclotriazane systems should contain remarkable dispersion interactions. However, the short-range dispersion cannot be reproduced thoroughly by DFT/B3LYP. A quantum cluster equilibrium model illustrates the neglected dispersion energies and the nonadditive energies can affect markedly the properties of the liquid consisting of the circular clusters.     

Bond length and the electron density at the bond critical point: X--X, Z--Z, and C--Z bonds (X = Li-F, Z = Na-Cl)

The aim was to investigate the relationship between the bond length and the electron density at the bond critical point in homonuclear X--X and Z--Z and heteronuclear C--Z bonds (X = Li-F, Z = Na-Cl). The d,rho(c) pairs were obtained from 472 target bonds in DFT-optimized (B3LYP/6-311+G(d,p)) small molecular species. These species were selected arbitrarily but with a view to maximize the range widths WR for each atom combination. It was found that (i) with one clear exception, the d(A - A) means (A = X or Z) correlate linearly with the bond lengths d(A(2)) of the respective diatomic molecules; (ii) the d(A - A) means correlate parabolically with n, the formal number of valence electrons in the atoms of the bond; and (iii) with increasing sample size N the ratio WR(rho(c))/WR(d) appears to converge toward a representation f [WR(rho(c))/WR(d)](N-->infinity) characteristic of A. Detailed analysis of the d,rho(c) relationship has shown that by and large simple power regression accounts best for the DFT data. The regression coefficients of d = arho(c) (-b) and rho(c) = alphad(-beta) (b, beta > 0) vary with n in a seemingly irregular manner but one that is consistent with simple chemical notions. The d(A(2)) can be approximated in terms of multilinear MO electron occupancies.     

The structure and chemical bonding in the N(2)-CuX and N(2)...XCu (X = F, Cl, Br) systems studied by means of the molecular orbital and Quantum Chemical Topology methods

Ab initio studies carried out at the MP2(full)/6-311+G(2df) and MP2(full)/aug-cc-pVTZ-PP computational levels reveals that dinitrogen (N(2)) and cuprous halides (CuX, X = F, Cl, Br) form three types of systems with the side-on and end-on coordination of N(2): N[triple bond]N-CuX (C(infinity v)), N(2)-CuX (C(2v)) stabilized by the donor-acceptor bonds and weak van der Waals complexes N(2)...XCu (C(2v)) with dominant dispersive forces. An electron density transfer between the N(2) and CuX depends on type of the N(2) coordination and a comparison of the NPA charges yields the [N[triple bond]N](delta+)-[CuX](delta-) and [N(2)](delta-)-[CuX](delta+) formula. According to the NBO analysis, the Cu-N coordinate bonds are governed by predominant LP(N2)-->sigma*(Cu-X) "2e-delocalization" in the most stable N[triple bond]N-CuX systems, meanwhile back donation LP(Cu)-->pi*(N-N) prevails in less stable N(2)-CuX molecules. A topological analysis of the electron density (AIM) presents single BCP between the Cu and N nuclei in the N[triple bond]N-CuX, two BCPs corresponding to two donor-acceptor Cu-N bonds in the N(2)-CuX and single BCP between electron density maximum of the N[triple bond]N bond and halogen nucleus in the van der Waals complexes N(2)...XCu. In all systems values of the Laplacian nabla(2)rho(r)(r(BCP)) are positive and they decrease following a trend of the complex stability i.e. N[triple bond]N-CuX (C(infinity v)) > N(2)-CuX (C(2v)) > N(2)...XCu (C(2v)). A topological analysis of the electron localization function (ELF) reveals strongly ionic bond in isolated CuF and a contribution of covalent character in the Cu-Cl and Cu-Br bonds. The donor-acceptor bonds Cu-N are characterized by bonding disynaptic basins V(Cu,N) with attractors localized at positions corresponding to slightly distorted lone pairs V(N) in isolated N(2). In the N[triple bond]N-CuX systems, there were no creation of any new bonding attractors in regions where classically the donor-acceptor bonds are expected and there is no sign of typical covalent bond Cu-N with the bonding pair. Calculations carried out for the N[triple bond]N-CuX reveal small polarization of the electron density in the N[triple bond]N bond, which is reflected by the bond polarity index being in range of 0.14 (F) to 0.11 (Cl).     

Effect of o-methyl group on rate, mechanism, and resonance contribution: aminolysis of Y-substituted phenyl X-substituted 2-methylbenzoates

Second-order rate constants have been determined spectrophotometrically for the reactions of 4-nitrophenyl X-substituted 2-methylbenzoates (2a-e) and Y-substituted phenyl 2-methylbenzoates (3a-e) with alicyclic secondary amines in 80 mol % H(2)O/20 mol % DMSO at 25.0 +/- 0.1 degrees C. The o-methyl group in the benzoyl moiety of 2a-e retards the reaction rate but does not influence the reaction mechanism. The Hammett plots for the reactions of 2a-e are nonlinear, while the corresponding Yukawa-Tsuno plots are linear with large r values (1.06-1.70). The linear Yukawa-Tsuno plots suggest that stabilization of the ground-state through resonance interaction between the electron donating substituent X and the carbonyl group is responsible for the nonlinear Hammett plots, while the large r values imply that the ground-state resonance interaction is significant. The reactions of 2a-e resulted in smaller rho(X) values but larger r values than the corresponding reactions of 4-nitrophenyl X-substituted benzoates (1a-e). The small rho(X) value for the reactions of 2a-e (e.g., rho(X) = 0.22) is suggested to be responsible for the large r value (e.g., r = 1.70). The reactions of 3a-e with piperidine are proposed to proceed in a stepwise manner with a change in the rate-determining step on the basis of the curved Br?nsted-type plot obtained. Microscopic rate constants associated with the reactions of 3a-e are also consistent with the proposed mechanism.     

Molecular macrocluster formation on silica surfaces in phenol-cyclohexane mixtures

The adsorption of phenol, an aromatic compound with a hydrogen-bonding group, onto a silica surface in cyclohexane was investigated by colloidal probe atomic force microscopy (AFM), attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR), and adsorption isotherm measurements. ATR-FTIR measurements on the silica surface indicated the formation of surface macroclusters of phenol through hydrogen bonding. The ATR-FTIR spectra were also measured on the H-terminated silicon surface to observe the effect of the silanol groups on the phenol adsorption. The comparison of the ATR-FTIR spectra for both the silicon oxide and H-terminated silicon surfaces proved that the silanol groups are necessary for the formation of phenol clusters on the surface. The surface force measurement using colloidal probe AFM showed a long-range attraction between the two silica surfaces in phenol-cyclohexane mixtures. This long-range attraction resulted from the contact of the adsorbed phenol layers for the phenol concentrations below 0.6 mol %, at which no significant phenol clusters formed in the bulk solution. The attraction started to decrease at 0.6 mol % phenol due to the exchange of the phenol molecules between the clusters in the bulk phase and on the surface. The surface density of phenol in the adsorbed layer was calculated on the basis of the long-range attraction and found to be much smaller than the liquid phenol density. The plausible structure of the adsorbed phenol layer was drawn by referring to the crystal structure of the bulk phenol and orientation of the phenol molecules on the surface, estimated by the dichroic analysis of ATR-FTIR spectroscopy. The investigation of the phenol adsorption on the silica surface in a nonpolar solvent using this novel approach demonstrated the effect of the aromatic ring on the surface packing density.     

DFT study on hydrogen-bonding adsorption mechanism of rutin onto macroporous adsorption resins functionalized with amino, hydroxyl, and carboxyl groups

The adsorption mechanism of a series of macroporous adsorption resins (p-(CH₃NH)PhL (L = NH₂, OH, COOH)) with rutin have been investigated using density functional theory calculations at B3LYP/6-31G(d,p) level of theory. Solvent effects on these species were explored using calculations that included a polarizable continuum model for the aqueous solvent. In this article, the geometry structure, interaction energies and the infrared spectra for the stable reactants and the adsorption complexes were obtained and analyzed. The results show that the hydrogen-bonding have been formed in the adsorption complexes. The higher interaction energy is calculated for the carboxyl group, while the resin with amino group has the highest adsorption capacity for rutin. The adsorption complexes become more and more stable as increasing the number of adsorbents. Our theoretical study is in good explanation for the experimental results.     

Water effect on the o-h dissociation enthalpy of para-substituted phenols: a DFT study

The effect of water on the O-H bond dissociation enthalpy (BDE) of para-substituted phenols has been investigated by means of DFT calculations. It is shown that the experimental BDE values are fairly well-reproduced by simple B3LYP/6-31G* calculations carried out on the phenol/phenoxyl-water complexes taking into account only hydrogen-bonding (HB) interactions of water molecules with molecular sites (HB model). On the contrary, the BDE values computed with the polarizable continuum model (PCM/B3LYP/6-31G*)8 are overestimated by about 3-4 kcal/mol. Discrepancy between theory and experiment increases using the PCM method in addition to the HB model. Calculations show that, in general, the HB interaction with water molecules decreases the BDE of phenols bearing electron-releasing groups while increasing the BDE of phenols bearing electron-withdrawing substituents. This opposite effect is explained by considering the resonance structures with charge separation both in phenols and in phenoxyl radicals. With electron donors, the phenoxyl radical is preferentially stabilized by the HB acceptor interaction with two water molecules, while with electron acceptors the phenol is preferentially stabilized by the HB donor interaction with one water molecule.     

Core-to-Rydberg band shift and broadening of hydrogen bonded ammonia clusters studied with nitrogen K-edge excitation spectroscopy

Nitrogen 1s (N 1s) core-to-Rydberg excitation spectra of hydrogen-bonded clusters of ammonia (AM) have been studied in the small cluster regime of beam conditions with time-of-flight (TOF) fragment-mass spectroscopy. By monitoring partial-ion-yield spectra of cluster-origin products, "cluster" specific excitation spectra could be recorded. Comparison of the "cluster" band with "monomer" band revealed that the first resonance bands of clusters corresponding to N 1s → 3sa(1)/3pe of AM monomer are considerably broadened. The changes of the experimental core-to-Rydberg transitions ΔFWHM (N 1s → 3sa(1)/3pe) = ~0.20/~0.50 eV compare well with the x ray absorption spectra of the clusters generated by using density functional theory (DFT) calculation. The broadening of the core-to-Rydberg bands in small clusters is interpreted as being primarily due to the splitting of non-equivalent core-hole N 1s states caused by both electrostatic core-hole and hydrogen-bonding (H(3)N···H-NH(2)) interactions upon dimerization. Under Cs dimer configuration, core-electron binding energy of H-N (H-donor) is significantly decreased by the intermolecular core-hole interaction and causes notable redshifts of core-excitation energies, whereas that of lone-pair nitrogen (H-acceptor) is slightly increased and results in appreciable blueshifts in the core-excitation bands. The result of the hydrogen-bonding interaction strongly appears in the n-σ* orbital correlation, destabilizing H-N donor Rydberg states in the direction opposite to the core-hole interaction, when excited N atom with H-N donor configuration strongly possesses the Rydberg component of anti-bonding σ* (N-H) character. Contributions of other cyclic H-bonded clusters (AM)(n) with n ≥ 3 to the spectral changes of the N 1s → 3sa(1)/3pe bands are also examined.     

Studies on the structure, stability, and spectral signatures of hydride ion-water clusters

The gas-phase structure, stability, spectra, and electron density topography of H(-)W(n) clusters (where n = 1-8) have been calculated using coupled-cluster CCSD(T) and M?ller-Plesset second-order perturbation (MP2) theory combined with complete basis set (CBS) approaches. The performance of various density functional theory (DFT) based methods such as B3LYP, M05-2X, M06, M06-L, and M06-2X using 6-311++G(d,p), and aug-cc-pVXZ (aVXZ, where X = D, T, and Q) basis sets has also been assessed by considering values calculated using CCSD(T)/CBS limit as reference. The performance of the functionals has been ranked based on the mean signed/unsigned error. The comparison of geometrical parameters elicits that the geometrical parameters predicted by B3LYP/aVTZ method are in good agreement with those values obtained at MP2/aVTZ level of theory. Results show that M05-2X functional outperform other functionals in predicting the energetics when compared to CCSD(T)/CBS value. On the other hand, values predicted by M06-2X, and M06 methods, are closer to those values obtained from MP2/CBS approach. It is evident from the calculations that H(-)W(n) (where n = 5-8) clusters adopt several interesting structural motifs such as pyramidal, prism, book, Clessidra, cubic, cage, and bag. The important role played by ion-water (O-H···H(-)) and water-water (O-H···O) interactions in determining the stability of the clusters has also been observed. Analysis of the results indicates that the most stable cluster is made up of minimum number of O-H···H(-) interaction in conjugation with the maximum number of O-H···O interactions. The Bader theory of atoms in molecules (AIM) and natural bond orbital (NBO) analyses has also been carried out to characterize the nature of interactions between hydride ion and water molecules. It can be observed from the vibrational spectra of H(-)W(n) clusters, the stretching frequencies involving ion-water interaction always exhibit larger redshift and intensities than that of water-water (inter solvent) interactions.