Understanding the molecular details of CO(2)-sorbent interactions is critical for the design of better carbon-capture systems. Here we report crystallographic resolution of CO(2) molecules and their binding domains in a metal-organic framework functionalized with amine groups. Accompanying computational studies that modeled the gas sorption isotherms, high heat of adsorption, and CO(2) lattice positions showed high agreement on all three fronts. The modeling apportioned specific binding interactions for each CO(2) molecule, including substantial cooperative binding effects among the guest molecules. The validation of the capacity of such simulations to accurately model molecular-scale binding bodes well for the theory-aided development of amine-based CO(2) sorbents. The analysis shows that the combination of appropriate pore size, strongly interacting amine functional groups, and the cooperative binding of CO(2) guest molecules is responsible for the low-pressure binding and large uptake of CO(2) in this sorbent material.

Molecular-dynamics simulations are used to study the stability of structure II hydrogen clathrates with different H2 guest occupancies. Simulations are done at pressures of 2.5 kbars and 1.013 bars and for temperatures ranging from 100 to 250 K. For a structure II unit cell with 136 water molecules, H2 guest molecule occupancies of 0-64 are studied with uniform occupancies among each type of cage. The simulations show that at 100 K and 2.5 kbars, the most stable configurations have single occupancy in the small cages and quadruple occupancy in the large cages. The optimum occupancy for the large cages decreases as the temperature is raised. Double occupancy in the small cages increases the energy of the structures and causes tetragonal distortion in the unit cell. The spatial distribution of the hydrogen guest molecules in the cages is determined by studying the guest-water and guest-guest radial distribution functions at various temperatures.

Position and orientation of water protons need to be specified when the molecular simulation studies are performed for clathrate hydrates. Positions of oxygen atoms in water are experimentally determined by X-ray diffraction analysis of clathrate hydrate structures, but positions of water hydrogen atoms in the lattice are disordered. This study reports a determination of the water proton coordinates in unit cell of structure I (sI), II (sII), and H (sH) clathrate hydrates that satisfy the ice rules, have the lowest potential energy configuration for the protons, and give a net zero dipole moment. Possible proton coordinates in the unit cell were chosen by analyzing the symmetry of protons on the hexagonal or pentagonal faces in the hydrate cages and generating all possible proton distributions which satisfy the ice rules. We found that in the sI and sII unit cells, proton distributions with small net dipole moments have fairly narrow potential energy spreads of about 1 kJ∕mol. The total Coulomb potential on a test unit charge placed in the cage center for the minimum energy∕minimum dipole unit cell configurations was calculated. In the sI small cages, the Coulomb potential energy spread in each class of cage is less than 0.1 kJ∕mol, while the potential energy spread increases to values up to 6 kJ∕mol in sH and 15 kJ∕mol in the sII cages. The guest environments inside the cages can therefore be substantially different in the sII case. Cartesian coordinates for oxygen and hydrogen atoms in the sI, sII, and sH unit cells are reported for reference.

Molecular dynamics simulations are used to compare microscopic structures and guest dynamics to macroscopic properties in structure II clathrate hydrates with cyclopentane, tetrahydrofuran (THF), 1,3-dioxolane, tetrahydropyran (THP), and p-dioxane as guests. Significant differences are observed between structural parameters and rotational dynamics for the different guests. The simulations show the formation of guest-host hydrogen bonds between the ether oxygen atoms of THF and THP and the cage water hydrogen atoms of the clathrate but the absence of similar hydrogen bonds in the clathrate hydrates of the other guests on the time scale of the calculations. This guest-host hydrogen bonding leads to the formation of Bjerrum L-defects in the clathrate water lattice where two adjacent water molecules have no covalently bonded hydrogen atom between them. Unlike Bjerrum defects of ice lattices, these guest-induced L-defects are not accompanied by the formation of a D-defect at an adjacent site in the water lattice. At the simulation temperature of 200 K, the guest-water hydrogen bonds in the THF clathrate are short lived (lifetime less than 1 ps) but in the THP they are longer lived (a minimum of 100 ps). A van't Hoff plot for the probability of defect formation in THF as a function of temperature gives an activation barrier of approximately 8.3 kJ/mol for guest-host defect formation in the THF clathrate. The consequences of the defect formation on the thermal expansivity, isothermal compressibility, dipole-dipole correlation function, and mechanical stability of the clathrate are discussed.

Molecular dynamics simulations are used to study the dynamics and transport properties of 12 room-temperature ionic liquids of the 1-alkyl-3-methylimidazolium [amim](+) (alkyl = methyl, ethyl, propyl, and butyl) family with PF(6)(-), NO(3)(-), and Cl(-) counterions. The explicit atom transferable force field of Canongia Lopes et al. [J. Phys. Chem. B 108, 2038 (2004)] is used in the simulations. In this first part, the dynamics of the ionic liquids are characterized by studying the mean-square displacement (MSD) and the velocity autocorrelation function (VACF) for the centers of mass of the ions at 400 K. Trajectory averaging was employed to evaluate the diffusion coefficients at two temperatures from the linear slope of MSD(t) functions in the range of 150-300 ps and from the integration of the VACF(t) functions at 400 K. Detailed comparisons are made between the diffusion results from the MSD and VACF methods. The diffusion coefficients from the integration of the VACFs are closer to experimental values than the diffusion coefficients calculated from the slope of MSDs. Both methods can show good agreement with experiment in predicting relative trends in the diffusion coefficients and determining the role of the cation and anion structures on the dynamical behavior of this family of ionic liquids. The MSD and self-diffusion of relatively heavier imidazolium cations are larger than those of the lighter anions from the Einstein results, except for the case of [bmim][Cl]. The cationic transference number generally decreases with temperature, in good agreement with experiments. For the same anion, the cationic transference numbers decrease with increasing length of the alkyl chain, and for the same cation, the trends in the cationic transference numbers are [NO(3)](-) < [Cl](-) < [PF(6)](-). The trends in the diffusion coefficient in the series of cations with identical anions are [emim](+) > [pmim](+) > [bmim](+) and those for anions with identical cations are [NO(3)](-) > [PF(6)](-) > [Cl](-). The [dmim](+) has a relatively low diffusion coefficient due to its symmetric structure and good packing in the liquid phase. The major factor for determining the magnitude of the self-diffusion is the geometric shape of the anion of the ionic liquid. Other important factors are the ion size and the charge delocalization in the anion.

There is interest in the role of ammonia on Saturn's moons Titan and Enceladus as the presence of water, methane, and ammonia under temperature and pressure conditions of the surface and interior make these moons rich environments for the study of phases formed by these materials. Ammonia is known to form solid hemi-, mono-, and dihydrate crystal phases under conditions consistent with the surface of Titan and Enceladus, but has also been assigned a role as water-ice antifreeze and methane hydrate inhibitor which is thought to contribute to the outgassing of methane clathrate hydrates into these moons' atmospheres. Here we show, through direct synthesis from solution and vapor deposition experiments under conditions consistent with extraterrestrial planetary atmospheres, that ammonia forms clathrate hydrates and participates synergistically in clathrate hydrate formation in the presence of methane gas at low temperatures. The binary structure II tetrahydrofuran + ammonia, structure I ammonia, and binary structure I ammonia + methane clathrate hydrate phases synthesized have been characterized by X-ray diffraction, molecular dynamics simulation, and Raman spectroscopy methods.ice | single crystal X-ray diffraction | hydrogen bonding | hydrate inhibitors | ethane A mmonia has long been seen as a key species in extraterrestrial space, both interstellar and on outer planets, moons, and comets and the interplay of ammonia, methane, and water has been the subject of a considerable number of studies and speculation (1-8). The main role assigned to ammonia has been that of an antifreeze for ice and clathrate hydrate formation, modifying the stability region of the solid ice and methane clathrate hydrate phases as a thermodynamic inhibitor (2, 5, 9). However, ammonia is a methane-sized molecule, thus based on size alone it has the potential for being a suitable guest for clathrate hydrate cages. Issues that may have prevented ammonia from being considered as a suitable clathrate guest molecule include the notion that guest species need to be hydrophobic in order to be incorporated into clathrates, and the observation of a number of stoichiometric nonclathrate phases of ammonia and water obtained upon cooling aqueous ammonia solutions (2). Previous experimental work on the water-ammonia and water-methaneammonia systems had not shown evidence for the enclathration of ammonia (5,(10)(11)(12)(13)(14)(15)(16)(17)(18)(19). Close inspection, however, shows that the low pressure ammonia dihydrate (18) and the high pressure phase II of ammonia monohydrate (19) have structural features in common with canonical clathrate and semiclathrate structures.Recent structural analysis and molecular simulations have shown that some guest molecules which form strong hydrogen bonds with the water framework of the clathrate hydrate lattice may nonetheless produce stable phases (20-23). It is therefore reasonable to consider that ammonia has potential as a clathrate guest molecule. Furthermore, previous experimental and computational studies hav...

A scanning-tunneling microscope has been used to induce efficient local desorption of benzene from Si(100) at low currents (<100 pA), sample biases (approximately -2.4 V) and temperatures (22 K). A theoretical model based upon first principles electronic structure calculations and quantum mechanical wave packet dynamics describes this process as occurring via transient ionization of a pi state of the adsorbed molecule. This model accounts for the unexpected efficiency and sharp threshold of the yield.

Molecular dynamics simulations are performed to determine the melting points of aluminum nanoparticles of 55-1000 atoms with the Streitz-Mintmire [Phys. Rev. B 1994, 50, 11996] variable-charge electrostatic plus potential. The melting of the nanoparticles is characterized by studying the temperature dependence of the potential energy and Lindemann index. Nanoparticles with less than 850 atoms show bistability between the solid and liquid phases over temperature ranges below the point of complete melting. The potential energy of a nanoparticle in the bistable region alternates between values corresponding to the solid and liquid phases. This bistability is characteristic of dynamic coexistence melting. At higher temperatures, only the liquid state is stable. Nanoparticles with more than 850 atoms undergo a sharp solid-liquid-phase transition characteristic of the bulk solid phase. The variation of the melting point with the effective nanoparticle radius is also determined.