The capacity of SX2 (X = F, Cl, and Br) to engage in different kinds of noncovalent bonds was investigated by ab initio calculations. SCl2 (SBr2) has two σ-holes upon extension of Cl (Br)?S bonds, and two σ-holes upon extension of S?Cl (Br) bonds. SF2 contains only two σ-holes upon extension of the F?S bond. Consequently, SCl2 and SBr2 form chalcogen and halogen bonds with the electron donor H2CO while SF2 forms only a chalcogen bond, i.e., no F···O halogen bond was found in the SF2:H2CO complex. The S···O chalcogen bond between SF2 and H2CO is the strongest, while the strongest halogen bond is Br···O between SBr2 and H2CO. The nature of these two types of noncovalent interaction was probed by a variety of methods, including molecular electrostatic potentials, QTAIM, energy decomposition, and electron density shift maps. Termolecular complexes X2S···H2CO···SX′2 (X = F, Cl, Br, and X′ = Cl, Br) were constructed to study the interplay between chalcogen bonds and halogen bonds. All these complexes contained S···O and Cl (Br)···O bonds, with longer intermolecular distances, smaller values of electron density, and more positive three-body interaction energies, indicating negative cooperativity between the chalcogen bond and the halogen bond. In addition, for all complexes studied, interactions involving chalcogen bonds were more favorable than those involving halogen bonds.
Based on the structure of MOF-808, different substituents were introduced to replace hydrogen atom on the phenyl ring of MOF-808. The GCMC method was used to study the effect of functional groups on the hydrogen storage properties of MOF-808-X (X?=??OH, ?NO2, ?CH3, ?CN, ?I). The H2 uptakes and isosteric heat of adsorption were simulated at 77 K. The results indicate that all these substituents have favorable impact on the hydrogen storage capacity, and –CN is found to be the most promising substituent to improve H2 uptake. These results may be helpful for the design of MOFs with higher hydrogen storage capacity.
Graphical abstract Atomistic structures of MOFs. (a) The structures of MOF-808-X. (b) Model of organic linker. Atom color scheme: C, gray; H, white; O, red; X, palegreen (X?=??OH, ?NO2, ?CH3, ?CN, ?I)
Detailed electrostatic potential (ESP) analyses were performed to compare the directionality of halogen bonds with those of hydrogen bonds and lithium bonds. To do this, the interactions of HOOOH with the molecules XF (X?=?Cl, Br, H, Li) were investigated. For each molecule, the percentage of the van der Waals (vdW) molecular surface that intersected with the ESP surface was used to roughly quantify the directionality of the halogen/hydrogen/lithium bond associated with the molecule. The size of the region of intersection was found to increase in the following order: ClF?<?BrF?<?HF?<?LiF. The maximum ESP in the region of intersection, VS, max, was observed to become more positive according to the sequence ClF?<?BrF?<?HF?<?LiF. For ClF and BrF, the positive electrostatic potential was concentrated in a very small region of the vdW molecular surface. On the other hand, for HF and LiF, the positive electrostatic potential was more diffusely scattered across the vdW surface than for ClF and BrF. Also, the optimized geometries of the dipolymers HOOOH···?XF (X?=?Cl, Br, H, Li) indicated that halogen bonds are more directional than hydrogen bonds and lithium bonds, consistent with the results of ESP analyses.
The σ-hole and π-hole of the protonated 2-halogenated imidazolium cation (XC3H4N2+; X = F, Cl, Br, I) were investigated and analyzed. The monomers of (CH3)3SiY(Y=F, Cl, Br, I), considered as the Lewis base, were combined with the σ-hole and π-hole of XC3H4N2+ to form the σ-hole and π-hole interactions in the bimolecular complexes (CH3)3SiY?·?·?·?XC3H4N2+ and (CH3)3SiY?·?·?·?C3(X)H4N2+(X/Y=F, Cl, Br, I), respectively. For both the σ-hole and π-hole interactions, the equilibrium geometries of complexes show regular changes according to the sequence of heavy sequence of the noncovalent interaction acceptors and donors. The electrostatic energy is the main contribution in the formation of both kinds of interactions, it has linear relations with the VS,max values of σ-hole and the V′S,max values of π-hole. Both the σ-hole and π-hole interactions belong to the closed-shell and noncovalent interactions. The π-hole interactions are stronger than the σ-hole interactions. For the π-hole interactions, the contribution percents of the dispersion energies are somewhat greater than those of the σ-hole interactions, while it is contrary for the polarization energy.
Quantum chemical calculations were performed to investigate the stability of the ternary complexes BeH2···XMH3···NH3 (X?=?F, Cl, and Br; M?=?C, Si, and Ge) and the corresponding binary complexes at the atomic level. Our results reveal that the stability of the XMH3···BeH2 complexes is mainly due to both a strong beryllium bond and a weak tetrel–hydride interaction, while the XMH3···NH3 complexes are stabilized by a tetrel bond. The beryllium bond with a halogen atom as the electron donor has many features in common with a beryllium bond with an O or N atom as the electron donor, although they do exhibit some different characteristics. The stability of the XMH3···NH3 complex is dominated by the electrostatic interaction, while the orbital interaction also makes an important contribution. Interestingly, as the identities of the X and M atoms are varied, the strength of the tetrel bond fluctuates in an irregular manner, which can explained by changes in electrostatic potentials and orbital interactions. In the ternary systems, both the beryllium bond and the tetrel bond are enhanced, which is mainly ascribed to increased electrostatic potentials on the corresponding atoms and charge transfer. In particular, when compared to the strengths of the tetrel and beryllium bonds in the binary systems, in the ternary systems the tetrel bond is enhanced to a greater degree than the beryllium bond.
Intramolecular and intermolecular hydrogen bonding in electronic excited states of calixarene building blocks bis(2-hydroxyphenyl)methane (2HDPM) monomer and hydrogen-bonded 2HDPM-H2O complex were studied theoretically using the time-dependent density functional theory (TDDFT). Twenty-four stable conformations (12 pairs of enantiomers) of 2HDPM monomer have been found in the ground state. From the calculation results, the conformations 1a and 1b which both have an intramolecular hydrogen bond are the most stable ones. The infrared spectra of 2HDPM monomer and 2HDPM-H2O complex in ground state and S1 state were calculated. The stretching vibrational absorption band of O2???H3 group in the monomer and complex disappeared in the S1 state. At the same time, a new strong absorption band appeared at the C=O stretching region. From the calculation of bond lengths, it indicates that the O2???H3 bond is significantly lengthened in the S1 state. However, the C1???O2 bond is drastically shortened upon electronic excitation to the S1 state and has the characteristics of C=O band. Furthermore, the intramolecular hydrogen bond O2???H3?·?·?·?O4 of the 2HDPM monomer and the intermolecular hydrogen bonds O2???H3?·?·?·?O7 and O7???H9?·?·?·?O4 of 2HDPM-H2O complex are all shortened and strengthened in the S1 state.
Figure
Intramolecular and intermolecular hydrogen bonding in electronic excited states of calixarene building blocks bis(2-hydroxyphenyl)methane (2HDPM) monomer and hydrogen-bonded 2HDPM-H2O complex were studied by TDDFT method 相似文献
Second group metal dimers can replace the carbon atom in benzene to form metallabenzene (C5H6M2) compounds. These complexes possess some aromatic character and promising hydrogen adsorption properties. In this study, we investigated the aromatic character of these compounds using aromaticity indices and molecular orbital analysis. To determine the nature of interactions between hydrogen and the metallic center, variation-perturbational decomposition of interaction energy was applied together with ETS-NOCV analysis. The results obtained suggest that the aromatic character comes from three π orbitals located mainly on the C5H5? fragment. The high hydrogen adsorption energy (up to 6.5 kcal mol?1) results from two types of interaction. In C5H6Be2, adsorption is controlled by interactions between the empty metal orbital and the σ orbital of the hydrogen molecule (Kubas interaction) together with corresponding back-donation interactions. Other C5H6M2 compounds adsorb H2 due to Kubas interactions enhanced by H2–π interactions.
Detailed ab initio molecular orbital calculations on the interactions of molecular hydrogen, H2, with various poly-aromatic hydrocarbons (PAHs) as a model system for graphene were carried out to accurately describe the physisorption phenomenon. The binding energies corrected for the basis set superposition error, ΔEbind(BSSE), were obtained using the optimized geometries at the MP2 level with a large basis set and were compared with the single point binding energies, denoted as ΔEbind(BSSE-s), using large basis sets on the geometries optimized at the small basis sets, such as SVP and TZVP. The calculations showed that the ΔEbind(BSSE-s) values were similar to those at the MP2 level with the large basis sets. The binding strength increased gradually with increasing size of the PAHs. The ΔEbind(BSSE-s) for an infinite graphene sheet was estimated to be ?1.70 kcal mol?1 using the non-linear curve fitting method. The present work could be expected to provide more useful and reliable information on H2 physisorption.
Graphical abstract Detailed ab initio molecular orbital calculations on the interactions of molecular hydrogen with various poly-aromatic hydrocarbons as a model system for graphene indicate that the perpendicular type A is the most favorable and the binding energy on an infinite graphene sheet is estimated to be ?1.70 kcal mol?1.
In the present work, in order to investigate the electronic excited-state intermolecular hydrogen bonding between the chromophore coumarin 153 (C153) and the room-temperature ionic liquid N,N-dimethylethanolammonium formate (DAF), both the geometric structures and the infrared spectra of the hydrogen-bonded complex C153–DAF+ in the excited state were studied by a time-dependent density functional theory (TDDFT) method. We theoretically demonstrated that the intermolecular hydrogen bond C1?=?O1···H1–O3 in the hydrogen-bonded C153–DAF+ complex is significantly strengthened in the S1 state by monitoring the spectral shifts of the C=O group and O–H group involved in the hydrogen bond C1?=?O1···H1–O3. Moreover, the length of the hydrogen bond C1?=?O1···H1–O3 between the oxygen atom and hydrogen atom decreased from 1.693 Å to 1.633 Å upon photoexcitation. This was also confirmed by the increase in the hydrogen-bond binding energy from 69.92 kJ mol?1 in the ground state to 90.17 kJ mol?1 in the excited state. Thus, the excited-state hydrogen-bond strengthening of the coumarin chromophore in an ionic liquid has been demonstrated theoretically for the first time. 相似文献
MP2(full)/aug-cc-pVDZ(-PP) computations predict that new triangular bonding complexes (where X? is a halide and H–C refers to a protic solvent molecule) consist of one halogen bond and two hydrogen bonds in the gas phase. Carbon tetrabromide acts as the donor in the halogen bond, while it acts as an acceptor in the hydrogen bond. The halide (which commonly acts as an acceptor) can interact with both carbon tetrabromide and solvent molecule (CH3CN, CH2Cl2, CHCl3) to form a halogen bond and a hydrogen bond, respectively. The strength of the halogen bond obeys the order CBr4???Cl? > CBr4???Br? > CBr4???I?. For the hydrogen bonds formed between various halides and the same solvent molecule, the strength of the hydrogen bond obeys the order C-H???Cl? > C-H???Br? > C-H???I?. For the hydrogen bonds formed between the same halide and various solvent molecules, the interaction strength is proportional to the acidity of the hydrogen in the solvent molecule. The diminutive effect is present between the hydrogen bonds and the halogen bond in chlorine and bromine triangular bonding complexes. Complexes containing iodide ion show weak cooperative effects.
Figure
The triangular bonding complexes consisting of halogen bond and hydrogen bonds were predict in the gas phase by computational quantum chemistry. 相似文献
To understand the chemical behavior of uranyl complexes in water, a bis-uranyl [(phen)(UO2)(μ2–F)(F)]2 (A; phen?=?phenanthroline, μ2?=?doubly bridged) and its hydrated form A?·?(H2O)n (n?=?2, 4 and 6) were examined using scalar relativistic density functional theory. The addition of water caused the phen ligands to deviate slightly from the U2(μ2–F)2 plane, and red-shifts the U–F-terminal and U?=?O stretching vibrations. Four types of hydrogen bonds are present in the optimized hydrated A?·?(H2O)n complexes; their energies were calculated to fall within the range 4.37–6.77 kcal mol-1, comparable to the typical values of 5.0 kcal mol-1 reported for hydrogen bonds. An aqueous environment simulated by explicit and/or implicit models lowers and re-arranges the orbitals of the bis-uranyl complex.
Figure
A bis(uranyl) complex [(phen)(UO2)(μ2–F)(F)]2 (A) and its solvated form A?·?(H2O)n were examined using scalar relativistic density functional theory. Hydrogen bonds cause the phen ligand to slightly deviate from the equatorial plane of the uranyl ion, resulting in a pronounced red-shift of the U–F-terminal and U?=?O asymmetric stretching vibrations. The calculated energies fall within 4.4?–6.8 kcal/mol. Explicit and/or implicit aqueous solvation re-arranges the molecular orbitals of the complex 相似文献
In this work, computations of density functional theory (DFT) were carried out to investigate the nature of interactions in solid 2,6-dibromo-4-nitroaniline (DBNA). This system was selected to mimic the hydrogen/halogen bonding found within crystal structures as well as within biological molecules. DFT (M06-2X/6-311++G**) calculations indicated that the binding energies for different of interactions lie in the range between ?1.66 and ?9.77 kcal mol?1. The quantum theory of atoms in molecules (QTAIM) was applied to provide more insight into the nature of these interactions. Symmetry-adapted perturbation theory (SAPT) analysis indicated that stability of the Br···Br halogen bonds is predicted to be attributable mainly to dispersion, while electrostatic forces, which have been widely believed to be responsible for these types of interactions, play a smaller role. Our results indicate that, for those nuclei participating in hydrogen/halogen bonding interactions, nuclear quadrupole resonance parameters exhibit considerable changes on going from the isolated molecule model to crystalline DBNA.
Figure
Electrostatic potential mapped on the surface of 2,6-dibromo-4-nitroaniline (DBNA) molecular electron density (0.001 e au?3). Color ranges for VS(r), in kcal?mol?1: red > 26.5, yellow 26.5–5.7, green 5.7– ?15.1, blue < ?15.1. Black circles Surface maxima, blue surface minima 相似文献
The character of the cooperativity between the HOX···OH/SH halogen bond (XB) and the Y―H···(H)OX hydrogen bond (HB) in OH/SH···HOX···HY (X = Cl, Br; Y = F, Cl, Br) complexes has been investigated by means of second-order Møller?Plesset perturbation theory (MP2) calculations and “quantum theory of atoms in molecules” (QTAIM) studies. The geometries of the complexes have been determined from the most negative electrostatic potentials (VS,min) and the most positive electrostatic potentials (VS,max) on the electron density contours of the individual species. The greater the VS,max values of HY, the larger the interaction energies of halogen-bonded HOX···OH/SH in the termolecular complexes, indicating that the ability of cooperative effect of hydrogen bond on halogen bond are determined by VS,max of HY. The interaction energies, binding distances, infrared vibrational frequencies, and electron densities ρ at the BCPs of the hydrogen bonds and halogen bonds prove that there is positive cooperativity between these bonds. The potentiation of hydrogen bonds on halogen bonds is greater than that of halogen bonds on hydrogen bonds. QTAIM studies have shown that the halogen bonds and hydrogen bonds are closed-shell noncovalent interactions, and both have greater electrostatic character in the termolecular species compared with the bimolecular species.
Figure
The character of the cooperativity between the X···O/S halogen bond (XB) and the Y―H···O hydrogen bond (HB) in OH/SH···HOX···HY (X=Cl, Br; Y=F, Cl, Br) complexes has been investigated by means of second-order Møller—Plesset perturbation theory (MP2) calculations and “quantum theory of atoms in molecules” (QTAIM) studies. 相似文献
This article analyzes the substitution effects on cooperativity between fluorin-centered halogen bonds in NCF?·?·?·?NCF?·?·?·?NCX and CNF?·?·?·?CNF?·?·?·?CNX complexes, where X?=?H, F, Cl, CN, OH, and NH2. These effects are investigated theoretically in terms of geometric and energetic features of the complexes, which are computed by ab initio methods. The topological analysis, based on the quantum theory of atoms in molecules (QTAIM), is used to characterize the interactions and analyze their enhancement with varying electron density at bond critical points. It is found that the complexes with electron-donating groups exhibit a strong cooperativity, while a much weaker cooperativity occurs in the NCF?·?·?·?NCF?·?·?·?NCCN and CNF?·?·?·?CNF?·?·?·?CNCN trimers. An excellent correlation is found between the cooperative energy in the ternary complexes and the calculated three-body interaction energies. The energy decomposition analysis (EDA) indicates that the electrostatic and dispersion effects play a main role in the cooperativity of fluorine-centered halogen bonding.
Figure
Structure of NCF···NCF···NCX and CNF···CNF···CNX complexes 相似文献
The interaction of external water molecules with hydrated pyrrole-2-carboxaldehyde PCL/(H2O)n complexes was investigated. The work was supported by both theoretical [DFT/TD-DFT methods using 6-311G++(d,p) basis set in the ground (S0) and excited (S1, S2, S3)states] and experimental [UV-Vis, FTIR and Raman] verification. The focus of the present work was on the weak intermolecular O–H?O, N–H?O–H hydrogen bonded interaction (IerHB) between PCL and external water molecules, and the influence of increasing the number of water molecules to form hydrated PCL/(H2O)n complexes. Effects were observed on different vibrational normal modes and on electronic transition levels. A hydrogen-bonded network of water induces a shift to higher energy in certain normal modes of PCL to form stable PCL/(H2O)n complexes by lowering the barrier energy. Potential energy distribution (PED) analysis indicates a significant charge transfer from PCL to water by creating a water bridge. Hydrogen bonding effects account for the substantial red shift and broadness in νNH, νCO vibrational modes. Water rearrangement turns out to be the main driving force for hydrated complex formation.
The electric dipole transitions between pure spin and mixed spin electronic states are calculated at the XMC-QDPT2 and MCSCF levels of theory, respectively, for different intermolecular distances of the C6H6 and O2 collisional complex. The magnetic dipole transition moment between the mixed-spin ground (“triplet”) and the first excited (“singlet”) states is calculated by quadratic response at MCSCF level of theory. The obtained results confirm the theory of intensity borrowing and increasing the intensity of electronic transitions in the C6H6?+?O2 collision. The calculation of magnetically induced current density is performed for benzene molecule being in contact with O2 at the distances from 3.5 to 4.5 Å. The calculation shows that the aromaticity of benzene is rising due to the conjugation of π-MOs of both molecules. The C6H6?+?O2 complex becomes nonaromatic at the short distances (r?<?3.5 Å). The computation of static polarizability in the excited electronic states of the C6H6?+?O2 collisional complex at various distances supports the theory of red solvatochromic shift of the a?→?X band.
The structures of the N-(hydroxymethyl)acetamide (model molecule of ceramide) dimers have been fully optimized at B3LYP/6–311++G** level. The intermolecular hydrogen bonding interaction energies have been calculated using the B3LYP/6–311++G**, B3LYP/6–311++G(2df,2p), MP2(full)/6–311++G** and MP2(full)/6–311++G(2df,2p) methods, respectively. The results show that the O–H···O, N–H···O, O–H···N, and C–H···O hydrogen bonding interactions could exist in N-(hydroxymethyl)acetamide dimers, and the O–H···O, N–H···O, and O–H···N hydrogen bonding interactions could be stronger than C–H···O. The three-dimensional network structure formed by ceramide molecules through intermolecular hydrogen bonding interactions may be the main reason why the stratum corneum of skin could prevent foreign substances from entering our body, as is in accordance with the experimental results. The stability of hydrogen-bonding interactions follow the order of (a)?>?(b)?≈?(c)?>?(d)?>?(e)?≈?(f)?>?(g)?>?(h). The analyses of the energy decomposition, frequency, atoms in molecules (AIM), natural bond orbital (NBO), and electron density shift are used to further reveal the nature of the complex formation. In the range of 263.0–328.0 K, the complex is formed via an exothermic reaction, and the solvent with lower temperature and dielectric constant is favorable to this process.
Graphical abstract The structures and the O–H···O=C, N–H···O=C and C–H···O=C H-bonding interactions in the N-(hydroxymethyl)acetamide (model molecule of ceramide) dimers were investigated using the B3LYP and MP2(full) methods.
The insertion reactions of the silylene H2Si with H2BXHn-1 (X?=?F, Cl, Br, O, N; n?=?1, 1, 1, 2, 3) have been studied by DFT and MP2 methods. The calculations show that the insertions occur in a concerted manner, forming H2Si(BH2)(XHn-1). The essences of H2Si insertions with H2BXHn-1 are the transfers of the σ electrons on the Si atom to the positive BH2 group and the electrons of X into the empty p orbital on the Si atom in H2Si. The order of reactivity in vacuum shows the barrier heights increase for the same-family element X from up to down and the same-row element X from right to left in the periodic table. The energies relating to the B-X bond in H2BXHn-1, and the bond energies of Si-X and Si-B in H2Si(BH2)(XHn-1) may determine the preference of insertions of H2Si into B-X bonds for the same-column element X or for the same-row element X. The insertion reactions in vacuum are similar to those in solvents, acetone, ether, and THF. The barriers in vacuum are lower than those in solvents and the larger polarities of solvents make the insertions more difficult to take place. Both in vacuum and in solvents, the silylene insertions are thermodynamically exothermic.
