The rich and complex potential energy surface of the ethanol dimer

Exploration of Free Energy Surface of the Au10 Nanocluster at Finite Temperature

The first step in comprehending the properties of Au10 clusters is understanding the lowest energy structure at low and high temperatures. Functional materials operate at finite temperatures; however, energy computations employing density functional theory (DFT) methodology are typically carried out at zero temperature, leaving many properties unexplored. This study explored the potential and free energy surface of the neutral Au10 nanocluster at a finite temperature, employing a genetic algorithm coupled with DFT and nanothermodynamics. Furthermore, we computed the thermal population and infrared Boltzmann spectrum at a finite temperature and compared it with the validated experimental data. Moreover, we performed the chemical bonding analysis using the quantum theory of atoms in molecules (QTAIM) approach and the adaptive natural density partitioning method (AdNDP) to shed light on the bonding of Au atoms in the low-energy structures. In the calculations, we take into consideration the relativistic effects through the zero-order regular approximation (ZORA), the dispersion through Grimme's dispersion with Becke-Johnson damping (D3BJ), and we employed nanothermodynamics to consider temperature contributions. Small Au clusters prefer the planar shape, and the transition from 2D to 3D could take place at atomic clusters consisting of ten atoms, which could be affected by temperature, relativistic effects, and dispersion. We analyzed the energetic ordering of structures calculated using DFT with ZORA and single-point energy calculation employing the DLPNO-CCSD(T) methodology. Our findings indicate that the planar lowest energy structure computed with DFT is not the lowest energy structure computed at the DLPN0-CCSD(T) level of theory. The computed thermal population indicates that the 2D elongated hexagon configuration strongly dominates at a temperature range of 50-800 K. Based on the thermal population, at a temperature of 100 K, the computed IR Boltzmann spectrum agrees with the experimental IR spectrum. The chemical bonding analysis on the lowest energy structure indicates that the cluster bond is due only to the electrons of the 6 s orbital, and the Au d orbitals do not participate in the bonding of this system.

Acid Dissociation in (HX)n (H2 O)n Clusters (X=F, Cl, Br, I; n=2, 3)

In this work, we analyze the interactions between two or three hydrogen halide molecules and the same number of water moieties through a systematic exploration of their potential energy surfaces. Our results indicate that the most stable HF and HCl aggregates do not experience dissociation of any of the acid fragments, even with three water molecules. In contrast, in the HBr and HI clusters, one of the acid fragments does dissociate. While the global minimum of (HBr)₃ (H₂ O)₃ is a hydrogen-bridged bihalide anion (BrHBr^- ), which is persistent at temperatures up to 203 K, the lowest energy structure of (HI)₃ (H₂ O)₃ has a separated ion pair, but the motif with a bihalide anion (IHI^- ) is only 0.2 kcal mol^-1 above the global minimum. Among the more stable structures is a broad spectrum of contacts, including water⋯water, HX⋯water, and HX⋯HX hydrogen bonds, halogen bonds, ionic and long-range X⋯H contacts.

Relative Populations and IR Spectra of Cu38 Cluster at Finite Temperature Based on DFT and Statistical Thermodynamics Calculations

The relative populations of Cu₃₈ isomers depend to a great extent on the temperature. Density functional theory and nanothermodynamics can be combined to compute the geometrical optimization of isomers and their spectroscopic properties in an approximate manner. In this article, we investigate entropy-driven isomer distributions of Cu₃₈ clusters and the effect of temperature on their IR spectra. An extensive, systematic global search is performed on the potential and free energy surfaces of Cu₃₈ using a two-stage strategy to identify the lowest-energy structure and its low-energy neighbors. The effects of temperature on the populations and IR spectra are considered via Boltzmann factors. The computed IR spectrum of each isomer is multiplied by its corresponding Boltzmann weight at finite temperature. Then, they are summed together to produce a final temperature-dependent, Boltzmann-weighted spectrum. Our results show that the disordered structure dominates at high temperatures and the overall Boltzmann-weighted spectrum is composed of a mixture of spectra from several individual isomers.

Effects of Temperature on Enantiomerization Energy and Distribution of Isomers in the Chiral Cu13 Cluster

In this study, we report the lowest energy structure of bare Cu₁₃ nanoclusters as a pair of enantiomers at room temperature. Moreover, we compute the enantiomerization energy for the interconversion from minus to plus structures in the chiral putative global minimum for temperatures ranging from 20 to 1300 K. Additionally, employing nanothermodynamics, we compute the probabilities of occurrence for each particular isomer as a function of temperature. To achieve that, we explore the free energy surface of the Cu₁₃ cluster, employing a genetic algorithm coupled with density functional theory. Moreover, we discuss the energetic ordering of isomers computed with various density functionals. Based on the computed thermal population, our results show that the chiral putative global minimum strongly dominates at room temperature.

Theoretical Prediction of Structures, Vibrational Circular Dichroism, and Infrared Spectra of Chiral Be4B8 Cluster at Different Temperatures

Lowest-energy structures, the distribution of isomers, and their molecular properties depend significantly on geometry and temperature. Total energy computations using DFT methodology are typically carried out at a temperature of zero K; thereby, entropic contributions to the total energy are neglected, even though functional materials work at finite temperatures. In the present study, the probability of the occurrence of one particular Be4B8 isomer at temperature T is estimated by employing Gibbs free energy computed within the framework of quantum statistical mechanics and nanothermodynamics. To identify a list of all possible low-energy chiral and achiral structures, an exhaustive and efficient exploration of the potential/free energy surfaces is carried out using a multi-level multistep global genetic algorithm search coupled with DFT. In addition, we discuss the energetic ordering of structures computed at the DFT level against single-point energy calculations at the CCSD(T) level of theory. The total VCD/IR spectra as a function of temperature are computed using each isomer's probability of occurrence in a Boltzmann-weighted superposition of each isomer's spectrum. Additionally, we present chemical bonding analysis using the adaptive natural density partitioning method in the chiral putative global minimum. The transition state structures and the enantiomer-enantiomer and enantiomer-achiral activation energies as a function of temperature evidence that a change from an endergonic to an exergonic type of reaction occurs at a temperature of 739 K.

Exploration of Free Energy Surface and Thermal Effects on Relative Population and Infrared Spectrum of the Be6B11- Flux-Ional Cluster

The starting point to understanding cluster properties is the putative global minimum and all the nearby local energy minima; however, locating them is computationally expensive and difficult. The relative populations and spectroscopic properties that are a function of temperature can be approximately computed by employing statistical thermodynamics. Here, we investigate entropy-driven isomers distribution on Be₆B₁₁⁻ clusters and the effect of temperature on their infrared spectroscopy and relative populations. We identify the vibration modes possessed by the cluster that significantly contribute to the zero-point energy. A couple of steps are considered for computing the temperature-dependent relative population: First, using a genetic algorithm coupled to density functional theory, we performed an extensive and systematic exploration of the potential/free energy surface of Be₆B₁₁⁻ clusters to locate the putative global minimum and elucidate the low-energy structures. Second, the relative populations’ temperature effects are determined by considering the thermodynamic properties and Boltzmann factors. The temperature-dependent relative populations show that the entropies and temperature are essential for determining the global minimum. We compute the temperature-dependent total infrared spectra employing the Boltzmann factor weighted sums of each isomer’s infrared spectrum and find that at finite temperature, the total infrared spectrum is composed of an admixture of infrared spectra that corresponds to the spectra of the lowest-energy structure and its isomers located at higher energies. The methodology and results describe the thermal effects in the relative population and the infrared spectra.

Activity coefficients of binary methanol alcohol mixtures from cluster weighting

The hydrogen bond network of different small alcohols is investigated via cluster analysis. Methanol/alcohol mixtures are studied with increasing chain length and branching of the molecule. Those changes can play an important role in different fields, including solvent and metal extraction. The extended tight binding method GFN2-xTB allows the evaluation and geometry optimization of thousands of clusters built via a genetic algorithm. Interaction energies and geometries are evaluated and discussed for the neat systems. Thermodynamic properties, such as vaporization enthalpies and activity coefficients, are calculated with the binary quantum cluster equilibrium (bQCE) approach using our in-house code peacemaker 2.8. Combined distribution functions of the distances against the angles of the hydrogen bonds are evaluated for neat and mixed clusters and weighted by the equilibrium populations achieved from bQCE calculations.

Growth Pattern, Stability, and Properties of Complexes of C2H5OH and nCO2 (n = 1-5) Molecules: A Theoretical Study

This work is dedicated to theoretically investigate the formation process of C2H5OH···nCO2 (n = 1-5) complexes and to shed light on the nature of interactions formed under the variation of CO2 concentration. It is found that CO2 molecules tend to locate around the polarized -OH group to interact with the lone pairs of the O atom. The interaction of ethanol with three CO2 molecules (C2H5OH···3CO2) induces the most stable structure in the sequence considered. The atoms in molecules (AIM), NCIplot, and natural bond orbital (NBO) analyses point out that the Oethanol···CCO2 tetrel bond overcomes hydrogen, chalcogen, and CO2···CO2 tetrel-bonded interactions and mainly contributes to the strength of C2H5OH···nCO2 (n = 1-5) complexes. All intermolecular interactions in the examined complexes are weakly noncovalent, and their positive cooperativity is evaluated to be slightly weaker than that of CO2 pure systems. SAPT2+ and molecular electrostatic potential (MEP) calculations indicate that the electrostatic force is the main factor underlying the attractive interplay in the complexes of C2H5OH and CO2.

Exploration of the potential energy surfaces of small ethanol clusters

The potential energy surfaces of small ethanol clusters, from dimer to pentamer, have been thoroughly explored using two different levels of theory. There is a clear relative energy gap between cyclic, linear and branched cyclic structures.

Fe3+chelating quinoline–hydrazone hybrids with proven cytotoxicity, leishmanicidal, and trypanocidal activities

Neglected tropical diseases cause great concern in developing countries where there are millions of reported infected humans. Our calculations support a direct relationship between biological activity and the Fe³⁺chelating ability of the shown set of quinoline–hydrazone hybrids.

Does H4SO5 exist?

The possible existence of H₄SO₅ in aqueous sulfuric acid is analyzed in detail. For bare H₄SO₅, the computed free energy barrier for the exergonic transformation of H₄SO₅ into the H₂SO₄H₂O complex is only 3.8 kcal mol^-1. The presence of water or sulfuric acid catalyzes the dehydration to such an extent that it becomes almost a barrierless process. In the gas phase, dehydration of H₄SO₅ is an autocatalytic reaction as the water molecule produced by the decomposition of one H₄SO₅ molecule induces further dissociation. Thus, in solution, the surrounding water molecules make the para-sulfuric acid a very vulnerable species to exist. The simulated Raman spectra also corroborate the absence of H₄SO₅ in solution.

Microsolvation of NO3−: structural exploration and bonding analysis

A rich and complex structural diversity is uncovered in the microsolvation of the nitrate anion.

How strong are the metallocene-metallocene interactions? Cases of ferrocene, ruthenocene, and osmocene

An exhaustive exploration of the potential energy surfaces of ferrocene, ruthenocene and osmocene dimers has been performed. Our computations involving dispersion show that only four different isomers are present in each metallocene dimer. The collective action of small interaction energies of dispersive nature leads to a dissociation energy of 7.5 kcal mol(-1) for the ferrocene dimer. Dispersion has strong effects on the geometrical parameters, reducing the M···M distances by almost 1 Å. Our results also reveal that inclusion of entropic factors modifies the relative stability of the complexes. The nature of bonding is examined using the energy decomposition analysis and the non-covalent interaction index. Both analyses indicate that dispersion is the major contributing factor in stabilizing a metallocene dimer.

Theoretical tools to distinguish O-ylides from O-ylidic complexes in carbene–solvent interactions

Orbital interactions and bond indices are among the theoretical tools suitable to distinguishO-ylides fromO-ylidic carbene–solvent complexes.