Effect of quantum nuclear motion on hydrogen bonding

This work considers how the properties of hydrogen bonded complexes, X-H center dot center dot center dot Y, are modified by the quantum motion of the shared proton. Using a simple two-diabatic state model Hamiltonian, the analysis of the symmetric case, where the donor (X) and acceptor (Y) have the same proton affinity, is carried out. For quantitative comparisons, a parametrization specific to the O-H center dot center dot center dot O complexes is used. The vibrational energy levels of the one-dimensional ground state adiabatic potential of the model are used to make quantitative comparisons with a vast body of condensed phase data, spanning a donor-acceptor separation (R) range of about 2.4-3.0 angstrom, i.e., from strong to weak hydrogen bonds. The position of the proton (which determines the X-H bond length) and its longitudinal vibrational frequency, along with the isotope effects in both are described quantitatively. An analysis of the secondary geometric isotope effect, using a simple extension of the two-state model, yields an improved agreement of the predicted variation with R of frequency isotope effects. The role of bending modes is also considered: their quantum effects compete with those of the stretching mode for weak to moderate H-bond strengths. In spite of the economy in the parametrization of the model used, it offers key insights into the defining features of H-bonds, and semi-quantitatively captures several trends. (C) 2014 AIP Publishing LLC.

Explosives Sensing by Using Electron-Rich Supramolecular Polymers: Role of Intermolecular Hydrogen Bonding in Significant Enhancement of Sensitivity

We demonstrate here that supramolecular interactions enhance the sensitivity towards detection of electron-deficient nitro-aromatic compounds (NACs) over discrete analogues. NACs are the most commonly used explosive ingredients and are common constituents of many unexploded landmines used during World WarII. In this study, we have synthesised a series of pyrene-based polycarboxylic acids along with their corresponding discrete esters. Due to the electron richness and the fluorescent behaviour of the pyrene moiety, all the compounds act as sensors for electron-deficient NACs through a fluorescence quenching mechanism. A Stern-Volmer quenching constant determination revealed that the carboxylic acids are more sensitive than the corresponding esters towards NACs in solution. The high sensitivity of the acids was attributed to supramolecular polymer formation through hydrogen bonding in the case of the acids, and the enhancement mechanism is based on an exciton energy migration upon excitation along the hydrogen-bond backbone. The presence of intermolecular hydrogen bonding in the acids in solution was established by solvent-dependent fluorescence studies and dynamic light scattering (DLS) experiments. In addition, the importance of intermolecular hydrogen bonds in solid-state sensing was further explored by scanning tunnelling microscopy (STM) experiments at the liquid-solid interface, in which structures of self-assembled monolayer of the acids and the corresponding esters were compared. The sensitivity tests revealed that these supramolecular sensors can even detect picric acid and trinitrotoluene in solution at levels as low as parts per trillion (ppt), which is much below the recommended permissible level of these constituents in drinking water.

Metallacyclocumulenes: A Theoretical Perspective on the Structure, Bonding, and Reactivity

CONSPECTUS: Transition metals help to stabilize highly strained organic fragments. Metallacycles, especially unsaturated ones, provide much variety in this area. We had a sustained interest in understanding new C-C bond formation reactions affected by binuclear transition metal fragments Cp2M. One such study led to the exploration of the bimetallic C-C cleavage and coupled complexes, where the acetylide ligands bridge two metal atoms. The underlying M-C interaction in these complexes inspired the synthesis of a five-membered cyclocumulene complex, which opened a new phase in organometallic chemistry. The metallacyclocumulene produces a variety of C-C cleavage and coupled products including a radialene complex. Group 4 metallocenes have thus unlocked a fascinating chemistry by stabilizing strained unsaturated C4 organic fragments in the form of five-membered metallacyclocumulenes, metallacyclopentynes, and metallacycloallenes. Over the years, we have carried out a comprehensive theoretical study to understand the unusual stability and reactivity of these metallacycles. The unique (M-C-beta) interaction of the internal carbon atoms with the metal atom is the reason for unusual stability of the metallacycles. We have also shown that there is a definite dependence of the C-C coupling and cleavage reactions on the metal of metallacyclocumulenes. It demonstrates unexpected reaction pathways for these reactions. Based on this understanding, we have predicted and unraveled the stabilization factors of an unusual four-membered metallacycloallene complex. Indeed, our prediction about a four-membered heterometallacycle has led to an interesting bonding situation, which is experimentally realized. This type of M-C bonding is intriguing from a fundamental perspective and has great relevance in synthesizing unusual structures with interesting properties. In this Account, we first give a short prologue of what led to the present study and describe the salient features of the structure and bonding of the metallacyclocumulenes. The unusual reaction pathway of this metallacycle is explored next. Similar features of the metallacyclopentynes and metallacycloallenes are briefly mentioned. Then, we discuss the exploitation of the unique M-C bonding to design some exotic molecules such as a four-membered metallacycloallene complex. Our efforts to build a conceptual framework to understand these metallacycles and to exploit their chemistry continue.

Hydrogen bonding, halogen bonding and lithium bonding: an atoms in molecules and natural bond orbital perspective towards conservation of total bond order, inter- and intra-molecular bonding

One hundred complexes have been investigated exhibiting D-X center dot center dot center dot A interactions, where X = H, Cl or Li and DX is the `X bond' donor and A is the acceptor. The optimized structures of all these complexes have been used to propose a generalized `Legon-Millen rule' for the angular geometry in all these interactions. A detailed Atoms in Molecules (AIM) theoretical analysis confirms an important conclusion, known in the literature: there is a strong correlation between the electron density at the X center dot center dot center dot A bond critical point (BCP) and the interaction energy for all these interactions. In addition, we show that extrapolation of the fitted line leads to the ionic bond for Li-bonding (electrostatic) while for hydrogen and chlorine bonding, it leads to the covalent bond. Further, we observe a strong correlation between the change in electron density at the D-X BCP and that at the X center dot center dot center dot A BCP, suggesting conservation of the bond order. The correlation found between penetration and electron density at BCP can be very useful for crystal structure analysis, which relies on arbitrary van der Waals radii for estimating penetration. Various criteria proposed for shared-and closed-shell interactions based on electron density topology have been tested for H/Cl/Li bonded complexes. Finally, using the natural bond orbital (NBO) analysis it is shown that the D-X bond weakens upon X bond formation, whether it is ionic (DLi) or covalent (DH/DCl) and the respective indices such as ionicity or covalent bond order decrease. Clearly, one can think of conservation of bond order that includes ionic and covalent contributions to both D-X and X center dot center dot center dot A bonds, for not only X = H/Cl/Li investigated here but also any atom involved in intermolecular bonding.

Dissociation Constants of Weak Acids from ab Initio Molecular Dynamics Using Metadynamics: Influence of the Inductive Effect and Hydrogen Bonding on pK(a) Values

The theoretical estimation of the dissociation constant, or pK(a), of weak acids continues to be a challenging field. Here, we show that ab initio CarParrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free-energy profile of the dissociation reaction provide reasonable estimates of the pK(a) value. Water molecules, sufficient to complete the three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. The free-energy profiles exhibit two distinct minima corresponding to the dissociated and neutral states of the acid, and the difference in their values provides the estimate for pK(a). We show for a series of organic acids that CPMD simulations in conjunction with metadynamics can provide reasonable estimates of pK(a) values. The acids investigated were aliphatic carboxylic acids, chlorine-substituted carboxylic acids, cis- and trans-butenedioic acid, and the isomers of hydroxybenzoic acid. These systems were chosen to highlight that the procedure could correctly account for the influence of the inductive effect as well as hydrogen bonding on pK(a) values of weak organic acids. In both situations, the CPMD metadynamics procedure faithfully reproduces the experimentally observed trend and the magnitudes of the pK(a) values.

Herringbone to cofacial solid state packing via H-bonding in diketopyrrolopyrrole (DPP) based molecular crystals: influence on charge transport

Themono-alkylation of DPP derivatives leads to cofacial pi-pi stacking via H-bonding unlike their di-alkylated counterparts, which exhibit a classical herringbone packing pattern. Single crystal organic field-effect transistor (OFET) measurements reveal a significant enhancement of charge carrier mobility for mono-hexyl DPP derivatives.

Stacking Interactions in RNA and DNA: Roll-Slide Energy Hyperspace for Ten Unique Dinucleotide Steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by B97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 134-147, 2015.

Experimental validation of `pnicogen bonding' in nitrogen by charge density analysis

The participation of a nitrogen atom acting as an electrophile in pnicogen bonding, a hitherto unexplored interaction has been established by experimental charge density analysis. QTAIM and NBO analyses ratify this observation.

Ionic Bonding between Artificial Atoms

Conventional solids are prepared from building blocks that are conceptually no larger than a hundred atoms. While van der Waals and dipole-dipole interactions also influence the formation of these materials, stronger interactions, referred to as chemical bonds, play a more decisive role in determining the structures of most solids. Chemical bonds that hold such materials together are said to be ionic, covalent, metallic, dative, or otherwise a combination of these. Solids that utilize semiconductor nanocrystal quantum dots as building units have been demonstrated to exist; however, the interparticle forces in such materials are decidedly not chemical. Here we demonstrate the formation of charge transfer states in a binary quantum dot mixture. Charge is observed to reside in quantum confined states of one of the participating quantum dots. These interactions lead to materials that may be regarded as the nanoscale analog of an ionic solid. The process by which these materials form has interesting parallels to chemical reactions in conventional chemistry.

Halogen Bonding Controls the Regioselectivity of the Deiodination of Thyroid Hormones and their Sulfate Analogues

The type1 iodothyronine deiodinase (1D-1) in liver and kidney converts the L-thyroxine (T4), a prohormone, by outer-ring (5) deiodination to biologically active 3,3,5-triiodothyronine (T3) or by inner-ring (5) deiodination to inactive 3,3,5-triiodothronine (rT3). Sulfate conjugation is an important step in the irreversible inactivation of thyroid hormones. While sulfate conjugation of the phenolic hydroxyl group stimulates the 5-deiodination of T4 and T3, it blocks the 5-deiodination of T4. We show that thyroxine sulfate (T4S) undergoes faster deiodination as compared to the parent thyroid hormone T4 by synthetic selenium compounds. It is also shown that ID-3 mimics, which are remarkably selective to the inner-ring deiodination of T4 and T3, changes the selectivity completely when T4S is used as a substrate. From the theoretical investigations, it is observed that the strength of halogen bonding increases upon sulfate conjugation, which leads to a change in the regioselectivity of ID-3 mimics towards the deiodination of T4S. It has been shown that these mimics perform both the 5- and 5-ring deiodinations by an identical mechanism.

Role of synergistic pi-pi stacking and X-H center dot center dot center dot Cl (X = C, N, O) H-bonding interactions in gelation and gel phase crystallization

The self-assembly of p-pyridyl-ended oligo-p-phenylenevinylenes (OPVs) in ethanol leads to the formation of either hollow or solid microrods. The corresponding protonated OPVs with n-butyl chains induce transparent gelation and also gel phase crystallization owing to various synergistic noncovalent interactions. The chloride ion-selective gelation, AIEE and stimuli responsiveness of the gel are also observed.

Effects of hydrogen bonding on amide-proton chemical shift anisotropy in a proline-containing model peptide

Longitudinal relaxation due to cross-correlation between dipolar ((HN-1H alpha)-H-1) and amide-proton chemical shift anisotropy (H-1(N) CSA) has been measured in a model tripeptide Piv-(L)Pro-(L)Pro-(L)Phe-OMe. The peptide bond across diproline segment is known to undergo cis/trans isomerization and only in the cis form does the lone Phe amide-proton become involved in intramolecular hydrogen bonding. The strength of the cross correlated relaxation interference is found to be significantly different between cis and trans forms, and this difference is shown as an influence of intramolecular hydrogen bonding on the amide-proton CSA. (C) 2015 Elsevier B.V. All rights reserved.

Conformational Stability and Intramolecular Hydrogen Bonding in 1,2-Ethanediol and 1,4-Butanediol

The gas-phase infrared spectra of 1,2-ED and 1,4-BD have been, recorded at three different temperatures using a multipass gas cell of 6 m optical path length. DFT calculation has also been carried out using 6-311++G** and aug-cc-pVDZ basis sets to look for the existence of intramolecular hydrogen bonding, in them from the red shift and infrared absorption intensity enhancement of the bonded O-H band compared to that of the free O-H hand. Equilibrium population analysis With 10 conformers of 1,2-ED and 1,4-BD at experimental temperatures were-carried out for the reconstruction of the Observed vibrational spectra at that temperature,using standard statistical relationships. The most abundant conformer at experimental temperatures, was identified. In 1,2-ED a red shift of 45 cm(-1) in the intramolecularly interacting O-H stretching vibrational band position and no significant intensity enhancement compared to that of-the free O-H have been observed. On the contrary, in one of the hydrogen bonded conformers of 1,4-BD, a 124 cm(-1) red shift in the O-H stretching frequency and 8.5 times-intensity enhancement for the ``bonded'' O-H compared to that of the ``free'' O-H is seen. On the basis of this comparative study, we have concluded that strong intramolecular hydrogen bonding exists in 1,4-BD. But there appears, to be weak intramolecular hydrogen bonding in 1,2-ED at temperatures of 303, 313, and 323 Km the gas phase We have found that most stable hydrogen-bonded conformers of 1,4-BD are less populated than some of the non-hydrogen-bonded conformers. Even for the 1,4-BD, the relative population of the g'GG'Gt conformer, which has a strong intramolecular,hydrogen bond, is less than what is predicted. Perhaps the intramolecular hydrogen bond plays a less Significant role in the relative stability,of the various Conformers than what has been predicted from calculations and prevails in the literature.

``Conformational Simulation'' of Sulfamethizole by Molecular Complexation and Insights from Charge Density Analysis: Role of Intramolecular S center dot center dot center dot O Chalcogen Bonding

Intramolecular S center dot center dot center dot O chalcogen bonding and its potential to lock molecular conformation have been examined in the crystal forms of sulfamethizole, a sulfonamide antibiotic. Molecular complexes of sulfamethizole, including salts and cocrystal, have been synthesized, and their crystal structures were analyzed in order to examine the possible conformational preferences of the molecule in various ionic states and supramolecular environments (neutral/cocrystal, anionic salt, and cationic salt forms). The electrostatic potential mapped on Hirshfeld surfaces generated for these crystal forms provides insights into the possible binding modes of the drug in different environments. Further, the observed conformation locking feature has been rationalized in terms of the experimental charge density features of the intramolecular S center dot center dot O chalcogen bonding in sulfamethizole. The study quantitatively illustrates and rationalizes an intriguing case of a local minimum of molecular conformation being exclusively preferred over the global minimum, as it facilitates more efficient intermolecular interactions in a supramolecular environment.

Halogen Bonding and Chalcogen Bonding in 4,7-Dibromo-5, 6-dinitro-2,1,3-benzothiadiazole

An organic solid, 4,7-dibromo-5,6-dinitro-2,1,3-benzothiadiazole, has been designed to serve as an illustrative example to quantitatively evaluate the relative merits of halogen and chalcogen bonding in terms of charge density features. The compound displays two polymorphic modifications, one crystallizing in a non-centrosymmetric space group (Z' = 1) and the other in a centrosymmetric space group with two molecules in the asymmetric unit (Z' = 2). Topological analysis based on QTAIM clearly brings out the dominance of the chalcogen bond over the halogen bond along with an indication that halogen bonds are more directional compared to chalcogen bonds. The cohesive energies calculated with the absence of both strong and weak hydrogen bonds as well as stacking interaction are indicative of the stabilities associated with the polymorphic forms.