#### Abstract

In order to understand the physics behind the surface properties and nano-scale phenomena, we are motivated first to investigate the inner bond strengths as well as the effect of number of neighboring atoms and their relative distance in addition to space positions (crystallography). Therefore, in order to study the effect of the nature of metallic bond on their physico-chemical properties, we first tried to investigate and introduce a mathematical model for transforming the bulk molar cohesion energy into microscopic bond strengths between atoms. Then an algorithm for estimating the nature of bond type including the materials properties and lattice scale “cutoff” has been proposed. This leads to a new fundamental energy scale free from the crystallography and number of atoms. The results of our model in case of fundamental energy scale of metals not only perfectly describe the inter relation between binding and melting phenomena but also adequately reproduce the bond strength for different bond types with respect to other estimations reported in literatures. The generalized algorithm and calculation methodology introduced here by us are suggested to be used for developing energy scale of bulk crystal materials to explain or predict any particular materials properties related to bond strengths of metallic elements.

#### 1. Introduction

Binding strength inside the material which is mainly expressed by cohesion energy is classified as a bulk thermodynamics property. Such energy scale is believed to rule most of the physicochemical and mechanical properties of material [1, 2]. However, even the fundamental researches at nanoscale and surface science (as the new area of materials science) have been interrelated strongly with the cohesion phenomenon [3, 4]. Therefore, any generalized hypothesis at bulk, nano-scale, or even surface of materials should have been granted on the well developed cohesive scale of bulk materials.

In researches including bulk metal properties not only the macroscopic properties such as cohesion energy but also the number of neighboring atoms and crystallography of material along with the effective length of interaction energies between atoms play dominant roles in fundamental theories [5–7].

Moreover, in recent years bond energy calculation has been applied in estimation of surface properties through (atomic scale) broken bond model; thus the bond strength should be known between different neighboring atoms. In latest works of Fu et al. in [8–10] these estimations have been mostly supported by complicated quantum chemistry information through empirical electron theory where the bond energy of neighboring atoms was evaluated from the value of bonding capability of covalent electron, screen factor upon the core electron, bond length, and the number of covalent electron pair. Obviously, lack of a simple and more engineering based procedure has been realized which leads to an inquiry for developing a generalized simple method.

Therefore, we are motivated to closely study the structural effect of different crystals on their cohesion energy and evaluate a structural free energy scale which defines the materials properties. In this report, we study the effect of crystallography and number of atomic neighbors on the cohesive energy in order to define a fundamental energy scale of material. This information not only enables us to evaluate the real interatomic bond strength inside the material but also can be used as a grand energy scale to investigate other fundamental properties of elements. This knowledge gives us a tool to evaluate the strength energy from a particular atom inside the material. Also we develop free software using which anyone can compute the bond strength scale at different distances as function of number of atoms and their relative distances for any crystal materials by knowing some simple input data without any adjusting parameters.

#### 2. Theory

Most problems in solid state science require detailed study of energies along with the atomistic structural information inside the material. A detail study of crystallography effect and its related features on inner cohesive strength of pure metals needs a fundamental model which not only evolves macroscopic cohesion properties but also includes the number of atoms (those which are affected by cohesive force of an arbitrary inner atom) and their relative distances inside the materials. In order to fulfill such conditions we need to develop a model for distribution of cohesion energy between the (effected) atoms as a function of their number and relative distances from an arbitrary inner atom.

##### 2.1. Modeling the Solid Materials

In most of theoretical investigation in materials science dealing with thermodynamics description of material, the dimension of cohesion energy is described as a molar quantity. Therefore, theoretically by applying that amount of energy to the material we can dissociate one mole of material into its atoms [1]. In order to investigate the amount of any particular bond strength we need to find an algorithm for transferring macroscopic (molar) cohesion energy into atomic (inner) bond strength scale.

Regarding Figure 1(a), as long as we consider all bonds between any two particular atoms inside the material, then, the total cohesive strength of atoms leads to of total bonds. However, as we count each bond twice, the total number of bonds inside material with atoms will be where is the total number of atoms and is the number of bonds in material.

Therefore, during the transformation of molar cohesive energy of material into the atomic bond strength scale, we have to apply a 1/2 coefficient in order to avoid counting of multiple bonds between two particular atoms. Considering a simple case of Figure 1(a) where we have only 1 type of inter atomic cohesive strength of between each pair of atoms of and , then regarding (1) we can write where is the cohesive strength of atoms, is the cohesive strength between 1 inner and all neighboring atoms, and is the number of neighbors (denoted by ccn (cohesion coordination number (in case of first ccn we reach the value of classical ))).

However, in reality for an arbitrary inner atom we have more than one type of cohesive strength as Figure 1(b) shows. Therefore, if we consider a homogeneous cohesive energy distribution (in all directions of crystal cell) among each group of atoms with equal distances from an arbitrary inner atoms (far enough from the surfaces), we can write where is the total cohesion energy of th group of atoms in [kJ/mol], is the total inner cohesion energy of material in [kJ/mol], and as the atomic group that represents the bond type.

We also suppose that the distribution of cohesive energy is homogenous between all atoms in each group ccn; therefore, regarding the coefficient of atomic bond scale in (1), the total cohesive energy in (3) can be rewritten as where is the average cohesive energy (bond strength) in th atomic group between each of the atom in that group and the particular inner atom from which we count the distances in [kJ/mol] or [kJ/mol no.] (see also Figure 2).

However, as we fundamentally subtract out the effect of number of neighbors by introducing the in (4), the effect of distance on the cohesive strength between all atoms and the arbitrary inner origin atom should be considered inside the function. Such contribution attributes through soft-sphere model [7] to most of the pair potentials theories [11–14], where the parameters with energy dimensions can be extracted out from the distance or lattice parameters terms. In aforementioned models, all pair potentials (energy terms) are supposed to be constant as a material property. Therefore, developing further (4) leads us to define a distance related dimensionless function along with a fundamental energy scale , as where is the relative atomic distance of th group of atoms from an arbitrary atom in (nm) and if . then representing the maximum distance for effectiveness of bond strength.

Further on, by introducing (5) into (4), the total cohesion energy of material can be written as

##### 2.2. Nature of Distance Function

In (5) is not a simple length related (metric) function, but it represents the fundamental nature (physics) of bonds among the atoms of a particular solid at different distances with respect to an arbitrary inner atom. By applying the soft-sphere model [7] and considering all potential function theories like embedded-atom model [12] EAM, Finnis-Sinclair model [11, 13] FS, or even the simple Lenard-Jones potential [15] LG, we consider an inverse power potential form for physics of bond strength. Therefore, such hypothesis would be a quite adequate approximation for defining the distance related cohesion energy term inside the crystal metals. Thus, we can write where defines the nature of particular binding inside material in a given crystal structure and is called power of potential function.

The function should include atomic diameter as a material property and the crystallographic parameter which is the geometrical position of each atom from an arbitrary inner atom. Therefore, we can write

If we measure the distances of each atom from an arbitrary inner atom (equicentral spheres), the first coordination distance can be defined at contact condition of two neighboring atoms . Also should remain dimensionless; therefore, by summarizing the aforementioned three constraints on we can conclude that

As it can be seen from (9), value of always represents a relative distance feature of cohesion effectiveness respect to the first contact coordination neighbor. Thus, if we have only the first neighboring coordination atoms, then each atom senses a maximum cohesive strength as , while by increasing the number of atoms and being placed far from the first coordination atoms (ccn = 1), the relative potential energy exponentially decreases [11–13, 15] by power as (7) shows it.

Therefore, if we would be able to construct the mathematical algorithms of all geometrical positions in fcc, bcc, and hcp crystallographic structures, then we would be able to extract out the from the function by applying the mathematical series for each atomic group distances by . Thus, (9) can be simplified to and by substituting the result of (9) into (6) we can write where is symbolizing a function which includes and values.

Equation (11) shows that macroscopic cohesion energy can be expressed as multiplication of a fundamental energy scale and function which is a complicated function of number of atoms and their relative distances. Obviously function should converge to a value in order to express the dependency of from by the effect of atomic group distances and geometry of crystal. Thus, by knowing the total as an overall thermodynamics quantity [1], the mathematical series which describe the equilibrium atomic distances (based on crystal structures) from a particular inner atom, and the number of neighboring atoms in each distances, along with knowing the value of (power of potential function: nature of interaction between the atoms), then we are able to subtract out the influences of distance, geometry, and number of atoms. Therefore, this leads to a fundamental energy scale which is free from any crystal structure information.

Before any analytical investigation, we have to study the crystallography of different structures and evaluate the required mathematical series which represent the equal atomic distances and number of neighboring atoms .

##### 2.3. Distance and Relative Number of Neighboring Atoms in fcc, bcc and hcp Lattices

There were already some attempts in literature to evaluate the effect of number of neighboring atoms and in some cases the effect of their relative atomic distances mostly for estimation of surface properties of pure metals in fcc, bcc, and hcp structures in [8, 9, 17, 18] and [10, 19, 20], respectively. The value of bond energy in this literature is originated from either empirical electron theory, dangling bond analysis methods, atomic potential simulation, or density function theory which mostly are based on parameters such as covalent electron pairs, bond length, number of electron, contribution of couple effect between spin and orbit, bond capability, screen effect, electron density, and much other quantum chemistry information.

In their reports Zhang et al*.* in [17] considered only the effect of 1st nearest neighbors for fcc metals, and Fu et al. in [8, 9] used up to 3rd nearest neighbors for fcc and bcc, respectively, and in [10] used up to 7th nearest neighbors for hcp structures, while Wu et al. in [18, 20] considered up to 12 slabs for fcc and hcp surfaces, respectively, and Matysina in [19] considered up to 3rd nearest neighbors effect for hcp crystal structures.

Unfortunately neither the application of above literature resources nor their mathematical formulations can provide a generalized mathematical series for estimating the number of neighboring atoms and their relative distances. However, Sloane and Teo in [21] reported a magnificent but complicated mathematical formulation for theta series and magic numbers in closed packed clusters (see also Appendix A) which yet few applications are known in which; one exactly can use this information in materials science researches.

However, not only our different viewpoint of such series but also our free-software program (based on information appearing in Appendix A) enables every researcher to evaluate and apply the effect of up to 50th neighboring atomic groups and their relative distances from an arbitrary inner atom in fcc, bcc, or hcp structure. For a simple representation Figure 3 shows examples up to 3rd nearest neighbors in fcc, bcc, and hcp cells.

**(a)**

**(b)**

**(c)**

Regarding (10) and (11) mathematical series presented in Appendix A enables us to evaluate the values of and at given th neighboring groups. Table 1 shows the result of first 15 neighboring atoms and their relative distances in fcc, bcc, and hcp crystals. Appendix A shows the generalized mathematical series which is able to reproduce these numbers. In addition, developed free-software computer code is able to evaluate these calculations.

##### 2.4. Flow Chart of Free Program and Its Algorithm

To facilitate calculation of based on our mathematical model, we have written a user friendly code named IBSE-Ver1 in visual basic. IBSE stands for *inner bond strength of solid elements*. In this software, first, appropriate value of should be chosen. Second, number of bond types, , and power in potential function are determined by user. Finally, user chooses crystal structure of the element, bcc, fcc, or hcp. To easily select an element, a periodic table is shown after starting the program which enables user to choose an element from it. Initial input values including and crystal structure are predefined in the software from [1, 16], or Table 6 for each element at 0 Kelvin. However, users can introduce their own values for or crystal structure. Also users can change and to see the effect of these parameters on . The algorithm of calculations is shown in Figure 4 as described completely in the text.

After clicking the calculate button, three kinds of information are shown as output of the software. First, the value of is shown in the text box depending on the value of input parameters (, crystal structure, , ). Second, is depicted against number of selected bond types, , in a graph. Third, there are buttons which by clicking them user can see some alternative information. More information about IBSE-Ver1 is available in the help of this software.

#### 3. Results and Discussions

##### 3.1. How to Select the Input Information

Regarding the aforementioned hypothesis in (11) in addition to we need to know the maximum number of atomic distances (cutoff length) and the nature of bonds (power of potential function) inside that particular material. For a given material at a particular environment conditions while knowing the value of at given temperature and pressure for a particular lattice structure (fcc, bcc, or hcp), Appendix A enables us to evaluate the number of neighboring atoms and their relative distances from an arbitrary inner atom. Obviously, this information can be achieved from our free-software program too. However, in order to estimate the value of we need to calculate the quantity of . In other words, by knowing the value of and we are able to calculate the function which is needed to estimate the value of .

One of the advantages of modeling the nature of bond strength by introducing the dimensionless function in (9) was the extraction of materials properties such as atomic diameter (see also (10)). Therefore, function depends on the number of bond types , which should be considered along with the number of atoms in each of these groups . This pure crystallographic information is one of the outcomes in Appendix A.

For sake of simplicity we will classify all metallic elements into three main lattice groups (fcc, bcc, and hcp) and try to evaluate the quantity of for each of these groups collectively. Therefore, we can study the mathematical behavior of in each three main crystallographic structures as function of and .

Obviously, by increasing the value of (increasing the effective bond length and number of bond types), the number of atoms sensing the chosen arbitrary atom from which these atomic groups are measured increases, but the bond energy that each atom experiences decreases too. Thus, there will be a length (cutoff) above which the addition of atomic groups is unrealistic. Therefore, regarding the geometrical and mathematical nature of function we expect a convergence feature in each structure. Figure 5 shows the results of such calculation for on the example of fcc crystals up to 18 atomic slabs for different values.

As it can be seen in Figure 5, by increasing the value of the decreasing rate of cohesive force (bond strength) increases which lead to overall bigger value of function. In addition, by increasing the value of above 8-9 a convergence feature is observable. Now that the convergence of has been proved, we try to find an algorithm which enables us to estimate the optimum value of (number of maximum slabs something similar to “cutoff”) and the nature of interaction which is represented by in (10).

In classical literatures, value of attributes to ion-ion, ion-dipole, and dipole-dipole interactions behavior, respectively, whereas in latter pair potential models like LG values vary around 6 [15] while in FS model [11, 14] the square root power of distance effect for Ni, Cu, Rh, Ag, Ir, and Al was found to be 6 while for Pd and Pb it was fitted by 7 and for Pt and Au it became 8. But in the following we propose a simple algorithm using which value of and can be evaluated more adequately in a simple manner.

Figure 6(a) shows the calculated results of for different series using different values in case of bcc crystal. Obviously lower amount of leads to divergent behavior of function (see , or in Figure 6).

Proposed algorithm: regarding the behaviors of simple metals, we can overview the following simple three boundary conditions (two plus one) using which the acceptable values of and can be estimated.(a)As represents the nature of bonds between the atoms, its value should be between 6 and 9 in order to satisfy the basics quantum chemistry requirements; thus . Obviously, for a given problem can be set differently.(b) represents the longest distance where an arbitrary atom can exerts its cohesive attraction force; therefore, it shows the so-called cutoff length of crystal. In simulation algorithms usually value of 2–2.5 a (a: lattice distance) is used [11]; therefore, regarding the values of Table 1, here we choose its corresponding value of in different crystals. Therefore, can either be estimated or the indirect experimental values of effective number of atomic layers can be transferred into the length scale.(c)In addition to the above two conditions, we consider a limitation of for the relative change of function in (11) while the values of , and vary. This condition shows the acceptable mathematical convergence limit for ; therefore, by varying we have .

By applying the condition (a), while searching with condition (c) (; see Figure 6(b), e.g., of bcc lattice), the following atomic distances (slabs) have been found:(i)in fcc by : 6, 7, 8, and 9 at : 23, 15, 11, and 8, respectively,(ii)in bcc by : 6, 7, 8, and 9 at : 17, 16, 11, and 10, respectively,(iii)in hcp by : 6, 7, 8, and 9 at : 19, 19, 15, and 13, respectively.

Therefore, considering the condition (b), using mathematical series in Appendix A in fcc, bcc, and hcp structures, respectively, the values of : 11, 10, and 13 (with cutoff values of almost *≈*2.3 a, *≈*2.5 a and *≈*2.7 a) could satisfy all three conditions. These values of : , , and in fcc, bcc, and hcp belong to : 8, 9, and , respectively, which will be used here for our calculation. Obviously, upon request applying similar algorithm using our free-software program a more reliable values can be generated for any particular condition.

##### 3.2. Structural Free Cohesive Scale of Elements (Computation and Verification)

By factorizing out the crystallographic aspect of total cohesive energy regarding the classical or corrected cohesion scales [1] along with separating the effect of atomic distances in (11), the achieved energy scale represents the pure material characteristics of each particular metal. By analyzing the aforementioned hypothesis, in (11), one can say that the nature of four physical quantity , , , and ccn can be expressed as follows

Obviously, regarding the temperature dependency of relation (12) we know that by increasing the temperature from 0 Kelvin the lattice expansion along with increasing vacancy effects also should be considered. However, as a macroscopic property at different temperature represents the real behavior of materials including the vacancies and lattice parameters; therefore, for the time being we neglect such contribution. The results of our calculation for at 0 K and at for fcc, bcc, and hcp metal materials are presented in Table 2. As it can be seen some metals melt in different lattice crystal than they are usually in solid state at 0 K. Applying the attached free-software, one can produce more reliable values for structural free energy scale separately for any metallic element.* (The attached free-software in this paper is available at our website *(http://www.fcrgroup.org).* At our website under the Achievements main menu goes to Setups and Software submenu and find IBSE-Ver1 software. The terms of service/use (TOS) for attached free-software ARE the correct citation to this paper.). *

is mainly affected by the nature of atomic interaction between each pair of the inner atoms. Moreover, considering the potential function phenomenon developed in literature and our aforementioned hypothesis, we can say that . Thus, by rearrangement of (11) we can write

###### 3.2.1. Verification with Melting Points

The bulk melting phenomenon in crystalline materials takes place when the bindings of crystallographic cells lose its structure; thus the best verification properties for cohesive energy scale are believed to be the melting point [1]. The value of cohesion energy has been defined in classical literature as a negative value of sublimation enthalpy while in the past decades there have been several attempts for correcting these values in order to fit the melting point more adequately [1, 2].

However, regarding the logic mentioned in (13) we expect that should verify more adequately the main physicochemical properties of materials such as melting with respect to classical cohesive energy scale. In the following this verification for three different main crystallographic forms has been shown (in Figure 7).

As our calculation results show, the correlation in three different crystallographic forms is perfectly represented by considering the boundary conditions of and . Therefore, we may introduce the concept as a simple structural free cohesive energy scale for pure crystal elements as a fundamental correlation scale regarding the physic-chemical properties of metals.

###### 3.2.2. Verification with Bond Strength for Different Neighboring Atoms

Applying the aforementioned computational algorithm in (11) and considering the relative atomic group distances in Table 1 along with the results of Appendix A, we are able to evaluate the bond strength inside the fcc, bcc, and hcp metals by in (4). simply expresses the bond strength at th atomic groups (bond type).

The latest reported about the applications of bond strengths were for estimating the surface properties of pure metals in works of Fu et al. in [8–10] which the bond strengths were calculated using the empirical electron theory. Fu et al. in [8, 9] used up to 3rd nearest neighbors (bond types: A, B, and C) for fcc and bcc, respectively, which are compared to our calculation results in Tables 3 and 4. Obviously, not only the simplicity of our model but also the fact that neither the number of atomic neighbors nor the nature of bonds is a limitation here could give extra advantages for our method. In addition to proposed value of and in Table 2, a set of calculation with (limitation in [8, 9] for fcc and bcc, resp.) and for power of potential function shows more similarity to values in literatures (see Tables 3 and 4).

Fu et al. in [10] reported up to 7th nearest neighbor (bond types: A, B, C, D, E, F, and G) bonds strength for hcp structures which is compared also to our results in Table 5. As our hcp crystal has considered to be a perfect cell, between type A and B and type E and F the nearest neighbor value in [10] has been considered for the first and forth type bond strength. As a result, we evaluate up to 5th nearest neighbor values from [10] to compare with our results.

In case of hcp metallic crystals we also evaluated a set of calculation with (limitation in [10] for hcp) while keeping for power of potential function with respect to Table 2 which shows more similarity to values in [10]. These calculations illustrate the ability and flexibility of model and software for fitting and adjustment. As it can be seen in Tables 3–5, not only the tendency of rate of decreasing the bond strengths is similar, but also the absolute estimated values adequately fall in a comparable energy zone. Moreover, using our free-software and generalized algorithms for any metals in any crystal structure at a particular temperature, one can regenerate a new set of data.

#### 4. Conclusions

A detailed investigation on the crystallography of fcc, bcc, and hcp lattices has been performed. This application of mathematical theta series gives us the information about the number of atoms and their relative equal distances from a central inner atom for up to 50th layers. Therefore, as long as the cutoff length (longest atomic distances from which an atom could sense the attraction of other one) is known, the maximum number of neighboring atoms and their distances from the central atom can be calculated.

Using these data, we introduced a model for transformations of bulk molar cohesion energy into inter atomic bond strengths inside the metallic elements. Obviously, in addition to number of neighboring atoms and their distances, the nature of bond type of particular material (as power in an inverse power potential) in given crystal structure should be known.

Therefore, an algorithm for estimating the value of neighboring atoms and distances from the central atom as function of crystal type and type of bond (nature of bond strength) has been proposed and for the sake of simplicity three collective results for all fcc, bcc, and hcp crystals have been evaluated and reported here.

The effect of number of neighboring atoms and their distance along with the physics of bond type inside the metallic elements has been investigated. This study leads to defining an energy scale free from crystallographic information. The aforementioned energy scale presented here is suggested to be used for scaling any fundamental properties of metallic elements which are interrelated to their inner cohesive feature.

We developed a free-software (IBSE-Ver1) which enables us to evaluate the number of neighboring atoms and their relative distance from a central atom in fcc, bcc, or hcp up to 50th layers. The free-software is able to calculate the structural free energy scale and interatomic bond strengths at different neighboring atomic distances depending on the total number of atoms, cutoff length, and the nature of bond type. The results of our calculation in case of structural free energy scale for the first time perfectly reproduced the expected classical hypothesis of linear dependency of bond energy and melting points. Also the results of bond strengths for different type of metallic elements adequately produced the tendencies reported in literature.

Obviously using the enclosed free-software more reliable values of structural free energy scale can be evaluated separately for each metallic element via considering better set of cutoff length and power of potential energy related to nature of bond types in particular crystal structure.

In next paper the mathematical series and other free-software for estimating the surface bond strength and related surface properties of all metallic elements at different crystal structures and crystal planes will be submitted.

#### Appendices

#### A. Application of Theta Series for Calculating Number of Neighboring Atoms and Their Relative Distance from a Central Atom in a Close Pack Structure

Assuming the origin of coordinates at the center of an arbitrary atom in a close pack structure, theta series are able to calculate number of neighbors and their distance from this central atom. These series are used for 2D (circular) or 3D (spherical) close pack structures. Let denote the number of atoms at distance from the origin. These atoms form a shell of radius . Central atom, its neighbors, and shells are shown in Figures 8(a) and 8(b) for 2D square lattice and 2D hexagonal lattice, respectively.

**(a)**

**(b)**

Theta series for a close packing structure () is given by

This equation is a power series in the variable . is the number of atoms (ccn in (4)), is the distance from the origin ( in (5)), and is norm of a vector .

For example in 2D square lattice shown in Figure 8(a), theta series is described by the following equation:

This power series predicts that there are 4, 4, 4, 8, 4,… neighbors at distance of , , , , and , respectively, from the central atom as shown in Figure 8(a).

Similarly, for hexagonal lattice shown in Figure 8(b), theta series is introduced as follows:

In the same way, all 2D and 3D close pack structures can be expressed using the following three special Jacobi theta series:

The Jacobi theta series satisfy large number of identities using which they can be simplified.

For some clusters the following series is also helpful: where is a nonzero number.

For fcc structure, theta series are given below which can predict up to 10 numbers of neighbors and their distance from central atom:

Similarly for hcp and in bcc

Sloane and Teo in [21] considered the (first coordination distance) by , , and during the construction of theta series for fcc, bcc, and hcp structure, respectively. Although in this paper we considered a more general algorithm of . Therefore, for transferring the series of distances for fcc, bcc, and hcp in (A.6), (A.8), and (A.7) into our values in Table 1, in above equations one should multiply the theta series values by , , and , respectively. For more information about theta series and magic numbers [21] is useful.

#### B. Selected Values of

For more details see Table 6.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.