• <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热这里只有是精品50| 亚洲狠狠婷婷综合久久图片| 9191精品国产免费久久| 97超级碰碰碰精品色视频在线观看| 嫁个100分男人电影在线观看| 国产精品九九99| 巨乳人妻的诱惑在线观看| 国产午夜福利久久久久久| 夜夜爽天天搞| 精品久久久久久久久久久久久| 日韩成人在线观看一区二区三区| 精品欧美一区二区三区在线| 国语自产精品视频在线第100页| 色播亚洲综合网| 欧美乱色亚洲激情| 国产成人影院久久av| 国产精品久久久久久亚洲av鲁大| 国产av一区在线观看免费| 长腿黑丝高跟| av免费在线观看网站| 99热6这里只有精品| 他把我摸到了高潮在线观看| 俺也久久电影网| 国产伦在线观看视频一区| 久久天堂一区二区三区四区| 淫妇啪啪啪对白视频| 18禁美女被吸乳视频| 日韩免费av在线播放| 真人做人爱边吃奶动态| xxxwww97欧美| 制服诱惑二区| 亚洲欧美一区二区三区黑人| 视频区欧美日本亚洲| 国产欧美日韩一区二区三| 国产精品香港三级国产av潘金莲| 曰老女人黄片| 麻豆av在线久日| 欧美黑人精品巨大| 欧美日韩瑟瑟在线播放| 久久久国产成人精品二区| cao死你这个sao货| 午夜免费观看网址| ponron亚洲| 麻豆一二三区av精品| 国产成人啪精品午夜网站| 久久午夜亚洲精品久久| 一二三四在线观看免费中文在| 久久精品国产99精品国产亚洲性色| 丰满人妻一区二区三区视频av | 熟女少妇亚洲综合色aaa.| 女警被强在线播放| 国产精品美女特级片免费视频播放器 | 黄色视频不卡| 亚洲 欧美 日韩 在线 免费| 人人妻人人看人人澡| 人妻丰满熟妇av一区二区三区| 老司机午夜福利在线观看视频| 久久久久久久精品吃奶| 成人午夜高清在线视频| 久久中文字幕人妻熟女| 久久中文字幕一级| 欧美日韩乱码在线| 18禁美女被吸乳视频| 免费电影在线观看免费观看| 亚洲av五月六月丁香网| 又爽又黄无遮挡网站| 久久精品亚洲精品国产色婷小说| 嫩草影院精品99| 九色国产91popny在线| 午夜影院日韩av| 两个人视频免费观看高清| 欧美最黄视频在线播放免费| 亚洲精品色激情综合| 搞女人的毛片| 女人被狂操c到高潮| 国产午夜精品久久久久久| 麻豆久久精品国产亚洲av| 国产精品亚洲一级av第二区| 久热爱精品视频在线9| 一级黄色大片毛片| svipshipincom国产片| 一级作爱视频免费观看| 19禁男女啪啪无遮挡网站| 身体一侧抽搐| 亚洲国产日韩欧美精品在线观看 | 巨乳人妻的诱惑在线观看| 日本黄色视频三级网站网址| 在线国产一区二区在线| 国产又黄又爽又无遮挡在线| 女人被狂操c到高潮| 99久久国产精品久久久| 免费看十八禁软件| 日本免费a在线| 久久精品综合一区二区三区| 国产成年人精品一区二区| 88av欧美| 国产91精品成人一区二区三区| 欧美日本亚洲视频在线播放| 精品一区二区三区视频在线观看免费| 亚洲国产高清在线一区二区三| 999精品在线视频| 小说图片视频综合网站| 亚洲成人免费电影在线观看| 在线观看66精品国产| 亚洲国产精品久久男人天堂| 热99re8久久精品国产| 麻豆av在线久日| 丝袜人妻中文字幕| 老鸭窝网址在线观看| 在线免费观看的www视频| 亚洲男人天堂网一区| 国产成年人精品一区二区| 久久久水蜜桃国产精品网| 午夜福利免费观看在线| 免费看美女性在线毛片视频| 国产伦人伦偷精品视频| cao死你这个sao货| 久久精品91无色码中文字幕| 亚洲性夜色夜夜综合| 国产三级在线视频| 99国产精品99久久久久| 久久久精品国产亚洲av高清涩受| 老鸭窝网址在线观看| av有码第一页| 欧美黄色片欧美黄色片| 免费看a级黄色片| 国产精品亚洲美女久久久| 久久久久国产精品人妻aⅴ院| 欧美大码av| 一区二区三区国产精品乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 国产亚洲av嫩草精品影院| 免费av毛片视频| 美女免费视频网站| 在线国产一区二区在线| 国产精品一区二区免费欧美| 欧美丝袜亚洲另类 | 久久久久国产精品人妻aⅴ院| 久久久久免费精品人妻一区二区| 国产激情偷乱视频一区二区| 成人欧美大片| 精品国内亚洲2022精品成人| 正在播放国产对白刺激| 非洲黑人性xxxx精品又粗又长| 午夜老司机福利片| 亚洲成av人片免费观看| 高清毛片免费观看视频网站| cao死你这个sao货| 男插女下体视频免费在线播放| 久久这里只有精品19| 国产精品久久电影中文字幕| 午夜精品久久久久久毛片777| 日本撒尿小便嘘嘘汇集6| 岛国视频午夜一区免费看| 亚洲乱码一区二区免费版| 人妻夜夜爽99麻豆av| 亚洲一区中文字幕在线| 亚洲精品在线美女| 亚洲精品一区av在线观看| 精品少妇一区二区三区视频日本电影| 日韩中文字幕欧美一区二区| 日韩免费av在线播放| 美女午夜性视频免费| 91麻豆精品激情在线观看国产| 国产一区二区在线av高清观看| 国产黄a三级三级三级人| 精品免费久久久久久久清纯| 国产成人影院久久av| 午夜激情福利司机影院| 亚洲无线在线观看| 国内毛片毛片毛片毛片毛片| 久久天堂一区二区三区四区| 人人妻人人澡欧美一区二区| 伊人久久大香线蕉亚洲五| 三级男女做爰猛烈吃奶摸视频| 国产精品98久久久久久宅男小说| 熟女电影av网| 女人高潮潮喷娇喘18禁视频| 国产精品,欧美在线| 俄罗斯特黄特色一大片| 一个人观看的视频www高清免费观看 | 国产亚洲av高清不卡| 欧美+亚洲+日韩+国产| 亚洲,欧美精品.| 亚洲av美国av| 一个人免费在线观看电影 | 十八禁网站免费在线| 最新美女视频免费是黄的| 欧美色欧美亚洲另类二区| 国产成人av激情在线播放| 亚洲成人免费电影在线观看| 12—13女人毛片做爰片一| 人人妻人人澡欧美一区二区| 国产免费av片在线观看野外av| 黄色视频,在线免费观看| 性色av乱码一区二区三区2| 成人手机av| 日韩大码丰满熟妇| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人精品无人区| 午夜福利视频1000在线观看| 欧美色视频一区免费| 欧美成狂野欧美在线观看| 亚洲七黄色美女视频| 欧美成人一区二区免费高清观看 | 欧美成人午夜精品| 老熟妇仑乱视频hdxx| 男女下面进入的视频免费午夜| 18禁观看日本| 每晚都被弄得嗷嗷叫到高潮| 亚洲av成人不卡在线观看播放网| 日韩欧美国产一区二区入口| 97人妻精品一区二区三区麻豆| 国产亚洲av嫩草精品影院| 在线观看免费午夜福利视频| 非洲黑人性xxxx精品又粗又长| 久久草成人影院| 高清在线国产一区| 亚洲人成网站高清观看| 在线国产一区二区在线| 天天一区二区日本电影三级| 亚洲欧美精品综合一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 99热6这里只有精品| 欧美丝袜亚洲另类 | 全区人妻精品视频| 人妻夜夜爽99麻豆av| 无遮挡黄片免费观看| 人人妻,人人澡人人爽秒播| 午夜福利18| 日本精品一区二区三区蜜桃| 亚洲七黄色美女视频| bbb黄色大片| 午夜久久久久精精品| 无限看片的www在线观看| 18禁观看日本| 日韩 欧美 亚洲 中文字幕| 欧美黑人精品巨大| 国产亚洲欧美在线一区二区| 久久久久国内视频| 别揉我奶头~嗯~啊~动态视频| 又大又爽又粗| 午夜福利视频1000在线观看| 精品熟女少妇八av免费久了| 免费人成视频x8x8入口观看| 欧美黑人巨大hd| 久久久国产精品麻豆| 精品久久久久久久毛片微露脸| 免费搜索国产男女视频| 久久久国产成人免费| 成人国语在线视频| 久久久久精品国产欧美久久久| 午夜精品一区二区三区免费看| 国产精品一区二区精品视频观看| 亚洲欧美日韩无卡精品| 亚洲美女视频黄频| av天堂在线播放| 国产免费av片在线观看野外av| 19禁男女啪啪无遮挡网站| 一个人免费在线观看电影 | 88av欧美| 香蕉丝袜av| 18禁裸乳无遮挡免费网站照片| 操出白浆在线播放| 免费看日本二区| 久久久久久久久中文| 777久久人妻少妇嫩草av网站| 人妻夜夜爽99麻豆av| 国产野战对白在线观看| 国产私拍福利视频在线观看| 舔av片在线| a级毛片在线看网站| 女警被强在线播放| 精品久久久久久成人av| 窝窝影院91人妻| 欧美不卡视频在线免费观看 | 日韩av在线大香蕉| 黄片大片在线免费观看| www国产在线视频色| 久久久精品大字幕| 精品福利观看| 亚洲精品在线观看二区| 亚洲国产欧洲综合997久久,| 亚洲黑人精品在线| 国产成人欧美在线观看| 一区二区三区国产精品乱码| 国产精品综合久久久久久久免费| 国产精品,欧美在线| 国产精品久久久久久人妻精品电影| 国产视频内射| 日本一区二区免费在线视频| 中出人妻视频一区二区| 久久 成人 亚洲| 一进一出抽搐gif免费好疼| 国产精品久久久久久精品电影| 午夜免费观看网址| 中文字幕熟女人妻在线| 男插女下体视频免费在线播放| 中文字幕精品亚洲无线码一区| 欧美大码av| 欧美日韩精品网址| 亚洲成人中文字幕在线播放| 国产成年人精品一区二区| 国产精品九九99| 99riav亚洲国产免费| 狂野欧美白嫩少妇大欣赏| 欧美成人免费av一区二区三区| 欧美乱色亚洲激情| 天天躁狠狠躁夜夜躁狠狠躁| av在线天堂中文字幕| a级毛片在线看网站| 亚洲免费av在线视频| 一级毛片女人18水好多| 亚洲一区中文字幕在线| 日韩欧美国产一区二区入口| 国产精品久久久人人做人人爽| 琪琪午夜伦伦电影理论片6080| 99riav亚洲国产免费| 岛国在线观看网站| 国产亚洲精品第一综合不卡| 国产精品一区二区三区四区久久| 日本免费a在线| 国内精品久久久久久久电影| 色综合婷婷激情| 男人舔奶头视频| 国产欧美日韩精品亚洲av| 国产亚洲av嫩草精品影院| 两个人看的免费小视频| 国产99白浆流出| av福利片在线观看| 小说图片视频综合网站| 国产精品香港三级国产av潘金莲| 亚洲av第一区精品v没综合| 欧美成人免费av一区二区三区| 成人永久免费在线观看视频| 国产男靠女视频免费网站| 免费在线观看完整版高清| 亚洲色图 男人天堂 中文字幕| 国产激情久久老熟女| 久热爱精品视频在线9| 午夜a级毛片| 国产黄色小视频在线观看| 国产精品美女特级片免费视频播放器 | 国产三级中文精品| 国产精品一区二区三区四区久久| tocl精华| 可以在线观看毛片的网站| 999久久久精品免费观看国产| av欧美777| 一夜夜www| 美女 人体艺术 gogo| 可以在线观看毛片的网站| 久久久久久免费高清国产稀缺| 精品国产美女av久久久久小说| 精品免费久久久久久久清纯| 亚洲欧美激情综合另类| 欧美最黄视频在线播放免费| 丁香六月欧美| 女同久久另类99精品国产91| 国产精品av视频在线免费观看| 亚洲av熟女| 欧美成人免费av一区二区三区| 久久精品亚洲精品国产色婷小说| 久久欧美精品欧美久久欧美| 国产精品一区二区三区四区久久| 嫁个100分男人电影在线观看| 国产片内射在线| 三级男女做爰猛烈吃奶摸视频| 日本 av在线| 大型av网站在线播放| 成人亚洲精品av一区二区| 一本大道久久a久久精品| 国产精品影院久久| 中文字幕熟女人妻在线| x7x7x7水蜜桃| 一区二区三区国产精品乱码| 一进一出抽搐gif免费好疼| 国产精品精品国产色婷婷| 国产欧美日韩一区二区精品| 99热这里只有是精品50| 黄色 视频免费看| 叶爱在线成人免费视频播放| 法律面前人人平等表现在哪些方面| 在线a可以看的网站| 亚洲男人天堂网一区| 欧美色视频一区免费| 看黄色毛片网站| 国产真人三级小视频在线观看| 国产高清有码在线观看视频 | 国产精品一区二区免费欧美| 国产一区二区三区视频了| 19禁男女啪啪无遮挡网站| av片东京热男人的天堂| 少妇熟女aⅴ在线视频| 久久精品aⅴ一区二区三区四区| 18禁裸乳无遮挡免费网站照片| 欧美大码av| 亚洲成人久久性| 天天添夜夜摸| 老司机午夜福利在线观看视频| 免费无遮挡裸体视频| 两个人免费观看高清视频| 成年人黄色毛片网站| 中文字幕高清在线视频| 日本五十路高清| 一个人观看的视频www高清免费观看 | 老司机靠b影院| av有码第一页| 久久精品国产99精品国产亚洲性色| 国产精品国产高清国产av| 亚洲国产欧洲综合997久久,| 午夜精品久久久久久毛片777| 男男h啪啪无遮挡| 亚洲 国产 在线| 国内毛片毛片毛片毛片毛片| 人妻久久中文字幕网| 一二三四在线观看免费中文在| 美女黄网站色视频| 麻豆国产av国片精品| 日日夜夜操网爽| 亚洲精品一区av在线观看| 亚洲精品中文字幕一二三四区| 法律面前人人平等表现在哪些方面| 国产精品亚洲一级av第二区| 亚洲av成人av| 久久久久性生活片| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美黑人欧美精品刺激| 成人av在线播放网站| 欧美乱码精品一区二区三区| 一级a爱片免费观看的视频| 亚洲人与动物交配视频| 亚洲欧美一区二区三区黑人| 99热这里只有精品一区 | 成人18禁在线播放| 99久久久亚洲精品蜜臀av| 香蕉丝袜av| 久久国产精品人妻蜜桃| 欧美日本亚洲视频在线播放| 中国美女看黄片| 熟妇人妻久久中文字幕3abv| 亚洲激情在线av| 亚洲人与动物交配视频| 婷婷亚洲欧美| 亚洲中文av在线| 亚洲avbb在线观看| 日韩精品青青久久久久久| 中亚洲国语对白在线视频| 看黄色毛片网站| 日韩 欧美 亚洲 中文字幕| 成年人黄色毛片网站| 女人爽到高潮嗷嗷叫在线视频| 18禁国产床啪视频网站| 很黄的视频免费| 777久久人妻少妇嫩草av网站| 国产一区二区在线av高清观看| 中出人妻视频一区二区| 老汉色∧v一级毛片| 九九热线精品视视频播放| 午夜成年电影在线免费观看| 少妇粗大呻吟视频| 老鸭窝网址在线观看| 看免费av毛片| 村上凉子中文字幕在线| 禁无遮挡网站| 久久久久国产精品人妻aⅴ院| 天堂√8在线中文| 精品久久久久久成人av| 一进一出抽搐gif免费好疼| 在线观看美女被高潮喷水网站 | 欧美日韩亚洲综合一区二区三区_| 国产精品 欧美亚洲| 午夜视频精品福利| 99久久无色码亚洲精品果冻| 国产1区2区3区精品| 国产av一区二区精品久久| 久久久久性生活片| 999久久久精品免费观看国产| 亚洲精品中文字幕一二三四区| 国产99久久九九免费精品| 国产精品香港三级国产av潘金莲| 亚洲人成电影免费在线| 麻豆成人av在线观看| 我要搜黄色片| 99国产精品99久久久久| 1024香蕉在线观看| 久久久久性生活片| 麻豆久久精品国产亚洲av| 老鸭窝网址在线观看| 在线观看舔阴道视频| 99热只有精品国产| 欧美日韩亚洲国产一区二区在线观看| 老熟妇仑乱视频hdxx| 日韩av在线大香蕉| 波多野结衣巨乳人妻| 国产免费av片在线观看野外av| 欧美黄色淫秽网站| 久久久久久九九精品二区国产 | 成人精品一区二区免费| 国产av一区二区精品久久| 亚洲免费av在线视频| 一区二区三区高清视频在线| 成人午夜高清在线视频| 日韩欧美在线乱码| 两个人免费观看高清视频| 欧美黄色淫秽网站| 999精品在线视频| 一个人观看的视频www高清免费观看 | 久久99热这里只有精品18| 99久久精品国产亚洲精品| 身体一侧抽搐| 在线观看免费日韩欧美大片| 小说图片视频综合网站| 国产成人精品无人区| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 操出白浆在线播放| 女同久久另类99精品国产91| 在线观看美女被高潮喷水网站 | 国产精品 国内视频| 在线观看舔阴道视频| 波多野结衣高清作品| 老司机福利观看| 波多野结衣高清无吗| 给我免费播放毛片高清在线观看| 又黄又爽又免费观看的视频| 露出奶头的视频| 午夜精品一区二区三区免费看| 日本 av在线| 久久久久久亚洲精品国产蜜桃av| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久久末码| 好看av亚洲va欧美ⅴa在| 九色成人免费人妻av| 淫妇啪啪啪对白视频| 宅男免费午夜| 亚洲精品国产一区二区精华液| 国产野战对白在线观看| 亚洲人与动物交配视频| 日韩欧美 国产精品| 国产亚洲av嫩草精品影院| av福利片在线观看| cao死你这个sao货| 国产精品乱码一区二三区的特点| 久久久久久大精品| 免费在线观看日本一区| 五月玫瑰六月丁香| 九色国产91popny在线| 在线免费观看的www视频| 久久久久精品国产欧美久久久| 欧美在线黄色| 精品欧美一区二区三区在线| 亚洲自拍偷在线| 亚洲成a人片在线一区二区| 欧美+亚洲+日韩+国产| 99精品久久久久人妻精品| 精品国产超薄肉色丝袜足j| 国产成人影院久久av| 久久国产精品影院| 人成视频在线观看免费观看| 亚洲自偷自拍图片 自拍| 国产精品精品国产色婷婷| 99久久99久久久精品蜜桃| 一个人免费在线观看的高清视频| 亚洲 欧美一区二区三区| 99精品欧美一区二区三区四区| 成人18禁在线播放| 男女下面进入的视频免费午夜| 99久久综合精品五月天人人| 国产成人系列免费观看| 别揉我奶头~嗯~啊~动态视频| 亚洲国产高清在线一区二区三| 亚洲精品久久国产高清桃花| 变态另类成人亚洲欧美熟女| 男女视频在线观看网站免费 | 国产成年人精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 18禁黄网站禁片免费观看直播| 可以免费在线观看a视频的电影网站| 2021天堂中文幕一二区在线观| 亚洲成a人片在线一区二区| av中文乱码字幕在线| 亚洲人成伊人成综合网2020| 一区二区三区高清视频在线| 后天国语完整版免费观看| 欧美日韩黄片免| 欧美乱妇无乱码| 妹子高潮喷水视频| 香蕉国产在线看| 日本黄大片高清| 亚洲在线自拍视频| 久久久久久久久久黄片| 久久国产精品人妻蜜桃| 久久中文看片网| 国产精品一区二区免费欧美| 丰满人妻一区二区三区视频av | 亚洲最大成人中文| 欧美日韩亚洲综合一区二区三区_| 99国产精品99久久久久| 欧美一区二区国产精品久久精品 | 全区人妻精品视频| 亚洲五月婷婷丁香| 日本黄大片高清| 欧美日韩一级在线毛片| 嫩草影视91久久| 久久久久性生活片| 国产精品久久久av美女十八| 一卡2卡三卡四卡精品乱码亚洲| 草草在线视频免费看| 制服人妻中文乱码| 长腿黑丝高跟| 男男h啪啪无遮挡| 制服丝袜大香蕉在线| 欧美性长视频在线观看| 夜夜爽天天搞|