How can we detect hydrogen bond local cooperativity in liquid water: A simulation study

A Classical Model for 3-body Interactions in Aqueous Ionic Systems

We present a classical electrostatic induction model to evaluate the 3-body Ion-Water-Water (I-W-W) and (W-W-W) interactions in aqueous ionic systems. The monatomic ions were described by a point charge and a dipole-dipole polarizability, while for the polyatomic ions distributed multipoles up to hexadecapole and distributed polarizabilities up to quadrupole-quadrupole were used. The accuracy of the classical model is benchmarked against an accurate dataset of 936 (I-W-W) and 2,184 (W-W-W) 3-body terms for 13 different monatomic and polyatomic cation and anion systems. The classical model shows excellent agreement with the reference MP2 and CCSD(T) 3-body energies. The Root-Mean-Square-Errors (RMSEs) for monatomic cations, monatomic anions, and polyatomic ions were 0.29 kcal/mol, 0.25 kcal/mol, and 0.12 kcal/mol, respectively. The corresponding RMSE for 1,744 CCSD(T)/aVTZ 3-body (W-W-W) energies, used to train MB-pol, was 0.12 kcal/mol. The accuracy of the classical model demonstrates that the 3-body term for aqueous ionic systems can be accurately modeled classically, without the need to fit to tens of thousands of high-level ab initio calculations. This approach provides a fast but accurate and efficient path towards modeling the 3-body effect in aqueous ionic systems that is fully transferable across systems with different ions.

Thermal conductivity temperature dependence of water confined in nanoporous silicon

Recently, it has been shown that high density nanoconfined water was the reason of the important enhancement of the effective thermal conductivity up to a factor of 50% of a nanoporous silicon filled with water. In this work, using molecular dynamics simulations, we further investigate the role of the temperatureT(from 285 to 360 K) on the thermal conductivity enhancement of nanohybrid porous silicon and water system. Furthermore, by studying and analysing several structural and dynamical parameters of the nanoconfined water, we give physical insights of the observed phenomena. Upon increasing the temperature of the system, the thermal conductivity of the hybrid system increases reaching a maximum forT= 300 K. With this article, we prove the existence of new heat flux channels between a solid matrix and a nanoconfined liquid, with clear signatures both in the radial distribution function, mean square displacements, water molecules orientation, hydrogen bond networks and phonon density of states.

Effects of molecular shape on alcohol aggregation and water hydrogen bond network behavior in butanol isomer solutions

Despite butanol isomers such as n-butanol, sec-butanol, isobutanol and tert-butanol having the same chemical formula, their liquid-liquid phase diagrams are distinct. That is, tert-butanol is miscible in water at all concentrations, while the other three butanol isomers are partially miscible under ambient conditions. The molecular shape of tert-butanol is close to globular and differs from the other three butanol molecules with a relatively long carbon chain. By performing molecular dynamics simulations and graph theoretical analysis of the four water-butanol isomer mixtures at varying concentrations, we show how distinct butanol aggregates are formed which depend upon the molecular shape and affect the water H-bond network structure and phase diagram in the binary liquid. The three butanol isomers of n-butanol, sec-butanol and isobutanol at concentrated solutions form chain-like alcohol aggregates, but tert-butanol forms small aggregates due to the distinct packing behavior caused by its globular molecular shape. By employing the graph theoretical analysis such as the degree distribution and the eigenvalue spectrum from the adjacency matrix in the graphical representation of the alcohol H-bond network, we show that the tert-butanol aggregates have a different morphological structure from that of the other three butanol isomers in aqueous solution. The graph theoretically distinct butanol aggregates are categorized into two groups, water-compatible and water-incompatible, depending upon the interaction between the alcohol and water molecules. Based upon our observations, we propose that the water-incompatible networks of n-butanol, sec-butanol and isobutanol aggregates do not change the water structure significantly, forming two separate liquid phases that are alcohol-rich and water-rich. However, the water-compatible network of tert-butanol aggregates has a considerable interaction with the water molecules and causes significant disruption of the water H-bond network, forming a homogeneous solution. Understanding the alcohol aggregation behavior and water structure in butanol-water mixtures provides a critical clue in appreciating fundamental issues such as miscibility and phase separation in aqueous solution systems.

Theoretical insight into the hybridization effect of donor and acceptor atoms on the cooperativity of C-H···N hydrogen bonds

In the present work, the influence of hybridization on cooperativity between C-H···N hydrogen bonds is theoretically investigated. Here, C₂H₆, C₂H₄, and C₂H₂ are considered as hydrogen bonding donor while NH₃, N₂H₄, N₂H₂, and N₂ act as the hydrogen bonding acceptor. The calculations are performed at MP2/aug-cc-pVTZ level. It is observed that the stability of systems is amplified as C(sp) > C(sp²) > C(sp³) and also N(sp³) > N(sp²) > N(sp). The role of interaction and deformation energies on the stability of the systems is examined. The results indicate the contribution of interaction energy is dominant in all complexes. The strength of C-H···N hydrogen bond is estimated using interaction energy. In agreement with cooperative energies, the C-H···N hydrogen bond is respectively weakened/strengthened in the triads containing C(sp) and C(sp²)/C(sp³) where two hydrogen bonds coexist. On the other hand, the C-H···N hydrogen bond is strengthened in the ternary systems including N(sp³) and N(sp²) while an opposite behavior is obtained in the triad having N(sp).

Mutual Influence of Pnicogen Bonds and Beryllium Bonds: Energies and Structures in the Spotlight

Pnicogen bonds, which are weak noncovalent interactions (NCIs), can be significantly modified by the presence of beryllium bonds, one of the strongest NCIs known. We demonstrate the importance of this influence by studying ternary complexes in which both NCIs are present, that is, the ternary complexes formed by a nitrogen base (NH₃, NHCH₂, and NCH), a phosphine (fluorophosphane, PH₂F) and a beryllium derivative (BeH₂, BeF₂, BeCl₂, BeCO₃, and BeSO₄). Energies, structures, and nature of the chemical bonding in these complexes are studied by means of ab initio computational methods. The pnicogen bond between the nitrogen base and the phosphine and the beryllium bond between the fluorine atom of fluorophosphane and the beryllium derivative show large cooperativity effects both on energies and geometries, with dissociation energies up to 296 kJ mol^-1 and cooperativity up to 104 kJ mol^-1 in the most strongly bound complex, CH₂HN:PH₂F:BeSO₄. In the complexes between the strongest nitrogen bases and the strongest beryllium donors, phosphorus-shared and phosphorus-transfer bonds are found.