The Lattice Compatibility Theory : Supports from the Generalized Simha-Somcynsky Chemical Physics-Related Theory

The earliest models used in the study of lattice structures are mean field theories, which do not contain structural dependence. The Lattice Compatibility Theory (LCT) proposes here a novel framework where the measure of the disorder is based on Urbach tailing features and lattice matching features between the host matrix and doping agent intrinsic structures. This study has been implemented on a particular compound (BTO:Co) and refers to the Simha-Somcynsky (SS) theory, a mean field theory where the measure of the disorder is stated as holes.


Introduction
The knowledge of doping agents behaviors within host lattice matrix is of considerable importance for the optimal design for applications such as semiconductor windows functional glasses, transparent electrodes in flat panel displays, buffer layers, and solar cells [1][2][3][4][5][6][7][8][9].Although such behaviors have been studied for many host-doping agent lattice systems, theoretical fundaments and updated principles are still important for predicting or correlating the stability behavior of many systems within a wide range of lattice shapes.
The first theories based on mean field theory and independent from the design of lattice structures failed in the statistical thermodynamics of branched macromolecules.Studies on branched structures served as attempts to mathematically correct the mean field theories.One of the most important of these studies is the Lattice-Cluster Theory (LCT), developed by Freed and Bawendi [1].
For complex lattice systems, other theories have been developed.Dee and Walsh [2,3] proposed the lattice theory as a tool for depicting the thermodynamic properties of heterogeneous structures.This theory was an enhancement of those of Flory [4,5] and Huggins [6] concerning chain structures.Lennard-Jones and Devonshire [7,8], and Prigogine et al. [9][10][11] developed the cell model, according to which molecules can move in neighbors-induced potential holes.Later, Simha et al. [12,13] noted that free volume theories, especially hole theories, could delineate the thermodynamic properties of solid structures, and they introduced the notion of free volume.The so-called Simha-Somcynsky (SS) model [12] defined accurately the holes in the lattice-like structure and determined the statistical behavior by finding out the combination of the holes formed within lattice intermolecular sites.
In the hole theory, a major change in volume is explained by the number of holes, and the change in cell size plays a minor role while in the cell theory [14][15][16], the changes in volume with changes in temperature and pressure can be explained only by the changes in cell size.Oppositely, the lattice theory stipulates that lattice size is fixed, and volume change is explained solely by the number of vacant sites [15,16].
Zhong et al. [17] indicated that the Simha-Somcynsky (SS) model could not generate sufficient free volume and still lacks fluidity.They modified the SS model by using the open-cell concept, and they induced a new empirical spatial parameter.

International Journal of Chemical Physics
In this paper, it is outlined that SS theory was independent from the macromolecular architecture, and it exhibited differences in behavior concerning some doped lattices.Hence new elements for explanation of these differences are presented and discussed in the framework of the Lattice Compatibility Theory (LCT).The paper is organized in the following way.In Section 2, we present some relevant theoretical fundaments to both Simha-Somcynsky Lattice-Hole and Lattice Compatibility Theories along with some patterns of the studied BTO lattice.In Section 3, we present comparative and conjoint analyses.Section 4 is the conclusion.

Theoretical Fundaments
2.1.Elements from the Generalized Simha-Somcynsky Lattice-Hole Theory.According to the generalized Simha-Somcynsky theory [12,13], any lattice can be considered as a succession of elementary molecules and holes (Figure 1).Each cell in the occupied fraction is either empty or contains the molecule van-der-Waals volume as well as an inherent free volume.
In this theory's framework, the most important aggregate which traduces entities repartition within the lattice is the volume fraction   .This parameter is the ratio of the occupied sites by a molecule  on all of lattice sites: where  is the number of molecules,  ℎ is the number of unoccupied cells, and   is the number of equivalent segments in a molecule .
The total occupied volume fraction is coupled with temperature through the minimization condition where   (,   ) is the molar Helmholtz free energy for a molecule  within the given structure at a temperature , and approximated by where  is universal gas constant (J/mol K),  is absolute temperature, and   is the number of nearest neighbors of the molecule .
The configurational partition function   (,   , ) only depends on the positional part of the degrees of freedom and the potential energy of the molecule.All molecules with the same mass will have identical contributions from their kinetic degrees of freedom (the momenta), but they may have distinct modes of interacting which will result in different contributions from the configurational part.Ordinarily, the kinetic and potential energy contributions to the canonical partition function are separable.The kinetic contribution can be evaluated analytically, giving, for small molecules, the thermal de Broglie wavelength.For complex molecules, and by developing good expansion of multidimensional integral over the 3D space and statistical mechanics-related perturbation techniques, the configurational partition function is represented as where (,   ) is the combinatorial factor, V  is the free volume,  0 is the internal energy when all of the segments are located at the center of the cell,  is Boltzmann constant, and  are the external degrees of freedom per molecule, as defined by Prigogine et al. [9][10][11] and Hirschfelder et al. [18].The combinatorial factor (,   ), arising from the mixing of molecules and holes, represents the occupation possibilities of the molecules into the lattice sites.It is obtained through the number of molecules  and volume fraction along with the number of nearest neighbors of the molecule as follows: The free volume V  is, as defined by Park and Kim [19], the difference between the volume per molecule and the hardcore volume of the molecule: where  is the cell volume,  * is the molecular characteristic volume, and  is the lattice geometric factor (i.e.,  −2 = 2, for face-centered cubic lattice).

Elements from the Lattice Compatibility Theory Analysis.
The Lattice Compatibility Theory, as mentioned in some recent studies [20][21][22][23], is based on the interaction of dopingelement lattice behavior versus host edifice.Preludes to this theory have been established by Boubaker [21] in the context of analysing Urbach tailing controversial behaviour in some nanocompounds as well as I-III-O 2 ternary oxide ceramics instability at low temperatures beside ternary tellurides and sulphides.They were also confirmed by Petkova and Boubaker [20] on the basis of investigation on some copperdoped bismuth sillenites (BSO) and germanates (BGO) compounds.

The stability of doping agents inside BTO host structures is favorized [sic] by geometrical compatibility, expressed in terms of matching patterns between doping agent intrinsic lattice and those of the host.
In the actual discussed case (cobalt-doped BTO lattice), the nature of the highest occupied band and the location of holes in BTO lattice structures have been demonstrated to be determinant.In this context, fundamental geometrical observations concerning the structure of BTO and the doping lattice were interpreted in terms of conventional orbitals patterns-linked geometry.
The Lattice Compatibility Theory tried to give a plausible understanding of the disparity concerning doping element incorporating dynamics starting from intrinsic doping element lattice properties in comparison to those of the host.In precedent study [20][21][22][23] materials, changes in Urbach energy have been associated with a cobalt-doping induced disorder in BTO matrices.In order to understand Urbach tailing alteration following doping agent insertion in host structures, Urbach energy   has been determined for doped and undoped samples through the equations where (ℎ]) represents, for each sample, the experimentally deduced optical absorption profile.Urbach energy   is a measure of the inhomogenous disorder and atomic scale dispersion inside structures as it indicates the width of the band tails of the localized states in presence of defects.Its analytical formulation deduced takes into account three components: structural disorder, carrierphonon interaction, and carrier-impurity, as follows:  where Ω 1 , Ω 2 , and Ω 3 are constants,   is Boltzmann constant,  is lattice strain related with the structural disorder,   is Debye temperature,   is Debye length,  * is carrier effective mass,  is impurity charge,  is electron charge,  is static dielectric permittivity, and ℎ is Planck's constant.
The width of the localized states (band tail energy or Urbach energy   ) has been estimated from the slopes of the plots of ln (]) versus energy ℎ] (Figure 2).In this context, fundamental geometrical observations concerning the structure of BTO and the doping lattice are interpreted in terms of comparison to cobalt intrinsic lattice parameter (Figure 2) as well as the fact that in oxygen-rich edifices, interstitial Co 2+ tetrahedral sites split the five degenerated atomic d orbitals into two groups, leading to 3 unpaired d electrons on Co 2+ , while all the d electrons of polyhedral Co 3+ are paired (with a null global magnetic moment), as illustrated by Figure 3.
Using a complete set of measurements, main lattice constants of cobalt intrinsic lattice have been compared to those of BTO and revealed an obvious compatibility with the cubic main metric and angular parameters of Co 2+ intrinsic lattice.It has been recorded that Co 2+ incorporation in BTO matrix was followed by a loss of the hopping motion of electrons which decreased the piling up of electrons at host matrix vacancies hence impeding the buildup of space charge polarization.The variation in lattice parameters with Co 3+ incorporation had been attributed to the increase of the unit cell volume of the host lattice with increasing Co 3+ content.In fact, since Co 3+ ions are trivalent and in order to maintain the charge balance during doping, incorporation had to be accompanied by an increase of unit cell volume either by a reduction of host cations valence or an oxygen content increase.
In the same context, and according to Muncaster et al. [24], Huang et al. [25], and Chmaissem et al. [26], cobalt ions intercalation dynamics in similar host structures with fixed coordination geometries were mainly governed by differences in ionic size and magnetic properties.Relative induced disorder in BTO:Co 3+ edifices show that it is fundamentally the matter of oxygen uptake.The amount of incorporated oxygen seems to be wholly dependent on the presence of additional locations where the bond lengths from oxygen to surrounding cobalt ions remain within an upper limit [26].

Lattice Compatibility Theory (LCT) and Simha-Somcynsky Conjoint Analysis
The most fundamental structure alteration which occurred in the host BTO lattice has been interpreted in terms of changes either in the configurational partition function   (,   , ) or in the lattice parameters, according to Simha-Somcynsky analysis and Lattice Compatibility Theory, respectively.In the Simha-Somcynsky analysis, the structure-dependent coefficients (  in (3)) increase in different ways during Co 2+ and Co 3+ incorporation (Figures 4 and 5).
In this context, for a molecule , the relative contribution of the doping agent is expressed in terms of relative variation of the molar Helmholtz free energy as follows: Under a first-order approximation, cobalt ions edifices presence induces a drastic decrease of Helmholtz free energy as well as a lower Urbach energy.It was recorded that these trends are not the same for Co 2+ and Co 3+ ions.The loss of stability which characterized the transition from Co(II) to Co(III) state and which was recorded by Amine et al. [27,28], Duncan et al. [29], Liu et al. [30], Santhanam and Rambabu [31], and Aklalouch et al. [32,33] has been satisfactorily explained in the Lattice Compatibility Theory (LCT) framework.The most probable explanation which can be provided through the Simha-Somcynsky theory lies in the expression of the Helmholtz free energy (see ( 3) and ( 7)).Indeed, the incorporation of cobalt ions contributes to the decrease of the Helmholtz free energy [34][35][36][37], but, as it is also equivalent to the increase of the number of occupied sites within the cell (Figures 4 and 5), it induces an increase of energy through the first two terms in (3).Moreover, the latter contribution of Co 3+ ions is more important due to its voluminous polyhedral configuration (Figure 5).

Conclusion
The presented work showcases some fundaments of the Lattice Compatibility Theory (LCT) framework in relation with the precedent Simha-Somcynsky theory-linked analyses.The most probable explanations which can be provided to some intriguing nanoscale unexpected alterations, through the LCT theory, have been discussed in terms of the well-known Simha-Somcynsky (SS) lattice hole theory.What this theory and more sophisticated frameworks have in common is that they incorporate a description for holes on the lattice, acting as free volumes.This approach, which became necessary since the lattice-based arrangement of molecules became the most encountered in the recent literature, was confronted to the results which have been recorded in an unupdated compound (BTO:Co).The results were as accurate as relevant.The LCT analyses used in lattice linear structures can also be used together with the SS model for other similar structures.
Figure 2: Plots of doped BTO optical coefficient versus energy (as guides to Urbach tailing quantification).