Two-dimensional interlocked pentagonal bilayer ice: how do water molecules form a hydrogen bonding network?

The plethora of ice structures observed both in bulk and under nanoscale confinement reflects the extraordinary ability of water molecules to form diverse forms of hydrogen bonding networks. An ideal hydrogen bonding network of water should satisfy three requirements: (1) four hydrogen bonds connected with every water molecule, (2) nearly linear hydrogen bonds, and (3) tetrahedral configuration for the four hydrogen bonds around an O atom. However, under nanoscale confinement, some of the three requirements have to be unmet, and the selection of the specific requirement(s) leads to different types of hydrogen bonding structures. According to molecular dynamics (MD) simulations for water confined between two smooth hydrophobic walls, we obtain a phase diagram of three two-dimensional (2D) crystalline structures and a bilayer liquid. A new 2D bilayer ice is found and named the interlocked pentagonal bilayer ice (IPBI), because its side view comprises interlocked pentagonal channels. The basic motif in the top view of IPBI is a large hexagon composed of four small pentagons, resembling the top view of a previously reported "coffin" bilayer ice [Johnston, et al., J. Chem. Phys., 2010, 133, 154516]. First-principles optimizations suggest that both bilayer ices are stable. However, there are fundamental differences between the two bilayer structures due to the difference in the selection among the three requirements. The IPBI sacrifices the linearity of hydrogen bonds to retain locally tetrahedral configurations of the hydrogen bonds, whereas the coffin structure does the opposite. The tradeoff between the conditions of an ideal hydrogen bonding network can serve as a generic guidance to understand the rich phase behaviors of nanoconfined water.

Highly confined water: two-dimensional ice, amorphous ice, and clathrate hydrates.

Understanding phase behavior of highly confined water, ice, amorphous ice, and clathrate hydrates (or gas hydrates), not only enriches our view of phase transitions and structures of quasi-two-dimensional (Q2D) solids not seen in the bulk phases but also has important implications for diverse phenomena at the intersection between physical chemistry, cell biology, chemical engineering, and nanoscience. Relevant examples include, among others, boundary lubrication in nanofluidic and lab-on-a-chip devices, synthesis of antifreeze proteins for ice-growth inhibition, rapid cooling of biological suspensions or quenching emulsified water under high pressure, and storage of H2 and CO2 in gas hydrates. Classical molecular simulation (MD) is an indispensable tool to explore states and properties of highly confined water and ice. It also has the advantage of precisely monitoring the time and spatial domains in the sub-picosecond and sub-nanometer scales, which are difficult to control in laboratory experiments, and yet allows relatively long simulation at the 10(2) ns time scale that is impractical with ab initio molecular dynamics simulations. In this Account, we present an overview of our MD simulation studies of the structures and phase behaviors of highly confined water, ice, amorphous ice, and clathrate, in slit graphene nanopores. We survey six crystalline phases of monolayer (ML) ice revealed from MD simulations, including one low-density, one mid-density, and four high-density ML ices. We show additional supporting evidence on the structural stabilities of the four high-density ML ices in the vacuum (without the graphene confinement), for the first time, through quantum density-functional theory optimization of their free-standing structures at zero temperature. In addition, we summarize various low-density, high-density, and very-high-density Q2D bilayer (BL) ice and amorphous ice structures revealed from MD simulations. These simulations reinforce the notion that the nanoscale confinement not only can disrupt the hydrogen bonding network in bulk water but also can allow satisfaction of the ice rule for low-density and high-density Q2D crystalline structures. Highly confined water can serve as a generic model system for understanding a variety of Q2D materials science phenomena, for example, liquid-solid, solid-solid, solid-amorphous, and amorphous-amorphous transitions in real time, as well as the Ostwald staging during these transitions. Our simulations also bring new molecular insights into the formation of gas hydrate from a gas and water mixture at low temperature.

Identification of alcohol conformers by Raman spectra in the C-H stretching region.

The spontaneous polarized Raman spectra of normal and deuterated alcohols (C2-C5) have been recorded in the C-H stretching region. In the isotropic Raman spectra, a doublet of -CαH stretching vibration is found for all alcohols at below 2900 cm(-1) and above 2950 cm(-1). By comparing the experimental and calculated spectra of various deuterated alcohols, the doublets are attributed to the -CαH stretching vibration of different conformers. For ethanol, the band observed at 2970 cm(-1) is assigned as the stretching vibration of -CαH in the Cα-O-H plane of the gauche-conformer, while the band at 2895 cm(-1) is contributed from both the -CαH2 symmetrical stretching vibration of the trans-conformer and the -CαH stretching vibration out of the Cα-O-H plane of the gauche-conformer. The population of gauche-conformer is estimated to be 54% in liquid ethanol. For the larger alcohols, the same assignments for the doublet are obtained, and the populations of gauche-conformers with plane carbon skeleton are found to be slightly larger than that of ethanol, which is consistent with results from molecular dynamics simulations.

Highly selective adsorption of methanol in carbon nanotubes immersed in methanol-water solution.

The systems of open-ended carbon nanotubes (CNTs) immersed in methanol-water solution are studied by molecular dynamics simulations. For the (6,6) CNT, nearly pure methanol is found to preferentially occupy interior space of the CNT. Even when the mass fraction (MF) of methanol in bulk solution is as low as 1%, the methanol MF within the CNT is still more than 90%. For CNTs with larger diameters, the methanol concentrations within CNTs are also much higher than those outside CNTs. The methanol selectivity decreases with increasing CNT diameter, but not monotonically. From microscopic structural analyses, we find that the primary reason for the high selectivity of methanol by CNTs lies on high preference of methanol in the first solvation shell near the inner wall of CNT, which stems from a synergy effect of the van der Waals interaction between CNT and the methyl groups of methanol, together with the hydrogen bonding interaction among the liquid molecules. This synergy effect may be of general significance and extended to other systems, such as ethanol aqueous solution and methanol/ethanol mixture. The selective adsorption of methanol over water in CNTs may find applications in separation of water and methanol, detection of methanol, and preservation of methanol purity in fuel cells.

Structural selection of graphene supramolecular assembly oriented by molecular conformation and alkyl chain.

Graphene molecules, hexafluorotribenzo[a,g,m]coronene with n-carbon alkyl chains (FTBC-Cn, n = 4, 6, 8, 12) and Janus-type "double-concave" conformation, are used to fabricate self-assembly on highly oriented pyrolytic graphite surface. The structural dependence of the self-assemblies with molecular conformation and alkyl chain is investigated by scanning tunneling microscopy and density functional theory calculation. An interesting reverse face "up-down" way is observed in FTBC-C4 assembly due to the existence of hydrogen bonds. With the increase of the alkyl chain length and consequently stronger van der Waals interaction, the molecules no longer take alternating "up-down" orientation in their self-assembly and organize into various adlayers with lamellar, hexagonal honeycomb, and pseudohoneycomb structures based on the balance between intermolecular and molecule-substrate interactions. The results demonstrate that the featured "double-concave" molecules are available block for designing graphene nanopattern. From the results of scanning tunneling spectroscopy measurement, it is found that the electronic property of the featured graphene molecules is preserved when they are adsorbed on solid surface.

Ab initio prediction of amorphous B84.

To explore the possible existence of boron clusters without carbon analogs, we study B(84) cluster as a prototypical system by ab initio calculations. Structures of several isomer forms of B(84) are optimized. Among these isomers, a group of amorphous (disordered) structures are found to be the most stable. Different from the high-symmetry isomers, the amorphous B(84) clusters are more stable than the fullerene B(80) in terms of cohesive energy per atom. These amorphous structures can be distinguished from other high-symmetry structures experimentally via, for example, infrared spectra. The radial and angular distribution functions of amorphous B(84) structures are more diffuse than those of high-symmetry structures. On the basis of these findings, we propose that amorphous structures may be generic for boron and dominate boron clusters in a range of cluster scale.

Room temperature electrofreezing of water yields a missing dense ice phase in the phase diagram.

Water can freeze into diverse ice polymorphs depending on the external conditions such as temperature (T) and pressure (P). Herein, molecular dynamics simulations show evidence of a high-density orthorhombic phase, termed ice χ, forming spontaneously from liquid water at room temperature under high-pressure and high external electric field. Using free-energy computations based on the Einstein molecule approach, we show that ice χ is an additional phase introduced to the state-of-the-art T-P phase diagram. The χ phase is the most stable structure in the high-pressure/low-temperature region, located between ice II and ice VI, and next to ice V exhibiting two triple points at 6.06 kbar/131.23 K and 9.45 kbar/144.24 K, respectively. A possible explanation for the missing ice phase in the T-P phase diagram is that ice χ is a rare polarized ferroelectric phase, whose nucleation/growth occurs only under very high electric fields.

First-principles study of electronic and magnetic properties of Co(n)Mn(m) and Co(n)V(m) (m + n < or = 6) clusters.

The electronic and magnetic properties of small Co(n)Mn(m) and Co(n)V(m) (m + n < or = 6) clusters are systematically studied using density functional theory. The results show that Co and V atoms prefer to aggregate in Co-Mn and Co-V clusters, respectively. Significant magnetic moment enhancement in Co-Mn clusters with Mn doping and reduction in Co-V clusters with V doping are found, consistent with experiment results for larger clusters [Phys. Rev. Lett. 2007, 98, 113401]. The results are discussed by analyzing the magnetic coupling type and local magnetic moment on each atoms. Density of states and vertical ionization potentials are calculated and show cluster size dependent behavior.

From isosuperatoms to isosupermolecules: new concepts in cluster science.

As an extension of the superatom concept, a new concept "isosuperatom" is proposed, reflecting the physical phenomenon that a superatom cluster can take multiple geometrical structures with their electronic structures topologically invariant. The icosahedral and cuboctahedral Au13(5+) units in the Au25(SCH2CH2Ph)18(-), Au23(SC6H11)16(-) and Au24(SAdm)16 nanoclusters are found to be examples of this concept. Furthermore, two isosuperatoms can combine to form a supermolecule. For example, the structure of the {Ag32(DPPE)5(SC6H4CF3)24}(2-) nanocluster can be understood well in terms of a Ag22(12+) supermolecule formed by two Ag13(8+) isosuperatoms. On the next level of complexity, various combinations of isosuperatoms can lead to supermolecules with different geometrical structures but similar electronic structures, i.e., "isosupermolecules". We take two synthesized nanoclusters Au20(PPhpy2)10Cl4(2+) and Au30S(StBu)18 to illustrate two Au20(6+) isosupermolecules. The proposed concepts of isosuperatom and isosupermolecule significantly enrich the superatom concept, give a new framework for understanding a wide range of nanoclusters, and open a new door for designing assembled materials.

Icosahedral B12-containing core-shell structures of B80.

Low-lying icosahedral (I(h)) B(12)-containing structures of B(80) are explored, and a number of core-shell isomers are found to have lower energy than the previous predicted B(80) fullerene. The structural transformation of boron clusters from tubular structure to core-shell structure may occur at a critical size less than B(80).

Electron impact ionization of water-doped superfluid helium nanodroplets: observation of He(H(2)O)(n)(+) clusters.

Electron impact mass spectra have been recorded for helium nanodroplets containing water clusters. In addition to identification of both H(+)(H(2)O)(n) and (H(2)O)(n)(+) ions in the gas phase, additional peaks are observed which are assigned to He(H(2)O)(n)(+) clusters for up to n=27. No clusters are detected with more than one helium atom attached. The interpretation of these findings is that quenching of (H(2)O)(n)(+) by the surrounding helium can cool the cluster to the point where not only is fragmentation to H(+)(H(2)O)(m) (where m < or = n-1) avoided, but also, in some cases, a helium atom can remain attached to the cluster ion as it escapes into the gas phase. Ab initio calculations suggest that the first step after ionization is the rapid formation of distinct H(3)O(+) and OH units within the (H(2)O)(n)(+) cluster. To explain the formation and survival of He(H(2)O)(n)(+) clusters through to detection, the H(3)O(+) is assumed to be located at the surface of the cluster with a dangling O-H bond to which a single helium atom can attach via a charge-induced dipole interaction. This study suggests that, like H(+)(H(2)O)(n) ions, the preferential location for the positive charge in large (H(2)O)(n)(+) clusters is on the surface rather than as a solvated ion in the interior of the cluster.

Low-temperature orientationally ordered structures of two-dimensional C60.

Orientationally ordered structures of two-dimensional (2D) C(60) at low temperature have been investigated theoretically and experimentally. Using total energy optimization with a phenomenological potential, we find the ground state is a close packed hexagonal lattice in which all the molecules have the same orientation. Several local minima of the potential energy surface are found to be associated with other 1 x 1 lattices as well as 2 x 2 lattices. The energies of the orientational domain boundaries of the 1x1 lattices are also computed, and two kinds of which yield negative values. A majority of these theoretical findings are confirmed by our low-temperature scanning tunneling microscopy study of a 2D C(60) array supported on a self-assembled monolayer.