MP2/aug-cc-pVTZ calculations were performed for complexes linked by hydrogen bonds. Three types of proton donating species were taken into account: H2O, CCl3H, and H3O+. These calculations are supported by the natural bond orbital (NBO) method and the quantum theory of atoms in molecules (QTAIM) approach. Numerous correlations between parameters of H-bonded systems were found. The most important are those which show the response of the system on the H-bond formation; for example, the increase of polarization of the A-H bond correlates with the strength of the hydrogen bond. Similar relationships were found for the σ-hole bonds while the π-hole bonds do not follow the trends known for the hydrogen bonds.
Protonation in the two-electron/two-proton reduction processes of 2,6-dichlorophenolindophenolate (DCIP) is investigated combining density functional theory (DFT) and molecular dynamics (MD) methods. DCIP (anion), DCIP?– (radical anion), and DCIP2? (dianion) are considered, including the electronic structure analysis from the prospective of quantum theory of atoms and molecules (QTAIM). It is shown that oxygen on the indophenolate moiety and nitrogen are the first and/or the second proton acceptor sites and their energetic order depends on the total charge of the system. MD simulations of differently charged species interacting with the solvent molecules have been performed for methanol, water, and oxonium cation (H3O+). Methanol and water molecules are found to form only hydrogen bonds with the solute irrespective of its charge. The calculated pKa values show that the imino group of DCIPH? is a weaker acid than water. While in the case of DCIP (and DCIP?–) plus oxonium cation, proton transfer from the solvent to the solute was evidenced for both aforementioned acceptor sites. In addition, MD simulations of bulks containing 15 and 43 molecules of water around the DCIP molecule have been performed, revealing the formation of 2–4 hydrogen bonds.
Graphical Abstract 2,6-Dichlorophenolindophenolate interacts with solvent molecules (water, oxonium cation and methanol). Hydrogen transfer and electronic structure are studied by DFT and molecular dynamics methods
The evolution of structural properties, thermodynamics and averaged (dynamic) total hardness values as a function of the composition of binary water–organic solvents, was rationalized in view of the intermolecular interactions. The organic solvents considered were ethanol, acetonitrile, and isopropanol at 0.25, 0.5, 0.75, and 1 mass fractions, and the results were obtained using molecular dynamics simulations. The site-to-site radial distribution functions reveal a well-defined peak for the first coordination shell in all solvents. A characteristic peak of the second coordination shell exists in aqueous mixtures of acetonitrile, whereas in the water–alcohol solvents, a second peak develops with the increase in alcohol content. From the computed coordination numbers, averaged hydrogen bonds and their lifetimes, we found that water mixed with acetonitrile largely preserves its structural features and promotes the acetonitrile structuring. Both the water and alcohol structures in their mixtures are disturbed and form hydrogen bonds between molecules of different kinds. The dynamic hardness values are obtained as the average over the total hardness values of 1200 snapshots per solvent type, extracted from the equilibrium dynamics. The dynamic hardness profile has a non-linear evolution with the liquid compositions, similarly to the thermodynamic properties of these non-ideal solvents.
Graphical abstract Computed dynamic total hardness, as a function of the cosolvent mass fraction for water–ethanol (EtOH), water–isopropanol (2PrOH) and water–acetonitrile (AN)
Using armchair-type single-walled carbon nanotubes (SWCNTs) of different sizes as model compounds for lignite, the effect of water molecule confinement on the water-holding capacity of lignite pores was investigated. Results indicated that the water-holding capacity of pores with diameters of <10 nm was eight times larger than that of pores with diameters of 100 nm. The configuration of the cluster of water molecules in each SWCNT and the binding energy between each SWCNT and the water molecules within it were calculated by means of density functional theory using a hybrid functional: M06-2X/6-311+G**, 6-31G*. The results prove that the configurations of the water molecules in the SWCNTs are very different to their configuration in the unconfined state. In vacuum, the cluster of three water molecules adopted a trimer configuration, while they presented a linear configuration in the 6.78 Å SWCNT. Similarly, in vacuum, the cluster of five water molecules formed a five-membered ring, while they favored a linear configuration in the 6.78 Å SWCNT, a zigzag configuration in the 8.14 Å SWCNT, and a trimer?+?1?+?1 configuration (i.e., a trimer plus two isolated water molecules) in the 9.49 Å, 10.85 Å, and 13.75 Å SWCNTs. There was found to be a degree of competition between the coupling energy of the water molecules with the SWCNT and the hydrogen bonding among the water molecules. When the diameter of the SWCNT was >1 nm, the hydrogen bonding among the water molecules dominated, while the coupling energy of the water molecules with the SWCNT amounted to only 30–40% of the total interaction energy of the water molecules.
Graphical Abstract Computed equilibrium structures of five water molecules confined in SWCNTs with diameters of 6.78 Å, 8.14 Å, 9.49 Å, 10.85 Å, and 13.75 Å, and in vacuum
The anti-hypertensive drugs amlodipine, atenolol and lisinopril, in ordinary and PEGylated forms, with different combined-ratios, were studied by molecular dynamics simulations using GROMACS software. Twenty simulation systems were designed to evaluate the interactions of drug mixtures with a dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane, in the presence of water molecules. In the course of simulations, various properties of the systems were investigated, including drug location, diffusion and mass distribution in the membrane; drug orientation; the lipid chain disorder as a result of drug penetration into the DMPC membrane; the number of hydrogen bonds; and drug surface area. According to the results obtained, combined drugs penetrate deeper into the DMPC lipid bilayer membrane, and the lipid chains remain ordered. Also, the combined PEGylated drugs, at a combination ratio of 1:1:1, enhance drug penetration into the DMPC membrane, reduce drug agglomeration, orient the drug in a proper angle for easy penetration into the membrane, and decrease undesirable lipotoxicity due to distorted membrane self-assembly and thickness.
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.
Hydrogen cyanide (HCN) and its tautomer hydrogen isocyanide (HNC) are relevant for extraterrestrial chemistry and possible relation to the origin of biomolecules. Several processes and reactions involving these molecules depend on their intermolecular interactions that can lead to aggregates and liquids especially due to the hydrogen bonds. It is thus important to comprehend, to describe, and to quantify their hydrogen bonds, mainly their nature and the cooperativity effects. A systematic study of all linear complexes up to pentamers of HCN and HNC is presented. CCSD(T)/CBS energy calculations, with and without basis set superposition error (BSSE) corrections for energies and geometries, provided a suitable set of benchmarks. Approximated methods based on the density functional theory (DFT) such as BP86, PBE, TPSS, B3LYP, CAM-B3LYP with and without dispersion corrections and long-range corrections, were assessed to describe the interaction energies and cooperativity effects. These assessments are relevant to select DFT functionals for liquid simulations. Energy decomposition analysis was performed at the PBE/STO-TZ2P level and provided insights into the nature of the hydrogen bonds, which are dominated by the electrostatic component. In addition, several linear relationships between the various energy components and the interaction energy were obtained. The cooperativity energy was also found to be practically linear with respect to the interaction energy, which may be relevant for designing and/or correcting empirical force fields.
Mechanisms for the activation of water, ammonia, and other small molecules by the PCcarbeneP nickel pincer complex were studied computationally with the aid of density functional theory. The calculation results indicate that the strongly donating, nucleophilic carbene center can engage in a variety of heterolytic splitting of E?H (E=H, C, N, O) bonds, some of which are reversible. The cleavage of E?H bonds across the Ni=C bond represents a new mode of bond activation by ligand cooperativity in nickel pincer complex. On the basis of the calculations, we also demonstrate that reversible H2 activation across the Ir=C bond via the PCcarbeneP iridium pincer complex was observed in the experiments, while other E?H (E=C, N, O) bonds were not activated. Our calculations are in good agreement with experimental observations and could provide new insights into ligand cooperativity in nickel pincer complexes.
Bond critical points (BCPs) in the quantum theory of atoms in molecules (QTAIM) are shown to be a consequence of the molecular topology, symmetry, and the Poincaré-Hopf relationship, which defines the numbers of critical points of different types in a scalar field. BCPs can be induced by a polarizing field or by addition of a single non-bonded atom to a molecule. BCPs and their associated bond paths are therefore suggested not to be a suitable means of identifying chemical bonds, or even attractive intermolecular interactions.
The correlation between the kinetic stability of molecules against temperature and variations in their geometric structure under optical excitation is investigated by the example of different organic pheromone molecules sensitive to temperature or ultraviolet radiation using the density functional theory. The kinetic stability is determined by the previously developed method based on the calculation of the probability of extension of any structural bond by a value exceeding the limit value Lмах corresponding to the breaking of the bond under temperature excitation. The kinetic stability calculation only requires the eigenfrequencies and vibrational mode vectors in the molecule ground state to be calculated, without determining the transition states. The weakest bonds in molecules determined by the kinetic stability method are compared with the bond length variations in molecules in the excited state upon absorption of light by a molecule. Good agreement between the results obtained is demonstrated and the difference between them is discussed. The universality of formulations within both approaches used to estimate the stability of different pheromone molecules containing strained cycles and conjugated, double, and single bonds allows these approaches to be applied for studying other molecules.
Molecular simulations were carried out to study the sodium dodecyl sulfate (SDS) surfactant with the interleukin 8 (IL8) protein as a model to investigate the influence of amphiphilic molecules on proteins. Simulations for an SDS micelle with an IL8 protein show that both aggregates, which were initially separated, eventually approach each other to form a single complex. The results showed that the protein was attached to the SDS micelle by the charged positive amino acids whereas less contacts were observed for the negatively charged amino acids. Structural protein properties, such as amino acid contacts and pair correlation functions were conducted between the micelle and the protein groups and they showed greater interactions between the surfactant headgroups and the positively charged residues in the protein. Moreover, hydrogen bonds were also calculated between both structures and a greater number of bonds among the SDS headgroups and the charged positive amino acids in the protein was found.
Graphical Abstract Attachment of the interleukin 8 protein with a sodium dodecyl sulfate (SDS) micelle is given by the charged positive amino acids as indicated by the interaction of those amino acids with the SDS headgroups.
Molecular dynamics simulations were performed to investigate the separation of trihalomethanes (THMs) from water using boron nitride nanosheets (BNNSs). The studied systems included THM molecules and a functionalized BNNS membrane immersed in an aqueous solution. An external pressure was applied to the z axis of the systems. Two functionalized BNNSs with large fluorinated-hydrogenated pore (F-H-pores) and small hydrogen-hydroxyl pore (H-OH-pores) were used. The pores of the BNNS membrane were obtained by passivating each nitrogen and boron atoms at the pore edges with fluorine and hydrogen atoms in the large pore or with hydroxyl and hydrogen atoms in the small pore. The results show that the BNNS with a small functionalized pore was impermeable to THM molecules, in contrast to the BNNS with a large functionalized pore. Using these membranes, water contaminants can be removed at lower cost.
Graphical Abstract A snapshot of the simulation system. The BNNS membrane with the large functionalized pore is located in the middle of the box. The size of the box is 3 × 3 × 5 nm3. Green chlorine, cyan carbon, red oxygen, white hydrogen
Density functional theory (DFT) calculations are performed to study the hydrogen-bonding in the DMSO-water and DMF-water complexes. Quantitative molecular electrostatic potential (MESP) and atoms-in-molecules (AIM) analysis are applied to quantify the relative complexation of DMSO and DMF with water molecules. The interaction energy of DMSO with water molecules was higher than in DMF-water complexes. The existence of cooperativity effect helps in the strong complex formation. A linear dependence was observed between the hydrogen bond energies EHB, and the total electron densities in the BCP’s of microsolvated complexes which supports the existence of cooperativity effect for the complexation process. Due to the stronger DMSO/DMF and water interaction, the water molecules in the formed complexes have a different structure than the isolated water clusters. NCI analysis shows that the steric area is more pronounced in DMF-water complex than the DMSO-water complex which accounts for the low stability of DMF-water complexes compared to the DMSO-water complex.
Graphical abstract NCI analysis shows that the steric area is more pronounced in DMF-water complex than the DMSO-water complex which accounts for the low stability of DMF-water complexes compared to the DMSO-water complex.
Thirteen X-ray crystal structures containing various non-covalent interactions such as halogen bonds, halogen–halogen contacts and hydrogen bonds (I?N, I?F, I?I, F?F, I?H and F?H) were considered and investigated using the DFT-D3 method (B97D/def2-QZVP). The interaction energies were calculated at MO62X/def2-QZVP and MP2/aug-cc-pvDZ level of theories. The higher interaction and dispersion energies (2nd crystal) of ?9.58 kcal mol?1 and ?7.10 kcal mol?1 observed for 1,4-di-iodotetrafluorobenzene bis [bis (2-phenylethyl) sulfoxide] structure indicates the most stable geometrical arrangement in the crystal packing. The electrostatic potential values calculated for all crystal structures have a positive σ-hole, which aids understanding of the nature of σ-hole bonds. The significance of the existence of halogen bonds in crystal packing environments was authenticated by replacing iodine atoms by bromine and chlorine atoms. Nucleus independent chemical shift analysis reported on the resonance contribution to the interaction energies of halogen bonds and halogen–halogen contacts. Hirshfeld surface analysis and topological analysis (atoms in molecules) were carried out to analyze the occurrence and strength of all non-covalent interactions. These analyses revealed that halogen bond interactions were more dominant than hydrogen bonding interactions in these crystal structures.
Graphical Abstract Molecluar structure of 1,4-Di-iodotetrafluorobenzene bis(thianthrene 5-oxide) moelcule and its corresponding molecular electrostatic potential map for the view of σ-hole.
Present molecular dynamics simulations indicate that the methanol component in a methanol/water mixture is more likely to be trapped in a cyclic peptide nanotube (CPNT), while water molecules tend to be present at the channel mouths as transient guests. Channel water resides mainly between methanol and the CPNT wall, resulting in a distinct decrease in the H-bond number per channel methanol. Six designed CPNTs with different channel diameters and outer surface characteristics all possess distinct selectivity to methanol over water. Of these, the amphipathic 8?×?(AQ)4-CPNT exhibits the best performance. Results in this study provide basic information for the application of a CPNT to enrich methanol from a methanol/water mixture.
The diffusion coefficients of 14 n-alkanes (ranging from methane to n-tetradecane) in liquid and supercritical methanol at infinite dilution (at a pressure of 10.5 MPa and at temperatures of 299 K and 515 K) were deduced via molecular dynamics simulations. Values for the radial distribution function, coordination number, and number of hydrogen bonds were then calculated to explore the local structure of each fluid. The flexibility of the n-alkane (as characterized by the computed dihedral distribution, end-to-end distance, and radius of gyration) was found to be a major influence and hydrogen bonding to be a minor influence on the local structure. Hydrogen bonding reduces the flexibility of the n-alkane, whereas increasing the temperature enhances its flexibility, with temperature having a greater effect than hydrogen bonding on flexibility.
Graphical abstract The flexibility of the alkane is a major influence and the hydrogen bonding is a minor influence on the first solvation shell; the coordination numbers of long-chain n-alkanes in the first solvation shell are rather low
A dispersion correction is introduced and tested for MNDO. The shift in electron density caused by the interaction between oscillating dipoles in the London picture of dispersion is mimicked by adding a small r?7-dependent attractive nucleus–electron potential to the core Hamiltonian. This potential results in a shift in electron density similar to that used by Feynman to explain dispersion. The resulting parameterized self-consistent and inherently multicenter treatment (MNDO-F) gives good results for CHNO compounds that do not exhibit hydrogen bonds, which MNDO cannot reproduce. This “Feynman” dispersion correction is also applicable to Hartree–Fock and density functional theory.
In this article, we explore the capacity of formed Schiff base complexes to trap metal atoms or ions, using their aromatic ends. The intrinsic geometry of each complex defines the process of substitution. Two cases were studied; one involving a trans Schiff base complex and the other considering how a salen ligand, with nickel systems traps chromium. We also assessed the nature of the new bonds and the frontier molecular orbitals.
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.
