Bare Surface Tension and Surface Fluctuations of Clusters with Long–Range Interaction

1. Bare Surface Tension and 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. We define three particle types: internaland surfaceparticlesand virtual chains.

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

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

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

7. 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

8. 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: at sufficiently high temperature (T = 0.95), when there is no non-negative solution for g0. This difficulty is removed in proposed theory.

9. 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.

10. 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:

11. Theoretical CF slice spectrum for different a

12. 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:

13. 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

14. Simulation cell: a cluster in equilibrium vapor environment

15. We isolate the surface 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.

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

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

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

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

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

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

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. 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.

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