Structures of water clusters in the solvent phase and relative stability compared to gas phase

Conformation Controlled Hydrogelation of Minimalistic α, γ Hybrid Peptide

A majority of short peptide (≤7 amino acids) hydrogels are primarily assembled via cross β-structure formation. In contrast to the natural trend, herein, we report the formation of supramolecular hydrogel from the ultrashort hybrid folded peptide composed of canonical α-amino acid and noncanonical γ-amino acid, Fmoc-γPhe-Phe-OH. The designed hybrid peptide hydrogel is composed of entangled fibers, has viscoelastic properties, exhibits proteolytic stability, and exhibits cytocompatibility with L929 fibroblast cells. Mutating the peptide sequence by altering the position of γPhe from the N-termini to C-termini transforms the self-assembly into crystalline aggregates. Combining FTIR, 2D NMR, and DFT calculations revealed that the hydrogel-forming peptide adopts a C9 H-bonded conformation, resembling the well-known γ-turn. However, the isomeric hybrid peptide adopts an extended structure. The present study highlights the importance of secondary structure in the higher order assembly of minimalist hybrid peptides and broadens the range of secondary structures to design short peptide-based hydrogels.

Solvation energies of the ferrous ion in water and in ammonia at various temperatures

CONTEXT

The solvation of metal ions is crucial to understanding relevant properties in physics, chemistry, or biology. Therefore, we present solvation enthalpies and solvation free energies of the ferrous ion in water and ammonia. Our results agree well with the experimental reports for the hydration free energy and hydration enthalpy. We obtained [Formula: see text] kJ mol[Formula: see text] for the hydration free energy and [Formula: see text] kJ mol[Formula: see text] for the hydration enthalpy of ferrous ion in water at room temperature. At ambient temperature, we obtained [Formula: see text] kJ mol[Formula: see text] as the [Formula: see text] ammoniation free energy and [Formula: see text] kJ mol[Formula: see text] for the ammoniation enthalpy. In addition, the free energy of solvation is deeply affected when the temperature increases. This pattern can be attributed to the rise of entropy when the temperature rises. Besides, the temperature does not affect the ammoniation enthalpies and the hydration enthalpy of the [Formula: see text] ion.

METHOD

All the geometry optimizations are performed at the MP2 methods associated with the 6-31++g(d,p) basis set of Pople. solvated phase structures of [Formula: see text] ion in water or in ammonia are performed using the PCM model. The [Formula: see text] program suite was used to perform all the calculations. The program TEMPO was also used to evaluate the temperature sensitivity of the different obtained geometries.

Hydration of [Formula: see text]aminobenzoic acid: structures and non-covalent bondings of aminobenzoic acid-water clusters

CONTEXT

Micro-hydration of the aminobenzoic acid is essential to understand its interaction with surrounding water molecules. Understanding the micro-hydration of the aminobenzoic acid is also essential to study its remediation from wastewater. Therefore, we explored the potential energy surfaces (PESs) of the para-aminobenzoic acid-water clusters, ABW[Formula: see text], [Formula: see text], to study the microsolvation of the aminobenzoic acid in water. In addition, we performed a quantum theory of atoms in molecules (QTAIM) analysis to identify the nature of non-covalent bondings in the aminobenzoic acid-water clusters. Furthermore, temperature effects on the stability of the located isomers have been examined. The located structures have been used to calculate the hydration free energy and the hydration enthalpy of the aminobenzoic acid using the cluster continuum solvation model. The hydration free energy and the hydration enthalpy of the aminobenzoic acid at room temperature are evaluated to be -7.0 kcal/mol and -18.1 kcal/mol, respectively. The hydration enthalpy is in perfect agreement with a previous experimental estimate. Besides, temperature effects on the calculated hydration enthalpy and free energy are reported. Finally, we calculated the gas phase binding energies of the most stable structures of the ABW[Formula: see text] clusters using twelve functionals of density functional theory (DFT), including empirical dispersion. The DFT functionals are benchmarked against the DLPNO-CCSD(T)/CBS. We have found that the three most suitable DFT functionals are classified in the following order: PW6B95D3 > MN15 > [Formula: see text]B97XD. Therefore, the PW6B95D3 functional is recommended for further study of the aminobenzoic acid-water clusters and similar systems.

METHODS

The exploration started with classical molecular dynamics simulations followed by complete optimization at the PW6B95D3/def2-TZVP level of theory. Optimizations are performed using Gaussian 16 suite of codes. QTAIM analysis is performed using the AIMAll program.

Structures of DMSO clusters and quantum cluster equilibrium (QCE)

Dimethylsulfoxide (DMSO) clusters are crucial for understanding processes in liquid DMSO. Despite its importance, DMSO clusters have received negligible attention due to the complexity of their potential energy surfaces (PESs). In this work, we explored the PESs of the DMSO clusters from dimer to decamer, starting with classical molecular dynamics, followed by full optimizations at the PW6B95-D3/def2-TZVP level of theory. In addition, the binding energies, the binding enthalpy per DMSO, and the quantum theory of atoms in molecules (QTAIM) analysis of the most stable isomers are reported. Temperature effects on the stability of the isomers have also been assessed. After thoroughly exploring the PESs of the DMSO clusters, 159 configurations have been used to apply the quantum cluster equilibrium (QCE) theory to liquid DMSO. The quantum cluster equilibrium theory has been applied to determine the liquid properties of DMSO from DMSO clusters. Thus, using the QCE, the population of the liquid DMSO, its infrared spectrum, and some thermodynamic properties of the liquid DMSO are predicted. The QCE results show that the population of the liquid DMSO is mainly dominated by the DMSO dimer and decamer, with the contribution in trace of the DMSO monomer, trimer, tetramer, pentamer, and octamer. More interestingly, the predicted infrared spectrum of liquid DMSO is in qualitative agreement with the experiment.

Solvation of Manganese(III) Ion in Water and in Ammonia

In this work, we have studied the solvation of manganese(III) ion in water and in ammonia using three levels of theory: MP2, MN15, and ωB97XD associated with the aug-cc-pVDZ basis set. The studied systems are constituted of Mn³⁺(H₂O)₆ and Mn³⁺(NH₃)₆ in gas and solvent phases as well as Mn³⁺(H₂O)₁₈ and Mn³⁺(NH₃)₁₈ in the gas phase. Four aspects of the solvation of manganese(III) ion have been examined for the aforementioned systems at the three levels of theory. First, we started by locating the Jahn-Teller elongated and compressed configuration in Mn³⁺(H₂O)₆ and Mn³⁺(NH₃)₆. Second, we calculated the spin state energies and the spin state free energies for temperatures ranging from 50 to 400 K to look at possible spin crossover in the studied systems. Third, we carried out a quantum theory of atoms in molecules (QTAIM) analysis, and we determined the ionic radii of manganese(III) ion in water and in ammonia. Fourth, we calculated the solvation free energies and the solvation enthalpies of manganese(III) ion in water and in ammonia using the cluster continuum solvation model. For these four aspects of the solvation of manganese(III) ion, most of the reported properties are provided in this work for the first time. We particularly found that the calculated solvation enthalpy of the manganese(III) ion in water is in good agreement with an experimental estimate.

OH-initiated degradation of 1,2,3-trimethylbenzene in the atmosphere and wastewater: Mechanisms, kinetics, and ecotoxicity

1,2,3-Trimethylbenzene (1,2,3-TMB) is an important volatile organic compound (VOC) present in petroleum wastewater and the atmosphere. This compound can be degraded by OH radicals via abstraction, addition and substitution mechanisms. Results show that the addition mechanism is dominant and H-abstraction is subdominant, while methyl abstraction and substitution mechanisms are negligible in the gas and aqueous phases. Moreover, H-abstraction products undergo further reactions with O₂, NO, NO₂, H₂O, and OH radicals in the atmosphere. Time-dependent density functional theory (TDDFT) calculations show that the degraded products, including 2,3,4-trimethylphenyl-nitroperoxoite, 1,2,3-trimethyl-4-nitrobenzene, 1,2,3-trimethyl-5-nitrobenzene, 2,6-dimethylbenzyl nitroperoxoite, 2,3-dimethylphenyl nitroperoxoite, 2,6-dimethylbenzaldehyde, and 2,3-dimethylbenzaldehyde, can photolyze under the sunlight. Kinetically, the calculated total rate constant is 5.57 × 10^-11 cm³ molecule^-1·s^-1 at 1 atm and 298 K, which is consistent with available experimental values measured in the atmosphere. In addition, the calculated total reaction rate constant in water is close to that in the gas phase. In terms of ecotoxicity, all degradation products are less toxic than the initial reactant to fish, green algae and daphnia. For mammals represented by rats, 1,2,3-TMB and its products are moderately toxic, except for 2,3-dimethylphenol and 2,6-dimethylphenol, which are slightly toxic.

Hydrogen bond networks of dimethylsulfoxide (DMSO) pentamer

Understanding of clusters of dimethylsulfoxide (DMSO) is important in several applications in Chemistry. Despite its importance, very few studies of DMSO clusters, (DMSO)_n, have been reported in comparison to systems such as water clusters or methanol clusters. In order to provide further understanding of DMSO clusters, we investigated the structures and non-covalent interactions of the (DMSO)_n, n=5. Therefore, the potential energy surface (PES) of the DMSO pentamer has been examined using classical molecular dynamics. The structures generated using classical molecular dynamics are further optimized at the PW6B95D3/aug-cc-pVDZ level of theory. To comprehend the non-covalent bondings in the DMSO pentamer, we carried out a quantum theory of atoms in molecule (QTAIM) analysis. In addition, the effects of temperature on the structural stability is investigated between 20 and 500K. It comes out that seven different kind of non-covalent bondings can be found in DMSO pentamers.

Adsorption free energy of phenol onto coronene: Solvent and temperature effects

Molecular modeling can considerably speed up the discovery of materials with high adsorption capacity for wastewater treatment. Despite considerable efforts in computational studies, the molecular modeling of adsorption processes has several limitations in reproducing experimental conditions. Handling the environmental effects (solvent effects) and the temperature effects are part of the important limitations in the literature. In this work, we address these two limitations using the adsorption of phenol onto coronene as case study. In the proposed model, for the solvent effects, we used a hybrid solvation model, with n explicit water molecules and implicit solvation. We increasingly used n=1 to n=12 explicit water molecules. To account for the temperature effects, we evaluated the adsorption efficiency using the adsorption free energy for temperatures varying from 200 to 400K. We generated initial configurations using classical molecular dynamics, before further optimisation at the ωB97XD/aug-cc-pVDZ level of theory. Polarisable continuum solvation model (PCM) is used for the implicit solvation. The adsorption free energy is evaluated to be -1.3kcal/mol at room temperature. It has been found that the adsorption free energy is more negative at low temperatures. Above 360K, the adsorption free energy is found to be positive.

Molecular simulations of the adsorption of aniline from waste-water

Molecular simulations of adsorption processes have received considerable attention. Despite the attention, exploration of the literature shows serious limitations, among which solvent and temperature effects are the most important. In this work, we propose a computational approach to study the adsorption of aniline (as an example of pollutant) from wastewater using coronene as adsorbent. We identified all possible adsorption sites using classical molecular dynamics for further optimization at the ωB97XD/aug-cc-pVDZ level of theory. Three different solvation schemes have been explored: implicit solvation of aniline + coronene, explicit solvation of aniline + coronene, and implicit-explicit solvation of aniline + coronene. For the explicit solvation, we used six water molecules, while the implicit solvation is performed using the PCM (polarizable continuum medium) solvation model. For each of the four cases (gas phase and the three solvation schemes), the adsorption free energy is evaluated as function of temperature from 200 K to 400 K. The results show that solvation has a considerable effects on the adsorption free energy. Furthermore, we noted that the adsorption free energy varies from -39.5 kJ mol^-1 at 200 K to 27.7 kJ mol^-1 at 400 K using the implicit-explicit solvation of aniline + coronene. This result highlights the importance of considering temperature effects in molecular simulations study of adsorption processes.

Conformational energies of microsolvated Na+ clusters with protic and aprotic solvents from GFNn-xTB methods

Performance of contemporary tight-binding semiempirical GFNn-xTB methods for the conformational energies of singly charged sodium clusters Na⁺ (S)_n (n = 4-8) with 3 protic and 8 aprotic solvents is examined against the reference RI-MP2/CBS method. The median Pearson correlation coefficients of ρ = 0.84 (GFN2-xTB) and ρ = 0.82 (GFN1-xTB) do not give the clear preference to any tested approach. GFN1-xTB method demonstrates more stable performance than its GFN2-xTB successor with the average mean absolute errors (MAEs)/mean signed errors (MSEs) of 1.2/0.2 and 2.3/1.6 kcal mol^-1 , respectively. Conformational energies produced by the computationally efficient DFT functional PBE and double-ζ basis set complemented with -D3(BJ) dispersion correction are suitable for the preliminary sampling (median ρ = 0.93), but should be used with a caution for the calculations of the average ensemble properties (MAE/MSE = 1.7/1.1 kcal mol^-1 ). Higher-ranking PBE0-D3(BJ) and ωB97M-V with triple-ζ basis sets yield significantly lower MAEs/MSEs of 0.55/0.20 and 0.51/0.23 kcal mol^-1 , respectively.

Dimethylsulfoxide (DMSO) Clusters Dataset: DFT Relative Energies, Non-Covalent Interactions, and Cartesian Coordinates

Theoretical understanding of dimethylsulfoxide (DMSO) liquid depends on the understanding of the DMSO clusters. In this work, we provide the structures and the energetics of the DMSO clusters. The structures have been generated using ABCluster and further optimized at the MP2/aug-cc-pVDZ level of theory. The final structures have been optimized at two different levels of theory: PW6B95D3/aug-cc-pVDZ and ωB97XD/aug-cc-pVDZ. The Cartesian coordinates of the structures optimized at the MP2/aug-cc-pVDZ level of theory are also reported. The relative energies of the structures can be used to locate the most favorable structures of the DMSO clusters. The Cartesian coordinates of the structures can be used for further investigations on DMSO clusters. In addition, we report the data related to the quantum theory of atoms in molecule (QTAIM) analysis of the investigated clusters. The QTAIM data reported in this work can be used to understand and determine the nature of non-covalent interactions in DMSO clusters. For further reading and discussion on the data reported here, please report to the original manuscript Malloum and Conradie (2022) [1].

Data to understand the nature of non-covalent interactions in the thiophene clusters

We have reported herein the data to understand the nature and number of non-covalent interactions that stabilize the structures of the thiophene clusters. In addition, we have also provided the optimized Cartesian coordinates of all the structures of the investigated thiophene clusters. Initially, the geometries have been generated using the ABCluster code which performs a global optimization to locate local and global minima structures of molecular clusters. The located geometries have been optimized at the MP2/aug-cc-pVDZ level of theory using Gaussian 16 suite of programs. To understand the nature of non-covalent interactions, we have performed a quantum theory of atoms in molecules (QTAIM) analysis on all the structures of the thiophene dimer. Furthermore, the QTAIM analysis has been performed also on the most stable structure of the thiophene trimer and tetramer. We have used the AIMAll program to perform the QTAIM analysis. The data reported in this paper contains the critical points, the bonds paths and their related properties, for each investigated structures. Besides, the data contains the optimized Cartesian coordinates of all the investigated structures of the thiophene clusters. This can be use for any further investigations involving thiophene clusters. For further information and analysis, the reader is referred to the original related research article (Malloum and Conradie, 2022).

QTAIM analysis dataset for non-covalent interactions in furan clusters

Furan clusters are very important to understand the dynamics and properties of the furan solvent. They can be used combined with quantum cluster equilibrium theory to theoretically determine the thermodynamics properties of the furan solvent. To understand the structures of the furan clusters, one needs to understand the non-covalent interactions that hold the furan molecules together. In this paper, we have provided the data necessary to understand the non-covalent interactions in furan clusters. Firstly, the structures of the furan clusters have been generated using classical molecular dynamics as implemented in the ABCluster code. Secondly, the generated structures have been fully optimized at the MP2/aug-cc-pVDZ level of theory. The optimized Cartesian coordinates of all the investigated structures are reported in this work to enable further investigations of the furan clusters. These Cartesian coordinates will save computational time for all further investigations involving the furan clusters. Thirdly, to understand the nature of the non-covalent interactions in furan clusters, we have performed a quantum theory of atoms in molecule (QTAIM) analysis using AIMAll program. Using QTAIM, we have provided the critical points, bond paths and their related properties for all the investigated structures. These data can be used to identify and classify the non-covalent interactions in furan clusters. The reader can refer to the original article for further information and discussion of the data provided herein Malloum and Conradie (2022) [1].

Structures, binding energies and non-covalent interactions of furan clusters

Understanding of the furan solvent is subjected to the knowledge of the structures of the furan clusters and interactions taking place therein. Although, furan clusters can be very important to determine the dynamics and the properties of the furan solvent, there has been only a few investigations reported on furan dimer. In this work, we have explored the potential energy surfaces (PESs) of the furan clusters using two incremental levels of theory. Structures have been initially generated using classical molecular dynamics followed by full optimization at the MP2/aug-cc-pVDZ level of theory. The results show that the most stable structure of the furan dimer has a stacking configuration while that of the trimer has a cyclic configuration. We have noted that the structures of the furan tetramer have no definite configurations. In addition, we have performed a quantum theory of atoms in molecule (QTAIM) analysis to identify all possible non-covalent interactions of the furan clusters. The results show that six different types of non-covalent interactions can be identified in furan clusters. We have noted that the CH⋯C and CH⋯O hydrogen bondings are the strongest non-covalent interactions while the H⋯H bonding interaction is found to be the weakest. Furthermore, we have assessed the performance of ten DFT functionals in calculating the binding energies of the furan clusters. The ten DFT functionals (M05, M05-2X, M06, M06-2X, M08HX, PBE0, ωB97XD, PW6B95D3, APFD and MN15) have been benchmarked to DLPNO-CCSD(T)/CBS. The functionals M05-2X and M06 are recommended for further affordable investigations of the furan clusters.