A Power Sequence Interaction Function for Liquid Phase Particles

: In this manuscript, a function is derived that allows the interactions between the atoms/molecules in nanoparticles, nanodrops, and macroscopic liquid phases to be modeled. One goal of molecular theories is the development of expressions to predict speciﬁc physical properties of liquids for which no experimental data are available. A big limitation of reliable applications of known expressions is that they are based on the interactions between pairs of molecules. There is no reason to suppose that the energy of interaction of three or more molecules is the sum of the pairwise interaction energies alone. Here, an interaction function with the limit value w = e 2 π /e is presented, which allows for the derivation of the atomic mass unit and acts as a bridge between properties of elementary particles and emergent properties of macroscopic systems. The following properties of liquids are presented using the introduced interaction function: melting temperatures of n-alkanes, nanocrystals of polyethylene, melting temperatures of metal nanoparticles, solid–liquid phase transition temperatures for water in nanopores, critical temperatures and critical pressures of n-alkanes, vapor pressures in liquids and liquid droplets, self-diffusion coefﬁcients of compounds in liquids, binary liquid diffusion coefﬁcients, diffusion coefﬁcients in liquids at inﬁnite dilution, diffusion in polymers, and viscosities in liquids.


Introduction
Macroscopic systems of atoms and molecules exist as gases, liquids, and solids. The starting point for the discussion of gases is the completely disordered distribution of the molecules of a perfect gas without a specific volume. The starting point for the discussion of solids is the ordered structure of a perfect crystal. Liquids lie between these two extremes. The potential energy of the particles in the liquid, responsible for the attractive interaction, is of the same order of magnitude as the kinetic energy. The consequence is a mobile structure, but with a specific volume, which considerably complicates the theoretical treatment of liquids. However, it is just these properties that determine the exceptional importance of liquids to the high variety of material systems, including living matter.
Macroscopic systems have emergent properties unknown for individual particles, like temperature, as an example. Expressions for bulk properties can be obtained from statistical mechanics. With these expressions, properties of macroscopic systems can be interrelated through a common link of intermolecular force laws. The Molecular Theory of Gases and Liquids, published by Hirschfelder et al. in 1954 [1], describes the interactions between two molecules in terms of functions based on quantum mechanics.
However, emergent properties are very difficult to reduce to just interactions between pairs of particles without consideration of the other nearby particles. An attempt to consider more than one pair resulted in the Axilrod-Teller equation [2], but without significant improvement.

The Interaction Function
A first physical assumption for an interaction model is a system of uniform particles with equal finite energies. This is also a basic assumption of quantum mechanics. The N particles of the system form n groups of x 1 , x 2 , . . . , x n particles. For modeling possible interactions between these particles, the physical system is correlated with the system of natural numbers. The n groups are correlated with the n natural numbers x 1 , x 2 , . . . , x n . With this correlation a connection with the fundamental laws of arithmetic is established.
A fundamental law of natural numbers is the relation g ≤ m between the arithmetical mean m = (x 1 + x 2 + . . . + x n )/n = N/n and the geometrical mean g = (x 1 ·x 2 · . . . ·x n ) 1/n . The general theorem states that the extreme value g = m is valid only if all the x i numbers are equal. This property is used as starting point for comparisons with physical systems consisting of uniform particles. A consequence of the maximum g = m is a system of uniform groups of unities. Such uniformity of elementary particles, atoms, and molecules of a specific structure is a characteristic property of matter.
However, in addition to the extreme g = m, the following relations with further extreme values result. If the natural number N can be written in different ways as a sum of n i identical prime numbers x i as a result of g = m, then the corresponding product y i = x i ni has a maximum value of x i = 3, as shown in the following examples where n i x i = N = 30. This property of grouping natural numbers leads to the function y = x q/x with q ≥ 1 and the maximum value y max = 3 q/3 for x = 3.
The functions y max and Y max in Equation (1) represent a form of collective organization in the system of natural numbers. The fundamental additive (extensive) and multiplicative (intensive) rules in arithmetic lead to microsystems made of three particles with a maximum relative intensity of interaction represented as y max , and to macrosystems made of x >> 1 particles, with a maximum relative intensity of interaction represented as Y max in Equation (1).
Although not decisive for the following, the x i interacting particles in Equation (1) can be treated mathematically as a permutation. The process can be understood formally as x! interacting steps. Each individual step is interpreted mathematically as one change of places between two numbers. The total number of such place exchanges is x!. The relative number p x of exchanges related to x! with no particle remaining in its starting position is then [4]: p x = 1/2! − 1/3! + . . . .+ (−1) x 1/x!, where lim p x = p e = 1/e for x >> 1. (2) Fluids 2021, 6, 354 3 of 31 The limit value p e = 1/e is designated as the maximum probability of a place exchange in the system. In this model it is assumed that the interaction between uniform particles in the system occurs with the maximum probability e q·pe = e q/e . Now a second fundamental assumption for the physical system is introduced with the number 2π. The first assumption of a finite energy says nothing about the dimension of the initial particles. Assuming a spherical particle with the relative radius r = 1, a circumference of 2π is the result. In conformity with quantum mechanics, any particle traveling with a linear momentum p should have a wavelength λ given in the de Broglie relation p = h/λ, with Planck's constant h. With a relative wavelength λ = 2π, a new initial value q = 2π is introduced into Equation (1) instead of q = 1, and the following number, w = e 2π/e , is the result for the maximum value Y max in macrosystems. In this way a combination of arithmetic and geometry occurs. The limit value pe = 1/e is designated as the maximum probability of a place exchange in the system. In this model it is assumed that the interaction between uniform particles in the system occurs with the maximum probability e q•pe = e q/e . Now a second fundamental assumption for the physical system is introduced with the number 2π. The first assumption of a finite energy says nothing about the dimension of the initial particles. Assuming a spherical particle with the relative radius r = 1, a circumference of 2π is the result. In conformity with quantum mechanics, any particle traveling with a linear momentum p should have a wavelength λ given in the de Broglie relation p = h/λ, with Planck´s constant h. With a relative wavelength λ = 2π, a new initial value q = 2π is introduced into Equation (1) instead of q = 1, and the following number, w = e 2π/e , is the result for the maximum value Ymax in macrosystems. In this way a combination of arithmetic and geometry occurs. The number w represents the limit value of the following power sequence [5,6]: with the limit value: The power sequence in Equation (3) is denoted in the following as the interaction function. The limit value w represents a relative energy density in the macrosystem.
A first application of the interaction function starts with the term w1 in Equation (3). This is the first term, where n = 1 in the power sequence The number w represents the limit value of the following power sequence [5,6]: with the limit value: lim n→∞ w n,e = w = e 2π/e = 10.08909 . . . ≈ 2π w 1 − w 1,e = 10.0898 . . . , where w 1 = (1 + 2π) 1/2 = 2.6987 . . . and w 1,e = (1 + 2π) 1/e = 2.076 . . . The power sequence in Equation (3) is denoted in the following as the interaction function. The limit value w represents a relative energy density in the macrosystem.
A first application of the interaction function starts with the term w 1 in Equation (3).
This is the first term, where n = 1 in the power sequence w n = 1 + 2π n n (1+1/n) n and with the same limit value w as in Equation (3).
In conformity with the Einstein equation E = mc 2 , a body at rest has a rest energy (mass energy) of E 0 = m 0 ·c 2 , with the speed of light c.
The value w 1 represents the relative energy of one elementary particle. In conformity with Equation (1), a maximum of interaction is the result with a combination of three In order to transform the relative energy value w 1 into a quantity expressed in IS units, the limit value w in Equation (3) must be related to the decimal system with the basic number 10. In this way w 1 ·w/10 represents the energy as J.
The interaction function is of fundamental importance for emergent properties in fluid phases and the first term w 1 forms a bridge between the fundamental atomic mass unit and macroscopic properties. This connection explains the importance of the relative mass values M of particles in macroscopic systems related to the n-alkanes with i carbon atoms, as reference homologous series for all organic substances, and the term (M − 2)/14 = i = n in w n,e = 1 + 2π n n e from Equation (3).

The Interaction Function and the Optimal Entropy
The next application of the interaction function refers to emergent properties of fluids. The starting point is a system with one mol of uniform particles with a relative molar energy w 1 = (1 + 2π) 0.5 . If these particles interact in conformity with Equation (3), the limit of relative molar energy w = e 2π/e is the result for a fluid phase. This is a reference value resulting from optimal interaction of uniform particles in a liquid phase as a consequence of Equation (1) where q = 2π. However, if no interactions occur between the n particles, then w 1,e = (1 + 2π) 1/e represents the fluid phase in this situation as a limit value for n >> 1. This can be understood as the critical state of a system without interactions and holds for a liquid as well as for a gas phase, with a compression factor Z = 1. However, the term w 1 = (1 + 2π) 1/2 in Equation (3) represents the relative molar energy of a perfect gas phase. That means that a supplementary relative internal molar kinetic energy E int related to the difference w 1 − w 1,e exists in the system. This energy determines a certain temperature T.
In order to establish a connection between the difference w 1 − w 1,e and SI units, start with n = 1000 mol of hypothetical particles with the atomic mass unit u = 1.66 × 10 −27 kg in a volume V = 1 m 3 and the Avogadro constant N A .
This system, where N = n·N A = 10 3 (mol) × 6.02214 × 10 23 (mol −1 ) particles, defines the mass unit of N·u = 1.000 kg. The total energy of the system is E = ∑n i ε i , where n i particles with energy ε i and the total number of particles is N = ∑n i . Using the Boltzmann distribution, E = (N/q) ∑ε i exp(−βε i ), where β = 1/kT and the translational partition function is q [7]. The following relationship between N/q and the two terms w 1 and w 1,e is now postulated: This takes into account the relative internal molar kinetic energy w 1 − w 1,e . With the translational partition function q = N/(w 1 − w 1,e ) = V/λ 3 and λ 3 = h 3 /(2πu 0 kT) 3/2 , a value of T = 2.98058 K is the result, using the previous values for N, V, u, the Planck constant h = 6.62607015 × 10 −34 Js, the Boltzmann constant k = 1.380649 × 10 −23 JK −1 , and the thermal de Broglie wavelength λ.
The self-diffusion of particles and the entropy of a system are both a result of random particle motion. With the Sackur-Tetrode equation, the molar entropy S m of the above system can be calculated at temperature T and pressure P: where λ 3 = 1.03407 × 10 −27 for T = 2.98058 K and P = 1 Pa, and the value of S m = 108.85 JK −1 mol −1 results from Equation (5), with the gas constant R = 8.31447 JK −1 mol −1 .
The same value results for T = 298.058 K and P = 10 5 Pa. The internal energy of the gas phase at T = 2.98058 K is RT = 24.78 J mol −1 .
The consequence of the above result is an optimum value of molar entropy for fluid systems at T = 298 K and P = 1 bar. At these conditions a maximum of the variety in the structure of matter (including living matter) occurs. As an example, the many enzyme reactions in water at room temperature can be mentioned. In the range of 273 < T < 373 K and 0.7 < P < 1.7 bar, deviations of the entropy from 108 JK −1 mol −1 are < 7%.
It can be supposed that the natural evolution occurs not only in the direction of increasing entropy, but also in the direction of increasing the variety of structures. Of course, the necessary conditions for such macroscopic systems are assumed, which occur with low probability in the known world.
The liquid phase is determined in the temperature range between the freezing point T f , the boiling point T b , and the critical temperature T c with the critical pressure P c , which means T f ≤ T ≤ T c .
In conformity with classic statistical mechanics, the molar volume heat capacity of monoatomic solids is C V,m = 3R JK −1 mol −1 (Dulong-Petit rule). This rule is valid only at high temperatures, as demonstrated with quantum mechanics (Einstein and Debye), where C V,m = 0 at T = 0. The molar entropy of melting, S m,S , is for many monoatomic solids < 3R. One reason is the existence of several phase transitions in the solid state [8]. Here a molar entropy of S m,S = 3R JK −1 mol −1 is assumed at T f . For the molar entropy of liquid evaporation, a value of S m,L = R·w = 83.9 K −1 mol −1 has been defined [5] (Troutons law). That means a total molar entropy results for a gas phase, S m,G = R(3 + w) = 108.83 JK −1 mol −1 , in the above optimal situation, in conformity with Equations (4) and (5) with the value 108.85 JK −1 mol −1 .
The molar entropy change for the isothermal expansion of a perfect gas is ∆S = R·ln(V 2 /V 1 ) JK −1 mol −1 . When V 2 = 2.479 × 10 −2 m 3 mol −1 for the molar volume of the perfect gas at 298.15 K, the value V 1 = 1.03 × 10 −6 m 3 mol −1 = 1.03 cm 3 mol −1 results for the molar volume of gas in a mol of liquid, independent of its composition [8].

Modeling of Liquid-Phase Properties with the Interaction Function
All prediction equations presented in the following are based on the interaction function of Equation (3). Its application fell into two categories [5]:

1.
Prediction of physical properties of nanoparticles or drops made of a limited number n of atoms or molecules as functions of n, with an asymptotical limit value for n >> 1.
With specific constants C and corresponding dimensions with basic IS units, the energy E n = Cw n of corresponding nanoparticles and E = Cw for macrosystems is the result.

2.
Prediction of physical properties of macroscopic systems with molecules from a homologous series, like n-alkanes, with i carbon atoms. When E i = C 0 w i,e w the energy of a macrosystem can be defined that contains molecules from a homologous series with i identical atomic groups in each molecule, then the dimensionless ratio E i /E ∞ = C 0 w i,e w/C 0 ww = w i,e /w is the result.
That means, even one molecule can be treated as a system of i interacting subparticles in conformity with the interaction function defined in Equation (3). The non-polar Methylene group -CH 2 -with two strong covalent bonded H-atoms is such a subparticle, with a certain individuality along the i carbon atoms in the chain. This is in contrast to the situation in a strong polar molecule, like water, with three atoms with a relative free mobility to each other.
Such predictions can be generalized for any organic compound with the relative molecular mass M and different polarity and structure, with an additional specific parameter for that molecule. The property of chain formation is characteristic for organic compounds and the n-alkanes are therefore the reference series for all organic compounds [5,6].
The structure of this series with i carbon atoms and relative molecular masses of M i = 14·i + 2 delivers a basis for the correlation of emergent properties of fluid systems with the interaction function of Equation (3). The relative molecular mass M of a compound can be related to the number i = (M − 2)/14, corresponding to a hypothetical alkane with i carbon atoms. When i = n the interaction function w n,e is used.
Compounds of the first category above, for example monoatomic systems, are correlated with the first term w 1,e of the interaction function.
The different emergent properties of liquid systems result as exponential functions. The specific polarity and specific structure of a compound is taken into account with an increment. This increment results often in the formation of a temperature-dependent factor f = a + bT or f = a + b/T. The relative simplicity of this temperature dependence results from the exponential form of most prediction equations, with the increment f in the exponent. This increment can be included in n = (M − 2 + f )/14 or can be used as a separate parameter in the prediction equation for the corresponding property. The uniform treatment of all these properties, based on w n,e , allows for a significant reduction in the use of additional empirical factors to a minimum, because the derived prediction equations have a real connection to the corresponding properties. The two empirical terms a and b in this increment term compensate for the approximations resulting from the simplifications used. The assignment of these increments to specific functional groups and structures is important. If experimental values for a property under consideration are available for the most important functional groups and specific structures, then a map of such values may be a helpful basis for estimations of values for unknown cases.
In the following, examples of emergent properties are discussed in connection with phase transitions, like melting or freezing points (T m and T f , respectively) and vapor pressures P vp , until the critical temperature T c at equilibrium, and with dynamic processes, such as diffusion coefficients D L and viscosity coefficients η L in liquids and diffusion coefficients D P of compounds in polymers.
It is obvious that two values obtained for two temperatures for a single specific parameter, related to the relative molecular mass M, can offer only an approximation of the requested property. However, if such simple results obtained for various functional groups and specific structures are in the same order of magnitude as results from sophisticated models that are based on several empirical constants with the need for high computational effort, it is more advantageous to use the simpler way.

Modeling Melting Temperatures of n-Alkane Macrosystems Using the Interaction Function
The melting or freezing point of a pure substance is one of the most important specific macroscopic properties for its characterization and identification. It is also a specific parameter for crystalline polymers. The stability of alkane crystalline structures results from two types of interactions: lateral association of chains with interactions in the x-and y-directions and interactions in the z direction across gaps between methyl end groups [9]. The melting point of a macroscopic alkane crystal sample is the result of the different orientations of the molecules along the coordination axis, but is determined mainly by the chain length of the molecules, e.g., by the number i of the methylene groups (particles), including the two methyl end groups, because the molecular interfaces are proportional to i. One can assume that the main interactions occur along the x-and y-directions and are proportional to an interaction function w i,e . However, the interaction between the two methyl groups from two neighbor molecules along the z-direction must also be taken into account. Therefore, a constant with a smaller value than w i,e must be used for the z-axis. The ratio w 1 /w i,e , where w 1 = (1 + 2π) 1/2 , was found to be the best approximation for the relative interaction in the z-direction. The ratio of the melting point of an alkane with i C-atoms T m,i and the limit value T m,∞ for i >> 1 can then be described by the following equation: In order to determine the limit value T m,∞ , accurate melting points of known alkanes are necessary.
From the experimental melting temperatures T m,i of the n-alkanes [10] the corresponding limit values T m,∞ were calculated with Equation (6). The first few members of the homologous series did not show a regular increase in their melting points with carbon number. This irregularity was caused by crystal structural differences between molecules with even and odd i-values as well as by the methyl end groups. To avoid these effects, alkane molecules from i > 23 were selected from [10].
In addition, the melt temperature of T m,192 = 400.65 K was determined [11] at the equilibrium between the molten and crystalline state for the synthesized straight-chain alkane C 192 H 386 . From all these data the main limit value T m,∞ = T m = 415.8 K was obtained for unfolded polymethylene with a standard deviation of s = 0.35.
Shorter chains have sharp melting points, which means that the transition from a solid to liquid state occurs within a small temperature interval, e.g., approximately 0.25 K for i = 44. A slower, more gradual melting process is observed for longer chains, which can lead to large measurement errors when not adequately taken into account.
In a model [9] that incorporates the possible crystal systems into which alkanes can crystallize (hexagonal, orthorhombic, monoclinic, or triclinic), the possible interactions between all atoms of an alkane were calculated and summed up using the known distances between atoms and bond angles between H-C and C-C. In this way, the results are expressed as a function of chain length for every crystal system. Using non-linear (parabolic) curve fitting onto the experimental melting point curves in the range from 26 ≤ i ≤ 100, the limit value for an infinitely long alkane chain was found to be T m,∞ = 415.14 K.

Modeling the Melting Points of Polyethylene Nanocrystals Using the Interaction Function
Many semi-crystalline polymers, like polyethylene (PE), form lamellae crystals that are 10-30 nm thick and are at least one order of magnitude larger in the lateral direction [13]. For the melting temperature T 0 m,i of an alkane crystal, with only one lamella made of the molecular chain with i carbon atoms but an unlimited lateral dimension, Equation (6) can be adapted in the following manner: No interaction between layers in the z-direction occurred in this case, and for i >> 1 and a macroscopic surface, an upper limit melting temperature T m 0 = 415.8·(2w/(2w + w 1 /w)) = 410.36 K was the result. For a microscopic monolayer lamella made of molecule chains with i Catoms and where j molecules are arranged in each lateral direction x and y, Equation (7) is used: As an example, a microcrystal lamella of 300 nm × 300 nm × 30 nm was considered, which is typical for PE. With the orthorhombic C-C lattice distance l c = 0.1273 nm in the z-direction [13] 30/l c = 236 = i and approximately j = 2400 molecules arranged in the x-and y-direction, T 0 m,i = 415.8·(2w i,e /(2w + w 1 /w)) = 398.14 K was the result with Equation (7) and T 0 m,j,i = 398.14·(w j,e /w) = 396.94 K with Equation (8). The measured melting temperatures collected from different sources for ultra-long n-alkanes [13] were between 395 and 404 K.
When alkane chains become long enough (i > 200) they start to fold. The layer thickness d of a lamella is then no longer proportional to the total chain length but rather proportional to the length of the chain segment between a CH 3 end group and a folding point or the segment between two folding points. Nevertheless, the contributions to the melting temperature in the x-and y-directions are assumed to depend on the whole chain, which means on the contribution of the w i,e value in Equation (7). This is in accordance with the situation in the critical state, with no distinctive direction in space [5]. The increase in the critical state temperature values T c,i occurs only in conformity with the w i,e -values. It is only important that the primary structure of the non-branched alkane chains with the covalent C-C bonds represent a linear arrangement. This is in accordance also with the finding that cyclic alkanes can have the same melting point as the linear alkanes with the same number of C atoms (i > 140) but with a thickness of the lamellae in the stack that is only half of that of the linear alkanes [14].
From a collection of available experimental data [13] on the melting temperatures of the ultra-long alkanes that exist in extended chains as well as folded-chain crystals with up to i = 390 carbon atoms, the melting temperatures are functions of i and not functions of the lamella thickness. This is also in line with the dynamic explanation of the melting process, where the mobility of the molecule is the dominant factor for the breakdown of the lattice [13]. However, the mobility depends on the number i of repeating units within the chain molecule. The melting transition is not as sharp as for pure one-component systems and covers a certain temperature range. In folded alkanes, it can also lie below the corresponding i value of the alkane. One explanation may be the fact that after folding the whole surface of the chain it is not available for interaction with the neighbor molecules, because the neighboring folds neutralize each other.
The limit value T m,∞ = 415.8 K and the molecular specific lattice distance l c are the only data needed in addition to the relative molecular mass M if the melting temperature for a nanocrystal with defined dimensions is calculated.

Modeling the Melting Temperatures of Metal Nanoparticles Using the Interaction Function
A solid rectangular material composed of similar n x ·n y ·n z = N spherical particles was assumed. The melting temperature of the material is the result of interactions between the particles arranged along the three-coordination axis. The relation T m,N /T m = (w nx,e + w ny,e + w nz,e )/(w x + w y + w z ) was Fluids 2021, 6, 354 9 of 31 considered in conformity with the interaction function between the melting temperature T m,N of the material with N particles and the melting temperature T m of a macroscopic sample. For a cube with n x = n y = n z , T m,N /T m = w n,e /w, where n = N 1/3 . That means a decrease in the melting temperature of material clusters occurs with a reduction in their size.
Various structures of nanoparticles are possible, depending on their mode of preparation. In a solid spherical cluster composed of approximately the same number N of spherical particles as in a cube, N 1/3 particles are arranged along the diameter d and a similar relation, T m,N /T m = w n,e /w as above for a cube, can be assumed. With the diameter d 0 of a particle, n = d/d 0 .
In a simple three-dimensional arrangement of particles, each atom is surrounded by six atoms at equal distances. However, spherical particles (metal atoms) crystallize preferably into one of the densest spherical packings, e.g., gold atoms with a face-centeredcubic internal structure, fcc (face centered cubic) [15]. In such a system, each atom is surrounded by 12 atoms at equal distances. The different arrangement possibilities in a crystal produce different coordination numbers. Depending on the number of direct neighbors, more than one of these can be affected by the interaction along one of the three coordination axes. The coordination number is taken into account by multiplying the particle number N 1/3 along the cluster diameter with a factor, for example a = 2 for the coordination number 12. The simplest way of modeling the melting temperature T m,d of a metal cluster with diameter d is to correlate the melting ratio T m,d /T m with the interaction function as T m,d /T m = w n,e /w, where n = a·d/d 0 and T m of the bulk. However, a consequence of the small radius of nanoparticles is the interaction of the surface atoms attracting inwards only, for example in the x-direction. Taking this into account for three-dimensional clusters with three coordination axes, w nx,e /3w < w n,e /w is the result for the surface layer, because w ny,e = w nz,e = 0. That means that a surface layer with the thickness of two atoms can be assumed to be liquid below the melting temperature T m,d .
If d − 4d 0 is used instead of d, the following equation for the melting temperature of a cluster is defined using the interaction function: The first suitable experimental data [16] to test this comes from the melting of gold clusters ranging in diameter from 2 to 24 nm. The number N of atoms in such largesized clusters lies between 300 and over 1/2 million. Even for such clusters, the melting points are significantly lower than the melting point of macroscopic gold samples where T m = 1336 K. In Figure 2, the values of measured melting-point data obtained for different cluster diameters d [16] are presented in comparison with the T m,d values calculated with Equation (9). For Au, the metal atom diameter [15] d 0 = 0.2884 nm for a coordination number 12 and a = 2 were used.
The pre-melting effect resulting from weaker bonded surface atoms has been mentioned in several studies [16][17][18]. Molecular dynamics were used for simulations of the melting temperature.
Further examples with reliable experimental data obtained from indium clusters [19], Pb [20,21], Sn [22], and Al [23] are in agreement with the above results for Au obtained with the interaction model.
For modeling the melting temperature of nanoparticles, molecular dynamic simulations [24] and analytical approaches [25] are described.
Equation (9) can be extended for melting temperatures T mw,d of nanowires with diameter d and T mf,d of nanofilms with thickness d. For nanowires with film thickness d. Equations (10) and (11) were verified with experimental data obtained for indium nanowires and nanofilms [26], using a = 2 like for fcc structures. The pre-melting effect resulting from weaker bonded surface atoms has been mentioned in several studies [16][17][18]. Molecular dynamics were used for simulations of the melting temperature.
Further examples with reliable experimental data obtained from indium clusters [19], Pb [20,21], Sn [22], and Al [23] are in agreement with the above results for Au obtained with the interaction model.
For modeling the melting temperature of nanoparticles, molecular dynamic simulations [24] and analytical approaches [25] are described.
Equation (9)  with n = aꞏ (11) with film thickness d. Equations (10) and (11) were verified with experimental data obtained for indium nanowires and nanofilms [26], using a = 2 like for fcc structures. All published models for the melting temperatures of nanoparticles are more or less based on simplified assumptions. The real properties of the clusters considered depend on many conditions found during their preparation and different modes of measurement [16]. The crystalline structure of very small clusters can be changed, for example, during the preparation. The pre-melting surface region of not-perfectly-spherical particles may be of different thicknesses at different positions, but a mean liquid layer of 1-2 atom diameters is a reliable assumption. All published models for the melting temperatures of nanoparticles are more or less based on simplified assumptions. The real properties of the clusters considered depend on many conditions found during their preparation and different modes of measurement [16]. The crystalline structure of very small clusters can be changed, for example, during the preparation. The pre-melting surface region of not-perfectly-spherical particles may be of different thicknesses at different positions, but a mean liquid layer of 1-2 atom diameters is a reliable assumption.
Here, only T m , the melting temperature of the bulk material, the atom diameter d 0 , and the corresponding coordination number are needed for a given nanoparticle with diameter d. In all published models more specific parameters are involved.

Modeling the Solid-Liquid Phase Transition Temperatures for Water in Nanopores Using the Interaction Function
The melting of ice and freezing of water in a series of mesoporous silica materials with pore diameters from 2.9 to 11.7 nm was studied using differential scanning calorimetry [27,28]. A lowering of the melting temperature up to 50 K was found for water in pores of radius r nm and was described with the equation T m − T m,r = C/(r − t), where C = 52.4 K nm, obtained from fitting experimental data and t = 0.6 nm. T m and T m,r are the bulk and pore phase transition temperatures, respectively. From the experimental investigation [27], a surface layer of non-frozen water needed to be considered in all pores, which reduced the effective pore radius from r to r − t. The thickness of this layer, t = 0.6 nm, corresponded to approximately two water molecules.
One cm 3 of H 2 O with a density of 0.93 ± 0.03 g cm −3 in the pores [27] contains N = 0.93 × 6.022 × 10 23 /18 = 3.111 × 10 22 water molecules, which means N 1/3 × 10 −7 = 3.145 molecules per nm, with a molecule diameter d 0 = 0.318 nm. With two molecules, t = 0.636 nm was the result. In contrast to the methylene groups in the n-alkanes with two non-polar covalent bonds between hydrogen and carbon, water, with its strong hydrogen bonds, was treated here in a first approximation as a system where each molecule had three interacting subparticles * in the crystal. Consequently, in a pore with the diameter d nm and the pore axis along the z-axis, a number n = 9.435(d − 2t) particles was assumed to interact along the pore diameter in the xand y-axis, respectively. Along the z-axis in a filled pore an unlimited number of N particles can be assumed. With these assumptions and taking into account the liquid layer of 0.636 nm, the following equation for the dimensionless ratio between the melting point T m,d and T m = 273.15 K was the result: This equation has the same structure as Equation (10) for nanowires. An explanation of the liquid layer t is similar to that of the metal nanoparticles above.
In Figure 3, the experimental values T m,d /T m [28] are compared to the theoretical curve obtained with Equation (12).  From the experimental results [27], a limiting pore diameter of 2.7 nm was determined. With this value and Equation (12), where n = 3•3.145•(d − 2•0.636), a maximum melting point depression of 61.76 K was the result, compared to the published values 61.6 K and 63 K [29]. This supports the validation of Equation (12).
If wn,e/w is used for nanoparticles of ice instead of Equation (12) similar to nanowires, a maximum depression of 92.6 K would result. This value is important for modeling vapor pressures of water nanodroplets.
In this approach only the melting temperature of ice and the molecular diameter d0 of water are needed.
The remarkable importance of the above study with water in nanopores lies in the homogeneous structures of the used mesoporous silica materials with constant diameters. This allows for precise measuring of the melting temperatures as function of the pore diameter and extrapolations to minimum values.

Modeling Vapor Pressures in Liquids Using the Interaction Function
An equation for vapor pressures of pure organic fluids and partial pressures in multicomponent systems was recently published [5]. A remarkable conclusion resulted from a comparison of the melting of metal clusters to the melting of ice in pores. In both cases, independent observations came to the necessary assumption of a thin liquid skin that surrounds the solid core. The liquid layer of water has been experimentally proven [29] and can be described as a uniform and continuous liquid mantel at the inner pore walls. For metal clusters with often different geometrical structures, the behavior of the liquid skin is less well established. In the case of metal clusters this skin is formed on the free surface, whereas the liquid layer in the pore is between the pore wall and the ice phase.
From the experimental results [27], a limiting pore diameter of 2.7 nm was determined. With this value and Equation (12), where n = 3·3.145·(d − 2·0.636), a maximum melting point depression of 61.76 K was the result, compared to the published values 61.6 K and 63 K [29]. This supports the validation of Equation (12).
If w n,e /w is used for nanoparticles of ice instead of Equation (12) similar to nanowires, a maximum depression of 92.6 K would result. This value is important for modeling vapor pressures of water nanodroplets.
In this approach only the melting temperature of ice and the molecular diameter d 0 of water are needed.
The remarkable importance of the above study with water in nanopores lies in the homogeneous structures of the used mesoporous silica materials with constant diameters. This allows for precise measuring of the melting temperatures as function of the pore diameter and extrapolations to minimum values.

Modeling Vapor Pressures in Liquids Using the Interaction Function
An equation for vapor pressures of pure organic fluids and partial pressures in multicomponent systems was recently published [5].
The exact correlation between the vapor pressure P vp of a liquid and the temperature T at equilibrium needs an equation of state that includes the gas and the liquid phase. The lack of exact equations of state leads to the critical state as a starting basis. Due to the large range of molecular masses, the model is developed in two steps: first for monoatomic systems and for the first small molecules in a homologous series, and second for homologous series beyond the first molecules.
The following critical reference state equation [5] for small molecules resulted: where the critical pressure P c,r = 2.375 × 10 6 Pa, the critical molar volume V c,r = 4.153 × 10 −5 m 3 mol −1 , the critical compression factor Z c,r = w 1 /w = 0.2675, and the critical temperature T c,r = 44.35 K, with the gas constant R = 8.31447·JK −1 mol −1 .
No additional specific constants are needed, only the number i of carbon atoms. The general vapor equation for organic fluids based on Equation (13) is: where the critical temperature T c,k = (w k,e /w)·T c,∞ of a hypothetical alkane with the relative molecular mass M + U and the limit value T c,∞ = 1036.5 K for polymethylene [5].
The general application of the vapor equation for organic fluids lies in the introduction of a dimensionless structure increment U with the same unit as the dimensionless molecular mass M. This increment U takes into account the polarity and structure difference of the compound in comparison to an n-alkane with the same molecular mass M. In this way a parameter k = (M − 2 + U)/14 is defined [5] and introduced into the interaction function w k,e = (1 + 2π/k) k/e instead of w i,e . This procedure is also used in the treatment of vapor pressures in liquid droplets.
In addition to the U-values for organic compounds, the U-values for water are predicted, because water stands for the compound with the upper limit of polar increments. With the two values of u 1 = 75.71 for T 1 = 298.15 K and u 2 = 68.06 for T 2 = T b = 373.14 K, the increment for water is U = U b = u 1 + [(u 1 − u 2 )/(T 1 − T 2 )]·(T − T 1 ) = 75.71 − 0.102·(T − 298.15). For water, the value M = 2 × 18.015 = 36.03 is used because Equation (14) does not give correct results for low M values. This example with water is given because it is used later for the prediction of vapor pressure in water nanodroplets.
For binary vapor-liquid equilibria between components 1 and 2, two increment differences, ∆U1,2 and ∆U2,1, are introduced into Equation (15) The increment differences are related to increments U1 and U2 and contain the additional empirical fractions δ1,2 and δ2,1. As an example, Figure 4 shows the vapor mole fractions Y 1 for heptane and Y 2 for toluene at the corresponding liquid mole fractions X 1 in the binary system heptane(1)-toluene (2)    With known specific fractions of the two U-values in binary systems, the vaporliquid equilibrium behavior of multicomponent systems can be predicted without additional empirical data as a good first approximation [5].
A big advantage of the above approach in comparison to other methods is the easy direct measurement of partial pressures with head space gas chromatography (HSGC) and the determination of the fraction values δ, with a few mixtures and molar fractions X1 between 0 and 1.

Modeling Vapor Pressure of Liquid Droplets with the Interaction Function
The vapor pressure in Equation (14) defined in [5] refers to macroscopic phases with n >> 1 molecules. For small liquid amounts, as in the case of droplets, smaller n values appeared and a decrease in interaction energy between the molecules occurred, which led to a higher vapor pressure.
In contrast to the solid phase, particles are mobile in all directions, and it was a priori difficult to decide the exact modus of the interactions. The ratio P/P * between the vapor pressure P of a liquid droplet with n molecules and the vapor pressure P * of the macroscopic liquid phase can be assumed, for example in a first approximation, to be propor- With known specific fractions of the two U-values in binary systems, the vapor-liquid equilibrium behavior of multicomponent systems can be predicted without additional empirical data as a good first approximation [5].
A big advantage of the above approach in comparison to other methods is the easy direct measurement of partial pressures with head space gas chromatography (HSGC) and the determination of the fraction values δ, with a few mixtures and molar fractions X 1 between 0 and 1.

Modeling Vapor Pressure of Liquid Droplets with the Interaction Function
The vapor pressure in Equation (14) defined in [5] refers to macroscopic phases with n >> 1 molecules. For small liquid amounts, as in the case of droplets, smaller n values appeared and a decrease in interaction energy between the molecules occurred, which led to a higher vapor pressure.
In contrast to the solid phase, particles are mobile in all directions, and it was a priori difficult to decide the exact modus of the interactions. The ratio P/P* between the vapor pressure P of a liquid droplet with n molecules and the vapor pressure P* of the macroscopic liquid phase can be assumed, for example in a first approximation, to be proportional to the ratio e −wn,e /e −w = exp(w − w n,e ), where w n,e = (1 + 2π/n) n/e from Equation (3). However, the correlation of the interaction function w n,e with the homologous series of n-alkanes with i carbon atoms suggests a stabilizing preference for interactions between the i particles along a string. This was a consequence of electromagnetic interaction.
If the attractive interaction between the particles in the three-dimensional liquid droplet, with mobility between the particles in all directions, occurred along a linear arrangement, then a maximum value was obtained only for spherical droplets. In a spherical droplet i = n 1/3 molecules from the n molecules of the droplet then interacted along the diameter. However, at the same time, all n molecules of the droplet could exchange their relative places according to their kinetic energy, in conformity with the structure of e 1/e from Equation (1), and in this way the droplet was a system of dynamically connected particles.
The relative interaction energy of a particle in the liquid phase is temperature dependent and assumed to be w k,e − w k,e (T m ·f v /T) (17) where w k,e = (1 + 2π/k) k/e and k from Equation (15). T m represents the melting/freezing temperature of the macroscopic liquid phase. At T < T m , no liquid phase exists. The factor f v in Equation (17) was the sole additional temperature-dependent specific parameter for the liquid stage. This led along the circumference 2π and with i particles along the droplet diameter to an amount defined as a k,i = 2π·(w k,e − w k,e ·(T m ·f v /T))·i (18) and the final relation P/P * = exp(w − w Ai,e ) with w Ai,e = 1 + 2π a k,i a k,i e (19) was the result. Checking the validity of the above relation started with mercury droplets. Liquid mercury with monoatomic particles represents the simplest case, for which w 1,e from Equation (3) instead of w k,e and its Tm = 234.316 K are used in Equation (18): The relation w − w Ai,e from Equation (19) can be compared to the results obtained with the well-known Kelvin Equation: where the surface tension is γ, V m is the molar volume of the liquid, and the droplet radius is r.
The number n of particles in a spherical droplet resulted from the following relation: where ρ (g cm −3 ) represents the liquid density at temperature T, M (g mol −1 ) is the molecular or atomic mass, r (nm) is the droplet radius, and N A = 6.022·10 23 (mol −1 ) is the Avogadro number. The number of atoms along a droplet diameter n 1/3 = i is listed in Table 2, for r = 10 nm. The values for the temperature-dependent factor f v were calculated using the assumed equation where a v and b v are fitting parameters from the experimental data. When f v = 0.295 + 0.002573·T the main error between the results obtained with Equations (19) and (20) was <1% for r = 10 nm ( Figure 5). This f v resulted from fitting the corresponding values for the five selected temperatures. Neglecting the factor f v = 1.062 needed for Hg at 25 • C, an error of e = −22% resulted for ln(P/P*) at a droplet radius of r = 10 nm.
where av and bv are fitting parameters from the experimental data. When fv = 0.295 + 0.002573•T the main error between the results obtained with Equations (19) and Equation (20) was < 1% for r = 10 nm ( Figure 5). This fv resulted from fitting the corresponding values for the five selected temperatures. Neglecting the factor fv = 1.062 needed for Hg at 25 °C, an error of e = −22% resulted for ln(P/P*) at a droplet radius of r = 10 nm. At r < 2 nm the results of the two methods deviated remarkably. At very small droplets, with r < 1 nm, the Kelvin equation was affected by thermodynamic limitations, as shown below.
As the next example, the n-alkane tetradecane was used, representing the class of compounds with no polarities. The following fv-value as a function of T K can be used in Equation (18)  At r < 2 nm the results of the two methods deviated remarkably. At very small droplets, with r < 1 nm, the Kelvin equation was affected by thermodynamic limitations, as shown below.
As the next example, the n-alkane tetradecane was used, representing the class of compounds with no polarities. The following f v -value as a function of T K can be used in Equation (18), with k = (M − 2)/14 for w k,e : f v = 0.41 + 0.001714T.

·100
The experimental measurement of vapor pressures for small liquid drops is extremely difficult. Of special interest in this connection is water, due to its importance in the form of fog, as an example. In a recently published study [30], a grand canonical screening approach was used to compute the vapor pressures of water nanodroplets from molecular dynamic simulations. This method demonstrated that the applicability of the Kelvin equation extends only down to small lengths of 1 nm. However, for clusters of 40 particles or less, the macroscopic thermodynamics and the molecular descriptions deviate from each other.
Equation (19) was not affected by thermodynamic limitations, as can be seen in the comparison of predicted vapor pressure of water nanodroplets to the results obtained with the grand canonical screening approach [30].
Relative vapor pressures P/P* for eight different water droplets composed by n molecules and the corresponding radius r at 278, 298, and 318 K are given in [30]. For 298 K these values were: For the Kelvin Equation (20), the surface tensions γ = 0.07199 (N/m) and the molar volumes V m = 1.807 × 10 −5 (m 3 /mol) for water at 298 K.
The value k = (M − 2 + U)/14 = 7.839 for water in Equation (15) resulted from [5], where U = U b for water, as shown above, and the corresponding w k,e = 5.46 for Equation (18) at 298 K. Instead of the product T m ·f v with the empirical factor f v , the minimum melting temperature, T m,min for ice nanoparticles was used in Equation (18). The value T m,min resulted from the above treatment of ice in nanopores, as in Section 4.4. Instead of (2w n,e + w)/3w in Equation (12), the ratio w n,e /w was used for ice nanoparticles and This result emphasizes the importance of reliable values obtained for melting temperatures of nanoparticles. No empirical f v -value is needed. Figure 6 shows the comparison of the logarithm of the relative vapor pressure as a function of the inverse radius at 298 K of the values from [30], and of the values obtained with Equation (19) and with the Kelvin Equation (20).
This result can be considered a validation of the interaction function as the basis for a correct model for vapor pressure predictions in liquid droplets. It explains the deviation between the results from Equation (19) and the Kelvin Equation (20), as shown in Figure 5 for Hg. This result can be considered a validation of the interaction function as the basis for a correct model for vapor pressure predictions in liquid droplets. It explains the deviation between the results from Equation (19) and the Kelvin Equation (20), as shown in Figure  5 for Hg.
Due to the correct applicability of the Kelvin Equation at r = 10 nm, a comparison of results obtained with Equation (20) and Equation (19) for w − wA,i at r = 10 nm allowed for the easy determination of the fv values for different structures from tabulated experimental surface tensions [10]. Essentially all estimation techniques for the surface tension were empirical [3] and it is therefore recommended to use the collection of experimental data as given in [10] or [12]. From these data a corresponding collection of fv-values as functions of specific structures can be derived, for which the corresponding U-values [5] are predictable. This makes prediction in the absence of experimental values possible, down to droplets with ≈ 8 particles.

Modeling Diffusion Constants in Liquids Using the Interaction Function
The interaction function can also be considered a starting point for dynamic processes like diffusion occurring in macroscopic systems. Although transport phenomena based on diffusion are determined from concentration gradients and therefore not equilibrium processes, self-diffusion also takes place at equilibrium. This means that the diffusion coefficient can be described as an exponential function of the interaction function in Equation (3) [6].
In conformity with the three fundamentally different states of aggregation, three different procedures for predicting diffusion coefficients were used. Here, the IS unit conversion factor Du = m 2 s −1 = 10 4 cm 2 s −1 was used. Diffusion in liquids has few similarities with diffusion in gases and solids. Due to the correct applicability of the Kelvin Equation at r = 10 nm, a comparison of results obtained with Equations (19) and (20) for w − w A,i at r = 10 nm allowed for the easy determination of the f v values for different structures from tabulated experimental surface tensions [10]. Essentially all estimation techniques for the surface tension were empirical [3] and it is therefore recommended to use the collection of experimental data as given in [10] or [12]. From these data a corresponding collection of f v -values as functions of specific structures can be derived, for which the corresponding U-values [5] are predictable. This makes prediction in the absence of experimental values possible, down to droplets with ≈8 particles.

Modeling Diffusion Constants in Liquids Using the Interaction Function
The interaction function can also be considered a starting point for dynamic processes like diffusion occurring in macroscopic systems. Although transport phenomena based on diffusion are determined from concentration gradients and therefore not equilibrium processes, self-diffusion also takes place at equilibrium. This means that the diffusion coefficient can be described as an exponential function of the interaction function in Equation (3) [6].
In conformity with the three fundamentally different states of aggregation, three different procedures for predicting diffusion coefficients were used. Here, the IS unit conversion factor Du = m 2 s −1 = 10 4 cm 2 s −1 was used. Diffusion in liquids has few similarities with diffusion in gases and solids.

Diffusion in Gases
At constant temperatures, the ratio D G,2 /D G,1 of the self-diffusion coefficients in an ideal gas at two states 1 and 2 equals the ratio V 2 /V 1 of the system volumes at these states. When V 1 = RT r /p u and V 2 = RT/p, a starting point for modeling diffusion coefficients is the relation [6] where w 1 is from Equation (3) for free moving monoatomic gases, p u = 1 Pa, and T r = 298.15 K.
At p = 1 bar = 10 5 Pa and T = 273 K, for example, D G = 1.36 cm 2 s −1 was the result with Equation (23), compared to 1.4 cm 2 s −1 measured for He [31]. At these conditions the gas phase can be considered ideal, with no further interactions between the gas particles. In this investigation no further treatment for the gas phase followed, since many accurate methods are already explained in the literature [32].

Diffusion in Solids
The other extreme case of diffusion is represented by a perfect solid state, in which no free moving particles exist. With the IS unit Du, one can assume the following expression for self-diffusion coefficients of atoms in a solid metal: The first term in the bracket represents the activation energy of diffusion as an interaction between a relative cross-sectional area of a monoatomic particle in a condensed state and the solid matrix with its melting temperature T m . The second term is a further brake on the diffusion of an atom in a solid matrix. Whereas the diffusing particles in the gas phase move through a matrix in which all particles are in the same degree of movement, the diffusing atoms in one solid metal move through a rigid matrix made of immovable particles. This different situation is taken into account with the limit value w.
A difference of 20 orders of magnitude between the diffusion constants in a solid matrix and a gas phase, as shown in the above examples, emphasizes the significant difference between these extremes of aggregate state [6].
The negative terms in the exponent of Equation (24) demonstrate the strong interaction between the particles in a solid state. A similar effect is the relative rest energy, E 1 3 defined above in Section 2 (The Interaction Function) as resulting from the three interacting particles w 1 .

Diffusion in Liquids
Situated between a state with completely free movement of the particles and a state with a practical absence of any movement, liquids represent extremely complex behavior. Not only do single particles in the liquid phase move through the macroscopic matrix, but more or less microscopic and macroscopic parts of the liquid itself also move through the system. An additional effect is the overlapping of the linear occurring interactions in conformity with electromagnetic interactions and the kinetic displacement of particle groups in the liquid. This may initiate turbulence in the liquid via twist effects.
Liquid state theories for calculating diffusion coefficients are quite idealized and none are satisfactory in providing relations for calculating D L [3].
Here the following expression for a self-diffusion coefficient of molecules in liquids was assumed: where: D L,r = w w 1 − w 1,e exp w 1,e − (w 1,e ·w) The reference factor D L,r refers to the relative kinetic energy of the liquid to the ratio between w and the initial free kinetic energy w 1 − w 1,e of the system (see above, Section 4). The expression in the bracket has the following meaning: As mentioned above, microscopic parts of the liquid moving through the matrix are represented by w 1,e . They are hindered by interaction with the matrix and the resistance from the matrix w 1,e − (w 1,e ·w) 1 2 − w. At the same time, macroscopic parts from the liquid, represented by w, are moving and are hindered from interacting with the surrounding liquid and its resistance w − (w·w) 1/2 − w = −w. In the above relation a simplified contribution of these two extremes is used.
The exponent in Equation (25) contains the specific diffusing particle, with its relative molecular mass M, represented as w n,e = (1 + 2π/n) n/e , where n = (M − 2)/14. Its freezing temperature is approximated by (w n,e /w)T m,∞ , with the asymptotic limit temperature T m,∞ = 416 K of polymethylene. This approximation represents a simplified form of Equation (6). The only empirical parameter is a fine-tuning constant, f D = a D + b D ·T, as a linear function of temperature, close to 1. This constant also compensates for the above simplifications in the expression.

Self-Diffusion Coefficients
Only a limited number of experimental data for self-diffusion coefficients are known. The validation procedure of the above relation for D L starts with the homologous series of n-alkanes.
For n-alkanes [34,35], the following values for the factor f D = a D + b D T in Equation (25) The factors f D1 and f D2 in Equation (27) were obtained by fitting experimental data in [34,35].
For a compound C n H 2n+2 only the number n of carbon atoms and the temperature were needed.
For heptane the percentage average absolute deviation AAD% = 6.3 for 7 D L -values in the temperature range 190 < K < 377 [35]. For dodecane and tetradecane the corresponding errors were AAD% < 5 [34]. The error (e%) between experimental and predicted D L -values was 0.1 for octadecane at 50 • C and −16 for C 32 H 66 at 100 • C [35]. Table 3 contains values for f D and AAD% values for 13 compounds at different temperatures. An increase in the a D -values and decrease in the b D -values resulted in increasing polarities of the compounds (Figure 7). This behavior allowed for a rough estimation of the D L -values for compounds with unknown f D -values. With more experimental data, closer correlations between f D -values and functional groups are possible. A rough estimation of self-diffusion values is also possible from viscosity coefficients, as shown below in Section 4.8. An increase in the aD-values and decrease in the bD-values resulted in increasing polarities of the compounds (Figure 7). This behavior allowed for a rough estimation of the DL-values for compounds with unknown fD-values. With more experimental data, closer correlations between fD-values and functional groups are possible. A rough estimation of self-diffusion values is also possible from viscosity coefficients, as shown below in Section 4.8.

Binary Liquid Diffusion Coefficients
A binary mixture of n-hexane (A) and n-dodecane (B) was considered as an additional example of liquid diffusion. The mutual diffusion of these two hydrocarbons increases as the mixture becomes richer in A. If the mole fraction X A → 1.0, the mutual diffusion coefficient D AB = D BA → D • BA , where this notation signifies that the limiting diffusivity represents the diffusion of B in a medium consisting essentially of A [3]. Similarly, D • AB is the diffusivity of A in essentially pure B. Self-diffusion coefficients differ from binarydiffusion coefficients, and there is no way to relate the two coefficients directly, without an additional specific parameter [3].
Liquid state theories for calculating binary liquid diffusion coefficients are not satisfactory in providing relations for calculating the mutual diffusion coefficient D AB . Different models are presented in [3]. Some models are based on the Stokes-Einstein equation, which strictly applies to macroscopic systems with large spherical molecules diffusing in a dilute solution.
The following equation is proposed for the mutual diffusion coefficient in a binary mixture of compounds with the relative molecular masses M A and M B as function of the mole fraction X B : f DA and f DB are for the pure compounds A and B, respectively, introduced above for the self-diffusion coefficients. δ D is an interaction factor of the two compounds in the binary system that is obtained from a comparison of the calculated data using Equation (28) with a few experimental points.
If A = B, w B,e = w A,e , f DA = f DB = f D , X A = X B = 1, and δ D = 0. That means that Equation (28) leads to Equation (25). For the binary system hexane(A)-dodecane(B), δ D = 0.86·X A + 0.105. The AAD% = 1.6 for D AB at 298.15 K [35] (Figure 8). These AAD% values from comparison with experimental data in [35] are significantly smaller than the corresponding predicted data obtained from [36,37].  [35] are significantly smaller than the corresponding predicted data obtained from [36,37] In addition to the molecular weights of solvent B and solute A and factor f DB = a DB + b DB T for the self-diffusion of B, factor ϕ B = c + d·w A,e is needed to predict the whole series of D • AB data in one solvent B. Factor ϕ B with its linear dependence from w A,e resulted from comparisons of data calculated with Equation (29) along with a few experimental data and allowed for the easy determination of c and d from two known D • AB -values.
In addition to the molecular weights of solvent B and solute A and factor fDB = aDB + bDBT for the self-diffusion of B, factor φB = c + d•wA,e is needed to predict the whole series of D°AB data in one solvent B. Factor φB with its linear dependence from wA,e resulted from comparisons of data calculated with Equation (29) along with a few experimental data and allowed for the easy determination of c and d from two known D°AB-values.
In Table 4 the experimental D°AB values for 6 compounds in decane are compared with the predicted values with equation (29).  In Table 4 the experimental D • AB values for 6 compounds in decane are compared with the predicted values with Equation (29).   Table 4 at 25 • C in n-hexane, n-heptane, n-octane, and n-decane were 2.3, 2.9, 2.1, and 0.84, respectively, in comparison to the experimental values [10].
In Table 5 the D • AB values for 26 compounds in water and the errors from comparison with predicted values obtained with Equation (29) are listed.

Modeling Viscosities of Liquids Using the Interaction Function
There is no satisfactory theoretical basis available yet for the estimation of liquid viscosities [3]. However, a strong connection between viscosity and liquid diffusion coefficients exists [1]. In this model, the following equation for the prediction of the liquid viscosity coefficient η L is proposed: where: η r = 1 10 8 D L,r = 4.3621 × 10 −4 Pa s and w n,e = (1 + 2π/n) n/e , where n = (M n − 2)/14. The viscosity reference factor η r results from the inverse value of the self-diffusion reference factor D L,r in Equation (26).
In the following, empirical factors f η = a η + b η T in Equation (30) obtained from fitting experimental data for a series of compounds and homologous series are given and the calculated viscosity coefficients are compared to the literature data [10] in the temperature The AAD% values of the viscosities for liquid n-alkanes with n carbon atoms from hexane to hexadecane, calculated with Equations (30) and (31) in the temperature range of −25 < t • C < 100, were < 4. For pentane AAD% = 12.6 and for octadecane, AAD% = 9.
The corresponding factors fη = a η + b η ·T for several other compounds are shown in Table 6. Table 6. The a η and b η values for the compounds with relative molecular masses M and the AAD% from a comparison of the calculated viscosity coefficients to the literature data [10] for N measurements in the temperature range of −25 < t < 100 • C. The fη values for alkylbenzenes, 2-ketones, and ethyl acetate between 25-50 • C were between 1.05 and 1.15. An increase in the a η -values and decrease in the b η -values resulted in increasing polarities of the compounds (Figure 10). This behavior is similar to that of the self-diffusion coefficients (Figure 7).
The Θ-values lie near 1 and for many organic compounds one mean value Θ = 1.25 can be used. A few examples of the AAD% values for DL estimated with Equation (32) using the corresponding ηL data are shown in Table 7. There are many experimental data for ηL values available in the literature for different structural groups [12]. From such collections the corresponding terms aη and bη can be easily determined. From aη and bη values for different structures, as exemplified in the above tables, estimates of DL-values are possible using Equation (32), for which no experimental data are available. An important relation between the self-diffusion coefficient and the viscosity coefficient in a liquid is the product Θ: The Θ-values lie near 1 and for many organic compounds one mean value Θ = 1.25 can be used. A few examples of the AAD% values for D L estimated with Equation (32) using the corresponding η L data are shown in Table 7. Table 7. Comparison of self-diffusion coefficients D L,Θ = 1.25/η L ·10 8 , D L,c calculated with Equation (25) and experimental values, and D L,ex for 4 organic compounds with different polarities from N data in the temperature range ∆T. For water, 21 data with D L,Θ = 2.2/η L ·10 8 are given.  There are many experimental data for η L values available in the literature for different structural groups [12]. From such collections the corresponding terms a η and b η can be easily determined. From a η and b η values for different structures, as exemplified in the above tables, estimates of D L -values are possible using Equation (32), for which no experimental data are available.

Modeling of Diffusion in Polymers Using the Interaction Function
The diffusion behavior of organic compounds solved in polymers lies between that of liquids and solids. As a consequence, the models for diffusion coefficient estimation are based on ideas drawn from diffusion in both liquids and solids [6].
With these values the activation energies for diffusion in HDPE and PET are represented in Figure 11a as functions of M and in Figure 11b as functions of 1/T. corresponding fP-values for Equations (33) and (34) were determined. They show a linear behavior and are represented by the two following functions. For HDPE: With these values the activation energies for diffusion in HDPE and PET are represented in Figure 11a  The behavior of the activation energy of diffusion in PET in conformity with Equations (34) and (35) is in good agreement with the experimental data from [42] and [43].
As an example, for diffusion modeling in HDPE in conformity with Equation (33), eight compounds with different structures used as additives in plastics with 100 < M < 1200 are listed in Table 8. Their diffusion coefficients DP in HDPE at 60 °C are shown in Figure 12.  The behavior of the activation energy of diffusion in PET in conformity with Equations (34) and (35) is in good agreement with the experimental data from [42,43].
As an example, for diffusion modeling in HDPE in conformity with Equation (33), eight compounds with different structures used as additives in plastics with 100 < M < 1200 are listed in Table 8. Their diffusion coefficients D P in HDPE at 60 • C are shown in Figure 12. Table 8. Additives with relative molecular masses M used in HDPE and their diffusion coefficients represented in Figure 12 [44].  For an unknown thermoplastic matrix, two experimental diffusion values for a compound with known M at two temperatures are needed for the determination of the linear function fP. For a glassy matrix in the interested temperature range, two additional measured diffusion values are needed in the glassy region.
The advantage of Equation (33) is the possibility of an easy prediction of diffusion coefficients over the whole mass and temperature range for a polymer, based on just two measured points for one compound with known M at two temperatures T.
To maintain the validity of Equation (33), the total concentrations of the diffusing components should remain below ≈5% in order to avoid softening the matrix, resulting in changes (increasing) to the diffusion rate. The same precaution must be taken for the contact media used when measuring the diffusion coefficient to avoid swelling of the polymer matrix. In order to take into account such influences, an additional parameter AP can be introduced into Equation (33), which must be determined experimentally with an additional measurement [6].  Table 8 [44]; ( ) n-alkanes with 20, 27, 32 and 40 carbon atoms [6]; ( ) water [5].
For an unknown thermoplastic matrix, two experimental diffusion values for a compound with known M at two temperatures are needed for the determination of the linear function f P . For a glassy matrix in the interested temperature range, two additional measured diffusion values are needed in the glassy region.
The advantage of Equation (33) is the possibility of an easy prediction of diffusion coefficients over the whole mass and temperature range for a polymer, based on just two measured points for one compound with known M at two temperatures T.
To maintain the validity of Equation (33), the total concentrations of the diffusing components should remain below ≈5% in order to avoid softening the matrix, resulting in changes (increasing) to the diffusion rate. The same precaution must be taken for the contact media used when measuring the diffusion coefficient to avoid swelling of the polymer matrix. In order to take into account such influences, an additional parameter A P can be introduced into Equation (33), which must be determined experimentally with an additional measurement [6].

Conclusions
The interaction function in Equation (3) for liquid particles is based on the fundamental properties of natural numbers. This function acts as a bridge between properties of elementary particles and emergent properties of macroscopic systems. The atomic mass unit and an internal molar kinetic energy term for a system with optimal molar entropy are derived as two examples. This leads to a temperature T 0 = 2.98 K in a macrosystem of particles with the atomic mass unit u and the molar entropy S m = 108.85 JK −1 mol −1 . The same entropy value results for the system at a reference temperature T r = 298.15 K and a pressure of 10 5 Pa. In such an environment, the maximum of structure variety, including life, is possible.
Even one organic molecule can be treated as a system of i interacting subparticles in conformity with the interaction function defined in Equation (3). The non polar Methylene group -CH 2 -with two strong covalent bonded H-atoms is such a subparticle, with a certain individuality along the i carbon atoms in a chain. The property of chain formation is characteristic for organic compounds and the n-alkanes are therefore the reference series for all organic compounds.
A variety of physical properties of liquids is described on the basis of the interaction function. All corresponding specific values for constants are based on Equation (3).
The important conclusion is the possibility to predict values for specific properties of liquids, starting from an uniform point of view, with the interaction function w n,e using a minimum of empirical constants as simple, essentially linear temperature functions in addition to the molecular mass M.
The correlation of the power sequence in Equation (3) with the homologous series of n-alkanes as a reference for all organic compounds is one main basis. All specific structure properties can be covered with additional empirical increments with the same dimensionless units as M. This allows, in theory, for the creation of corresponding maps of such structure increments for all important functional groups, similar to charts of infrared spectra, as an example. Such data can then be easily used by introducing them directly into the prediction equations.
The correlation of emergent properties of fluid systems with the interaction function in Equation (3) provides a reliable basis for applications with a minimum of empirical data and computational effort. That way it offers an attractive alternative to many sophisticated prediction methods. Practical applications for diffusion modeling in plastic materials for food and pharmaceutical packaging and for modeling interactions with the environment are typical examples of the usefulness of such simple approaches. The interaction function represents the real basic phenomenon of the liquid state, in contrast to all models starting alone from pair interactions.