Abstract. In the paper, the particle dynamics in two-dimensional (2D) Yukawa systems were studied numerically. New data on a density of the thermal fluctuations of pair-interaction forces have been obtained for non-ideal systems in a wide range of their parameters. Comparisons of these thermal fluctuations with the internal energy density were performed. For the first time we have found, that for strongly correlated fluid 2D Yukawa systems the dielectric constant is inversely proportional to the second derivative of a pair potential at the mean inter-particle distance.The study of physical properties in non-ideal systems is of significant interest from the basic point of view and, at present, it is a subject of an intensive theoretical and experimental research in various areas of physics [1][2][3][4][5][6]. The main problem involved in the studies of these systems is associated with the absence of analytical theory of liquid. To predict the physical properties of non-ideal systems the various semi-empirical approaches and computer simulations with the help of numerical techniques (Monte-Carlo or molecular dynamics methods) are commonly used [1][2][3][4][5][6]. So, for an analysis of the thermodynamic properties of liquids (such, as the thermal coefficients of pressure, the heat capacities, the isothermal compressibility, etc.) the equations of state can be used [5][6][7][8]. In case of isotropic pair interactions (with the interaction energy φ(r) the physical properties of non-ideal systems, such as the energy density, and the pressure are determined by the temperature T , the concentration, n, and the pair correlation function, g(r), which can be measured experimentally or may be found from the computer simulations [5][6][7][8]. So, the internal energy density U per a particle of a system (calorific equation of state) may be written ashere m = 2, 3 is the number of dimensions in the system and n = r −m p , where r p is the mean inter-particle distance. In the scope of a classical electrostatics the energy density, U , in the nonideal quasi-equilibrium three-dimensional systems (with an accuracy of constant not dependent on T ) may be presented as [9,10]