The statistical and thermodynamic description of a system of interacting particles is one of the basic problems in condensed matter physics. In this direction, significant successes have been achieved  -  . The special method used in these researches is based on a density functional theory  -  . Within the framework of this approach, the phase transitions in a system of interacting particles from a gas phase to liquid and also from liquid to solid are described  -  . Various models of the phase formation and separation in such systems are considered depending on the strength and nature of inter-particle interaction    .
The classical density-functional theory has become increasingly popular in the study of liquid-solid transition     . This is partially due to the success in providing an accurate description of the fluid to close-packed solid transition for a system of hard spheres   . On the basis of these results, it seems that the density-functional theory works well at least for systems of monatomic particles with hard interactions for which the equilibrium solid phase is close packed    . The question is whether the success of the density-functional theories with hard interactions extends to systems with softer interactions where the thermodynamically stable solid phase is not necessarily closely packed. The density-functional theory applied to nonuniform classical liquids has been able to depict a wide range of physical properties of simple solid systems  -  . It is well known that the hard sphere model is used as the “reference” system, but to make the system more realistic, we need to include some attractive interactions. Accordingly, a considerable effort has gone into developing a free energy function describing nonuniform hard sphere systems and at present there exists a quite better functional approach for these systems. Unfortunately no equation of state is known which allow for hard disk and the hard sphere fluids in the close-packing limit  . A new problem appears when the density-functional theory is applied to less packed crystal structures such as the bbc lattice, which can be at most a metastable state of the system within the density functional approach. It should be noted that the density func- tional theory gives better results for higher densities of solid state matter    . Thus to study the fluid state, it seems useful to start from the solid state density func- tional.
The problem of the solid state formation has not yet been finally solved. From the point of the density functional theory, this means that there is no common density functional that describes the system of interacting particles at all densities from gas to solid state. Instead we have one density functional well adapted for gas phase and another one for solid or liquid phase. Another problem is that the system of interacting particles crystallizes after passing through many intermediate states such as a super- cooled liquid. The description of a supercooled liquid state encounters serious dif- ficulties both in computer modelling and in theoretical approach, when it remains the problem of reduction of possible states of the system in its statistical description. We always can describe criteria of instability in a system but we can not specify a final state. A considerable fraction of the elements, when compressed, undergo a first order phase transition from a disordered state (a gas, fluid, or liquid) to a face-centered cubic (fcc) or hexagonal close packed (hcp) crystal. Computer simulations as well as experiments show that this freezing transition is already apparent in the simple model of interacting particles     . Theory can explain only the transition “afterwards”. A complete theory of hard sphere freezing should therefore identify this property of the partition function that leads to observed symmetry breaking. Thus, it seems that one obtain better description when the final state is given a priori and the conditions under which this state exists are determined from the density functional theory. In this sense the description of simple model system of interacting spheres with a finite energy of contact repulsion and long-range attraction of a gravitational type is offered. It is known that particle systems with purely attractive gravitational interaction exhibit collapse, sometimes called a zero-order phase transition   . If no short-range cutoff are introduced, then the discontinuous jump is infinite. This makes all normal noncollapsed states of the self-attractive system metastable with respect to such a collapse. On the other hand, if the same form of short-range cutoffs are introduced, the entropy and free energy jumps are finite. In this case, in the result of a collapse, the system goes into a nonsingular state with a dense core, the precise nature of which depends on the details of the short-range behavior of the potential. This model system of interacting particles can be used to describe the final state of the mutually gravitating particles with a firm nucleus or the colloidal particles in liquids and liquid crystals.
The second aspect of the solid state formation is that for systems with long-range attraction and short-range repulsion the system crystallizes passing through states with different packing (for example, cubic structure can transform to hexagonal). In the proposed approach we can find not only the criteria of formation of the particular crystal structure but we can also describe the transformations of the system from one packing to another. The phase transitions with change of packing due to the change of particle shape resulted from strong interparticle attraction (for example, to become a deformed spherical colloidal particle) are also can be described within the developed approach. It should be emphasized that the contact repulsion between the particles implies the use of the Fermi statistics even for classical particles and finite energy of contact repulsion that is because particles can not penetrate each other. This feature was not taken into account in the statistical description of phase transitions into a crystal state  .
Let us consider a model system of particles on a lattice interacting via the following potential
where the first term accounts for the long-range Newtonian attraction between particles and the second one corresponds to the short-range repulsion with being the particle radius, here is the Heaviside step function. Within the density-functional theory the free energy for this system can be written in the form  
where is the one-particle distribution function (the probability to find a particle in point) and the summation is over the lattice. We assume hard sphere model when each site of the lattice can be occupied by only one particle. This feature is accounted by Fermi type of the entropy part in (2).
By minimizing the free energy (2) we obtain some particles distribution which corresponds to some equilibrium state of the system. Depending on the tem- perature and interaction strength the distribution can be homogeneous or inhomo- geneous. If the equilibrium state is disordered then the one-particle distribution function is homogeneous, that is, where is the relative particle concentration. Any inhomogeneity in particles distribution can be described by the additional term in the expression that in the Fourier space can be expressed as, where the sum is now over re- ciprocal lattice vector. In the density-functional theories the solid density is given
a priori by a Gaussian parametrization, where the
sum is over the lattice vector and depends upon the particular crystal structure, and is the order parameter. The value of parameter being different for gas, liquid, and solid phases determines the phase of the system. In our case, where is the packing factor. It should be noticed that the multiplier is deter- mined by the particular arrangement of atoms in a unit cell of the crystal. The value corresponds to a cubic cell of linear size equal to the particle’s diameter. The small values of the order parameter correspond to the situation typical for solids, when atoms are strongly localized around the ideal lattice sites, and the free energy can be reduced to
where. Then the equation of state is
In our case, thus (4) reduces to
Note that is the whole system volume including the volume occupied by par- ticles. In the limit of low packing factor, , we obtain the equation of state of lattice gas with interaction between particles:
This equation is like a Van der Waals equation of state for hard spheres with the second part of this equation representing the interaction in the system. The formula (6) may serve as a rather good approximation for an equation of state for systems with a repulsive hard core and a newtonian attractive potentials. The thermodynamic pro- perties of elastic hard-sphere systems depend on the temperature in a trivial manner. Two methods can be used to calculate the pressure in the system. The first is based on calculating the rate from which the pressure can be deduced by the formula
, where is the low-density collision rate
for large packing factors of hard spheres gas, is the mean square velocity and is proportional to. The second method utilizes the following relation between the equation of state of hard spheres gas and the value of the radial distribution function
at:. These two methods give close values but the
second one requires some care because the radial distribution function varies rapidly in the vicinity of  . Another equation of state of the hard spheres gas,
, where a particle volume is excluded, has quit different form
though asymptotically is similar to above ones. When and thus, where is the size of the system, the collapse in system is possible. The corresponding packing factor is. The con- dition of collapse satisfy to zero value of the pressure, this yields
In the case when
Thus a collapse is possible when the mean distance between particles is less than
To study the effect of inhomogeneity of particles distribution we substitute into (2) and make use of the long-wavelength expansion of the concentration. The result is
For the hexagonal packing we have  :
and when this reduces to
If then the inhomogeneity length becomes smaller than the contact distance between particles:. In this case the attraction obviously prevails repulsion resulting in deformation of a particle with a soft nucleus. Thus the form of a cell can change from the initial cubic cell to the hexagonal. The origin of such a transformation can be explained in the following way. An attraction between particles causes a collapse in the system. The particles begin to press upon one another that results in the deformation of their shapes. A final equilibrium state in such a system is when the system volume is equal to the total volume of all the particles. But this does not change the volume of a separate particle but only changes the form of the particle. For liquid particles it is possible to estimate pressure created by a gravitational attraction and to compare it to the Laplacian pressure in a separate liquid particle. We
can estimate it from the condition, where is the surface energy and
is the number of possible nearest neighbors. The number of nearest neighbors in the close packing. When particles are deformed the distance between their centers
must be less than. Thus the condition corresponds to the condition of
deformation for such a structural unit with a soft nucleus. If then the formula (10) gives the period of particles distribution to be
. In this case the equilibrium structure of the system is
determined by only the attractive and kinetic energy and is not influenced by the short-range repulsion. In a system of emulsion particles in a fluid media with the screened Coulomb-like interaction it is possible to observe structural phase transitions with change of the shape of the droplets.
Thus the behavior of interacting particles with long-range attraction and short-range repulsion can be described in all equilibrium phases  . Starting from a gas phase the attraction results in the condensation of the system. In the condensed matter phase with lowering the temperature the attraction induces transition into solid state with the packing factor of cubic structure. At some ratio of the attraction and repulsion in- tensities there can be a new structural transition with change of the packing factor. The cubic lattice transforms into hexagonal structure. Thus the form of the structure unit can change without change of its volume, for example, sphere can transform into some dodecahedron (or disk into hexagonal cell). The possible volume of the system will be equal to the sum of individual volumes of separate structural units. Thus even a collapse of the system produced by a strong attraction between particles will be finished by the creation of a new structure with the structural factor up to unit, and this state is final for systems with strong long-range attraction and short-range repulsion of a soft type.
This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF-2014R1A1A2058029), and supported by Brain Pool Pro- gram through the Korean Federation of Science and Technology Societies (KOFST) funded by the Ministry of Science, ICT and Future Planning (161S-1-3-1617).