Exploring Hydration Patterns of Aldehydes and Amides: Ab Initio Investigations

Microhydration of a Carbonyl Group: How does the Molecular Electrostatic Potential (MESP) Impact the Formation of (H2O)n:(R2C═O)Complexes?

The presence of a carbonyl group in a molecule usually leads to the identification of a π-hole on the molecular electrostatic potential (MESP) of the species. How does this electrophilic site influence the formation of microhydrated complexes? To address this point, a panel of R₂CO solutes with various MESPs was selected, and we identified the structures and properties of several complexes containing one, two, three and six water molecules. The following solutes were considered in the present study: H₂CO, F₂CO, Cl₂CO,(NC)₂CO and H₂C═CO. Geometry optimizations and frequency calculations were carried out at the LC-ωPBE/6-311++G(d,p) level, with the GD3BJ empirical correction for dispersion. For a number of n water molecules around the R₂CO solute, the structure and the features of the most stable (H₂O)_n:(R₂CO) complexes are highly dependent on the MESP of the isolated R₂CO solute. The formation of pi-hole bondings appears to play a decisive role in the initiation of a three-dimensional organization of water molecules around the solute.

Molecular Hydration of Carbonic Acid: Ab Initio Quantum Chemical and Density Functional Theory Investigation

Molecular hydration of carbonic acid (H₂CO₃) is investigated in terms of bonding patterns in H₂CO₃···(H₂O) _n [ n = 1-4] hydrogen-bonded clusters within ab initio quantum chemical and density functional theory (DFT) frameworks. Successive addition of water molecules to H₂CO₃···H₂O entails elongation of O-H (hydroxyl) bond as well as contraction of specific intermolecular hydrogen bonds signifying hydration of carbonic acid; these structural features get markedly enhanced under the continuum solvation framework. A comparison between the structurally similar clusters H₂CO₃···(H₂O) _n and HCOOH···(H₂O) _n [ n = 1-3] brings out the structural stability of the former. The present investigation in conjunction with the binding energy behavior of approaching water molecule(s) should serve as a precursor for pathways exploring aqueous dissociation of H₂CO₃ for larger clusters, as well as development of force-field potentials for acid dissociation process.

Modelling water diffusion in plasticizers: development and optimization of a force field for 2,4-dinitroethylbenzene and 2,4,6-trinitroethylbenzene

A classical all-atom force field has been developed for 2,4,6-trinitroethylbenzene and 2,4-dinitroethylbenzene and applied in molecular dynamics simulations of the two pure and two mixed plasticizer systems. Bonding parameters and partial charges were derived through electronic and geometry optimization of the single molecules. The other required parameters were derived from values already available in the literature for generic nitro aromatic compounds, which were adjusted to reproduce to a high level of accuracy the densities of 2,4-dinitroethylbenzene, 2,4,6-trinitroethylbenzene and the energetic plasticizers K10 and R8002. This force field has been applied to both K10 and R8002, which when used as plasticizers form an energetic binder with nitrocellulose. Nitrocellulose decomposes in storage, under varying conditions, but in particular where it may become increasingly dry. Following the derivation of the force field, we have therefore applied it to calculate water diffusion coefficients for each of the different materials at 298 K and 338 K, thereby providing a starting point for understanding water behaviour in a nitrocellulose binder.

A classical all-atom force field has been developed for the plasticizer molecules 2,4,6-trinitroethylbenzene and 2,4-dinitroethylbenzene which can be used to investigate properties and energetic output of nitrocellulose-based propellants and binders.

Two faces of triel bonds in boron trihalide complexes

The N⋅⋅⋅B triel bonds in complexes of boron trihalides, BX₃ (X = F, Cl, Br, and I), with species acting as Lewis bases through the nitrogen center, NH₃ , N₂ , and HCN, are analyzed theoretically (MP2/aug-cc-pVTZ calculations). It is confirmed that stronger Lewis acid properties of the boron center are observed for the BCl₃ moiety than for the BF₃ one in complexes with the strong Lewis base (NH₃ ); while the opposite order is observed for complexes with the weak Lewis base (N₂ ). The BX₃ NCH complexes (for X = Cl, Br, and I) are characterized by two tautomeric forms and by two corresponding N⋅⋅⋅B distances, the shorter one possesses characteristics of the covalent bond. In a case of the BF₃ NCH complex one energetic minimum is observed. Ab initio calculations are supported by an analysis of molecular electrostatic potentials (EPs) and electron density distributions. The quantum theory of 'atoms in molecules' and the decomposition of the energy of interaction are applied. The aforementioned acidity orders as well as the existence of two tautomers for some of complexes result partly from the electrostatic interactions' balance; the EP distribution is different for the BF₃ species than for the other BX₃ species where X = Cl, Br, and I. © 2017 Wiley Periodicals, Inc.

Theoretical investigations of structures and properties of Au and S substituted poly (2,3-di-(4′-hydroxyphenyl)-1,4-phenylene ethynyl) at different electric field intensities

Theoretical investigations on Au and S substituted poly (2,3-di-(4′-hydroxyphenyl)-1,4-phenylene ethynyl) as a molecular wire have been carried out by incorporating the external electric field (EF). The results demonstrate that both geometric and electronic structures of conjugated molecular wire are sensitive to EF. The Z component of the dipole moment and total dipole moment increases continuously since EF polarizes the molecule. Static EF modifies the frontline molecular orbitals and HLG. The I-V curve indicates the symmetrical trend respect to the original point.

Systematic Study of the Tautomerism of Uracil Induced by Proton Transfer. Exploration of Water Stabilization and Mutagenicity

To systematically investigate all the possible tautomerisms from uracil (U) and its enol form (U) induced by proton transfer, we describe a study of structural tautomer interconversion in the gas phase, in a continuum solvent, and in a microhydrated environment with 1 or 2 explicit water molecules, using density functional theory (DFT) calculations by means of the B3LYP exchange and correlation functions. A total of 62 geometries including 25 transition states were optimized, and the geometrical parameters have been discussed. Some rules of the configuration variation in tautomerization were summarized. The relative stabilities of all the tautomers were established. When a proton transfers from the di-keto form to the keto-enol form, water molecules in different regions show absolutely opposite effects: some assist, whereas others hinder the tautomerization. However, when a proton transfers from the keto-enol form to the di-enol form, water molecules in different regions show similar effects: the Gibbs free energy always increases and the activation energy always decreases. Additionally, some important factors that obviously affect the activation energy and Gibbs free energy were found and discussed in detail. The reasons that water molecules can assist or prevent the proton transfer were given. Furthermore, on the basis of our calculated results, we explain why it is hard to detect the di-enol form of uracil in general experiments.

Hydrogen Bonding in Phenol, Water, and Phenol−Water Clusters

Structure, stability, and hydrogen-bonding interaction in phenol, water, and phenol-water clusters have been investigated using ab initio and density functional theoretical (DFT) methods and using various topological features of electron density. Calculated interaction energies at MP2/6-31G level for clusters with similar hydrogen-bonding pattern reveal that intermolecular interaction in phenol clusters is slightly stronger than in water clusters. However, fusion of phenol and water clusters leads to stability that is akin to that of H(2)O clusters. The presence of hydrogen bond critical points (HBCP) and the values of rho(r(c)) and nabla(2)rho(r(c)) at the HBCPs provide an insight into the nature of closed shell interaction in hydrogen-bonded clusters. It is shown that the calculated values of total rho(r(c)) and nabla(2)rho(r(c)) of all the clusters vary linearly with the interaction energy.

Many-body interaction analysis: Algorithm development and application to large molecular clusters

A completely automated algorithm for performing many-body interaction energy analysis of clusters (MBAC) [M. J. Elrodt and R. J. Saykally, Chem. Rev. 94, 1975 (1994); S. S. Xantheas, J. Chem. Phys. 104, 8821 (1996)] at restricted Hartree-Fock (RHF)/MA Plesset 2nd order perturbation theory (MP2)/density functional theory (DFT) level of theory is reported. Use of superior guess density matrices (DM's) for smaller fragments generated from DM of the parent system and elimination of energetically insignificant higher-body combinations, leads to a more efficient performance (speed-up up to 2) compared to the conventional procedure. MBAC approach has been tested out on several large-sized weakly bound molecular clusters such as (H(2)O)(n), n=8, 12, 16, 20 and hydrated clusters of amides and aldehydes. The MBAC results indicate that the amides interact more strongly with water than aldehydes in these clusters. It also reconfirms minimization of the basis set superposition error for large cluster on using superior quality basis set. In case of larger weakly bound clusters, the contributions higher than four body are found to be repulsive in nature and smaller in magnitude. The reason for this may be attributed to the increased random orientations of the interacting molecules separated from each other by large distances.