Search

Molecular dynamics simulations have been performed to investigate the solvation characteristics of neutral fullerene (C(60)) and charged fulleride anion (C(60)(5-)) in liquid ammonia. Potassium ions are present as counterions in the system containing fulleride ion. In addition to solvation characteristics, dynamical properties of solvation shells are also found out for both the neutral and anionic solutes. Our results reveal the presence of a rather large solvation shell of ammonia molecules around the C(60)(5-) ion. It is found that the ammonia molecules are more closely packed in the first solvation shell of C(60)(5-) than that of C(60). The distributions of ammonia molecules in the solvation shells of C(60) and C(60)(5-) solutes together with hydrogen bonding characteristics of the solvent in different solvation shells are investigated. It is found that the solvation of the small counterions (K(+)) in liquid ammonia is affected very little by the presence of the large C(60)(5-) anion. Regarding the dynamics of ammonia in solvation shells, it is found that the residence, translational and rotational dynamics of ammonia molecules differ significantly between the solvation shells of the neutral and charged fullerene solutes, especially in the first solvation shells. The average lifetimes of ammonia-ammonia hydrogen bonds are calculated from both continuous and intermittent hydrogen bond correlation functions. The calculations of binding energies reveal that the hydrogen bonds are weaker, hence short lived in the solvation shell of C(60)(5-) compared to those in the solvation shell of neutral C(60) and also in bulk liquid ammonia.

We have investigated the solvation structure and migration kinetics of an excess proton and a hydroxide ion in water-ammonia mixed liquids of varying composition by means of ab initio molecular dynamics simulations. The excess proton is always found to be attached to an ammonia molecule to form the ammonium ion. Migration of the excess proton is found to occur very occasionally from one ammonia to the other but no proton transfer to a water molecule is observed during the entire simulations. Also, when the ammonium ion is solvated in water only, its hydrogen bond dynamics and rotation are found to occur at a faster rate than those in water-ammonia mixtures. For water-ammonia mixtures containing a proton less, the defect is found to stay like the hydroxide ion. For these systems, occasional proton transfer is found to occur only through the hydrogen bonded chains of water molecules in these water-ammonia mixtures. No proton transfer is found to take place from an ammonia molecule. The presence of ammonia molecules makes the realization of proper presolvated state of the hydroxide ion to accept a proton a more difficult process and, as a result, the rate of proton transfer and migration kinetics of the hydroxide ion in water-ammonia mixtures are found to be slower than that in liquid water and these rates are found to slow down further with increase of ammonia concentration., (© 2012 American Institute of Physics)

First principles molecular dynamics simulations are carried out to investigate the solvation of an excess electron and a lithium atom in mixed water-ammonia cluster (H(2)O)(5)NH(3) at a finite temperature of 150 K. Both [(H(2)O)(5)NH(3)](-) and Li(H(2)O)(5)NH(3) clusters are seen to display substantial hydrogen bond dynamics due to thermal motion leading to many different isomeric structures. Also, the structures of these two clusters are found to be very different from each other and also very different from the corresponding neutral cluster without any excess electron or the metal atom. Spontaneous ionization of Li atom occurs in the case of Li(H(2)O)(5)NH(3). The spatial distribution of the singly occupied molecular orbital shows where and how the excess (or free) electron is primarily localized in these clusters. The populations of single acceptor (A), double acceptor (AA), and free (NIL) type water and ammonia molecules are found to be significantly high. The dangling hydrogens of these type of water or ammonia molecules are found to primarily capture the free electron. It is also found that the free electron binding motifs evolve with time due to thermal fluctuations and the vertical detachment energy of [(H(2)O)(5)NH(3)](-) and vertical ionization energy of Li(H(2)O)(5)NH(3) also change with time along the simulation trajectories. Assignments of the observed peaks in the vibrational power spectra are done and we found a one to one correlation between the time-averaged populations of water and ammonia molecules at different H-bonding sites with the various peaks of power spectra. The frequency-time correlation functions of OH stretch vibrational frequencies of these clusters are also calculated and their decay profiles are analyzed in terms of the dynamics of hydrogen bonded and dangling OH modes. It is found that the hydrogen bond lifetimes in these clusters are almost five to six times longer than that of pure liquid water at room temperature.

The preferential solvation of solutes in mixed solvent systems is an interesting phenomenon that plays important roles in solubility and kinetics. In the present study, solvation of a lithium atom in aqueous ammonia solution has been investigated from first principles molecular dynamics simulations. Solvation of alkali metal atoms, like lithium, in aqueous and ammonia media is particularly interesting because the alkali metal atoms release their valence electrons in these media so as to produce solvated electrons and metal counterions. In the present work, first principles simulations are performed employing the Car-Parrinello molecular dynamics method. Spontaneous ionization of the Li atom is found to occur in the mixed solvent system. From the radial distribution functions, it is found that the Li(+) ion is preferentially solvated by water and the coordination number is mostly four in its first solvation shell and exchange of water molecules between the first and second solvation shells is essentially negligible in the time scale of our simulations. The Li(+) ion and the unbound electron are well separated and screened by the polar solvent molecules. Also the unbound electron is primarily captured by the hydrogens of water molecules. The diffusion rates of Li(+) ion and water molecules in its first solvation shell are found to be rather slow. In the bulk phase, the diffusion of water is found to be slower than that of ammonia molecules because of strong ammonia-water hydrogen bonds that participate in solvating ammonia molecules in the mixture. The ratio of first and second rank orientational correlation functions deviate from 3, which suggests a deviation from the ideal Debye-type orientational diffusion. It is found that the hydrogen bond lifetimes of ammonia-ammonia pairs is very short. However, ammonia-water H-bonds are found to be quite strong when ammonia acts as an acceptor and these hydrogen bonds are found to live longer than even water-water hydrogen bonds.

First principles molecular dynamics simulations are carried out to investigate the solvation of an excess electron and a lithium atom in mixed water-ammonia cluster (H2O)5NH3 at a finite temperature of 150 K. Both [(H2O)5NH3]- and Li(H2O)5NH3 clusters are seen to display substantial hydrogen bond dynamics due to thermal motion leading to many different isomeric structures. Also, the structures of these two clusters are found to be very different from each other and also very different from the corresponding neutral cluster without any excess electron or the metal atom. Spontaneous ionization of Li atom occurs in the case of Li(H2O)5NH3. The spatial distribution of the singly occupied molecular orbital shows where and how the excess (or free) electron is primarily localized in these clusters. The populations of single acceptor (A), double acceptor (AA), and free (NIL) type water and ammonia molecules are found to be significantly high. The dangling hydrogens of these type of water or ammonia molecules are found to primarily capture the free electron. It is also found that the free electron binding motifs evolve with time due to thermal fluctuations and the vertical detachment energy of [(H2O)5NH3]- and vertical ionization energy of Li(H2O)5NH3 also change with time along the simulation trajectories. Assignments of the observed peaks in the vibrational power spectra are done and we found a one to one correlation between the time-averaged populations of water and ammonia molecules at different H-bonding sites with the various peaks of power spectra. The frequency-time correlation functions of OH stretch vibrational frequencies of these clusters are also calculated and their decay profiles are analyzed in terms of the dynamics of hydrogen bonded and dangling OH modes. It is found that the hydrogen bond lifetimes in these clusters are almost five to six times longer than that of pure liquid water at room temperature. [ABSTRACT FROM AUTHOR]

The structure and energetics of water-ammonia mixed clusters with an excess electron, [(H2O)n(NH3)m]- with m=1, n=2–6 and m=2, n=2, and also the corresponding neutral clusters are investigated in detail by means of ab initio quantum chemical calculations. The authors focus on the localization structure of the excess electron with respect to its surface versus interiorlike states, its binding to ammonia versus water molecules, the spatial and orientational arrangement of solvent molecules around the excess electron, the changes of the overall hydrogen-bonded structure of the clusters as compared to those of the neutral ones and associated dipole moment changes, vertical detachment energies of the anionic clusters, and also the vertical attachment energies of the neutral clusters. It is found that the hydrogen-bonded structure of the anionic clusters are very different from those of the neutral clusters unlike the case of water-ammonia dimer anion, and these changes in structural arrangements lead to drastically different dipole moments of the anionic and the neutral clusters. The spatial distribution of the singly occupied molecular orbital holding the excess electron shows only surface states for the smaller clusters. However, for n=5 and 6, both surface and interiorlike binding states are found to exist for the excess electron. For the surface states, the excess electron can be bound to the dangling hydrogens of either an ammonia or a water molecule with different degrees of stability and vertical detachment energies. The interiorlike states, wherever they exist, are found to have a higher vertical detachment energy than any of the surface states of the same cluster. Also, for interiorlike states, the ammonia molecule with its dangling hydrogens is always found to stay on top or on a far side of the charge density of the excess electron without participating in the hydrogen bond network of the cluster; the intermolecular hydrogen bonds are formed by the water molecules only which add to the overall stability of these anionic clusters. [ABSTRACT FROM AUTHOR]

The ammonia dimer (NH[sub 3])[sub 2] has been investigated using high-level ab initio quantum chemistry methods and density functional theory. The structure and energetics of important isomers are obtained to unprecedented accuracy without resorting to experiment. The global minimum of eclipsed C[sub s] symmetry is characterized by a significantly bent hydrogen bond which deviates from linearity by as much as ≈20°. In addition, the so-called cyclic C[sub 2h] structure, resulting from further bending which leads to two equivalent “hydrogen bonding contacts,” is extremely close in energy on an overall flat potential energy surface. It is demonstrated that none of the currently available [generalized gradient approximation (GGA), meta-GGA, and hybrid] density functionals satisfactorily describe the structure and relative energies of this nonlinear hydrogen bond. We present a novel density functional, HCTH/407+, which is designed to describe this sort of hydrogen bond quantitatively on the level of the dimer, contrary to, e.g., the widely used BLYP functional. This improved generalized gradient approximation functional is employed in Car–Parrinello ab initio molecular dynamics simulations of liquid ammonia to judge its performance in describing the associated liquid. Both the HCTH407+ and BLYP functionals describe the properties of the liquid well as judged by analysis of radial distribution functions, hydrogen bonding structure and dynamics, translational diffusion, and orientational relaxation processes. It is demonstrated that the solvation shell of the ammonia molecule in the liquid phase is dominated by steric packing effects and not so much by directional hydrogen bonding interactions. In addition, the propensity of ammonia molecules to form bifurcated and multifurcated hydrogen bonds in the liquid phase is found to be negligibly small. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We have carried out ab initio molecular dynamics simulations of a liquid-vapor interfacial system consisting of a mixture of water and ammonia molecules. We have made a detailed analysis of the structural and dynamical properties of the bulk and interfacial regions of the mixture. Among structural properties, we have looked at the inhomogeneous density profiles of water and ammonia molecules, hydrogen bond distributions, orientational profiles, and also vibrational frequency distributions of bulk and interfacial molecules. It is found that the interfacial molecules show preference for specific orientations so as to form water-ammonia hydrogen bonds at the interface with ammonia as the acceptor. The structure of the system is also investigated in terms of inter-atomic voids present in the system. Among the dynamical properties, we have calculated the diffusion, orientational relaxation, hydrogen bond dynamics, and vibrational spectral diffusion in bulk and interfacial regions. It is found that the diffusion and orientation relaxation of the interfacial molecules are faster than those of the bulk. However, the hydrogen bond lifetimes are longer at the interface which can be correlated with the time scales found from the decay of frequency time correlations. [ABSTRACT FROM AUTHOR]