180 likes | 297 Vues
This study investigates the relationship between nucleation rates and the microscopic energy of small cluster formation within a model Lennard-Jones system. By calculating growth and decay rate constants at varying temperatures, we derive nucleation rates from kinetic formalism and analyze the scaling behavior through the comparison of logJ versus lnS. Using Monte Carlo simulations, we explore the dependence of nucleation on supersaturation ratios, revealing evidence of scaling emerging from the temperature-dependent free energy differences. Our findings contribute to the understanding of non-ideal gas behavior in cluster formation.
E N D
A Study of Scaled Nucleationin a Model Lennard-Jones System Barbara Hale and Tom Mahler Physics Department University of Missouri – RollaJerry Kiefer Physics Department St. Bonaventure University
Motivation To understand how scaling of the nucleation rate is related to the microscopic energy of formation of small clusters.
Scaling:Wölk and Strey Water DataCo = [Tc/240-1]3/2 ; Tc = 647.3 K B. Hale, J. Chem. Phys. 122, 204509 (2005)
Schmitt et al Toluene (C7H8) data Co = [Tc /240-1]3/2 ; Tc = 591.8K
Study of Scaling in LJ System • calculate rate constants for growth and decay of model Lennard-Jones clusters at three temperatures; • determine model nucleation rates from kinetic nucleation rate formalism; • compare logJ vs lnS and logJ vs lnS/[Tc/T-1]3/2
Model Lennard-Jones System Law of mass action dilute vapor system with volume, V; non-interacting mixture of ideal gases; each n-cluster size is ideal gas of Nn particles; full atom-atom LJ interaction potential; separable classical Hamiltonian
Law of Mass Action Nn/[Nn-1N1] = Q(n)/[Q(n-1)Q(1)n] Q(n) = n-cluster canonical configurational partition function
Relation to Growth & Decay (n-1)Nn-1N1= (n)Nn we calculate: Q(n)/[Q(n-1)Q(1)n]= Nn/[Nn-1N1] = (n-1)/(n)
Kinetic Nucleation Rate Formalism 1/J = n=1,M 1/Jn ; M large Jn = (n)(N1S)2j=2,n [N1S(j-1)/(j)] S = N1exp/N1
Free Energy Differences - f(n) = ln [Q(n)/(Q(n-1)Q(1))]calculated = ln [ (ρoliq/ρovap)(j-1)/(j) ] Use Monte Carlo Bennett technique.
Monte Carlo Simulations Ensemble A: (n -1) cluster plus monomer probe interactions turned off Ensemble B: n cluster with normal probe interactions Calculate f(n) =[F(n)-F(n-1)]/kT
The nucleation rate can be calculated for a range of supersaturation ratios, S. 1/J = n=1,M 1/Jn ; M large Jn = (n)(N1S)2j=2,n [N1S(j-1)/(j)] S = N1exp/N1
Comments & Conclusions • Experimental data indicate that Jexp is a function of lnS/[Tc/T-1]3/2. • No “first principles” derivation of scaling exists. • Monte Carlo LJ cluster simulations show evidence of scaling. • Scaling appears to emerge from [Tc/T-1] dependence of the f(n).