On the Correlation between Bond-Length Change and Vibrational Frequency Shift in Hydrogen-Bonded Complexes: A Computational Study of Y···HCl Dimers (Y = N2, CO, BF)

Factors contributing to halogen bond strength and stretch or contraction of internal covalent bond

The halogen bond formed by a series of Lewis acids TF₃X (T = C, Si, Ge, Sn, Pb; X = Cl, Br, I) with NH₃ is studied by quantum chemical calculations. The interaction energy is closely mimicked by the depth of the σ-hole on the X atom as well as the full electrostatic energy. There is a first trend by which the hole is deepened if the T atom to which X is attached becomes more electron-withdrawing: C > Si > Ge > Sn > Pb. On the other hand, larger more polarizable T atoms are better able to transmit the electron-withdrawing power of the F substituents. The combination of these two opposing factors leaves PbF₃X forming the strongest XBs, followed by CF₃X, with SiF₃X engaging in the weakest bonds. The charge transfer from the NH₃ lone pair into the σ*(TX) antibonding orbital tends to elongate the covalent TX bond, and this force is largest for the heavier X and T atoms. On the other hand, the contraction of this bond deepens the σ-hole at the X atom, which would enhance both the electrostatic component and the full interaction energy. This bond-shortening effect is greatest for the lighter X atoms. The combination of these two opposing forces leaves the T-X bond contracting for X = Cl and Br, but lengthening for I.

Intermolecular interactions in the crystalline structure of some polyhalogenated Di- And triamino Pyridines: Spectroscopical perspectives

Supramolecular synthon is identified as a unit and provides important structural and energetic information in the study of organic crystals. However, the direct estimation of the supramolecular interaction remains challenging. In the present work six polyhalogenated di- or triamino pyridines were synthesised, their crystalline structure was characterised, and corresponding supramolecular synthons were studied using a combination of quantum mechanical calculations and FT-IR and Raman spectroscopy. Some distinctive features were identified especially for three vibrational normal modes (RNMs) related to the pyridine ring (viz. RNM₁, RNM₃ and RNM₇) in the vibrational spectra (FT-IR and Raman) of the solid samples, which are due to the supramolecular interactions, hydrogen bond (hb) in particular, according to the quantum mechanical calculations. The comparison between the IR and Raman spectra of experimental and simulated results indicates that the adjacent intermolecular hydrogen bonds between two same molecules extensively exist in the solid samples. Moreover, some quantitative correlation was established among the dimerisation energies for hb dimers (hb1 dimers for compounds 1 and 2), the ring structure defined by the distribution of the substituents and quantitative characteristics of the vibrational spectra, for instance, the splitting magnitudes for RNM₃(2) in IR spectra and the peak gap between RNM₁ and RNM₂ in Raman spectra.

Theoretical study on L-H+-L with identical donors: short strong hydrogen bond or not?

Short strong hydrogen bonds (SSHBs) play crucial role in many chemical processes. Recently, as the representative of SSHBs, [F-H-F]- was experimentally observed. [F-H-F]- has a symmetric structure, which can be described as a H+ acid shared by two terminal F- donors (F--H+-F-). To explore whether two identical donors are bound to result in SSHBs, we performed theoretical studies on a series of compounds (L-H+-L) with two identical electron donors (L corresponds to donors containing group 14, 15, 16 and 17 elements). The results show that identical donors do not definitely lead to SSHBs. Instead, typical hydrogen bonds also exist. We found that both electronegativity and basicity contribute to the patterns of hydrogen bonds, where more electronegative and weaker donors benefit to SSHBs. Besides, it was found that zero-point energies also respond to the hydrogen bonding systems. This systemic work is expected to provide more insights into SSHBs.

Theoretical investigation on the improper hydrogen bond in κ-carrabiose⋯Y (Y = HF, HCl, HBr, NH3, H2O, and H2S) complexes

The nature of H-bonds in κ-carrabiose⋯Y (Y = HF, HCl, HBr, NH₃, H₂O, and H₂S) complexes was studied. For this aim, the structure of isolated κ-carrabiose was optimized using three global hybrids functional: B3LYP, PBE0, and M06-2X combined with 6-311G** basis set. Subsequently, the κ-carrabiose in the presence of HF, HCl, HBr, NH₃, H₂O, and H₂S was optimized using the CBS-4 M method. NBO analyses were then carried out at the MP2/6-311G** level of theory. A particular interest was focused on C₍₁₈₎-H₍₃₄₎⋯Y bond. The results reveal that the C₍₁₈₎-H₍₃₄₎⋯Y bond is an improper H-bond since a significant contraction of C₍₁₈₎-H₍₃₄₎ was observed during the complexation leading to a significant blueshifted stretching frequency. The NBO analyses have shown that the formation of the improper H-bonds C₍₁₈₎-H₍₃₄₎⋯Y (Y = F, Cl, Br, N, O, and S) is principally due to the increase of the s-character of the hybrid orbital in carbon atom (rehybridization) in κ-carrabiose⋯Y complexes. Regarding the polarization, it was proved that more the H-bond center (carbon in C₍₁₈₎-H₍₃₄₎⋯Y) becomes less positive, the hydrogen more positive, and Y more negative; more the contraction of the C₍₁₈₎-H₍₃₄₎ bond is important. It was also confirmed for intramolecular H-bonds in κ-carrabiose⋯Y complexes that the rehybridization is responsible for H-bonds nature either proper or improper.

Interaction-Deletion: A Composite Energy Method for the Optimization of Molecular Systems Selectively Removing Specific Nonbonded Interactions

The complex interactions between different portions of a large molecule can be challenging to analyze through traditional electronic structure calculations. Moreover, standard methods cannot easily quantify the physical consequences of individual pairwise interactions inside a molecule. By creating a set of molecular fragments, we propose a composite energy method to explore changes in a molecule caused by removing selected nonbonded interactions between different molecular portions. Energies and forces are easily obtained with this composite approach, allowing geometry optimizations that lead to chemically meaningful structures that describe how the omitted interactions contribute to changes in the local geometrical minima. We illustrate the application of our new hybrid scheme by computing the influence of intramolecular hydrogen-bonding interactions in two small molecules: 1,6-(tG⁺G⁺TG⁺G⁺g^-)-hexanediol and a cyclic analogue, cis-1,4-cyclohexanediol. The resulting structural and energetic changes are interpreted to yield key physical insights and quantify concepts such as "preparation energy" or "reorganization energy". We demonstrate that the composite method can be extended to larger molecular systems by showing its application on a Si(100) surface model containing interactions between dissociated ammonia molecules on adjacent surface dimers. The scheme's efficacy is also tested by applying it to systems having multiple intramolecular interactions, viz., 3₁₀-polyglycine and H⁺GPGG. Furthermore, the cooperative nature of intramolecular hydrogen bonds is explored by using interaction-deletion in 2-nitrobenzene-1,3-diol.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Impact of donor acidity and acceptor anharmonicity on vc=o spectral shifts in O-H···O=C H-bonded ketone-alcohol complexes: An IR spectroscopic investigation

O-H···O=C Hydrogen bonding (H-bonding) results in spectral shifts in both ν_O-H and ν_C=O modes. A large number of investigations exist in literature that focuses on how the spectral shifts vary with certain properties of the donors and acceptors. However information on how the magnitude of spectral shift is dictated individually by the donor and acceptor is not yet clear to us. Here, IR spectroscopy in room temperature CCl₄ solution has been used to investigate how ν_C=O spectral shifts of ketones in H-bonded complexes with alcohols are influenced by change in donor and acceptor properties. For this purpose 25 number of O-H···O=C H-bonded complexes formed by 5 different ketones and 5 alcohols have been studied. The magnitude of red-shifts shown by the ν_C=O bands were found to show systematic trends with changing donor and acceptor species. It was found that for a particular ketone, the magnitude of shift increases monotonically with acidity of the alcohols. Spectral shifts were found to increase in a linear fashion with decreasing p^K_a of the donor alcohols. On the other hand, when the alcohol were kept fixed, the spectral shift was found to be dependent on the CO bond strength of the ketones. We found spectral shifts to linearly increase with increasing anharmonicity constant and decreasing dissociation energy of CO bond. Finally, it has been shown that there exist concomitant correlations of the spectral shifts with donor p^K_a and acceptor ν_C=O anharmonicity constant/bond dissociation energy. The relations have been validated for H-bonded complexes of 1,4-cyclohexanedione with the above mentioned five alcohols.

Catalyzed Atomic Layer Deposition of Silicon Oxide at Ultralow Temperature Using Alkylamine

We report the catalyzed atomic layer deposition (ALD) of silicon oxide using Si₂Cl₆, H₂O, and various alkylamines. The density functional theory (DFT) calculations using the periodic slab model of the SiO₂ surface were performed for the selection of alternative Lewis base catalysts with high catalytic activities. During the first half-reaction, the catalysts with less steric hindrance such as pyridine would be more effective than bulky alkylamines despite lower nucleophilicity. On the other hand, during the second half-reaction, the catalysts with a high nucleophilicity such as triethylamine (Et₃N) would be more efficient because the steric hindrance is less critical. The in situ process monitoring shows that the calculated atomic charge is a good indicator for expecting the catalyst activity in the ALD reaction. The use of Et₃N in the second half-reaction was essential to improving the growth rate as well as the step coverage of the film because the Et₃N-catalyzed process deposited a SiO₂ film with a step coverage of 98% that is better than 93% of the pyridine-catalyzed process. The adsorption of pyridine, ammonia (NH₃), or trimethylamine (Me₃N) salts was more favorable than that of Et₃N, n-Pr₃N, or ⁱPr₃N salts. Therefore, Et₃N was expected to incorporate less amine salts in the film as compared to pyridine, and the compositional analyses confirmed that the concentrations of Cl and N by the Et₃N-catalyzed process were significantly lower than those by the pyridine-catalyzed process.

Quantum Interference Contribution to the Dipole Moment of Diatomic Molecules

The interference energy partitioning analysis method developed by our group and used to study the nature of the chemical bond was extended to partition the electric dipole moment in quasi-classical and interference contributions. Our results show that interference participates in charge displacement in polar molecules, providing, directly or indirectly, a relevant contribution for the total dipole moment. A linear correlation was found between the interference contribution of the dipole moment from the bond electron group, μ_INT(bond), and the difference of electronegativity of the atoms which form the bond, ΔX_AB. This interesting result reinforces the fact that electronegativity is not a property of an atom alone, but rather a property of the atom in the molecule and that ΔX_AB can only be associated with that part of the total charge displacement resulting from the formation of the chemical bond. The partitioning of the total dipole moment into quasi-classical and interference contributions provides new insights about the reasons for the failure of the ΔX_AB criterion in predicting the correct orientation of the dipole moment in several molecules. The results of the present work also bring additional evidence for the previously proposed mechanism of formation of polar bonds.

Solvent-Independent Anharmonicity for Carbonyl Oscillators

The physical origins of vibrational frequency shifts have been extensively studied in order to understand noncovalent intermolecular interactions in the condensed phase. In the case of carbonyls, vibrational solvatochromism, MD simulations, and vibrational Stark spectroscopy suggest that the frequency shifts observed in simple solvents arise predominately from the environment's electric field due to the vibrational Stark effect. This is contrary to many previously invoked descriptions of vibrational frequency shifts, such as bond polarization, whereby the bond's force constant and/or partial nuclear charges are altered due to the environment, often illustrated in terms of favored resonance structures. Here we test these hypotheses using vibrational solvatochromism as measured using 2D IR to assess the solvent dependence of the bond anharmonicity. These results indicate that the carbonyl bond's anharmonicity is independent of solvent as tested using hexanes, DMSO, and D₂O and is supported by simulated 2D spectra. In support of the linear vibrational Stark effect, these 2D IR measurements are consistent with the assertion that the Stark tuning rate is unperturbed by the electric field generated by both hydrogen and non-hydrogen bonding environments and further extends the general applicability of carbonyl probes for studying intermolecular interactions.

Fingerprints of inter- and intramolecular hydrogen bonding in saligenin–water clusters revealed by mid- and far-infrared spectroscopy

Saligenin (2-(hydroxymethyl)phenol) exhibits both strong and weak intramolecular electrostatic interactions.

Theoretical investigation of gas-phase molecular complex formation between 2-hydroxy thiophenol and a water molecule

A systematic study of the interaction energies and hydrogen bonding interaction of a gas-phase molecular complex between 2-hydroxy thiophenol and a water molecule.

Vibrational solvatochromism. II. A first-principle theory of solvation-induced vibrational frequency shift based on effective fragment potential method

Vibrational solvatochromism is a solvation-induced effect on fundamental vibrational frequencies of molecules in solutions. Here we present a detailed first-principle coarse-grained theory of vibrational solvatochromism, which is an extension of our previous work [B. Błasiak, H. Lee, and M. Cho, J. Chem. Phys. 139(4), 044111 (2013)] by taking into account electrostatic, exchange-repulsion, polarization, and charge-transfer interactions. By applying our theory to the model N-methylacetamide-water clusters, solute-solvent interaction-induced effects on amide I vibrational frequency are fully elucidated at Hartree-Fock level. Although the electrostatic interaction between distributed multipole moments of solute and solvent molecules plays the dominant role, the contributions from exchange repulsion and induced dipole-electric field interactions are found to be of comparable importance in short distance range, whereas the charge-transfer effect is negligible. The overall frequency shifts calculated by taking into account the contributions of electrostatics, exchange-repulsion, and polarization terms are in quantitative agreement with ab initio results obtained at the Hartree-Fock level of theory.

Nature of the chemical bond and origin of the inverted dipole moment in boron fluoride: a generalized valence bond approach

The generalized product function energy partitioning (GPF-EP) method has been applied to investigate the nature of the chemical bond and the origin of the inverted dipole moment of the BF molecule. The calculations were carried out with GPF wave functions treating all of the core electrons as a single Hartree-Fock group and the valence electrons at the generalized valence bond perfect-pairing (GVB-PP) or full GVB levels, with the cc-pVTZ basis set. The results show that the chemical structure of both X (1)Σ(+) and a (3)Π states is composed of a single bond. The lower dissociation energy of the excited state is attributed to a stabilizing intraatomic singlet coupling involving the B 2sp-like lobe orbitals after bond dissociation. An increase of electron density on the B atom caused by the reorientation of the boron 2sp-like lobe orbitals is identified as the main responsible effect for the electric dipole inversion in the ground state of BF. Finally, it is shown that π back-bonding from fluorine to boron plays a minor role in the electron density displacement to the bonding region in both states. Moreover, this effect is associated with changes in the quasi-classical component of the electron density only and does not contribute to covalency in either of the states. Therefore, at least for the case of the BF molecule, the term back-bonding is misleading, since it does not contribute to the bond formation.

Evidence of blue-shifting N–H⋯N hydrogen bonding despite elongation of the N–H bond

N–H⋯N hydrogen bonding between pyrrole and N₂ results blue shift of ν_N–H accompanied by an increase in the N–H bond length.

Vibrational solvatochromism: towards systematic approach to modeling solvation phenomena

Vibrational solvatochromic frequency shift of IR probe is an effect of interaction between local electric field and IR probe in condensed phases. Despite prolonged efforts to develop empirical maps for vibrational frequency shifts and transition dipoles of IR probes, a systematic approach to ab initio calculation of vibrational solvatochromic charges and multipoles has not been developed. Here, we report on density functional theory (DFT) calculations of N-methylacetamide (NMA) frequency shifts using implicit and coarse-grained models. The solvatochromic infrared spectral shifts are estimated based on the distributed multipole analysis of electronic densities calculated for gas-phase equilibrium structure of NMA. Thus obtained distributed solvatochromic multipole parameters are used to calculate the amide I vibrational frequency shifts of NMA in water clusters that mimic the instantaneous configurations of the liquid water. Our results indicate that the spectral shifts are primarily electrostatic in nature and can be quantitatively reproduced using the proposed model with semi-quantitative accuracy when compared to the corresponding DFT results.

Ultrafast vibrational frequency shifts induced by electronic excitations: naphthols in low dielectric media

We study the solvent-induced frequency shifts of the OH-stretching mode of 1-naphthol and 2-naphthol in nonpolar/weakly polar solvents, subject to electronic excitation, with ultrafast UV/mid-infrared pump-probe spectroscopy and theoretical modeling based on Pullin's perturbative treatment of vibrational solvatochromic effects. The model is parametrized at the density functional theory (DFT) level, including the B3LYP/TZVP and TD-B3LYP/TZVP descriptions, for the naphthol chromophores in the S(0)- and (1)L(b)-states and accounts for both the static and the optical dielectric response of the solvent on time scales comparable to that of the OH-stretching vibrational motions. The favorable comparison between experimental and theoretical values of the solvent-induced vibrational frequency shifts suggests that the ultrafast dielectric response of the solvent contributes predominantly to the solvatochromic shifts in solvents of moderate polarity where specific solute-solvent interactions are absent.

The O-H stretching mode of a prototypical photoacid as a local dielectric probe

We investigate the OH stretch vibrational frequency shifts of a prototype photoacid, 2-naphthol (2N), when dissolved in solvents of low polarity. We combine femtosecond mid-infrared spectroscopy and a theoretical model based on the Pullin-van der Zwan-Hynes perturbative approach to explore vibrational solvatochromic effects in the ground S(0) and the first electronically excited (1)L(b) states. The model is parametrized using density functional theory (DFT), at the B3LYP/TZVP and TD-B3LYP/TZVP levels for the 2N chromophore in the S(0) and (1)L(b) states, respectively. From the agreement between experiment and theory we conclude that vibrational solvatochromic effects are dominated by the instantaneous dielectric response of the solvent, while time-dependent nuclear rearrangements are of secondary importance.

Weak H-bonds. Comparisons of CH···O to NH···O in proteins and PH···N to direct P···N interactions

Whereas CH···O H-bonds are usually weaker than interpeptide NH···O H-bonds, this is not necessarily the case within proteins. The nominally weaker CH···O are surprisingly strong, comparable to, and in some cases stronger than, the NH···O H-bonds in the context of the forces that hold together the adjacent strands in protein β-sheets. The peptide NH is greatly weakened as proton donor in certain conformations of the protein backbone, particularly extended structures, and forms correspondingly weaker H-bonds. The PH group is a weak proton donor, but will form PH···N H-bonds. However, there is a stronger interaction in which P can engage, in which the P atom, not the H, directly approaches the N electron donor to establish a direct P···N interaction. This approach is stabilized by the same sort of electron transfer from the N lone pair to the P-H σ* antibond that characterizes the PH···N H-bond.