Unexpected Strengthening of the H-Bond Complexes in a Polar Solvent Due to a More Efficient Solvation of the Complex Compared to Isolated Monomers

Solvation-Enhanced Salt Bridges

Salt bridges formed by amidines and carboxylic acids represent an important class of noncovalent interaction in biomolecular and supramolecular systems. Isothermal titration calorimetry was used to study the relationships between the strength of the interaction, the chemical structures of the components, and the nature of the solvent. The stability of the 1:1 complex formed in chloroform changed by 2 orders of magnitude depending on the basicity of the amidine and the acidity of the acid, which is consistent with proton transfer in the complex. Polar solvents reduce the stabilities of salt bridges formed with N,N'-dialkylamidines by up to 3 orders of magnitude, but this dependence on solvent polarity can be eliminated if the alkyl groups are replaced by protons in the parent amidine. The enhanced stability of the complex formed by benzamidine is due to solvation of the NH sites not directly involved in salt bridge formation, which become significantly more polar when proton transfer takes place, leading to more favorable interactions with polar solvents in the bound state. Calculation of H-bond parameters using density functional theory was used to predict solvent effects on the stabilities of salt bridges to within 1 kJ mol-1. While H-bonding interactions are strong in nonpolar solvents, and solvophobic interactions are strong in polar protic solvents, these interactions are weak in polar aprotic solvents. In contrast, amidinium-carboxylate salt bridges are stable in both polar and nonpolar aprotic solvents, which is attractive for the design of supramolecular systems that operate in different solvent environments.

The Stability of Hydrogen-Bonded Ion-Pair Complex Unexpectedly Increases with Increasing Solvent Polarity

The generally observed decrease of the electrostatic energy in the complex with increasing solvent polarity has led to the assumption that the stability of the complexes with ion-pair hydrogen bonds decreases with increasing solvent polarity. Besides, the smaller solvent-accessible surface area (SASA) of the complex in comparison with the isolated subsystems results in a smaller solvation energy of the latter, leading to a destabilization of the complex in the solvent compared to the gas phase. In our study, which combines Nuclear Magnetic Resonance, Infrared Spectroscopy experiments, quantum chemical calculations, and molecular dynamics (MD) simulations, we question the general validity of this statement. We demonstrate that the binding free energy of the ion-pair hydrogen-bonded complex between 2-fluoropropionic acid and n-butylamine (CH3CHFCOO-…NH3But+) increases with increased solvent polarity. This phenomenon is rationalized by a substantial charge transfer between the subsystems that constitute the ion-pair hydrogen-bonded complex. This unexpected finding introduces a new perspective to our understanding of solvation dynamics, emphasizing the interplay between solvent polarity and molecular stability within hydrogen-bonded systems.

Impact of dielectric constant of solvent on the formation of transition metal-ammine complexes

The DFT-level computational investigations into Gibbs free energies (ΔG) demonstrate that as the dielectric constant of the solvent increases, the stabilities of [M(NH3 )n ]2+/3+ (n = 4, 6; M = selected 3d transition metals) complexes decrease. However, there is no observed correlation between the stability of the complex and the solvent donor number. Analysis of the charge transfer and Wiberg bond indices indicates a dative-bond character in all the complexes. The solvent effect assessed through solvation energy is determined by the change in the solvent accessible surface area (SASA) and the change in the charge distribution that occurs during complex formation. It has been observed that the SASA and charge transfer are different in the different coordination numbers, resulting in a variation in the solvent effect on complex stability in different solvents. This ultimately leads to a change between the relative stability of complexes with different coordination numbers while increasing the solvent polarity for a few complexes. Moreover, the findings indicate a direct relationship between ΔΔG (∆Gsolvent -∆Ggas ) and ΔEsolv , which enables the computation of ΔG for the compounds in a particular solvent using only ΔGgas and ΔEsolv . This approach is less computationally expensive.

Trends in the stability of covalent dative bonds with variable solvent polarity depend on the charge transfer in the Lewis electron-pair system

In general, the stability of neutral complexes with dative bonds increases as the polarity of the solvent increases. This is based on the fact that the dipole moment of the complex increases as the charge transferred from the donor to the acceptor increases. As a result, the solvation energy of the complex becomes greater than that of subsystems, causing an increase in the stabilization energy with increasing solvent polarity. Our research confirms this assumption, but only when the charge transfer is sufficiently large. If it is below a certain threshold, the increase in the complex's dipole moment is insufficient to result in a higher solvation energy than subsystems. Thus, the magnitude of the charge transfer in the Lewis electron-pair system determines the stability trends of dative bonds with varying solvent polarity. We used molecular dynamics (MD) simulations based on an explicit solvent model, which is considered more reliable, to verify the results obtained with a continuous solvent model.

Exploring the Non-Covalent Bonding in Water Clusters

QTAIM and source function analysis were used to explore the non-covalent bonding in twelve different water clusters (H2O)n obtained by considering n = 2–7 and various geometrical arrangements. A total of seventy-seven O−H⋯O hydrogen bonds (HBs) were identified in the systems under consideration, and the examination of the electron density at the bond critical point (BCP) of these HBs revealed the existence of a great diversity of O−H⋯O interactions. Furthermore, the analysis of quantities, such as |V(r)|/G(r) and H(r), allowed a further description of the nature of analogous O−H⋯O interactions within each cluster. In the case of 2-D cyclic clusters, the HBs are nearly equivalent between them. However, significant differences among the O−H⋯O interactions were observed in 3-D clusters. The assessment of the source function (SF) confirmed these findings. Finally, the ability of SF to decompose the electron density (ρ) into atomic contributions allowed the evaluation of the localized or delocalized character of these contributions to ρ at the BCP associated to the different HBs, revealing that weak O−H⋯O interactions have a significant spread of the atomic contributions, whereas strong interactions have more localized atomic contributions. These observations suggest that the nature of the O−H⋯O hydrogen bond in water clusters is determined by the inductive effects originated by the different spatial arrangements of the water molecules in the studied clusters.