• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Chemical and physical studies of metallic alloy-based old Indian coins with LIBS coupled with multivariate analysis

    2023-10-08 08:20:56VikasGUPTAAbhishekhKumarRAITejmaniKUMARRajendharJUNJURIManojKUMARandRAI
    Plasma Science and Technology 2023年9期

    Vikas GUPTA,Abhishekh Kumar RAI,Tejmani KUMAR,Rajendhar JUNJURI,G Manoj KUMAR and A K RAI,?

    1 Laser Spectroscopy Research Laboratory,Department of Physics,University of Allahabad,Prayagraj 211002,India

    2Department of Earthand Planetary Sciences,University of Allahabad,Prayagraj 211002,India

    3AdvancedCentre ofResearch in High Energy Materials,University of Hyderabad,Hyderabad 500046 India

    4 Leibniz Institute of Photonic Technology,Jena D-07745,Germany

    Abstract The present work aims to demonstrate the capabilities of Laser-induced Breakdown Spectroscopy(LIBS) coupled with a multivariate technique for rapid quantification and classification of old Indian coins made of various alloys.Thirteen old Indian coins in different years of circulation,(1922-1986) were selected for the study.The concentrations were determined by Calibration free LIBS(CF-LIBS) method.The concentration of cuprum(Cu) is negligible,and aluminum(Al)is maximum in the first five coins,and vice-versa in the remaining eight coins.Two different multivariate methods,Principal Component Analysis(PCA)and Soft Independent Modelling of Class Analogy(SIMCA) have been used to classify and identify the coins.PCA classified all thirteen samples into four main alloy categories.The discernment of unknown samples to their probable class membership of alloy was performed using SIMCA.The surface hardness(Brinell hardness number) is linearly correlated with the plasma temperature and LIBS intensity ratios.The sample surface of the first and fifth coin belongs to Al-alloy,having the least surface hardness,and it became harder for Cu-Ni alloy,Ni-brass alloy,and bronze alloy.The hardness of the surface is more for bronze sample twelve.It is also observed that the plasma temperature increases monotonically with the Brinell hardness number.This analysis provides valuable information on fabrication methodology and explains large diversification in the elementary composition of old coins.

    Keywords:LIBS,metallic alloy,stoichiometric ablation,PCA,SIMCA,old indian coins,surface hardness

    1.Introduction

    We observed extensive diversification in the elementary composition and fabrication methodology in Mughals/British coins and Indian rupee coins during India’s independence[1-4].During that period,both decimal and non-decimal coins are in circulation.Before India’s independence,British Indian coins were categorized/named as non-decimal coins.After India’s independence,British Indian coins were used as a frozen currency until India became a republic in 1950.The first rupee coins of the Republic of India were minted in 1950.These coins are referred to as Anna series or pre-decimal coinage,which includes a half rupee,quarter rupee,two annas,one anna,half anna and one pice coins.In the Annaseries,one rupee was divided into 16 annas or 64 pice,with each anna equal to 4 pice.In 1957,India shifted to the decimal system.After adopting the decimal system,the coins minted between 1957 and 1964 were printed with the legend‘Naya Paisa’(New Paisa).With the Indian Coinage(Amendment) Act of 1955,India shifted to the ‘Decimal System’,according to which a new ‘decimal series’ was introduced in 1957.The rupee was now divided into 100 paise instead of 16 annas.These newly introduced coins were initially struck with the words ‘Naya Paisa’ to distinguish them from the previous coins,but the word‘Naya’was dropped in 1964.A steady increase in pricing gradually led these coins to be minted only in Al from 1964,and these new Al series dropped the word‘Naya’from it.Later a new metal was introduced,and we soon had Ferritic Stainless Steel coins of 10,25,and 50 paise introduced in 1988 and one rupee in 1992.Due to a considerable decline in demand for coins of lower denominations and an increase in metal costs,circulation of lower denomination coins was reduced,and the Government of India,in consultation with the Reserve Bank of India,decided to mint all denomination coins in ferritic stainless steel in 2006 [1].These coins,categorized as Anna series(NDC),Naya Paisa(Frozen Currency/DC),and Paisa(NDC) were mainly made of the composition of metals/alloys like copper-nickel(cupronickel),copper-nickel-zinc(nickel-brass),bronze and aluminium-magnesium [1,2].

    The composition of metals/alloys,in ancient times as well as during India’s independence was very expensive and highly valued materials.They had a face value at the minimum size and weight,that is why many coins were made of gold or silver or copper,or iron in that period.Therefore,the first and most important is that metals/alloys are more resistant and can stand long travels and use.The second is that in ancient times the value of coins was their face value.There were no central banks to back up what they said their worth was.The coins that are planned for circulation should have some special properties/requirements like the ease of stamping,anti-corrosion properties,and high durability so that they can be used for a long time.Achieving this goal has tended to necessitate the use of alloys.Most of the coins during India’s independence are made of mainly aluminummagnesium alloy and copper alloy(cupro-nickel,nickelbrass,and bronze).Metals,like copper,silver,gold,and iron were also used in ancient times to make coins because of their hardness and high durability.Iron may have been replaced because of its usage in war equipment and its tendency to get corroded easily.The different compositions of metals or alloys have been widely used to make coins.Other metals like gold and silver were replaced because of the variation of their face value during different periods.Most of the coins during India’s independence were made up of copper.Brass,the composition of copper and zinc with a 1:3 ratio may be chosen because of its attractive gold color with high durability and anti-corrosion properties[5].Also,it is relatively cheaper and stronger than copper.The coins made by the melting copper along with nickel or tin were also used to make coins known as copper-nickel and bronze.The use of copper in most of the coins is because of being a naturally antibacterial material.The durability of copper-nickel coins is more than that of other alloys,that is why it was used to make coins.Bronze has also been used to make coins.Bronze is harder than copper due to alloying that metal with tin or other metals.Bronze is also more fusible(i.e.more readily melted) and is hence easier to cast.It is harder than pure copper and far more resistant to corrosion.We observed that ancient coins were made of gold,silver,copper,or alloys containing a high percentage of these three elements.The value of a coin was the market value of the metal in it.Modern coins represent a promise from the issuing government rather than the coin’s actual value.It becomes necessary to change to cheaper metals when the metal in a coin costs more than the actual coin.At that time,Al was also a suitable metal to make coins that had low cast and lightweight.Also,we could make more coins than copper or zinc from the same process with the same amount of material with the least cast.Also,lighterweight of coins could make them easier to carry.Metals/alloy were also used as a coin because it is hard and difficult to break.It has also been observed that in ancient periods/during Indian independence,bronze was a hard metal than copper,which has been widely used to make coins and artifacts [6-10].

    To get an idea about the dynasty of the different periods(Mugal,British,and Indian),it is of utmost necessity to identify and quantify the major and minor constituents of the coins.Therefore,we aim to present an example for quick identification,quantification,and classification of various currencies/coins using Laser-induced Breakdown Spectroscopy(LIBS) coupled with multivariate methods.

    Different analytical techniques like inductively coupled plasma-atomic emission spectrometry(ICP-AES),x-ray fluorescence spectroscopy(XRF),inductively coupled plasma-mass spectrometry(ICP-MS),and Energy Dispersive x-ray Analysis(EDX) are used for the elemental analysis of metallic objects.These techniques are more promising and give precise discernment of different elements/alloys used to prepare the coins [11-19].However,these techniques need sample preparation or involve solvent/dissolution.Sometimes it is impossible to take the sample to the laboratory due to their size or their surrounding conditions.Also,sometimes the things are so precious that they need analysis without damaging and breaking them.The coins associated with the surrounding objects on the wall needin situcharacterization.Therefore,we need a profound technique to overcome all these drawbacks in the above-discussed analytical tools for rapidly identifying and classifying metallic alloy-based objects.Recently,LIBS has emerged as a more effective tool for compositional analysis and identification/classification of materials/alloys used in producing different objects/coins[20-30].Because of its ability ofin situ,quasi-non-destructive,and no/minimal sample preparation,we can use this technique to identify and classify metallic alloys in different fields like archaeology and precious metallic alloys/coins.In addition to that,the LIBS can establish the correlation of mechanically measured surface hardness of the different coins with plasma parameters and intensity ratio of ionic and atomic lines [22,26,27].The conventional techniques such as theVickers hardness test,the Brinell hardness test,etc are destructive and time-consuming compared to LIBS.In the present work,we have measured the surface hardness of some coins with LIBS technique and compared it with the result found by the Brinell hardness number(BHN).

    Table 1.Obverse and reverse view of old Indian coins with their composition.

    The combination of LIBS with multivariate methods gives a quick and simple methodology for identifying and classifying metallic alloys,which is cost-effective and instant.Principle component analysis(PCA) and Soft Independent Modelling of Class Analogies(SIMCA) are well-known multivariate methods used for the classification and identification of samples [31-33].It extracts valuable information from a large LIBS data set.Multivariate analysis methods are used to evaluate and collect statistical data to clarify and explain relationships between different variables associated with this data.The above technique of LIBS is performed not only as a spectrochemical analytical technique but also as an easy way to measure the surface hardness of solid targets(metallic coins) simultaneously.

    2.Material and methods

    Thirteen old Indian coins have been selected for the study made up of two main categories of alloys,Cu alloy and Al-Mg alloy.The Al-Mg alloy is composed of mainly aluminum and 2%-3% magnesium where five coins belong to this category(alloy 1).Another five coins belong to cupro-nickel alloy(alloy 2).In the remaining three coins,two coins are composed of nickel-brass alloy(alloy 3) and one is bronze(alloy 4).Generally,cupro-nickel composed of Cu and Ni as main matrix elements with a ratio of 3:1 yields white colour.Similarly,in nickel-brass,the composition of Cu-Ni-Zn is 75%,20%,and 5%,respectively,yielding golden yellow colour.The appearance of cupro-nickel is nearly the same as with Al-Mg alloy [5].The composition of the selected coins was given on the RBI website [1].

    The photographs and their composition of the 13 Indian coins are presented in table 1.Generally,these coins belong to different alloys like Al-Mg alloy,Cu-Ni alloy,Ni-brass,and bronze alloy,etc.Also,they were classified as Decimal(DC)and Non-decimal(NDC)coins.All the coins except c6,c7,c8,and c11 belong to NDC coins,and the remaining belong to DC.

    2.1.LIBS instrumental setup

    The LIBS spectra were recorded by accumulating five laser shots to get one spectrum.Thirty-five sets of LIBS spectra of each sample have been acquired by using a Q-switched Nd:YAG laser(Innolas split light 1200).The schematic diagram of the LIBS experimental setup is shown in figure 1.A planoconvex quartz lens with a 150 mm focal length is used to focus a pulsed laser beam at the sample’s surface.The laser delivered energy of 45 mJ,having an irradiance of nearly 1012W cm-2at the focal spot,generating a fugacious plasma.The emitted signal from plasma is recorded with the help of a Mechelle spectrometer(ME-5000) coupled to ICCD(Andor,iSTAR DH-734),having a 0.1 nm spectral resolution at 500 nm.The motion in theX-Ydirection of the sample stage is controlled with a motion controller(EPM 2000),which provides a fresh surface for each laser exposure.The surface contamination is removed by using defocused laser shots on the surface of the coins by moving the coins in a controlled manner using a suitable translational stage.Then the spectral signature was taken at the same fresh spot position by retracing.The other experimental parameters like the focal length of the lens,energy of the laser beam,lens to sample distance,repetition rate,etc are optimized to get better LIBS spectra.These experimental parameters are maintained the same during the whole experiment.

    Figure 1.Schematic diagram of LIBS instrumental setup.

    2.2.Estimation of surface hardness

    In this section,we have used the Brinell hardness test to estimate the surface hardness of coins in terms of Brinell Hardness Number(HBN) using the formula given in equation(1).It is commonly used to determine the surface hardness of metals and alloys.In this method,we have indented a hardened steel or carbide ball under a load into the sample surface.The HBN is obtained by dividing the load,in kilograms,by the spherical area of the indentation in square millimeters;this area is a function of the ball diameter and the depth of the indentation

    whereF-applied load(kg),D-indenter diameter(mm),andd-indentation diameter(mm).Here,we have taken 187.5 kg equivalent load for all samples,and 2.5 mm of indenter has been used.The impression of the indenter on the coins surfaces c1,c9,c11,and c12 has been shown in figure 2,and the hardness of the four coins c1,c9,c11,and c12 is listed in table 2.

    3.Results and discussion

    3.1.Qualitative analysis of the samples

    The LIBS spectra of the coins in the UV-NIR spectral region are shown in figure 3 in the stacked form,providing the identification and discrimination of the samples based on the appearance of the spectral signature of significant elements.We have used the atomic spectra database(ASD) of the National Institute of Standards and Technology(NIST) for the assignment of the spectral lines of various elements in the LIBS spectra of the coins [34].

    Figure 2.Impression of indenter on the surface of the coins c1,c9,c11,and c12.

    Table 2.The samples with their respective BHN.

    Figure 3 shows the LIBS spectra of thirteen old Indian coins(c1-c13) in the spectral range of 240.0-595.0 nm.As shown in figure 3(a),a similar spectral signature in the LIBS spectra is observed for the samples c1-c5,which belong to only decimal coins(DC).Another group of coins(c6-c10)combination of decimal and non-decimal coins(Anna series)shows nearly similar LIBS spectra shown in figure 3(b).The LIBS spectra of the remaining three samples(c11-c13)show minor dissimilarity in spectral signature presented in figure 3(c).The spectral lines of Zn near 330.2 and 472.3 nm were not identified in sample c12,while it is present in samples c11 and c13.In contrast,the spectral signature of Sn near 283.9 nm is identified in sample c12,while it is absent in samples c11 and c13.Thus,it would look simple to classify the various coins based on the spectral signature of the major elements,as shown in figure 3.After analyzing the whole spectral range(220.0-875.0 nm) of LIBS spectra,we have seen that most of the major peaks were observed in the spectral region of 240.0-595.0 nm.Various major peaks of significant elements responsible for categorizing the samples in a different group of alloys have been observed in this spectral region.We have observed mainly the spectral lines of Al at 256.8,265.0,305.0,308.2,and 394.4 nm and those of Al II at 263.1,281.6 nm,and Mg at 279.1,and 517.1 nm in the LIBS spectra of the first five samples(c1-c5) with higher spectral intensity.The LIBS spectra of the remaining ten coins(c6-c13)clearly show the presence of spectral lines of Cu as significant elements with minor differences in spectral intensities of different elements.Copper is a common element that is found in all samples except c1-c5.We have identified the persistent lines of copper at 324.7,327.3,510.5,515.3,and 521.8 nm and the strong lines of copper as major peaks at 570.0 and 578.2 nm.We have observed the Ni and Zn spectral lines as major peaks at 251.1,310.18,341.50,352.4,361.9,334.5,472.2,and 481.0 nm,respectively.Spectral lines of Ni with Cu were also identified as significant peaks in the coins c6-c10.The coins c11 and c13reveal the presence of Ni and Zn,along with Cu as major peaks.While in coin c12,spectral lines of Zn and Sn were observed as major and minor peaks,respectively,along with Cu.The spectral signature of Sn is present in sample c12,shown separately in figure 3(d),which is not observed in samples c11 and c13.Al,Ca,Mg,and Na were found in all samples with different spectral intensities as major peaks in the samples c1-c5,while they were seen as minor peaks in all remaining samples.The spectral signatures of Ca and Na were observed as minor peaks in all coins.Therefore,we may conclude that different alloys were used to make these coins.The wavelengths of the spectral lines of various elements found in a different group of samples are summarized in table 3.

    Figure 3.LIBS spectra of different coins(c1-c13) in the spectral range of 240.0-595.0 nm.

    Figure 4.Comparison of the intensity of the spectral lines of the elements in different samples.

    Figure 5.Comparison of spectral intensity of elements(Cu,Ni,Zn,and Sn) present in the sample.

    Figure 6.Comparison of the spectral intensities of the sample c10,c11 and c12(a)for Cu at 578.3 nm(b)for Ni at 341.5 and 342.4 nm and(c) for Zn at 472.3 and 481.1 nm.

    Figure 7.The Boltzmann plot for the Cu I and Al I.

    Figure 8.The variation of concentration of Cu and Al in coins c1-c13.

    After analyzing the spectral lines of the elements listed in table 3,we may conclude that all samples can mainly be divided into four groups of alloys.The first five coins belong to Al-Mg alloy,and the remaining coins belong to Cu alloy.Further classification among the different Cu alloys is challenging without comparing the spectral intensity of the elements present in the samples.Therefore,a comparison of spectral intensity using a bar plot has been shown in figure 4.Single spectra have been taken from each sample to compare spectral intensities.

    It is clear from the bar plot shown in figure 4 that the relative concentration of Al and Mg is higher in the first five samples belonging to alloy one compared to other samples.For other samples belonging to Cu-alloy(alloy 2),the relative concentration of Cu and Ni is approximately the same in samples c6-c10.For further classification of three coins(c10,c11,and c12) we have compared the relative intensity of spectral lines of Cu,Ni,Zn,and Sn,shown in figure 5.

    Table 3.Wavelengths of the different spectral lines of major and minor elements present in samples.

    Based on the relative intensity of these elements,we can conclude that sample c10 belongs to the copper-nickel alloy and sample c11 belongs to the Ni brass alloy.The sample c12 is different from c10 and c11 due to the presence of Sn.Therefore,it is clear that sample c12 belongs to bronze alloy.The relative spectral intensities of Cu and Ni at wavelengths 578.3 nm,341.5 nm,and 343.4 nm,respectively,are shown in figures 6(a)-(b).Also,the relative spectral intensity of Zn is shown in figure 6(c).

    3.2.Quantitative analysis

    To verify the result of qualitative analysis,we have calculated the concentration of each species in the samples using the CFLIBS [35-37].The CF-LIBS method does not need any reference materials.Therefore CF-LIBS has been used in this work for quantitative analysis.For quantitative analysis,we have looked over Plasma’s temporal evolution.The integrated intensity of the spectral lines of the elements in the LIBS spectra plays a vital role,as the concentration of the elements present in the sample is directly proportional to the intensity of the spectral line given in equation(2).The following equation(2) is the line intensityInm,of a transition in LIBS spectra.The suffix inInmindicates the transition from upper energy levelEnto lower energy levelEm

    whereAnmis the transition probability for the transitionEntoEmandgnis the degeneracy parameter forEn,kBis the Boltzmann constant,CsandUs(T) are the relative concentration and partition function of that species,respectively,andFis an experimental factor.But before using the intensities of spectral lines for the analytical purpose,the laserinduced plasma should fulfil the following three conditions.

    (i) Laser ablation should be stoichiometric.

    (ii) The plasma should be optically thin.

    (iii) Plasma should be in local thermal equilibrium(LTE).

    To satisfy the first condition,the irradiance of the laser shot at the focus spot at the target sample(coins) should be higher than 1012W cm-2[27,35].This is a sufficient condition for stoichiometric ablation.Irradiance is dependent on the radiant flux of Nd:YAG laser pulse and the energy received by the area of the focal spot having a diameter(D)given in the equation below

    where,λanddare the wavelength and diameter of the laser beam,andfis the focal length of the convex lens used to focus the laser beam on the target sample.We calculated irradiance for the present experiment and found near 5.0×1012W cm-2,which is much greater than 109W cm-2.Therefore we have achieved the criteria of the stoichiometric nature of the ablated part to proceed with further analysis.

    For optically thin plasma,the theoretical value of the intensity ratio of two atomic lines of the same species should equal the experimental intensity ratio.These intensity ratios of two emission lines corresponding to the same species and having the nearly same upper energy level can be representedby the following equation

    whereInmandIijare the integrated spectra line intensities i.e.the experimental values corresponding to the different transitions of the same species.

    For the spectral lines of Cu at 510.5 nm and 570.0 nm having the nearly same upper energy level,we have found that the experimentally estimated intensity ratio(9.23)is quite near to the theoretically calculated intensity ratio(9.30) at gate delay near 1.5μs.A similar calculation was performed for all coins and found that experimentally determined intensity ratios are in good agreement with the theoretically calculated intensity ratio for the spectra recorded between 1.0 and 2.0μs gate delay.Thus,the delay time of 1.0-2.0μs is the best temporal window for quantitative analysis.Therefore,we have used this window for quantitative estimation of the elements.

    In the present experiment,the McWhirter criterion has been used to fulfill LTE as it is an essential necessary condition for LTE.A detailed description has been described in earlier publications [10,22,25,38].

    3.2.1.Elemental composition of the coins using CF-LIBS.After satisfying the three conditions,we applied the CF-LIBS algorithm [36,37].We have estimated the concentration of each element present in the sample with the help of the CFLIBS algorithm and the results are summarized in table 4.The Boltzmann plot for the species(Cu I and Al I) is shown in figures 7(a) and(b).Al is present in the first five coins(c1)-(c5) while Cu is present in last eight coins(c6)-(c13).

    3.2.2.Variation of the constituents in coins based on designation.The samples named c1-c5,c9,c10,c12,and c13 belong to decimal coins(DC),and the remaining samples,c6-c8,and c11 belong to non-decimal coins(NDC).Therefore,out of 13 coins,nine coins belong to DC.The first five coins are made of mainly Al and Mg.The concentration of Al varies from(95.75±0.23)% to(97.11±0.19)%.From table 4,we observed that copper alloys like bronze,Ni-brass,and cupro-nickel had also been used to make DC.For four coins(c6-c8 and c11)belonging to NDC,two different alloys Cu-Ni alloys(c6-c8) and Ni-brass alloy(c11),have been used to make these coins.The variation of Cu in Cu-Ni alloy varies from(68.70 ± 0.10)% to(76.23 ±0.03)% and that of Ni from(21.27 ± 0.10)% to(27.58 ±0.13)%.The concentration of Cu is higher in c11 as compared to c6-c8 coins belonging to Cu-Ni alloy.Therefore,we observed no significant pattern in the manufacturing process of decimal and non-decimal coins.Although,a stronger negative correlation(negligible concentration) was observed in the group of coins c1-c5 and c6-c13 based on the concentration of copper and aluminum,respectively shown in figure 8.

    3.2.3.Variation of the constituents in the coins based on alloy.After analyzing the composition of the coins,weconclude that these coins belong to mainly aluminum and copper alloys.The samples of the decimal coins from c1 to c5 belong to Al-Mg alloy,while the coins from c6-c10 and c11,c13,and c12 belong to different alloys like copper-nickel,brass,and bronze,respectively because of the relative concentration and composition as discussed in section 3.2 and also from results of CF-LIBS mentioned in table 4.

    Table 4.Concentration(%) of the elements in the coins obtained by CF-LIBS.

    3.2.4.Year-wise variation.Looking over the manufacturing year,we observe that these coins had been minted in different years mentioned in table 1.The coin c8 is the oldest coin out of all thirteen coins,while c4 is the latest one.Figure 9(a)shows the variation of Cu in various coins manufactured in different years.We observed that the coins minted up to the year 1962 contain a higher amount of Cu except coin c2.While the concentration of aluminum was increased in the coins minted after 1962.After that,copper was not used to minting coins.This indicates the change in the economic condition of India during that period.It seems that the manufacturing of coins shifted from copper to aluminum after almost 1962[1].It is clear from figures 9(a)and(b)that most of the decimal coins are made of mainly Al and belong to Al alloy.In contrast,most of the coins in the category of nondecimal coins were made of primarily copper.The manufacturing procedure of coins belonging to Cu-Ni alloy looks the same as the composition of decimal coins c9 and c10,and non-decimal coins c6,c7,and c8 are nearly the same.However,the issuing dates of these two groups of coins are very different.

    Figure 9.The year-wise variation of concentration of Cu and Al in coins c1-c13.

    Figure 10.The intensity ratio of ionic and atomic spectral lines versus the BHN for different alloy samples.

    Figure 11.Variation of plasma temperature with BHN.

    Figure 12.The scores plot(a) PC1 against PC2 and loading plot((b),(c)) of PCA results.

    3.3.Surface hardness measurement

    We have chosen four coins from each group of alloys to determine the surface hardness.We have used the Brinell hardness test to measure the hardness of the coin discussed in the material and method section,and the results are summarized in table 2.We have also compared the above results with the result obtained by the LIBS discussed in the following section.

    3.3.1.Correlation of intensity ratio of the of ionic and atomic spectral lines with the surface hardness of the material.We have tried to find a correlation between surface hardness measured by the LIBS(measuring the intensity ratio of ionic to atomic spectral lines) using Ca and Cu lines with the mechanically mustered surface hardness(Brinell hardness number).In this section,we have correlated the surface hardness of four samples belonging to Al and Cu alloy by measuring the intensity ratio of the ionic to atomic spectral lines.Measurement of the intensity ratio of the spectral lines requires a proper selection of well-resolved and intense spectral lines.For this purpose,we have identified ionic spectral lines of Ca II at 393.4 nm and atomic spectral lines of Ca I at 422.7 nm because these spectral lines are present in allsamples and well resolved.In addition to this,we have also used the ionic and atomic spectral lines of Cu at 237.0 nm and 578.7 nm,respectively for the samples belonging to the copper alloy for the determination of the surface hardness,which gives the best correlation with the mechanically measured hardness as compared to Ca lines.

    We observed from the graph shown in figure 10 that both have a positive correlation.A best linear fit was found for the copper lines shown in figure 10(b) with an RMSE of nearly 0.93.

    3.3.2.Correlation of the plasma temperature with the surface hardness.In this section,we have estimated the plasma temperature in LTE.With the help of equation(2),we can easily determine the plasma temperature for any species like Cu I,Al I by taking the logarithm of equation(2)

    whereInmis the integrated intensity of that species for different transitions and others are spectroscopic parameters corresponding to each transition.Csis the concentration of species andUs(T) is the partition function [36,37].By plotting the left-hand side versusEn,the slope(-1/kBT) of the obtained line gives plasma temperature.The emission lines of Cu(I) and Al(I) have been used to determine temperature under the condition of LTE for the coins c1 and coins c6,c9,c11,and c12 respectively.According to these requirements,the Al I and Cu I lines selected to determine the plasma temperature are listed in table 5 with their required parameters [34].We found that the plasma temperature varies from 11050 to 12450 K of coins c1,c9,c11,and c12,respectively,belonging to a different alloy.

    The linear relation between plasma temperature and hardness for the four coins of different alloy samples has been shown in figure 11.It is noticed here that the plasma temperature increases as the hardness of the metallic alloy increases.

    The probable explanation for the increase in surface hardness with the rise in plasma temperature is given below.

    (1) The breakdown threshold of the materials decreases with the increase of the surface hardness[26,27].Therefore,a large amount of the preceding laser beam pulse is utilized to heat the generated plasma,which increases plasma temperature.

    (2) This increase is attributed to the increase in the collision force.

    (3) Cowpeet al[38]have studied bioceramics samples using the LIBS to see the relationship between plasmatemperature and the sample hardness,and they also demonstrated a linear relationship between sample surface hardness and plasma temperature.

    Therefore,we conclude from analysis that the increasing order of hardness belongs to the material/alloys used to make coins.The sample surface of c1 and c5 belongs to Al-alloy,having the least surface hardness,and large for Cu-Ni alloy,Ni-brass alloy,and bronze alloy.The hardness of the surface is more for bronze sample c12.

    3.4.Classification using multivariate methods

    It is time-consuming to classify/discriminate samples/coins using the qualitative analysis of the spectral lines of the LIBS spectra.However,the identification and classification may also be made without any extended spectral analysis of each coin.Also,there is no need for a comparative study of spectral line intensity for the classification.With the help of the multivariate analysis of the LIBS data set,we can easily classify and identify the alloy used to manufacture a particular coin.A detailed analysis is given in the following section.

    3.4.1.PCA analysis.We have recorded 35 spectra for each coin,which account for a total of 455 sets for thirteen samples.LIBS spectra cover the spectra range of 220-875 nm,a large number of variables as the intensity at a different wavelength.Reducing the number of variables is the first step before performing the multivariate methods[31,32].It is possible by observing raw data in a simple graph.For statistical analysis of LIBS data,spectra obtained from different coins are arranged,and various data sets are prepared in the form of a matrix.We have created matrices of dimensions(455×145 32)made up of variables with 14532 features as columns.The 35 spectra taken from each sample have been represented by rows having dimensions 455 for 13 samples(c1×c13).The Unscrambler-X software(CAMO Software India Pvt.Ltd) is used to pre-process and solve the data matrix,which firstly makes covariance matrices and then further gives the scores and loading plots.

    The sample classification on the PCA score plot is shown in figure 12(a).We have found classified data on the PC1 versus PC2 graph.The PC1 discriminates the data with 92% of the total variance and the PC2 with 5% of the total variance.There is no need for a higher PC as we have covered significant variance for classification.The score plot of PCA visually represents the presence of different clusters,which implies that different alloys have been used to make coins.Each point on the PCA score plot represents the single spectra of one sample.Therefore,we found spectral data of all samples clustered in a different region,with their groups describing the similar composition of the coins.The PC1 component reveals the significant dissimilarity between two groups of coins,mainly coins made of Al-Mg and Cu alloys,which were classified on the positive and negative sides of the PC1.The responsible transitions of Al(3p2P1/2,3/2-4s2S1/2)and Mg(3s3p3P012-3s4s3S1) at wavelength 394.4 nm,396.1 nm for Al and Mg at 516.7 nm,517.2 nm,518.3 nm for the classification of these two groups were easily identified with a loading plot of PC1 shown in figure 12(b).

    The LIBS spectra of the coins(c1-c5) belonging to Al-Mg alloy are clustered on the positive side of the PC1 axis.All other LIBS spectra of the coins(c6-c13)belonging to Cu alloy clustered on the negative side of the PC1 axis.Further,the LIBS spectra of the coins made of the copper alloy have been classified on the PC2 axis in two different groups having major similarities with minor differences.In the third quaordinant on the PCA score plot,we observed the clustering of LIBS spectra of the coins c6-c10 while the LIBS spectra of coins c11-c13 clustered in the first quaordinant of the PCA score plot i.e.; on the positive side of PC2.However,we observed two groups of clustering in the first quaordinant on the score plot of PCA.The loading plot of PC2 identified the spectral lines which are responsible for the further discrimination of the coins belonging to Cualloy.Some of the major responsibilities,strong transition lines of different elements like Cu(4p2P3/2,1/2-4s22D3/2),Zn(4s4p3P1,2-4s5s3S1)were easily identified on the loading plot of PC2 at wavelengths of 570.0 nm,578.2 nm for Cu and that of Zn at 472.2 nm and 481.0 nm shown in figure 12(c).The coins belonging to Alloy 2 and Alloy 3 seem to be very similar on the PCA score plot.It was clustered into different groups.

    3.4.2.SIMCA analysis.We have performed SIMCA for the supervised classification of data and identification of a new sample with their group of alloy[33].For this,We have taken the LIBS data set of four coins(c1,c10,c11,and c12),one from each alloy,to create a model.Out of 35 LIBS data sets of each coin,we have selected 34 as a trainee and one LIBSdata set has been kept out for a test for each alloy.For each of them,four different PCA models were created for each group of alloys by selecting the different regions on the score plot of PCA.This multivariate method calculates the principle components which are the linear combinations of the variables that are associated with most of the valuable information.Table 5 shows the distance of each PCA model for different groups of alloys(alloys 1-4).The similarity and dissimilarity have been easily observed by analyzing the relative distance of all PCA models.There is a significant distance between alloy 1 and with other three classes of alloy corresponding to Cu alloy.However,the model distance of others is approximately small and belongs to nearly the same class of Cu alloy.Table 5 of model distance shows that there is much more similarity in alloy 3 and alloy 4 as they have most of the spectral signatures common.However,the dissimilarity is revealed with the model distance of alloy 2 to those of alloy 3 and alloy 4.Now,we supposed the ‘test’LIBS data set of the same coins as four unknown samples which have been selected for identification/classification named s1,s2,s3,and s4.Another LIBS data set of a sample belonging to Nordic gold has been intentionally chosen for cross-validation named s5.The projection of a new sample corresponding to their class member has been listed in table 6.All samples have been correctly classified except sample s5,a Nordic gold alloy,which has been intensionally chosen for the cross-validation of SIMCA classification.It reveals dissimilar samples belonging to none of the sample class alloys in this model.

    Table 5.PCA model distance.

    Table 6.Predicted results of a classification in group.

    4.Conclusions

    LIBS can simultaneously analyze,classify and measure the surface hardness of different samples,especially metallic alloys.Emission spectra of thirteen old Indian coins made up of various alloys were recorded without any sample pretreatment.The presence of the characteristic spectral lines of the various elements in the LIBS spectra of the coins reveals that the studied coins are made up of four different alloys.The intensity of the spectral lines of different elements was measured to compare the relative concentration of the elements/constituents in each coin.CF-LIBS approach is used to simultaneously determine the concentration of all constituents/elements of each coin studied in the present work.The concentration of Al in coins c1-c5 was estimated and found to vary from(95.75±0.23)% to(97.11±0.19)%.We found a higher concentration of Cu((91.18 ± 0.16)%) in bronze sample c12,defined as one Naya paisa(DC).The measurement of the surface hardness was demonstrated by correlating the surface hardness with the LIBS intensity ratio.Additionally,the plasma excitation temperature has been estimated for each coin and correlated with the surface hardness of different samples.We found one Naya paisa coin(c12) belonging to bronze alloy is harder than the others.Principle Component Analysis(PCA)has visually shown the classification of the coins into four groups.Soft Independent Modelling of Class Analogies(SIMCA) was used to identify the categories of new samples based on the model developed in the present work.We observed precise classification for all four classes of alloy by which the coins were made off.Finally,a suitable multivariate model for the category of a group of a sample belonging to four different alloys is developed.

    Acknowledgments

    Mr Vikas Gupta is thankful to UGC-CSIR for providing financial assistance as a junior research fellow for the research work.

    久久人妻福利社区极品人妻图片| 如日韩欧美国产精品一区二区三区| 日本vs欧美在线观看视频| 老司机午夜福利在线观看视频| 法律面前人人平等表现在哪些方面| 精品国内亚洲2022精品成人| 国产精品一区二区精品视频观看| 国产av在哪里看| 少妇被粗大的猛进出69影院| 国产精品av久久久久免费| 真人一进一出gif抽搐免费| 国产精品久久电影中文字幕| 久久精品亚洲熟妇少妇任你| 97碰自拍视频| 91精品三级在线观看| av视频在线观看入口| 国产成人一区二区三区免费视频网站| 中文字幕av电影在线播放| 黑人巨大精品欧美一区二区mp4| 在线视频色国产色| 黄频高清免费视频| 精品久久蜜臀av无| 两性午夜刺激爽爽歪歪视频在线观看 | 两个人免费观看高清视频| 久久人妻福利社区极品人妻图片| x7x7x7水蜜桃| 午夜日韩欧美国产| av天堂在线播放| 欧美日本中文国产一区发布| 国产色视频综合| 国产精品 欧美亚洲| 国产免费男女视频| 天天躁狠狠躁夜夜躁狠狠躁| 极品教师在线免费播放| 日本vs欧美在线观看视频| 中文字幕av电影在线播放| 妹子高潮喷水视频| 国产私拍福利视频在线观看| 99riav亚洲国产免费| 精品久久蜜臀av无| 成年版毛片免费区| 久久精品aⅴ一区二区三区四区| 精品国产乱码久久久久久男人| 人人妻人人澡人人看| 9色porny在线观看| 每晚都被弄得嗷嗷叫到高潮| 99国产精品免费福利视频| 麻豆av在线久日| 午夜福利成人在线免费观看| 久久精品亚洲熟妇少妇任你| 免费少妇av软件| 叶爱在线成人免费视频播放| 十分钟在线观看高清视频www| 波多野结衣巨乳人妻| 亚洲精品国产区一区二| 国产一区二区在线av高清观看| 欧美乱妇无乱码| 午夜成年电影在线免费观看| 一级黄色大片毛片| 久久中文看片网| 亚洲五月天丁香| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 美女免费视频网站| 国产成人精品久久二区二区91| 欧美色欧美亚洲另类二区 | 51午夜福利影视在线观看| 亚洲人成伊人成综合网2020| 久久久久久久久中文| 久久久久久人人人人人| 亚洲av电影在线进入| 久99久视频精品免费| 脱女人内裤的视频| 亚洲精品在线观看二区| 51午夜福利影视在线观看| 国产单亲对白刺激| 成年版毛片免费区| 在线免费观看的www视频| 十八禁人妻一区二区| 国产精品永久免费网站| 国产99白浆流出| 亚洲成av片中文字幕在线观看| 18禁裸乳无遮挡免费网站照片 | 男女做爰动态图高潮gif福利片 | 欧美乱码精品一区二区三区| 国产三级在线视频| 大型av网站在线播放| 国产亚洲欧美98| 搡老熟女国产l中国老女人| 女人高潮潮喷娇喘18禁视频| 免费看美女性在线毛片视频| 可以免费在线观看a视频的电影网站| 在线观看舔阴道视频| 久久欧美精品欧美久久欧美| 亚洲人成77777在线视频| 国产麻豆成人av免费视频| 亚洲精品国产色婷婷电影| 国产精品一区二区精品视频观看| 亚洲成人国产一区在线观看| 制服诱惑二区| 黑人操中国人逼视频| 亚洲av成人不卡在线观看播放网| 久久人妻福利社区极品人妻图片| 99久久综合精品五月天人人| 免费无遮挡裸体视频| 亚洲国产精品sss在线观看| 国产一区二区激情短视频| 精品欧美国产一区二区三| 久久人妻av系列| 99国产精品一区二区三区| 亚洲一区二区三区色噜噜| 99国产精品免费福利视频| 99精品久久久久人妻精品| 国产精品电影一区二区三区| 国产精品久久久久久亚洲av鲁大| 久久中文字幕人妻熟女| 国产精品久久电影中文字幕| 一级a爱视频在线免费观看| 男女做爰动态图高潮gif福利片 | 老司机在亚洲福利影院| 男女床上黄色一级片免费看| 禁无遮挡网站| 91在线观看av| 欧美在线黄色| 欧美乱妇无乱码| 免费观看精品视频网站| 午夜福利18| 99国产综合亚洲精品| 久9热在线精品视频| 国产1区2区3区精品| 天天躁夜夜躁狠狠躁躁| 叶爱在线成人免费视频播放| 十分钟在线观看高清视频www| 久久精品亚洲熟妇少妇任你| av视频免费观看在线观看| 丝袜人妻中文字幕| 黑人操中国人逼视频| 日本欧美视频一区| 香蕉国产在线看| 中文字幕高清在线视频| 色在线成人网| 麻豆av在线久日| 91老司机精品| 久久国产乱子伦精品免费另类| 亚洲av第一区精品v没综合| 日韩欧美在线二视频| 夜夜看夜夜爽夜夜摸| 天堂√8在线中文| 国产精品美女特级片免费视频播放器 | 精品人妻在线不人妻| 这个男人来自地球电影免费观看| 欧美一级a爱片免费观看看 | 久久久久久免费高清国产稀缺| 老汉色av国产亚洲站长工具| 日本vs欧美在线观看视频| www日本在线高清视频| 国产高清激情床上av| 1024视频免费在线观看| 一区二区三区国产精品乱码| 搡老岳熟女国产| 久久天躁狠狠躁夜夜2o2o| 人人澡人人妻人| 伊人久久大香线蕉亚洲五| 大陆偷拍与自拍| 狂野欧美激情性xxxx| 757午夜福利合集在线观看| 精品国产美女av久久久久小说| 狠狠狠狠99中文字幕| 国产激情久久老熟女| 国产精品99久久99久久久不卡| 日本精品一区二区三区蜜桃| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品合色在线| 一二三四社区在线视频社区8| 青草久久国产| 欧美乱色亚洲激情| 久久精品影院6| 日韩一卡2卡3卡4卡2021年| 搞女人的毛片| 亚洲精品粉嫩美女一区| 午夜成年电影在线免费观看| 女人被躁到高潮嗷嗷叫费观| 免费女性裸体啪啪无遮挡网站| 欧美日韩亚洲国产一区二区在线观看| 老熟妇仑乱视频hdxx| 国产精品av久久久久免费| 国产精品 欧美亚洲| 精品一区二区三区视频在线观看免费| 久久中文字幕人妻熟女| 国产精品一区二区在线不卡| 精品电影一区二区在线| 一个人观看的视频www高清免费观看 | 亚洲专区国产一区二区| 99在线人妻在线中文字幕| 国产精品美女特级片免费视频播放器 | 男女床上黄色一级片免费看| 久久精品国产99精品国产亚洲性色 | 久久久国产欧美日韩av| 夜夜爽天天搞| 欧美久久黑人一区二区| 国产真人三级小视频在线观看| a级毛片在线看网站| 久久青草综合色| 无遮挡黄片免费观看| 99riav亚洲国产免费| 热re99久久国产66热| 精品国产一区二区三区四区第35| 后天国语完整版免费观看| 精品久久久久久久人妻蜜臀av | 91国产中文字幕| 99精品在免费线老司机午夜| 在线视频色国产色| 亚洲一区中文字幕在线| 婷婷六月久久综合丁香| 女警被强在线播放| 久久国产精品影院| 国内毛片毛片毛片毛片毛片| 人人澡人人妻人| 午夜福利18| 大香蕉久久成人网| 国产免费男女视频| 亚洲精品国产精品久久久不卡| 精品久久蜜臀av无| 免费在线观看影片大全网站| 最近最新中文字幕大全免费视频| 久久人妻av系列| 亚洲精品一区av在线观看| 欧美在线一区亚洲| 亚洲av成人一区二区三| 中文字幕人妻丝袜一区二区| 亚洲欧美精品综合久久99| 精品一区二区三区视频在线观看免费| 淫秽高清视频在线观看| 国产精品一区二区免费欧美| 美女午夜性视频免费| 久久午夜亚洲精品久久| 老司机深夜福利视频在线观看| 美女高潮喷水抽搐中文字幕| 亚洲欧美精品综合久久99| 黄色 视频免费看| 91av网站免费观看| 黄色毛片三级朝国网站| 亚洲激情在线av| 禁无遮挡网站| 手机成人av网站| 亚洲av成人av| 欧美另类亚洲清纯唯美| 亚洲精品在线观看二区| 这个男人来自地球电影免费观看| 一本综合久久免费| 19禁男女啪啪无遮挡网站| 精品日产1卡2卡| 一进一出抽搐gif免费好疼| 国产成人精品久久二区二区91| 亚洲 欧美一区二区三区| 久久人人精品亚洲av| 给我免费播放毛片高清在线观看| 免费无遮挡裸体视频| 亚洲无线在线观看| 国内精品久久久久精免费| 亚洲国产看品久久| 桃色一区二区三区在线观看| 99久久久亚洲精品蜜臀av| 最近最新中文字幕大全免费视频| 真人做人爱边吃奶动态| 精品国产美女av久久久久小说| a在线观看视频网站| 91麻豆av在线| 一进一出抽搐动态| 18禁国产床啪视频网站| 男人操女人黄网站| 十分钟在线观看高清视频www| 看片在线看免费视频| 午夜精品国产一区二区电影| av在线天堂中文字幕| 国产精品久久电影中文字幕| 欧美精品亚洲一区二区| 精品一品国产午夜福利视频| 久久久久久久精品吃奶| 国产成人av教育| 午夜免费鲁丝| 欧美亚洲日本最大视频资源| 国产午夜精品久久久久久| 在线观看免费午夜福利视频| 精品久久久久久久人妻蜜臀av | 午夜福利,免费看| 国产熟女午夜一区二区三区| 中文字幕高清在线视频| 欧美日本亚洲视频在线播放| 日本免费a在线| 成人三级做爰电影| 老熟妇乱子伦视频在线观看| 国产av一区二区精品久久| 在线观看免费日韩欧美大片| 怎么达到女性高潮| 村上凉子中文字幕在线| 午夜福利视频1000在线观看 | 一个人观看的视频www高清免费观看 | 亚洲精品在线观看二区| 国产麻豆成人av免费视频| 精品一区二区三区视频在线观看免费| 一区二区三区精品91| 国语自产精品视频在线第100页| 国产av精品麻豆| 久久 成人 亚洲| 男女床上黄色一级片免费看| 国产在线精品亚洲第一网站| 亚洲男人天堂网一区| 欧美日韩亚洲国产一区二区在线观看| 精品国产一区二区三区四区第35| 色综合欧美亚洲国产小说| 乱人伦中国视频| tocl精华| 色哟哟哟哟哟哟| a级毛片在线看网站| 国产精品乱码一区二三区的特点 | 亚洲色图av天堂| 欧美日韩亚洲综合一区二区三区_| 真人做人爱边吃奶动态| 亚洲av日韩精品久久久久久密| 夜夜躁狠狠躁天天躁| 国产精品免费一区二区三区在线| 在线播放国产精品三级| 国产欧美日韩一区二区三| 欧美日韩亚洲国产一区二区在线观看| 国产又色又爽无遮挡免费看| 熟女少妇亚洲综合色aaa.| 午夜福利18| aaaaa片日本免费| 午夜福利高清视频| 好男人在线观看高清免费视频 | 国产精品一区二区三区四区久久 | 黑人操中国人逼视频| av有码第一页| 国产成人av激情在线播放| 精品乱码久久久久久99久播| 国产精品一区二区免费欧美| 亚洲第一青青草原| 国产精品 欧美亚洲| av在线播放免费不卡| 中出人妻视频一区二区| 亚洲第一av免费看| 熟女少妇亚洲综合色aaa.| 久久国产精品男人的天堂亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 九色亚洲精品在线播放| 亚洲熟妇熟女久久| 欧美丝袜亚洲另类 | 亚洲专区国产一区二区| 麻豆av在线久日| 曰老女人黄片| 变态另类丝袜制服| 一级毛片精品| 久久久久久久午夜电影| 高潮久久久久久久久久久不卡| 村上凉子中文字幕在线| 在线av久久热| 午夜两性在线视频| tocl精华| 中文字幕色久视频| 日韩欧美国产一区二区入口| 日韩有码中文字幕| 女同久久另类99精品国产91| 久久久久国内视频| 在线视频色国产色| 两人在一起打扑克的视频| 国产成人啪精品午夜网站| 亚洲国产精品999在线| 亚洲国产欧美一区二区综合| 成年版毛片免费区| 国产一卡二卡三卡精品| xxx96com| 巨乳人妻的诱惑在线观看| 亚洲av日韩精品久久久久久密| 亚洲一区高清亚洲精品| 欧美亚洲日本最大视频资源| 日韩视频一区二区在线观看| 久久久久精品国产欧美久久久| 精品少妇一区二区三区视频日本电影| 国产黄a三级三级三级人| 一区二区三区精品91| 午夜成年电影在线免费观看| 一本综合久久免费| av在线播放免费不卡| 国产成人精品久久二区二区91| 日本撒尿小便嘘嘘汇集6| 国产av一区在线观看免费| 国产欧美日韩一区二区三区在线| 黄片小视频在线播放| a在线观看视频网站| 欧美国产精品va在线观看不卡| 精品人妻在线不人妻| 国产极品粉嫩免费观看在线| 亚洲一区中文字幕在线| 国产精品影院久久| 88av欧美| 狂野欧美激情性xxxx| 亚洲五月色婷婷综合| 久久中文字幕一级| 欧美日韩一级在线毛片| 亚洲国产日韩欧美精品在线观看 | 亚洲av成人一区二区三| 国产一区二区三区视频了| 嫩草影院精品99| 无人区码免费观看不卡| 9191精品国产免费久久| 天天躁夜夜躁狠狠躁躁| 一进一出好大好爽视频| 99国产精品一区二区蜜桃av| 国产激情久久老熟女| 中文亚洲av片在线观看爽| 久久精品国产亚洲av香蕉五月| 久久亚洲真实| 欧美中文综合在线视频| 一级a爱视频在线免费观看| 亚洲第一青青草原| 亚洲第一欧美日韩一区二区三区| 亚洲自偷自拍图片 自拍| 国产一区在线观看成人免费| 久久亚洲真实| 国产日韩一区二区三区精品不卡| 亚洲精品中文字幕一二三四区| 我的亚洲天堂| 国产xxxxx性猛交| 国产熟女午夜一区二区三区| 999久久久国产精品视频| 一本久久中文字幕| 欧美日韩黄片免| av超薄肉色丝袜交足视频| 午夜两性在线视频| 亚洲国产精品成人综合色| 在线十欧美十亚洲十日本专区| 一级片免费观看大全| 亚洲专区字幕在线| 亚洲黑人精品在线| 成人永久免费在线观看视频| 亚洲精品中文字幕一二三四区| 国产成人av激情在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 欧美中文日本在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 人人妻人人澡人人看| 亚洲欧美日韩无卡精品| aaaaa片日本免费| 欧美激情久久久久久爽电影 | 国产人伦9x9x在线观看| 女人高潮潮喷娇喘18禁视频| 成熟少妇高潮喷水视频| 人人妻人人澡欧美一区二区 | 午夜激情av网站| 精品国产美女av久久久久小说| 黄色a级毛片大全视频| 色综合亚洲欧美另类图片| 日韩大尺度精品在线看网址 | 成在线人永久免费视频| 少妇熟女aⅴ在线视频| 91成年电影在线观看| 女人被狂操c到高潮| 久热这里只有精品99| 精品欧美一区二区三区在线| 宅男免费午夜| svipshipincom国产片| 午夜福利,免费看| 欧美中文综合在线视频| 亚洲国产毛片av蜜桃av| 国产精品永久免费网站| 色综合亚洲欧美另类图片| 国产99久久九九免费精品| 欧美黄色片欧美黄色片| 啦啦啦韩国在线观看视频| 国产私拍福利视频在线观看| 免费久久久久久久精品成人欧美视频| 天堂动漫精品| 精品一区二区三区av网在线观看| 亚洲性夜色夜夜综合| 亚洲色图综合在线观看| 国产精品精品国产色婷婷| 国产精品电影一区二区三区| 欧美激情 高清一区二区三区| 欧美最黄视频在线播放免费| 一二三四社区在线视频社区8| 91九色精品人成在线观看| 日韩三级视频一区二区三区| 国产精品乱码一区二三区的特点 | 美女扒开内裤让男人捅视频| 国产成年人精品一区二区| 一卡2卡三卡四卡精品乱码亚洲| 成人国产一区最新在线观看| 欧美成人免费av一区二区三区| www.自偷自拍.com| √禁漫天堂资源中文www| 99精品欧美一区二区三区四区| 亚洲av电影不卡..在线观看| 亚洲熟妇中文字幕五十中出| 女人爽到高潮嗷嗷叫在线视频| 亚洲黑人精品在线| 欧美不卡视频在线免费观看 | 他把我摸到了高潮在线观看| 美女扒开内裤让男人捅视频| 在线天堂中文资源库| 精品国产乱子伦一区二区三区| 欧美日韩亚洲综合一区二区三区_| 窝窝影院91人妻| 91av网站免费观看| 韩国av一区二区三区四区| 久久天躁狠狠躁夜夜2o2o| 好男人在线观看高清免费视频 | 欧美日韩中文字幕国产精品一区二区三区 | 午夜激情av网站| 在线观看午夜福利视频| 午夜激情av网站| 久99久视频精品免费| 一进一出抽搐gif免费好疼| 91麻豆av在线| 久久国产乱子伦精品免费另类| 国产国语露脸激情在线看| 国产精品久久久av美女十八| 99久久综合精品五月天人人| www.999成人在线观看| 两个人看的免费小视频| 免费观看人在逋| 欧美日本视频| 久久精品成人免费网站| 国产精品久久久久久亚洲av鲁大| 很黄的视频免费| 国产主播在线观看一区二区| 亚洲国产精品久久男人天堂| 国产精品日韩av在线免费观看 | 亚洲少妇的诱惑av| √禁漫天堂资源中文www| 久久精品aⅴ一区二区三区四区| 久久 成人 亚洲| 精品国产美女av久久久久小说| 波多野结衣巨乳人妻| 色综合站精品国产| 国产精品98久久久久久宅男小说| 色播亚洲综合网| 成人精品一区二区免费| 久久久精品欧美日韩精品| 亚洲 欧美一区二区三区| 搡老妇女老女人老熟妇| 午夜免费鲁丝| 亚洲欧美一区二区三区黑人| 欧美成人午夜精品| 国产麻豆69| 国产av在哪里看| 久久国产精品人妻蜜桃| 国产伦人伦偷精品视频| 国产1区2区3区精品| 成人亚洲精品一区在线观看| 又紧又爽又黄一区二区| 在线观看午夜福利视频| 亚洲美女黄片视频| 成人18禁高潮啪啪吃奶动态图| 亚洲男人天堂网一区| 免费在线观看视频国产中文字幕亚洲| 久久性视频一级片| 国产一级毛片七仙女欲春2 | 国产精品 欧美亚洲| 欧美日韩中文字幕国产精品一区二区三区 | 午夜福利影视在线免费观看| 午夜免费激情av| 黑人巨大精品欧美一区二区mp4| aaaaa片日本免费| 亚洲,欧美精品.| 精品久久久久久久久久免费视频| 制服丝袜大香蕉在线| 女人爽到高潮嗷嗷叫在线视频| 亚洲五月天丁香| 国产乱人伦免费视频| 国产麻豆69| 变态另类成人亚洲欧美熟女 | 精品卡一卡二卡四卡免费| 18美女黄网站色大片免费观看| 国产精品一区二区精品视频观看| 69精品国产乱码久久久| 久久人人97超碰香蕉20202| 午夜精品在线福利| 亚洲精品在线观看二区| 中文字幕人妻熟女乱码| 黄色成人免费大全| 国产色视频综合| 在线观看66精品国产| aaaaa片日本免费| 91字幕亚洲| 日韩一卡2卡3卡4卡2021年| 身体一侧抽搐| 法律面前人人平等表现在哪些方面| 日本 av在线| 国产97色在线日韩免费| 男人舔女人下体高潮全视频| 亚洲男人的天堂狠狠| 成人亚洲精品一区在线观看| 久久精品亚洲精品国产色婷小说| 一区在线观看完整版| 日本免费一区二区三区高清不卡 | 亚洲中文字幕一区二区三区有码在线看 | 色播亚洲综合网| 亚洲五月色婷婷综合| 狂野欧美激情性xxxx| 黄色女人牲交| 老司机深夜福利视频在线观看| 日本 欧美在线| 亚洲成人精品中文字幕电影| 欧洲精品卡2卡3卡4卡5卡区| 久久中文字幕人妻熟女| 老熟妇乱子伦视频在线观看| 97超级碰碰碰精品色视频在线观看| 大型黄色视频在线免费观看| 成人18禁高潮啪啪吃奶动态图| 久久精品国产亚洲av高清一级| 老司机在亚洲福利影院| 女人精品久久久久毛片| 波多野结衣av一区二区av| 国产激情欧美一区二区|