Observation of two-step nucleation in methane hydrates

In this work we show that homogeneous nucleation of methane hydrate can, under appropriate conditions, be a very rapid process, achieved within tens of nanoseconds. In agreement with recent experimental results on different systems, we find that the nucleation of a gas hydrate crystal appears as a two-step process. It starts with the formation of disordered solid-like structures, which will then spontaneously evolve to more recognizable crystalline forms. This previously elusive first-stage state is confirmed to be post-critical in the nucleation process, and is characterized as processing reasonable short-range structure but essentially no long-range order. Its energy, molecular diffusion and local structure reflect a solid-like character, although it does exhibit mobility over longer (tens of ns) timescales. We provide insights into the controversial issue of memory effects in methane hydrates. We show that areas locally richer in methane will nucleate much more readily, and no 'memory' of the crystal is required for fast re-crystallization. We anticipate that much richer polycrystallinity and novel methane hydrate phases could be possible.     

Heterogeneous crystal growth of methane hydrate on its sII [001] crystallographic face

This paper presents a systematic molecular simulation study of the heterogeneous crystal growth of methane hydrate sII from supersaturated aqueous methane solutions. The growth of sII hydrate on the [001] crystallographic face is achieved through utilization of a recently proposed methodology, and rates of crystal growth of 1 A/ns were sustained for the molecular models and specific conditions employed in this work. Characteristics of the crystals grown as well as properties and structure of the interface are examined. Water cages with a 5(12)6(3) arrangement, which are improper to both sI and sII structures, are identified during the heterogeneous growth of sII methane hydrate. We show that the growth of a [001] face of sII hydrate can produce an sI crystalline structure, confirming that cross-nucleation of methane hydrate structures is possible. Defects consisting of two methane molecules trapped in large 5(12)6(4) cages and water molecules trapped in small and large cages are observed, where in one instance we have found a large 5(12)6(4) cage containing three water molecules.     

Unusual crystalline and polycrystalline structures in methane hydrates

We present previously unreported crystalline and polycrystalline structures for methane hydrates obtained at relatively low pressures. These structures, which contain unusual cages with 12 pentagonal faces and 3 hexagonal faces, were observed during the atomistic simulations of the crystal growth of sI and sII methane hydrates. These 51263 cages have a significant impact on the structure of the resulting crystal and could explain several experimental observations regarding in-situ transformations between sI and sII hydrates. We document a previously unidentified structure of methane hydrates which we designate structure sK. Additionally, we predict a polycrystalline structure consisting of this new hydrate and sI and suggest a mechanism for the formation of a polycrystalline structure consisting of sequences of sI and sII hydrates.     

Molecular insights into the heterogeneous crystal growth of si methane hydrate

In this paper we report a successful molecular simulation study exploring the heterogeneous crystal growth of sI methane hydrate along its [001] crystallographic face. The molecular modeling of the crystal growth of methane hydrate has proven in the past to be very challenging, and a reasonable framework to overcome the difficulties related to the simulation of such systems is presented. Both the microscopic mechanisms of heterogeneous crystal growth as well as interfacial properties of methane hydrate are probed. In the presence of the appropriate crystal template, a strong tendency for water molecules to organize into cages around methane at the growing interface is observed; the interface also demonstrates a strong affinity for methane molecules. The maximum growth rate measured for a hydrate crystal is about 4 times higher than the value previously determined for ice I in a similar framework (Gulam Razul, M. S.; Hendry, J. G.; Kusalik, P. G. J. Chem. Phys. 2005, 123, 204722).     

Molecular dynamics methodology to investigate steady-state heterogeneous crystal growth

In this paper a new molecular dynamics simulation methodology to investigate steady-state heterogeneous crystal growth from a supercooled liquid is presented. The method is tested on pure component systems such as Lennard-Jonesium and water/ice, as well as multicomponent systems such as methane hydrate crystals. The setup uses periodicity in all three directions and two interfaces; at one interface, crystallization occurs, while at the other, melting is enforced by locally heating the crystal only near that interface. Steady-state conditions are achieved when the crystal is melted at the same rate as the growth occurs. A self-adaptive scheme that automatically modifies the rate of melting to match the rate of growth, crucial for establishing steady-state conditions, is described. In contrast with the recently developed method of Razul et al. [Mol. Phys. 103, 1929 (2005)], where the rates of growth (melting) were constant and the temperatures determined, the present approach fixes the supercooling temperature at the growing interface and identifies the corresponding steady-state crystal growth rate that corresponds to the thermodynamic force provided. The static properties of the interface (e.g., the interfacial widths) and the kinetics of the crystal growth are found to reproduce well previous findings. The importance of establishing steady-state conditions in such investigations is also briefly discussed.     

Molecular insights into the potential and temperature dependences of the differential capacitance of a room-temperature ionic liquid at graphite electrodes

Molecular dynamics simulation studies of the structure and the differential capacitance (DC) for the ionic liquid (IL) N-methyl-N-propylpyrrolidinium bis(trifluoromethane)sulfonyl imide ([pyr(13)][TFSI]) near a graphite electrode have been performed as a function temperature and electrode potential. The IL exhibits a multilayer structure that extends 20-30 ? from the electrode surface. The composition and ion orientation in the innermost layer were found to be strongly dependent on the electrode potential. While at potentials near the potential of zero charge (PZC), both cations and anions adjacent to the surface are oriented primarily perpendicular to the surface, the counterions in first layer orient increasingly parallel to the surface with increasing electrode potential. A minimum in DC observed around -1 V(RPZC) (potential relative to the PZC) corresponds to the point of highest density of perpendicularly aligned TFSI near the electrode. Maxima in the DC observed around +1.5 and -2.5 V(RPZC) are associated with the onset of "saturation", or crowding, of the interfacial layer. The asymmetry of DC versus electrode polarity is the result of strong interactions between the fluorine of TFSI and the surface, the relatively large footprint of TFSI compared to pyr(13), and the tendency of the propyl tails of pyr(13) to remain adsorbed on the surface even at high positive potentials. Finally, an observed decreased DC and the disappearance of the minimum in DC near the PZC with increasing temperature are likely due to the increasing importance of entropic/excluded volume effects (interfacial crowding) with increasing temperature.     

Evaluation of site-site bridge diagrams for molecular fluids

The presence of bridge functions in formally exact integral equation theories is the primary obstacle preventing the extraction of exact fluid structure from these theories. The bridge functions are typically neglected but in many fluids their impact may be significant. Each bridge function can be subdivided into bridge diagrams, which are well defined but difficult to evaluate. The calculation of bridge diagrams for the Chandler-Silbey-Ladanyi (CSL) integral equation theory is the subject of this paper. In particular, we evaluate the diagrams required to yield an exact theory up to the first power in density [O(rho(1))] and provide algorithms that remain feasible for any molecule. Further, the bridge diagrams are evaluated and compared with the f-bond and h-bond formulations. Exact bridge diagrams are numerically evaluated for several chiral molecules, for two polar dimers, and for SPC/E water. The quality of the diagrams is assessed in two ways: First, the predicted interatomic distributions are compared with those obtained from Monte Carlo simulations. Second, the connectivity constraints are evaluated and the errors in satisfying these exact relationships are compared for the f-bond and h-bond formulations. For apolar fluids, a clear improvement in CSL theory is evident with the inclusion of O(rho(0)) and O(rho(1)) diagrams. In contrast, for polar fluids, the inclusion of bridge diagrams does not lead to improvement in the structural predictions.     

Inhibiting manganese(Ⅱ)from catalyzing electrolyte decomposition in lithium-ion batteries

A once overlooked source of electrolyte degradation incurred by dissolved manganese(Ⅱ)species in lithium-ion batteries has been identified recently.In order to deactivate the catalytic activity of such manganese(II)ion,1-aza-12-crown-4-ether(A12C4)with cavity size well matched manganese(Ⅱ)ion is used in this work as electrolyte additive.Theoretical and experimental results show that stable complex forms between A12C4 and manganese(II)ions in the electrolyte,which does not affect the solvation of Li ions.The strong binding effect of A12C4 additive reduces the charge density of manganese(II)ion and inhibits its destruction of the PF₆^-
structure in the electrolyte,leading to greatly improved thermal stability of manganese(II)ions-containing electrolyte.In addition to bulk electrolyte,A12C4 additive also shows capability in preventing Mn²⁺ from degrading SEI on graphite surface.Such bulk and interphasial stability introduced by A12C4 leads to significantly improved cycling performance of LIBs.