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

Microhydration of small protonated polyaromatic hydrocarbons: a first principles study

Using first principles methodology, we investigate the microsolvation of protonated benzene (BzH+), protonated coronene (CorH+) and protonated dodecabenzocoronene (DbcH+). Gas phase complexes of these small protonated polyaromatic hydrocarbons (H+PAHs) with mono-, di-, and tri-hydrated water molecules are considered. Their most stable forms are presented, where we discuss their structural, energetic aromaticity and IR and UV spectral features. In particular, we focus on the analysis of the bonding and various non-bonded interactions between these protonated aromatics and water clusters. The strength of non-bonded interactions is quantified and correlated with their electron density profiles. Furthermore, insights into the interfacial interactions and stability of these complexes were obtained through non-covalent index and symmetry-adapted perturbation theory (SAPT0) analyses. We also discuss the effects of the extension of the π aromatic cloud on the water solvation of these protonated aromatics. In particular, we extended our predictions for the S0 → S1 and S0 → T1 wavelength transitions of micro hydrated H+PAHs to deduce those of these species solvated in aqueous solution. The present findings should be useful for understanding, at the microscopic level, the effects of water interacting with H+PAHs, which are relevant for organic chemistry, astrochemistry, atmospheric chemistry, combustion and materials science.

Formation of Eigen or Zundel Features at Protonated Water Cluster-Aromatic Interfaces

Interfacial interactions of protonated water clusters adsorbed at aromatic surfaces play an important role in biology, and in atmospheric, chemical and materials sciences. Here, we investigate the interaction of protonated water clusters ((H+ H2 O)n (where n=1-3)) with benzene (Bz), coronene (Cor) and dodecabenzocoronene (Dbc)). To study the structure, stability and spectral features of these complexes, computations are done using DFT-PBE0(+D3) and SAPT0 methods. These interactions are probed by AIM electron density topography and non-covalent interactions index (NCI) analyses. We suggest that the excess proton plays a crucial role in the stability of these model interfaces through strong inductive effects and the formation of Eigen or Zundel features. Also, computations reveal that the extension of the π-aromatic system and the increase of the number of water molecules in the H-bounded water network led to a strengthening of the interactions between the corresponding aromatic compound and protonated water molecules, except when a Zundel ion is formed. The present findings may serve to understand in-depth the role of proton localized at aqueous medium interacting with large aromatic surfaces such as graphene interacting with acidic liquid water. Besides, we give the IR and UV-Vis spectra of these complexes, which may help for their identification in laboratory.

Design of Ascorbic Acid Eutectic Mixtures With Sugars to Inhibit Oxidative Degradation

L-Ascorbic acid (ASC), commonly known as vitamin C, acts as an anti-oxidant in the biological system. It is extensively used as an excipient in pharmaceutical industry, food supplements in fruit juices, and food materials due to its free radicals scavenging activity. Main drawback of ASC is its poor aqueous stability owing to the presence of lactone moiety that is easily oxidized to dehydroascorbic acid and further degraded. To improve aqueous stability and inhibit oxidative degradation, ASC was co-crystallized to constitute binary eutectic compositions with mono and di-saccharides such as glucose, sucrose, lactose, and mannitol. The eutectics were confirmed by their (single) lower melting endotherm compared to ASC and sugars, although Powder X-ray diffraction (PXRD) and Fourier transform Infrared spectroscopy (FT-IR) data confirmed the characteristics of their physical mixture. Scanning electron microscope (SEM) images of the binary eutectics confirmed their irregular morphology. The ASC eutectics exhibited improved shelf-life by 2–5-fold in weakly acidic (pH 5) and neutral (pH 7) aqueous buffer medium, whereas the eutectic with glucose enhanced shelf-life only by 1.1–1.2-fold in acidic medium (pH 3.3 and 4). Notably, stabilizing effect of the sugar eutectics decreased with increasing acidity of the medium. In addition, higher binding energy of the disaccharide eutectics partially supports the aqueous stability order of ASC in the neutral pH medium due to more number of non-bonded interactions than that of monosaccharides.

A Benchmark Protocol for DFT Approaches and Data-Driven Models for Halide-Water Clusters

Dissolved ions in aqueous media are ubiquitous in many physicochemical processes, with a direct impact on research fields, such as chemistry, climate, biology, and industry. Ions play a crucial role in the structure of the surrounding network of water molecules as they can either weaken or strengthen it. Gaining a thorough understanding of the underlying forces from small clusters to bulk solutions is still challenging, which motivates further investigations. Through a systematic analysis of the interaction energies obtained from high-level electronic structure methodologies, we assessed various dispersion-corrected density functional approaches, as well as ab initio-based data-driven potential models for halide ion-water clusters. We introduced an active learning scheme to automate the generation of optimally weighted datasets, required for the development of efficient bottom-up anion-water models. Using an evolutionary programming procedure, we determined optimized and reference configurations for such polarizable and first-principles-based representation of the potentials, and we analyzed their structural characteristics and energetics in comparison with estimates from DF-MP2 and DFT+D quantum chemistry computations. Moreover, we presented new benchmark datasets, considering both equilibrium and non-equilibrium configurations of higher-order species with an increasing number of water molecules up to 54 for each F, Cl, Br, and I anions, and we proposed a validation protocol to cross-check methods and approaches. In this way, we aim to improve the predictive ability of future molecular computer simulations for determining the ongoing conflicting distribution of different ions in aqueous environments, as well as the transition from nanoscale clusters to macroscopic condensed phases.

Twisted Eigen Can Induce Proton Transfer at a Hydrophobic-Hydrophilic Interface

The investigation of proton localization at a hydrophobic-hydrophilic interface is an important problem in chemical and materials sciences. In this study, protonated benzene (i.e., benzenium ion) and water clusters [BZH⁺W_n (where n = 1-6)] are selected as prototype models to understand the interfacial interactions and proton transfer mechanism between a carbonaceous surface and water molecules. The excess protons can localize in the vicinity of the hydrophobic-hydrophilic interface, and these clusters are stabilized by various kinds of noncovalent interactions. Calculations are carried out using ab initio (MP2) and density functional theory B3LYP methods to shed more light on geometries, energetics, and spectral signatures of the protonated species [H⁺(H₂O)_n] at the interfaces. These calculations revealed few low-lying isomers, which have not been reported earlier. Scrutiny of the results reveals that proton localization in the hydrophilic environment is more stable than the hydrophobic benzene π-cloud. Furthermore, the occurrence of an O-H⁺···π hydrogen bond significantly influences the O-H⁺···O interactions in the water clusters and also intensively affects the vibrational modes of the Eigen cation. Thus, the aromatic π-clouds can stabilize the Eigen cation and at the same time, a twisted form of Eigen (one O-H⁺···π → two O-H⁺···π) can enhance the proton transfer through the water chain via a Grotthuss-type mechanism. The vibrational spectra of these clusters reveal that there is a large red-shifted frequency for the O-H⁺···O, O-H⁺···π, and O-H···π modes of interaction. The energetic values and vibrational frequencies obtained from the B3LYP method are in close agreement with the MP2 level and experimental values, respectively.

Combined DFT and MD simulation studies of protein stability on imidazolium–water (ImH+Wn) clusters with aromatic amino acids

The hydrated clusters of protonated imidazole (ImH⁺) can induce protein denaturation through various kinds of monovalent interactions such as cation···π (stacking), N–H⋯π (T-shaped) and water-mediated O–H⋯O H-bonds.

Multicomponent solid forms of the uric acid reabsorption inhibitor lesinurad and cocrystal polymorphs with urea: DFT simulation and solubility study

Lesinurad (systematic name: 2-{[5-bromo-4-(4-cyclopropylnaphthalen-1-yl)-4H-1,2,4-triazol-3-yl]sulfanyl}acetic acid, C₁₇H₁₄BrN₃O₂S) is a selective uric acid reabsorption inhibitor related to gout, which exhibits poor aqueous solubility. High-throughput solid-form screening was performed to screen for new solid forms with improved pharmaceutically relevant properties. During polymorph screening, we obtained two solvates with methanol (CH₃OH) and ethanol (C₂H₅OH). Binary systems with caffeine (systematic name: 3,7-dihydro-1,3,7-trimethyl-1H-purine-2,6-dione, C₈H₁₀N₄O₂) and nicotinamide (C₆H₆N₂O), polymorphs with urea (CH₄N₂O) and eutectics with similar drugs, like allopurinol and febuxostat, were prepared using the crystal engineering approach. All these novel solid forms were confirmed by XRD, DSC and FT-IR. The crystal structures were solved by single-crystal and powder X-ray diffraction. The crystal structures indicate that the lesinurad molecule is highly flexible and the triazole moiety, along with the rotatable thioacetic acid (side chain) and cyclopropane ring, is almost perpendicular to the planar naphthalene moiety. The carboxylic acid-triazole heterosynthon in the drug is interrupted by the presence of methanol and ethanol molecules in their crystal structures and forms intermolecular macrocyclic rings. The caffeine cocrystal maintains the consistency of the acid-triazole heterosynthons as in the drug and, in addition, they are bound by several auxiliary interactions. In the binary system of nicotinamide and urea, the acid-triazole heterosynthon is replaced by an acid-amide synthon. Among the urea cocrystal polymorphs, Form I (P-1, 1:1) consists of an acid-amide (urea) heterodimer, whereas in Form II (P2₁/c, 2:2), both acid-amide heterosynthons and urea-urea dimers co-exist. Density functional theory (DFT) calculations further support the experimentally observed synthon hierarchies in the cocrystals. Aqueous solubility experiments of lesinurad and its binary solids in pH 5 acetate buffer medium indicate the apparent solubility order lesinurad-urea Form I (43-fold) > lesinurad-caffeine (20-fold) > lesinurad-allopurinol (12-fold) ≃ lesinurad-nicotinamide (11-fold) > lesinurad, and this order is correlated with the crystal structures.

Theoretical Insights into the Role of Water Molecules in the Guanidinium-Based Protein Denaturation Process in Specific to Aromatic Amino Acids

Noncovalent interactions between the guanidinium cation (Gdm⁺) and aromatic amino acids (AAs) in the water molecules have been studied using quantum chemical calculation and molecular dynamics (MD) simulations. Our studies show that there are two different modes of interactions between Gdm⁺ and AAs with and without water molecules. It is observed that nonhydrated Gdm⁺ interacts with AAs through N-H···π interactions, whereas hydrated clusters of Gdm⁺ are stabilized by stacking interactions with the help of the water-mediated hydrogen bond. Thus, different hydration patterns have significant effects on the predominant cation···π interactions in AAs-Gdm⁺ complexes. Findings from MD simulation elicit that the interaction pattern of Gdm⁺ with AAs varies as Phe < Tyr < Trp. Both the QM and MD calculations show a similar trend in the interaction of AAs with Gdm⁺. Moreover, the interaction of AAs with Gdm⁺ depends on the spatial orientation of AAs in the protein and the concomitant local structure, that is, the AAs present in the unstructured region of protein such as coils and bends exhibit higher binding for Gdm⁺ when compared to the AAs present in the structured region of the protein such as the α-helix and the β-sheet. Our study clearly reveals that H-bonded water molecules and the hydration pattern of Gdm⁺ as well as the positional presence of these AAs in the protein structure context play determining roles in the denaturation of protein by the Gdm⁺ cation.

Insights on the interaction of Zn2+ cation with triazoles: Structures, bonding, electronic excitation and applications

At present, we investigate the structures, the stability, the bonding and the spectroscopy of the Zn²⁺-triazole complexes (Zn²⁺-Tz), which are subunits of triazolate based porous materials and Zn-enzymes. This theoretical work is performed using ab initio methods and density functional theory (DFT) where dispersion correction is included. Through these benchmarks, we establish the ability and reliability of M05-2X+D3 and PBE0+D3 functionals for the correct description of Zn²⁺-Tz bond since these DFTs lead to close agreement with post Hartree-Fock methods. Therefore, M05-2X+D3 and PBE0+D3 functionals are recommended for the characterization of larger organometallic complexes formed by Zn and N-rich linkers. For Zn²⁺-Tz, we found two stable σ-type complexes: (i) a planar structure where Zn²⁺ links to unprotonated nitrogen and (ii) an out-of-plane cluster where carbon interacts with Zn²⁺. The most stable isomers consist on a coordinated covalent bond between the lone pair of unprotonated nitrogen and the vacant 4s orbital of Zn²⁺. The roles of covalent interactions within these complexes are discussed after vibrational, NBO, NPA charges and orbital analyses. The bonding is dominated by charge transfer from Zn²⁺ to Tz and intramolecular charge transfer, which plays a vital role for the catalytic activity of these complexes. These findings are important to understand, at the microscopic level, the structure and the bonding within triazolate based macromolecular porous materials and Zn-enzymes.

Ab initio and density functional theory (DFT) studies on triflic acid with water and protonated water clusters

The structure, stability and infrared spectral signatures of triflic acid (TA) with water clusters (W_n) and protonated water clusters (TAH⁺W_n, n = 1 - 6) were computed using DFT and MP2 methods. Our calculations show that a minimum of three water molecules are necessary to stabilize the dissociated zwitterionic form of TA. It can be seen from the results that there is no significant movement of protons in smaller (n = 1 and 2) and linear (n = 1 - 6) types of water clusters. Further, the geometries of TAW_n clusters first form a neutral pair (NP) to contact ion pair (CIP), then form a solvent separated ion pair (SSIP) in a water hexamer. These findings reveal that proton transfer may take place through NP to CIP and then CIP to SSIP. The calculated binding energies (BEs) of ion pair clusters is always higher than that of NP clusters (i.e., more stable than the NP). Existing excess proton linear chain clusters transfer a proton to adjacent water molecules via a Grotthuss mechanism, whereas the same isomers in the branched motifs do not conduct protons. Examination of geometrical parameters and infrared frequencies reveals hydronium ion (H₃O⁺ also called Eigen cation) formation in both TAW_n and protonated TAW_n clusters. The stability of Eigen water clusters is three times higher than that of other non-Eigen water clusters. Our study shows clearly that formation of ion pairs in TAW_n and TAH⁺W_n clusters greatly favors proton transfer to neighboring water molecules and also enhances the stability of these complexes.

Microscopic investigations of site and functional selectivity of triazole for CO2 capture and catalytic applications

Ab initio and DFT studies on CO₂ interacting with different tautomers and isomers of triazole (TZ) are carried out to understand the adsorption mechanism, site selectivity and their mutual preferential attracting sites.

Structure, Spectroscopy, and Bonding within the Zn(q+)-Imidazole(n) (q = 0, 1, 2; n = 1-4) Clusters and Implications for Zeolitic Imidazolate Frameworks and Zn-Enzymes

Using density functional theory (DFT) with dispersion correction and ab initio post Hartree-Fock methods, we treat the bonding, the structure, the stability, and the spectroscopy of the complexes between Zn(q+) and imidazole (Im), Zn(q+)Imn (where q = 0, 1 and 2; n = 1-4). These entities are subunits of zeolitic imidazolate frameworks (ZIFs) and Zn-enzymes, which possess relevant roles in industrial and biological domains, respectively. We also investigate the Imn (n = 2-4) clusters for comparison. For each species, we determine several new structures that were not found previously. Our calculations show a competition between atomic metal solvation, by either σ-type interactions or π-stacking type interaction, and proton transfer through hydrogen bonding (H-bonding) in charged species. This results in several geometrical environments around the metal. These are connected with structural properties and the functional role of Zn cation within ZIFs and Zn-enzymes. Moreover, we show that the Zn(2+)Imn subunits do not absorb in the visible domain, which may be related to the photostability of ZIFs. Our findings are important for the development of new applications of ZIFs and metalloenzymes.

Carbon dioxide interaction with isolated imidazole or attached on gold clusters and surface: competition between σ H-bond and π stacking interaction

Using first principle methodologies, we investigate the subtle competition between σ H-bond and π stacking interaction between CO2 and imidazole either isolated, adsorbed on a gold cluster or adsorbed on a gold surface. These computations are performed using MP2 as well as dispersion corrected density functional theory (DFT) techniques. Our results show that the CO2 interaction goes from π-type stacking into σ-type when CO2 interacts with isolated imidazole and Au clusters or surface. The balance between both types of interactions is found when an imidazole is attached to a Au20 gold cluster. Thus, the present study has great significance in understanding and controlling the structures of weakly-bound molecular systems and materials, where hydrogen bonding and van der Waals interactions are competing. The applications are in the fields of the control of CO2 capture and scattering, catalysis and bio- and nanotechnologies.

Characterization of Znq+–imidazole (q = 0, 1, 2) organometallic complexes: DFT methods vs. standard and explicitly correlated post-Hartree–Fock methods

Benchmarking DFts for the characterization of the Zn^q+–imidazole (q= 0, 1, 2) complexes.

Role of size and shape selectivity in interaction between gold nanoclusters and imidazole: a theoretical study

We present a theoretical study on the structure, stability, spectra and electronic properties of imidazole (Im) adsorbed on gold nanoclusters (Aun, n = 2, 4, 6, 8, 10, and 20). These computations were performed using various density functional theories with and without inclusion of Grimme's (D3) dispersion correction. For small clusters, we also carried out wavefunction-based ab initio (MP2 and SCS-MP2) computations for comparison. Vibrational, atoms in molecules (AIM) and natural bond orbital (NBO) analyses clearly reveal the occurrence of charge transfer (CT) through covalent (N1-Au) and noncovalent interactions that play important roles in the stability of the Im@Aun complexes with anchor assisted H-bonds (Cα-H · Au). Therefore, gold clusters can act as H-bond acceptors with biomolecules for development of new materials and applications. Our study establishes also the ability and reliability of PBE0 and M05-2X functionals compared to B3LYP and PBE for an accurate description of covalent and noncovalent interactions between Im and gold clusters since they lead to close agreement with MP2. Finally, we show that the Au8 cluster may be viewed as large enough to mimic the 3D gold surface.