Thermal Surface Fluctuations of Clusters with Long-Range Interaction

1. Thermal Surface Fluctuations of Clusters withLong-Range Interaction D.I. Zhukhovitskii Joint Institute for High Temperatures, RAS

2. Liquid―vapor interface structure: smooth or stratified ? (van der Waals) (Gibbs) Gas Gas Intermediate phase Liquid Liquid

3. Aim of research: 1. Working out a proper method for MD simulation of such clusters in vapor environment. 2. Calculation of slice spectra. 3. Estimation of fission threshold. 4. Development of a theory of surface fluctuations for clusters with long–range interaction.

4. System under consideration Cluster particles are assumed to interact via the pair additive potential where and the long–range component

5. Molecular dynamics simulation Systems with multiple length and time scales require special integrators to prevent enormous energy drift. In the force rotation approach, an artificial torque of the long–range force components Fi arising from cluster rotation is removed by rotation of these forces. We impose the condition y1, y2, andy3are the Euler angles. They are solutions of equation set

6. Simulation cell: a cluster in equilibrium vapor environment

7. Definition of a cluster: a particle belongs to the cluster if it has at least one neighbor particle at the distance less than rb, which belongs to this cluster. The problem is, how to definerb.

9. We define three particle types: internaland surfaceparticlesand virtual chains. Particle 1 with the radius vector r1 that belongs to the cluster will be called internal if there exists at least one particle 2 with radius vector r2 belonging to the same cluster that forms more than four bonds such that the conditions are satisfied.

11. Average configurations yield smooth density distribution inside the transitional region: < >

12. We isolate the surface particles (pivot particles) situated between two parallel planes. The particle polar coordinates are the values of a continuous function The slice spectrum are defined as the averages both over configurations and over the Euler cluster rotation angles: The total spectrum is a sum of the capillary fluctuations (CF) and bulk fluctuations (BF) spectra.

13. Bulk fluctuations arise from randomicity of particles location. Hence, they coincide with fluctuations of the surface particles of a cluster truncated by a sphere. The total spectral density is

14. Different components of the surface fluctuations spectral amplitudes for a cluster comprising 30000 particles at a= 0 and the temperature of 0.75interparticle potential well depth. (1),bulk fluctuations, Sk= Rk; (2), capillary fluctuations, Sk = Qk;(3), the total spectral amplitude, Sk = Sk;(4), the total spectral amplitude without isolation of the virtual chains.

15. CF spectral amplitudes for clusters comprising 20000 particles at a= 445, T = 0.955: theory,simulation.BF amplitudes are shown for comparison

16. CF spectral amplitudes for clusters comprising 20000 particles at a= 10, T = 0.75: theory,simulation.BF amplitudes are shown for comparison

17. CF spectral amplitudes for clusters comprising 20000 particles: theory,simulation.BF amplitudes are shown for comparison

18. CF spectral amplitudes for clusters comprising 20000 particles at a = –4.96, T = 0.75: theory,simulation.BF amplitudes are shown for comparison

19. Second spectral amplitude for clusters comprising 20000 particles as a function of a

20. Deformation parameters of clusters comprising 20000 particles , x = (c/a)2/3 – 1, at T= 0.75

21. Precursor stage of a supercritical cluster

22. Fission of a supercritical cluster

23. Ratios of the second slice spectral amplitudes calculated in three reciprocally perpendicular planes, the plane of a maximum amplitude and the planes of intermediate and minimum amplitude, as a function of time for a supercritical cluster

24. Autocorrelation function and correlation decay time for the second slice spectral amplitude for different a

25. Gibbs and Smoluchowski treated the liquid–vapor interface as a more or less abrupt change of the density and predicted that this interface is perturbed by thermal fluctuations. Mandelstam (1913) and Buff, Lovett, and Stillinger (1965) obtained where is the bare surface tension. The interface width diverges due to short-wavelength fluctuations. A simple cutoff at the interparticle distance leads to the critical point paradox. A way to overcome it is introduction of the bend rigidity (Helfrich, 1973). This yields the wave vector dependent bare surface tension where is the bend rigidity. Unfortunately, results obtained by different researchers are inconsistent. Thus, Mecke (1999) obtained a decreasing dependence ; some derived more complicated dependences.

26. Theory of cluster capillary fluctuations Probability of cluster fluctuation is defined by corresponding change in the Gibbs free energy where Assuming small fluctuation amplitudes we have derived where g0 is the bare surface tension.Based on the equipartition theorem we arrive at the amplitudes of fluctuation modes

27. Limitation of the maximum surface curvature by formation of a virtual chain

28. Formation of virtual chains limits the local curvature of the fluctuation surface: This allows one to write and to find the spectrum cutoff number and otherwise. If we introduced a common cutoff then we would arrive at failure of the capillary wave theory (critical point paradox): at sufficiently high temperature (T = 0.95), when there is no non-negative solution for g0. This difficulty is removed in proposed theory.

29. By definition, the bare surface tension g0 refers to a flat (nonperturbed) interface. Due to the parachor considerations, it depends on the surface density, which is independent on the field strength (field pressure vanishes on the surface). Therefore, g0 is field independent. The quantity is also field independent by definition. Due to the relation the ordinary surface tension g proved to be field independentas well. Bulk fluctuations Bulk fluctuations are characterized by the radial distribution of surface particles and the distribution of their number

30. 1. The case a= 0. The interface variance and proportional interface width diverge with cluster size. 2. The casea > 0 (pseudogravitation). The maximum of spectral slice amplitude is reached atkmax=(aL02/8)1/4. Divergence of interface variance at R→ ∞ is removed: In the case of gravitational attraction, the interface variance vanishes with the increase in R:

31. Theoretical CF slice spectrum for different a

32. Surface variance s2 as a function of cluster size at a = 0

33. 3. The case a < 0(Coulomb-like repuilsion).The surface variance is The maximum value a = –10 corresponds to singularity of The cluster becomes unstable with respect to fission. The classical fission threshold [Bohr and Wheeler (1939), Frenkel (1939)] supposes greater charge:

34. Conclusions • A leading order theory of surface fluctuations is proposed for clusters with a long–range particles interaction. • CF are damped by the attractive long–range interaction; the surface tension is independent of the field strength. • For the repulsive interaction, the fission threshold is defined by the bare rather than ordinary surface tension. • A nonlinear theory of large fluctuations is required.

35. Thank you for the attension! For more details, visit http://oivtran.ru/dmr