Thermodynamic Models for the Size-dependent Melting of Nanoparticles: Different Hypotheses

A careful comparison of the experimental results reported in the literature reveals different variations of the melting temperature even for the same materials. Though there are different theoretical models, thermodynamic model has been extensively used to understand different variations of size-dependent melting of nanoparticles. There are different hypotheses such as homogeneous melting (HMH), liquid nucleation and growth (LNG) and liquid skin melting (LSM) to resolve different variations of melting temperature as reported in the literature. HMH and LNG account for the linear variation where as LSM is applied to understand the nonlinear behaviour in the plot of melting temperature against reciprocal of particle size. However, a bird's eye view reveals that either HMH or LSM has been extensively used by experimentalists. It has also been observed that not a single hypothesis can explain the size-dependent melting in the complete range. Therefore we describe an approach which can predict the plausible hypothesis for a given data set of the size-dependent melting temperature. A variety of data have been analyzed to ascertain the hypothesis and to test the approach.

Simulation of Influence of Bilayer Melting on Dynamics and Thermodynamics of Interfacial Water

We investigate the effect of bilayer melting transition on thermodynamics and dynamics of interfacial water using molecular dynamics simulation with the two-phase thermodynamic model. We show that the diffusivity of interface water depicts a dynamic crossover at the chain melting transition following an Arrhenius behavior until the transition temperature. The corresponding change in the diffusion coefficient from the bulk to the interface water is comparable with experimental observations found recently for water near 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) vesicles Phys. Chem. Chem. Phys. 13, 7732 (2011)]. The entropy and potential energy of interfacial water show distinct changes at the bilayer melting transition, indicating a strong correlation in the thermodynamic state of water and the accompanying first-order phase transition of the bilayer membrane. DOI: 10.1103/PhysRevLett.110.018303

Anisotropy induced crossover from weakly to strongly first order melting of two dimensional solids

Melting and freezing transitions in two dimensional (2D) systems are known to show highly unusual characteristics. Most of the earlier studies considered atomic systems: the melting of 2D molecular solids is still largely unexplored. In order to understand the role of anisotropy as well as multiple energy and length scales present in molecular systems, here we report computer simulation studies of melting of 2D molecular systems. We computed a limited portion of the solid-liquid phase diagram. We find that the interplay between the strength of isotropic and anisotropic interactions can give rise to rich phase diagram consisting of isotropic liquid and two crystalline phases-honeycomb and oblique. The nature of the transition depends on the relative strength of the anisotropic interaction and a strongly first order melting turns into a weakly first order transition on increasing the strength of the isotropic interaction. This crossover can be attributed to an increase in stiffness of the solid phase free energy minimum on increasing the strength of the anisotropic interaction. The defects involved in melting of molecular systems are quite different from those known for the atomic systems.

Remote sensing study of snowline altitude at the end of melting season, Chandra-Bhaga basin, Himachal Pradesh, 1980-2007

Glaciers have a direct relation with climate change. The equilibrium line altitude (ELA) is the most useful parameter to study the effect of climate change on glaciers. The ELA is a theoretical snowline at which the glacier mass balance is zero. Snowline altitude (SLA) at the end of melting season is generally regarded as the ELA. Glaciers of Chandra-Bhaga basin in Lahaul-Spiti district of Himachal Pradesh were chosen to study the ELA, using satellite images from 1980 to 2007. A total of 19 glaciers from the Chandra-Bhaga basin were identified and selected to carry out the study of ELA variation over 27 years. This study reveals that the mean SLA of the sub-basin has increased from 5009 +/- 61m to 5401 +/- 21m from 1980 to 2007. The study is in agreement with the reported increase in the temperature and decrease in winter snowfall of North-West Himalaya in the last century.

Melting and solidification behavior of Pb-Sn embedded alloy nano-particles

Nano-sized bimetallic dispersoids consisting of (Pb) and beta-(Sn) phases of eutectic composition (Pb26.1Sn73.9) embedded in aluminum and Al-Cu-Fe quasicrystalline matrices have been prepared by rapid solidification processing. The two phases, face centered cubic (Pb) and body center tetragonal, beta-(Sn) solid solution co-exist in all the embedded nanoparticles at room temperature. The phases bear crystallographic orientation relationship with the matrix. In situ TEM study has been carried out for the alloy particles to study the melting and the solidification behavior. The detailed microscopic observations indicate formation of a single-phase metastable fcc (Pb) in the nano-particles prior to the melting during heating. Solidification of these particles begins with nucleation of fcc (Pb), which phase separates into fcc (Pb) and beta-(Sn) lamellae in the solid state. In situ X-ray diffraction study is carried out to obtain lattice parameter of metastable fcc (Pb) and thereby an estimate of amount of Sn dissolved in the metastable (Pb) prior to the melting. The results are discussed in terms of a metastable phase diagram between fcc Pb and fcc Sn and invoking the size effect on the metastable phase diagram. The size factor is found to play a critical role in deciding the pathway of phase transformation as well as the extension of solid solubility of Sn in fcc (Pb) in the nano-particles.

Energy funneling and macromolecular conformational dynamics: a 2D Raman correlation study of PEG melting

Cascading energy landscapes through funneling has been postulated as a mechanistic route for achieving the lowest energy configuration of a macromolecular system (such as proteins and polymers). In particular, understanding the molecular mechanism for the melting and crystallization of polymers is a challenging fundamental question. The structural modifications that lead to the melting of poly(ethylene glycol) (PEG) are investigated here. Specific Raman bands corresponding to different configurations of the PEG chain have been identified, and the molecular structural dynamics of PEG melting have been addressed using a combination of Raman spectroscopy, 2D Raman correlation and density functional theory (DFT) calculations. The melting dynamics of PEG have been unambiguously explained along the C-O bond rotation coordinate.

Comparative Study of Protein Unfolding in Aqueous Urea and Dimethyl Sulfoxide Solutions: Surface Polarity, Solvent Specificity, and Sequence of Secondary Structure Melting

Elucidation of possible pathways between folded (native) and unfolded states of a protein is a challenging task, as the intermediates are often hard to detect. Here, we alter the solvent environment in a controlled manner by choosing two different cosolvents of water, urea, and dimethyl sulfoxide (DMSO) and study unfolding of four different proteins to understand the respective sequence of melting by computer simulation methods. We indeed find interesting differences in the sequence of melting of alpha helices and beta sheets in these two solvents. For example, in 8 M urea solution, beta-sheet parts of a protein are found to unfold preferentially, followed by the unfolding of alpha helices. In contrast, 8 M DMSO solution unfolds alpha helices first, followed by the separation of beta sheets for the majority of proteins. Sequence of unfolding events in four different alpha/beta proteins and also in chicken villin head piece (HP-36) both in urea and DMSO solutions demonstrate that the unfolding pathways are determined jointly by relative exposure of polar and nonpolar residues of a protein and the mode of molecular action of a solvent on that protein.

Nonequilibrium ``Melting'' of a Charge Density Wave Insulator via an Ultrafast Laser Pulse

We employ an exact solution of the simplest model for pump-probe time-resolved photoemission spectroscopy in charge-density-wave systems to show how, in nonequilibrium, the gap in the density of states disappears while the charge density remains modulated, and then the gap reforms after the pulse has passed. This nonequilibrium scenario qualitatively describes the common short-time experimental features in TaS2 and TbTe3, indicating a quasiuniversality for nonequilibrium ``melting'' with qualitative features that can be easily understood within a simple picture.

Succinate esters: odd-even effects in melting points

Dialkyl succinates show a pattern of alternating behavior in their melting points, as the number of C atoms in the alkane side chain increases, unlike in the dialkyl oxalates Joseph et al. (2011). Acta Cryst. B67, 525-534]. Dialkyl succinates with odd numbers of C atoms in the alkyl side chain show higher melting points than the immediately adjacent analogues with even numbers. The crystal structures and their molecular packing have been analyzed for a series of dialkyl succinates with 1 - 4 C atoms in the alkyl side chain. The energy difference (Delta E) between the optimized and observed molecular conformations, density, Kitaigorodskii packing index (KPI) and C-H center dot center dot center dot O interactions are considered to rationalize this behavior. In contrast to the dialkyl oxalates where a larger number of moderately strong C-H center dot center dot center dot O interactions were characteristic of oxalates with elevated melting points, here the molecular packing and the density play a major role in raising the melting point. On moving from oxalate to succinate esters the introduction of the C2 spacer adds two activated H atoms to the asymmetric unit, resulting in the formation of stronger C-H center dot center dot center dot O hydrogen bonds in all succinates. As a result the crystallinity of long-chain alkyl substituted esters improves enormously in the presence of hydrogen bonds from activated donors.

Use of polydispersity index as control parameter to study melting/freezing of Lennard-Jones system: Comparison among predictions of bifurcation theory with Lindemann criterion, inherent structure analysis and Hansen-Verlet rule

Using polydispersity index as an additional order parameter we investigate freezing/melting transition of Lennard-Jones polydisperse systems (with Gaussian polydispersity in size), especially to gain insight into the origin of the terminal polydispersity. The average inherent structure (IS) energy and root mean square displacement (RMSD) of the solid before melting both exhibit quite similar polydispersity dependence including a discontinuity at solid-liquid transition point. Lindemann ratio, obtained from RMSD, is found to be dependent on temperature. At a given number density, there exists a value of polydispersity index (delta (P)) above which no crystalline solid is stable. This transition value of polydispersity(termed as transition polydispersity, delta (P) ) is found to depend strongly on temperature, a feature missed in hard sphere model systems. Additionally, for a particular temperature when number density is increased, delta (P) shifts to higher values. This temperature and number density dependent value of delta (P) saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (delta (TP)). This value (delta (TP) similar to 0.11) is in excellent agreement with the experimental value of 0.12, but differs from hard sphere transition where this limiting value is only 0.048. Terminal polydispersity (delta (TP)) thus has a quasiuniversal character. Interestingly, the bifurcation diagram obtained from non-linear integral equation theories of freezing seems to provide an explanation of the existence of unique terminal polydispersity in polydisperse systems. Global bond orientational order parameter is calculated to obtain further insights into mechanism for melting.

On the issues of resolving a low melting combination as a definite eutectic or an elusive cocrystal: A critical evaluation

Cocrystals and eutectics are different yet related crystalline multi-component adducts with diverse applications in pharmaceutical and materials fields. Recently, they were shown to be alternate products of cocrystallization experiments. Whereas a cocrystal shows distinct diffraction, spectroscopic and thermal signatures as compared to parent components, the hallmark of a eutectic is its low melting nature. However, in certain cases, there can be a problem when one resorts to design a cocrystal and assess its formation vis-A -vis a eutectic. In the absence of a gold standard method to make a cocrystal, it is often difficult to judge how exhaustive should the cocrystallization trials be to ensure the accomplishment of a desired/putative cocrystal. Further, a cocrystal can manifest with intermolecular interactions and/or crystal structure similar to that of its parent compounds such that the conventional diffraction and spectroscopic techniques will be of little help to conclusively infer the formation of cocrystal in the lack of single crystals. Such situations combined with low melting behavior of a combination brings the complication of resolving the combination as a cocrystal or eutectic since now both the adducts share common features. Based on the curious case of Caffeine-Benzoic acid combination, this study aims to unfold the intricate issues related to the design, formation and characterization of cocrystals and eutectics for a way forward. The utility of heteronuclear seeding methodology in establishing a given combination as a cocrystal-forming one or a eutectic-forming one in four known systems is appraised.

Influence of Local Magnetic Fields on P-Doped Si Floating Zone Melting Crystal Growth in Microgravity

A two-dimensional axisymmetric numerical model is presented to study the influence of local magnetic fields on P-doped Si floating zone melting crystal growth in microgravity. The model is developed based on the finite difference method in a boundary-fitted curvilinear coordinate system. Extensive numerical simulations are carried out, and parameters studied include the curved growth interface shape and the magnetic field configurations. Computed results show that the local magnetic field is more effective in reducing the impurity concentration nonuniformity at the growth interface in comparison with the longitudinal magnetic field. Moreover, the curved growth interface causes more serious impurity concentration nonuniformity at the growth interface than the case with a planar growth interface.

Application of enthalpy method for moving boundary problem in laser surface melting

A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization. The enthalpy method was applied to solve this two-phase axisymmetrical melting problem Computational results of temperature fields were obtained, which provide useful information to practical laser treatment processing. The validity of enthalpy method in solving such problems is presented.

Influence of grain boundary on melting

The temperature behaviour of an Al bicrystal with surfaces consisting of (110) and (111) crystals is simulated using molecular dynamics. The result shows that the (110) crystal losses its crystalline order at 820K, whereas the disorder does not propagate through the (111) crystal at this temperature. Instead, some disordered atoms are recrystallized into the (111) crystal and the initial grain boundary changes into a stable order-disorder interface. Thus, it was discovered that at a temperature near its melting point, the (111) crystal grew and obstructed the propagation of disorder. Such an obstruction is helpful for understanding melting.

Numerical Simulation of Transient Thermal Field and Residual Stress in Laser Melting Process

以激光熔凝表面强韧化处理为背景,应用空间弹塑性有限单元和高精度数值算法同时考虑材料组织性能的变化模拟工件的温度场及残余应力,研究激光熔凝加工中瞬时温度场及残余应力数值模拟,同时考虑相变潜热及相变塑性的影响,用算例验证了模型的正确性,给出了不同时刻温度场分布及残余应力分布。