Structures and spectroscopy of medium size protonated ammonia clusters at different temperatures, H+(NH3)10–16

Clusters of solvated ferrous ion in water-ammonia mixture: Structures and noncovalent interactions

The behavior of metal ions is commonly studied in pure solvent although, in our daily life, these metals are involved in mixtures of solvents. In the present study, we investigated structures, relative stabilities and temperature dependance of solvated ferrous ion in water-ammonia mixture solvent at 0K and at various temperatures ranging from 25K to 400K. All the calculations are performed at the MN15 level of theory associated with the aug-cc-pVDZ basis set. For deep understanding of binding patterns in solvated ferrous ion in water-ammonia mixture solvent, noncovalent interactions are presented based on the QTAIM analysis using AIMAll. Our results prove that the ferrous ion is more stable when it is solvated by ammonia instead of water. In addition, hydrogen bonds are weakened by the presence of ammonia molecules. The temperature dependence of the different obtained geometries indicates that from s=6 (s is the sum of water and ammonia molecules around the ferrous ion), when the number of water molecules is almost equal to that of ammonia, the structures with coordination number 5 are dominant. However, the coordination number is six when there are a maximum water molecules (rich water solution) or maximum ammonia molecules (rich ammonia solution) around the ferrous ion (for s≥6). The QTAIM analysis shows that there are two coordination bondings and four hydrogen bondings. Furthermore, it is found that the Fe2+⋯N coordination bondings are stronger than the Fe2+⋯O confirming that the ferrous ion prefers to be solvated by ammonia instead of water.

Hydrogen bond network structures of protonated dimethylamine clusters H+(DMA)n (n = 3-7)

Infrared spectroscopy of protonated dimethylamine clusters, H+(DMA)n, (n = 3-7), and their Ar-tagged clusters was performed in the NH and CH stretching vibrational region to explore their hydrogen bond network structures. A stable isomer search and vibrational spectral simulations of the clusters were also carried out to support the interpretations of the observed spectra. Weakly hydrogen-bonded NH stretching vibrational bands, which are characteristic of cyclic structures of small-sized protonated clusters, are observed in the spectra of the Ar-tagged clusters of n ≥ 5, while only linear chain type structures are suggested for the Ar-tagged clusters of n = 3-4 and the bare clusters of all the sizes. These results demonstrate that the size and temperature dependence of the hydrogen bond network structures of the protonated dimethylamine clusters is analogous to that of protonated monohydric alcohol clusters.

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.

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.

Mg(II) and Ca(II) Microsolvation by Ammonia: Born-Oppenheimer Molecular Dynamics Studies

We report the structural and energetic features of the Mg²⁺ and Ca²⁺ cations in ammonia microsolvation environments. Born-Oppenhemier molecular dynamics studies are carried out for [Mg(NH₃)_n]²⁺ and [Ca(NH₃)_n]²⁺ clusters with n = 2, 3, 4, 6, 8, 20, and 27 at 300 K based on hybrid density functional theory calculations. We determine binding energies per ammonia molecule and the metal cation solvation patterns as a function of the number of molecules. The general trend for Mg²⁺ is that the Mg-N distances increase as a function of n until the first solvation shell is populated by six ammonia molecules, and then the distances slightly decrease while CN = 6 does not change. For Ca²⁺, the first solvation shell at room temperature is populated by eight ammonia molecules for clusters with more than one solvation shell, leading to a different structure from that of [Ca(NH₃)₆]²⁺ hexamine. The evaporation of NH₃ molecules was found at 300 K only for Mg²⁺ clusters with n ≥ 10; this was not the case for Ca²⁺ clusters. Vibrational spectra are obtained for all of the clusters, and the evolution of the main features is discussed. EXAFS spectra are also presented for the [Mg(NH₃)₂₇(NH₃)₂₇]²⁺ and [Ca(NH₃)₂₇]²⁺ clusters, which yield valuable data to be compared with experimental data in the liquid phase, as previously done for the aqueous solvation of these dications.

Structures and relative stability of hydrated ferrous ion clusters and temperature effects

Structures of solvated ferrous ion clusters have been investigated in the singlet and quintet spin states of the ferrous ion. Relative stabilities of isomers are also discussed at different temperatures.

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.

Intramolecular Reactions in Ionized Ammonia Clusters: A Direct Ab Initio Molecular Dynamics Study

Ammonia cluster cations are a chemical species that has recently attracted considerable research attention as an ion-molecule reaction species in the planetary atmosphere, surface reaction species in materials chemistry, and super-alkali species. Reactions of the radical cation of an ammonia cluster, [(NH₃)_n]⁺ (n = 2-6), following the ionization of the parent neutral cluster, were investigated using direct ab initio molecular dynamics to elucidate the reactions of the ammonia cluster cation under astrochemical conditions. The calculations showed that two competing reaction channels-proton transfer (PT) channel and complex formation channel-operate after the ionization of neutral clusters. In the PT channel, a proton of NH₃⁺ was transferred to a neighboring ammonia molecule. The PT channel was found in all clusters (n = 2-6). Reaction via the PT channel became faster with increasing cluster size and saturated around n = 5-6. In the complex formation channel, a face-to-face complex having a H₃N-NH₃⁺ structure (with a N-N bond) was formed. This channel was found only in larger clusters (n = 5-6). Time scales of PT and complex formation channels were calculated to be 20-30 and 40-50 fs, respectively. The reaction mechanism was discussed based on the results of theoretical calculations.

Structures of the solvated copper(ii) ion in ammonia at various temperatures

We investigated theoretically the structures and relative stabilities of the solvated copper(ii) ion in ammonia, Cu²⁺(NH₃)_n, n = 1–10.

Global and local minima of protonated acetonitrile clusters

Potential energy surfaces of protonated acetonitrile clusters have been explored to locate global and local minima energy structures. The structures are stabilized by strong hydrogen bonds, anti-parallel dimers, dipole–dipole and CH⋯N interactions.

Structures, relative stability and binding energies of neutral water clusters, (H2O)2–30

We have revised the structures of neutral water clusters, (H₂O)_n=2–30, with the affordable M06-2X functional, presenting up to 25 isomers for each cluster size.

Structures of solvated ferrous ion clusters in ammonia and spin-crossover at various temperatures

The solvated ferrous ion in ammonia is hexa-coordinated, irrespective of the temperature.

Ionization of Ammonia Nanoices with Adsorbed Methanol Molecules

Large ammonia clusters represent a model system of ices that are omnipresent throughout the space. The interaction of ammonia ices with other hydrogen-boding molecules such as methanol or water and their behavior upon an ionization are thus relevant in the astrochemical context. In this study, ammonia clusters (NH₃) _N with the mean size N̅ ≈ 230 were prepared in molecular beams and passed through a pickup cell in which methanol molecules were adsorbed. At the highest exploited pickup pressures, the average composition of (NH₃) _N(CH₃OH) _M clusters was estimated to be N: M ≈ 210:10. On the other hand, the electron ionization of these clusters yielded about 75% of methanol-containing fragments (NH₃) _n(CH₃OH) _mH⁺ compared to 25% contribution of pure ammonia (NH₃) _nH⁺ ions. On the basis of this substantial disproportion, we propose the following ionization mechanism: The prevailing ammonia is ionized in most cases, resulting in NH₄⁺ core solvated most likely with four ammonia molecules, yielding the well-known "magic number" structure (NH₃)₄NH₄⁺. The methanol molecules exhibit a strong propensity for sticking to the fragment ion. We have also considered mechanisms of intracluster reactions. In most cases, proton transfer between ammonia units take place. The theoretical calculations suggested the proton transfer either from the methyl group or from the hydroxyl group of the ionized methanol molecule to ammonia to be the energetically open channels. However, the experiments with selectively deuterated methanols did not show any evidence for the D⁺ transfer from the CD₃ group. The proton transfer from the hydroxyl group could not be excluded entirely or confirmed unambiguously by the experiment.

Spectroscopy of Proton Coordination with Ethylenediamine

Protonated ethylenediamine monomer, dimer, and trimer were produced in the gas phase by an electrical discharge/supersonic expansion of argon seeded with ethylenediamine (C₂H₈N₂, en) vapor. Infrared spectra of these ions were measured in the region from 1000 to 4000 cm^-1 using laser photodissociation and argon tagging. Computations at the CBS-QB3 level were performed to explore possible isomers and understand the infrared spectra. The protonated monomer exhibits a gauche conformation and an intramolecular hydrogen bond. Its parallel shared proton vibration occurs as a broad band around 2785 cm^-1, despite the formally equivalent proton affinities of the two amino groups involved, which usually leads to low frequency bands. The barrier to intramolecular proton transfer is 2.2 kcal mol^-1 and does not vanish upon addition of the zero-point energy, unlike the related protonated ammonia dimer. The structure of the dimer is formed by chelation of the monomer's NH₃⁺ group, thereby localizing the excess proton and increasing the frequency of the intramolecular shared proton vibration to 3157 cm^-1. Other highly fluxional dimer structures with facile intermolecular proton transfer and concomitant structural reorganization were computed to lie within 2 kcal mol^-1 of the experimentally observed structure. The spectrum of the trimer is rather diffuse, and a clear assignment is not possible. However, an isomer with an intramolecular proton transfer like that of the monomer is most consistent with the experimental spectrum.

Solvation energies of the proton in methanol revisited and temperature effects

Various functionals assessing solvation free energies and enthalpies of the proton in methanol.

Quantitative Temperature Dependence of the Microscopic Hydration Structures Investigated by Ultraviolet Photodissociation Spectroscopy of Hydrated Phenol Cations

To discuss the temperature effect on microscopic hydration structures in clusters, relative populations of the isomers having different hydration structures at well-defined temperatures are quite important. In the present study, we measured ultraviolet photodissociation spectra of the temperature-controlled hydrated phenol cation [PhOH(H₂O)₅]⁺ trapped in the 22-pole ion trap. Two isomers having a distinct hydration motif with each other are identified in the spectra, and a clear change in the relative populations is observed in the temperature range from 30 to 150 K. This behavior is quantitatively interpreted by statistical mechanical estimation based on density functional theory calculations. A ring with tail-type hydration motif is dominant in cold conditions, whereas a chain-like motif is dominant in hot conditions. The present study provides very quantitative information about the temperature effect on the microscopic hydration structures.