^{1}

^{1}

^{2}

^{*}

In this paper, a mathematical relation was found between interatomic Hooke’s force constant and both the bulk modulus and interatomic distance in solid crystals, considering that the forces which have effect on an atom are only those resulted from the neighboring atoms, and the forces are subject to Hooke’s law as the deflections of atoms from their equilibrium positions are very small. This work has been applied on some solid semiconducting crystals of diatomic primitive cell, including crystals of mono-atomic primitive cell automatically, by using linear statistical fitting with computer programming and, then, using mathematical analysis, proceeding from the vibrational dispersion relation of solid linear lattice, these two methods have been used in the process in order to support each other and for the result to be satisfying and reasonable. This is a contribution to the process of using computer programming in physics to facilitate mathematical analyses and obtain the required relations and functions by designing and developing appropriate computer programs in line with the macro and micro natures of materials. The importance of this is in enhancing our understanding of the interatomic actions in cells and of the crystal structure of materials in general and semiconductors in particular, as it is a step of the initial steps to facilitate the process of calculating energies and extracting mathematical relations between correlation energy and temperature as well as between sub-fusion and fusion energies with temperature.

Obtaining a mathematical relation of interatomic Hooke’s force constant to bulk modulus and interatomic distance in solids, especially in semiconductors, is very important as a relation connecting macroscopic and microscopic quantities that helps in understanding the interaction between atoms and calculating energies and other quantities. Dedicated bond force constant and bulk modulus of Cn fullerenes (n = 20, 28, 36, 50, 60) are computed using density functional theory (DFT) and finite element analysis (FEA) by Peon et al. The bond force constants predicted by DFT are then used as an input for finite element analysis (FEA) of the fullerene [^{nd} degree F = − β x + γ x 2 [

In this paper, it has been assumed that Hooke’s force is the only affecting interatomic force under nearest-neighbor approximation ( F = − β U ), which was used to find a mathematical expression for solid crystals of diatomic primitive cell, including those of mono-atomic one automatically, relating interatomic Hooke’s force constant (β) to bulk modulus and interatomic distance by using two methods, linear statistical fitting with computer programming and mathematical analysis. We used the vibrational dispersion relation of a linear lattice of diatomic primitive cell and Data of some semiconductors (C (diamond), Si, Ge, SiC, ZnS, ZnTe, CdS and CdTe).

In spite of existing other forces than Hooke’s one and which may be of other atoms than the neighbors, but that do not affect much the process of calculating energies and other quantities as was shown in Einstein and Debye theories for specific heat in solids [

In order to get the mathematical relation that we seek, we considered that the interatomic forces in a linear lattice are subject to Hooke’s law with small atomic deflections as ( F = − β U ), where F is the Hooke’s force, U is small deflection and (β) is the force constant, with neglecting other effects to facilitate the process mathematically for scientific purposes [

As we know, the two functions of the dispersion relation which describe linear vibrations of a diatomic primitive cell as shown in

The acoustic branch

ω a c 2 = β [ ( 1 M + 1 m ) − ( 1 M + 1 m ) 2 − 4 m M sin 2 ( q a 2 ) ] , (1)

The optical branch

ω o p 2 = β [ ( 1 M + 1 m ) + ( 1 M + 1 m ) 2 − 4 m M sin 2 ( q a 2 ) ] , (2)

where m and M are the masses of both atoms in the primitive cell ( M ≥ m ) and a is the interatomic distance. The possible deflections of two atoms are

U ( l ) = U m e i ( q ⋅ l a − ω t ) (for the bigger mass), (3)

u ( l ) = u m e i ( q ⋅ l a − ω t ) (for the smaller mass), (4)

where U_{m} and u_{m} are the maximum deflections for both atoms, ω is the angular frequency and t is the time [

Where the allowed values of q (wave number) are given as follows:

q = l 2 π N a , N is the number of primitive cells in the crystal, l = ± 1 , ± 2 , ± 3 , ⋯

This means that all the solution points are within the first two Brillouin zones, so, if the number of atoms in each primitive cell is three, then all the solutions are within the first three Brillouin zones [

We proceeded from the acoustic branch of the dispersion relation in a linear solid lattice of diatomic cell (Equation (1)) and the fact that where q approaches zero, it approximately takes the following form [

ω = s q , (5)

where s is the wave (sound) speed via the lattice, as shown in

We should make some mathematical simplifications in Equation (1) (the acoustic branch) to remove the function dependence on the physical characteristics as follows.

Equation (1) may be rewritten under the following form:

ω 2 = β m [ ( m M + 1 ) − m 2 M 2 + 1 + 2 m M cos a q ] , (6)

or

m β ω = ( m M + 1 ) − m 2 M 2 + 1 + 2 m M cos a q , (7)

assuming that x = a q , p = M m ( M ≥ m , then p ≥ 1 )and ϖ = m β ω , then

ϖ = ( 1 + 1 p ) − 1 + 1 p 2 + 2 p cos x . (8)

By using a VISUAL BASIC program that generates 50 ( ϖ , x ) points beside x = 0 via substituting with values of x in Equation (8) and uses them for linear statistical fitting as ϖ = η 0 + η 1 x after calculating each of η 0 and η 1 , we can reach the sought relation.

When running the program which needs to enter the values of p = M/m for each matter, the η 1 values are as

η 1 | p = M m | M ( × 10 26 ) [ | m ( × 10 26 ) [ | Crystal |
---|---|---|---|---|

0.495478885554644 | 1 | 1.994022912 | 1.994022912 | C-C |

0.495478885554644 | 1 | 4.663788810 | 4.663788810 | Si-Si |

0.495478885554644 | 1 | 12.05213349 | 12.05213349 | Ge-Ge |

0.381690077840471 | 2.338884263 | 4.663788810 | 1.994022912 | Si-C |

0.400556205681922 | 2.039301310 | 10.85505562 | 5.322928773 | Zn-S |

0.406612439957611 | 1.951667177 | 21.18545600 | 10.85505562 | Zn-Te |

0.327164494602450 | 3.505926388 | 18.66179645 | 5.322928773 | Cd-S |

0.479475162402413 | 1.135231317 | 28.76019326 | 18.66179645 | Cd-Te |

Because that the curve ϖ ( x ) passes through the origin (0, 0), so η 0 ≅ 0 and η 1 is the tangent slope of linear fitting function ϖ = η 1 x besides x = 0 .

ϖ = η 1 x ⇒ m β _ ω _ = η 1 a q , (9)

where ω _ and β _ are calculated by linear fitting process, then

ω _ = a η 1 β _ m q , (10)

which we compare with Equation (5), so

ω _ = a η 1 β _ m q ≅ s q , then (11)

β = m s 2 a 2 η 1 2 . (12)

Proceeding from Equation (8), which is

ϖ = ( 1 + 1 p ) − 1 + 1 p 2 + 2 p cos x .

A first estimation of the force constant β may be derived from simple consideration. Indeed, as x approaches zero, cos θ ≅ 1 − x 2 2 , so, Equation (8) can be rewritten under the following form

ϖ 2 = ( 1 + 1 p ) − 1 + 1 p 2 + 2 p ( 1 − x 2 2 ) , around x = 0 (13)

or ϖ 2 = ( 1 + 1 p ) [ 1 − 1 − x 2 p ( 1 + 1 / p ) 2 ] , (14)

for values of x around zero, we have x 2 p ( 1 + 1 / p ) 2 ≪ 1 , and, consequently

(as ( 1 + y ) 1 / 2 ≅ 1 + 1 2 y as y ≅ 0 ), Equation (14) reduces to the simple form

ϖ 2 = x 2 2 ( p + 1 ) , (15)

as x = a q , p = M m ( M ≥ m , then p ≥ 1 )and ϖ = m β ω , then from Equation 15, one can obtain the following expression of ω .

ω = ( β 2 ( M + m ) a ) q (16)

According to Equation (5), we get

s = ( β 2 ( M + m ) a ) , (17)

then β = 4 m ¯ s 2 a 2 , (18)

where m ¯ = m + M 2 .

Comparing the estimated force constant and the analytically extracted one, we get the results shown in

We can now introduce the interatomic force constant, β, via the following relation

β = 4 m ¯ s 2 a 2 ,

β = 4 m ¯ s 2 a 2 Extracted | β _ = m s 2 a 2 η 1 Estimated | S (m/s) | B [^{11}) | a [^{−10} m) | Crystal |
---|---|---|---|---|---|

78.31346728 | 79.7487862 | 11155.092 | 4.4000 | 3.56 | C-C |

26.90453009 | 27.39779461 | 6520.984 | 0.9908 | 5.43 | Si-Si |

24.05782553 | 24.49891544 | 3998.367 | 0.8501 | 5.66 | Ge-Ge |

45.98786306 | 47.27084074 | 8084.071 | 2.1142 | 4.35 | Si-C |

13.61016036 | 13.95506009 | 3508.744 | 0.5032 | 5.41 | Zn-S |

12.30523188 | 12.60025426 | 2673.068 | 0.4046 | 6.10 | Zn-Te |

13.34362773 | 13.84632183 | 2942.99 | 0.4783 | 5.82 | Cd-S |

9.939901797 | 10.12465468 | 2288.505 | 0.3068 | 6.48 | Cd-Te |

s = B ρ = B v c m + M ^{,} ρ is the density [

v c is the primitive cell volume, such as v c = a 3 n _ [

where n _ is the number of atoms in a primary cell;

Substituting from Equation (20) in Equation (19) then in Equation (18), we get

β = 2 n _ B a , (21)

which is the relation we are looking for.

As an example, for F. C. C. Structure, as n _ = 4 , β = 1 2 B a .

It is worth noting that a similar relation has been shown to be valid for solids with a variety of crystal structures [

Comparing Equation (12) and Equation (18), we find

m s 2 a 2 η 1 2 = 4 m ¯ s 2 a 2 , (22)

then m s 2 a 2 η 1 2 = 4 m ¯ s 2 a 2 , (23)

which gives us the exact values of the slope η 1 for the matters used in this work as shown in

As direct results from the obtained force constant relation, we can calculate the maximum acoustic angular frequency ( ω a c ) max , the minimum optical angular frequency ( ω o p ) min and the maximum optical angular frequency ( ω o p ) max , which can easily be measured experimentally.

From

β = 1 2 M ( ω a c ) max 2 = 1 2 m ( ω o p ) min 2 = 1 2 μ ( ω o p ) max 2 , (24)

where μ is the reduced mass ( 1 μ = 1 M + 1 m ) .

By substituting with β = 2 ( M + m ) s 2 a 2 , we can deduce that

Crystal | p = M / m | η 1 Extracted | η 1 Estimated |
---|---|---|---|

C-C | 1 | 0.5 | 0.495478885554644 |

Si-Si | 1 | 0.5 | 0.495478885554644 |

Ge-Ge | 1 | 0.5 | 0.495478885554644 |

Si-C | 2.338884263 | 0.386976257 | 0.381690077840471 |

Zn-S | 2.039301310 | 0.405600166 | 0.400556205681922 |

Zn-Te | 1.951667177 | 0.411577201 | 0.406612439957611 |

Cd-S | 3.505926388 | 0.333114054 | 0.327164494602450 |

Cd-Te | 1.135231317 | 0.483907705 | 0.479475162402413 |

( ω o p ) max (×10^{13}) | ( ω o p ) min (×10^{13}) | ( ω a c ) max (×10^{13}) | S (m/s) | p = M m | a (×10^{−10}) [ | Crystal |
---|---|---|---|---|---|---|

12.534 | 8.863 | 8.863 | 11155.092 | 1 | 3.56 | C-C |

4.804 | 3.397 | 3.397 | 6520.984 | 1 | 5.43 | Si-Si |

2.826 | 1.998 | 1.998 | 3998.367 | 1 | 5.66 | Ge-Ge |

8.115 | 6.792 | 4.441 | 8084.071 | 2.338884263 | 4.35 | Si-C |

2.761 | 2.262 | 1.584 | 3508.744 | 2.039301310 | 5.41 | Zn-S |

1.852 | 1.506 | 1.078 | 2673.068 | 1.951667177 | 6.10 | Zn-Te |

2.435 | 2.147 | 1.147 | 2942.99 | 3.505926388 | 5.82 | Cd-S |

1.214 | 0.885 | 0.831 | 2288.505 | 1.135231317 | 6.48 | Cd-Te |

( ω a c ) max = 2 ( 1 + 1 p ) ⋅ s a , (25)

( ω o p ) min = 2 ( 1 + p ) ⋅ s a , (26)

and ( ω o p ) max = 2 ( 1 + p ) p ⋅ s a , (27)

which

The simple approach, detailed in this work, based on linear fitting by computer programming and mathematical analysis, permitted computation of the constant force (CF) as a linear function of the bulk modulus and interatomic distance parameters, under the assumption that only interactions between first nearest neighbors are considered. The two derived relations from the two above-mentioned methods have led to similar results for some mono- and diatomic molecules of the semiconductor type. Furthermore, these two relations enabled us to express the interatomic force constant as proportional to the bulk modulus and the interatomic distance.

The importance of deriving this relation is that many physical quantities are linked to the force constant, the bulk modulus and the interatomic distance and so, it could be a contribution to increase crystal structure understanding and may be used in many other theoretical works.

Moreover, this resulted mathematical relation makes it possible to calculate interatomic Hooke’s force constant for solid matters in terms of bulk modulus and interatomic distance, and so to calculate bulk modulus in terms of the others, which can help understanding the relationship between sub-fusion and fusion energies with fusion temperature.

The results of this work will be used in an attempt to extract mathematical relations between correlation energy, sub-fusion energy and fusion energy with temperature, which will clarify many physical issues and provide broader research possibilities.

The authors declare no conflicts of interest regarding the publication of this paper.

Joghlaf, M., Ababou, Y. and Sayouri, S. (2021) A Simple Approach to Compute Interatomic Force Constant for Mono and Diatomic Semiconductors. Journal of Applied Mathematics and Physics, 9, 11-20. https://doi.org/10.4236/jamp.2021.91002