• <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.

    欧美成人一区二区免费高清观看| 精品久久久久久久久久免费视频| 赤兔流量卡办理| 国产精品亚洲一级av第二区| 99久久精品热视频| 久久久久国内视频| 精品99又大又爽又粗少妇毛片| 精品久久久久久久人妻蜜臀av| 日本-黄色视频高清免费观看| 婷婷六月久久综合丁香| 欧美xxxx黑人xx丫x性爽| 欧美成人精品欧美一级黄| 午夜视频国产福利| 一区二区三区高清视频在线| 亚洲熟妇熟女久久| 国产精品美女特级片免费视频播放器| 欧美潮喷喷水| 老熟妇仑乱视频hdxx| 亚洲av免费在线观看| 夜夜爽天天搞| 日本在线视频免费播放| 亚洲色图av天堂| 国产 一区精品| 国产一区二区三区在线臀色熟女| 在线免费十八禁| 老女人水多毛片| 一个人看的www免费观看视频| av国产免费在线观看| 97在线视频观看| 日韩高清综合在线| av在线老鸭窝| 日韩欧美 国产精品| 97超视频在线观看视频| 禁无遮挡网站| 亚洲三级黄色毛片| 国产精品综合久久久久久久免费| 长腿黑丝高跟| 亚洲最大成人手机在线| 免费观看精品视频网站| 小说图片视频综合网站| 国产亚洲91精品色在线| 给我免费播放毛片高清在线观看| 精品一区二区三区av网在线观看| 亚洲av五月六月丁香网| www.色视频.com| 日韩av在线大香蕉| 亚洲图色成人| 在线观看66精品国产| 午夜a级毛片| 校园人妻丝袜中文字幕| 国产麻豆成人av免费视频| 天堂网av新在线| 一级黄色大片毛片| 少妇被粗大猛烈的视频| 亚洲av成人精品一区久久| 色综合色国产| 搡老岳熟女国产| 日韩强制内射视频| 在线免费十八禁| 亚洲在线观看片| 国产毛片a区久久久久| av天堂在线播放| 国产黄a三级三级三级人| 国产精华一区二区三区| 一本精品99久久精品77| 精品福利观看| 一级a爱片免费观看的视频| 嫩草影院精品99| 在线免费观看不下载黄p国产| 国产在线精品亚洲第一网站| 国产精品亚洲一级av第二区| 色综合亚洲欧美另类图片| 亚洲精品影视一区二区三区av| 99热这里只有是精品在线观看| 赤兔流量卡办理| 日本五十路高清| 俄罗斯特黄特色一大片| 99热6这里只有精品| 亚洲人成网站在线观看播放| 天天躁夜夜躁狠狠久久av| 亚洲无线在线观看| 亚洲成人久久爱视频| 精品国内亚洲2022精品成人| 精品熟女少妇av免费看| 久久久久久伊人网av| 美女内射精品一级片tv| 亚洲美女搞黄在线观看 | 国产精品,欧美在线| 男女视频在线观看网站免费| 久久精品国产99精品国产亚洲性色| 欧美高清成人免费视频www| 欧美一区二区精品小视频在线| 国产精品人妻久久久影院| 亚洲va在线va天堂va国产| 精品久久久久久成人av| 免费无遮挡裸体视频| 亚洲精品日韩av片在线观看| 久久久精品欧美日韩精品| 秋霞在线观看毛片| a级一级毛片免费在线观看| 天天躁夜夜躁狠狠久久av| 亚洲精品粉嫩美女一区| 亚洲四区av| 久久人妻av系列| 99视频精品全部免费 在线| 我的女老师完整版在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美高清性xxxxhd video| 日韩精品中文字幕看吧| 婷婷色综合大香蕉| 日韩成人av中文字幕在线观看 | 国产伦精品一区二区三区视频9| 一夜夜www| 国产精品伦人一区二区| 国产精品99久久久久久久久| 无遮挡黄片免费观看| 超碰av人人做人人爽久久| 日本黄色片子视频| 亚洲国产日韩欧美精品在线观看| 久久久午夜欧美精品| 中国国产av一级| 女人被狂操c到高潮| 亚洲,欧美,日韩| 亚洲国产精品sss在线观看| 日韩成人av中文字幕在线观看 | 亚州av有码| 最近手机中文字幕大全| 日本一本二区三区精品| 午夜福利在线观看免费完整高清在 | 久久久午夜欧美精品| 中文字幕熟女人妻在线| 亚洲一区高清亚洲精品| 大香蕉久久网| 久久久久久久久大av| 免费观看人在逋| 在现免费观看毛片| 亚洲自拍偷在线| 有码 亚洲区| 国产成人a∨麻豆精品| 成年免费大片在线观看| 久久人妻av系列| 网址你懂的国产日韩在线| 亚洲自拍偷在线| a级毛片a级免费在线| 国产亚洲av嫩草精品影院| 国产又黄又爽又无遮挡在线| 国产精品亚洲一级av第二区| 97超视频在线观看视频| 久久久久国内视频| 国产麻豆成人av免费视频| 日本在线视频免费播放| 国产激情偷乱视频一区二区| 精品少妇黑人巨大在线播放 | av视频在线观看入口| 一本一本综合久久| 99热这里只有精品一区| 18+在线观看网站| 你懂的网址亚洲精品在线观看 | 小蜜桃在线观看免费完整版高清| 亚洲中文字幕日韩| 男人舔女人下体高潮全视频| 露出奶头的视频| 久久久色成人| 久久久a久久爽久久v久久| 国产精品1区2区在线观看.| 久久精品国产自在天天线| 欧美一区二区亚洲| 精品久久久久久成人av| 国内精品美女久久久久久| 成人特级av手机在线观看| 美女高潮的动态| 亚洲av五月六月丁香网| 亚洲国产欧洲综合997久久,| 在线免费观看不下载黄p国产| 亚洲欧美日韩东京热| 日日撸夜夜添| 午夜激情欧美在线| 国产亚洲精品久久久久久毛片| 在线免费观看的www视频| 人妻少妇偷人精品九色| 国产伦一二天堂av在线观看| 日韩制服骚丝袜av| 男人和女人高潮做爰伦理| 在线国产一区二区在线| 亚洲精品成人久久久久久| 一边摸一边抽搐一进一小说| 赤兔流量卡办理| 国产一区二区在线观看日韩| 午夜福利视频1000在线观看| 久久久欧美国产精品| 国产不卡一卡二| 在现免费观看毛片| 色5月婷婷丁香| 亚洲美女视频黄频| 日韩亚洲欧美综合| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人91sexporn| 国产人妻一区二区三区在| 一个人看视频在线观看www免费| 久久久久久伊人网av| 少妇熟女欧美另类| 欧美不卡视频在线免费观看| 五月玫瑰六月丁香| 国产亚洲欧美98| 国产老妇女一区| 欧美性猛交黑人性爽| 晚上一个人看的免费电影| 成人特级黄色片久久久久久久| 欧美性猛交黑人性爽| 免费黄网站久久成人精品| 12—13女人毛片做爰片一| 日本爱情动作片www.在线观看 | 日韩强制内射视频| 欧洲精品卡2卡3卡4卡5卡区| 1000部很黄的大片| 精品免费久久久久久久清纯| 老师上课跳d突然被开到最大视频| av天堂在线播放| 一边摸一边抽搐一进一小说| 国产精品人妻久久久久久| 一区福利在线观看| 国产黄色视频一区二区在线观看 | 91在线观看av| 国产伦精品一区二区三区视频9| 简卡轻食公司| 亚洲欧美日韩东京热| 亚洲熟妇中文字幕五十中出| 欧美极品一区二区三区四区| 我的老师免费观看完整版| 男女之事视频高清在线观看| 真实男女啪啪啪动态图| 亚洲成人精品中文字幕电影| 久久久久九九精品影院| 91在线观看av| 国产私拍福利视频在线观看| 色视频www国产| 九九在线视频观看精品| 黄色视频,在线免费观看| 99久久成人亚洲精品观看| 午夜福利在线在线| 最后的刺客免费高清国语| 中国美白少妇内射xxxbb| 亚洲一级一片aⅴ在线观看| 日本黄色片子视频| 久久精品91蜜桃| 精品免费久久久久久久清纯| 国产探花在线观看一区二区| 午夜a级毛片| 麻豆乱淫一区二区| 在线天堂最新版资源| 我要看日韩黄色一级片| 天堂网av新在线| 3wmmmm亚洲av在线观看| 熟女人妻精品中文字幕| 亚洲av五月六月丁香网| 91久久精品国产一区二区三区| 九九爱精品视频在线观看| 亚洲av中文av极速乱| 久久午夜福利片| 国产亚洲精品综合一区在线观看| 国产精品,欧美在线| 久久精品国产99精品国产亚洲性色| 日本黄大片高清| 丝袜美腿在线中文| 麻豆一二三区av精品| 香蕉av资源在线| 色视频www国产| 身体一侧抽搐| 国产成人a∨麻豆精品| 18+在线观看网站| 国产精品免费一区二区三区在线| 老司机福利观看| 91午夜精品亚洲一区二区三区| 久久久久九九精品影院| 伊人久久精品亚洲午夜| 国产一区二区激情短视频| 在线播放国产精品三级| 国产三级在线视频| 最近在线观看免费完整版| 三级经典国产精品| 美女xxoo啪啪120秒动态图| 成人美女网站在线观看视频| 一区福利在线观看| 99热这里只有是精品在线观看| 成年女人看的毛片在线观看| 99国产精品一区二区蜜桃av| 午夜免费男女啪啪视频观看 | 国产伦一二天堂av在线观看| 哪里可以看免费的av片| 99久久成人亚洲精品观看| 日韩欧美精品v在线| 国产精品一区www在线观看| 亚洲精品一卡2卡三卡4卡5卡| 欧美性猛交╳xxx乱大交人| 欧美潮喷喷水| 高清日韩中文字幕在线| 成年女人毛片免费观看观看9| 亚洲精品亚洲一区二区| 99久国产av精品国产电影| 看免费成人av毛片| 亚洲最大成人av| 天堂网av新在线| 香蕉av资源在线| 色av中文字幕| 人人妻,人人澡人人爽秒播| 男人和女人高潮做爰伦理| h日本视频在线播放| 日本爱情动作片www.在线观看 | 一个人观看的视频www高清免费观看| 中文字幕av在线有码专区| 亚州av有码| 在线播放国产精品三级| 五月伊人婷婷丁香| a级毛色黄片| 日韩欧美 国产精品| 俺也久久电影网| 欧美不卡视频在线免费观看| 偷拍熟女少妇极品色| 成人欧美大片| 国产一区亚洲一区在线观看| 色尼玛亚洲综合影院| 97人妻精品一区二区三区麻豆| 久久久久性生活片| 日韩欧美 国产精品| 亚洲四区av| 在线观看66精品国产| 亚洲av熟女| 久久久久久久久久成人| 国产一区二区亚洲精品在线观看| 日本一本二区三区精品| 国产成人91sexporn| 欧美极品一区二区三区四区| 五月伊人婷婷丁香| 最新在线观看一区二区三区| av在线观看视频网站免费| 国产精品一区www在线观看| 免费观看在线日韩| 最近手机中文字幕大全| av天堂中文字幕网| 日韩av在线大香蕉| 日韩欧美 国产精品| 国产亚洲精品久久久com| 国产精品亚洲一级av第二区| 天天躁日日操中文字幕| 婷婷亚洲欧美| 欧美日韩乱码在线| 18+在线观看网站| 欧美性猛交黑人性爽| 啦啦啦韩国在线观看视频| 日韩三级伦理在线观看| 蜜桃久久精品国产亚洲av| 国产真实伦视频高清在线观看| 国产精品一二三区在线看| 亚洲精品亚洲一区二区| 美女 人体艺术 gogo| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看| 夜夜看夜夜爽夜夜摸| 欧美成人a在线观看| 老熟妇乱子伦视频在线观看| 校园人妻丝袜中文字幕| 国产精品亚洲美女久久久| 黄色视频,在线免费观看| 亚洲人成网站在线观看播放| 91久久精品国产一区二区三区| 欧美潮喷喷水| 国国产精品蜜臀av免费| 成人欧美大片| 亚洲精品粉嫩美女一区| 国产精品野战在线观看| АⅤ资源中文在线天堂| 精品久久久久久久久久免费视频| 亚洲自偷自拍三级| 99热这里只有是精品50| 亚洲精品色激情综合| av黄色大香蕉| 久久久久久伊人网av| 中文在线观看免费www的网站| 久久精品国产亚洲av香蕉五月| 别揉我奶头 嗯啊视频| 看非洲黑人一级黄片| 女的被弄到高潮叫床怎么办| 国产不卡一卡二| 亚洲av电影不卡..在线观看| 少妇被粗大猛烈的视频| 亚洲av电影不卡..在线观看| 一级a爱片免费观看的视频| 99久国产av精品| 人人妻人人澡欧美一区二区| 美女高潮的动态| 午夜福利18| 91精品国产九色| 精华霜和精华液先用哪个| 国产v大片淫在线免费观看| 五月伊人婷婷丁香| 成人亚洲精品av一区二区| 搡女人真爽免费视频火全软件 | 无遮挡黄片免费观看| 国产69精品久久久久777片| 久久久午夜欧美精品| 夜夜夜夜夜久久久久| av天堂在线播放| 久久99热6这里只有精品| 人人妻人人澡欧美一区二区| 欧美丝袜亚洲另类| 亚洲精品在线观看二区| 久久精品国产亚洲网站| 日产精品乱码卡一卡2卡三| 麻豆av噜噜一区二区三区| 亚洲成人精品中文字幕电影| 午夜福利成人在线免费观看| 黄片wwwwww| 欧美极品一区二区三区四区| 热99在线观看视频| 国产久久久一区二区三区| 色综合亚洲欧美另类图片| 日韩欧美国产在线观看| 亚洲熟妇中文字幕五十中出| 可以在线观看毛片的网站| 成人特级av手机在线观看| 国产亚洲精品综合一区在线观看| 俄罗斯特黄特色一大片| 亚洲国产精品合色在线| 一个人观看的视频www高清免费观看| 欧美最新免费一区二区三区| 成人无遮挡网站| 欧美日韩一区二区视频在线观看视频在线 | 美女xxoo啪啪120秒动态图| 国产精品亚洲美女久久久| 亚洲av第一区精品v没综合| 日本熟妇午夜| 欧美zozozo另类| 三级经典国产精品| 最后的刺客免费高清国语| 成年女人看的毛片在线观看| 亚洲国产色片| 免费看a级黄色片| 久久久国产成人精品二区| 美女大奶头视频| 日本在线视频免费播放| 丰满人妻一区二区三区视频av| 十八禁国产超污无遮挡网站| 亚洲精品粉嫩美女一区| 中文字幕av在线有码专区| 天美传媒精品一区二区| 免费av不卡在线播放| 毛片女人毛片| 男女啪啪激烈高潮av片| 亚洲精品456在线播放app| 深夜a级毛片| 直男gayav资源| 99热只有精品国产| 国产精品野战在线观看| 欧美日韩乱码在线| 一级黄色大片毛片| 能在线免费观看的黄片| 日韩欧美精品v在线| 日本三级黄在线观看| 亚洲一级一片aⅴ在线观看| 国产精品野战在线观看| 男女之事视频高清在线观看| 欧美日韩综合久久久久久| 在线免费十八禁| 村上凉子中文字幕在线| 久久久久久久久久久丰满| 国产一区亚洲一区在线观看| 国产亚洲精品综合一区在线观看| 欧美激情在线99| 国产精品久久久久久久电影| 热99re8久久精品国产| 国产成人91sexporn| 国产精品一及| 成年av动漫网址| 欧美日本亚洲视频在线播放| 熟女人妻精品中文字幕| 此物有八面人人有两片| 成人性生交大片免费视频hd| 成人亚洲欧美一区二区av| 中国美女看黄片| 亚洲人成网站在线播放欧美日韩| 精品人妻熟女av久视频| 免费在线观看影片大全网站| 春色校园在线视频观看| 全区人妻精品视频| 久久午夜亚洲精品久久| 色哟哟哟哟哟哟| 精品人妻熟女av久视频| 亚洲精品亚洲一区二区| 国产69精品久久久久777片| 亚洲自拍偷在线| 亚洲人与动物交配视频| av.在线天堂| 午夜免费激情av| 欧美+亚洲+日韩+国产| 国产v大片淫在线免费观看| 香蕉av资源在线| 亚洲人成网站在线观看播放| 久久6这里有精品| 亚洲第一区二区三区不卡| 中文字幕精品亚洲无线码一区| 亚洲精品日韩av片在线观看| 亚洲av成人精品一区久久| 淫秽高清视频在线观看| 国内精品一区二区在线观看| 一本一本综合久久| 欧美日本亚洲视频在线播放| 国产又黄又爽又无遮挡在线| 露出奶头的视频| 99久久九九国产精品国产免费| 亚洲欧美日韩高清专用| 国产午夜精品久久久久久一区二区三区 | 精品一区二区三区av网在线观看| 听说在线观看完整版免费高清| 少妇高潮的动态图| 1024手机看黄色片| 深夜a级毛片| 久久久久久久亚洲中文字幕| 偷拍熟女少妇极品色| 三级男女做爰猛烈吃奶摸视频| 国产黄色视频一区二区在线观看 | 人人妻,人人澡人人爽秒播| 中文字幕免费在线视频6| 深夜精品福利| 不卡视频在线观看欧美| 给我免费播放毛片高清在线观看| 国产真实伦视频高清在线观看| 天堂影院成人在线观看| 我要看日韩黄色一级片| 三级男女做爰猛烈吃奶摸视频| 可以在线观看的亚洲视频| www日本黄色视频网| 国产私拍福利视频在线观看| 亚洲成av人片在线播放无| 天天躁日日操中文字幕| 欧美一区二区精品小视频在线| 成熟少妇高潮喷水视频| 国产av不卡久久| 一个人看视频在线观看www免费| 禁无遮挡网站| 日本与韩国留学比较| 天堂影院成人在线观看| 一级毛片久久久久久久久女| 99视频精品全部免费 在线| 一级黄片播放器| 久久久久久伊人网av| 在线国产一区二区在线| 欧美成人一区二区免费高清观看| 日本五十路高清| av黄色大香蕉| 我要看日韩黄色一级片| 夜夜看夜夜爽夜夜摸| 久久99热6这里只有精品| 亚洲一级一片aⅴ在线观看| 中国美白少妇内射xxxbb| 亚洲精品国产av成人精品 | 欧美高清性xxxxhd video| 亚洲内射少妇av| 国产黄片美女视频| 热99re8久久精品国产| 国产成人91sexporn| 18+在线观看网站| 国产成人影院久久av| 亚洲av中文av极速乱| 乱码一卡2卡4卡精品| 免费观看人在逋| 日本黄色视频三级网站网址| 久久草成人影院| 国产欧美日韩精品亚洲av| 伦理电影大哥的女人| av女优亚洲男人天堂| 亚洲美女视频黄频| 在线国产一区二区在线| 日本一二三区视频观看| 亚洲成人精品中文字幕电影| 国产探花在线观看一区二区| 亚洲aⅴ乱码一区二区在线播放| 中文字幕av成人在线电影| 久久精品国产亚洲av天美| 最近视频中文字幕2019在线8| 亚洲成人av在线免费| 麻豆精品久久久久久蜜桃| 高清日韩中文字幕在线| 成人永久免费在线观看视频| 一本一本综合久久| 国产不卡一卡二| 久久久欧美国产精品| 18禁裸乳无遮挡免费网站照片| 国产精品福利在线免费观看| 男女那种视频在线观看| 女人十人毛片免费观看3o分钟| 欧美一区二区精品小视频在线| 久久精品国产鲁丝片午夜精品| 国产精品精品国产色婷婷| 久久精品91蜜桃| 美女高潮的动态| 99久久精品热视频| 欧美潮喷喷水| 麻豆国产97在线/欧美| 久久99热这里只有精品18| 亚洲丝袜综合中文字幕| 欧美性感艳星| 久久久久国产精品人妻aⅴ院| 国产精品久久久久久av不卡| 国产成人精品久久久久久| 一进一出抽搐动态| 亚洲精品日韩av片在线观看| 久久精品国产亚洲av天美| 亚洲av美国av| 精品国内亚洲2022精品成人| 色综合站精品国产| 伦精品一区二区三区| 嫩草影院入口| 国产高潮美女av| 国内久久婷婷六月综合欲色啪| 久久精品影院6|