Temperature dependence of the dynamic viscosity of liquid mercury

The dependence of the dynamic viscosity of mercury on temperature is calculated and expressed in terms of the cluster associate model, based on the Boltzmann distribution and with normalization at the melting point. The resulting refined equation for mercury viscosity adequately reflects this dependence over the range of the liquid state, including the critical point.     

Melting and freezing characteristics and structural properties of supported and unsupported gold nanoclusters

Molecular dynamics simulations in conjunction with MEAM potential models have been used to study the melting and freezing behavior and structural properties of both supported and unsupported Au nanoclusters within a size range of 2 to 5 nm. In contrast to results from previous simulations regarding the melting of free Au nanoclusters, we observed a structural transformation from the initial FCC configuration to an icosahedral structure at elevated temperatures followed by a transition to a quasimolten state in the vicinity of the melting point. During the freezing of Au liquid clusters, the quasimolten state reappeared in the vicinity of the freezing point, playing the role of a transitional region between the liquid and solid phases. In essence, the melting and freezing processes involved the same structural changes which may suggest that the formation of icosahedral structures at high temperatures is intrinsic to the thermodynamics of the clusters, rather than reflecting a kinetic phenomenon. When Au nanoclusters were deposited on a silica surface, they transformed into icosahedral structures at high temperatures, slightly deformed due to stress arising from the Au-silica interface. Unlike free Au nanoclusters, an icosahedral solid-liquid coexistence state was found in the vicinity of the melting point, where the cluster consisted of coexisting solid and liquid fractions but retained an icosahedral shape at all times. These results demonstrated that the structural stability in the structures of small Au nanoclusters can be enhanced through interaction with the substrate. Supported Au nanoclusters demonstrated a structural transformation from decahedral to icosahedral motifs during Au island growth, in contrast to the predictions of the minimum-energy growth sequence: icosahedral structures appear first at very small cluster sizes, followed by decahedral structures, and finally FCC structures recovered at very large cluster sizes. The simulations also showed that island shapes are strongly influenced by the substrate, more specifically, the structural characteristic of a Au island is not only a function of size, but also depends on the contact area with the surface, which is controlled by the wetting of the cluster to the substrate.     

Molecular dynamics simulations of the melting of aluminum nanoparticles

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.     

An investigation of melting/freezing characteristics of nanoparticle-enhanced phase change materials

The aim of this study is to investigate the melting/freezing characteristics of paraffin by adding Cu nanoparticles. Cu/paraffin composite phase change materials (PCMs) were prepared by a two-step method. The effects of Cu nanoparticles on the thermal conductivity and the phase change heat transfer of PCMs were investigated by the Hot Disk thermal constants analyzer and infrared monitoring methods, respectively. The maximum thermal conductivity enhancements up to 14.2% in solid state and 18.1% in liquid state are observed at the 2?wt% Cu/paraffin. The photographs of infrared monitoring suggest that the melting and freezing rates of Cu/paraffin are enhanced. For 1?wt% Cu/paraffin, the melting and freezing times can be saved by about 33.3 and 31.6%, respectively. The results provide that adding nanoparticles is an efficient way to enhance the phase change heat transfer of PCMs.     

Transferable Pair Potentials for Liquid Iron,Cobalt and Nickel

Effective inter-ionic pair potentials for liquid iron, cobalt and nickel are derived from second order pseudopotential perturbation theory with the transferable electron-ion potential of Fiolhais and co-workers which was originally developed for the solid state. The liquid structure is obtained by the thermodynamically self consistent liquid state theory, the variational modified hypernetted chain (VMHNC) approximation. It has been found that the calculated structural properties at the investigated thermodynamic states just above the melting point agree well with experiment. In this work we have determined the parameters of liquids Fe, Co and Ni for the universal choice of the evanescent core pseudopotentials of Fiolhais'. We have shown that this pseudopotential is transferable to the liquid state if used this parameterisation.     

Metal clusters that freeze into high energy geometries

Heat capacities measured for isolated aluminum clusters show peaks due to melting. For some clusters with around 60 and 80 atoms there is a dip in the heat capacities at a slightly lower temperature than the peak. The dips have been attributed to structural transitions. Here we report studies where the clusters are annealed before the heat capacity is measured. The dips disappear for some clusters, but in many cases they persist, even when the clusters are annealed to well above their melting temperature. This indicates that the dips do not result from badly formed clusters generated during cluster growth, as originally suggested. We develop a simple kinetic model of melting and freezing in a system consisting of one liquidlike and two solidlike states with different melting temperatures and latent heats. Using this model we are able to reproduce the experimental results including the dependence on the annealing conditions. The dips result from freezing into a high energy geometry and then annealing into the thermodynamically preferred solid. The thermodynamically preferred solid has the higher freezing temperature. However, the liquid can bypass freezing into the thermodynamically preferred solid (at high cooling rates) if the higher energy geometry has a larger freezing rate.     

Simulations of the solid, liquid, and melting of 1-n-butyl-4-amino-1,2,4-triazolium bromide

Molecular dynamics simulations are used to study the solid and liquid properties and to predict the melting point of 1-n-propyl-4-amino-1,2,4-triazolium bromide ([patr][Br]) using a force field based on the one developed by Canongia Lopes et al. (J. Phys. Chem. B 2004, 108, 2038) for dialkyl substituted imidazolium salts, which was modified by including terms from the general AMBER force field. Electrostatic charges for the intermolecular interactions were determined from gas-phase ab initio electron structure calculations of the triazolium cation. Simulations of the solid state at 100 K reproduced the experimental density to within 4%. Simulations from 100 K to the melting point and the liquid from 333 to 500 K were performed to determine the temperature dependence of the densities of the two phases. The structures of the solid and liquid phases are characterized with radial distribution functions, which show that there are strong spatial correlations among neighboring ion pairs in liquid [patr][Br]. The dynamic behavior of the ions in the liquid state is also studied by computing velocity autocorrelation functions and the mean-square displacements between the ions. The melting point is determined by simulating void-induced melting. Changes in the density, intermolecular energy, and Lindemann index are used as indicators of the melting transition. The computed melting point is 360 +/- 10 K, which is within 10% of the experimental value 333 K.     

Crystallization, melting, and structure of water nanoparticles at atmospherically relevant temperatures

Water nanoparticles play an important role in atmospheric processes, yet their equilibrium and nonequilibrium liquid-ice phase transitions and the structures they form on freezing are not yet fully elucidated. Here we use molecular dynamics simulations with the mW water model to investigate the nonequilibrium freezing and equilibrium melting of water nanoparticles with radii R between 1 and 4.7 nm and the structure of the ice formed by crystallization at temperatures between 150 and 200 K. The ice crystallized in the particles is a hybrid form of ice I with stacked layers of the cubic and hexagonal ice polymorphs in a ratio approximately 2:1. The ratio of cubic ice to hexagonal ice is insensitive to the radius of the water particle and is comparable to that found in simulations of bulk water around the same temperature. Heating frozen particles that contain multiple crystallites leads to Ostwald ripening and annealing of the ice structures, accompanied by an increase in the amount of ice at the expense of the liquid water, before the particles finally melt from the hybrid ice I to liquid, without a transition to hexagonal ice. The melting temperatures T(m) of the nanoparticles are not affected by the ratio of cubic to hexagonal layers in the crystal. T(m) of the ice particles decreases from 255 to 170 K with the particle size and is well described by the Gibbs-Thomson equation, T(m)(R) = T(m)(bulk) - K(GT)/(R - d), with constant K(GT) = 82 ± 5 K·nm and a premelted liquid of width d = 0.26 ± 0.05 nm, about one monolayer. The freezing temperatures also decrease with the particles' radii. These results are important for understanding the composition, freezing, and melting properties of ice and liquid water particles under atmospheric conditions.     

Ionic microgels as model systems for colloids with an ultrasoft electrosteric repulsion: structure and thermodynamics

We present a theoretical analysis of the structural properties and phase behavior of spherical, loosely cross-linked ionic microgels that possess a low monomer concentration. The analysis is based on the recently derived effective interaction potential between such particles [A. R. Denton, Phys. Rev. E 67, 011804 (2003)]. By employing standard tools from the theory of the liquid state, we quantitatively analyze the pair correlations in the fluid and find anomalous behavior above the overlap concentration, similar to the cases of star-branched neutral and charged polymers. We also employ an evolutionary algorithm in order to predict the crystalline phases of the system without any a priori assumptions regarding their symmetry class. A very rich phase diagram is obtained, featuring two reentrant melting transitions and a number of unusual crystal structures. At high densities, both the Hansen-Verlet freezing criterion [J.-P. Hansen and L. Verlet, Phys. Rev. 184, 151 (1969)] and the Lindemann melting criterion [F. A. Lindemann, Phys. Z. 11, 609 (1910)] lose their validity. The topology of the phase diagram is altered when the steric interactions between the polymer segments become strong enough, in which case the lower-density reentrant melting disappears and the region of stability of the fluid is split into two disconnected domains, separated by intervening fcc and bcc regions.     

Anomalous freezing and melting behaviour of capillary confined CO2

In this paper we present our recent positron annihilation study of the liquid»solid phase boundary for CO₂ confined in nanometer pores of VYCOR glass. We find that CO₂ remains liquid in the pores far below the bulk freezing temperature and there is pronounced hysteresis between freezing and melting compared to that seen at the gas-liquid boundary in the pores. On freezing we see evidence of open space created in the pores. This leads to complex melting behaviour possibly involving the formation of gas-liquid interfaces. We see that frezing in the pores is totally irreversible, so that any solid which forms (no matter how small) remains stable up to the higher melting temperature. In contrast melting is more reversible (possibly indicating nucleation centres which permit immediate re-freezing). Finally, the pre-frozen state in the pores is different to the post-melted state.     

Equation of state for mercury: revisited

An analytical equation of state by Song and Mason is developed to calculate the PVT properties of mercury. The equation of state is based on the statistical-mechanical perturbation theory of hard convex bodies and can be written as a fifth-order polynomial in the density. There exists three temperature-dependent parameters in the equation of state; the second virial coefficient, an effective molecular volume, and a scaling factor for the average contact pair distribution function of hard convex bodies. The temperature-dependant parameters have been calculated using corresponding-states correlations based on the surface tension and the liquid density at the normal boiling point. Employing the present equation of state, we have calculated the PVT properties of mercury over a wide range of temperatures and pressures. The average absolute deviation for the calculated density of mercury in the saturation and compressed states is 1.09%.     

Thermodynamic properties for liquid mercury using GMA equation of state

The important known regularities and thermodynamic properties of liquid mercury have been studied based on the average potential energy. Recognised regularities, the linearity of Zeno contour, bulk modulus and secant bulk modulus as functions of temperature, isochors of pressure versus temperature and near linearity of the inverse isobaric expansion coefficient have been investigated, all evaluated using the Goharshadi–Morsali–Abbaspour equation of state. The validity of the equation of state in predicting thermophysical properties is confirmed by a statistical parameter, absolute average deviation, with a maximum value of 0.41, showing excellent agreement with the experiment at temperatures between 293.15 and 323.15?K from low to high pressures.     

Observation of water transport in the micro-porous layer of a polymer electrolyte fuel cell with a freezing method and cryo-scanning electron microscope

Micro-porous layers (MPLs) play an important role in the water management of polymer electrolyte fuel cells (PEFCs), however, the detailed mechanism of how the produced water is drained from these layers is not well understood. This paper observed the cross-sectional distribution of liquid water inside the cathode MPL to elucidate details of the phase state of the water transported through the MPL. The freezing method and cryo-scanning electron microscope (cryo-SEM) are used for the observations; the freezing method enables immobilization of the liquid water in the cell as ice forms by the freezing, and the cryo-SEM can visualize the water distribution in the vicinity of the MPL at high resolution without the ice melting. It was shown that no liquid water accumulates inside the MPL in operation at 35 °C, while the pores of the MPL are filled with liquid water under very low cell temperature operation, at 5 °C. These results indicate that the produced water passes through the MPL not as a liquid but in the vapor state in usual PEFC operation. Additionally, liquid water at the interface between the MPL and a catalyst layer (CL) was identified, and the effect of the interfacial contact on the water distribution was examined.     

Surface-induced shift of melting/freezing temperatures at interfaces between two semi-infinite media

Many physical properties and the structure of liquids at interfaces differ from those of the same liquids in the bulk state. In particular, the freezing temperatures of numerous liquid systems at interfaces between two fluids can be either increased or decreased compared to the temperature of corresponding bulk phase transition. The approach proposed in this paper allows one to predict the surface-induced shift of the melting/freezing temperatures in various systems. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 14–18, January, 2007.     

Glass transition and phase state of organic compounds: dependency on molecular properties and implications for secondary organic aerosols in the atmosphere

Recently, it has been proposed that organic aerosol particles in the atmosphere can exist in an amorphous semi-solid or solid (i.e. glassy) state. In this perspective, we analyse and discuss the formation and properties of amorphous semi-solids and glasses from organic liquids. Based on a systematic survey of a wide range of organic compounds, we present estimates for the glass forming properties of atmospheric secondary organic aerosol (SOA). In particular we investigate the dependence of the glass transition temperature T(g) upon various molecular properties such as the compounds' melting temperature, their molar mass, and their atomic oxygen-to-carbon ratios (O:C ratios). Also the effects of mixing different compounds and the effects of hygroscopic water uptake depending on ambient relative humidity are investigated. In addition to the effects of temperature, we suggest that molar mass and water content are much more important than the O:C ratio for characterizing whether an organic aerosol particle is in a liquid, semi-solid, or glassy state. Moreover, we show how the viscosity in liquid, semi-solid and glassy states affect the diffusivity of those molecules constituting the organic matrix as well as that of guest molecules such as water or oxidants, and we discuss the implications for atmospheric multi-phase processes. Finally, we assess the current state of knowledge and the level of scientific understanding, and we propose avenues for future studies to resolve existing uncertainties.     

Dynamics of transition-metal clusters

The atomic structure and thermodynamic properties of transition-metal 6- and 7-atom clusters are investigated using the molecular dynamics method, where Gupta's potential taking into account many-body interaction is employed. The caloric and the structural fluctuations are studied. The “fluctuating state” is found forN=6 in the region of the temperature near and below the melting point, where clusters undergo structural transition from one isomer to others without making any topological change. The fluctuating state differs from the “coexistence state” found in Ar clusters [1] i.e. the former involves no liquid state. In the liquid state the motion of atom-permutation occurs besides the breathing motion. On the other hand, the fluctuating state is not found forN=7 but only the motion of atom-permutation in the liquid state. The coexistence state is found in both cases of 6- and 7-atom clusters. We also discuss a possibility of larger clusters displaying the fluctuating state.     

Melting of monatomic glass with free surfaces

Melting of monatomic glass with free surfaces has been studied by molecular dynamics simulations in models with Lennard-Jones-Gauss interatomic potential. Models have been heated up from a glassy state toward a normal liquid state. Atomic mechanism of melting has been analyzed via monitoring spatio-temporal arrangements of liquid-like atoms occurred during heating process. Liquid-like atoms are detected via the Lindemann criterion of melting. It is clear that the transition from glass into supercooled liquid of our "ordinary" glass with free surfaces exhibits a non-heterogeneous behavior, i.e., although liquid-like atoms initiate/grow mainly in the surface shell, significant amount of liquid-like atoms also initiates/grows simultaneously in the interior during heating process. We found three characteristic temperatures of melting of glass with a free surface. Temperature dependence of structure and various thermodynamic quantities of the system upon heating is also presented and discussed.     

Accurate freezing and melting equations for the Lennard-Jones system

Analyzing three approximate methods to locate liquid-solid coexistence in simple systems, an observation is made that all of them predict the same functional dependence of the temperature on density at freezing and melting of the conventional Lennard-Jones (LJ) system. The emerging equations can be written as T=Aρ(4)+Bρ(2) in normalized units. We suggest to determine the values of the coefficients A at freezing and melting from the high-temperature limit, governed by the inverse 12th power repulsive potential. The coefficients B can be determined from the triple point parameters of the LJ fluid. This produces freezing and melting equations which are exact in the high-temperature limit and at the triple point and show remarkably good agreement with numerical simulation data in the intermediate region.     

Aligned/unaligned conducting polymer cryogels with three-dimensional macroporous architectures from ice-segregation-induced self-assembly of PEDOT-PSS

Porous conducting polymers are of great interest because of the huge potential to combine high surface areas in the dry state with physical properties relevant to organic electronics. Aligned or unaligned conducting polymer cryogels with 3D macroporous architectures have been prepared using the ice-segregation-induced self-assembly (ISISA) of different poly(3,4-ethylenedioxythiophene)-poly(styrenesulfonate) (PEDOT-PSS) freezing precursors as a dispersion or a formed hydrogel. The chemical composition and molecular structure of the resulting conducting polymer cryogels have been investigated by X-ray photoelectron spectroscopy and Raman spectroscopy, respectively. The morphologies of the PEDOT-PSS cryogels, together with their textural structures, have been revealed by scanning electron microscopy, mercury porosimetry, and nitrogen sorption tests. Processing PEDOT-PSS via ISISA endows the conducting polymers with novel properties, as demonstrated by a series of X-ray diffraction, differential scanning calorimetry, and electrical conductivity tests. These conducting polymer cryogels with aligned/unaligned macroporous architectures suggest the potential in the development of electronic components, tissue engineering, and next-generation catalytic and separation supports.