We study melting of a face-centered crystalline solid consisting of polydisperse Lennard-Jones spheres with Gaussian polydispersity in size. The phase diagram reproduces the existence of a nearly temperature invariant terminal polydispersity (δ(t) =/~ 0.11), with no signature of reentrant melting. The absence of reentrant melting can be attributed to the influence of the attractive part of the potential upon melting. We find that at terminal polydispersity the fractional density change approaches zero, which seems to arise from vanishingly small compressibility of the disordered phase. At constant temperature and volume fraction the system undergoes a sharp transition from crystalline solid to the disordered amorphous or fluid state with increasing polydispersity. This has been quantified by second- and third-order rotational invariant bond orientational order, as well as by the average inherent structure energy. The translational order parameter also indicates a similar sharp structural change at δ =/~ 0.09 in case of T(*) = 1.0, φ = 0.58. The free energy calculation further supports the sharp nature of the transition. The third-order rotationally invariant bond order shows that with increasing polydispersity, the local cluster favors a more icosahedral arrangement and the system loses its local crystalline symmetry. Interestingly, the value of structure factor S(k) of the amorphous phase at δ =/~ 0.10 (just beyond the solid-liquid transition density at T(*) = 1) becomes 2.75, which is below the value of 2.85 required for freezing given by the empirical Hansen-Verlet rule of crystallization, well known in the theory of freezing.

No abstract

Aqueous binary mixtures often exhibit dramatic departure from the predicted hydrodynamic behavior when transport properties are plotted against composition. We show by inherent structure (IS) analysis that this sharp composition dependent breakdown of the Stokes–Einstein relation can be attributed to the non-monotonic variation in the average inherent structure energy of these mixtures. Further IS analysis reveals the existence of a unique ground state, stabilized by both the formation of an optimum number of H-bonds and a favorable hydrophobic interaction at this composition. The surprisingly sharp turnaround behavior observed in the effective hydrodynamic radius also owes its origin to the same combination of these two factors. Interestingly, the temperature dependence of isothermal compressibility shows a minimum at the particular composition. Extensive studies on water–dimethyl sulfoxide and water–ethanol mixtures using two different force-fields of water reveal many features that are nearly universal. A justification of this quasi-universal behavior is provided in terms of a mode-coupling theory (MCT) of viscosity, which can serve as the starting point of a remarkable correlation observed with the nearest neighbor structure, as captured by the first peaks of the radial distribution function, and the slowdown in the intermediate scattering function at intermediate wavenumbers. Therefore, the formation of the local structure captured through IS analysis can be correlated with the MCT.

No abstract

An energy landscape view of phase separation and nonideality in binary mixtures is developed by exploring their potential energy landscape (PEL) as functions of temperature and composition. We employ molecular dynamics simulations to study a model that promotes structure breaking in the solute-solvent parent binary liquid, at low temperatures. The PEL of the system captures the potential energy distribution of the inherent structures (IS) of the system and is obtained by removing the kinetic energy (including that of intermolecular vibrations). The broader distribution of the inherent structure energy for structure breaking liquid than that of the structure making liquid demonstrates the larger role of entropy in stabilizing the parent liquid of the structure breaking type of binary mixtures. At high temperature, although the parent structure of the structure breaking binary mixture is homogenous, the corresponding inherent structure is found to be always phase separated, with a density pattern that exhibits marked correlation with the energy of its inherent structure. Over a broad range of intermediate inherent structure energy, bicontinuous phase separation prevails with interpenetrating stripes as signatures of spinodal decomposition. At low inherent structure energy, the structure is largely phase separated with one interface where as at high inherent structure energy we find nucleation type growth. Interestingly, at low temperature, the average inherent structure energy () exhibits a drop with temperature which signals the onset of crystallization in one of the phases while the other remains in the liquid state. The nonideal composition dependence of viscosity is anticorrelated with average inherent structure energy.

A systematic procedure is outlined for scaling analysis of momentum and heat transfer in gas tungsten arc weld pools. With suitable selections of non-dimensionalised parameters, the governing equations coupled with appropriate boundary conditions are first scaled, and the relative significance of various terms appearing in them is analysed accordingly. The analysis is then used to predict the orders of magnitude of some important quantities, such as the velocity scale at the top surface, velocity boundary layer thickness, maximum temperature increase in the pool, and time required for initiation of melting. Some of the quantities predicted from the scaling analysis can also be used for optimised selection of appropriate grid size and time steps for full numerical simulation of the process. The scaling predictions are finally assessed by comparison with numerical results quoted in the literature, and a good qualitative agreement is observed.

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 (δ P) above which no crystalline solid is stable. This transition value of polydispersity (termed as transition polydispersity, δ 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, δ P shifts to higher values. This temperature and number density dependent value of δ P saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (δ TP). This value (δ TP ∼ 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 (δ 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.

We explore the potential energy landscape of structure breaking binary mixtures (SBBM) where two constituents dislike each other, yet remain macroscopically homogeneous at intermediate to high temperatures. Interestingly, we find that the origin of strong composition dependent non-ideal behaviour lies in its phase separated inherent structure. The inherent structure (IS) of SBBM exhibits bi-continuous phase as is usually formed during spinodal decomposition. We draw analogy of this correlation between non-ideality and phase separation in IS to explain observation of non-ideality in real aqueous mixtures of small amphiphilic solutes, containing both hydrophilic and hydrophobic groups. Although we have not been able to obtain IS of these liquids, we find that even at room temperature these liquids sustain formation of fluctuating, transient bicontinuous phase, with limited lifetime (τ 20 ps). While in the model (A, B) binary mixture, the non-ideal composition dependence can be considered as a fluctuation from a phase separated state, a similar scenario is expected to be responsible for the unusually strong non-ideality in these aqueous binary mixtures.