Accelerated oxygen evolution kinetics on Ir-doped SrTiO3 perovskite by NH3 plasma treatment

2022-11-21 09:30LiLiDeng邓丽丽XiaoPingMa马晓萍ManTingLu卢曼婷YiHe何弈RongLeiFan范荣磊andYuXin辛煜
Chinese Physics B 2022年11期
关键词:丽丽

Li-Li Deng(邓丽丽), Xiao-Ping Ma(马晓萍), Man-Ting Lu(卢曼婷),Yi He(何弈), Rong-Lei Fan(范荣磊), and Yu Xin(辛煜)

School of Physical Science and Technology,Jiangsu Key Laboratory of Thin Films,Soochow University,Suzhou 215006,China

Exploring low-cost and high-performance catalysts for oxygen evolution reaction(OER)remains to be a great challenge. Iridium-based perovskite oxide has large potential in OER because of its intrinsic activity and outstanding physicochemical properties. In this study, iridium-doped strontium titanate (Ir-STO) solution is brushed on a Ti sheet by the traditional method to obtain the Ir-STO/Ti electrodes after being calcined at a high temperature. The microstructure and electrocatalysis properties of the Ir-STO are further modified by a facile and scalable NH3-plasma strategy. In addition to the doping of Ir, the NH3 plasma treatment further results in N-doping into Ir-STO, which enriches active species and causes oxygen vacancies near doped sites. The resulting N,Ir-STO/Ti electrode reveals excellent acidic OER activity with the lowest overpotential of 390 mV at 10 mA/cm2 and the smallest Tafel slope of 140 mV/dec after 10-min plasma treatment. Therefore,the great potential of activated N,Ir-STO/Ti is regarded as a catalyst for the OER,and thus making a new opportunity for developing other perovskite catalysts via NH3 plasma treatment.

Keywords: electrocatalysts,NH3 plasma,oxygen vacancies,N doping

1. Introduction

Electrochemical splitting of water to produce hydrogen driven by renewable feedstocks (generated from wind or solar,etc.) is considered to be a promising approach to mass production of hydrogen in the future.[1–3]However, the anodic semi-reaction-oxygen evolution reaction (OER) process of water splitting requires an inevitable dynamic overpotential,resulting in high energy consumption.[4,5]Noble metal oxides such as IrO2and RuO2are known as the most effective OER catalysts,but their low reserves and high cost greatly limit the large-scale commercial applications.[6–8]Therefore, it is urgent to develop highly active electrocatalysts with low noble metal content.

A great deal of research has been done on transition metal-based materials.[9]Among them, perovskite oxides have attracted considerable attention because of their low price, good stability and adjustable physicochemical properties.[10–12]The general formula for perovskite is ABO3,

where transition metal B plays a key role in the catalytic process. Seitzet al.[13]demonstrated the role of B site in perovskite by the epitaxial growth of SrIrO3films on the surface of SrTiO3catalyst. Metal-ions doping is recognized as an efficient strategy to improve intrinsic electronic conductivity and modulate electrocatalytic activity.[14]Actually,on the premise of ensuring the structural stability of perovskite, introducing metal ions with different charges or sizes into the perovskite crystal lattice can make the lattice strain,thus adjusting the oxidation state.[15]Lianget al.[16]synthesized solid solution perovskite Ir-STO by polymerization complex method.Ir doping activated the inert STO matrix to generate Ti sites with catalytic activity. They found the synergistic effect of Ir and Ti to enhance the oxygen adsorption on the surface in the case of reduced iridium consumption.

Notably, when metal cations are doped into perovskite,oxygen vacancies will be induced in perovskite crystals due to the electroneutrality principle.[17]Meffordet al.[18]and She et al.[19]have confirmed the change of the oxygen vacancy concentration in La1-xSrxCoO3-δand La1-xSrxFeO3-δperovskites with substitution of La3+by Sr2+,respectively.Other studies showed that oxygen vacancies in metal oxides can improve the catalytic activity of OER oxides by reducing the adsorption energy of water.[20,21]Plasma treatment can effectively carve oxides,selectively remove oxygen on the surface,and produce oxygen vacancies.[22]Under the bombardment of high-energy ions (Ar+, N2+, H2O+), oxygen vacancies can be created on metal oxide surfaces. Luet al.[23]replaced the La site with Sr in LaCoO3, and then treated the sample with Ar plasma, resulting in the increase of the oxygen vacancies and a large number of active sites with high OER activity and stability.

Recent studies have shown that the performance of perovskite-based catalysts can be improved by introducing non-metallic elements (such as P, N, S,etc.) into the perovskite lattice to change the surface properties and electron/phase structure of perovskite.[24–26]Owing to different properties of nitrogen and oxygen(such as polarizability,electronegativity and charge),the introduction of N into perovskite oxide can change the physical and chemical properties.[26]Specifically,since N is less electronegative than O,N-doping can reduce the bandgap, thus improving electron transport in perovskite oxides.[27,28]The N partially replaces the O in perovskite oxides,forming perovskite with abundant oxygen vacancies and high oxidation states of B-site metal cations.[29]Zhanget al.[30]doped LaNiO3with N,which significantly increases the Ni3+content and the surface oxygen vacancies and promotes the anodic oxygen evolution reaction. In our experiment,we have found that the concentration of oxygen vacancies can be optimized by N,Ir doping into the STO matrix.

In this paper, in order to improve OER electrocatalytic performance, two strategies are employed for SrTiO3perovskites to generate oxygen vacancies. First, we partially replace Ti with Ir to change the electrocatalytic performance of SrTiO3. Then, N doping into the Ir-STO coatings by using NH3plasma treatment to produce oxygen vacancies. The Ir-STO coatings will be synthesized by high-temperature thermal decomposition on titanium sheets. The atomic ratio of Ir to Ti is 1:2, significantly reducing the amount of iridium used compared with the industrial proportion of 70%such as in the IrO2/RuO2system.[31]And then, capacitively coupled NH3plasma is used to treat the above samples and thus to adjust oxygen vacancies in the bulk. Herein, a 13.56-MHzdriven capacitively coupled NH3plasma with input power of 300 W is chosen with the different treatment times for the Ir-STO/Ti samples. The samples are characterized by using field emission-scanning electron microscopy (FE-SEM),x-ray diffractometer(XRD),x-ray photoelectron spectroscopy(XPS),and electrochemical workstation.

2. Experiment

2.1. Synthesis of Ir-STO/Ti

A Ti sheet (0.5-mm thick) with an exposed surface area of 1 cm×1 cm was used. It was sandblasted, ultrasonically cleaned for 10 min in a 15-vol% Na2CO3solution, and then etched for 30 min in oxalic acid (15 vol%). Sr(NO3)2,C16H36O4Ti, and H2IrCl6·6H2O were dissolved in a mixture of glycol and citric acid(the ratio of glycol to citric acid was 9:1)to form the precursor solution,sequentially. The Ti sheet was brushed with the precursor solution and then annealed in a muffle furnace at 300°C for 10 min. The above process was repeated three times,and finally,the samples were annealed at 600°C for 4 h to fabricate Ir-STO/Ti electrodes.

The fabricated Ir-STO/Ti electrodes were then treated by capacitively coupled NH3plasma to dope N atoms in the Ir-STO/Ti. The annealed Ir-STO/Ti was transferred into the capacitively coupled plasma(CCP)chamber with the discharge power of 300 W and the NH3flow rate of 80 sccm, and the samples were placed about 1 cm away from the lower plate to conduct double-sided treatment.The working pressure was set to be 5 Pa and the bias voltage to be 452 V nominally. In order to avoid large damage to Ir-STO/Ti due to plasma-induced degradation (PID), the treatment time of Ir-STO/Ti was controlled within 20 min. For controlling the nitridation degree,the treatment was performed for different durations of 5 min,10 min, 15 min, and 20 min, and the corresponding samples were denoted as N,Ir-STO/Ti-5 min,N,Ir-STO/Ti-10 min,N,Ir-STO/Ti-15 min,and N,Ir-STO/Ti-20 min,respectively.

2.2. Characterizations of structure

The crystal structure of the prepared samples were analyzed by using an XRD with bruker D8 advance diffractometer at a scanning range of 20°–80°.The surface morphology of the prepared Ir-STO/Ti and N,Ir-STO/Ti electrodes were characterized via FESEM with a Regulus Su8100 microscope. The chemical composition information of the samples was tested by XPS with a Thermo ESCALAB250XI spectrometer. The information about oxygen vacancies in Ir-STO/Ti and N, Ir-STO/Ti coatings were measured by electron spin resonance spectrometer(EPR,Model JES-X320,Japan).

2.3. Electrochemical measurements

Electrochemical tests were performed on a standard three-electrode system of electrochemical workstation. A three-electrode system consists of a working electrode with electrocatalyst of Ir-STO/Ti and N,Ir-STO/Ti,a platinum plate as the counter electrode and a saturated calomel electrode(SCE) as a reference electrode. From the following formula calculated are all potentials of the reversible hydrogen electrode(RHE)in this paper:

The linear scanning voltammetry (LSV) polarization curves were measured in a 0.5-M H2SO4solution from 1.0 V to 1.8 V(versusRHE)in potential steps of 5 mV/s. The cyclic voltammogram (CV) curves were made between 0.1 V and 1.5 V(versusRHE) at a scan rate of 50 mV/s through experiments in 0.5-M H2SO4solution. Electrochemical impedance spectroscopy (EIS) was used at a voltage of 1.7 V (versusRHE)and AC sinusoidal signal amplitude of 5.0 mV, and in a frequency range of 100 kHz–0.1 Hz. ZSimpDemo software is used to fit EIS data in the corresponding equivalent circuit.

3. Results and discussion

Interaction of plasma with surface of materials plays a crucial role in the modification of semiconductor fabrication.[32]However,the heavy elements for the Ir-STO/Ti may not be easily etched by NH3plasma treatment even under the ion bombardment due to ineffective momentum transfer in between. The SEM results(not shown here)demonstrate that the morphologies of the samples do not change significantly,indicating little effect of the plasma treatment on the morphology of perovskites.[23,33]Under the action of ion bombardment, the elements with small mass O in the Ir-STO coating surface can be easily etched under the bombardment of high energy ions of N+2,NH+,etc. Some of the neutral groups such as NHx,N in the plasma will diffuse onto the surface and fill in these vacancies with great probability, although they are also bombarded by ions. Therefore,the displacement reaction will occur in between the active species such as NHx, N,etc. in the NH3plasma[34]and lattice O atom in the Ir-STO coatings.This is an important reason for N atoms to dope into the perovskite matrix under the plasma action,confirmed by the XRD and XPS characterization.

Fig.1. XRD patterns of(a)Ir-STO/Ti and STO/Ti,with inset showing magnified part of patterns between 31.8°–32.8°;(b)XRD patterns of Ir-STO/Ti and N,Ir-STO/Ti treated by NH3 plasma for 5 min,10 min,15 min,20 min,respectively; (c)magnified XRD patterns with 2θ between 31.8°–32.8°;(d)error bar of typical SrTiO3 (110)diffraction peak.

The XRD results, depicted in Fig. 1(a), confirm that the STO/Ti sample and Ir-STO/Ti sample each present a perovskite-type phase. In comparison, no evident new peaks are observed in the XRD patterns but the diffraction peak(SrTiO3) gradually shifted toward small diffraction angles due to the doping Ir into STO/Ti as indicated in the inset of Fig. 1(a). Considering that the measurement error of this instrument is 0.01°, we can confirm that the slight shift(0.09°±0.01°, 2θ) originates from a high ionic radius iridium dopant into STO/Ti(Ir4+=62.5 pm,Ti4+=60.5 pm).[16]Figure 1(b) shows the crystal structure of samples after NH3plasma treatment, the bulk crystal structure does not change much in comparison with Ir-STO/Ti matrix. From the detailed observation of XRD patterns, it is found that the typical SrTiO3(110) diffraction peak of N, Ir-STO/Ti shown in Figs. 1(c)–1(d), shifts toward high diffraction angles in comparison with pristine Ir-STO/Ti. The shifting of 2θvalue can be attributed to two main factors. On the one hand, at the beginning of NH3plasma treatment,numerous energetic particles rapidly attack and react with surface lattice oxygen to form oxygen vacancies.[35]The charge neutrality after the generation of oxygen vacancies is compensated for by the reduction of Ti4+(0.605 ˚A) into Ti3+(0.670 ˚A).[36]On the other hand, the N atoms can easily occupy the oxygen vacancies,thus replacing the O atom positions and dope into substitutional sites, resulting in the modification of the electronic states.[37]As mentioned above,it seems that N and Ir doping may adjust the lattice parameter of STO/Ti.

The most direct evidence to determine whether N doping exists in Ir-STO/Ti sample comes from XPS analysis. The XPS measurement spectra indicate that besides the presence of Ir, Ti, Sr, and O, the information from N of Ir-STO/Ti samples treated by the NH3plasma also appears as shown in Fig. 2(a). It can be seen from Fig. 2(b) that N 1s chemical state changes greatly with the increase of NH3plasma treatment time. The N 1s XPS spectrum as shown in the figure can be well fitted to three characteristic peaks: oxynitride (Ir/Ti–O–N) species (399.8 eV), Ir–N bond (397.7 eV), and Ti–N bond (396.2 eV), respectively.[38–40]Various peaks of N 1s binding energy between 395.8 eV and 397.8 eV can be found elsewhere,[41]which are generally considered as being from other nitrogen source. These findings suggest that N may be introduced into Ir-STO/Ti lattice as the alternative dopant by NH3plasma treatment. Compared with Ir-STO/Ti, the peaks of Ir–N and Ti–N for the N,Ir-STO/Ti samples increase gradually with the increase of plasma treatment time,indicating an effective N-doping in the Ir-STO/Ti matrix.

Fig.2. XPS spectra of Ir-STO/Ti and N,Ir-STO/Ti with different NH3 plasma treatment times: (a)survey spectra,(b)high resolution spectra of N 1s,(c)Ir 4f,and(d)O 1s.

Table 1. XPS data fitting results of different Ir and O species of Ir-STO/Ti and N,Ir-STO/Ti with different NH3 plasma treatment times.

Owing to charge balancing in the N, Ir-STO/Ti, the incorporation of Ir and N atoms into the STO matrix may cause the valence state of the iridiumon the B-site to change. Figure 2(c)shows that the Ir 4f spectra of the samples reveal two doublets corresponding to the spin–orbit components Ir 4f7/2and Ir 4f5/2, respectively. The doublets are divided into Ir3+4f7/2(61.40 eV),Ir4+4f7/2(63.01 eV),Ir3+4f5/2(64.30 eV),and Ir4+4f5/2(66.01 eV).[42,43]Therefore,by integrating the peak area,the relative values of content of these Ir species are calculated as shown in Table 1.The Ir content decreases possibly due to its relatively high sputtering yield. The Ir4+content of the samples treated by NH3plasma increases significantly.This may be due to the charge compensation caused by N and Ir co-doping.[29]Interestingly, the Ir4+state may activate the inert STO matrix,generate efficient Ti catalytic sites,thus improving the performance of the catalyst.[16]

The incorporation of Ir and N into the STO matrix causes a slight change both in the lattice parameters as mentioned before and in the chemical state of O changes. The O 1s spectra are collected from the Ir-STO/Ti and N, Ir-STO/Ti samples as shown in Fig. 2(d). Typical peaks at 529.20 eV, 529.9 eV,and 531.2 eV correspond to the lattice oxygen (O2-), highly oxidative oxygen species (O2-/O-), and hydroxyl groups or surface-adsorbed oxygen species (–OH/O2), including –OH,O-,and O2-2 ,respectively.[40,44]In Table 1 listed are the relative values of content of each oxygen-related species for these N, Ir-STO/Ti coatings, which are obtained from the area ratio of corresponding speak. It should be pointed out that the relative value of content of O2-/O-species is closely related to surface oxygen vacancies.[45]It can be seen that the N, Ir-STO/Ti coating treated with plasma for 10 min shows the highest content of O2-/O-(0.42)in all the samples,which means that the quantity of oxygen vacancies is the highest. Note that the increasing of NH3plasma treatment time leads the O2-/Ocontent in N,Ir-STO/Ti coatings to decrease slightly.It may be due to the fact that excessive N atoms occupy oxygen vacancies,thereby replacing O sites and doping into the substitution sites.[46]Additionally,the increased–OH/O2content of the N,Ir-STO/Ti samples compared with that of the Ir-STO/Ti indicates the enhanced hydrophilicity of the catalyst, which can increase the adsorption of oxygen species on the catalyst surface,thus improving the catalytic activities.[47]

The changes of oxygen vacancies for the N, Ir-STO/Ti coatings are further probed by the low-temperature electron paramagnetic resonance(EPR)spectrum. As shown in Fig.3,Ir-STO/Ti and N,Ir-STO/Ti coatings each present a response at ag-value of 2.003, indicating the existence of isolated oxygen vacancies. It can also be seen from Fig. 3 that the peak value of EPR increases gradually when the plasma treatment time is less than 10 min and decreases gradually when the plasma treatment time is more than 10 min. For the plasma treatment time of 10 min, the peak value of EPR is the strongest, indicating the highest oxygen vacancy content in N, Ir-STO/Ti coating. The results from EPR are consistent with the XPS analysis ones obtained before. The ions from the bulk plasma will bombard the sample surface, oxygen bonded to either Ti or Ir on the surface of matrix will be sputtered or etched by ions or the other active species such as NHxand N.The presence of N dopant facilitates the formation of oxygen vacancies. As a result, suitable plasma treatment of 10 min can improve the vacancy density. With excessive plasma treatment, more oxygen vacancies including oxygen vacancies caused by Ir doping will be replaced by N atoms.That is to say,the incorporation of both N and Ir into the STO matrix may modulate the oxygen vacancies in the bulk.[48]

Fig. 3. EPR spectra of Ir-STO/Ti and N, Ir-STO/Ti with different NH3 plasma treatment times.

The LSV tests are conducted to measure some parameters of Ir-STO/Ti and N,Ir-STO/Ti electrodes at 5-mV/s scanning rate in 0.5-M H2SO4solution. As shown in Fig.4(a),the OER catalytic activity of Ir-STO/Ti is remarkable with an overpotential of 410 mV at a current density of 10 mA/cm2, while the OER activity of STO/Ti is negligible. This result indicates that Ir doping enhances the OER activity of STO/Ti,which is consistent with the result from other research.[16]Figure 4(b)shows the LSV curves of Ir-STO/Ti treated by NH3plasma for different treatment times. It is found that the N,Ir-STO/Ti-10-min catalyst exhibits the best activity with the smallest overpotential of 390 mV at a current density of 10 mA/cm2. Basically,the N,Ir-STO/Ti coatings show better activity than Ir-STO/Ti except the one with 20 min. Owing to the saturation of oxygen vacancy and the decrease of Ir content, the electrocatalytic activity of N,Ir-STO/Ti-20-min electrode is lower than that of Ir-STO/Ti. Among them,the lowest Tafel slope of N,Ir-STO/Ti-10 min(140 mV/dec)is also found as shown in Fig.4(c).

Figure 4(d) shows the voltammetry curves of various samples recorded at a scanning rate of 50 mV/s in a 0.5-M H2SO4solution. The pair of peaks at 1 V–1.2 V(versusRHE)are attributed to the surface redox transition of Ir(III)/Ir(IV),which is a pair of typical redox peaks in the electrochemical process of Ir-based electrocatalysts.[49,50]Within this 0.4 V–1.44 V potential range, hydrogen evolution reaction or OER will not occur on the electrode surface. Under this condition, the current belongs to the electric double layer charging current.[51,52]Integrating the CV curve in this potential range can obtain the voltammetric charge capacity(q*),which represents the number of protons exchanged by the electrode with the solution.[53]Theq*is proportional to the electrochemically active surface area of the oxide electrode (EASA),[54]and can reflect the number of electrochemically active centers on the surface of the electrode. Figure 4(e)shows the highest EASA of about 200 mC/cm2for the N,Ir-STO/Ti-10 min sample. Meanwhile,the electrochemical double-layer capacitance(Cdl) for the sample in Fig. 4(f) also shows the highest value of 6.18 mF/cm2.

Figure 5 shows the impedance spectra of the Ir-STO/Ti and N,Ir-STO/Ti samples in a 0.5-M H2SO4solution. The resistance of the H2SO4solution is characterized byRs, which is connected in series with the coating system unit.[55]TheRctis the resistance of charge transfer on the electrode solution interface, and its value reflects the electrocatalytic activity of OER.[56]TheQdlvalue represents the number of active points on the inner surface of the coating.[57,58]The fitted EIS data for all the electrical parameters are listed in Table 2. TheQdlvalue of the sample treated by NH3plasma for 10 min is significantly higher than those of other samples, indicating that the NH3plasma treatment for 10 min can maximize the number of active points on the coating surface. Moreover,theRct(11.10 Ω/cm2) value of NH3plasma-treated samples for 10 min is the smallest,implying a better conductivity.This test result is consistent with the polarization measurement mentioned above.

Fig.4. (a)Polarization curves of pristine STO/Ti and the Ir-STO/Ti. (b)Polarization curves of Ir-STO/Ti and N,Ir-STO/Ti,with inset showing chronopotentiometry test at 10 mA/cm2. (c)Corresponding Tafel plots. (d)CV curves of Ir-STO/Ti and N,Ir-STO/Ti. (e)Voltammetric charges obtained by integrating the curves. (f)Value of Cdl for the Ir-STO/Ti and N,Ir-STO/Ti with different NH3 plasma treatment times.

Fig.5. Nyquist plots of N,Ir-STO/Ti prepared by different times,with inset showing equivalent circuit model.

Table 2. Fitting EIS parameters of Ir-STO/Ti and N,Ir-STO/Ti samples.

For 13.56-MHz-driven capacitively coupled NH3plasma,active species such as NH2,H,and N can be observed by optical emission spectra.[59]As these species transport into the material surface, interaction in between will occur. Among them,extreme reducible H atom causes the metal oxides to reduce to metal by hydrogen reacting with oxygen atoms on the surface. For STO, the large bond energies of Sr–O and Ti–O(426.3 kJ/mol,666.5 kJ/mol)[60]become a little bit difficult to reduce H atoms. It is noteworthy that the treated samples are placed on the powered electrode surface, the self-bias on the sheath approaches to-360 V as 300-W radio frequency power is applied. For the treatment condition, discharge pressure is set to be 5.0 Pa,that is,ion acceleration during transportation through sheath will endure collision with neutral particles less than 4 times averagely due to ion mean free path of 1.2 mm and the calculated sheath width of about 3.7 mm.Ion energy to the STO surface on the powered electrode is enough to cause Ti–O or Sr–O bonds to break. Thus, in the combination of ion bombardment assistance, the reduction of H atom on the STO surface becomes relatively easy. As O atom is devoured by the H atom,the left sites will be occupied by N-containing species with a significant probability. A shift of STO/Ti lattice (110) diffraction peak to a lower angle in XRD and the appearance of N element and its change of chemical state in XPS have confirmed the binding of N into Ir-STO/Ti matrix.Additionally,from the diffraction angle in the opposite direction due to N and Ir doping in Fig.1,it can be inferred that the two kinds of elements can modulate the lattice parameters of N,Ir-STO/Ti. From Fig.4,N and Ir doping may optimize the bulk oxygen vacancies with a maximal vacancy concentration at 10-min NH3plasma treatment. Therefore, the best electrocatalytic performance of N,Ir-STO/Ti coatings is observed with NH3plasma treatment time increasing up to 10 min.

4. Conclusions

In summary, we have successfully created more O vacancies in SrTiO3perovskite by Ir doping and NH3plasma treatment. By combining two techniques,we achieve the outstanding OER performance of N,Ir-STO/Ti compared with the pristine STO/Ti. In addition,we further explore the influence of different plasma treatment time on the electrocatalytic performance of N, Ir-STO/Ti. It is found that suitable plasma treatment will improve the cacancy density of N, Ir-STO/Ti,which can improve the charge transport properties, resulting in superior OER performance. The N,Ir-STO/Ti electrode reveals excellent acidic OER activity with the lowest overpotential of 390 mV at 10 mA/cm2and the smallest Tafel slope of 140 mV/dec after 10-min NH3plasma treatment. With excessive plasma treatment, the redundant N atoms can easily occupy the oxygen vacancies,leading the vacancy concentration and mobility to decrease. The results reveal that a moderate plasma treatment can balance the generation of vacancies and N doping to obtain the best OER characteristics. Our studies make a new opportunity for designing electrocatalysts with significantly low noble metal content and providing a generic method to realize controllable doping of other perovskite catalysts via plasma treatment.

Acknowledgements

Project supported by the Priority Academic Program Development (PAPD) Program of Jiangsu Higher Education Institutions, Jiangsu Province, China and the National Natural Science Foundation of China(Grant No.11675117).

猜你喜欢
丽丽
ASYMPTOTIC STABILITY OF A BOUNDARY LAYER AND RAREFACTION WAVE FOR THE OUTFLOW PROBLEM OF THE HEAT-CONDUCTIVE IEAL GAS WH SY*
Westward Movement
画一画
《咔嚓!老田就爱高丽丽》
漂亮的生日礼物
丽丽学画画
Green product development
On the ecological concept in design
I love my family
丽丽的一天