Mal-Soon Lee

High-pressure polymeric phases of carbon dioxide

Understanding the structural transformations of solid CO2 from a molecular solid characterized by weak intermolecular bonding to a 3-dimensional network solid at high pressure has challenged researchers for the past decade. We employ the recently developed metadynamics method combined with ab initio calculations to provide fundamental insight into recent experimental reports on carbon dioxide in the 60–80 GPa pressure region. Pressure-induced polymeric phases and their transformation mechanisms are found. Metadynamics simulations starting from the CO2-II (P42/mnm) at 60 GPa and 600 K proceed via an intermediate, partially polymerized phase, and finally yield a fully tetrahedral, layered structure (P-4m2). Based on the agreement between calculated and experimental Raman and X-ray patterns, the recently identified phase VI [Iota V, et al. (2007) Sixfold coordinated carbon dioxide VI. Nature Mat 6:34–38], assumed to be disordered stishovite-like, is instead interpreted as the result of an incomplete transformation of the molecular phase into a final layered structure. In addition, an α-cristobalite-like structure (P41212), is predicted to be formed from CO2-III (Cmca) via an intermediate Pbca structure at 80 GPa and low temperatures (<300 K). Defects in the crystals are frequently observed in the calculations at 300 K whereas at 500 to 700 K, CO2-III transforms to an amorphous form, consistent with experiment [Santoro M, et al. (2006) Amorphous silica-like carbon dioxide. Nature 441:857–860], but the simulation yields additional structural details for this disordered solid.

First-principles investigation of finite-temperature behavior in small sodium clusters.

A systematic and detailed investigation of the finite-temperature behavior of small sodium clusters, Na(n), in the size range of n=8-50 are carried out. The simulations are performed using density-functional molecular dynamics with ultrasoft pseudopotentials. A number of thermodynamic indicators such as specific heat, caloric curve, root-mean-square bond-length fluctuation, deviation energy, etc., are calculated for each of the clusters. Size dependence of these indicators reveals several interesting features. The smallest clusters with n=8 and 10 do not show any signature of melting transition. With the increase in size, broad peak in the specific heat is developed, which alternately for larger clusters evolves into a sharper one, indicating a solidlike to liquidlike transition. The melting temperatures show an irregular pattern similar to the experimentally observed one for larger clusters [Schmidt et al., Nature (London) 393, 238 (1998)]. The present calculations also reveal a remarkable size-sensitive effect in the size range of n=40-55. While Na(40) and Na(55) show well-developed peaks in the specific-heat curve, Na(50) cluster exhibits a rather broad peak, indicating a poorly defined melting transition. Such a feature has been experimentally observed for gallium and aluminum clusters [Breaux et al., J. Am. Chem. Soc. 126, 8628 (2004); Breaux et al., Phys. Rev. Lett. 94, 173401 (2005)].

Aluminum Zintl anion moieties within sodium aluminum clusters.

Through a synergetic combination of anion photoelectron spectroscopy and density functional theory based calculations, we have established that aluminum moieties within selected sodium-aluminum clusters are Zintl anions. Sodium-aluminum cluster anions, Na(m)Al(n)(-), were generated in a pulsed arc discharge source. After mass selection, their photoelectron spectra were measured by a magnetic bottle, electron energy analyzer. Calculations on a select sub-set of stoichiometries provided geometric structures and full charge analyses for both cluster anions and their neutral cluster counterparts, as well as photodetachment transition energies (stick spectra), and fragment molecular orbital based correlation diagrams.

Far-infrared absorption of water clusters by first-principles molecular dynamics.

Based on first-principle molecular dynamic simulations, we calculate the far-infrared spectra of small water clusters (H(2)O)(n) (n = 2, 4, 6) at frequencies below 1000 cm(-1) and at 80 K and at atmospheric temperature (T>200 K). We find that cluster size and temperature affect the spectra significantly. The effect of the cluster size is similar to the one reported for confined water. Temperature changes not only the shape of the spectra but also the total strength of the absorption, a consequence of the complete anharmonic nature of the classical dynamics at high temperature. In particular, we find that in the frequency region up to 320 cm(-1), the absorption strength per molecule of the water dimer at 220 K is significantly larger than that of bulk liquid water, while tetramer and hexamer show bulklike strengths. However, the absorption strength of the dimer throughout the far-infrared region is too small to explain the measured vapor absorption continuum, which must therefore be dominated by other mechanisms.

The effects of electronic structure and charged state on thermodynamic properties : An ab initio molecular dynamics investigations on neutral and charged clusters of Na39, Na40, and Na41

In this paper we explore the effects of the electronic structure, the charge state, and the nature of energy distribution of isomers on the thermodynamic properties of sodium clusters. The focus of the work is to isolate the effects of these ingredients on thermodynamic behavior by choosing specific clusters. Toward this end we investigate Na(39) (-), Na(40), and Na(41) (+), which are the electronic closed shell systems which differ in number of atoms and charge state. We also examine Na(39), Na(39) (+), Na(40) (+), and Na(41) clusters having different charges of these clusters. Our density functional molecular dynamics simulations show that all electronic shell-closing clusters have similar melting temperature of approximately 310 K. Remarkably, it is observed that an addition of even one electron to Na(39) increases the melting temperature by about 40 K and makes the specific heat curve sharper. All the cationic clusters show broadened specific heat curves.

Electronic structures, equilibrium geometries, and finite-temperature properties of Na n (n=39-55) from first principles

Density-functional theory has been applied to investigate systematics of sodium clusters Na n in the size range of n=39-55. A clear evolutionary trend in the growth of their ground-state geometries emerges. The clusters at the beginning of the series (n=39-43) are symmetric and have partial icosahedral (two-shell) structure. The growth then goes through a series of disordered clusters (n=44-52) where the icosahedral core is lost. However, for n≥53, a three-shell icosahedral structure emerges. This change in the nature of the geometry is abrupt. In addition, density-functional molecular dynamics has been used to calculate the specific heat curves for the representative sizes n=43, 45, 48, and 52. These results along with already available thermodynamic calculations for n=40, 50, and 55 enable us to carry out a detailed analysis of the heat capacity curves and their relationship with respective geometries for the entire series. Our results clearly bring out strong correlation between the evolution of the geometries and the nature of the shape of the heat capacities. The results also firmly establish the size-sensitive nature of the heat capacities in sodium clusters.

Effects of geometric and electronic structure on the finite temperature behavior of Na 58 , Na 57 , and Na 55 cluster

We report the equilibrium geometries and the electronic structures of Na n clusters in the size range of n=55-62 using density-functional method. An analysis of the evolutionary trends in their ground state geometries reveals that Na 58 has a spherical shape which is driven by the closed-shell nature of the electronic structure. This structure shows a significant large network connected by short bonds among the surface atoms as well as between core and surface atoms, which affects its finite-temperature behavior. By employing ab initio density-functional molecular dynamics, we calculate the specific heat of Na 58 and Na 57 . We observe two distinct features in their specific-heat curves as compared to that of Na 55 : (1) Both clusters show very broad melting transition. (2) The calculated melting temperature of Na 58 is ~375 K, the highest one studied so far, and that of Na 57 is also relatively high (~350 K). Thus, when a cluster has a (nearly) geometric closed-shell structure as well as a (nearly) electronic closed-shell one, it shows a high melting temperature. Our calculations clearly bring out the size-sensitive nature of the specific-heat curve in sodium clusters.

Mixtures of planetary ices at extreme conditions

The interiors of Neptune and Uranus are believed to be primarily composed of a fluid mixture of methane and water. The mixture is subjected to pressures up to several hundred gigapascal, causing the ionization of water. Laboratory and simulation studies so far have focused on the properties of the individual components. Here we show, using first-principle molecular dynamic simulations, that the properties of the mixed fluid are qualitatively different with respect to those of its components at the same conditions. We observe a pressure-induced softening of the methane-water intermolecular repulsion that points to an enhancement of mixing under extreme conditions. Ionized water causes the progressive ionization of methane and the mixture becomes electronically conductive at milder conditions than pure water, indicating that the planetary magnetic field of Uranus and Neptune may originate at shallower depths than currently assumed.