• 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

    • Fulltext


        Click here to view fulltext PDF

      Permanent link:

    • Keywords


      polydispersity index; Lindemann ratio; Hansen-Verlet rule; bifurcation diagram; terminal polydispersity; inherent structure energy

    • Abstract


      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.

    • Author Affiliations


      Sarmistha Sarkar1 2 Rajib Biswas1 Partha Pratim Ray2 Biman Bagchi1

      1. Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560 012, India
      2. Department of Physics, Jadavpur University, Kolkata 700 032, India
    • Dates

  • Journal of Chemical Sciences | News

© 2021-2022 Indian Academy of Sciences, Bengaluru.