Premelting Transition Research Articles

1/f2 noise as a precursor of structural reconstructions near the melting point of crystalline materials with different types of chemical bonds

Transition phenomena near the melting point (premelting effect) are fundamental processes occurring in materials with different types of chemical bonding. At T ≥ 0.8Tm, there is a fluctuating heat release. The analysis of fluctuation processes allows obtaining information about dynamic reconstructions in various subsystems and interconnections between them within the system. The purpose of this work was to study the spectral characteristics of heat fluctuations in stationary modes of premelting materials with ionic, covalent, and metallic chemical bonds (KCl, Ge, Cu, and Sb) and to determine the type of the fluctuation process using the Hurst parameter. The spectral characteristics of heat fluctuations in stationary modes of pre-melting KCl, Ge, Cu, and Sb at T* ~ 0.9Tm were determined by wavelet analysis. This method allows analyzing the behavior of complex systems at critical points in order to identify certain correlations and development trends in them. The study showed that in the premelting region, the frequency spectrum of heat fluctuations was characteristic of 1/f2 noise or nonlinear Brownian noise, which is a precursor of structural reconstructions during phase transitions. The type of fluctuation processes in the premelting region of KCl, Ge, Cu, and Sb was determined using the Hurst parameter (H). It was shown that in stationary modes of premelting H > 0.5. Consequently, the previous trend of the process dynamics was very likely to develop in the same direction. However, with a decrease in the energy of the chemical bond H Æ 0.5, which indicated a decrease in the stability of the system and a likely change in the development trend for the structural reconstructions inthe premelting transition region. Thus, near the melting point, there are unstable dynamic states, which are precursors to structural changes in the system, which have certain developmental trends. This should be taken into account when calculating the stability and reliability of materials and systems

Disordering complexion transition of grain boundaries in bcc metals: Insights from atomistic simulations

Disordering complexion transition of grain boundaries in bcc metals: Insights from atomistic simulations

Prediction of superionic state in LiH2 at conditions enroute to nuclear fusion

Hydrogen and lithium, along with their compounds, are crucial materials for nuclear fusion research. High-pressure studies have revealed intricate structural transitions in all these materials. However, research on lithium hydrides beyond LiH has mostly focused on the low-temperature regime. Here, we use density functional theory and ab initio molecular dynamics simulations to investigate the behavior of LiH2, a hydrogen-rich compound, near its melting point. Our study is particularly relevant to the low-pressure region of the compression pathway of lithium hydrides toward fusion. We discovered a premelting superionic phase transition in LiH2 that has significant implications for its mass transportation, elastic properties, and sound velocity. The theoretical boundary for the superionic transition and melting temperature was then determined. In contrast, we also found that the primary compound of lithium hydrides, LiH, does not exhibit a superionic transition. These findings have important implications for optimizing the compression path to achieve the ignition condition in inertial confinement fusion research, especially when lithium tritium-deuteride (LiTD) is used as the fuel.

The Heat Capacity of PuO2 at High Temperature: Molecular Dynamics Calculations

Abstract A new generation of fast breeder reactors (FBRs) is under development with the objective of making nuclear energy more sustainable. Most promising reactor designs are loaded, at least during their early phase of deployment, with UO2–PuO2 mixed oxide fuel (MOX). Concentrations of plutonium dioxide that are foreseen for FBRs range up to 30 mol%. This highlights the need for a sound and deep knowledge of the thermophysical properties of PuO2. This statement is valid in the case of heat capacity, as evaluations on MOX fuel are usually carried out by using the Neumann–Kopp rule. Heat capacity is relevant for thermal conductivity and performance under transient conditions. However, measurements on the heat capacity of plutonium dioxide are scarce or even lacking at high temperature. Numerical methodologies such as molecular dynamics (MD) calculations have been employed to overcome the difficulties encountered in experimental measurements. Besides numerical also theoretical models have been applied as valuable tools for interpretation of enthalpy measurements. Nevertheless, due to the mentioned lack of experimental measurement issues such as the existence of the Bredig transition and the formation of defects at high temperatures are still debated in nuclear fuel research. Excess enthalpy seen in measurements of actinides oxides has been explained by means of either electronic disorder or anion disorder. In the case of plutonium dioxide, a common consensus has been reached on the hypothesis that anion disorder leads to a significant increase in heat capacity at high temperature. Konings and Beneš have developed a model that accounts for this phenomenon. Their correlation has been often included in models of heat capacity and employed for recommendations. However, in the high-temperature region, MD calculations showed an underestimation of model predictions that was not compensated by the presence of a peak of heat capacity that has been interpreted as the Bredig transition. Based on these observations, this paper presents MD evaluations on the heat capacity of PuO2 at high temperature that are mostly focused on the formation energy of oxygen Frenkel pairs (OFPs) and its correlation with the model proposed by Konings and Beneš. Besides an interatomic potential published in the open literature and developed in compliance with the experimental thermal expansion of PuO2, a second interatomic potential has been applied in calculations. This latter is featured by a lower formation energy of OFP. The contribution due to defects formation was calculated by means of a simplified theoretical model of heat capacity. Results of calculations in the very high-temperature domain showed an increase in the contribution due to OFP defects consistent with the model by Konings and Beneš. Predictions suggest the onset of a premelting transition around 85% of melting temperature without the presence of a peak of heat capacity. Major deviations from the recommended model have been noted in the intermediate temperature region where the effect of clustering of defects should play a significant role. Therefore, the value of formation energy of OFP proposed by Konings and Beneš could be interpreted as an effective value that accounts for the two processes (defects clustering and premelting transition) that could contribute, according to our results, to the heat capacity of plutonium dioxide at high temperature. This conclusion is consistent with the numerical evaluations of OFP formation energy that are in general higher than proposed by Konings and Beneš.

Emergence of dynamical disorder and phase metastability in carbon nanobowls

The condensed phases of the carbon-based nanobowl corannulene are herein investigated in the temperature range 200–573 K, combining differential scanning calorimetry, synchrotron X-ray diffraction, and quasielastic neutron-scattering. For the first time, we identify the presence of a well-defined thermal event at 382 K, a figure well below a melting point of 540 K. Contrary to naïve expectation, detailed analysis of the neutron-scattering data below and above this pre-melting transition signals the emergence of correlated stochastic dynamics within both thermodynamically stable solid phases of the material. We find that these are progressively responsible for the suppression of molecular and supramolecular order over mesoscopic length scales, and are associated with the formation of high-symmetry rotor-like states exhibiting localized stochastic motions. Upon cooling from the melt, we have also discovered a robust hysteresis associated with the existence of hitherto-unknown metastable liquid (deep-supercooled) and disordered-solid phases. This behaviour is markedly different from that observed in the quintessential Buckminsterfullerene C 60 or other chemically substituted fullerene adducts studied to date at this level of detail. These results evince new and yet-to-tapped opportunities for the use of the stable and metastable phases of carbon-based nanobowls in novel applications exploiting the emergence of dynamical disorder at the nanoscale.

Base-specific pre-melting and melting transitions of DNA in presence of ionic liquids probed by synchrotron-based UV resonance Raman scattering

Hydrated ionic liquids (ILs) have been identified as solvent media able to enhance the structural stability of deoxyribonucleic acid (DNA). In this work, we investigate the molecular interaction between imidazolium-based ILs and DNA during its thermal unfolding pathway, by exploiting synchrotron-based UV Resonance Raman scattering (UVRR) experiments. This technique gives a selective focus on the thermal responses of specific nucleobases in the structure of DNA, providing the experimental sensitivity to both cooperative and local structural changes occurring along the complex unfolding process of DNA. UVRR measurements probe two distinct temperature-dependent phenomena occurring in the DNA double-helix, i.e. a non-cooperative pre-melting mainly involving adenine bases and a cooperative melting transition primarily localized on guanine tracts. The analysis of Raman spectra reveals that both the cation and anion of the ionic liquids strongly interact with the structure of DNA, thus affecting the melting process but not perturbing the pre-melting transition that precedes the complete separation of the strands of DNA. Overall these results suggest that the dominant interaction occurs between the imidazolium cation and the bases of guanine and thymine in the structure of DNA, in agreement with previous results of molecular dynamics simulations. • Thermal stability of DNA in water/ionic liquids is explored by UV Resonance Raman • Thermal responses of specific nucleobases are selectively detected • Melting and pre-melting are specifically probed using proper excitation wavelength • Both cation and anion of ionic liquids strongly interact with DNA structure • The dominant interaction with ionic liquid involves guanine and thymine

A quantitative model of a cooperative two-state equilibrium in DNA: experimental tests, insights, and predictions - CORRIGENDUM.

Quantitative parameters for a two-state cooperative transition in duplex DNAs were finally obtained during the last 5 years. After a brief discussion of observations pertaining to the existence of the two-state equilibrium per se, the lengths, torsion, and bending elastic constants of the two states involved and the cooperativity parameter of the model are simply stated. Experimental tests of model predictions for the responses of DNA to small applied stretching, twisting, and bending stresses, and changes in temperature, ionic conditions, and sequence are described. The mechanism and significance of the large cooperativity, which enables significant DNA responses to such small perturbations, are also noted. The capacity of the model to resolve a number of long-standing and sometimes interconnected puzzles in the extant literature, including the origin of the broad pre-melting transition studied by numerous workers in the 1960s and 1970s, is demonstrated. Under certain conditions, the model predicts significant long-range attractive or repulsive interactions between hypothetical proteins with strong preferences for one or the other state that are bound to well-separated sites on the same DNA. A scenario is proposed for the activation of the ilvPG promoter on a supercoiled DNA by integration host factor.

A quantitative model of a cooperative two-state equilibrium in DNA: experimental tests, insights, and predictions

Quantitative parameters for a two-state cooperative transition in duplex DNAs were finally obtained during the last 5 years. After a brief discussion of observations pertaining to the existence of the two-state equilibrium per se, the lengths, torsion, and bending elastic constants of the two states involved and the cooperativity parameter of the model are simply stated. Experimental tests of model predictions for the responses of DNA to small applied stretching, twisting, and bending stresses, and changes in temperature, ionic conditions, and sequence are described. The mechanism and significance of the large cooperativity, which enables significant DNA responses to such small perturbations, are also noted. The capacity of the model to resolve a number of long-standing and sometimes interconnected puzzles in the extant literature, including the origin of the broad pre-melting transition studied by numerous workers in the 1960s and 1970s, is demonstrated. Under certain conditions, the model predicts significant long-range attractive or repulsive interactions between hypothetical proteins with strong preferences for one or the other state that are bound to well-separated sites on the same DNA. A scenario is proposed for the activation of the ilvPG promoter on a supercoiled DNA by integration host factor.

Temperature-dependence of the bending elastic constant of DNA and extension of the two-state model. Tests and new insights

Review and analyses of the experimental data indicate that in nearly all cases bending elastic constants of the effective springs between bp of DNA actually undergo a net increase with increasing T from 278 to 315 K. The exceptions to this rule are bending elastic constants obtained from equilibrium topoisomer distributions of a 2686 bp pUC19 DNA by assuming a fixed T-independent value of the torsion elastic constant. When the same data are analyzed using measured T-dependent values of the torsion elastic constant, which decline with increasing T, a modest increase in bending elastic constant with increasing T is obtained. After revising the torsion elastic constants of the previously formulated two-state cooperative transition model to account for additional data, that model is fitted to the bending elastic constants reckoned from the aforementioned topoisomer distributions to determine the best-fit values for each state. The rather good fit implies a strong negative linear correlation between the inverse bending and inverse torsion elastic constants as T is varied. Predictions of the resulting two-state model, wherein each state has fixed bending and torsion elastic constants, agree surprisingly well with single-molecule relative extension and torque data. The same model also yields good agreement with numerous other experimental data. With increasing T the equilibrium is shifted from the (longer, torsionally stiffer, flexurally softer) b-state toward the (shorter, torsionally softer, flexurally stiffer) a-state. This transition is suggested to be the origin of the so-called broad pre-melting transition exhibited by many, but not all, DNAs.

Grain boundary premelting of monolayer ices in 2D nano-channels

ABSTRACTWe employ molecular-dynamics (MD) simulations to study grain boundary (GB) premelting in ices confined in two-dimensional hydrophobic nano-channels. Premelting transitions are observed in symmetric tilt GBs in monolayer ices and involve the formation of a premelting band of liquid phase water with a width that grows logarithmically as the melting temperature is approached from below, consistent with the existing theory of GB premelting. The premelted GB is found rough for a broad range of temperature below the melting temperature, the two ice-premelt interfaces bounding the melted GB are engaged with long wave-length parallel coupled fluctuations. Based on current MD simulation study, one may expect GB premelting transitions exist over a wide range of low dimensional phases of confined ice and shows important consequences for crystal growth of low dimensional ices.

Molecular dynamics simulation of charged colloids confined between hard walls: pre-melting and pre-freezing across the BCC–fluid coexistence

Molecular dynamics (MD) computer simulations are used to study the structure of hard-core Yukawa systems confined between two parallel hard walls. States around the coexistence between a fluid and a body-centered cubic (BCC) crystal are considered. In all cases a pronounced layering in the vicinity of the walls is observed. Using a thermodynamic integration scheme, we determine the wall-fluid interfacial free energy γ which is negative and monotonically decreasing with increasing bulk density of the fluid. In the case of the fluid, the layers next to the walls undergo a transition from a fluid to a hexagonal structure. This pre-freezing transition occurs well below the coexistence bulk density of the fluid. The confined BCC crystal in (111) orientation shows melted regions between crystalline face-centered cubic (FCC) layers close to the wall and the BCC bulk region.

Measurement and interpretation of the thermo-physical properties of UO2 at high temperatures: The viral effect of oxygen defects

Abstract Values are reported of specific heat, thermal conductivity and thermal diffusivity of UO2 from 1500 K to 2900 K based on laser flash measurements. Experiment is complemented by the development of solid state physics models that aid in the interpretation of the results. Specific heat is shown to exhibit a smooth maximum at 2715 K ± 100 K, consistent with a competition between two processes - oxygen defect interactions (net attraction) and saturation of oxygen interstitial sites. The specific heat model and measurements show, for the first time that a gradual pre-melting transition is consistent with high temperature literature values – enthalpy increment measurements and independently measured high temperature oxygen defect concentrations. Thermal conductivity exhibits a minimum consistent with: 1) an increase in electronic thermal conductivity via polaron production and mobilization and 2) degradation in lattice thermal conductivity due to phonon - phonon scattering and phonon - defect scattering. It is predicted that the high concentration of oxygen defects should contribute significantly to electrical conductivity and thermal expansion at high temperatures.

Grain boundary complexions and pseudopartial wetting

Abstract The important class of grain boundary (GB) complexions includes the few nanometer thick layers having composition which strongly differs from that of the abutting grains. Such GB complexions are frequently called intergranular films (IGFs) and can be observed close to the lines of wetting, prewetting and premelting complexion transitions in the bulk phase diagrams. In the majority of systems, the direct transition between complete and partial GB wetting takes place (by changing temperature, pressure, etc.). However, in certain conditions the so-called pseudopartial (or pseudoincomplete, or frustrated complete) GB wetting appears in a phase diagram between complete and partial wetting. In case of pseudopartial GB wetting, the thin GB layer of a complexion (IGF or 2-D interfacial phase) can coexist with large droplets (or particles) of the wetting phase with a non-zero dihedral (contact) angle. Thus, such IGFs can be observed in the two-phase (or multiphase) fields of bulk phase diagrams, in the broad intervals of concentrations, temperature and/or pressure. The IGFs driven by the pseudopartial GB wetting can drastically modify the properties of polycrystals. In this review, we discuss this phenomenon for the technologically important Fe–Nd–B-based hard magnetic alloys, WC–Co cemented carbides and Al-based light alloys.

Modes of surface premelting in colloidal crystals composed of attractive particles.

Crystal surfaces typically melt into a thin liquid layer at temperatures slightly below the melting point of the crystal. Such surface premelting is prevalent in all classes of solids and is important in a variety of metallurgical, geological and meteorological phenomena. Premelting has been studied using X-ray diffraction and differential scanning calorimetry, but the lack of single-particle resolution makes it hard to elucidate the underlying mechanisms. Colloids are good model systems for studying phase transitions because the thermal motions of individual micrometre-sized particles can be tracked directly using optical microscopy. Here we use colloidal spheres with tunable attractions to form equilibrium crystal-vapour interfaces, and study their surface premelting behaviour at the single-particle level. We find that monolayer colloidal crystals exhibit incomplete premelting at their perimeter, with a constant liquid-layer thickness. In contrast, two- and three-layer crystals exhibit conventional complete melting, with the thickness of the surface liquid diverging as the melting point is approached. The microstructures of the surface liquids differ in certain aspects from what would be predicted by conventional premelting theories. Incomplete premelting in the monolayer crystals is triggered by a bulk isostructural solid-solid transition and truncated by a mechanical instability that separately induces homogeneous melting within the bulk. This finding is in contrast to the conventional assumption that two-dimensional crystals melt heterogeneously from their free surfaces (that is, at the solid-vapour interface). The unexpected bulk melting that we observe for the monolayer crystals is accompanied by the formation of grain boundaries, which supports a previously proposed grain-boundary-mediated two-dimensional melting theory. The observed interplay between surface premelting, bulk melting and solid-solid transitions challenges existing theories of surface premelting and two-dimensional melting.

Melting point trends and solid phase behaviors of model salts with ion size asymmetry and distributed cation charge.

The melting point trends of model salts composed of coarse grain ions are examined using NPT molecular dynamics simulations. The model salts incorporate ion size asymmetry and distributed cation charge, which are two common features in ionic liquids. A series of single-phase and two-phase simulations are done at set temperatures with 50 K intervals for each salt, and the normal melting point is estimated within 50 K. The melting point trends are then established relative to a charge-centered, size symmetric salt with a normal melting point between 1250 K and 1300 K. We consider two sets of size asymmetric salts with size ratios up to 3:1; the melting point trends are different in each set. The lowest melting point we find is between 450 K and 500 K, which is a reduction of over 60% from the charge-centered, size symmetric case. In both sets, we find diversity in the solid phase structures. For all size ratios with small cation charge displacements, the salts crystallize with orientationally disordered cations. When the partial cation charge is far enough off-center in salts with ion size ratios near 1:1, the salts can become trapped in glassy states and have underlying crystal structures that are orientationally ordered. At ion size ratios near 3:1, the salts with large cation charge displacements show premelting transitions at temperatures as low as 300 K. After the premelting transition, these salts exist either as fast ion conductors, where the smaller anions move through a face centered cubic (fcc) cation lattice, or as plastic crystals, where ion pairs rotate on a fcc lattice.

Solid state 13C-NMR, infrared, X-ray powder diffraction and differential thermal studies of the homologous series of some mono-valent metal (Li, Na, K, Ag) n-alkanoates: A comparative study

A comparative study of the molecular packing, lattice structures and phase behaviors of the homologous series of some mono-valent metal carboxylates (Li, Na, K and Ag) is carried out via solid state FT-infrared and (13)C-NMR spectroscopes, X-rays powder diffraction, density measurements, differential scanning calorimetry, polarizing light microscopy and variable temperature infrared spectroscopy. It is proposed that, for lithium, sodium and potassium carboxylates, metal-carboxyl coordination is via asymmetric chelating bidentate bonding with extensive intermolecular interactions to form tetrahedral metal centers, irrespective of chain length. However, for silver n-alkanoates, carboxyl moieties are bound to silver ions via syn-syn type bridging bidentate coordination to form dimeric units held together by extensive head group inter-molecular interactions. Furthermore, the fully extended hydrocarbon chains which are crystallized in the all-trans conformation are tilted at ca. 30°, 27°, 15° and 31° with respect to a normal to the metal plane, for lithium, sodium, silver and potassium carboxylates, respectively. All compounds are packed as lamellar bilayer structures, however, lithium compounds are crystallized in a triclinic crystal system whilst silver, sodium and potassium n-alkanoates are all monoclinic with possible P1 bravais lattice. Odd-even alternation observed in various physical features is associated with different inter-planar spacing between closely packed layers in the bilayer which are not in the same plane; a phenomenon controlled by lattice packing symmetry requirements. All compounds, except silver carboxylates, show partially reversibly first order pre-melting transitions; the number of which increases with increasing chain length. These transitions are associated, for the most part, with lamellar collapse followed by increased gauche-trans isomerism in the methylene group assembly, irrespective of chain length. It is proposed that the absence of mesomorphic transitions in their phase sequences is due to a lack of sufficient balance between attractive and repulsive electrostatic and van der Waals forces during phase change. The evidence presented in this study shows that phase behaviors of mono-valent metal carboxylates are controlled, mainly, by head group bonding.

Quantum Size Effects in the Size–Temperature Phase Diagram of Gallium: Structural Characterization of Shape‐Shifting Clusters

Finite temperature analysis of cluster structures is used to identify signatures of the low-temperature polymorphs of gallium, based on the results of first-principle Born-Oppenheimer molecular dynamics simulations. Pre-melting structural transitions proceed from either the β- and/or the δ-phase to the γ- or δ-phase, with a size- dependent phase progression. We relate the stability of each isomer to the electronic structures of the different phases, giving new insight into the origin of polymorphism in this complicated element.

Thermal conductivity and energetic recoils in UO2 using a many-body potential model

Classical molecular dynamics simulations have been performed on uranium dioxide (UO2) employing a recently developed many-body potential model. Thermal conductivities are computed for a defect free UO2 lattice and a radiation-damaged, defect containing lattice at 300 K, 1000 K and 1500 K. Defects significantly degrade the thermal conductivity of UO2 as does the presence of amorphous UO2, which has a largely temperature independent thermal conductivity of ∼1.4 Wm−1 K−1. The model yields a pre-melting superionic transition temperature at 2600 K, very close to the experimental value and the mechanical melting temperature of 3600 K, slightly lower than those generated with other empirical potentials. The average threshold displacement energy was calculated to be 37 eV. Although the spatial extent of a 1 keV U cascade is very similar to those generated with other empirical potentials and the number of Frenkel pairs generated is close to that from the Basak potential, the vacancy and interstitial cluster distribution is different.

Theoretical investigation of structural and thermo-mechanical properties of thoria up to 3300 K temperature

Abstract We investigated thermo-mechanical properties of thoria up to a very high temperature (3300 K). We demonstrate that, using first-principles molecular dynamics, it is possible to predict thermal expansion of thoria in agreement with experiment. The new generalized gradient approximation functional within the density functional theory predicts, in agreement with experiment, not only the relative thermal expansion but also the absolute values of the lattice constant as a function of temperature. The molecular dynamics approach has an advantage over the previously used quasi-harmonic method, because it can be used even at temperatures (above 2700 K) where the longitudinal optical mode breaks in thoria. The calculated phonon dispersion agrees well with the experimental relation, measured using inelastic neutron scattering. The temperature, at which the negative frequency in the optical mode appears, coincides with the λ-type pre-melting transition reported in thoria.

Phase-field modeling of grain-boundary premelting using obstacle potentials.

We investigate the multiorder parameter phase field model of Steinbach and Pezzolla [Physica D 134, 385 (1999)] concerning its ability to describe grain boundary premelting. For a single order parameter situation solid-melt interfaces are always attractive, which allows us to have (unstable) equilibrium solid-melt-solid coexistence above the bulk melting point. The temperature-dependent melt layer thickness and the disjoining potential, which describe the interface interaction, are affected by the choice of the thermal coupling function and the measure to define the amount of the liquid phase. Due to the strictly finite interface thickness the interaction range also is finite. For a multiorder parameter model we find either purely attractive or purely repulsive finite-ranged interactions. The premelting transition is then directly linked to the ratio of the grain boundary and solid-melt interfacial energy.