Neutral and Charged 1-Butyl-3-methylimidazolium Triflate Clusters: Equilibrium Concentration in the Vapor Phase and Thermal Properties of Nanometric Droplets

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P ~ 0.01 Pa) and temperatures (T 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0 T 700 K range. A near-continuous transition between a liquidlike phase at high T and a solidlike phase at low T takes place at T ~ 190 K in close correspondence with the bulk glass point Tg ~ 200 K. Solidification is accompanied by a transition in the dielectric properties of the droplet, giving rise to a small permanent dipole embedded into the solid cluster. The simulation results highlight the molecular precursors of several macroscopic properties and phenomena and point to the close competition of Coulomb and dispersion forces as their common origin.

Local and semilocal density functional computations for crystals of 1-alkyl-3-methyl-imidazolium salts

The accuracy and reliability of popular density functional approximations for the compounds giving origin to room temperature ionic liquids have been assessed by computing the T=0 K crystal structure of several 1-alkyl-3-methyl-imidazolium salts. Two prototypical exchange-correlation approximations have been considered, i.e., the local density approximation (LDA) and one gradient corrected scheme [PBE-GGA, Phys. Rev. Lett. 77, 3865 (1996)]. Comparison with low-temperature x-ray diffraction data shows that the equilibrium volume predicted by either approximations is affected by large errors, nearly equal in magnitude (~10%), and of opposite sign. In both cases the error can be traced to a poor description of the intermolecular interactions, while the intramolecular structure is fairly well reproduced by LDA and PBE-GGA. The PBE-GGA optimization of atomic positions within the experimental unit cell provides results in good agreement with the x-ray structure. The correct system volume can also be restored by supplementing PBE-GGA with empirical dispersion terms reproducing the r-6 attractive tail of the van der Waals interactions.

Ion Association in [bmim][PF6]/Naphthalene Mixtures: An Experimental and Computational Study

Mixtures of room temperature ionic liquids (IL) with neutral organic molecules provide a valuable testing ground to investigate the interplay of the ionic and molecular-dipolar state in dense Coulomb systems at near ambient conditions. In the present study, the viscosity eta and the ionic conductivity a of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6])/naphthalene mixtures at T = 80 degrees C have been measured at 10 stoichiometries spanning the composition range from pure naphthalene to pure [bmim][PF6]. The viscosity grows nearly monotonically with increasing IL mole fraction (x), whereas the conductivity per ion displays a clear peak at x approximate to 15%. The origin of this maximum has been investigated using molecular dynamics simulations based on a classical force field. Snapshots of the simulated samples show that the conductivity maximum is due to the gradual transition in the IL component from an ionic state at high x to a dipolar fluid made of neutral ion pairs at low x. At concentrations x <0.20 the ion pairs condense into molecular-thin filaments bound by dipolar forces and extending in between nanometric droplets of IL. These results are confirmed and complemented by the computation of dynamic and transport properties in [bmim][PF6]/naphthalene mixtures at low IL concentration.

Spontaneous spin polarization and electron localization in constrained geometries: The Wigner transition in nanowires

The Wigner transition in a jellium model of cylindrical nanowires has been investigated by density-functional computations using the local spin-density approximation. A wide range of background densities rho(b) has been explored from the nearly ideal metallic regime (r(s)=[3/4 pi rho(b)](1/3)=1) to the high correlation limit (r(s)=100). Computations have been performed using an unconstrained plane wave expansion for the Kohn-Sham orbitals and a large simulation cell with up to 480 electrons. The electron and spin distributions retain the cylindrical symmetry of the Hamiltonian at high density, while electron localization and spin polarization arise nearly simultaneously in low-density wires (r(s)similar to 30). At sufficiently low density (r(s)>= 40), the ground-state electron distribution is the superposition of well defined and nearly disjoint droplets, whose charge and spin densities integrate almost exactly to one electron and 1/2 mu(B), respectively. Droplets are arranged on radial shells and define a distorted lattice whose structure is intermediate between bcc and fcc. Dislocations and grain boundaries are apparent in the droplets' configuration found by our simulations. Our computations aim at modeling the behavior of experimental low-carried density systems made of lightly doped semiconductor nanostructures or conducting polymers.

Interaction of Room Temperature Ionic Liquid Solutions with a Cholesterol Bilayer

Water solutions of representative (IC(4)mim][Cl] and [C(4)mim][Tf2N] room temperature ionic liquids (ILs) in contact with a neutral lipid bilayer made of cholesterol molecules has been investigated by molecular dynamics simulations based on an empirical force field model. The results show that both ILs display selective adsorption at the water-cholesterol interface, with partial inclusion of ions into the bilayer. In the case Of [C(4)mim][Cl], the adsorption of ions at the water-cholesterol interface is limited by a sizable bulk solubility of the IL, driven by the high water affinity of [Cl](-). The relatively low Solubility Of [C(4)mim][Tf2N], instead, gives rise to a nearly complete segregation of the IL component on the bilayer, altering its volume, compressibility, and electrostatic environment. The computational results display important similarities to the results of recent experimental measurements for ILs in contact with phospholipid model membranes (see Evans, K. O. Int. J. Mol. Sci. 2008, 9, 498-511 and references therein).

Structural motifs of cholesterol nanoparticles

The growth sequence of gas-phase cholesterol clusters (Ch(N)) with up to N=36 molecules has been investigated by atomistic simulation based on an empirical force field model. The results of long annealings from high temperature show that the geometric motifs characterizing the structure of pure cholesterol crystals already appear in nanometric aggregates. In all clusters molecules tend to align along a common direction. For cluster sizes above the smallest ones, dispersion interactions among the hydrocarbon body and tails of cholesterol cooperate with hydrogen bonding to give rise to a bilayer structure. Analysis of snapshots from the annealing shows that the condensation of hydrogen bonds into a connected network of rings and chains is an important step in the self-organization of cholesterol clusters. The effect of solvation on the equilibrium properties of medium-size aggregates is investigated by short molecular dynamics simulations for the N=30 and N=40 clusters in water at near ambient conditions and in supercritical carbon dioxide at T=400 K.

Structural changes in quasi-one-dimensional many-electron systems: From linear to zigzag and beyond

Many-electron systems confined to a quasi-one-dimensional geometry by a cylindrical distribution of positive charge have been investigated by density functional computations in the unrestricted local spin density approximation. Our investigations have been focused on the low-density regime, in which electrons are localized. The results reveal a wide variety of different charge and spin configurations, including linear and zig-zag chains, single-and double-strand helices, and twisted chains of dimers. The spin-spin coupling turns from weakly antiferromagnetic at relatively high density, to weakly ferromagnetic at the lowest densities considered in our computations. The stability of linear chains of localized charge has been investigated by analyzing the radial dependence of the self-consistent potential and by computing the dispersion relation of low-energy harmonic excitations.

Mesophases in Nearly 2D Room-Temperature Ionic Liquids

Computer simulations of (i) a [C(12)mim][Tf2N] film of nanometric thickness squeezed at kbar pressure by a piecewise parabolic confining potential reveal a mesoscopic in-plane density and composition modulation reminiscent of mesophases seen in 3D samples of the same room-temperature ionic liquid (RTIL). Near 2D confinement, enforced by a high normal load, as well as relatively long aliphatic chains are strictly required for the mesophase formation, as confirmed by computations for two related systems made of (ii) the same [C(12)mim][Tf2N] adsorbed at a neutral solid surface and (iii) a shorter-chain RTIL ([C(4)mim][Tf2N]) trapped in the potential well of part i. No in-plane modulation is seen for ii and iii. In case ii, the optimal arrangement of charge and neutral tails is achieved by layering parallel to the surface, while, in case iii, weaker dispersion and packing interactions are unable to bring aliphatic tails together into mesoscopic islands, against overwhelming entropy and Coulomb forces. The onset of in-plane mesophases could greatly affect the properties of long-chain RTILs used as lubricants.

Melting of a tetrahedral network model of silica

Thermal properties of an idealised tetrahedral network model of silica are investigated by Monte Carlo simulations. The interatomic potential consists of anharmonic stretching and bending terms, plus a short range repulsion. The model includes a bond interchange rule similar to the well known Wooten, Winer and Weaire (WWW) algorithm (see Phys. Rev. Lett., 1985, 54, 1392). Simulations reveal an apparent first order melting transition at T = 2200 K. The computed changes in the local coordination upon melting are consistent with experimental and ab initio data.

Spontaneous spin polarization and charge localization in metal nanowires: the role of a geometric constriction

An idealized jellium model of conducting nanowires with a geometric constriction is investigated by density functional theory (DFT) in the local spin density (LSD) approximation. The results reveal a fascinating variety of spin and charge patterns arising in wires of sufficiently low (r(s) >= 15) average electron density, pinned at the indentation by an apparent attractive interaction with the constriction. The spin-resolved frequency-dependent conductivity shows a marked asymmetry in the two spin channels, reflecting the spontaneous spin polarization around the wire neck. The relevance of the computational results is discussed in relation to the so-called 0.7 anomaly found by experiments in the low-frequency conductivity of nanowires at near-breaking conditions (see 2008 J. Phys.: Condens Matter 20, special issue on the 0.7 anomaly). Although our mean-field approach cannot account for the intrinsic many-body effects underlying the 0.7 anomaly, it still provides a diagnostic tool to predict impending transitions in the electronic structure.

Large-scale density functional calculations of the surface properties of the Wigner crystal

The surface properties of the jellium model have been investigated by large supercell computations in the density functional theory-local spin-density (DFT-LSD) approach for planar slabs with up to 1000 electrons. A wide interval of densities has been explored, extending into the stability range of the Wigner crystal. Most computations have been carried out on nominally paramagnetic samples with an equal number of spin-up and spin-down electrons. The results show that within DFT-LSD spontaneous spin polarization and charge localization start nearly simultaneously at the surface for r(s) similar to 20, then, with decreasing density, they progress toward the center of the slab. Electrons are fully localized and spin polarized at r(s) = 30. At this density the charge distribution is the superposition of disjoint charge blobs, each corresponding to one electron. The distribution of blobs displays both regularities and disorder, the first being represented by well-defined planes and simple in-plane geometries, and the latter by a variety of surface defects. The surface energy, surface dipole, electric polarisability, and magnetization pattern have been determined as a function of density. All these quantities display characteristic anomalies at the density of the localization transition. The analysis of the low-frequency electric conductivity shows that in the fluid paramagnetic regime the in-plane current preferentially flows in the central region of the slab and the two spin channels are equally conducting. In the charge localized, spin-polarized regime, conductivity is primarily a surface effect, and an apparent asymmetry is observed in the two spin currents.

Thermodynamic properties and atomistic structure of the dry amorphous silica surface from a reactive force field model

A force field model of the Keating type supplemented by rules to break, form, and interchange bonds is applied to investigate thermodynamic and structural properties of the amorphous SiO2 surface. A simulated quench from the liquid phase has been carried out for a silica sample made of 3888 silicon and 7776 oxygen atoms arranged on a slab similar to 40 angstrom thick, periodically repeated along two directions. The quench results into an amorphous sample, exposing two parallel square surfaces of similar to 42 nm(2) area each. Thermal averages computed during the quench allow us to determine the surface thermodynamic properties as a function of temperature. The surface tension turns out to be gamma=310 +/- 20 erg/cm(2) at room temperature and gamma=270 +/- 30 at T=2000 K, in fair agreement with available experimental estimates. The entropy contribution Ts-s to the surface tension is relatively low at all temperatures, representing at most similar to 20% of the surface energy. Almost without exceptions, Si atoms are fourfold coordinated and oxygen atoms are twofold coordinated. Twofold and threefold rings appear only at low concentration and are preferentially found in proximity of the surface. Above the glass temperature T-g=1660 +/- 50 K, the mobility of surface atoms is, as expected, slightly higher than that of bulk atoms. The computation of the height-height correlation function shows that the silica surface is rough in the equilibrium and undercooled liquid phase, becoming smooth below the glass temperature T-g.