Scientific Research

An Academic Publisher

Method to Measure Indoor Radon Concentration in an Open Volume with Geiger-Mueller Counters: Analysis from First Principles

**Author(s)**Leave a comment

^{218}Po,

^{214}Po) and the conversion of the alpha count rate into a radon concentration. A unique feature of this method, in comparison with standard methodologies to measure radon concentration, is the absence of a fixed control volume. Advantages afforded by the reported GM method include: 1) it provides a direct in-situ value of radon level, thereby eliminating the need to send samples to an external testing laboratory; 2) it can be applied to monitoring radon levels exhibiting wide short-term variability; 3) it can yield short-term measurements of comparable accuracy and equivalent or higher precision than a commercial radon monitor sampling by passive diffusion; 4) it yields long-term measurements statistically equivalent to commercial radon monitors; 5) it uses the most commonly employed, overall least expensive, and most easily operated type of nuclear instrumentation. As such, the method is par-ticularly suitable for use by researchers, public health personnel, and home dwellers who prefer to monitor indoor radon levels themselves. The results of a consecutive 30-day sequence of 24 hour mean radon measurements by the proposed GM method and a commercial state-of-the-art radon monitor certified for radon testing are compared.

KEYWORDS

1. Introduction

Radon is a colorless, odorless, chemically inert radioactive gas produced in the decay series of uranium-238 (), which leads to radon-222 () with half-life 3.82 days, of thorium-232 (), which leads to “thoron” or radon-220 () with half-life 55.6 seconds, and of uranium-235 (), which leads to “actinon” or radon-219 () with half-life 3.96 seconds [1] . Because comprises a minute fraction (~0.72%) of terrestrial uranium compared with (99.274%), and given the very short half-lives of and compared with, the primary radon isotope of concern in this paper is.

In the decay series leading from to a stable isotope of lead (), the immediate progenitor of is an isotope of radium (), which decays by alpha emission to. Since is found virtually everywhere throughout the Earth’s crust, and are present ubiquitously in soil, rocks, and water [2] . As a consequence, radon is the largest single source of natural environmental radiation [3] . The mean radon dose of 2 mSv/year is greater than 50% of the mean background rate (3.6 mSv/year) [4] , and inhalation of radon is the leading cause of lung cancer among non-smokers [5] . Thus, the determination of radon concentration in homes and workplaces is not only of physical interest, but also a matter of public health concern.

Adverse health effects from inhalation come primarily, not from the gaseous radon itself, but from its radioactive progeny―isotopes of polonium (,), lead (), and bismuth (), formed in the chain of decays shown in the top row below:

(1)

Since the half-life of is 22.3 years, an effective secular equilibrium among the nuclides of the top row can be established within a few hours. Over decades, however, decays, as shown in the bottom row, to the truly stable isotope.

In response to increasing concern over radon exposure, numerous investigations have been made during the past 50 years or more regarding methods to measure indoor radon concentration. This is a vast subject, but, in brief, contemporary radon methodologies ordinarily employ one of the following categories of detectors [6] : 1) alpha track, 2) activated charcoal, 3) electret ion chamber, 4) electronic integrator, 5) continuous monitor. Although these five methodologies differ significantly in their modes of operation, characteristic uncertainties, sampling periods, and costs, they all share an important common feature, viz. the radon gas is collected in a fixed specified control volume. Without a built-in control volume in the design of each of these devices, there would be no way to convert the actual measurements made per unit of time―e.g. a number of particle counts, or number of particle tracks, or accumulated ion charge― into a concentration (or density), i.e. a number of radon atoms per unit of volume. From the known sampling time and calculated particle density, one can ultimately calculate the radioactivity expressed in Bq/m^{3} or (particularly in the US) in pCi/L, where 1 becquerel (Bq) = 1 disintegration/second, and 1 pico-curie per liter (pCi/L) = 37 Bq/m^{3}.

From the standpoint of the convenience of end users of radon activity information, who in the majority of cases are not nuclear scientists or engineers with sophisticated instrumentation to perform radon activity measurements themselves, each of the preceding numbered methodologies, with perhaps exception of (5), requires that the collected sample bearing a fixed volume of radon be sent away to a testing laboratory for analysis. The delay time of at least a few days is not merely an inconvenience, but can significantly affect the uncertainty of the estimated radon activity as a result of sample decay throughout the transit period and waiting time to measurement. (Recall that the half-life of is about 3.8 days, and the half-lives of the progeny in the chain of transmutations (1) leading to (quasi-stable) is much shorter than that.)

Commercially available continuous radon monitors (methodology 5 above), pose other kinds of problems, particularly to researchers who need to monitor the short- term radon concentration in their laboratories. University physics labs are often located in the basements of buildings with concrete floors and cinder-block walls, i.e. materials through which radon, originating in the underlying soil, can diffuse and which, themselves, give rise to radon exhalation [7] . As a gas denser than air, the radon subsequently accumulates close to ground-level regions of a room where apparatus is located and researchers work. Whereas lay persons (e.g. home residents) might take the readings of their radon monitors at face value, professional researchers, especially physicists, often need to know in detail what the device measures and how it calculates the displayed radon activity. These critical details are regarded as commercial proprietary information, which is not revealed even when requested for academic research purposes [8] .

Another drawback is that continuous radon monitors usually record a signal that is proportional to the absolute radon concentration, but for which the proportionality “constant” is an empirical number not deducible from first principles. It must be determined by calibration with a primary standard, can vary among the units sold, and is built into the device algorithm which, as mentioned, is a proprietary secret. Thus, apart from maintaining their own continuously operating nuclear counting equipment, which is (a) costly, (b) ties up apparatus needed for other purposes, and (c) ordinarily requires expertise in alpha, beta, and/or gamma spectroscopy (for examples of such setups, see [9] [10] ), interested researchers cannot themselves verify the radon concentration displayed by a commercial monitor or even ascertain how the engineers or technicians who programmed the monitor calculated the result.

Last, but by no means insignificant, is the relatively high statistical uncertainty in the values of radon activity concentrations measured by commercially available continuous monitors over the short term. The literature of one such state-of-the-art product, claimed by the manufacturer as certified for official radon testing everywhere (except in the US and Canada for which they sell a more expensive model) assigns a relative uncertainty (ratio of standard deviation to mean) of 20% to the displayed 7-day mean activity concentration [11] . However, researchers in laboratories or residents in homes with wide short-term variation in radon levels may need to keep track of the mean daily radon activity. Assuming Gaussian statistics (either as an approximation to Poisson statistics or by virtue of the Central Limit Theorem [12] ), the corresponding relative uncertainty in the mean 24-hour measurement of the commercial monitor would be larger by, yielding a total relative uncertainty of about 53%.

This paper describes a simple, novel method to determine radon concentration directly by means of alpha particle counting in an open volume (rather than fixed control volume) with the most commonly available, overall least expensive, and most easily operated type of nuclear instrumentation: a standard pancake Geiger-Mueller (GM) radiation counter. Although there have been previous attempts to use GM counters to detect the presence of radon [13] , this paper presents a methodology and theoretical analysis to measure radon activity concentration quantitatively and accurately. Among the advantages of this method are:

・ Measurement results are obtained directly after sampling; no external testing laboratory is required to convert the raw data into an absolute radon activity concentration.

・ A relatively short measurement period of 24 hours suffices to provide adequate precision (relative uncertainty of about 10% at 100 Bq/m^{3}) for most personal health- related purposes, such as to ascertain whether a local radon concentration exceeds the published standards set by the US Environmental Protection Agency [5] at 4 pCi/L, or the average concentration reference level of 100 Bq/m^{3} set by the World Health Organization [14] .

・ The operational relations derived in this paper, based on fundamental physical principles and known properties of alpha particles and detector materials, are easy to implement and straightforward to interpret.

The simplicity of the measurement method together with the transparency of the analysis makes it possible for a broad demographic such as researchers, public health officials, and home residents, particularly in developing countries, to estimate local radon activity concentration in real time without the need of specialized and expensive nuclear equipment and without having to rely on public or private testing laboratories or on the readings of commercial devices whose operation and programming are shrouded in secrecy.

2. Experimental Protocol and Analysis

2.1. Apparatus and Methodology

The method described in this paper calls for use of two pancake-style GM radiation counters: one with a detector window unobstructed so as to record incident, , and radiations; the other whose detector window is covered by a sheet of paper to block incident radiation but let pass and radiations. Note that the detector aperture itself is covered by an ultra-thin layer of material such as mica that permits passage of alpha particles. The two GM counters are placed at the same horizontal level (cm above the floor) adjacent to one another at a distance of meter from walls, so as not to be affected by alpha emission from radionuclides in the building materials. The two counters are begun at the same time and set to count for a specified period, e.g. 24 hours. Upon termination of sampling, the difference in readings of the two counters can be converted into a radon-generated alpha count per unit of time from which the radon activity concentration can be calculated by means of the theory developed in this paper. Although there is no artificial control volume (apart from the boundaries of the entire room), the region of alpha particle detectability is actually tightly constrained by the physics of radon diffusion in air and alpha particle interactions in matter.

As the only naturally occurring radioactive gas to emanate from common building materials or percolate into a room from the underlying soil, radon is the only alpha emitter and progenitor of alpha emitters (isotopes of polonium) that can contribute to the GM alpha particle count rate under the conditions described above. Alpha particles issuing from radioactive materials in the walls, floor, and ceiling are not detected because of the intrinsically short range of alpha particles in matter. As a relatively heavy ion (in comparison to the much lighter beta particle), an alpha particle emitted in radioactive decay loses energy almost entirely by creating electron-ion pairs in its passage through matter [15] . Each collision in air, for example, absorbs about 35 eV of energy. Thus, a 5.5 MeV alpha particle emitted in the decay of to, creates about electron-ion pairs in air before abruptly coming to rest.

Apart from slight variations due to fluctuations in air density, an alpha track is effectively straight with a range closely correlated with initial energy E according to the relation [16] [17]

(2)

where range is in cm and energy E is in MeV for. For example, from Equation (2) the range in air of the 5.5 MeV alpha emitted by is 4.03 cm.

For passage through condensed matter, the range of an alpha particle can be calculated from its range in air by means of the Bragg-Kleeman rule [16] [17] , one form of which is

(3)

where signifies mass density and A the effective atomic mass number. The effective atomic mass number for a compound of molar mass M comprising atoms of type i of elemental mass number is given by the relation

(4)

with weight

. (5)

For a mixture of molecules (such as air) of effective mass M, is to be interpreted as the mass fraction (number of atoms × volume fraction) of molecular constituent i.

For purposes of this paper, the molar (or volume) composition of air is 78% N_{2}, 21% O_{2}, and 1% Ar [18] . Hence, the effective molar mass of air is. From Equation (5) it then follows that

, (6)

leading to the effective atomic mass number of air

. (7)

Given the density of air (at 1 atm and 20˚C) and from Equation (7), the Bragg-Kleeman relation (3) reduces to

. (8)

As applied to a GM window, which is here taken to be Muscovite mica (KAl_{3}Si_{3}O_{10}(OH)_{2}) of density [19] and effective atomic mass number (from Equations ((4) and (5))), Equation (8) leads to a range of 20.4 mm for a 5.5 MeV alpha particle.

The above examples serve to illustrate the very short range of alpha particles produced by radon and its progeny and provide a basis for the effective alpha range function and alpha transmission function to be constructed shortly.

2.2. Relation of Radon Alpha Particle Count to Radon Activity Concentration

It is assumed in this paper that all radionuclides decay exponentially in time in accordance with standard nuclear physics, although this assumption has been challenged in the literature and is a matter of current research. For example, see [20] [21] for experimental tests of the standard theory of radioactive decay. According to standard nuclear physics, the probability of a nuclear transmutation within short time interval dt is, where is the associated decay rate constant. The activity or number of decays per unit of time, is then given by

(9)

where is the number of radioactive nuclei in the sample at time t. Solving Equation (9) results in the familiar exponential decay law

(10)

from which follows the relation between half-life and decay rate:.

If, as in the case of, the progeny themselves decay, leading to a chain of K unstable daughter products, then the sequence of decays, as denoted in relation (1) is described by the Bateman rate equations [22]

(11)

where is the number of atoms. If the sequence is terminated at, then.

A general method to solve coupled linear differential equations like (11) is given in [23] . The solution to (11) can also be expressed in iterative form [24]

(12)

which is particularly useful for numerical solution by computer. Of relevance to the present work is the limiting case known as secular equilibrium, a steady-state condition of equal activities between a long-lived parent and much shorter-lived daughter radionuclides. Mathematically, the criteria for secular equilibrium are

. (13)

In light of the preceding background, consider a GM detector with circular window of radius over whose surface slowly diffuses a current density of atoms of steady-state number density. Assuming azimuthal symmetry of the detector and steady-state radon density in the vicinity of the detector throughout the time of measurement, the differential number of alpha counts in sampling time T from radon atoms in a differential volume element located at polar coordinates relative to the center of the window (Figure 1) is given by the relation

(14)

in which the various factors are defined as follows:

= mean number of alpha particles from radon and its progeny contributing to the count for each decay.

= decay rate of = [25] .

= alpha particle range function in air = fraction of alpha particles reaching a distance r in air from their source (Figure 2).

= solid angle of the GM window subtended at the point of alpha emission (Figure 1, Figure 3).

= alpha transmission function = surface area of the GM window through which an alpha particle can penetrate to the interior fill gas (Figures 4(a)-(c)).

= intrinsic alpha particle detection efficiency (for energy ³ 4 Mev). Under normal operating conditions for a GM tube, even a single ion pair formed within the fill gas can trigger a full electric discharge. Thus, the intrinsic efficiency for detecting alpha particles is [26] . From this point on, it will be assumed that.

Figure 1. Azimuthally symmetric detection geometry. A radon atom located at coordinates decays, producing an alpha particle that reaches the detector window of radius if emitted within solid angle.

Figure 2. 5.5 MeV alpha particle range function (in air), Equation (41).

Figure 3. Fractional solid angle of detection, Equation (46), as a function of polar angle for radial coordinates (in cm): 0.5 (red), 1.0 (blue), 2.0 (violet), a = 2.25 (black), 3.0 (orange), 4.0 (green).

(a)(b)(c)

Figure 4. (a) Geometry of transmission of alpha particle through mica window of radius and thickness d. The alpha, emitted at point relative to the vertical symmetry axis (z) of the window, must impact the window at angle to the surface into order to pass through the mica along a path £ range; (b) a line at constant rotated about the vertical axis through generates a circle of radius; (c) the lenticular intersection (dark gray) of this circle with the circular mica window of radius is the region of alpha transmission.

Explicit expressions for the functions, , and will be given in due course. Once known, however, they can be used to calculate the following useful quantities:

Effective detection volume (in the half space above the GM window):

. (15)

Effective cross section of radon flow through surface normal to the GM window surface:

. (16)

The upper limit to the radial integrals (15) and (16) is given as, although in principle it can not exceed the dimensions of the room. Operationally, the function, which takes the form illustrated in Figure 2, constrains the upper limit to a size comparable to the alpha range in air.

The overall efficiency of the detection process is given by the range-weighted volume average of the product of the fractional quantities and

. (17)

Thus, the volume average of Equation (14) takes the much simplified form

. (18)

The number of radon atoms diffusing over the detector window through cross section in time dt is given by

(19)

where the radon current density

(20)

is the product of the diffusion velocity and radon atomic density. Equation (20) and explicit expressions for

(21)

and a characteristic diffusion length

(22)

are obtained from solution of the diffusion equation based on Fick’s law [27] with implementation of appropriate boundary conditions. This analysis is given in Appendix 1 for the one-dimensional case in which radon diffuses toward the detector from one wall. The experimental radon diffusion coefficient in air is [28] [29]

. (23)

The diffusion of radon over the detector results in an effective residence time in the detection volume

(24)

analogous to the residence time in radioactive flow measurements [30] such as employed in high-performance radioactive liquid chromatography [31] . The principal differences between standard methods of radioactive flow measurement and the novel GM method of radon measurement reported here are: (A) in the former the detection volume is fixed by the design of the apparatus, the flow rate is a deterministic variable adjustable by the experimenter, and a finite sample of radioactive material is used, whereas (B) in the latter the detection volume is determined by the energy-dependent range of the alpha particle, the flow rate is a stochastic variable determined by the diffusion of radon in air, and the quantity of radioactive material is likewise a stochastic variable determined by the duration of sampling.

Substitution of Equations ((21) and (24)) into Equation (18) leads to the compact expression

(25)

where is the observed alpha count rate, and the radon activity per unit of volume (activity concentration) is

. (26)

Inversion of Equation (26) therefore yields the sought-for radon activity concentration (e.g. in Bq/m^{3})

. (27)

The bracketed product in the numerator of Equation (25) yields a dimensionless quantity, the number of radon disintegrations that occur in a time interval given by the residence time in the denominator. In other words, under conditions of 100% detection efficiency and alpha emission from radon only, the number of alphas counted during time T would equal the actual number of radon disintegrations in the detection volume during time.

2.3. Alpha Emission from Radon Progeny

Derivation of the activity concentration relation (27) assumed (a) effectively equivalent ranges for the three alphas produced by radon and polonium decays and (b) secular equilibrium among radon and its progeny. Examination of these two assumptions will lead to an extension of relation (27).

Consider assumption (a) first. The energy of an alpha particle determines its range and therefore the effective range function, detection volume (15), flux cross section (16), GM window transmission function, and detection efficiency (17). Assuming a stationary-state condition, although not necessarily secular equilibrium (to be examined next), among radon and its progeny, the analysis leading to Equation (25) can be extended to yield

(28)

where indices respectively refer to alpha particles from, ,. The functions are calculated according to the expressions given in the previous section (and sections to follow) with insertion of the appropriate alpha range for air and for mica in the associated range function and transmission function. The fractions, which sum to unity, give the relative proportions and effective alpha number

(29)

of the activities and in the equilibrium mixture within the detection volume. Note that the subscript i on the activity symbol corresponds to the progeny number (and not the alpha particle label j) in the sequence (1) of daughter nuclides of.

Table 1 summarizes the energies and ranges of the three alphas that potentially contribute to the GM count from decays of, , and.

The ranges in air and mica, calculated from Equations ((2) and (8)) respectively, increase approximately as the 3/2 power of alpha energy.

It is important to bear in mind that the three residence times

(30)

still depend on the diffusion velocity of. Table 2 summarizes the diffusion properties of the three alpha-emitting nuclides with diffusion coefficients D taken from Reference [32] . Only radon has a macroscopic characteristic diffusion length on the order of meters compared with the diffusion lengths of the two polonium progeny, which are on the order of a few centimeters and nanometers, respectively. Thus, in the flow Equation (19) only radon, and not its progeny, reach the detector from points of origin well outside the alpha detection volume, a distance of at least 1 m from the walls. (Note: a number in parentheses in Table 2 signifies a power of 10.)

Consider next assumption (b) regarding secular equilibrium, expressed by relation (13). Solution (12) of the Bateman equation (11) for a closed system of radon and its

Table 1. Alpha particle energy and mean range.

Table 2. Diffusion properties of radon and alpha emitting progeny.

progeny leads to the results displayed graphically in Figure 5. As seen in the figure, secular equilibrium is established within about 4 hours, except for whose activity remains close to 0 over sampling times short in comparison to its half-life of 22 years. Under conditions of secular equilibrium, the three alpha-emitting nuclides in the decay sequence from to have the same activity, and it follows from (29) that

(SECULAR EQUILIBRIUM). (31)

Depending on environmental conditions (e.g. relative humidity, temperature, aerosol content, and other factors), the radon progeny, whether initially charged or neutral, are subject shortly after generation to physical and chemical interactions with the surfaces of a room as well as with particles and molecules within the room’s atmosphere. Because of these interactions, the daughter nuclides―in particular, ,―can deposit (“plate out”) on solid surfaces or become part of airborne molecular clusters [33] . The resulting equilibrium activities and equilibrium alpha number can therefore differ from the values (31) for secular equilibrium. Note, however, that, since decays nearly instantly, it has the same activity as its progenitor.

In place of the Bateman equations, the new equilibrium conditions can be determined from solution of the Jacobi equations [34] [35] , which in their simplest form for a closed, unventilated room with clean air (corresponding to experimental conditions in this paper), become

(32)

Figure 5. Evolution in time of the ratio of activities of and progeny to the activity of at: (red), (blue), (green), and (black), (brown). Secular equilibrium is achieved in about 4 hours, except for whose half-life is 22.3 y. The activity of, whose half-life is 165 μs, is practically identical to that of its progenitor,.

where

(33)

is the loss rate, is termed the deposition velocity, and is the surface to volume ratio of the room. The equilibrium progeny activities relative to the activity of radon obtained from the steady-state solution to Equation (32) are given by

. (34)

In the experimental section of this paper indoor radon activity concentrations are reported for measurements taken in an unventilated basement room with cement block walls and very low dust and aerosol content in the air. Progeny resulting from radon decay are more likely to become part of air molecular clusters than to plate out. Experimental values for the deposition velocity of attached radon progeny have been found experimentally to range from about 0.03 to 0.20 m/h [36] . To estimate the most unbiased distribution of a physical variable consistent with known information, statistical physicists employ the principle of maximum entropy (PME) [37] [38] . Given only the range of the deposition velocity, the PME leads to a uniform distribution with respective mean and variance

(35)

which is in excellent accord with the value 0.0117 m/h obtained from a recent statistical fit to data taken in indoor dwellings [39] .

Substitution of and the value (for the room in which measurements were made) into Equations ((33) and (34)) leads to equilibrium activity ratios

(36)

from which follow the equilibrium proportions and number of alpha particles per radon decay

(JACOBI EQUILIBRIUM). (37)

Inversion of Equation (28) then yields the general operational relation for radon activity concentration measured with two GM counters in sampling time T

. (38)

With the values given in Equations ((31) and (37)), the constant ratio of radon activity to alpha count rate from Equation (38) is

.(39)

To good approximation, the standard deviation of activity (38) is then

(40)

in which the uncertainty in alpha count rate is governed by Poisson-Gauss statistics. (Note: The characteristic feature of Poisson statistics is that the variance equals the mean. For a mean count per sampling interval, the Poisson distribution is adequately represented by a Gaussian distribution with.)

It is to be noted explicitly that relation (39) does not depend on a non-calculable empirical proportionality constant that must be determined independently by calibration against a primary standard. Rather, the physical quantities whose numerical values enter Equation (38) are measurable properties of particles (e.g. alphas), materials (e.g. mica), and local environment (e.g. air). As such, the more accurately these physical quantities are known for the conditions under which measurements are made, the better will be the resulting estimate of radon activity concentration. Nevertheless, the numerical values in relation (39) show that, for moderate levels of radon in an unventilated room with clean air, the analysis based on Equation (28) with equilibrium values (37) gives results comparable to the analysis which assumed secular equilibrium values (31).

2.4. Alpha Range Function Q(r)

Figure 2 shows a nearly identical theoretical replica of the experimental transmissivity―i.e. fraction of particles that pass through an absorbing medium as a function of thickness―of a 5.5 MeV alpha particle in air [40] . The probability of transmission is essentially 100% until the alpha particle has slowed sufficiently to generate a very large number of ion pairs and thereby rapidly come to rest. The narrow variation of alpha ranges about the mean range (defined at 50% transmission) follows a normal distribution [41] with relative uncertainty of about 5% [42] . The overall shape of the empirical range function closely resembles the occupation probability function of a Fermi-Dirac particle [43] , and can be accurately represented by an expression of the form

(41)

where is the alpha mean range parameter calculable from Equation (2) for air (and Equation (3) for mica), and is a fall-off parameter that fits the experimental transmissivity curve. The ranges of the three alpha particles are summarized in Table 1. For all three range functions in air, the value provides a satisfactory match.

2.5. Solid Angle Subtended by the Detector Window

The solid angle subtended by a detector in an experiment to measure nuclear radiations is ordinarily defined only for configurations with a point-like or planar source symmetrically located on or about the symmetry axis normal to the detecting surface [44] . Such configurations do not apply to the measurement of alpha particles emitted by radon and polonium atoms located randomly with respect to the detector in a three- dimensional volume. For this configuration a more general relation is needed.

Figure 1 shows a radon atom at arbitrary radial distance r from the center of the azimuthally symmetric GM window and at polar angle to the vertical symmetry axis of the window. Of the steradians into which the atom can emit an alpha particle, only the shaded portion of the unit sphere centered on the atom permits the alpha to arrive at the window surface. From basic geometry, the fraction of the surface of the unit sphere is given by

(42)

where, from Figure 1, the conical half-angle of subtended at the center of the unit sphere is

. (43)

The problem, then, is to express angles and in terms of the integration variables r and.

The first step is to obtain the lengths

(44)

from the law of cosines. The second step is to express in terms of the sides of a triangle from the relation

. (45)

After some algebraic manipulation, it then follows that

. (46)

As a quick plausibility check of Equation (46), set (source on the normal symmetry axis), whereupon one obtains the expected result [45]

. (47)

Figure 3 shows the variation in fractional solid angle (46) as a function of polar angle for different values of radial coordinate. Examination of Figure 1 shows immediately why equals the maximum value 1/2 for radial coordinate and equals 0 for. The solid black line in the figure denotes the threshold value.

2.6. Alpha Particle Transmission Function

Upon arrival at the surface of the GM mica window, an alpha particle will fail to pass through the mica and ionize the interior fill gas if its path length exceeds the range in mica, as given in Table 2. This constraint limits the minimum incident angle between the path of the particle and the window surface to

(48)

where d is the window thickness, as shown in Figure 4(a).

Examination of Figure 4(a) and Figure 4(b) shows that the portion of the window surface through which an alpha can be transmitted cannot exceed the area of the circle of radius from the center (directly below the point of alpha emission) to the point of incidence at angle to the surface. However, depending on the location of the alpha at the instant of emission, a portion of this circle can extend beyond the area of the detector window (Figure 4(b)), which is, itself, a circle of radius with center at O. Thus, the relevant transmission area is the lenticular intersection of these two circles, shown as the dark shaded area in Figure 4(c).

The area of intersection of two circles of respective radii and, whose centers are separated by the distance, is readily derivable from integral calculus or from theorems of plane geometry [46] and is given by

(49)

where, in the context of the geometry of the GM window,

. (50)

The operation Re in Equation (49) to “take the real part” of the bracketed expression is necessary to include the case of non-overlapping circles. In that case, expression (49) returns 0; otherwise it can return an imaginary number if. In all other cases, the actual area of intersection is a real positive number. The operation Re assures that the integral (17) leads to a correct real value for the overall detection efficiency.

Figure 6 shows the variation in for the 5.5 MeV alpha as a function of for different values of r. As expected from the geometry of Figure 4(a), the alpha transmissivity of the mica window at fixed r is greatest for close to 0 and decreases toward 0 as approaches. For r greater than about 0.84 cm, there is a portion of the angular range for which the entire detection window is transmissible.

Figure 6. Alpha transmissivity, Equation (49), in mica as a function of polar angle for different values of the radial coordinate (in cm) of the emission point: 2.0 (red), 1.0 (blue), 0.8 (green), 0.6 (violet), 0.4 (orange), 0.2 (plum).

Figure 7. Radial dependence of factors contributing to detection efficiency of 5.5 MeV alpha emitted from a point source on the normal symmetry axis. (a) alpha transmissivity in mica (blue); (b) solid angle (red); (c) efficiency (dashed black); (d) integrated efficiency (violet).

2.7. Alpha Particle Detection Efficiency

Figure 7 shows the variation with distance r from the center of the mica window of various functions contributing to detection efficiency of a 5.5 MeV alpha emitted at r on the vertical symmetry axis of the window (i.e. standard configuration for many beam experiments). Plot 7a (solid blue) depicts the transmissivity. Alphas emitted at a distance r exceeding about 0.8 cm from the center of the window are transmissible through the entire window surface. For decreasing distances below about 0.8 cm, the area of transmission rapidly falls toward 0. This effect is purely geometric. The circular area of alpha transmissibility may be thought of as the base of a cone of height r and fixed apex angle complementary to, Equation (48), as shown in Figure 4(a) & Figure 4(b). As the point of alpha emission (i.e. the apex of the cone) approaches the window surface (i.e.), the area of the base decreases as. Plot 7b (solid red) depicts the solid angle subtended at the detector. The behavior is effectively opposite that of the transmissivity, achieving its maximum value of 0.5 for alphas emitted directly at the window surface and monotonically dropping to 0 with increasing separation. The product of the two functions yields the detection efficiency of the alpha from a point source, which is seen in plot 7c (dashed black) to reach a maximum of about 32% for a source separation of approximately 0.8 cm. Finally, plot 7d (solid violet) shows the directionally integrated efficiency

, (51)

which peaks at about 15% at a source separation of 1.13 cm.

The overall detection efficiencies of the individual 5.5, 6.0 and 7.7 MeV alpha particles arriving at the detector from anywhere within the half-space above the GM window are obtained by substitution in Equation (17) of the respective range functions to yield

. (52)

The mean alpha detection efficiencies under conditions of secular equilibrium (31) and Jacobi equilibrium (37) are

, (53)

which agree closely with the empirically determined alpha efficiency 4.2% provided by the manufacturer of the author’s GM counters [47] .

3. Short-Term and Long-Term Measurements of Indoor Radon Concentration

A comparison of the GM method with a commercial radon monitor to measure short- term and long-term indoor radon concentrations was made by means of two Inspector Radiation Alert detectors (to be referred to as GM1 and GM2) manufactured by S.E. International Inc (SEI), each with a halogen-quenched, pancake-style GM tube with mica window of areal density ~2.0 mg/cm^{2} and effective window diameter 4.5 cm. The manufacturer-specified accuracy (in counts per minute) is, tested by the US National Institutes of Standards and Technology (NIST). The detector is sensitive to alpha particles down to 2 MeV, beta particles down to 0.16 MeV, and gamma photons down to 10 KeV. A prior (bg)-background count of duration 192 hours showed that GM2 was slightly more sensitive than GM1 by, an asymmetry taken into account in the subsequent data analysis. The GM display shows the mean count for the previous 30-second time interval, updated at intervals of 3 seconds.

The two GM counters were placed, as described in Section 2.1, on a raised (0.5 m) platform in a basement room with concrete walls, at a distance of 1.5 m from the one wall that was lower than the outside ground level due to the exterior sloping terrain. From data of radon penetration and mobility [48] , it can be reasonably concluded that the main source of radon in the room was the radon flux emerging primarily from this wall. Distances to the other walls well exceeded 2 m.

On the same platform and at the same distance from the subterranean wall as the GM counters, were placed two Corentium Home radon monitors [49] to serve for comparison. The Corentium monitor performs hourly samplings of radon through a passive diffusion chamber by means of alpha spectrometry with a silicon photodiode. The sensitivity of the instrument claimed by the manufacturer is 5.55 counts per hour (cph) per pCi/L, or, equivalently, a relative uncertainty of 20% at 100 Bq/m^{3} after 7 days. As pointed out in Section 1, the relative uncertainty (ratio of standard deviation to mean) of a 24-hour measurement would then be about 53%. The monitor displays three readings, the mean radon activity concentration of the preceding (a) 24 hours, (b) 7 days, and (c) period up to 1 year. The measurement durations of relevance to the experiment reported here are 24 hours and 30 days.

Figure 8 displays graphically the sequence of 24-hour radon measurements made over a period of 30 days by (a) the method using two GM counters (red) with radon activity calculated by Equations ((37) and (38)), (b) Corentium monitor No. 1 referred to as COR1 (blue), (c) Corentium monitor No. 2 referred to as COR2 (green), and the mean value of the two Corentium readings referred to as ACOR (black). For radon activities of about 100 Bq/m^{3} (or 2.7 pCi/L), GM1 recorded about 62,300 counts and GM2 recorded about 67,800 counts in 24 h, resulting in a net alpha count rate of approximately in basic agreement with relation (39). The algorithms for calculating the alpha count rate (cpm) and associated uncertainty from the raw data are given in Appendix 2. A complete tabulation of the GM and COR data is given in Table 3.

As Figure 8 shows, the radon concentrations measured by GM fall close to the center of the wide ACOR error bars in nearly every sample except the first when COR1 and COR2 were initially put into service and significantly underestimated the radon level in the room. This is a feature acknowledged by the manufacturer, who markets the Corentium monitor as useful primarily for sample periods well in excess of 48 h, and preferably at least 7 days. It is to be noted, therefore, that the 24 h measurements of COR1 and COR2 often varied to a greater extent between each other than each did with the GM measurement. The relative uncertainty of a 24 h GM-measured radon concentration is approximately at 100 Bq/m^{3}, in comparison with for

Figure 8. Sequence of 24-hour measurements of radon activity by GM method (red), commercial monitor COR1 (blue), commercial monitor COR2 (green), mean value ACOR of the two commercial monitors (black). The standard deviation of radon activity concentrations by (a) GM is approximately 0.32 pCi/L, (b) commercial monitor ranges from about 0.90 to 1.44.

the mean reading of COR1 and COR2.

The method to measure radon concentration with GM counters, although especially useful for short-term (24 h) sample periods, also provides accurate long-term radon levels with an uncertainty that decreases with the square root of sampling time. From the data in Table 3 one can calculate the mean radon concentration over the 30-day sampling period for each modality (GM vs commercial monitor):

. (54)

The results (54) show that, with respect to the Corentium monitors which are marketed as a standard for radon testing, the proposed GM-based method yielded a statistically equivalent mean long-term value of radon concentration. Assuming that the two sets of measurements are samples from normal distributions, the statistic z for testing the equivalence of two means [50]

(55)

is itself normally distributed as―i.e. as a standard normal variate of 0 mean and unit variance. It then follows that the probability that a subsequent series of GM measurements would yield a value greater or lesser than twice the observed value in (55) (that is, two standard deviations from the mean difference of 0) is

. (56)

The conventional statistical threshold for judging whether an event could have occurred through pure chance is 5%. Since probability (56) well exceeds the 5% threshold, the difference between the two means in (54) is regarded as statistically insignificant. It

Table 3. GM counts and radon activity (24 h sampling over 30 days).

SE = secular equilibrium; JE = Jacobi equilibrium

is worth emphasizing, however, that the difference of the two mean long-term measurements in (54) would be statistically significant if those two concentrations were both measured by the GM method. In that case, the relevant statistic would be

(57)

with associated p-value

(58)

which is considerably below the 5% threshold of significance. The reason for the difference in outcomes between (56) and (58) is due to the much lower standard error of the GM measurements in comparison to measurements by the commercial monitor. This seminal point is elucidated further in Figure 9.

Figure 9 shows a plot of the ACOR (24 h) radon concentrations (black points) as a function of the corresponding GM-measured alpha count rate. Superposed on the plot is the linear relation (38) for conditions of secular equilibrium (blue) and Jacobi equilibrium (red). The slopes of the two lines, given by expression (39), fall within the error bars (standard deviation) of the commercial monitor. Also superposed is the maximum-likelihood (ML) line of regression [51] (green) to the ACOR measurements with slope, which signifies statistically a flat line, although a slight rise of the data points with increasing alpha cpm values is marginally visible if the points near 3.6 cpm are disregarded. Assuming that the measurements of the commercial monitor comprise a linear trend with normally distributed random noise, one can interpret Figure 9 in the following heuristic way. In the limit of numerous 24 h measurements, the accumulation of black points along each vertical line of constant alpha cpm would

Figure 9. Plot of ACOR radon activity concentration (black) against corresponding GM alpha count rate. The three superposed lines show the GM activities based on Jacobi equilibrium (red) and secular equilibrium (blue), as well as the ML line of regression (green) to the ACOR points.

be distributed with a Gaussian density centered on the intersecting theoretical red or blue line (depending on equilibrium conditions) given by Equation (38). One can infer from the figure that many such points would be required before the ML line of regression to the radon concentration estimated by commercial monitor would be statistically equivalent to theoretical relation (38) of the GM method.

The GM counters and commercial radon monitors used in this experimental comparison both sampled environmental radon through a process of passive diffusion. One may inquire, therefore, why it is that the GM method can make consistent short-term measurements with much lower statistical uncertainty, whereas short-term measurements of the commercial monitor show wide scatter and insensitivity to the total alpha count rate. An explanation, at least in part (since the programming of the commercial monitor is a proprietary secret), is that the built-in control volume according to Corentium is 24 cm^{3}, whereas the effective detection volumes of the 5.5, 6.0, and 7.7 MeV alpha particles in open air are, according to Equation (15), about 139, 210, and 700 cm^{3} respectively.

4. Conclusions

In the comprehensive description of methodologies in the World Health Organization Handbook [6] for measuring indoor radon activity (summarized in Section 1), the use of GM counters is not included. Indeed, one can find explicit statements from manufacturers of GM counters [52] that their utility for detecting (let alone actually measuring) radon is debatable. Reasons given or implied include (a) GM tubes do not detect radon directly; (b) GM tubes cannot determine whether a detected alpha came from radon or its progeny; (c) GM efficiency for alpha detection is low; (d) there is no fixed volume within which a radon sample can be confined and its concentration measured.

While the preceding statements are not altogether incorrect, they do not constitute valid arguments against the use of GM counters as an effective methodology to measure―and not merely detect―indoor concentration of radon. The statements presuppose the use of only a single GM counter and assume the necessity of an initially determined volume of radon. As shown in this paper, neither supposition is justified.

In this paper I describe and demonstrate a method to measure indoor radon using two GM counters, which, together with the theoretical analysis derived here for converting an alpha count rate into a radon concentration, provides a number of advantages over alternative methodologies:

・ The method is especially suitable for daily monitoring of radon since it provides an activity concentration after just 24 h in situ sampling that eliminates the need to send a prepared sample to a testing laboratory. For continuous monitoring, this saves users much time and considerable expense. Other methodologies (e.g. alpha-track counting) may have higher accuracy, but are not a convenient or economical option for researchers or home dwellers who require frequent or continuous radon monitoring.

・ The GM method yields more precise mean short-term and long-term radon concentrations than do commercial radon monitors employing alpha detection in a passive diffusion chamber.

・ The GM method yields radon concentrations of comparable accuracy to that of commercial monitors employing sampling by passive diffusion.

・ The method employs the simplest, most versatile, and overall least expensive radiation detection instrumentation (simply two pancake GM counters) for non-nuclear researchers and home dwellers who are not trained in nuclear spectrometry or do not have their own nuclear spectrometry instrumentation. As such, the GM method can be especially useful in geographical areas where access to specialized radon testing equipment or testing laboratories is inadequate. Or, in fact, useful in any milieu for people who want to monitor indoor radon levels themselves.

The principles behind the experimental method and theoretical analysis described in this paper are physically verifiable and fully transparent (in contrast to the operation and programming of commercial monitors which are patent-protected intellectual property). The use of two monitors, so placed as to receive alpha particles only from radon and its polonium progeny, makes it possible to separate the alpha signal from the beta and gamma background. Whereas a commercial radon monitor may use energy selection to count alphas from radon only, the method described here counts alphas from both sources and uses the laws of physics and properties of materials to calculate the mean number of alphas detected per radon decay.

An essentially novel feature to the method described here is the means by which the alpha count rate is converted into a radon concentration even though there is no built-in control volume. This is accomplished by calculating an effective detection volume (15), detection cross section (16), and residence time (24), which are determined by the physical laws governing alpha particle interactions in matter and the diffusion of atoms in air. Moreover, the theoretical analysis presented here can be applied to different environmental conditions of indoor radon measurement, once the effects of these conditions on the diffusion of radon gas and the equilibrium of radon and its progeny have been ascertained.

Appendix 1: Radon Diffusion

Mass conservation of nuclei in a process of one-dimensional diffusive flow with nuclear decay is expressed by the differential equation

(59)

in which is the particle density, is the particle current density, is the constant decay rate, and the relation between current density and concentration gradient is expressed by Fick’s law

(60)

with diffusion constant D. In the case of a steady-state flow, the derivative vanishes, and Equation (59) reduces to the form

(61)

with characteristic diffusion length

. (62)

The solution to Equation (61) and the associated current density take the general form

(63)

with diffusion velocity

. (64)

Coefficients are determined by initial conditions and/or boundary conditions. Taking the origin to be the center of the detector window where the particle concentration and current density are measured to be and leads to coefficients

. (65)

For radon measured in an open volume, the condition that and vanish requires that and. It then follows from Equation (65) that

(66)

and from Equation (63) that

. (67)

Appendix 2: Algorithm for Calculating the Alpha Count Rate in cpm

Define the following measurable quantities:

= GM1 24 hour count

= GM2 24 hour count

= (GM2 − GM1) background offset rate (192 h sample time)

= standard deviation of background offset rate

The alpha count rate (cpm) and corresponding standard deviation is then determined from the expressions

(68)

. (69)

Submit or recommend next manuscript to SCIRP and we will provide best service for you:

Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.

A wide selection of journals (inclusive of 9 subjects, more than 200 journals)

Providing 24-hour high-quality service

User-friendly online submission system

Fair and swift peer-review system

Efficient typesetting and proofreading procedure

Display of the result of downloads and visits, as well as the number of cited articles

Maximum dissemination of your research work

Submit your manuscript at: http://papersubmission.scirp.org/

Or contact wjnst@scirp.org

Conflicts of Interest

The authors declare no conflicts of interest.

Cite this paper

*World Journal of Nuclear Science and Technology*,

**6**, 232-260. doi: 10.4236/wjnst.2016.64024.

[1] | L’Annunziata, M.F. (1998) Handbook of Radioactivity Analysis. Academic Press, New York, 749. |

[2] |
De Felice, P. (2007) Primary Standards of Radon. Metrologia, 44, S82-S86. http://dx.doi.org/10.1088/0026-1394/44/4/S11 |

[3] | Bodansky, D. (1996) Nuclear Energy: Principles, Practices, and Prospects. Springer, New York, 33. |

[4] | Hafemeister, D. (2007) Physics of Societal Issues: Calculations on National Security, Environment, and Energy. Springer, New York, 190. |

[5] |
A Citizen’s Guide to Radon (2012) US Environmental Protection Agency EPA-402/K-12/002. http://www.epa.gov/radon |

[6] | Zeeb, H. and Shannoun, F. (Eds.) (2009) WHO Handbook on Indoor Radon. World Health Organization, 20. |

[7] |
Kumar, A. and Chauhan, R.P. (2014) Active and Passive Measurements of Radon Diffusion Coefficient from Building Construction Materials. Environmental Earth Sciences, 72, 251-257. http://dx.doi.org/10.1007/s12665-013-2951-5 |

[8] | Silverman, M.P. (2015) Personal Correspondence with the Corentium AS (Oslo, Norway), Manufacturer of the Corentium Digital Radon Monitor. |

[9] |
Silverman, M.P., Strange, W., Silverman, C.R. and Lipscombe, T.C. (1999) Tests of Alpha-, Beta-, and Electron Capture Decays for Randomness. Physics Letters A, 262, 265-273. http://dx.doi.org/10.1016/S0375-9601(99)00668-4 |

[10] |
Silverman, M.P. and Strange, W. (2009) Search for Correlated Fluctuations in the β + Decay of Na-22. Europhysics Letters, 87, Article ID: 32001. http://dx.doi.org/10.1209/0295-5075/87/32001 |

[11] |
Corentium FAQ. https://corentium.com/ca-en/product/home/ |

[12] | Kendall, M.G. and Stuart, A. (1963) The Advanced Theory of Statistics: Vol 1: Distribution Theory. Hafner, New York, 193-195. |

[13] |
Riley, J.A. (1996) Student Measurement of Radon Gas Concentrations. American Journal of Physics, 64, 72-77. http://dx.doi.org/10.1119/1.18295 |

[14] |
WHO Radon and Health Fact Sheet (2011) http://www.who.int/mediacentre/factsheets/fs291/en/ |

[15] |
Bichsel, H., Groom, D.E. and Klein, S.R. (2015) Chapter 32: Passage of Particles through Matter. Particle Data Group. http://pdg.lbl.gov/2015/reviews/rpp2015-rev-passage-particles-matt er.pdf |

[16] | Lapp, R.E. and Andrews, H.L. (1972) Nuclear Radiation Physics. 4th Edition, Prentice-Hall, Englewood Cliffs, 198-203. |

[17] | Tsoulfanidis, N. and Landsberger, S. (1995) Measurement & Detection of Radiation. 2nd Edition, Taylor & Francis, Washington DC, 132-136. |

[18] |
Schlatter, T.W. (2010) Atmospheric Composition. Encyclopedia of Aerospace Engineering. Wiley, Hoboken. http://dx.doi.org/10.1002/9780470686652.eae319 |

[19] |
Canberra GM Pancake Detectors. http://www.canberra.com/products/detectors/pdf/Pancake-GM-SS- C37516.pdf |

[20] |
Silverman, M.P. (2015) Search for Non-Standard Radioactive Decay Based on Distribution of Activities. Europhysics Letters, 110, Article ID: 52001. http://dx.doi.org/10.1209/0295-5075/110/52001 |

[21] |
Silverman, M.P. (2016) Search for Anomalies in the Decay of Radioactive Mn-54. Europhysics Letters, 114, Article ID: 62001. http://dx.doi.org/10.1209/0295-5075/114/62001 |

[22] | Friedlander, G. and Kennedy, J.W. (1955) Nuclear and Radiochemistry. Wiley, New York, 129-137. |

[23] |
Silverman, M.P. (2014) A Certain Uncertainty: Nature’s Random Ways. Cambridge University Press, New York, 315-318. http://dx.doi.org/10.1017/cbo9781139507370 |

[24] |
Silverman, M.P. and Strange, W. (2004) The Distribution of Composite Measurements: How to Be Certain of the Uncertainties in What We Measure. American Journal of Physics, 72, 1-14. http://dx.doi.org/10.1119/1.1738426 |

[25] |
National Nuclear Data Center. Brookhaven National Laboratory. http://www.nndc.bnl.gov/nudat2/chartNuc.jsp |

[26] | Knoll, G.F. (2010) Radiation Detection and Measurement. 4th Edition, Wiley, Hoboken, 217. |

[27] |
Ghez, R. (2001) Diffusion Phenomena. Springer, New York, 4-9. http://dx.doi.org/10.1007/978-1-4757-3361-7 |

[28] | Nazaroff, W.W. and Nero Jr., A.V. (Eds.) (1988) Radon and Its Decay Products in Indoor Air. Wiley, New York, 65-69. |

[29] |
Environmental Science Division, Argonne National Laboratory. Effective Radon Diffusion Coefficient. http://web.ead.anl.gov/resrad/datacoll/radon.htm |

[30] |
Hull, D.E. (1958) The Total-Count Technique: A New Principle in Flow Measurement. International Journal of Applied Radiation and Isotopes, 4, 1-15. http://dx.doi.org/10.1016/0020-708X(58)90020-6 |

[31] |
Reeve, D.R. and Crozier, A. (1977) Radioactivity Monitor for High-Performance Liquid Chromatography. Journal of Chromatography, 137, 271-282. http://dx.doi.org/10.1016/S0021-9673(00)81350-3 |

[32] | Nazaroff, W.W. and Nero, A.V. (1988) op. cit., pp 66, 206. |

[33] | Nazaroff, W.W. and Nero, A.V. (1988) op. cit pp 206-211 |

[34] |
Jacobi, W. (1972) Activity and Potential Alpha Energy of 222Radon and 220Radon Daughters in Different Air Atmospheres. Health Physics, 22, 441-450. http://dx.doi.org/10.1097/00004032-197205000-00002 |

[35] | Nazaroff, W.W. and Nero, A.V. (1988) op. cit., pp 161, 178-182. |

[36] | Nazaroff, W.W. and Nero, A.V. (1988) op. cit., p 178. |

[37] |
Jaynes, E.T. (1957) Information Theory and Statistical Mechanics. Physical Review, 106, 620-630. http://dx.doi.org/10.1103/PhysRev.106.620 |

[38] | Silverman, M.P. (2014) A Certain Uncertainty, op. cit., pp 49-61. |

[39] |
Mishra, R., Mayya, Y.S. and Kushwaha, H.S. (2009) Measurement of 220Rn/222Rn Progeny Deposition Velocities on Surfaces and Their Comparison with Theoretical Models. Journal of Aerosol Science, 40, 1-15. http://dx.doi.org/10.1016/j.jaerosci.2008.08.001 |

[40] | Knoll, G.F. (2010) op. cit., p. 34. |

[41] | Lapp, R.E. and Andrews, H.L. (1972) op. cit., p 199. |

[42] | Tsoulfanidis, N. and Landsberger, S. (1995) op. cit., p 106. |

[43] | Callen, H.B. (1985) Thermodynamics and an Introduction to Thermostatistics. Wiley, New York, 396-397. |

[44] | Tsoulfanidis, N. and Landsberger, S. (1995) op. cit., pp 218-226. |

[45] | Knoll, G.F. (2010) op. cit., p. 120. |

[46] |
Wolfram Math World. Circle-Circle Intersection. http://mathworld.wolfram.com/Circle-CircleIntersection.html |

[47] | SE International Inc (2015) Inspector USB Table of Built-In Efficiencies. Personal Communication from Chris Thomas, Project Manager, S.E.I. Inc. |

[48] |
Porstendoerfer, J. (1994) Properties and Behaviour of Radon and Thoron and Their Decay Products in the Air. Journal of Aerosol Science, 25, 219-263. http://dx.doi.org/10.1016/0021-8502(94)90077-9 |

[49] |
Corentium (Oslo, Norway). https://corentium.com/product/home/ |

[50] | Ross, S.M. (2004) Probability and Statistics for Scientists and Engineers. Elsevier, New York, 312-314. |

[51] |
Hald, A. (1999) On the History of Maximum Likelihood in Relation to Inverse Probability and Least Squares. Statistical Science, 14, 214-222. http://dx.doi.org/10.1214/ss/1009212248 |

[52] | https://medcom.com/support/frequently-asked-questions/ |

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.