the physics involved ranges from subnanosecond submillimeter electron avalanches governed by atomic and molecular ionization and attachment cross sections to centimeterscale propagating space charge waves ( streamers ) to microsecondscale meterlong channels of ionized gas described by nonequilibrium plasma and electrodynamic behavior . on a larger scale ,
spark behavior in nature ranges from 100  km lightning channels that last seconds and span conditions from the upper reaches of a thundercloud to the microsecondscale processes by which such a channel attaches to the ground . in such a diverse range of behavior one can find many puzzling phenomena .
the process by which a discharge initiates in the cloud is not understood nor is the process by which lightning channels ( leaders ) extend and branch .
the extension process sometimes proceeds as a series of microsecondscale steps separated by 10  s of microseconds of relative quiescence .
more puzzling , the steps seem to be immediately preceded by formation of a hot conductive channel displaced from the end of the existing channel that rapidly connects with the main channel , though the mechanism by which this  space leader 
the transition from streamer to leader is the focus of much research as are the dynamics of the streamer itself and the role of energetic particles in the process .
this role of energetic particles is especially interesting in the light of observations of xray production by lightning and the hypothesis that energetic particles produced by lightning contribute to the production of terrestrial gamma ray flashes .
xrays as produced by lightning seem to occur coincident with stepwise extensions of the channel [ howard et al . , 2010 ;
dwyer et al . , 2011 ] and have energies around a few 100  kev [ moore et al . , 2001 ; dwyer et al . , 2003 ;
these xrays imply the existence of a much larger population of higherenergy electrons . while lowenergy electrons encounter high dynamic friction , relativistic electrons encounter much lower friction and can  run away  to very high energies when driven by electric fields , a fact first recognized by wilson   and modeled in the context of electric fields near lightning channels in carlson [ 2009 , chapters 5 and 6 ] . 
seed  energetic electrons necessary to initiate such a process can come from energetic background radiation ( i.e. , cosmic rays [ carlson et al . , 2008 ] ) , lowenergy electrons in local electric fields strong enough to overcome the maximum friction force [ e.g. , moss et al .
, 2006 ; li et al . , 2009 ; chanrion and neubert , 2010 ] , or by feedback from prior generations of energetic electrons [ dwyer , 2003 , 2007 ] .
regardless , electric fields may drive avalanche growth of populations of such runaway electrons[gurevich et al .
unfortunately , the relative importance of these processes for phenomena like xray production by lightning or terrestrial gamma ray flashes is not well understood . in the context of these puzzles , any additional information about the processes involved
, laboratory studies can shed light on the detailed dynamics of leader extension , and while the energy scales are necessarily smaller ( of order 1  mv compared to the 10  s of megavolts of natural lightning ) , lab sparks even produce xrays .
dwyer [ 2005a ] report the first detection of xrays associated with such sparks in a study of 14 discharges of 1.5  mv , with total energy deposited in a detector up to mev scale .
variation of signals among xray detectors with different attenuators suggest xray energies from 30  kev to 150  kev piling up to produce the mevscale observations .
report similar results with total energy deposited up to several mev with 1  mv discharge .
( elaborated in nguyen ) also report similar results for 0.881  mv sparks with a much larger data set and with a variety of detectors , attenuators , and positions of detector .
nguyen et al .   show intensity variations with highvoltage ( hv ) polarity and with distance between detector and spark gap , suggesting that such intensity variations could result from xray production by positive streamers from metal structures near the detectors .
report a further study of 241 discharges of 1  mv with various polarities and gap distances and push the maximum energy deposited in a single detector up to 50  mev with varied attenuators providing evidence for individual photon energies exceeding 300  kev and statistical evidence for an average photon energy at most 230  kev .
dwyer et al .   also report collimator experiments supporting the production of xrays by a diffuse source in the gap between the electrodes and that bursts of xrays produced late in the discharge come from elsewhere .
march and montany report that hv pulses with rapid risetime tend to produce more xrays , and march and montany   examines the effect of electrode geometry .
kochkin et al . , kochkin et al . , and kochkin et al .  
add nanosecondresolution photography as a useful tool , mapping the streamer clouds produced by positive [ kochkin et al .
, 2012 ] and negative [ kochkin et al . , 2015 ] hv discharges .
in such experiments the hv risetime is much longer than the time for streamers to propagate across the gap , and in experiments with negative high voltage the negative streamers tend to appear and grow in bursts [ see , e.g. , kochkin et al . , 2014 , figure 6 ] .
demonstrate for negative hv that the xrays tend to be emitted during these bursts of streamer development but are emitted on a much shorter timescale than the burst , supporting the suggestion that the large and transient electric fields that result from interaction of negative and positive streamer fronts may play a role in energetic electron and thus xray production .
[ 2012 , 2015 ] also describe that attenuator experiments they claim are consistent with 200  kev characteristic xray energy .
these estimates are based simply on comparison of registration rate ( the fraction of sparks for which xrays are observed ) for various attenuators with the expected attenuation in number of photons incident on the detector , assuming all sparks emit the same number of photons and that pileup does not affect their observations , though they acknowledge that pileup has a significant effect .
together , these studies provide a reasonably complete picture of xray emissions from sparks : free electrons are pushed to overcome friction in the highfield region ahead of a negative streamer ( possibly briefly enhanced by a nearby positive streamer ) , gain at least enough energy to run away in the lower fields surrounding the streamer , then undergo bremsstrahlung in the surrounding air to emit xrays .
these intense streamer fields exist only as the discharge develops and are unlikely to exist during the discharge itself , though the intense electrodynamic environment during discharge may drive discharges elsewhere in the lab that also produce detectable xrays .
average  energy , the properties of the photon spectrum are unknown and difficult to judge due to pileup , and the fact that some sparks produce copious xrays while others produce no measurable xrays is not addressed . the present paper attempts to shed light on these topics by examining the energies measured in xray detectors placed near a 1  m spark gap at the technical university of eindhoven during 921 shots as described in section  2 .
these shots give us the opportunity to statistically analyze the distributions of spark xray fluence , number of photons detected by multiple detectors , and photon energy .
these distributions are described in section  3 and constrained by comparison to observations in section  4 .
the experiment was conducted in the highvoltage lab at the technical university of eindhoven with a haefly 2  mv marx generator configured to produce 1  mv pulses with 1  s risetime .
all tests were carried out with a 1  m pointpoint gap with negative hv polarity .
the spark voltage , ground electrode current , and highvoltage electrode current were recorded by lecroy fourchannel storage oscilloscopes configured to record 2  s of data at 10  ghz sample frequency when triggered by the marx generator , as were signals direct from the output of photomultiplier tubes monitoring a variety of energetic radiation detectors .
the radiation detectors used in the experiment included scintillating plastic optical fibers placed at a variety of locations around the spark to monitor energetic electrons .
analysis of these data has been presented [ ostgaard et al . , 2014 ] and will be published separately . due to the need to run many shots with the scintillating fiber detectors at various locations ,
a total of 950 shots were carried out from which we recorded usable data for 921 . during all shots ,
two additional xray detectors were running , providing a wealth of data about the photon population .
these two detectors are each composed of a 1.5  inch long 1.5  inch diameter labr3(ce ) scintillator monitored by a photomultiplier tube ( pmt ) .
these detectors are placed next to each other ( separation 10  cm ) in an electromagnetic compatibility ( emc ) cabinet to shield the detectors from the electromagnetic noise produced by the spark .
the location of the detectors relative to the spark corresponds to location h as shown in figure  1 of kochkin et al . .
the scintillator material is separated from the spark by roughly 0.5  mm of aluminum in the scintillator housing and emc cabinet wall and approximately 2  m of air .
the signals as recorded by the oscilloscopes are in volts and here have been converted to mev by use of cs137 and co60 gamma ray emitters as sources of known energy .
the lowest energy detectable by this setup is 20  kev , comparable to the minimum photon energies transmitted through the 0.5  mm aluminum shielding , while the maximum depends on the settings of the oscilloscope .
signal from detector 1 clips at just above 5  mev , while detector 2 clips at just above 3.5  mev .
the highvoltage trace shown is very consistent from one spark to the next , as is the highvoltage electrode current .
the pulselike features on the highvoltage electrode current occur when bursts of corona and streamer activity carry charge away from the electrode .
see kochkin et al .   for a detailed discussion of such processes as seen in highspeed camera imagery .
( top left ) the voltage of the highvoltage electrode with respect to the ground electrode .
the pulses seen in the oscilloscope traces indicate deposition of energy in the scintillator , but a single pulse from the pmts may be produced by multiple xray photons entering the scintillator .
these photons may arrive at slightly different times and produce a visibly altered pulse shape , but typically no time structure is visible .
the photon spectrum is therefore not directly measurable , and pulse pileup must be treated statistically . to compile data for such statistics , we search through the data for pulses . to search , we first smooth the data by lowpass filtering
once pulses are identified , we collect a variety of data including pulse time , height , integral , and duration .
pulse integral is a more robust measure of deposited energy in case multiple pulses arrive near simultaneously , so our analysis here reports energy as determined by pulse integral .
while pulse integral is somewhat resistant to saturation when oscilloscope signals clip , the energy of a saturated pulse can not reliably be determined , so we enforce a maximum energy for each detector corresponding to the energy of the largest unclipped pulse . in cases when multiple pulses are observed in a single detector , we add their energies together to better capture variation from one spark to the next , though note that the majority of sparks have zero or only one burst so this combination of bursts only affects a small fraction of our data .
the data set we consider here is therefore a set of pairs of numbers , each pair associated with a single spark and each number with the pair giving the total energy deposited in a detector during the given spark .
a scatterplot showing the energies deposited in each detector in each shot is shown in figure  2 .
most sparks ( 57% ) produce no detectable signals in either detector , indicating less than 20  kev deposited , while approximately 3% of sparks saturate both detectors , indicating at least 35  mev of energy deposited .
each point represents a spark , where e 
 1 was deposited in detector 1 , while e 
 2 was deposited in detector 2 .
the dotted line shows e 
 1 = e 
 2 . at e 
 1 = e 
 2 = 0 , 523 points overlap , and 20 points
the points clearly cluster in the lower left , but a significant number of points appear in the middle and upper regions of the plot . attempting to explain this distribution most simply , one might assume all sparks are identical and that all photon energies are equal and treat the observed distribution as solely due to poisson statistical fluctuations in the number of photons observed . in this case , the high number of points in the lower left ( 57% of events undetected ) implies a poisson mean of 0.562 . with this mean , observation of events that saturated both detectors is then incredibly unlikely as deposition of more than 8  mev ( 3  mev in one and 5  mev in the other ) requires at least 40 photons for the mean energy 200  kev consistent with the results quoted above .
observation of at least 40 photons from a poisson distribution with mean 0.562 has a probability around 10 , much less than the 3% of observed sparks that saturate both detectors . adding the complication of a photon energy distribution is logical , but does not help enough , since the maximum photon energy must be less than 1  mev as the photons will be much less energetic than the 1  mv maximum spark voltage .
assuming the same poisson mean as before and attempting to explain the events that saturated both detectors as due to an extreme fluctuation in photon energy such that all photons carried 1  mev of energy means that we must have at least eight photons instead of at least 40 , but the poisson probability of observing at least eight photons given a mean of 0.562 is roughly 10 , still much smaller than the 3% observed .
clearly , there must be a significant variability in spark xray fluence at the location of the detectors from one spark to the next .
this variability could come from intrinsic variability in spark xray luminosity or it could come from some variability in the geometry if xray emissions are not isotropic . regardless , a full explanation of the distribution of points in figure  2 must include the effect of the distribution of spark xray fluence , the distribution of numbers of photons detected by each detector ( given the xray fluence ) , and the distribution of photon energies . these distributions and their properties are the focus of this paper .
in attempting to build a statistical model of the distribution of points in figure  2 , we need to know the form of the distributions of physical properties relevant to the point locations .
there are three main distributions at work : the distribution of spark xray fluence , the distribution of numbers of photons incident on each detector ( given spark fluence ) , and the distribution of photon energies . the distribution of spark fluence is determined in principle by the processes at work in electron acceleration .
one spark has an intrinsically higher luminosity than another by having more energetic electrons , perhaps because random streamer branching events happened to lead to more negative streamers or perhaps fewer but more intense or more rapidly growing negative streamers or perhaps more interactions between negative streamers and positive streamers .
fluence variability may also result from anisotropy in the directionality of xray emissions : perhaps all sparks emit many xrays , but the emissions are not always beamed toward the detectors .
these processes may be chaotic or perhaps intrinsically random , but either way there will be some variability in xray fluence . to represent this variability , we treat the xray fluence as the expected value of the number of photons hitting a detector ,  , and assume some probability distribution for the occurrence of various values of . as there is yet no theoretical expectation for the form of this distribution , we must make some assumption .
examining the clustering of points in the lower left of figure  2 , we expect lower fluences to be more common than higher fluences , but the existence of points in the upper right suggests that the tail of the distribution to high fluences is quite strong . to represent distributions that follow this general trend of decreasing with a potentially long tail , we assume the distribution for  is a power law with index  : 
 ( 1)p()=(+1)max+1min+1 where  
 min and  
 max are lower and upper limits on  necessary to ensure that the distribution is normalizable for 1 and the leading constant ensures normalization . by examining the clustering in the lower left in figure  2
, we expect  < 0 , but the number of points saturating both detectors suggests that  is not too negative and that the distribution of  is in some sense quite hard .
given the fluence of the spark ( i.e. , given  ) , the number of photons hitting a detector , n , should be poisson distributed with mean . normalized , 
 ( 2)p(n)=nen ! 
 
this assumes that there is a very large number of photons produced by the spark , each with a very small probability of hitting the detector , and that the two detectors have the same probability of catching each photon .
this implicitly assumes that xray emissions from the spark are uniform on the scale of the 10  cm separation of the detectors , but this is reasonable as the detectors are roughly 2  m from the spark and bremsstrahlung from electrons with only a few hundred kev of energy is not strongly beamed .
finally , given a number of xrays hitting each detector , one can calculate the distribution of deposited energy by applying a distribution of xray energies .
this unfortunately is unknown as we do not know the distribution of electron energies responsible for the bremsstrahlung emissions [ chanrion and neubert , 2010 ] .
figure 2 presents calculations of energy distributions of free electrons as produced in constant electric fields that appear linear on a semi  log plot ( inset ) and that show the effects of an exponential cutoff up to the threshold of runaway electron behavior .
the resulting photon distribution will share some of these characteristics , but in principle connecting the photon and electron distributions requires detailed treatment of bremsstrahlung , a topic beyond the scope of the current work .
as such , we simply assume that the photon energy 
 e is also exponentially distributed with mean photon energy  : 
 ( 3)p(e)=ee/ where the division by  ensures normalization for 
 0<e<. note that for simplicity we are assuming that the photon energy distribution does not depend on xray fluence .
it is reasonable to expect that  depends on  somewhat , but including this dependence would require still more assumptions and additional free parameters . assuming an exponential photon energy distribution is especially convenient since the distribution of energies deposited in a detector e , a sum of the energies of n photons ( 
 e=i=1nei ) will be distributed as the sum of n independent exponentially distributed random variables .
the distribution of such a sum of exponential random variables can be expressed in closed form and is known as the erlang distribution : 
 ( 4)pn,(e)=en1ee/n(n1 ) ! 
 
note , however , that if n = 0 , no photons strike the detector , the energy deposited is exactly zero , and p 
 n = 0,(e ) = (e ) , a dirac delta function , instead of the continuous distribution above . if the fluence of the spark (  ) and the number of xrays hitting each detector ( n 
 1 , n 
 2 ) were known , the joint distribution of energy depositions would simply be given by a product of two erlang distributions . unfortunately , the fluence of the spark and the number of xrays hitting each detector are not known . as such
, we must compute the marginal distribution of observed energies by integrating ( i.e. , taking an weighted average ) over the possible values of the unknown quantities , weighted by the probability distributions associated with those quantities .
for one detector , 
 ( 5)p,(e)=minmaxdp()n=0p(n)pn,(e ) and for the joint distribution of simultaneous observations with two detectors , 
 ( 6)p,(e1,e2)=minmaxdp()n1=0n2=0p(n1)pn1,(e1)p(n2)pn2,(e2 ) where the integral over  covers the range of possible spark xray fluence and the sums over n or n 
 1 and n 
 2 cover the possible numbers of photons observed by each detector . the resulting joint probability distribution of e 
 1 and e 
 2 should be able to reproduce the observations in figure  2
there is a finite probability that both detectors observe exactly zero energy , a discrete probability represented by the n 
 1 = n 
 2 = 0 term in the sums for which 
 pn1, and 
 pn2, become delta functions as described above .
likewise , the n 
 1 = 0 , n 
 2>0 , and n 
 2=0 ,
n 
 1>0 terms in the probability distribution contain a single delta function . only the terms for which both n 
 1 and n 
 2 are greater than zero lack delta function divergence behavior .
these delta functions make the probability distribution difficult to work with directly but are capable of representing the clustering of points at the origin and along the axes in figure  2 .
these probability distributions form the basis of our model , and we seek to infer the parameters  ,  ,  
 min , and  
 max by comparison of our model to data .
comparison of a probability distribution to data can be done in many ways , for example , the kolmogorovsmirnov test or the andersondarling test .
such tests work on the basis of the cumulative distribution function ( cdf ) f , which for a probability distribution p(x ) in one dimension is defined as 
 f(x)=xp()d. since here we work in two dimensions ( e 
 1,e 
 2 ) , we need a twodimensional analog of the cumulative distribution function . here
we follow fasano and franceschini   in taking our set of data points { p 
 i } = { ( e 
 1i , e 
 2i ) } and for each point , assigning four cdflike values : 
 fdati < < , 
 fdati< , 
 fdati < , and 
 fdati , each giving the fraction of the observed data set in the corresponding quadrant relative to the data point in question .
for example , 
 fdati < < gives the fraction of data { p 
 j } for which e 
 1j < e 
 1i and e 
 2j < e 
 2i .
these 
 { fdati  } can be plotted if an order is assigned to the data points for use as an abscissa .
here it is convenient to order each set of 
 fdati  separately such that 
 fdati  is monotonically increasing , i.e. , sort each set of 
 fdati  from smallest to largest and use its place in the resulting list as its abscissa .
for example , the grey , curve starts at 0.03 since 3% of sparks saturate both detectors : the data point that has the fewest data points with energies greater than or equal to its energy must be one of those doubly saturated events , and thus , 3% of the data satisfies the , condition . moving to the right along the , curve , larger cdf values correspond to data points that are below and/or to the left of the upper right corner , i.e. , data points for which larger fractions of the data satisfy the , condition .
the large jump from about 0.35 to 1.0 at cdf rank 400 occurs when incrementally increasing the fraction of the data set that satisfies the , condition requires looking at a point where e 
 1 = e 
 2 = 0 , which thus includes all other such points ( 57% of all data ) as well as points for which e 
 1 = 0 or e 
 2 = 0 .
the cdf therefore jumps up to 1 since all of the data has e 
 10 and e 
 20 .
the direction of the inequality is related to the shape of the distribution , with  and < < related to the lower left / upper right balance and width while the <  and  < versions related to upper left / lower right balance and width .
one can imagine other ways of constructing such cdfs , but this is a straightforward technique that as we will see shortly is quite effective .
each horizontal point represents a point in the ( 2d ) data set , for which the vertical coordinate represents the fraction of the data set with energies related to the energies of the given data point by the inequalities associated with the given curve ( i.e. , the 
 fdati  described in the text ) .
each curve represents a different inequality , and the data have been reordered for each curve to make the curve in question monotonically increasing .
two sets of curves are shown , thick curves from the data and thin curves predicted given the best fit parameters discussed in the text , but the curves overlap so much any deviations are difficult to judge .
such curves as in figure  3 can be predicted on the basis of the probability density ( equation  ( 6 ) ) , assigning a predicted 
 fpredi  to each data point by integrating the expected distributions over the region relevant to the inequalities in question .
for example , 
 ( 7)fpredi<<=e1ide1e2ide2p,(e1,e2 ) 
  similar calculations can be made for predicted 
 fpredi < , 
 fpredi< , 
 fpredi by changing the limits on the integrals to go from the point in question to 
 +  as appropriate .
regardless , the exponentials and powers in the integrand can be manipulated to convert the integrals here and in equation  ( 6 ) into upper and lower incomplete gamma function evaluations ( e.g. , 
 (k , x)=xxk1exdx ) , and the summations in equation  ( 6 ) can be computed numerically and truncated without significant loss of accuracy . for a given  ,  ,  
 min , and  
 max ,
the result of this exercise is four sets of numbers ( 
 { fpredi < < } , 
 { fpredi < } , 
 { fpredi< } , and 
 { fpredi } ) that can be compared to the corresponding sets associated with the data ( 
 { fdati < < } , 
 { fdati < } , 
 { fdati< } , and 
 { fdati } ) . treating these numbers as analogous to the cumulative distribution function
, the kolmogorovsmirnov test statistic would be the maximum deviation between any two corresponding 
 f  i  .
the kolmogorovsmirnov test is sometimes criticized as not very powerful [ e.g. , razali and wah , 2011 ] , so instead , we calculate an analog of the andersondarling test statistic : 
 ( 8)s=i(fdati<<fpredi<<)2fpredi<<(1fpredi<<)+(fdati<fpredi<)2fpredi<(1fpredi<)+(fdati<fpredi<)2fpredi<(1fpredi<)+(fdatifpredi)2fpredi(1fpredi ) i.e. , the squared deviation between the predicted cdf and the observed cdf , weighted as in the andersondarling statistic , summed over all four cdf types ( inequalities ) , and summed over all points in the observed data set .
some numerical problems arise when both numerator and denominator are exactly zero , i.e. , one of the 
 fpredi  is exactly 0 or 1 , as occurs , for example , for 
 fpredi < < for a point , where e 
 1i = e 
 2i = 0
. since such points involve exact match between predicted and observed cdf and therefore should not contribute to the test statistic , we add a small factor ( 0.001 ) to the denominator of each term to ensure the division of zero by zero results in zero .
this does not significantly affect the other terms in the test statistic or the overall results .
we then fit our predicted distribution to the observed distribution  minimizing s  by varying the parameters of our calculated cdf with the downhill simplex algorithm ( nelder and mead   as implemented in johnson ) .
the optimal values are  =  1.29 ,  = 86  kev ,  
 min = 0.022 , and  
 max = 113 . since our multidimensional cdf and modified andersondarling test statistic are so far removed from their original application , we only use the statistic in judging the quality of the fit in the optimization process described above and make no attempt to apply the distributions typically associated with the andersondarling test statistic to assess the uncertainty in our fit results .
instead , we judge uncertainty by bootstrap , repeating our fit process many times with alternative data sets produced from our original data by sampling with replacement .
the fit results are each approximately normally distributed , with mean and standard deviation given as follows :  =  1.29  0.04 ,  = 86  7  kev ,  
 min = 0.022  0.006 , and  
 max = 110  25 .
plotting these bootstrapped fit results , one parameter versus another , shows very little correlation between results , with the exception of  and  
 min which are negatively correlated : more negative  is associated with higher  
 min .
this makes sense given that more negative  is associated with more events with low xray fluence , so raising the minimum fluence is necessary to retain the balance between events with and without detectable signal .
the cdfs computed from the best fit parameters are also shown in figure  3 as thin curves , but they overlap with the data so much that no systematic deviations are evident . a residuals plot showing the differences between the two sets of curves is shown in figure  4 .
some deviations can be seen but are well within the typical size of the fluctuations due to random distribution of points in the data set .
a positive deviation indicates a larger data cdf than computed from the best fit parameters .
monte carlo simulations of data sets drawn from the distributions described above show no obvious deviations from the distribution of the data when plotted as in figure  2 .
it is also worth noting that since our fitting procedure captures the joint distribution , it also captures the distribution within each detector separately . in this paper , we have not attempted to compensate for saturation to construct a true energy distribution , but kochkin et al .
, using essentially the same experimental setup as used here , do present such a spectrum [ kochkin et al . , 2015 , figure 12 ] constructed by a sophisticated procedure of
fitting pulse shapes to observed oscilloscope records . though the analysis in kochkin et al .  
includes the data set used here , the analysis procedures are completely different , and kochkin et al .  
results are presented for energy deposition per burst not per spark , but the fact that relatively few sparks have multiple bursts means that it is still useful to directly compare the results from kochkin et al .   to equation  ( 5 ) .
evaluation of equation  ( 5 ) with our fit result values of  ,  ,  
 min , and  
 max , normalized and overlaid with the data from kochkin et al .
, is shown in figure  5 , demonstrating much better agreement than the single exponential distribution used in kochkin et al . .
inferred a 200  kev characteristic burst energy ( which is asserted to be roughly equal to characteristic photon energy ) but showed an extremely poor fit at high burst energies that they attribute to pileup .
properly accounting for pileup as we do here with the poisson distribution of observed counts ( p 
 (n ) ) and for fluence variability with the assumption of a power law ( p 
 ( ) ) gives us a much better fit with a much lower mean photon energy ( 86  kev ) , suggesting that pileup and fluence variation is essential to consider in interpretation of such data .
comparison of the predicted total energy deposition spectrum from this paper with observations from figure 12 of kochkin et al . .
the red squares show the observed spectrum , while the black curves show the predictions based on our distributions .
the thick curve is the prediction , while the thin curves represent the 
 n standard deviation expected for the spread in the data . as a final comparison to data ,
while we ran no attenuator experiments in the data set described above , kochkin et al .  
they report results for lead attenuators of varying thickness , determining the fraction of sparks that produce observable signals in a single detector by running 50 sparks with each attenuator .
the distributions determined above can easily be used to construct results for attenuated scenarios for comparison with the observations in kochkin et al .   by monte carlo simulation .
first , we draw a random spark fluence , then draw a random number of photons from such a spark , then draw that many random photon energies from the photon energy spectrum , then attenuate those photons probabilistically by use of the mass attenuation coefficients from hubbell and seltzer , and then finally classify the event as containing detected xrays or not based on whether or not any photons made it through the lead shield .
we repeat this procedure 10 times , compute the overall fraction of events with detected xrays , and repeat for each attenuator thickness .
comparison with data requires an estimate of the uncertainty in the data , here estimated by calculating a 68% confidence interval based on the binomial distribution .
sixtyeight percent confidence corresponds to a 1 error bar , and all but one of our predictions are consistent with the observations at this confidence level .
the distributions determined in this paper are thus also consistent with attenuator experiments , though further attenuator observations would be useful to narrow the uncertainties in the observations .
predictions based on the distributions determined in this papera 
  the 68% confidence intervals correspond to 1 error bars and are calculated based on binomial statistics .
to summarize , we have studied observations of xray emission by sparks and attempted to understand those observations by modeling the process as a combination of a distribution of spark fluence , poisson statistics , and an xray energy spectrum .
these distributions successfully reproduce not only the data used here but also independent analysis of similar experiments on overall energy distribution and signal attenuation in kochkin et al . .
the two main results are a welldefined photon mean energy and the distribution of overall spark fluence .
the 86  7  kev mean photon energy determined above is consistent with most earlier estimates but somewhat lower than some .
we feel this mean energy is quite reasonable given that electrons accelerated in the environment of a 1  mv spark realistically must have at most a few hundred kev of energy and bremsstrahlung photons produced must have a lower energy than that of the source electrons . as a crude estimate of electron energy ,
assume bremsstrahlung produced by the relatively low energy electrons here follows approximately the dn / de 
   1/e 
  distribution seen for bremsstrahlung produced by highenergy electrons .
this 1/e 
  distribution spans the range from the energy of the electron as a maximum to a minimum energy at most 20  kev . using 20  kev as the minimum energy and requiring that the average photon energy is 86  kev require the electron energy to be roughly 200  kev , while using 5  kev as the minimum energy gives nearly 400  kev as the electron energy .
this 200400  kev energy range is not inconsistent with the potential available for electron acceleration at the time of xray production as seen in figure  1 , but this is a crude estimate in need of refinement by modeling of electron energy distributions and the resulting bremsstrahlung before it can be taken very seriously .
the exponential xray spectrum motivated by the expectation of an at least approximately exponential distribution in electron energies very closely reproduces the observed distributions , though is likely not the only such distribution to do so , and as noted above we have assumed that fluence (  ) and mean energy are independent which may or may not be true . while the average energy should be close to correct , the other distributions at work in this system complicate the analysis
. the more important factor at work is the very broad distribution of observed spark xray fluence (  ) , which , as mentioned earlier , we assume captures both intrinsic variability in spark luminosity and variability in the spatial distribution of xray emissions .
the results obtained above , which the distribution of  is very hard ( power law index 1.29 ) and covers a broad range , from 0.02 photons to over 100 photons for the experiment conducted here , pose a challenge to the idea that the strong field near the head of a negative streamer is all that is required to accelerate electrons into the energy regime needed for 100  kev xray production .
all sparks involve extensive streamer production [ see , e.g. , kochkin et al . , 2014 , figure 2 ] , so the number or distribution of streamer production can not directly explain the broad distribution of spark xray fluence . as discussed in kochkin et al .
, the very short duration and appearance of multiple bursts of xray emissions suggests that xray emission occurs during some fast transient process like streamer collision , but it is not clear how streamer collision frequency would be distributed and whether such a distribution could reproduce the characteristics observed here .
while we have not proved that xray fluence truly follows a power law , only that our assumption of a power law is consistent with our data , it is perhaps not unreasonable that a power law distribution might arise here in the context of dielectric breakdown since power laws arise in studies of systems that similarly approach breakdown and then collapse [ bak et al . ,
regardless of the true form of the distribution , we hope our analysis here is useful for evaluation of xray production mechanisms , which must explain the frequency of both dim and bright sparks in roughly the balance described here . as for the relative contributions of spark luminosity variation and nonuniform xray spatial distribution to the observed fluence variability
the observed fluence variability should correspond roughly to the degree of spatial variability in xray emissions .
one can set an upper limit on the spatial variability present in xray emissions in an experiment such as this by examining the directional distribution of bremsstrahlung from a unidirectional beam of 500  kev electrons .
[ 1979 , figure  5 ] and khn and ebert [ 2014 , figure  7b ] .
both figures present plots of the intensity of emissions versus angle relative to the direction of the electron beam , with tseng et al .   showing a factor of around 100 difference from highest to lowest fluence , while khn and ebert   show a factor of around 300 . that the best fit distribution constructed here requires a factor of  
 max/ 
 min = 100/0.02
= 5000 from highest to lowest fluence suggests that even this extremely aggressive assumption of unidirectional 500  kev electrons falls short of explaining the observed fluence variability , but this analysis is admittedly crude .
a slightly less crude assessment can be made by monte carlo : first , draw a random beam direction uniformly from the highvoltage electrode in a hemisphere facing the ground electrode , then determine the angle from the center of that beam to the location of the detector and evaluate the relative fluence by extracting the distributions shown in tseng et al .   and khn and ebert   for a relatively high energy photon .
though this ignores the xray energy dependence of the bremsstrahlung angular distribution , it should reasonably account for the degree of variability . repeat this procedure many times , compiling a histogram estimating the fluence distribution , then scale these histograms such that they are normalized and adjust the relative fluence to correspond to  by shifting the relative fluence until the histograms correctly predict the fraction of events that would fall at intensities too low to be detected .
histograms constructed in this manner are compared to the best fit  power law in figure  6 .
predicted fluence distributions assuming only spatial variability due to beaming of bremsstrahlung from unidirectional 500  kev electrons .
the curve marked  tseng  is derived from figure  5 of tseng et al . , 
the similarity of slopes is striking , especially for the khn and ebert   case ; but it is clear that fluence variability due to beaming can not fully account for the observations . in the context of the 86  kev average photon energy determined here , the maximum fluence of the bremsstrahlungderived cases ,   10 , would result in 800  kev deposited on average , far from the 35  mev needed to saturate the detectors .
put another way , if we attempt to account for the large fraction of lowfluence events simply as emissions beamed away from the detectors , the maximum fluence would then be too low to account for events that saturate the detectors .
the scale of variability in the distributions derived from bremsstrahlung , even for the extreme case of unidirectional 500  kev electrons , is therefore too low to fully account for the variability observed , providing clear evidence for strong sparktospark variability . crudely comparing the spread of the khn and ebert   curve to our power law in figure  6 ,
there is at least another order of magnitude variability in fluence that can not be explained as spatial variability and that is a very conservative estimate since in a more realistic scenario there will be directional dispersion of electrons , lower energy electrons , and possibly multiple beams , all of which lead to reductions in the spatial variability due to beaming alone , thus requiring even more sparktospark variation .
we hope the results described above will provide a basis for comparison to results of other spark experiments to determine if the distributions of mean photon energy and spark xray fluence depend on the experimental setup .
we further hope that theoretical work can find a basis for such distributions to confirm or refute our assumptions and that such theories can shed light on the process of runaway electron acceleration in sparks as suggested to be relevant to other electrical phenomena such as lightning or terrestrial gamma ray flashes .
 abstractxray emission by sparks implies bremsstrahlung from a population of energetic electrons , but the details of this process remain a mystery . 
 we present detailed statistical analysis of xray spectra detected by multiple detectors during sparks produced by 1  mv negative highvoltage pulses with 1  s risetime . with over 900 shots , 
 we statistically analyze the signals , assuming that the distribution of spark xray fluence behaves as a power law and that the energy spectrum of xrays detectable after traversing 2  m of air and a thin aluminum shield is exponential . 
 we then determine the parameters of those distributions by fitting cumulative distribution functions to the observations . 
 the fit results match the observations very well if the mean of the exponential xray energy distribution is 86  7  kev and the spark xray fluence power law distribution has index 1.29  0.04 and spans at least 3 orders of magnitude in fluence . 
Introduction
Experimental Setup
Statistical Model
Comparison to Data
Discussion