Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access August 26, 2020

Multivariate analysis for the classification of copper–lead and copper–zinc glasses

  • Dimitar Dimitrov , Miroslava Nedyalkova EMAIL logo , Sergio Madurga , Ludmila Naneva and Vasil Simeonov
From the journal Open Chemistry

Abstract

The similarity patterns in the physicochemical properties of copper–lead and copper–zinc borate glasses were identified by means of finding similarity within the objects of study using multivariate statistical analysis. As exploratory methods of multivariate analysis, cluster analysis, principal components analysis, and two-way clustering were applied for a set of copper–lead and copper–zinc borate glasses. Specific correlations among the physicochemical properties of copper glasses were interpreted. In particular, the effect of Pb and Zn doping metal ion in copper glasses in the structural and mechanical properties is identified. Interestingly, the degree of lead content determines two kinds of glasses with specific physicochemical properties.

1 Introduction

In the design of glass compositions, transition metals oxides being a constitutional part of the glass systems could play a dual role related to their physicochemical and structural properties (coordination number, different redox state at specific conditions, etc.) [1,2,3,4,5,6,7,8].

Oxide glasses exhibit interesting properties such as density, elastic moduli, hardness, glass transition temperature, and liquid fragility, and the changes in chemical composition cause significant alternations in all measured physical properties. It was also found that the change of copper content in copper–lead borate glasses results in a slow disproportionation of B-containing groups due to change in the coordination number of Cu+ when compared with Cu2+ ions. It causes changes in the physicochemical properties such as hardness, thermal stability, density, thermal expansion coefficient, and chemical durability. Zinc metaborate glasses are subject to an enhanced disproportionation of metaborate in other structural units since the lead metaborate series exhibit different structural alternations [9,10,11,12,13,14,15,16].

Specific changes in the physicochemical parameters of the borate glasses due to changes in their chemical composition could be subject to another experimental data approach, namely, intelligent data analysis (exploratory data mining) [17,18]. This option is rarely used in the analysis of the properties of noncrystalline materials. Multivariate statistics methods such as cluster analysis (CA) or principal components analysis (PCA) offer a new opportunity for data modeling, classification, and interpretation. They are superior with respect to their information content to the traditional correlation analysis if relationships between physicochemical properties and chemical composition are sought.

Research objectives in the current work may be linked with a validation of the efficiency of exploratory data analyses to study and interpret a set of experimental data from the literature source including values for different physicochemical parameters in different copper–zinc and copper–lead metaborate glass series. In addition, relationships for chemical composition versus physicochemical parameters are evaluated by regression analysis.

2 Methods

2.1 Multivariate statistical methods

In the present study, two multivariate statistical methods were used: CA and PCA.

CA is a well-known approach for searching similarity patterns (called clusters) within the data approach to reveal groups of similarity (clusters) within large data sets. The identification of the similarity groups makes it possible to better interpret specific relations between the objects of the study (glass types) or among the parameters characterizing them (physicochemical variables).

Three different modes of clustering were applied: hierarchical cluster, nonhierarchical cluster, and two-way CA.

Hierarchical CA is a nonsupervised method for linkage of objects in clusters that follows a simple algorithm:

  • normalization of the raw experimental data (to avoid differences in dimensionality which could deteriorate the clustering; usually the z-transform procedure is applied);

  • introduction of a similarity measure (Euclidean distances or squared Euclidean distances are often calculated as a measure of object similarity);

  • choice of linkage procedure (among many options we have chosen the Ward’s method of linkage);

  • selection of criteria for determination of the cluster significance, as the Sneath’s criterion based on the Dmax distance in the similarity matrix.

The output plot of the analysis is called a dendrogram, which allows the visual determination of significant clusters for later interpretation. No preliminary conditions are introduced for hierarchical clustering [19].

Nonhierarchical clustering belongs to supervised pattern recognition methods. It requires a preliminary determination of the number of clusters to which the objects (or variables) should belong. Therefore, it demands a priori hypothesis for the number of clusters for the proper preselection of the number of clusters based on expert opinion or preliminary information. The K-means clustering was conducted, which could be considered both as supervised and unsupervised cluster approaches. It aims to separate the objects into a preliminary chosen number of clusters so that each object belongs to the cluster with the nearest distance that represents the prototype of the cluster. This procedure is iterative.

Two-way clustering is also used in the present study. It is useful in cases where both objects and variables are expected to contribute simultaneously to discover meaningful patterns for the clusters. The resulting structure (clusters) is by nature not homogeneous. However, it is recognized by some researchers that this method offers a powerful tool for exploratory data analysis since it makes possible (like other multivariate statistical approaches such as correspondence analysis) to reveal relationships among the clusters of objects and the clusters of variables.

PCA (very similar to factor analysis) is a linear dimensionality reduction strategy, which introduces new directions in the original data space, making it possible to project the data in a lower-dimensionality space formed by new linearly uncorrelated variables called principal components or latent factors [20]. It is worth mentioning that the main application of this analysis is to reduce the number of variables and to detect structure in the data set, which allows classifying objects or variables. Therefore, PCA is applied as a data reduction method [19,20].

  1. Ethical approval: The conducted research is not related to either human or animal use.

3 Results

3.1 Hierarchical CA

The input data set consists of 16 borate glass samples (eight of them with Pb doping and eight with Zn doping) described by 11 parameters (physicochemical variables) [21].

The goal is to reveal patterns of similarity among the different glass samples or among their physicochemical descriptors.

In Figure 1, the hierarchical dendrogram for linkage of 11 variables is presented (z-transformed raw data; squared Euclidean distances as similarity measures; Ward’s method of linkage, and Sneath’s test of cluster significance).

Figure 1 Hierarchical dendrogram for clustering of 11 physicochemical parameters.
Figure 1

Hierarchical dendrogram for clustering of 11 physicochemical parameters.

Three clusters could be interpreted:

C1:Cg – packing density and CuO (%mol).

The essential role in glass properties is related to the atomic packing density (Cg). The atomic packing density is described as the ratio between the minimum theoretical volume occupied by the constituting ions and the effective volume of the glass. The correlation of Cg with the copper content is readily observed on the dendrogram. One could conclude that this tendency is general for the Cu-containing glass network. For CuO levels above 20 mol%, Cg tends to level-off, and no further accumulative behavior was detected (Supporting Information).

For the case of the glass system with lead, the packing density decreases for above 20 mol% of PbO. This observation is a consequence of the structural properties and the increasing significance of the packed structure pattern in the packing density.

C2: glass transition temperature Tg, K: elastic modulus, glass, 〈U0/V0〉 volume density of energy (equation (4) from SI); 〈U0/V0〉 volume density of energy (equation (5) from SI); E, Young’s modulus, G, Shear modulus, Poisson ration.

The similarity between glass transition temperature Tg and elastic modulus, glass K, is depicted on the dendrogram (Figure 1). The K is directly related to the external force required to compress or extend interatomic distances in opposition to the internal forces that seek to establish equilibrium interatomic distance.

The K is directly related to the external force required to compress or extend interatomic distances in opposition to the internal forces that seek to establish equilibrium interatomic distance.

The elastic properties of glasses define many of their mechanical properties. Besides, the thermal or pressure change of the bulk (K), shear (G), and Young (E) moduli offers information on the variations in the atomic glass structure. This cluster reveals a strong correlation among these physical properties. The values of the constants K, G, E as well as of the Poisson’s ratio could be measured experimentally by various methods for solids, including mechanical deformation or sound wave propagation techniques. For liquids, the elastic constants can be obtained from velocity measurements of high-frequency sound waves. These explanations are directly supported by the results of the hierarchical CA.

In the last cluster, C3, the similarity pattern is formed between density and coefficient of thermal expansion (CTE). This linkage could be used to model new Cu glasses with particular thermal characteristics. C3 is a very homogeneous cluster since C1 and C2 could be considered as one bigger cluster, which could be divided into two subclusters.

Figure 2 represents the hierarchical clustering of the glass samples (objects of the study). There is a well-expressed separation between Pb- and Zn-glass types. It could be stated that higher concentrations of Pb (PCB 20 and PCB 30) form an intermediate cluster between both classes of samples and are closer to the Zn type of glass rather than to Pb type of glass. This should be explained by specific relationships among the physicochemical parameters.

Figure 2 Hierarchical dendrogram for clustering of 16 glass samples.
Figure 2

Hierarchical dendrogram for clustering of 16 glass samples.

3.2 PCA

The z-transformed dataset was also subject to PCA. It was found that three principal components explain over 95% of the total variance of the system. Figure 3 shows the projection of the variables on the plane of factors 1 and 2, and Figure 4 shows the projection on the plane of the objects.

Figure 3 Projection of 11 variables on the plane described by factor 1 vs factor 2.
Figure 3

Projection of 11 variables on the plane described by factor 1 vs factor 2.

Figure 4 Projection of 16 cases on the plane described by factor 1 vs factor 2.
Figure 4

Projection of 16 cases on the plane described by factor 1 vs factor 2.

In general, most of the groupings observed by hierarchical CA are confirmed by PCA. It is readily seen that CTE and density variables form a well-expressed group since all other variables belong to another class.

The same conclusion is valid for the projection of the objects, that is, well-expressed separation of Pb and Zn cases with an intermediate position of PCB 20 and PCB 30).

The relationship between glass samples and physicochemical characteristics is demonstrated by a two-way clustering approach. It reveals how close the clusters of cases (glass systems) are connected with the clusters of variables (physicochemical parameters). Thus, discriminant parameters for each glass cluster could be extracted.

The foremost object for separation of the glass samples into two classes is the variances in the physicochemical parameters. Pb-glasses are characterized by the highest values of density and CTE since Zn-glasses show the lowest level of these two characteristics. Again, PCB 20 and PCB 30 differ from the Pb-glass pattern having low values of CTE and density.

All other physicochemical parameters reveal their higher levels for the Zn-glasses class and lower ones for the Pb-glasses class.

4 Conclusions

The present study has used intelligent data analysis to elucidate the relationship between the copper–lead and copper–zinc borate glasses with different copper contents. Several correlations among glasses and their physicochemical properties based on multivariate analysis techniques as CA, PCA, and two-way clustering were obtained. This allows for identifying the effects of Pb and Zn doping metal ions in copper glasses in the structural and mechanical properties. One cluster is identified for Zn-doped glasses, whereas two clusters with differentiated properties are identified with copper glasses doped with Pb.


tel: +359-81-61-799

Acknowledgments

The authors D. D. and M. M. are grateful for the financial support for this article by the Scientific Research Department (NIS) of Sofia University – 80-10-158/24.04.2020/NIS – Sofia University projects found 2020. The author M. N. is grateful for the additional support by the program “Information and Communication Technologies for a Single Digital Market in Science, Education and Security” of the Scientific Research Center, NIS-3317.

  1. Conflict of interest: The authors declare no conflict of interest.

References

[1] Metwalli E. Copper redox behavior, structure and properties of copper lead borate glasses. J Non-Cryst Solids. 2003;317:221–30.10.1016/S0022-3093(02)01853-7Search in Google Scholar

[2] Bae B-S, Weinberg MC. Crystallization of copper metaphosphate glass. J Am Ceram Soc. 1993;76:1395–400.10.1111/j.1151-2916.1993.tb03917.xSearch in Google Scholar

[3] Vasantharani P, Sangeetha N. Characterization of lead based binary and ternary glass systems using spectroscopic methods. Int J Res Pure Appl Phys. 2013;3:1–6.Search in Google Scholar

[4] Gautam C, Yadav AK, Singh AK. A review on infrared spectroscopy of borate glasses with effects of different additives. ISRN Ceram. 2012;2012:1–17.10.5402/2012/428497Search in Google Scholar

[5] Rajyasree C, Teja PMV, Murthy KVR, Rao DK. Optical and other spectroscopic studies of lead, zinc bismuth borate glasses doped with CuO. Phys B Condens Matter. 2011;406:4366–72.10.1016/j.physb.2011.08.082Search in Google Scholar

[6] Pal M, Roy B, Pal M. Structural characterization of borate glasses containing zinc and manganese oxides. J Mod Phys. 2011;2:1062.10.4236/jmp.2011.29129Search in Google Scholar

[7] Rao TR, Krishna CR, Thampy UU, Reddy CV, Reddy YP, Rao PS, et al. Physical and spectral investigations of Cu2+-doped alkali zinc borate glasses. Appl Magn Reson. 2011;40:339–50.10.1007/s00723-011-0208-6Search in Google Scholar

[8] Thulasiramudu A, Buddhudu S. Optical characterization of Cu2+ ion-doped zinc lead borate glasses. J Quant Spectrosc Radiat Transf. 2006;97:181–94.10.1016/j.jqsrt.2005.04.006Search in Google Scholar

[9] Ehrt D. Zinc and manganese borate glasses–phase separation, crystallisation, photoluminescence and structure. Phys Chem Glas-Eur J Glas Sci Technol Part B. 2013;54:65–75.Search in Google Scholar

[10] Krogh-Moe J. The infrared spectra and the structure of some anhydrous zinc borates. Z Für Krist - Cryst Mater. 1962;117:166–70.10.1524/zkri.1962.117.2-3.166Search in Google Scholar

[11] Harrison DE, Hummel FA. Phase equilibria and fluorescence in the system zinc oxide-boric oxide. J Electrochem Soc. 1956;103:491–8.10.1149/1.2430391Search in Google Scholar

[12] Möncke D, Kamitsos EI, Herrmann A, Ehrt D, Friedrich M. Bonding and ion–ion interactions of Mn2+ ions in fluoride-phosphate and boro-silicate glasses probed by EPR and fluorescence spectroscopy. J Non-Cryst Solids. 2011;357:2542–51.10.1016/j.jnoncrysol.2011.02.017Search in Google Scholar

[13] Ehrt D. Photoluminescence in the UV-VIS region of polyvalent ions in glasses. J Non-Cryst Solids. 2004;348:22–9.10.1016/j.jnoncrysol.2004.08.121Search in Google Scholar

[14] Hana NA, Taibi M, Haddad M, Aride J, Benkhouja K. Role of copper ions in dielectric and structural investigations of lead–borate glasses. Phase Transit. 2017;90:534–47.10.1080/01411594.2016.1227986Search in Google Scholar

[15] Kashif I, Ratep A. Effect of copper addition on BO4, H2O groups and optical properties of lithium lead borate glass. Opt Quantum Electron. 2017;49:231.10.1007/s11082-017-1067-7Search in Google Scholar

[16] Stefan R, Culea E, Pascuta P. The effect of copper ions addition on structural and optical properties of zinc borate glasses. J Non-Cryst Solids. 2012;358:839–46.10.1016/j.jnoncrysol.2011.12.079Search in Google Scholar

[17] Hristov H, Nedyalkova M, Madurga S, Simeonov V. Boron oxide glasses and nanocomposites: synthetic, structural and statistical approach. J Mater Sci Technol. 2017;33:535–40.10.1016/j.jmst.2016.07.016Search in Google Scholar

[18] Nedyalkova M, Hristov H, Simeonov V. Statistical approach to study of lithium magnesium metaborate glasses. Open Chem. 2017;15:61–6.10.1515/chem-2017-0008Search in Google Scholar

[19] Massart DL, Kaufman L. The Interpretation of Analytical Chemical Data by the Use of Cluster Analysis. New York: Wiley; 1983.Search in Google Scholar

[20] Massart DL, Vandeginste BG, Buydens LMC, De Jong S, Lewi PJ, Smeyers-Verbeke J, et al. Handbook of chemometrics and qualimetrics: part A. Appl Spectrosc. 1998;52:302A.10.1366/0003702981944896Search in Google Scholar

[21] Yao ZY, Möncke D, Kamitsos EI, Houizot P, Célarié F, Rouxel T, et al. Structure and mechanical properties of copper–lead and copper–zinc borate glasses. J Non-Cryst Solids. 2016;435:55–68.10.1016/j.jnoncrysol.2015.12.005Search in Google Scholar

Received: 2019-11-18
Revised: 2020-05-11
Accepted: 2020-06-04
Published Online: 2020-08-26

© 2020 Dimitar Dimitrov et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 28.3.2024 from https://www.degruyter.com/document/doi/10.1515/chem-2020-0140/html
Scroll to top button