The potential energy landscape contribution to the dynamic heat capacity

The dynamic heat capacity of a simple polymeric, model glassformer was computed using molecular dynamics simulations by sinusoidally driving the temperature and recording the resultant energy. The underlying potential energy landscape of the system was probed by taking a time series of particle positions and quenching them. The resulting dynamic heat capacity demonstrates that the long time relaxation is the direct result of dynamics resulting from the potential energy landscape. Moreover, the equilibrium (low frequency) portion of the potential energy landscape contribution to the heat capacity is found to increase rapidly at low temperatures and at high packing fractions. This increase in the heat capacity is explained by a statistical mechanical model based on the distribution of minima in the potential energy landscape.     

Dynamical behavior of one-dimensional water molecule chains in zeolites: nanosecond time-scale molecular dynamics simulations of bikitaite

Nanosecond scale molecular dynamics simulations of the behavior of the one-dimensional water molecule chains adsorbed in the parallel nanochannels of bikitaite, a rare lithium containing zeolite, were performed at different temperatures and for the fully and partially hydrated material. New empirical potential functions have been developed for representing lithium-water interactions. The structure and the vibrational spectrum of bikitaite were in agreement both with experimental data and Car-Parrinello molecular dynamics results. Classical molecular dynamics simulations were extended to the nanosecond time scale in order to study the flip motion of water molecules around the hydrogen bonds connecting adjacent molecules in the chains, which has been observed by NMR experiments, and the dehydration mechanism at high temperature. Computed relaxation times of the flip motion follow the Arrhenius behavior found experimentally, but the activation energy of the simulated system is slightly underestimated. Based on the results of the simulations, it may be suggested that the dehydration proceeds by a defect-driven stepwise diffusion. The diffusive mechanism appears as a single-file motion: the molecules never pass one another, even at temperatures as high as about 1000 K, nor can they switch between different channels. However, the mean square displacement (MSD) of the molecules, computed with respect to the center of mass of the simulated system, shows an irregular trend from which the single-file diffusion cannot be clearly evidenced. If the MSDs are evaluated with respect to the center of mass of the molecules hosted in each channel, the expected dependence on the square root of time finally appears.     

四氯合镉酸正十一烷铵的合成、晶体结构及低温热容

合成了四氯合镉酸正十一烷铵配合物(C₁₁H₂₃NH₃)₂CdCl₄(s)[简写: C₁₁Cd(s)]. 用X 射线单晶衍射技术、化学分析和元素分析确定其晶体结构和化学组成. 利用其晶体学数据计算出晶格能为: U_POT＝908.18 kJ·mol^-1. 利用精密自动绝热热量计测定了它在78～395 K 温区的低温热容, 结果表明, 该配合物在此温区出现两次连续的固-固相转变, 计算出两次相变的峰温、摩尔焓及摩尔熵分别为: T_trs,1＝(321.88±0.07) K, Δ_trsH_m,1＝(37.59±0.17) kJ·mol^-1, Δ_trsS_m,1＝(117.24±0.12) J·K^-1·mol^-1, T_trs,2＝(323.81±0.30) K, Δ_trsH_m,2＝(12.42±0.02) kJ·mol^-1 和Δ_trsS_m,2＝(38.36±0.09) J·K^-1·mol^-1. 用最小二乘法将实验摩尔热容对温度进行拟合, 得到热容随温度变化的多项式方程. 用此方程进行数值积分,得到此温区每隔5 K 的舒平热容值和相对于298.15 K 时的热力学函数值.     

Reversible and irreversible heat capacity of poly(trimethylene terephthalate) analyzed by temperature‐modulated differential scanning calorimetry

The heat capacity of poly(trimethylene terephthalate) (PTT) has been analyzed using temperature‐modulated differential scanning calorimetry (TMDSC) and compared with results obtained earlier from adiabatic calorimetry and standard differential scanning calorimetry (DSC). Using quasi‐isothermal TMDSC, the apparent reversing and nonreversing heat capacities were determined from 220 to 540 K, including glass and melting transitions. Truly reversible and time‐dependent irreversible heat effects were separated. The extrapolated vibrational heat capacity of the solid and the total heat capacity of the liquid served as baselines for the analysis. As one approaches the melting region from lower temperature, semicrystalline PTT shows a reversing heat capacity, which is larger than that of the liquid, an observation that is common also for other polymers. This higher heat capacity is interpreted as a reversible surface or bulk melting and crystallization, which does not need to undergo molecular nucleation. Additional time‐dependent, reversing contributions, dominating at temperatures even closer to the melting peak, are linked to reorganization and recrystallization (annealing), while the major melting is fully irreversible (nonreversing contribution). © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 622–631, 2000     

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.     

Quasielastic neutron scattering of poly(methyl phenyl siloxane) in the bulk and under severe confinement

Quasielastic neutron scattering was utilized to investigate the influence of confinement on polymer dynamics. Poly(methyl phenyl siloxane) chains were studied in the bulk as well as severely confined within the approximately 1-2 nm interlayer spacing of intercalated polymer/layered organosilicate nanohybrids. The temperature dependence of the energy resolved elastic scattering measurements for the homopolymer and the nanocomposites exhibit two distinct relaxation steps: one due to the methyl group rotation and one that corresponds to the phenyl ring flip and the segmental motion. Quasielastic incoherent measurements show that the very local process of methyl rotation is insensitive to the polymer glass transition temperature and exhibits a wave-vector independent relaxation time and a low activation energy, whereas it is not affected at all by the confinement. At temperatures just above the calorimetric glass transition temperature, the observed motion is the phenyl ring motion, whereas the segmental motion is clearly identified for temperatures about 60 K higher than the glass transition temperature. For the nanohybrid, the segmental motion is found to be strongly coupled to the motion of the surfactant chains for temperatures above the calorimetric glass transition temperature of the bulk polymer. However, the mean square displacement data show that the segmental motion in confinement is faster than that of the bulk polymer even after the contribution of the surfactant chains is taken into consideration.     

Heat capacity and physical transitions in collagen and solubility of water in it

The heat capacities of cattle and poultry collagen samples were measured by adiabatic calorimetry in the range 80–330 K. The temperatures of physical transitions in these samples were determined. Endothermic relaxation transitions (γ₁, γ₂, and β) were observed in the heat capacity curves of the studied collagen samples. The solubility of water in collagen samples obtained from various sources was determined by calorimetry from the enthalpy of melting of the “free” water phase.     

Computer simulation study of thermodynamic scaling of dynamics of 2Ca(NO3)2·3KNO3

Molecular dynamics (MD) simulations of the glass-former 2Ca(NO(3))(2·3KNO(3), CKN, were performed as a function of temperature at pressures 0.1 MPa, 0.5 GPa, 1.0 GPa, and 2.0 GPa. Diffusion coefficient, relaxation time of the intermediate scattering function, and anion reorientational time were obtained as a function of temperature and densitiy ρ. These dynamical properties of CKN scale as ρ(γ)∕T with a common value γ = 1.8 ± 0.1. The scaling parameter γ is consistent with the exponent of the repulsive part of an effective intermolecular potential for the repulsion between the atoms at shortest distance in the equilibrium structure of liquid CKN, Ca(2+), and oxygen atoms of NO(3)(-). Correlation between potential energy and virial is obeyed for the short-range terms of the potential function, but not for the whole potential including coulombic interactions. Decoupling of diffusion coefficient and reorientational relaxation time from relaxation time take place at a given ρ(γ)∕T value, i.e., breakdown of Stokes-Einstein and Debye-Stokes-Einstein equations result from combined thermal and volume effects. The MD results agree with correlations proposed between long-time relaxation and short-time dynamics, lnτ ∝ 1∕, where the mean square displacement concerns a time window of 10.0 ps. It has been found that scales as ρ(γ)∕T above and below the glass transition temperature, so that thermodynamic scaling of liquid dynamics can be thought as a consequence of theories relating short- and long-time dynamics, and the more fundamental scaling concerns short-time dynamical properties.     

金属铜升温熔化过程的分子动力学模拟

采用分子动力学方法模拟了金属铜的升温熔化过程.原子间作用势采用FS (Finnis-Sinclair)势，结构分析采用双体分布函数(PCF)、均方位移(MSD)等方法.计算结果表明,在连续升温过程中，金属铜在1444 K熔化，在该熔化点的扩散系数为4.31×10－9 m2•s－1.上述结论与实验值相当接近，并且比之采用EAM镶嵌原子势所作模拟得到的结果更佳，说明FS势可以用来处理象液铜这样较复杂的无序体系.本文指出了升温速率在金属熔化过程中所起的作用.     

Simulations of Metal Cu in Heating Process

Introduction In order to study the short-range order as well as thermodynamic properties, two distinct techniques of computer simulation, namely, the molecular dynamics and Monte Carlo methods, are most frequently em-ployed. In both techniques, the interaction potential is the primary input for computation. Mitra and co-workers1 have used a two-body model with Coulombic interac-tions and a power-law repulsion, fitted to the short- range structure and melting temperature of cristobalite. Thre…     

Melting behavior of water in cylindrical pores: carbon nanotubes and silica glasses

We report a study of the effects of confinement in multi-walled carbon nanotubes and mesoporous silica glasses (SBA-15) on the solid structure and melting of both H(2)O and D(2)O ice, using differential scanning calorimetry, dielectric relaxation spectroscopy, and neutron diffraction. Multi-walled nanotubes of 2.4, 3.9 and 10 nm are studied, and the SBA-15 studied has pores of mean diameter 3.9 nm; temperatures ranging from approximately 110 to 290 K were studied. We find that the melting point is depressed relative to the bulk water for all systems studied, with the depression being greater in the case of the silica mesopores. These results are shown to be consistent with molecular simulation studies of freezing in silica and carbon materials. The neutron diffraction data show that the cubic phase of ice is stabilized by the confinement in carbon nanotubes, as well as in silica mesopores, and persists up to temperatures of about 240 K, above which there is a transition to the hexagonal ice structure.     

Modeling the molecular dynamics of liquid metals at high pressures: Liquid potassium

A way of building potassium models by molecular dynamics using the embedded atom model (EAM) is developed. The contribution from pairwise interaction is presented as power series at an interpartial distance. Embedded potential parameters are determined by the experimental dependence of pressure on volume for a static compression of potassium at 300 K to a pressure of 53 GPa (potential A). By using potential A to describe shock compression and choosing the appropriate temperature at given degree of compression (up to 40000 K at compression to 0.29 of initial volume) it is shown that the model pressure can be made equal to the pressure indicated by the Rankine-Hugoniot relations. The model energy is lower than the actual energy determined by the relations, and the difference in energies increases with temperature almost linearly; such growth corresponds to an excess average heat capacity of about 11.6 J/(mol K), compared to the model heat capacity. It is established that the reasons for this divergence are the inability of the EAM potential to describe the temperature dependency of metal properties precisely, and the appearance of an energy contribution upon heating that is dependent on temperature but not on atom coordinates. Adding another summand to the potential energy (which is dependent on temperature only) allows us to match the heat capacities of real potassium and the models. The dependence of potassium’s melting temperature on pressure is calculated. The calculated melting temperature at 41.2 GPa is 1231 K. Additional data (e.g., the actual temperature on the Rankine-Hugoniot curve and precise quantum mechanics calculations of heat capacity at extreme conditions) is required to eliminate potential ambiguity.     

The use of the embedded atom model for liquid metals: Liquid potassium

The method for calculations the embedded atom potential for liquid metals based on the diffraction data on the structure close to the melting temperature was applied to potassium. The embedded atom potential parameters were adjusted using the data on the structure of potassium at 343, 473, and 723 K and the thermodynamic properties of potassium at temperatures up to 37240 K. The use of the molecular dynamics method and the embedded atom potential gave close agreement with the experimental data on the structure, density, and potential energy of liquid metal along the p ? 0 isobar at temperatures up to 2200 K and along the shock adiabat up to a pressure of ～85 GPa and 37240 K. The calculated bulk compression modulus at 343 K was close to its actual value, and the self-diffusion coefficients increased under isobaric heating conditions following a power law with an exponent of 1.6478. The melting temperature of body-centered potassium with the embedded atom potential was (319 ± 1) K, which was close to the actual melting temperature. The potential obtained incorrectly described crystalline potassium.     

Heat capacity and glass transition of ethylene oxide clathrate hydrate

The heat capacity of structure I ethylene oxide clathrate hydrate EO-6.86 H₂O was measured in the temperature range 6–300 K with an adiabatic calorimeter. The temperature and enthalpy of congruent melting were determined to be (284.11 ± 0.02) K and 48.26 kJ mol^–1, respectively. A glass transition related to the proton configurational mode in the hydrogen-bonded host was observed around 90 K. This glass transition was similar to the one observed previously for the structure II tetrahydrofuran hydrate but showed a wider distribution of relaxation times. The anomalous heat capacity and activation enthalpy associated with the glass transition were almost the same as those for THF-hydrate.Dedicated to Dr D. W. Davidson in honor of his great contributions to the sciences of inclusion phenomena.Author for correspondence.     

Dielectric Relaxation Studies of Binary Mixture of α‐Picoline and Methanol Using Time Domain Reflectometry at Different Temperatures

Complex permittivity spectra of binary mixtures of varying concentrations of α‐picoline and methanol (MeOH) were obtained using time domain reflectometry (TDR) technique over frequency range 10 MHz to 25 GHz at 283.15 K, 288.15 K, 293.15 K and 298.15 K temperatures. The dielectric relaxation parame‐ ters namely static permittivity (σ₀), high frequency limit permittivity (σ_oo1) and the relaxation time (ρ) were determined by fitting complex permittivity data to the single Debye/Cole‐Davidson model. Complex non linear least square (CNLS) fitting procedure was carried out using LEVMW software. The excess static permittivity (σ₀^E) and the excess inverse relaxation time (1/ρ)^E which contains information regarding mo‐ lecular structure and interaction between polar — polar liquids, were also determined. From the experimental data, effective Kirkwood correlation factor (g^eff) and Bruggeman factor (f_B) were calculated. Excess parameters were fitted to the Redlich‐Kister polynomial equation. The values of static permittivity and relaxation time increase non‐linearly with increase in the mol fraction of MeOH at all temperatures. The values of excess static permittivity (σ₀^E) and the excess inverse relaxation time (1/ρ)^E are negative for the studied α‐picoline — MeOH system at all temperatures.     

Quasi-isothermal measurement of frequency dependent heat capacity of semicrystalline polyethylene at the melting temperature using light heating modulated temperature DSC

Frequency dependence of the heat capacity of semicrystalline polyethylene was measured using light heating modulated temperature DSC (LMDSC). The quasi-isothermal measurement was carried out in the melting temperature range of the crystal. Two decades of the frequency, from 0.01 to 1 Hz, was covered by the LMDSC instrument constructed in the author's laboratory. It was found that in the melting temperature range polyethylene exhibited Debye relaxation with the relaxation time of 14 s and had excess heat capacity independent of the kinetics. The excess heat capacity and the relaxation strength could be attributed to the disordered crystals generated during rapid cooling from the molten state and/or its surrounding amorphous region instead of the stable crystals reorganized during the quasi-isothermal measurement.     

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     

Ultrafast laser studies of the photothermal properties of gold nanocages

Au nanocages were synthesized via a galvanic replacement reaction. The extinction peak of these hollow structured particles is shifted into the near-IR compared with the Ag nanocube templates. Energy transfer from the Au nanocages into the surrounding environment (water) as well as the coherently excited vibrational modes of the nanocages were studied by femtosecond pump-probe spectroscopy. The time scale for energy relaxation was found to increase with the size of the particles, with the relaxation time being independent of the laser intensity. The time scales for relaxation are comparable to those for solid spherical gold particles and are consistent with energy relaxation being controlled by heat dissipation in the solvent. The period of the coherently excited vibrational mode is proportional to the dimensions of the nanocages. Intensity-dependent measurements show that in solution the nanocages maintain their integrity up to lattice temperatures of 1100 +/- 100 K.     

How anesthetics, neurotransmitters, and antibiotics influence the relaxation processes in lipid membranes

We used pressure perturbation calorimetry to investigate the relaxation time scale after a jump into the melting transition regime of artificial lipid membranes. This time is equivalent to the characteristic rate of domain growth. The studies were performed on single-component large unilamellar and multilamellar vesicle systems with and without the addition of small molecules such as general anesthetics, neurotransmitters, and antibiotics. These drugs interact with membranes and affect melting points and profiles. In all systems, we found that heat capacity and relaxation times are linearly related to each other in a simple manner, and we outline the theoretical origin of this finding. Thus, the influence of a drug on the time scale of domain formation processes can be understood on the basis of their influence on the heat capacity profile. This allows estimations of the characteristic relaxation time scales in biological membranes.