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

    Direct measurement of an energy-dependent single-event-upset cross-section with time-of-flight method at CSNS

    2023-03-13 09:18:20BiaoPei裴標(biāo)ZhixinTan譚志新YongningHe賀永寧XiaolongZhao趙小龍andRuiruiFan樊瑞睿
    Chinese Physics B 2023年2期
    關(guān)鍵詞:永寧小龍

    Biao Pei(裴標(biāo)) Zhixin Tan(譚志新) Yongning He(賀永寧)Xiaolong Zhao(趙小龍) and Ruirui Fan(樊瑞睿)

    1Spallation Neutron Source Science Center,Dongguan 523803,China

    2Instituteof High Energy Physics,Chinese Academy of Sciences(CAS),Beijing 100049,China

    3School of Microelectronics,Xi’an Jiaotong University,Xi’an 710049,China

    4State Key Laboratory of Particle Detection and Electronics,Bejing 100049,China

    Keywords: static random-access memory,soft error rate,neutron SEU cross-section,time-of-flight

    1.Introduction

    Single event effects (SEEs) are running errors in semiconductor devices caused by a single,energetic particle which originates from cosmic rays and/or their showers.SEEs take place anytime and anywhere,from the galaxy space or the near space where satellites fly-by to the territorial ground where mankind lives.[1]In order to emulate their impacts to semiconductor devices, many particle facilities (for protons, neutrons and also ions)were built to provide test beams with a resemble spectrum.Special efforts were put up for neutrons as they are hard to generate and control.In particular, atmospheric neutron beam lines from spallation neutron sources are in the center of attention,such as FP30L at LANSCE,since neutrons are the dominant irradiation source in the territorial environment.However,atmospheric neutron beams only address a group of common applications in the SEE tests.As mankind explores the space,[2,3]more demands with diverse neutron spectra will be encountered.In these cases, previous method with a test beam of a specific spectrum is no longer applicable.The requirements for a diversity of applications demand a new way with a profound comprehension in SEEs,in particular,a study on their cross-section with a variety of neutron energies.

    There are kinds of effects caused by neutrons.The most common one is the single event upset(SEU),which is a state change in logic or bit-flips in memory.SEUs are ‘soft’ and non-destructive.It has been widely investigated.Although SEU research with single energy points had been collected for energy-dependence evaluation years ago, the summary is limited by the data points.Furthermore,the results in Johansson’s work[4]and a later report from Dyer[5]show a contradiction in the general trends range from 1 MeV to 200 MeV,which may be resulted from different devices and experimental conditions.In this perspective, an atmospheric-like neutron beam line with a wide spectrum is very suitable for the task of the cross-section measurements.In 2020, a pioneer work in energy dependence of SEU cross-section was reported by Iwashitaet al.[6]They exploited the hardware router and logic comparator within the chip to construct an energy-resolved measurement of SEU induced by the pulsed neutron source with tolerances in nanoseconds.Specifically,they focused on the soft errors that take place in the configuration random access memory(CRAM),which is a configurable static storage resource inside the chip.Circumstantial information indicates that commercial Xilinx FPGAs were adopted in the research.[7]Iwashita’s results suggest that the SEU cross-sections are steadily increased from 1 MeV to 20 MeV and show a plateau when the neutron energy is greater than 20 MeV.Apparently,these results are more elaborate than previous reports.However,the experimental design in Iwashita’s report is sophisticated.As the total size of CRAM in programming is out of control,authors introduced another cyclic redundancy check board to evaluate the total upsets, which makes the test complicated.

    The Back-n beam line is a white neutron test beam with a wide neutron spectrum range from 1 eV to 200 MeV at the China Spallation Neutron Source(CSNS).[8,9]The most probable energy of Back-n beam line is about 1 MeV and the portion of neutron energy greater than 1 MeV is about 48.2%.[10]Many SEE tests were accomplished at the Back-n neutron beam line and it has been a good alternative for semiconductor device tests before the atmospheric beam line (BL11) of CSNS completes its construction.Driven by 1.6 GeV protons,the atmospheric neutron irradiation spectrometer(ANIS)of CSNS is a promising test station to deliver a full spectrum neutron beam from MeV to GeV, which will be even more comparable to the territorial neutron environment than ICE House@LANSCE.Except the atmospheric beam line(BL11@CSNS)at 41°,Back-n beam line is the only fast neutron beam line that directs to the spallation target and outputs MeV neutrons without moderators.Back-n is initially designed for nuclear data measurements with a maximum flight path of 78.2 m which result in a good energy resolution in the range of resonance neutron spectrum.[10]The long flight path endows Back-n beam line an advantage in energy-resolved applications and research, such as energy-resolved neutron imaging.[11-13]As an atmospheric-like beam line,a comparison of Back-n beam line with other fast neutron beam lines for SEE irradiation in spectra may help to explain the situation.[14]Figure 1 lists the fluence of five fast neutron beam lines in the world and the JEDEC, IEC references.The label ‘Y’ means the sum of fluence in a year.As shown,the IEC*Y may serve as the upper limit among these lines while the JEDEC*Y is almost the prescribed minimum.FP30L@LANSCE has the most similar spectra shape to that of the JEDEC standard,which makes it the most favorable atmospheric neutron test site.The section in 10-200 MeV of TRIUMF NIF is the largest one,which is suggested to be helpful in the SEE tests.ANITA is a neutron facility driven by 180 MeV protons for accelerated single event effect tests,and is operated at the Svedberg Laboratory in Uppsala, Sweden.The spectra of VESUVIO@ISIS and Back-n@CSNS drop rapidly when the neutron energy is greater than 100 MeV.They are described as ‘soft’beam lines.In fact,the portion of Back-n flux from 1 MeV to 20 MeV is significantly higher than that of other white neutron sources.This may be helpful in the investigation of the threshold energy of SEU,and it also suggests that the Back-n beam line has a chance to get a better cross-section result at the energy of several MeV.In addition,Fig.1 exposes a fact that all five spectra are different from the realistic standard, and are just varying in degrees.

    To keep the test simple and compact, we prefer to measure directly the SEU bits in the block random access memory(BRAM)induced by neutrons.BRAM is a block memory that embedded throughout the FPGA which is often used for large data storage.They are continuous bits with excellent time performance.This feature is very useful when constructing modules for the time-of-flight measurement.The next step after determining the method and the objects is the technical scheme,including the chip selection and the circuit design.As a test in nanoseconds, the running clock of the circuit is important and will greatly influence the consuming time.A chip with high running frequency means a faster reading time and a smaller time shift error in SEU identification.Thanks to the advances in the semiconductor industry, commercial FPGAs with a frequency in hundreds MHz is very common nowadays,such as Kintex series in 28 nm from Xilinx Corp.Except the running frequency, the bit size of RAM for SEU tests is also related to the time consumed in the reading, which makes a great impact on the timing error in the time-of-flight (TOF)measurements.There are two effects caused by the bit size.On the one hand, under the same condition, the SEU occurrence is proportional to the used bit size.On the other hand,although a hardware parallel power could be exploited by circuit redesign, the size and depth of the RAM will increase the circuit complexity and eventually prolong the measurement time since the discrimination of the soft error consumes a certain time of locating and reading.Thus,an FPGA circuit with proper design is critical for the TOF measurements.In this test,we prefer to read a small memory of 400 kbit with a Xilinx chip running at a clock of 500 MHz.

    Fig.1.Comparison of the spectra of the CSNS Back-n white neutron source,the other white neutron sources and the JEDEC,IEC standards.The label‘Y’represents an accumulation of fluence in a year.

    2.Experimental measurement

    The experiment was performed at the CSNS Back-n beam line end station#1,which is 57.2 m far from the spallation center,as illustrated in the upright of Fig.2.The schematic in the left bottom of Fig.2 presents the experimental setup.Most of the devices are arranged in the underground tunnel and the experimental instructions are operated remotely in the control room.Fast neutrons are travelling from the left to the right,as indicated by the yellow horizontal arrow.The device under test(DUT,Kintex-7 chip)was installed in the center of the incident flange of a vacuum chamber,and been collimated by the laser markers.In order to timing the soft errors, another Spartan-6 FPGA communication board was fixed outside the beam spot.A local computer was setup in charge of the operation for the FPGA boards.The time-of-flight measurement was triggered by the signal from the rapid cycling synchrotron(RCS) kicker, as the red line shown in Fig.2.As shown, the kicker signal was extracted from the head of the ring target beam transport line and traveled near light speed with a short path.In this way,we constructed the trigger signal before the neutrons for TOF measurements.

    As a station for nuclear data measurement, the spectrum of Back-n beam line was measured accurately with a fission cross section measurement (FIXM) detector.[8]Started from its commissioning, several long-term spectral measurements demonstrate that the spectrum of Back-n beam line is very stable.Thus,no fission ion chamber was arranged after the DUT in this experiment and the spectrum data was inherited from its former measurements.[10]

    Fig.2.Schematic of the experimental setup at the CSNS Back-n white neutron source.Neutron beam transfer from the left to the right and the DUT was install in the center of the beam spot.

    2.1.SEU cross-section

    As mentioned,to estimate the DUT reliability,the soft error rate(SER)in a specific neutron irradiation is defined as an integration of SER within each energy section[6]

    whereσSEU(En)is the SEU cross-section at each energy sectionEn,Nbitsis the number of bits in the sample under test,andΦ(En)is the spectral fluence of neutron in units of n/MeV/cm2which designates the radiation fluence of the device in the test.Thus,the SEU cross-section is written as

    whereNSEU(En) is the number of soft error induced by neutrons within energy section ΔEnin the experiment.As illustrated in Fig.2,the soft error induced by neutrons will be registered by DUT output while its corresponding energy is measured by a timing process via time-of-flight method.

    2.2.Time-of-flight method

    Neutron TOF technique is a general method to determine the kinetic energy of a traveling neutron, by measuring the time it takes to fly from the spallation target surface to the circuit detector.The neutron kinetic energy is obtained from the time Δtand the distanceLby the following equation:

    wherem0is the rest mass of a neutron,cis the speed of light andv=L/Δtrepresents the velocity of the neutron,nonrelativistic in usual.

    In simple, we need to counter the SEUs and record its flight time Δt=tstop-tstart.However,the crux of a TOF measurement is to determine the start time and the end time of the neutron flight.Figure 3(a)shows the time relation within the measurements.As shown,it comes naturally that the timeTseu,when the soft error was registered,is set as the end time of the neutron flight.Since the maximum neutron energy of 200 MeV at the Back-n beam line is much lower than that in LANSCE, and the main section of Back-n spectrum is along the low energy section,as illustrated in Fig.1,the result is that the soft error counter histogram in time-space will have a slow increase,which suggests a large error if we determine the start time by the fastest neutrons.Therefore, in our measurements the start time cannot be treated like that in LANSCE, and a calibration of the start time(also written asT0)is demanded.

    Figure 3(b)presents an oscilloscope shot with annotations for theT0calibration.A vertical dotted line in purple represents the start timeT0when neutrons begin to fly.However,it is unknown and we need to figure it out.As mentioned,the CSNS beam profiling system provides bunch signals from the ring kicker, annotate asTk.Since the rising edge of kicker signals (Tk) is very stable, the kicker signal is often used to synchronize other devices after it.In our test,each TOF measurement was triggered by theTksignal,as the green line with legend CH2 plotted in Fig.3.When the test was being carried out, a C6D6 gamma detector was arranged in the end-station#2 to record the gamma flash signal, as shown in the red line with legend CH0.The gamma flash is a phenomenon of burst gamma emission which is concurrent with the spallation and the neutron flight.The blue line in Fig.3(b)is another trigger signal from a fast-current-transformer(FCT)which is not used for this analysis since it is later than the gamma flash.

    Fig.3.(a)The experimental configuration and time relation within the measurements.(b) T0 calibration for TOF measurements.The green line Tk is the kicker signal input to the FPGA circuit to trigger the TOF measurement.The prompt gamma burst of the first bunch is 3492 ns later behind Tk and the time for a photon to travel 76.2 m is 254 ns.

    As shown, the peak of the gamma flash caused by the first proton bunch is ΔTd=3492 ns behind theTksignal.In addition, since the detector is located 76.2 m far from the spallation target, a correction of the transportation time of gamma ΔTg= 254 ns from the target to the detector position is demanded.So, the start time of the neutron flight isT0=Tk+ΔTd-ΔTg.The flight time of a neutron soft error is written as

    whereTkis the kicker signal injected into the communication board to trigger the TOF measurements, andTseu-Tkis the output from the communication board.

    2.3.Circuit design for high-speed SEU detection

    Figure 4 presents the schematic diagram of the register transfer level (RTL) circuits designed for the cross-section measurements.As shown in Fig.4(a),BRAM modules are the core units in our design.40 modules are constructed in the test and each module outputs only one bit to flag whether a soft error occurs.The module bit size used in the test is 320×32 bits.In details, each module includes 5 BRAM blocks.They are configured as 512×72(four blocks)and 512×36(one block),respectively.So the total bit number is 512×(72×4+36).Therefore,partial bits(320 out of 324)and addresses(32 out of 512) are used to construct the soft error detection logic.Similar to Iwashita’s report,[6]the soft error probe is completed by logic comparators.As shown in the magnified drawing,all 320 bits are split into 20 groups.Initially,every bit is initialized in‘0’.To start with soft error detection,each group of 16 bits is input into a NOR gate.Then, the output from NOR signals will pass a two-level AND operation.Finally, a flag bit,whose change represents a upset occurs,is output for further process.The direct and parallel logic connection in the module results in a real time state detection.So whenever a bit upset within the BRAM module is detected,an error signal is output as a TOF stop which designates a soft error occured.In addition,the above process is implemented in a finite state machine(FSM),which describes the sequential logic.

    The circuit of the TOF board is presented in Fig.4(b).It is a Spartan-6 chip running at a frequency of 200 MHz.The kicker signalTkinstructs the start of the TOF measurements.An automatic counter realtime reg is cleared by the rising edge of theTksignal every 40 ms, indicating the start of a new round.Whenever the rising/falling edge of the error signal occurs,the real time register will latch the value of the realtime reg counter to get the neutron TOF time induced by the soft error.Finally, the time difference (Tseu-Tk) as the timing result is sent to PC via an UART module.Figure 4(c)analyzes the time sequence during a whole SEU detection.In fact, there are various time delays in the signal processes which should been calculated carefully.Suppose that the two asynchronous signals (kicker and error) have the same time difference compared to the clock edge.As shown in Fig.4(b),the regtime reg counter is cleared by the kicker signal with two cycles delay,which means 10 ns will be consumed at the beginning of the TOF.However,before the SEU event is registered,another consumed time should be evaluated.The SEU flag signal requires a delay of 3 cycles in the Kintex-7 chip and 2 cycles in the S6 board before the value of the reltime reg counter is latched.This process takes 16 ns before a SEU event is registered, as shown in Fig.4(c), which means that the neutron TOF has been increased by 16 ns.Therefore, the time delay in the circuits is 16-10=6 ns, which should be subtracted in the final result.In the end,the circuit was tested before the real experiment.Figure 4(d)presents our HDL simulation result.As shown,after the emulated upset event is injected,each process in Kintex-7(FSM,comparator and DFF)takes 2 ns and an error signal is successful raised after that.However, the time difference between the two asynchronous signals relative to the clock edge cannot always be equal, so there is a maximum time shift of±5 ns in our formula,which will eventually be expressed as an energy uncertainty in the following calculation.

    Fig.4.Schematic diagram of RTL for high-speed SEU detection and TOF measurement circuits.(a)Kintex-7 FPGA,it is the device under test(DUT),which operates at 500 MHz and outputs a signal that triggers the TOF calculation when a soft error occurs.(b)Spartan-6 chip,running at a frequency of 200 MHz,and outputs the neutron TOF time.(c)Analysis of FPGAs internal signal delay.(d)HDL simulation.

    3.Result analysis and discussion

    The experiment took 104.5 hours with 908 effective SEU records in total.The time-of-flight count histogram shows a ragged edge due to the scarce events, as shown in Fig.5(a).As predicted,the soft error counts are increasing much slower than that in Iwashita’s report.[6]The count curve reaches its peak of 44 at a time about 1000 ns,and then they are decreasing steadily.

    Fig.5.Results of the time-of-flight spectra for the SEU counts and the unfold process.The total soft errors in DUT are 908 counts.

    The CSNS is a pulsed neutron source running at a frequency of 25 Hz.Meanwhile, the time structure of its driven proton consists of two bunches, with a separation of 410 ns.It is a common design for modern accelerator and spallation station to increase its yield for hot/cold neutron applications.However,this time structure makes a big obstacle for the TOF measurements,especially for fast neutron.So,the most tricky work is the unfolding of the data since it is hard to accumulate enough counts(such as millions of events within 100 bins per degree for flux measurements) in the SEE experiment in order to satisfy the statistics requirement for the process.[15]Therefore, we drop the histogram concept and the minimal unit of the unfolding falls to single event.However,neither it is scarce enough to be completed with eyes.Thus, we developed a new method to tackle the separation process.The basic hypothesis of the unfolding is that the two hunches will have the same physical consequences except a time separation.The fastest SEU event must be attributed to the first bunch while its corresponding event in the second bunch should be 410 ns later after it.This process is called a pairing.In practice, the nearest event to that time tick will be assigned to the second bunch.Therefore,the first unchecked event is assigned to the first bunch, while the nearest event after 410 ns to the second bunch,and then the pairing is repeated.Apparently,there will be a difference between the real position and the expected position when looking for the nearest event from the second bunch.However,that difference is indicative.We kept that as a residue which will be balanced in the next searching process.The unfolded result is presented in Fig.5(b).The blue and red lines are the SEUs attributed to the first bunch and the second bunch,respectively,with a time correction.As shown,the first and the second bunches show a good overlapping, especially in the first peak region.The line in Fig.5(c) represents the ratio of the time residue with respect to the separation interval.The up and down fluctuations around thex-axis suggest that the balance technique helps to control the residue during the pairing.As shown, there are three fluctuation regions after the unfolding process.The most ragged position is near the 400th pair, which corresponds to the second peak region in Fig.5(b).This is consistent with the data in Fig.5(b).The result demonstrates that our unfolding is coherent and correct.Figure 5(d)shows the final result after unfolding.The unfolding process makes a big influence on the distribution,turning it from a relative smooth distribution into an ragged one.

    Fig.6.Measured neutron spectral fluence Φ(En).

    Figure 6 presents the fluence during the test,where a maximum fluence over 1×1012n/cm2/MeV was achieved.The bit SEU cross-section ranging from 0.5 MeV to 200 MeV is presented with the red dots and line in Fig.7,including the error bar.As shown,the large time uncertainty results in a wide energy uncertainty,which is concerned in high energy region.However,an overall trend is upward as a function of the neutron energy.For clarity, a magnified picture is inserted in the bottom center to show the details in the energy section from 0 to 10 MeV.The longest TOF time is 5955-410=5545 ns,which corresponds to a neutron single upset event with energy 557 keV (second bunch), while its pairing event in the first bunch is 5350 ns/598 keV.The result confirms that the SEU occurred much below the threshold energy of the nuclear cross-section of silicon(2.75 MeV),which is consistent with Iwashita’s report.Such a low energy is beyond the conventional recognition, which raises the hypothesis of its relationship to the manufacturing scale.Further investigations in both experimental tests and Monte Carlo simulations are demanded.In addition, the time-of-flight in the energy section near 0.5 MeV is about 5850 ns, where the energy resolution is little than 5% at the Back-n beam line.[10]Thus, it is also demonstrated that the measurement with a pulsed spallation neutron source is a reliable method for the critical energy investigation.Since the fluence data are accumulated in a longterm measurement, the contribution from fluence can be ignored,and the statistical errors of the cross section are mainly determined by the event counts.The difficulty in this SEU test is the little event count, which is caused by the small cross section.Since only one DUT was installed,the calculation of the margin of error relies upon the sample size, which is the count numbers in the bins.As shown in the unfolded ToF result in Fig.5(d),the count number in the high energy region is near 80,which results in an error of about 11%.Thus the high energy region shows a small statistical error while a relatively obvious error bar in the low energy region.For an average energy of 5.6 MeV,however,the unfolding process dramatically reduces the count number, which illustrates a remarkable error in the insert.The influence of the unfolding process will be discussed later.In addition,the highest data points display an obvious deviation,which may be attributed to the division process since the flux drops very quickly when the neutron energy is close to 200 MeV.It will be excluded in the following fitting process.

    One detail that is easily overlooked in Iwashita’s report is the start energy in the cross-section plot which was fixed as 1 MeV,this is incoherent with their TOF histogram.As illustrated,the interval of proton bunch in WNR station is fixed as 1.8μm, and the soft error event is not ended before 1800 ns.The maximum time interval of 1800 ns in the TOF histogram designates a neutron energy of 650 keV.Thus, as a matter of fact, Iwashitaet al.had measured a threshold energy, which is very close to our result.However,the discussion in Ref.[6]was short and soft errors lower than 1 MeV were ignored.A possible explanation is that these data were regarded as unreliable for the statistics reason.Nevertheless,the critical energy for neutron SEU events and its relation with manufacturing scale are very important questions.Obviously, a conclusion that the threshold energy less than 2.75 MeV is insufficient.The real reason why they gave up on the discussion may be that they noticed the threshold energy is out of their capability,and the discussion on the critical energy can not be completed.Thus, the Back-n beam line has the chance in this subject.In fact, the long flight path of the Back-n beam line meets the TOF measurement down to epicadmium neutron, which means it is possible to catch the soft error caused by even low energy neutrons.In addition,the interval of 1.8μm discussed above is originated from micro-pulse setting at WNR beam line,which is a counterpart to our double bunch mode.However,the superposition of events was not mentioned.

    To evaluate the result, we have to compare it with other reports.Since the design rule has a great impact on the SEU cross-section and the manufacturing node has made rapid progress over the years, we have to focus on the literature of the last five years.Obviously, the report from Iwashitaet al.cannot be ignored.As shown, the cross-section data from Iwashitaet al.plotted in blue asterisks are about an order of magnitude smaller than our results.Another SRAM neutron SEU measurement with mono-energy neutron 2.5 MeV and 14 MeV from Chenet al.is considered here.[16]Their DUTs were manufactured with the latest design rule of 40 nm,which attracts our interest for the comparison and other details are ignored.As shown, their SEU cross-sections are much higher than our result in the same energy.In addition,soft error test data are important indices in order to help customers design products for the high reliability market,and also many chip manufacturers publicly release their reliability data results.Xilinx Corp.has published its soft error data in their chip user guide, including both CRAM and BRAM.[17]As listed in Table 20 of Ref.[17], the BRAM in Kintex-7 chip has a cross-section of 5.57×10-15cm2per bit with an error of±18% according to the LANSCE test, while the CRAM cross-section is 5.69×10-15cm2/bit,very close to that of the BRAM.Although these results are only coarse cross-sections in total spectrum, they may serve as an average of the SEU cross-sections for our comparison.The manufacturer BRAM cross-section is plotted as the horizontal purple dash-dotted line in Fig.7.As shown, the average cross-section line from the Xilinx is higher than the blue asterisks from Iwashitaet al.A common sense is that an average value of many sections should be in the middle of the whole data.Thus,considering the authority of the manufacturer data and their sufficient test in beam time,the result from Iwashitaet al.is doubtful if they did use the same Xilinx series chip in the tests.In our opinion, the complex experimental design should account for the situation.Thus, although there are some fluctuations, the red dotted result in our test is a reasonable result.

    Another valuable topic is the subject of test data equivalence between the Back-n beam line and the FP30L@LANSCE.Some former report also touched the issue with a coarse estimate.Here we may repeat it in order to have an extensive comprehension.In our case,a coarse estimate of SEU cross-section with respect to the Back-n beam line will result inσseu=N/(Φ·tbits)=908/(1.62×107×η×T×400 kb)=7.55×10-16cm2/bit,whereη=48.2%represents the ratio of neutron with energy great than 1 MeV.[10]Since the spectrum of the Back-n neutron beam line is much different from that of LANSCE, these two experiments cannot be compared directly.Nevertheless, a previous experimental report may help in the issue.[18]Wanget al.had compared the SEUs result from the Back-n beam line and the JEDEC atmospheric environment.They concluded that the‘soft’characteristics of Back-n beam line will exaggerate the beam fluence which decreases the SEU cross-section.As listed in Table 6,[18]for HM62V8100 SRAM and considering a threshold energy of 0.6 MeV, a correction factor of 5.56 was suggested by the authors.Therefore, if the difference between LANSCE and the JEDEC standard is ignored, the total SEU cross-section of 4.2×10-15cm2/bit in our test is close to the result(about 5.6×10-15cm2/bit)provided by the chip manufacturer.In short,importing a correction factor to convert the SEU cross-section between different test stations is an effective method for a coarse data equivalence.The above estimate also suggests that no systematic error was included in our result.In addition, the measured threshold energy of 0.5 MeV also confirms the speculative energy value in report.[18]

    As suggested,[19,20]the SER cross-section is usually formulated in the 4-parameters Weibull function

    whereσsrepresents the saturated cross-section andEcis the critical energy.The fitted line is plotted in black(dashed line)and the parameters are labeled upon the line in the upper center of Fig.7.The coefficient of determinationR-square=0.9386 suggests the fitting is fine, except for some fluctuations.The critical energyEc=0.545 MeV is consistent with our experimental observation.

    Although the fit shows a good result,there are some fluctuations in the bit cross-section curve which are related to the valleys in the counter histogram of Fig.5(d).Since the coarse estimate is consistent with the manufacturer’s data, the total SEU counts are as expected.Therefore,an explanation is that events in this energy section were counted in its neighbor bins due to the TOF timing error.With this in mind, we explore the time shift in neutron energy.There are two sources of time uncertainty in Eq.(4), including both the start timeT0and the end timeTseu.The uncertainty in the start time should be attributed to the time structure of protons at CSNS.Except the separation of the two proton bunches,the pulse width of a proton bunch is about 80 ns at the running power of 100 kW.The spallation neutron bunch is suggested to have an equivalent pulse length for any energy section.Therefore, since it is impossible to distinguish whether the neutron comes from the proton bunch head or its tail,the TOF start has a time uncertainty of 80 ns.The second time uncertainty comes from the reading process.Since the target BRAM is 32 bit in depth,the reading process will consume 32 cycles which means a time uncertainty of 64 ns.In addition, a time shift of 5 ns in the asynchronous case of the Spartan-6 board had been mentioned.Therefore,the maximum time uncertainty in total will be 80+64+5=149 ns.The error bar inx-axis in Fig.7 presents the corresponding energy uncertainty.As shown,the fixed time uncertainty expresses a large energy uncertainty in the high energy region while a much shorter error bar is shown in the range of 1-10 MeV.However,the insert also shows two points deviate the fit curve,which is corresponding to the valleys in Fig.5,and we will explain its origin below.

    Fig.7.Result of neutron SEU cross-section for the SRAM-based FPGA.The red dotted line is the bit SEU cross-section measured in this test.The blue asterisks are the result from Iwashita’s report in LANSCE.[6] The green diamonds are the mono-energy cross-sections results in Ref.[16].

    Since it is impossible to get the accurate time, the event may occur in every tick in a time window with a width similar to the total uncertainty.So,the measurement uncertainty may act like a moving window which smooths the counts in time series.In common, the stacking under double-bunches mode will result in an obvious shoulder shape in its time histogram.But, the moving effect of the time uncertainty smooths the curve in to a relative gentle one.Unfortunately,the consequential double-bunches unfolding process breaks the smoothed count histogram into a ragged state with peaks and valleys.In short, both the time structure of double bunches and the moving window should be responsible for the fluctuations.Although the CSNS accelerator has a single bunch mode,the cost is that the neutron flux at each terminal will be reduced by half.Obviously, it is hard for other beam lines to accept the single bunch mode as a long term running state.Although the single bunch beam time is very scarce, we wish to make a test with single bunch mode in the future.

    In addition, although from the perspective of a single event, the bunch interval of 410 ns is like a big time uncertainty, it is should be regarded as a ‘fog’ before the truth and the underlying physics could be unveiled with the unfolding process.This has been proven by many experiments, such as time-of-flight spectrum test and nuclear data measurement.[8,10,15]The only difference is the size of statistics.Our test had collected 908 events in nearly 5 days and the size number is decent.Therefore,although the data set is not good enough, the result is acceptable for a SEE test.Meanwhile,for a small cross section and/or little fluence,the small number is inevitable at the very head or tail.The small number of soft error in low energy also has important physical meaning,otherwise,the investigate and discuss of the threshold energy will fall in vain.As to the reading depth of 32 bit,since we have a 80 ns proton bunch, a reading time uncertainty of 64 ns is acceptable.As an attempt,we prefer to a small number of BRAM blocks with large reading depth to increase the bit size,which helps us in the RTL timing constraint.

    Another important information in Fig.7 is the crossing point of the fitted line and the‘a(chǎn)verage’line.The correspondingx-tic of the crossing point suggests an ‘effective’ neutron energy at about 30.9 MeV, where the SEU cross-section induced by mono-energy neutron source will be equivalent to the one from ICE House at LANSCE.This is an important information to evaluate the similarity between a beam line and the atmospheric neutron environment if we regard the FP30L@LANSCE as the standard beam line for semiconductor device irradiation.However, the results from Chenet al.indicate that the effective neutron energy is below 14 MeV for their DUTs.In our opinion, it is unlikely that the effective neutron energy is as small as that,while some dozens of MeV is more probably in our guess.This should be checked in the future experiments.Further SEU cross-section measurements with mono-energy neutron higher than 30.9 MeV will help to verify the statement.If more measurements confirm our result with the concept of the effective neutron energy, it will be a new index to evaluate the characteristics of a beam line whether it is‘hard’as the atmospheric test beam line.

    4.Conclusions

    In summary,an energy-resolved single-event-upset crosssection measurement was performed with 28 nm BRAMs using TOF techniques at the Back-n beam line at CSNS.The result confirms that the soft error threshold energy is about 0.5 MeV,which is consistent with previous report.The comparison suggests that the direct measurement is a reliable method under this circumstance and the measured cross section is consistent with the manufacturer reliability data.In addition, an ’effective’ neutron energy is proposed as an index in order to evaluate the characteristics of an atmospheric test beam line.

    Acknowledgements

    Project supported by the National Natural Science Foundation of China (Grant Nos.2032165 and 62004158), the National Key Scientific Instrument and Equipment Development Project of China (Grant No.52127817), the State Key Laboratory of Particle Detection and Electronics(Grant Nos.SKLPDE-ZZ-201801 and SKLPDE-ZZ-202008),and the Special Funds for Science and Technology Innovation Strategy of Guangdong Province, China (Grant No.2018A0303130030).

    猜你喜歡
    永寧小龍
    El regreso del dragón
    辛永寧:慢性乙肝患者隨訪復(fù)查的那些事兒
    肝博士(2022年3期)2022-06-30 02:48:24
    永寧之門(mén)
    小小小小龍
    劉小龍
    讓“數(shù)”“形”結(jié)合更暢通
    我愛(ài)冬天
    生物胺與食品安全
    風(fēng)中的祈禱詞
    詩(shī)選刊(2015年4期)2015-10-26 08:45:28
    可愛(ài)的小龍人
    欧美区成人在线视频| 蜜臀久久99精品久久宅男| 久久久久久久精品精品| 男女无遮挡免费网站观看| 一二三四中文在线观看免费高清| 亚洲第一av免费看| 视频中文字幕在线观看| 男的添女的下面高潮视频| 亚洲不卡免费看| 激情 狠狠 欧美| 成年美女黄网站色视频大全免费 | 久久精品夜色国产| 少妇被粗大猛烈的视频| 欧美一级a爱片免费观看看| 青春草国产在线视频| 色哟哟·www| 欧美日韩视频精品一区| 午夜免费男女啪啪视频观看| 一级av片app| 我要看黄色一级片免费的| 亚洲精品色激情综合| 国产精品国产三级专区第一集| 欧美高清成人免费视频www| 欧美成人精品欧美一级黄| 2018国产大陆天天弄谢| 国产黄色免费在线视频| 久久久久国产网址| av在线app专区| 国产成人aa在线观看| 99久久精品热视频| tube8黄色片| 午夜福利视频精品| 成人国产麻豆网| 久久国产精品大桥未久av | 人妻制服诱惑在线中文字幕| 2021少妇久久久久久久久久久| 热re99久久精品国产66热6| 成人黄色视频免费在线看| 蜜臀久久99精品久久宅男| 身体一侧抽搐| 日本vs欧美在线观看视频 | 日韩强制内射视频| 2022亚洲国产成人精品| 下体分泌物呈黄色| 久久久久人妻精品一区果冻| 天天躁夜夜躁狠狠久久av| 日本vs欧美在线观看视频 | 我要看日韩黄色一级片| 一级毛片我不卡| 国产精品久久久久久精品电影小说 | 国产 精品1| av黄色大香蕉| 亚洲性久久影院| 免费久久久久久久精品成人欧美视频 | 精品人妻熟女av久视频| 欧美三级亚洲精品| 热99国产精品久久久久久7| 乱码一卡2卡4卡精品| 免费观看的影片在线观看| 久久影院123| 最近手机中文字幕大全| 亚洲av成人精品一区久久| 欧美3d第一页| av黄色大香蕉| 国产有黄有色有爽视频| 新久久久久国产一级毛片| 亚洲aⅴ乱码一区二区在线播放| 国产高潮美女av| 一区二区三区乱码不卡18| 亚洲欧美一区二区三区国产| 亚洲精品视频女| 成人美女网站在线观看视频| 嫩草影院入口| 国产精品一区二区性色av| 美女视频免费永久观看网站| 欧美xxxx性猛交bbbb| 中国国产av一级| 激情五月婷婷亚洲| 激情 狠狠 欧美| 天美传媒精品一区二区| 亚洲精品久久午夜乱码| 欧美成人午夜免费资源| 一本一本综合久久| 日产精品乱码卡一卡2卡三| 91久久精品国产一区二区成人| 久久99热这里只频精品6学生| 中文字幕亚洲精品专区| 中文字幕亚洲精品专区| 午夜日本视频在线| 成人无遮挡网站| 多毛熟女@视频| 婷婷色综合www| 亚洲精品一区蜜桃| 两个人的视频大全免费| 日韩av免费高清视频| 99久国产av精品国产电影| 中文字幕人妻熟人妻熟丝袜美| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品乱久久久久久| 国产大屁股一区二区在线视频| 一本久久精品| av免费在线看不卡| 久久精品人妻少妇| 男女啪啪激烈高潮av片| 久久热精品热| kizo精华| av在线app专区| 国产伦精品一区二区三区视频9| 99国产精品免费福利视频| 亚洲欧美中文字幕日韩二区| 成人亚洲欧美一区二区av| 国产精品三级大全| 日本欧美视频一区| 日韩一区二区视频免费看| 在线免费观看不下载黄p国产| 精品人妻一区二区三区麻豆| 国产av码专区亚洲av| 亚洲伊人久久精品综合| 男人和女人高潮做爰伦理| 亚洲人成网站高清观看| 男的添女的下面高潮视频| 97在线视频观看| 男人爽女人下面视频在线观看| 久久国产精品男人的天堂亚洲 | 大片免费播放器 马上看| 高清在线视频一区二区三区| 夜夜看夜夜爽夜夜摸| 精品午夜福利在线看| 日本欧美国产在线视频| 五月天丁香电影| 国产成人freesex在线| 国产成人a区在线观看| 国产黄色视频一区二区在线观看| 中文乱码字字幕精品一区二区三区| 久久99热6这里只有精品| 亚洲,欧美,日韩| 成人黄色视频免费在线看| 天天躁日日操中文字幕| 少妇精品久久久久久久| 亚洲av国产av综合av卡| 久久国产亚洲av麻豆专区| 精品亚洲乱码少妇综合久久| 成年美女黄网站色视频大全免费 | 在线天堂最新版资源| 国产永久视频网站| 国产黄频视频在线观看| 1000部很黄的大片| 自拍偷自拍亚洲精品老妇| 亚洲av福利一区| 久久人人爽av亚洲精品天堂 | 国产毛片在线视频| 女人久久www免费人成看片| 色综合色国产| 国产一区二区在线观看日韩| 性色avwww在线观看| 夜夜爽夜夜爽视频| 欧美日韩精品成人综合77777| 欧美日韩国产mv在线观看视频 | 中国国产av一级| 纯流量卡能插随身wifi吗| 国产男女超爽视频在线观看| 国产精品人妻久久久影院| 黑人猛操日本美女一级片| 日日摸夜夜添夜夜添av毛片| 亚洲丝袜综合中文字幕| 国产淫语在线视频| 看免费成人av毛片| 成人亚洲欧美一区二区av| 不卡视频在线观看欧美| 你懂的网址亚洲精品在线观看| 国产永久视频网站| 亚洲av福利一区| 国产精品一区二区性色av| 亚洲欧美日韩无卡精品| 成人特级av手机在线观看| 久久久久精品性色| 久久久a久久爽久久v久久| 夜夜看夜夜爽夜夜摸| 中国三级夫妇交换| 一本久久精品| 精品人妻一区二区三区麻豆| 久久鲁丝午夜福利片| 亚洲欧美一区二区三区黑人 | 国产老妇伦熟女老妇高清| 少妇 在线观看| 九草在线视频观看| 一级毛片电影观看| 国产在线一区二区三区精| 国产精品一区二区性色av| 国产综合精华液| 中文在线观看免费www的网站| 女人久久www免费人成看片| 免费高清在线观看视频在线观看| 这个男人来自地球电影免费观看 | 亚洲国产精品专区欧美| 国产老妇伦熟女老妇高清| 黄色欧美视频在线观看| 人妻 亚洲 视频| 小蜜桃在线观看免费完整版高清| 精品国产一区二区三区久久久樱花 | 国产精品久久久久久久久免| videossex国产| 毛片女人毛片| 色5月婷婷丁香| 国产成人91sexporn| 少妇人妻精品综合一区二区| 久久国产精品大桥未久av | 97精品久久久久久久久久精品| 久久鲁丝午夜福利片| 一级片'在线观看视频| 如何舔出高潮| 亚洲精品亚洲一区二区| 夜夜骑夜夜射夜夜干| 亚洲精品一二三| 九九久久精品国产亚洲av麻豆| 97精品久久久久久久久久精品| 波野结衣二区三区在线| 国产亚洲91精品色在线| 人妻少妇偷人精品九色| 最近最新中文字幕免费大全7| 黑人高潮一二区| 精品一区二区三区视频在线| 欧美激情极品国产一区二区三区 | 亚洲人成网站高清观看| 一区二区av电影网| 亚洲精品456在线播放app| 久久人妻熟女aⅴ| 天美传媒精品一区二区| 久久久久久久久久久丰满| 国产成人精品一,二区| 亚洲国产精品专区欧美| xxx大片免费视频| 久久久精品94久久精品| 青春草亚洲视频在线观看| 国产高清有码在线观看视频| 国产精品久久久久久久久免| 22中文网久久字幕| 91精品伊人久久大香线蕉| 日韩欧美一区视频在线观看 | 91精品国产九色| 一级二级三级毛片免费看| 久久国内精品自在自线图片| 亚洲精品一区蜜桃| 亚洲精品aⅴ在线观看| 久久 成人 亚洲| 久久久欧美国产精品| 精品久久久久久电影网| 国产精品一二三区在线看| 一级毛片黄色毛片免费观看视频| 午夜免费鲁丝| 97超碰精品成人国产| 亚洲自偷自拍三级| 中文在线观看免费www的网站| 国产黄频视频在线观看| 欧美精品亚洲一区二区| 插逼视频在线观看| 在线观看免费高清a一片| 国产在线男女| 精品久久国产蜜桃| 久久久久久久久久久丰满| 街头女战士在线观看网站| 国产淫片久久久久久久久| 伦理电影大哥的女人| 中国国产av一级| 精品久久久噜噜| 欧美bdsm另类| 国产无遮挡羞羞视频在线观看| 边亲边吃奶的免费视频| 亚洲av免费高清在线观看| 免费大片黄手机在线观看| 亚洲欧洲日产国产| 免费少妇av软件| 少妇丰满av| 亚洲精品第二区| av.在线天堂| 在线观看免费高清a一片| 日日摸夜夜添夜夜添av毛片| 少妇的逼好多水| 尤物成人国产欧美一区二区三区| 久久国产乱子免费精品| 超碰97精品在线观看| 国产高潮美女av| 18禁裸乳无遮挡免费网站照片| 国产欧美日韩精品一区二区| 精品国产一区二区三区久久久樱花 | 校园人妻丝袜中文字幕| 一区二区三区乱码不卡18| 男人狂女人下面高潮的视频| 成人美女网站在线观看视频| 中文资源天堂在线| 日韩人妻高清精品专区| 777米奇影视久久| 亚洲成人手机| 狂野欧美激情性bbbbbb| 国产日韩欧美亚洲二区| 少妇丰满av| 久久精品久久久久久久性| 国产在线一区二区三区精| 婷婷色av中文字幕| 建设人人有责人人尽责人人享有的 | 欧美另类一区| 丰满迷人的少妇在线观看| 18禁动态无遮挡网站| 日韩av在线免费看完整版不卡| 欧美性感艳星| 欧美国产精品一级二级三级 | 亚洲欧洲国产日韩| 免费不卡的大黄色大毛片视频在线观看| 久久综合国产亚洲精品| a级毛色黄片| 欧美三级亚洲精品| 春色校园在线视频观看| 夜夜看夜夜爽夜夜摸| 欧美成人a在线观看| 99re6热这里在线精品视频| 狂野欧美激情性xxxx在线观看| 欧美97在线视频| 亚洲国产精品专区欧美| 国产人妻一区二区三区在| 中文资源天堂在线| av网站免费在线观看视频| 搡女人真爽免费视频火全软件| 国产精品精品国产色婷婷| 国模一区二区三区四区视频| 日本黄色片子视频| 亚洲不卡免费看| 免费大片黄手机在线观看| 国产精品精品国产色婷婷| 国产日韩欧美在线精品| 久久久午夜欧美精品| 九色成人免费人妻av| 国产在视频线精品| 国产v大片淫在线免费观看| 美女中出高潮动态图| 亚洲欧美精品自产自拍| 美女主播在线视频| 中国国产av一级| 精品视频人人做人人爽| 亚洲国产av新网站| 狂野欧美激情性bbbbbb| 欧美老熟妇乱子伦牲交| 制服丝袜香蕉在线| 3wmmmm亚洲av在线观看| 青春草亚洲视频在线观看| 一级黄片播放器| 亚洲三级黄色毛片| 欧美xxⅹ黑人| 性色avwww在线观看| 国产精品99久久久久久久久| 久久精品人妻少妇| 一个人看的www免费观看视频| 中文天堂在线官网| 日韩在线高清观看一区二区三区| 成年美女黄网站色视频大全免费 | 91精品国产国语对白视频| 日本黄色片子视频| 国产黄片美女视频| 国产精品免费大片| 国产人妻一区二区三区在| 老女人水多毛片| 精品熟女少妇av免费看| 日本猛色少妇xxxxx猛交久久| 久久人人爽人人爽人人片va| 亚洲国产精品一区三区| 亚洲三级黄色毛片| 又大又黄又爽视频免费| 国产片特级美女逼逼视频| 欧美激情国产日韩精品一区| 一本色道久久久久久精品综合| 中文字幕免费在线视频6| 国产永久视频网站| 一级毛片我不卡| 亚洲人成网站在线观看播放| 菩萨蛮人人尽说江南好唐韦庄| 国产成人aa在线观看| 亚洲久久久国产精品| 97超视频在线观看视频| 久久人人爽av亚洲精品天堂 | 一区二区三区乱码不卡18| 亚洲av综合色区一区| 黄色一级大片看看| 99久久精品热视频| 97热精品久久久久久| 大又大粗又爽又黄少妇毛片口| 九草在线视频观看| 国产精品国产三级国产av玫瑰| 最近2019中文字幕mv第一页| 一个人看的www免费观看视频| 男男h啪啪无遮挡| 黄色日韩在线| 午夜福利影视在线免费观看| 久久人妻熟女aⅴ| 国产老妇伦熟女老妇高清| 亚洲精品aⅴ在线观看| 夜夜骑夜夜射夜夜干| 看免费成人av毛片| 多毛熟女@视频| 最近手机中文字幕大全| 一级爰片在线观看| 日本免费在线观看一区| 国产伦精品一区二区三区视频9| 久久影院123| 极品教师在线视频| 亚洲第一区二区三区不卡| 国产精品嫩草影院av在线观看| 永久免费av网站大全| videossex国产| 国产探花极品一区二区| 亚洲国产欧美人成| 建设人人有责人人尽责人人享有的 | 国产精品不卡视频一区二区| 精品人妻视频免费看| 久久久久网色| 日韩,欧美,国产一区二区三区| 亚洲精品日韩在线中文字幕| 成人漫画全彩无遮挡| 99热这里只有是精品在线观看| 91精品一卡2卡3卡4卡| 亚洲欧美成人综合另类久久久| 久久韩国三级中文字幕| 99久久人妻综合| 最新中文字幕久久久久| 新久久久久国产一级毛片| 五月玫瑰六月丁香| 久热这里只有精品99| 91久久精品国产一区二区三区| 久久婷婷青草| 一级毛片 在线播放| 狂野欧美激情性bbbbbb| 秋霞伦理黄片| 久久久午夜欧美精品| 一级a做视频免费观看| 色婷婷久久久亚洲欧美| 国产成人精品福利久久| 丰满迷人的少妇在线观看| 国产乱来视频区| 岛国毛片在线播放| 美女xxoo啪啪120秒动态图| 久久亚洲国产成人精品v| 深爱激情五月婷婷| 黄色配什么色好看| 亚洲婷婷狠狠爱综合网| 日韩国内少妇激情av| 欧美国产精品一级二级三级 | 在线观看免费高清a一片| 在线亚洲精品国产二区图片欧美 | 七月丁香在线播放| 亚洲,一卡二卡三卡| 人人妻人人看人人澡| 久久影院123| 亚洲精品一二三| 一级片'在线观看视频| 午夜精品国产一区二区电影| av.在线天堂| 99九九线精品视频在线观看视频| 午夜激情福利司机影院| 日韩欧美一区视频在线观看 | 亚洲精品aⅴ在线观看| 色哟哟·www| 99久久综合免费| 日本欧美视频一区| 草草在线视频免费看| 婷婷色av中文字幕| 少妇人妻久久综合中文| 日韩一区二区视频免费看| 少妇裸体淫交视频免费看高清| 精品酒店卫生间| 我要看黄色一级片免费的| 久久精品人妻少妇| 久久久久久久国产电影| 欧美精品国产亚洲| 欧美激情国产日韩精品一区| 日韩欧美精品免费久久| 亚洲精品视频女| 中文字幕制服av| 久久久久视频综合| 免费看不卡的av| 久久综合国产亚洲精品| 少妇的逼水好多| 亚洲av男天堂| 韩国av在线不卡| 亚洲电影在线观看av| 日韩成人av中文字幕在线观看| 18禁动态无遮挡网站| 亚洲电影在线观看av| 久久青草综合色| 蜜桃在线观看..| 国产精品一区二区在线不卡| 黄色日韩在线| 乱系列少妇在线播放| 哪个播放器可以免费观看大片| 老师上课跳d突然被开到最大视频| 一级毛片电影观看| 久久久久久久久久久免费av| 乱系列少妇在线播放| av国产免费在线观看| 亚洲人成网站高清观看| 青青草视频在线视频观看| 国产av精品麻豆| 免费黄色在线免费观看| 亚洲怡红院男人天堂| 国产色婷婷99| 免费观看性生交大片5| 天堂中文最新版在线下载| 青春草国产在线视频| 久久久成人免费电影| 人妻 亚洲 视频| 一二三四中文在线观看免费高清| 永久免费av网站大全| 如何舔出高潮| 国产亚洲欧美精品永久| 99re6热这里在线精品视频| 亚洲综合精品二区| 日日啪夜夜撸| 国产日韩欧美亚洲二区| 在线观看国产h片| h日本视频在线播放| 老司机影院成人| 一区在线观看完整版| 成人美女网站在线观看视频| 日韩成人av中文字幕在线观看| 欧美日韩一区二区视频在线观看视频在线| 国产精品一二三区在线看| 久久久久久久久大av| 97超碰精品成人国产| 欧美日韩国产mv在线观看视频 | 韩国av在线不卡| 日本色播在线视频| 中文字幕亚洲精品专区| 久久国产乱子免费精品| 麻豆乱淫一区二区| 午夜福利高清视频| 亚洲精品日本国产第一区| 精品国产一区二区三区久久久樱花 | 国产成人91sexporn| 国产欧美日韩一区二区三区在线 | 男的添女的下面高潮视频| 亚洲国产精品成人久久小说| 午夜福利网站1000一区二区三区| 日韩免费高清中文字幕av| 久久久久久久久久人人人人人人| 久久97久久精品| 人人妻人人看人人澡| 中文精品一卡2卡3卡4更新| 亚洲不卡免费看| 国产 一区精品| 一级黄片播放器| 亚洲精品亚洲一区二区| 日韩 亚洲 欧美在线| 日本午夜av视频| 日韩电影二区| 天天躁日日操中文字幕| 汤姆久久久久久久影院中文字幕| 国产亚洲av片在线观看秒播厂| 国产精品一区二区在线不卡| 国产av码专区亚洲av| 久久精品久久久久久久性| 精品午夜福利在线看| 免费看不卡的av| 校园人妻丝袜中文字幕| 大香蕉久久网| 日本av手机在线免费观看| 人人妻人人爽人人添夜夜欢视频 | 欧美日韩综合久久久久久| 18禁在线无遮挡免费观看视频| 网址你懂的国产日韩在线| av国产精品久久久久影院| 伊人久久精品亚洲午夜| 老司机影院成人| 日韩精品有码人妻一区| 国产精品人妻久久久久久| 高清不卡的av网站| 久久国产精品大桥未久av | 性色avwww在线观看| 国产免费一级a男人的天堂| 人妻少妇偷人精品九色| 国产极品天堂在线| 91精品国产九色| 国产伦精品一区二区三区视频9| 亚洲美女搞黄在线观看| 亚洲熟女精品中文字幕| 91精品国产九色| 亚洲av中文av极速乱| 在线 av 中文字幕| 欧美精品一区二区大全| 久久精品国产亚洲网站| 国产亚洲5aaaaa淫片| 亚洲人与动物交配视频| 日日啪夜夜撸| 国产欧美日韩一区二区三区在线 | 免费观看在线日韩| 在线亚洲精品国产二区图片欧美 | 欧美97在线视频| 国产成人精品婷婷| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久av| 日韩一区二区视频免费看| 精品一区在线观看国产| 亚洲不卡免费看| 日本vs欧美在线观看视频 | 在线观看人妻少妇| 国产男女内射视频| 日日摸夜夜添夜夜添av毛片| 校园人妻丝袜中文字幕| 国产精品一区二区在线观看99| 精品一区二区三区视频在线| 亚洲天堂av无毛| 国产午夜精品一二区理论片| 18禁裸乳无遮挡免费网站照片| 美女主播在线视频| 日韩在线高清观看一区二区三区| 国产精品嫩草影院av在线观看| 亚洲怡红院男人天堂| 久久精品国产亚洲av涩爱| a 毛片基地|