Molecular simulation of gas hydrate nanoclusters in water shell: Structure and phase transitions

Gas hydrate nanoclusters surrounded by water shells are studied by the molecular dynamic method. Hydrates of methane (sI structures) and krypton (sII structures), as well an ice nanocluster in a supercooled water shell, are considered. The main attention was focused on studying the local structure and phase transitions. Variations in local partial densities with an increase in temperature are monitored. Melting points of nanosized samples of gas hydrates are determined using caloric curves. Additional information on the behavior of the considered systems is obtained from the temperature dependences of diffusion coefficients and the Lindemann criterion. Two-phase transitions are revealed for gas hydrate nanoclusters. The first phase transition at 210 K can be assigned to the melting of the ice shell. The second transition at 230–235 K is identified as the phase transition in the hydrate core. The melting of ice cluster is observed at 215 K, which corresponds to the melting point of bulk crystal upon the use of the SPC/E water model.     

Ion calorimetry: Using mass spectrometry to measure melting points

Calorimetry measurements have been used to probe the melting of aluminum cluster cations with 63 to 83 atoms. Heat capacities were determined as a function of temperature (from 150 to 1050 K) for size-selected cluster ions using an approach based on multicollision-induced dissociation. The experimental method is described in detail and the assumptions are critically evaluated. Most of the aluminum clusters in the size range examined here show a distinct peak in their heat capacities that is attributed to a melting transition (the peak is due to the latent heat). The melting temperatures are below the bulk melting point and show enormous fluctuations as a function of cluster size. Some clusters (for example, n = 64, 68, and 69) do not show peaks in their heat capacities. This behavior is probably due to the clusters having a disordered solid-like phase, so that melting occurs without a latent heat.     

Melting of Water Molecule Clusters in a Strong Electric Field under the Conditions Modeling Arctic Stratosphere

The melting of (H₂O)₄₀ice microparticles under the conditions of subarctic stratosphere was numerically simulated using the Monte-Carlo method. The melting point of the microparticles was found to be 60 K lower than the melting temperature of bulk ice. The melting was detected by the behavior of the internal energy of microparticles, their heat capacity and electric susceptibility. The melting was accompanied by a qualitative change in the molecular orientation order in a cluster. Abrupt changes in the molecular arrangement in the cluster were not found. An electric field destroys the molecular orientation order in the cluster, and the clear-cut phase transition disappears. An electric field increases the rotational mobility of molecules.     

Neon scattering off Sodium clusters at finite temperatures

We present findings from computer simulations of collisions of neon atomic beams with Na20 atomic clusters at different internal temperatures. A functional form for the double differential cross section is determined, and no simple signature of a phase transition is seen, even though the clusters undergo a melting phase transition in the temperature range investigated (100 K–400 K). However, such experiments can be used effectively to measure the internal cluster temperature.     

Melting and glass transition for Ni clusters

The melting of NiN clusters (N = 29, 50-150) has been investigated by using molecular dynamics (MD) simulations with a quantum corrected Sutton-Chen (Q-SC) many-body potential. Surface melting for Ni147, direct melting for Ni79, and the glass transition for Ni29 have been found, and those melting points are equal to 540, 680, and 940 K, respectively. It shows that the melting temperatures are not only size-dependent but also a symmetrical structure effect; in the neighborhood of the clusters, the cluster with higher symmetry has a higher melting point. From the reciprocal slopes of the caloric curves, the specific heats are obtained as 4.1 kB per atom for the liquid and 3.1 kB per atom for the solid; these values are not influenced by the cluster size apart in the transition region. The calculated results also show that latent heat of fusion is the dominant effect on the melting temperatures (Tm), and the relationship between S and L is given.     

A Hundred‐Year‐Old Experiment Re‐evaluated: Accurate Ab Initio Monte Carlo Simulations of the Melting of Radon

State‐of‐the‐art relativistic coupled‐cluster theory is used to construct many‐body potentials for the noble‐gas element radon to determine its bulk properties including the solid‐to‐liquid phase transition from parallel tempering Monte Carlo simulations through either direct sampling of the bulk or from a finite cluster approach. The calculated melting temperature are 200(3) K and 200(6) K from bulk simulations and from extrapolation of finite cluster values, respectively. This is in excellent agreement with the often debated (but widely cited) and only available value of 202 K, dating back to measurements by Gray and Ramsay in 1909.     

Ab initio molecular dynamical investigation of the finite temperature behavior of the tetrahedral Au19 and Au20 clusters

Density functional molecular dynamics simulations have been carried out to understand the finite temperature behavior of Au19 and Au20 clusters. Au20 has been reported to be a unique molecule having tetrahedral geometry, a large HOMO-LUMO energy gap, and an atomic packing similar to that of the bulk gold (Li, J.; et al. Science 2003, 299, 864). Our results show that the geometry of Au19 is exactly identical with that of Au20 with one missing corner atom (called a vacancy). Surprisingly, our calculated heat capacities for this nearly identical pair of gold clusters exhibit dramatic differences. Au20 undergoes a clear and distinct solid-like to liquid-like transition with a sharp peak in the heat capacity curve around 770 K. On the other hand, Au19 has a broad and flat heat capacity curve with continuous melting transition. This continuous melting transition turns out to be a consequence of a process involving a series of atomic rearrangements along the surface to fill in the missing corner atom. This results in a restricted diffusive motion of atoms along the surface of Au19 between 650 to 900 K during which the shape of the ground state geometry is retained. In contrast, the tetrahedral structure of Au20 is destroyed around 800 K, and the cluster is clearly in a liquid-like state above 1000 K. Thus, this work clearly demonstrates that (i) the gold clusters exhibit size sensitive variations in the heat capacity curves and (ii) the broad and continuous melting transition in a cluster, a feature that has so far been attributed to the disorder or absence of symmetry in the system, can also be a consequence of a defect (absence of a cap atom) in the structure.     

Thermodynamic characteristics and the temperatures of relaxation transitions of poly(methacrylic acid)

The heat capacity of poly(methacrylic acid) containing 2.5 wt % water was measured in a vacuum adiabatic calorimeter at temperatures between 80 and 325 K. The heat capacity of anhydrous poly(methacrylic acid) was calculated, and its standard enthalpies of combustion and formation were determined. On the basis of the enthalpy of melting of the “free”-water phase, the limit of water solubility in the polymer was found calorimetrically at 273 K. The temperatures of relaxation transitions (the glass transition and the β and γ transitions) of poly(methacrylic acid) mixtures with water were determined via differential thermal analysis in the region 80–550 K. In addition, the determination of the temperatures of transitions of anhydrous poly(methacrylic acid) was performed via extrapolation to zero water content of the concentration dependences of the relaxation-transition temperatures.     

Phase transitions in mesophase macromolecules. III. The transitions in poly(ethylene terephthalate-co-p-oxybenzoate)

The glass and melting transitions of poly(ethylene terephthalate-co-p-oxybenzoate)s have been studied by differential scanning calorimetry. Despite the higher glass transition expected for polyoxybenzoate, there is almost no change in the glass transition temperature up to 63 mol % oxybenzoate (353 ± 4 K). At high oxybenzoate concentration there seems to be a discontinuous jump to a glass transition of 450 K. This high glass transition has been observed in two-phase materials down to 30 mol % oxybenzoate. The melting transition shows signs of isodimorphism with a minimum in melting temperature at about 60–70 mol % oxybenzoate, 60 K below the melting temperature of poly(ethylene terephthalate). The thermal properties are little affected by the change of the noncrystalline parts of the molecules to a mesophase structure. The transition to the isotropic phase could not be analyzed because of prior decomposition.     

Orientational melting and reorientational motion in a cubane molecular crystal: a molecular simulation study

Detailed molecular simulations are carried out to investigate the effect of temperature on orientational order in cubane molecular crystal. We report a transition from an orientationally ordered to an orientationally disordered plastic crystalline phase in the temperature range 425-450 K. This is similar to the experimentally reported transition at 395 K. The nature of this transition is first order and is associated with a 4.8% increase in unit cell volume that is comparable to the experimentally reported unit cell volume change of 5.4% (Phys. Rev. Lett. 1997, 78, 4938). An orientational order parameter, eta(T), has been defined in terms of average angle of libration of a molecular 3-fold axis and the orientational melting has been characterized by using eta(T). The orientational melting is associated with an anomaly in specific heat at constant pressure (C(P)) and compressibility (kappa). The enthalpy of transition and entropy of transition associated with this orientational melting are 20.8 J mol(-1) and 0.046 J mol(-1) K(-1), respectively. The structure of crystalline as well as plastic crystalline phases is characterized by using various radial distribution functions and orientational distribution functions. The coefficient of thermal expansion of the plastic crystalline phase is more than twice that of the crystalline phase.     

Calorimetric study of the relationship between molecular structure and liquid crystallinity of rod-like mesogens II. Heat capacities and phase transitions of 4-(trans-4-propylcyclohexyl)benzonitrile

The molar heat capacity of the rod-like compound 4-(trans-4-propylcyclohexyl)benzonitrile (3-CBCN), purity of 99.8mol%, has been measured with an adiabatic calorimeter at temperatures between 15 and 385K. 3-CBCN is a nematogenic mesogen, whose melting and clearing points are 316.33 and 319.09 K, respectively. The enthalpy and entropy gained at fusion are 20.4 kJmol -1 and 64.4 J K -1 mol -1, respectively; those for the nematic-to-isotropic transition are 1.1 kJmol -1 and 3.5 J K -1 mol -1 respectively. 3-CBCN exhibits a supercooled nematic phase, whose molar heat capacities have been measured from 25 K below the melting point. The molar and transition entropies of 3-CBCN are discussed in relation to those of 4-propylbiphenyl-4-carbonitrile (3-BBCN) and trans,trans-4'-propylbicyclohexyl- 4-carbonitrile (3-CCCN). There seems to exist a correlation between these values and mesophase stability. Finally, Eidenschink's theoretical model for the nematic-to-isotropic transition has been applied to 3-CBCN; the transition enthalpy estimated according to this model agrees well with the observed value.     

The properties of ion-water clusters. II. Solvation structures of Na+, Cl-, and H+ clusters as a function of temperature

Ion-water-cluster properties are investigated both through the multistate empirical valence bond potential and a polarizable model. Equilibrium properties of the ion-water clusters H+(H2O)100, Na+(H2O)100, Na+(H2O)20, and Cl-(H2O)17 in the temperature region 100-450 K are explored using a hybrid parallel basin-hopping and tempering algorithm. The effect of the solid-liquid phase transition in both caloric curves and structural distribution functions is investigated. It is found that sodium and chloride ions largely reside on the surface of water clusters below the cluster melting temperature but are solvated into the interior of the cluster above the melting temperature, while the solvated proton was found to have significant propensity to reside on or near the surface in both the liquid- and solid-state clusters.     

Crystallization and melting of a branched polyethylene with precisely controlled chemical structure

The heat capacity of a linear polyethylene with dimethyl branches, at every 21st backbone atom was analyzed by differential scanning calorimetry (DSC) and quasi-isothermal temperature-modulated DSC. This novel copolyethylene (PE2M) is relatively difficult to crystallize from the melt. On subsequent heating, a first, sharp melting peak is followed by a sharp cold-crystallization and crystal perfection and a smaller endotherm, before reaching the main melting at 315–320 K, close to the melting temperatures of eicosane and tetracontane. The low-temperature melting is sensitive to the cooling rate and disappears below 1.0 K min⁻¹. The cold crystallization can be avoided by heating with rates faster than 80 K min⁻¹. The PE2M exhibits some reversing and reversible melting, which is typical for chain-folded polymers. The glass transition of semicrystalline PE2M is broadened and reaches its upper limit at about 260 K (midpoint at about 0.355 K). Above this temperature, the crystals seem to have a heat capacity similar to that of the liquid. A hypothesis is that the melting transition can be explained by changes in crystal perfection without major alteration of the crystal structure and the lamellar morphology. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3461–3474, 2006     

Isomeric transitions in size-selected methanol hexamers probed by OH-stretch spectroscopy

We have measured the isomeric transition between the energetically lowest lying isomers of S6 and C2-symmetry of (CH3OH)6. The clusters are size-selected by deflection in collisions with He, and the isomers are identified by their infrared spectra of the OH-stretching vibration. The measurements are carried out at three source temperatures 253, 300 and 373 K which correspond to the cluster temperatures 93, 106 and 135 K. The latter ones are estimated by a relaxation model that accounts for the cluster formation and the energy released by the condensation. The transition takes place at a cluster temperature of about 102 K which is in agreement with the Molecular Dynamics simulation of such a transition at about 117 K using a realistic model potential.     

三甲醇丙烷—季戊四醇固体溶液的热容和相变焓

前曾报道季戊四醇在461.60K有一个固-固相变,其相交含为41.37kJmol~(-1)。但是,作为低温储能材料,该物质的相变温度偏高。从有关固-固相变储能材料的热力学研究中,我们发现三甲醇丙烷-季戊四醇固体溶液具有较低的固-固相变温度。本文测定了三甲醇丙烷-季戊四醇(摩尔比60:40)固体溶液的相变热参数和热容。     

Transitions in poly(diethyl siloxane)

An adiabatic heat capacity study of poly(diethylsiloxane) confirms that it has a single glass transition occurring at 130°K, the lowest glass transition reported to date for a high molecular weight polymer. The two previously reported glass transitions are first-order thermodynamic peaks whose location is dependent upon prior thermal history. Combination of these data with low-temperature x-ray diffraction indicates that the transitions in this temperature range are related to a crystal–crystal transformation. A crystal melting transition is observed near 270°K. In addition an anomalous rise in heat capacity near 60°K suggests a sub-glass transition of unknown origin.     

New Structure Formation on γ-irradiation of Poly(chlorotrifluoroethylene)

The radiation chemistry of PCTFE at different temperatures has been studied. The polymer was irradiated under vacuum to absorbed doses of up to 1500 kGy. Three irradiation temperatures were chosen. These included ambient temperature, a temperature well above the T_g and a temperature above the crystalline melting temperature. These were 298, 423 and 493 K, respectively. The formation of new structures was identified by solid-state FTIR and ¹⁹F NMR. No branching was observed below the melting temperature, but branches were observed above the melting temperature. G-values for chain-end formation were 1.5 and 2.4 at room temperature and 423 K, respectively and the G-value for the formation of double bonds was found to be <0.1. For the irradiations at 493 K, the G-values for the formation of chain ends, double bonds and branching points were 3.6, 0.2 and 0.5, respectively. The presence of long-chain branches within the polymer structure could not be proven for radiolysis at 493 K, but scission predominates and network formation does not occur upon irradiation. DSC studies of the polymers irradiated at ambient temperature were consistent with chain scission leading to an increase in the percentage crystallinity, as observed for other fluoropolymers.     

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.     

Proton transfer-induced conformational changes and melting in designed peptides in the gas phase

The conformations of protonated RA15K, RA20K and RA15H (R = arginine, A = alanine, K = lysine, and H = histidine) have been examined in the gas phase as a function of temperature. These peptides were designed so that intramolecular proton transfer will trigger conformational changes between a helix (proton sequestered at the C-terminus) and globule (proton sequestered at the N-terminus). Kinetically controlled structural transitions occur below 400 K (from helix to globule for RA15H, and from globule to helix for RA15K and RA20K). As the temperature is raised, the compact globule found at room temperature expands, accesses more configurations, and becomes entropically favored. At around 500 K, the RA15K and RA20K helices undergo a melting transition. The transition is broad, as expected for a phase transition in a finite system, and becomes narrower as the peptide size increases. In the helical conformation, the two basic residues are well separated; as a result, the proton transfer necessary to drive the melting transition probably involves a mobile proton. For doubly protonated RA15K, a dumbbell-like conformation (resulting from repulsion between the two protonated basic residues) is found at high temperature.     

Specific heat capacity and melting/crystallization characteristics of polytridecanolactone

An investigation of the thermodynamical properties of polytridecanolactone (PTDL) was made with the aid of a differential scanning calorimeter (DSC). PTDL is a linear polyester and belongs to the polylactones, which have been poorly investigated. In this paper we contribute with specific heat capacity in the range 180-400 K, and melting and glass transition characteristics. Further, we present unique results corresponding to the effect of different cooling rates on crystallization temperatures and crystallization energies. PTDL has a melting temperature of 350 K, and a glass transition at about 237 K. The crystallization results show that PTDL crystallizes easily, with a crystallization degree of about 80%. In addition, the crystallization energy decreases with increasing cooling rate, and levels out at a constant value at higher cooling rates. The crystallization temperature, on the other hand, shows an increasing sensitivity of cooling rate, where the supercooling is increasing more rapidly at higher cooling rates. © 1994 John Wiley & Sons, Inc.