Radionuclide Dating Is Rigorous

[Calilasseia]

Since one of the perennial canards that defenders of valid science have to endure on rationalist forums, a canard that is practically a masturbatory obsession with creationists, is the “radionuclide dating is based upon assumptions” canard, I thought it apposite to produce this post, for the specific purpose of destroying this canard once and for all. This post was, needless to say, inspired by another individual’s post elsewhere on the rigour of the Earth sciences, and I owe him a debt for inspiring me to get off my backside and post this.

Bear in mind that I’m able to provide this dissertation, not only because I learned the relevant theoretical work, but because I performed actual data measurements of the sort presented below in my physics classes - indeed, one of the eyebrow raising aspects of that time, that I delight in sharing, is that I was permitted to work with real live plutonium in the physics laboratory of my school, at the age of 15.

This latest abridgement of the post contains a few items not specified in older versions. Those familiar with the older versions of this post in various forums, will notice the few but relevant changes where they occur. Also, this forum limits posts to 32,000 characters, so the exposition presented here will be spread across two posts.

Given the nature of the material being covered, this will be a long post. Consequently, it will probably be spread over more than one post. However, I promise that persevering with it will be worth the effort.

Radionuclide Dating Is Rigorous

In order to address this topic at the proper level of detail, something that creationists prefer to avoid at all costs, I shall first begin with a discourse on the underlying physics of radionuclide decay, the precise mathematical law that this process obeys, and how that law is derived, both empirically and theoretically. Note that the decay law was first derived empirically, courtesy of a large body of work by scientists such as Henri Becquerel, Marie Curie, Ernest Rutherford and others. Indeed, the SI unit of activity was named the Becquerel in recognition of that scientist’s contribution to the early days of the study of radionuclide decay, and 1 Bq equals one transformation (decay) per second within a sample of radionuclide.

However, the underlying physics had to wait until the advent of detailed and rigorous quantum theories before it could be elucidated, and is based upon the fact that the nuclei of radionuclides are in an excited state with respect to the sum total of the quantum energy states of the constituent particles (which, being fermions, obey Fermi-Dirac statistics, and consequently, Pauli’s exclusion principle applies). In order for the system to move to a lower energy state, and settle upon a stable set of quantum numbers, various transformations need to take place, and these transformations result in the nucleus undergoing specific and well-defined structural changes, involving the emission of one or more particles. As well as the most familiar modes of decay, namely α and β^- decay, other decay modes exist, and a full treatment of the various decay modes possible, along with the underlying quantum physics, is beyond the scope of this exposition, as it requires a detailed understanding of the behaviour of the appropriate quantum operators, and as a corollary, a detailed understanding of the behaviour of Hilbert spaces, a level of knowledge that is, sadly, not widespread. With this limitation in mind, however, it is still possible to deduce a number of salient facts about radionuclide decay, which I shall now present.

Empirical Determination Of The Decay Law

Initially, the determination of the decay law was performed empirically, by observing the decay of various radionuclides in the laboratory, taking measurements of the number of decay events, and plotting these graphically, with time along the x-axis, and counts along the y-axis. Upon performing this task, the data for many radionuclides is seen to lie upon a curve, and determination of the nature of that curve requires a little mathematical understanding.

To determine the nature of a curve, various transformations can be performed upon the data, which when applied, produce plots that allow linear regression techniques to be applied to the transformed data, to allow the nature of the curve function to be determined. This standard technique was devised by mathematicians back in the 19th century - the first uses of this technique can be traced back to Legendre in 1805, and Gauss in 1809, followed by the work of Quetelet expanding the use of linear regression to the social sciences. Other mathematicians integrated logarithmic transformations with linear regression, as illustrated below, to expand the remit of the technique to elucidating non-linear relationships between data.

The result of each of these transformations is as follows:

[1] Plot log_e(y) against x - if the result is a straight line, then the relationship is of the form:

log_e(y) = kx + C (where C is some constant, in particular, the y-intercept of the straight line)

which can be rewritten:

log_e(y) - log_e(C₀) = kx (where C = log_e(C₀))

which rearranges to:

y = C₀e^kx

where C₀ is derived from the y-intercept of the straight line produced by the transformed plotting, and k is the gradient of the transformed line.

[2] Plot y against log_e(x) - if the result is a straight line, then the relationship is of the form:

y = k log_e(x) + C, where k is the gradient of the line, and C is the y-intercept of the straight line.

[3] Plot log_e(y) against log_e(x) - if the result is a straight line, then the relationship is of the form:

log_e(y) = k log_e(x) + C (where C is the y-intercept of the straight line thus produced)

This rearranges to:

log_e(y) - log_e(C₀) = k log_e(x) (where C = log_e(C₀))

Which in turn rearranges to:

log_e(y/C₀) = log_e(x^k)

Which finally gives us the relationship:

y = C₀x^k, where k is the gradient of the straight line produced by the transformed data, and C₀ is derived from the y-intercept of the straight line produced by the transformed data.

The above procedures allow us to determine the nature of the mathematical relationships governing large bodies of real world data, when those bodies of real world data yield curves as raw plots of y against x. By applying the relevant transformations to radonuclide decay data, it was found that transformation [1] transformed the data into a straight line plot (within the limits of experimental error, of course), and consequently, this informed the scientists examining the data that the decay law was of the form:

N = C₀e^kt

where C₀ and k were constants to be determined from the plot, and which were regarded as being dependent upon the particular radionuclide in question.

Now, if we are start with a known amount of radionuclide, and observe it decaying, then each decay event we detect with a Geiger counter represents one nucleus undergoing the requisite decay transformation. Since the process is random, over a long period of time, decaying nuclei will emit α or β particles in all directions with equal frequency, so we don’t need to surround the material with Geiger counters in order to obtain measurements allowing a good first approximation to the decay rate.

Obviously if we’re engaged in precise work, we do set up our experiments to do this, especially with long-lived nuclei, because the decay events for long-lived nuclei are infrequent, and we need to be able to capture as many of them as possible in order to determine the decay rate with precision. Let’s assume that we’re dealing with a relatively short-lived radionuclide which produces a steady stream of decay events at a reasonably fast rate, in which case we can simply point a single Geiger counter at it, and work out what proportion of these events we are actually capturing, because that proportion will be the ratio of the solid angle subtended by your Geiger counter, divided by the solid angle of an entire sphere (this latter value being 4π). When we have computed this ratio (let’s call it R), which will necessarily be a number less than 1 unless we have surrounded your sample with a spherical shell of Geiger counters, we then start collecting count data, say once per second, and plotting that data.

In a modern setup we’d use a computer to collect this mass of data (a facility that wasn’t available to the likes of Henri Becquerel, Röntgen and the Curies when they were engaged in their work), in order to have as large a body of data as possible to work with. Before working with the raw data, we transform it by taking each of the data points and dividing it by R to obtain the true count.

Once the data has been collected, transformed and plotted, the end result should be a nice curve. At this point, we’re interested in knowing what sort of curve we have, and there are two ways we can determine this. One way is to take the transformed data set, comprising count values c₁, c₂, c₃, … , c_n, where n is the number of data points collected, compute the following values:

r₁ = c₂ - c₁
r₂ = c₃ - c₂
r₃ = c₄ - c₃
…
r_n-1 = c_n - c_n-1

and then plot a graph with r_k on the vertical axis, and c_k on the horizontal axis. This should give a reasonable approximation to a straight line, and the slope of that straight line, obtained via regression analysis, will give the first approximation to the decay constant k. At this point, we know we are dealing with a relationship of the form dN/dt = -kN, and you can then apply the integral calculus to that equation (see below). Technically, what we are doing here is approximating the derivative by computing first differences, another standard mathematical technique from the discipline of numerical analysis - textbooks on this topic are numerous, I might add here.

However, as a double check, we can also perform a logarithmic regression on the data, plotting log_e(c_k) against time, which should also reveal a straight line, and again, the slope of that line will give you the value of k, which should be in good agreement with the value obtained earlier using the more laborious plot of r_k against c_k. In other words, applying the transformation [1] above to the data set, and extracting an exponential relationship from the data. Since we now know that the data is of the form:

log_eN = -kt

we can then derive the exponential form and check that it tallies with the integral calculus result.

Once we have that function coupling the decay rate to time, we can then work backwards, and feed in the values of the known starting mass and the experimentally obtained decay constant k, and see if the function obtained reproduces the transformed data points. If the result agrees with observation to a very good fit, we’re home and dry.

This is, essentially, how the process was done when the decay law was first derived - lots of data points were collected from observation of real radionuclide decay, and the above processes applied to that data, to derive the exponential decay law. When this was done for multiple radionuclides, it was found that they all obeyed the same basic law, namely:

N = N₀e^-kt

where N₀ is your initial amount of radionuclide, N is the amount remaining after time t, and k is the decay constant for the specific radionuclide.

Now, having determined this decay law empirically, it’s time to fire up some calculus, and develop a theoretical derivation of the decay law. Which I shall now proceed to do.

Theoretical Derivation Of The Decay Law And Comparison With The Above Empirical Result

Upon noting, using the calculation of first differences in the empirical determination above, that the rate of change of material with time, plotted against the material remaining, is constant, this immediately leads us to conclude that the decay law is governed by a differential equation. An appropriate differential equation is therefore:

dN/dt = -kN

which states that the amount of material undergoing decay is a linear function of the amount of material present (and furthermore, the minus sign indicates that the process results in a reduction of material remaining). Rearranging this differential equation, we have:

dN/N = -k dt

Integrating this, we have:

∫dN/N = - ∫ k dt

Our limits of integration are, for the left hand integral, the initial amount at t=0, which we call N₀, and the amount remaining after time t, which we call N_t. Our limits of integration for the right hand integral are t=0 and t=t_p, the present time.

Thus, we end up with:

log_eN -log_eN₀ = -kt_p

By an elementary theorem of logarithms, this becomes:

log_e(N/N₀) = -kt_p

Therefore, exponentiating both sides, we have:

N/N₀ = e^-kt

or, the final form:

N = N₀e^-kt
The half-life of a radionuclide is defined as the amount of time required for half the initial amount of material to decay, and is called T_½. Therefore, feeding this into the equation for the decay law,

½N₀ = N₀e^-kt

Cancelling N₀ on both sides, we have:

½ = e^-kt

log_e½ = -kt

By an elementary theorem of logarithms, we have:

log_e2 = kt

Therefore T_½ = log_e2/k

Alternatively, if the half-life is known, but the decay constant k is unknown, then k can be computed by rearranging the above to give:

k = log_e2/T_½

Which allows us to move seamlessly from one system of constants (half-lives) to another (decay constants) and back again.

If the initial amount of substance N₀ is known (e.g., we have a fresh sample of radionuclide prepared from a nuclear reactor), and we observe the decay over a time period t, then measure the amount of substance remaining, we can determine the decay constant empirically as follows:

N = N₀e^-kt

N/N₀ = e^-kt

log_e(N/N₀) = -kt

Therefore:

(1/t) log_e(N₀/N) = k

On the left hand side, the initial amount N₀, the remaining amount N and the elapsed time t are all known, therefore k can be computed using the empirically observed data.

Once again, this agrees with the empirical data from which the law was derived in the earlier exposition above, and consequently, we can be confident that we have alighted upon a correct result.

Once we have the decay law in place, it simply remains for appropriate values of k to be determined, which will be unique to each radionuclide. This work has been performed by scientists, and as a result of decades of intense labour in this vein in physics laboratories around the world, vast bodies of radionuclide data are now available.

Kaye & Laby’s Tables of Physical & Chemical Constants, devised and maintained by the National Physical Laboratory in the UK, contains among the voluminous sets of data produced by the precise laboratory work of various scientists a complete table of the nuclides, which due to its huge size, is split into sections to make it more manageable, in which data such as half-life, major emissions, emission energies and other useful data are included. The sections are:

[1] Hydrogen to Flourine (H¹ to F²⁴)

[2] Neon to Potassium (Ne¹⁷ to K⁵⁴)

[3] Calcium to Copper (Ca³⁵ to Cu⁷⁵)

[4] Zinc to Yttrium (Zn⁵⁷ to Y¹⁰¹)

[5] Zirconium to Indium (Zr⁸¹ to In¹³³)

[6] Tin to Praesodymium (Sn¹⁰³ to Pr¹⁵⁴)

[7] Neodymium to Thulium (Nd¹²⁹ to Tm¹⁷⁷)

[8] Ytterbium to gold (Yb¹⁵¹ to Au²⁰⁴)

[9] Mercury to Actinium (Hg¹⁷⁵ to Ac²³³)

[10] Thorium to Einsteinium (Th²¹² to Es²⁵⁶)

[11] Fermium to Oganesson (names for elements 112 onwards not officially recognised by IUPAC at the time of publication) (Fm²⁴² to Og²⁹⁴)

Now, the above exhaustively compiled data gives rise to yet more data, in the form of the tables covering the major decay series. These arise from the observation of which radionuclides decay into which other radionuclides (or in the case of certain radionuclides, which stable elements are formed after decay), and all of these decay events follow specific rules, according to whether α decay, β^- decay, or one of the other possible decay modes for certain interesting radionuclides, takes place. Again, data is supplied in the above tables with respect to all of this.

Now, we come to the question of how this data is pressed into service. Since the above work couples radionuclide decay to time, via a precise mathematical law, we can use this data to provide information on the age of any material that contains radionuclides. This can be performed by performing precise quantitative measurements of parent radionuclides and daughter products, all of which is well within the remit of inorganic chemists (since the chemistry of the relevant elements has been studied in detail, in some cases for over 200 years) and of course, modern gas chromatograph mass spectrometry can be brought to bear upon the process, yielding results with an accuracy that past chemists reliant upon earlier techniques could only dream of. Consequently, it is now time to cover the business of dating itself see the following post…

[Calilasseia]

Welcome to Part 2.

Radionuclide Dating - The Basics

With the data obtained above, it becomes possible to trace the decay of suitably long-lived elements in geological strata, locate specific isotopes, determine by precise quantitative analysis the amounts present in a given sample, and compare these with calculations for known decay observations in the laboratory, whence the time taken for the observed isotope composition of the sample can be determined. Given that several isotopes have extremely long half-lives, for example, U²³⁸ has a half-life of 4,500,000,000 years, and Th²³² has a half-life of 14,050,000,000 years, and several of the daughter isotopes also have usefully long half-lives, one can determine the age of a rock sample, where multiple isotopes are present, by relating them to the correct decay series and utilising the observed empirically determined half-lives of laboratory samples to determine the age of a particular rock sample, cross correlating using multiple isotopes where these are present and enable such cross correlation to be performed. Thus, errors can be eliminated in age determinations by the use of multiple decay series and the presence of multiple long-lived isotopes - any errors arising in one series will yield a figure different from that in another series, and the calculations can thus be cross-checked to ensure that they are consilient.

Referring to the data tables above, I have selected a number of isotopes of interest. These are isotopes whose half-lives have been determined to lie within a specific range, and which moreover are not known to be produced in the Earth’s crust by any major synthesis processes (except for the various Technetium isotopes, which can arise if Molybdenum isotopes are coincident with Uranium isotopes in certain rocks, but this exception is rare and well documented). The isotopes in question, in increasing atomic mass order, are:

Al²⁶ : 740,000 years
Cl³⁶ : 301,000 years
Ca⁴¹ : 103,000 years
Mn⁵³ : 3,740,000 years
Fe⁶⁰ : 1,500,000 years
Kr⁸¹ : 213,000 years
Zr⁹³ : 1,530,000 years
Nb⁹² : 34,700,000 years
Tc⁹⁷ : 2,600,000 years
Tc⁹⁸ : 4,200,000 years
Tc⁹⁹ : 211,000 years
Pd¹⁰⁷ : 6,500,000 years
Sn¹²⁶ : 100,000 years
I¹²⁹ : 15,700,000 years
Cs¹³⁵ : 2,300,000 years
Sm¹⁴⁶ : 103,000,000 years
Gd¹⁵⁰ : 1,790,000 years
Dy¹⁵⁴ : 3,000,000 years
Hf¹⁸²: 9,000,000 years
Re^186m : 200,000 years
Pb²⁰⁵ : 15,200,000 years
Bi²⁰⁸ : 368,000 years
Bi^210m : 3,040,000 years
Np²³⁶ : 154,000 years
Np²³⁷ : 2,140,000 years
Pu²⁴² : 373,300 years
Pu²⁴⁴ : 81,700,000 years
Cm²⁴⁷ : 15,600,000 years
Cm²⁴⁸ : 340,000 years

The reason I have chosen these isotopes is very simple. Namely, that they would all be present in measurable quantities in the Earth’s crust, and detectable by modern mass spectrometry among other techniques, if the planet was, say, only 6,000 years old, as various enthusiasts for mythology continue to assert. This is because because the half-lives of all these radionuclides are a good deal longer than 6,000 years. So, what do we find when we search for these isotopes in Earth rocks?

NONE of them are present in measurable quantities.

Now, one can safely assume that at the end of 20 half-lives, any measurable amount of a particular radionuclide has effectively vanished - the amount left is ½²⁰, or just 0.000095367% of the original mass that was present originally. So even for isotopes of common elements, this fraction represents a vanishingly small amount of material that would test even the world’s best mass spectrometer labs to detect in a sample. So, what does the observation of no measurable quantity of the above isotopes mean? It means that at least 20 half-lives of the requisite isotopes must have elapsed for those isotopes to disappear. Taking each isotope in turn, this means that:

[1] Sn¹²⁶, being absent, must have disappeared over a period of 20 half lives = 20 × 100,000 years = 2,000,000 years. Therefore the Earth must be at least 2,000,000 years old for all the Sn¹²⁶ to have disappeared.

[2] Ca⁴¹, being absent, must have disappeared over a period of 20 half lives = 20 × 103,000 years = 2,060,000 years. Therefore the Earth must be at least 2,060,000 years old for all the Ca⁴¹ to have disappeared.

[3] Np²³⁶, being absent, must have disappeared over a period of 20 half lives = 20 × 154,000 years = 3,080,000 years. Therefore the Earth must be at least 3,080,000 years old for all the Np²³⁶ to have disappeared.

[4] Re^186m, being absent, must have disappeared over a period of 20 half lives = 20 × 200,000 years = 4,000,000 years. Therefore the Earth must be at least 4,000,000 years old for all the Re^186m to have disappeared.

[5] Tc⁹⁹, being absent, must have disappeared over a period of 20 half lives = 20 × 211,000 years = 4,220,000 years. Therefore the Earth must be at least 4,220,000 years old for all the Tc⁹⁹ to have disappeared.

[6] Kr⁸¹, being absent, must have disappeared over a period of 20 half lives = 20 × 213,000 years = 4,260,000 years. Therefore the Earth must be at least 4,260,000 years old for all the Kr⁸¹ to have disappeared.

[7] Cl³⁶, being absent, must have disappeared over a period of 20 half lives = 20 × 301,000 years = 6,020,000 years. Therefore the Earth must be at least 6,020,000 years old for all the Cl³⁶ to have disappeared.

[8] Cm²⁴⁸, being absent, must have disappeared over a period of 20 half lives = 20 × 340,000 years = 6,800,000 years. Therefore the Earth must be at least 6,800,000 years old for all the Cm²⁴⁸ to have disappeared.

[9] Bi²⁰⁸, being absent, must have disappeared over a period of 20 half lives = 20 × 368,000 years = 7,360,000 years. Therefore the Earth must be at least 7,360,000 years old for all the Bi²⁰⁸ to have disappeared.

[10] Pu²⁴², being absent, must have disappeared over a period of 20 half lives = 20 × 373,000 years = 7,460,000 years. Therefore the Earth must be at least 7,460,000 years old for all the Pu²⁴² to have disappeared.

[11] Al²⁶, being absent, must have disappeared over a period of 20 half lives = 20 × 740,000 years = 14,800,000 years. Therefore the Earth must be at least 14,800,000 years old for all the Al²⁶ to have disappeared.

[12] Fe⁶⁰, being absent, must have disappeared over a period of 20 half lives = 20 × 1,500,000 years = 30,000,000 years. Therefore the Earth must be at least 30,000,000 years old for all the Fe⁶⁰ to have disappeared.

[13] Zr⁹³, being absent, must have disappeared over a period of 20 half lives = 20 × 1,530,000 years = 30,600,000 years. Therefore the Earth must be at least 30,600,000 years old for all the Zr⁹³ to have disappeared.

[14] Gd¹⁵⁰, being absent, must have disappeared over a period of 20 half lives = 20 × 1,790,000 years = 35,800,000 years. Therefore the Earth must be at least 35,800,000 years old for all the Gd¹⁵⁰ to have disappeared.

[15] Np²³⁷, being absent, must have disappeared over a period of 20 half lives = 20 × 2,140,000 years = 42,400,000 years. Therefore the Earth must be at least 42,400,000 years old for all the Np²³⁷ to have disappeared.

[16] Cs¹³⁵, being absent, must have disappeared over a period of 20 half lives = 20 × 2,300,000 years = 46,000,000 years. Therefore the Earth must be at least 46,000,000 years old for all the Cs¹³⁵ to have disappeared.

[17] Tc⁹⁷, being absent, must have disappeared over a period of 20 half lives = 20 × 2,600,000 years = 52,000,000 years. Therefore the Earth must be at least 52,000,000 years old for all the Tc⁹⁷ to have disappeared.

[18] Dy¹⁵⁴, being absent, must have disappeared over a period of 20 half lives = 20 × 3,000,000 years = 60,000,000 years. Therefore the Earth must be at least 60,000,000 years old for all the Dy¹⁵⁴ to have disappeared.

[19] Bi^210m, being absent, must have disappeared over a period of 20 half lives = 20 × 3,040,000 years = 60,800,000 years. Therefore the Earth must be at least 60,800,000 years old for all the Bi^210m to have disappeared.

[20] Mn⁵³, being absent, must have disappeared over a period of 20 half lives = 20 × 3,740,000 years = 74,800,000 years. Therefore the Earth must be at least 74,800,000 years old for all the Mn⁵³ to have disappeared.

[21] Tc⁹⁸, being absent, must have disappeared over a period of 20 half lives = 20 × 4,200,000 years = 84,000,000 years. Therefore the Earth must be at least 84,000,000 years old for all the Tc⁹⁸ to have disappeared.

[22] Pd¹⁰⁷, being absent, must have disappeared over a period of 20 half lives = 20 × 6,500,000 years = 130,000,000 years. Therefore the Earth must be at least 130,000,000 years old for all the Pd¹⁰⁷ to have disappeared.

[23] Hf¹⁸², being absent, must have disappeared over a period of 20 half lives = 20 × 9,000,000 years = 180,000,000 years. Therefore the Earth must be at least 180,000,000 years old for all the Hf¹⁸² to have disappeared.

[24] Pb²⁰⁵, being absent, must have disappeared over a period of 20 half lives = 20 × 15,200,000 years = 304,000,000 years. Therefore the Earth must be at least 304,000,000 years old for all the Pb²⁰⁵ to have disappeared.

[25] Cm²⁴⁷, being absent, must have disappeared over a period of 20 half lives = 20 × 15,600,000 years = 312,000,000 years. Therefore the Earth must be at least 312,000,000 years old for all the Cm²⁴⁷ to have disappeared.

[26] I¹²⁹, being absent, must have disappeared over a period of 20 half lives = 20 × 15,700,000 years = 314,000,000 years. Therefore the Earth must be at least 314,000,000 years old for all the I¹²⁹ to have disappeared.

[27] Nb⁹², being absent, must have disappeared over a period of 20 half lives = 20 × 34,700,000 years = 694,000,000 years. Therefore the Earth must be at least 694,000,000 years old for all the Nb⁹² to have disappeared.

[28] Pu²⁴⁴, being absent, must have disappeared over a period of 20 half lives = 20 × 81,700,000 years = 1,634,000,000 years. Therefore the Earth must be at least 1,634,000,000 years old for all the Pu²⁴⁴ to have disappeared.

[29] Sm¹⁴⁶, being absent, must have disappeared over a period of 20 half lives = 20 × 103,000,000 years = 2,060,000,000 years. Therefore the Earth must be at least 2,060,000,000 years old for all the Sm¹⁴⁶ to have disappeared.

This is an inescapable conclusion from observational reality, given that these isotopes are not found in measurable quantities in the Earth and would be found in measurable quantities if the Earth was only 6,000 years old, indeed, hardly any of the Sm¹⁴⁶ would have disappeared in just 6,000 years, and it would form a significant measurable percentage of the naturally occurring Samarium that is present in crustal rocks. The fact that NO Sm¹⁴⁶ is found places a minimum limit on the age of the earth of 2,060,000,000 years - over two billion years - and of course, dating using other isotopes with longer half lives that can be measured precisely has established that the age of the Earth is approximately 4.5 billion years.

Now since the decay of these isotopes obeys a precise mathematical law as derived above, and this law has been established through decades of observation of material of known starting composition originating from nuclear reactors specifically for the purpose of determining precise half-lives, which is one of the tasks that the UK National Physical Laboratory (whose data I cited above) performs on a continuous basis in order to maintain scientific databases, the provenance of all of this is beyond question. The tables I have linked to above are the result of something like half a century of continuous work establishing half-lives for hundreds upon hundreds of radionuclides, and not ONE of them has EVER been observed to violate that precise mathematical law which I opened this post with under the kind of conditions in which those materials would exist on Earth if they were present.

The majority of those isotopes are nowadays ONLY obtained by synthesis within nuclear reactors, and observation of known samples of these materials confirms again and again that not only does the precise mathematical law governing radionuclide decay apply universally to all of these isotopes, but that the half-lives obtained are valid as a consequence. The laws of nuclear physics would have to be rewritten wholesale for any other scenario to be even remotely valid, and that rewriting of the laws of nuclear physics would impact upon the very existence of stable isotopes including stable isotopes of the elements that make up each and every one of us, none of which would exist if the various wacky scenarios vomited forth on creationist websites to try and escape this were ever a reality.

Now, having covered this, it’s time to deal with a topic of importance, and one that constitutes another of those epistemological bones creationists love to chew upon, without realising that the questions surrounding this topic have already been answered. That question is …

Transport - And How Scientists Deal With It

Even in the very early days of the development of radionuclide dating techniques, scientists were aware that chemical elements in rock samples might be subject to various chemical processes resulting in the transport of material in or out of a given rock sample. After all, inorganic chemists had been performing detailed investigations of the aqueous chemistry of geologically significant elements for the best part of a century before the development of radionuclide dating. Scientists interested in utilising radionuclides for dating purposes, therefore had a vast body of knowledge from inorganic chemistry to draw upon, in order to determine whether or not transport was likely to take place in a given rock stratum, and given the known behaviour of the relevant elements and their various salts in aqueous conditions, they were able to devise means of testing whether transport was likely to be a significant factor in a given rock stratum of interest. Which they promptly set about doing.

For example, a geologist can reasonably conclude that virtually zero transport has taken place in a sample where the radionuclides are found as sulphide salts, because the only sulphides that are soluble in water are those of the alkali metals - ALL other sulphides are insoluble, with solubility constants that are extremely low. Lead (II) sulphide, for example (known geologically as the mineral galena), has a K_sp of 3.4 × 10^-28, mercury (II) sulphide (cinnabar) has a K_sp of around 4 × 10^-53, and silver sulphide (acanthite or argentite) has a K_sp of 1.6 × 10^-49.

Indeed, thanks to the large body of knowledge bestowed upon geologist by inorganic chemists, all manner of tests can be performed in order to determine if transport is likely to be a problem before a sample is submitted for dating analysis, so that this can be taken into account and proper corrections applied to the material in question. Indeed, scientists have spent time devising a technique, known as isochron dating, whose purpose is specifically that of checking whether the isotopes used for dating have been subject to transport, and developing appropriate corrections to the dates that would have been obtained without such checking. This technique takes advantage of the fact that many radionuclides produce a stable decay product that is isotopically different from a typical sample of the daughter element arising non-radiogenically, and appropriate comparisons of the amounts of non-radiogenic daughter element with that of the radiogenic daughter element can be used to determine if transport took place.

Let’s see how isochron dating works, shall we?

When an igneous stratum is formed (and dating techniques tend to focus upon igneous strata, because the deep mantle of the Earth is the primary source of the radionuclides required), molten material solidifies, and as it does so, the chemical compounds of various elements crystallise out of the melt differentially, resulting in the formation of well-defined mineral species. Moreover, some melts possess a greater chemical affinity for some elements than others, and consequently, there will be further differentiation, based upon those chemical affinities. Zircons, for example, have a chemical affinity for uranium, and incorporate uranium salts into their crystals far more readily than, say, lead salts. Likewise, some other crystalline formations will incorporate other elements preferentially, including elements that are of utility value in dating.

A typical element that is of use is rubidium, with the Rb⁸⁷ isotope having a half-life of 4.75 × 10¹⁰ years. This decays via β^- decay into Sr⁸⁷, an isotope of strontium, and this decay product is stable. However, a non-radiogenic isotope of strontium, namely Sr⁸⁶, also exists in geological strata. The utility of this will become apparent shortly.

Now, any melt starting off with a quantity of Rb⁸⁷ will, after sufficient time has elapsed, start to acquire quantities of Sr⁸⁷. This will take time, given the long half-life of Rb⁸⁷, and indeed, detectable quantities of Sr⁸⁷ amenable to mass spectrometry will only start to appear after 0.001 half-lives have elapsed. In the case of Rb⁸⁷, this is 45 million years, so a sample that has an Rb⁸⁷ age that is indistinguishable from zero could be as much as 40 million years old, which means that this is no indication of an allegedly “young” Earth. The moment any detectable traces of Sr⁸⁷ appear in a sample, however, then we’re dealing with a rock that is at least 45 million years old, which flushes blind assertions about the Earth being only 6,000 years old down the toilet to begin with.

But, this isn’t the point: the point is, that any mineral that acquires quantities of strontium upon formation will acquire a specific ratio of the two isotopes Sr⁸⁶ and Sr⁸⁷, and that ratio can be used to determine the initial amount of Rb⁸⁷ that was acquired during formation as well. Which means that no “assumptions” about initial material present are needed.

Now, since all the minerals that acquire strontium will acquire the same ratio of Sr⁸⁶ to Sr⁸⁷ at the start, we can use that ratio as the y-axis for a plot. However, different minerals within the sample will acquire different quantities of Rb⁸⁷, and we can use the ratio of Rb⁸⁷ to Sr⁸⁶ to form the x-axis of the plot. If our initial values are:

P = amount of Rb⁸⁷ at the time of stratum formation
D₁ = amount of Sr⁸⁶ at the time of stratum formation
D₂ = amount of Sr⁸⁷ at the time of stratum formation

then we plot an x-y plot consisting of:

x = P/D₁, y = D₂/D₁

the axes therefore correspond to:

x = increasing enrichment of Rb⁸⁷ in the sample with increasing value
y = increasing enrichment of Sr⁸⁷ in the sample with increasing value

Now, the global composition of the melt, from which the stratum eventually forms, will have a given point value on this plot. As the melt cools, and minerals crystallise out, different minerals will migrate along a straight line in this plot, as all the minerals will inherit the same value of y (= D₂/D₁), but inherit different values of x (=P/D₁). The result will be, at zero age, a horizontal line connecting the points for those minerals in the plot.

Now, as the Rb⁸⁷ decays in each mineral, it will produce Sr⁸⁷. Therefore, over geological time, the data points will move upwards and to the left.

Now, because decay occurs in a proportional manner, courtesy of the decay law derived above, when 20% of the Rb⁸⁷ has decayed in one mineral, then 20% of the Rb⁸⁷ will have decayed in all the minerals present in the stratum. This means that all the points will move upwards and to the left of the plot. This means that those points corresponding to the minerals with the greatest initial Rb⁸⁷ concentration will move the farthest, and the movement will be such that all of the points will remain on a straight line. Indeed, those familiar with Monte Carlo simulation methods can produce a computer simulation that reproduces this exact result. As a consequence, any stratum containing the appropriate minerals will yield, for a given age, a line of points whose slope increases with time, and the slope of that line can be used to determine the age of the sample.

Note that we don’t need to know the initial amounts of any of the elements present in the sample in order for this to work. All we need to know is the present-day ratios of those elements. And, from that data, we can reconstruct the original composition of the melt.

Now, here is the fun part. If the sample suffered any transport, then this will become immediately apparent because the points will deviate significantly from a straight line plot. This is because the minerals containing the various elements will undergo transport differentially. For example, if the minerals lost Rb⁸⁷ due to various processes of aqueous chemistry, then this would shift all of the points to the left, but because different minerals are involved in the process, each with a different chemistry, and a different response to the aqueous processes that could remove Rb⁸⁷ from them, the amount of x-shift will be different for each mineral, and as a consequence, straight line correlation will be destroyed. Immediately scientists see this, they know that something unusual has happened to the sample after formation, and that dates obtained from it are unlikely to be reliable.

Indeed, if the Earth were only 6,000 years old, not only would ALL isochron plots involving minerals with long half-lives be flat, horizontal lines, but random loss or gain of the parent nuclide would not affect those results. If the Earth were only 6,000 years old, any “contamination” that creationists assert would be present that would render radionuclide dating methods invalid, would not even be noticeable on a purported “young Earth”. Contamination would only ever be noticeable if the Earth was old. Which means that far from “invalidating” radionuclide dating, detectable contamination in a sample via the isochron method actually reinforces an old age for the Earth.

Worse still, from the creationist standpoint, the known possibilities for systematic alteration of an isochron plot, that results in modification of the data points such that they remain collinear, result in an underestimate of the age of the sample - in other words, the isochron plot records them as being younger than they actually are, not older. Complete homogenisation of the stratum with respect to the content of Sr⁸⁷ will reset the isochron ‘clock’ to zero, and partial homogenisation will result in a line of lower slope than previously, causing the technique to underestimate the true age.

Once such a plot is complete, any strong correlation to a straight line (as yielded by regression analysis, which will give us an indication how much we can trust a date from the plot, and what value of error to apply to the plot), will yield the following values:

[1] the y-intercept of the line is the value of D₂/D₁ at solidification;

[2] the slope of the line is equal to ΔP/(P - ΔP), where ΔP is the amount of Rb⁸⁷ lost to decay, and thus allows us to determine the age of the sample, given the known half-life of Rb⁸⁷ (we can back-calculate using the decay law to determine what value of P was initially present to produce the observed result).

There is a body of technical literature on the subject of isochron dating, and, needless to say, it requires a fair amount of work to plough through, but for those who have acquired the relevant scientific and mathematical background, it is well worth exerting this effort, not least in order to prepare oneself for the inevitable quote mining of the scientific papers that will arise in creationist circles.

And now, I’d like to turn attention to this:

Purported “Exceptions” To The Decay Law

One of the more interesting findings of recent years was that rhenium-187, a radionuclide that undergoes decay by electron capture, can have its half-life altered under laboratory conditions. However, the mechanism of alteration offers no hope to anyone entertaining fantasies about a 6,000 year old Earth, for the following reasons:

[1] The mechanism only applies to isotopes that decay via electron capture, and both α decaying and β^- decaying isotopes are unaffected by this mechanism;

[2] The mechanism only has a small effect upon the half-life of Re¹⁸⁷ unless extreme conditions are applied (and the scientists had to resort to extreme measures to determine this);

[3] The conditions required to produce this mechanism in Re¹⁸⁷ atoms on Earth are physically unreal.

What the relevant scientists did, was to take Re¹⁸⁷ atoms, and subject them to ever more extreme levels of ionisation. In order to do this, they had to resort to a particle accelerator, and use large amounts of energy to ionise the Re¹⁸⁷ atoms to the required amount. As any student of basic physical chemistry understands, it takes a certain amount of energy to remove one electron from a neutral atom, but, once that electron is removed, it requires more energy to remove a second electron, because that electron is now being removed from a positively charged ion, and more energy is needed to overcome the additional electrostatic attraction now present. Consequently, as one removes yet more electrons, more energy is needed for each subsequent electron.

In order to speed up the decay process of Re¹⁸⁷, scientists stripped all of the electrons from the Re¹⁸⁷ atoms, making them extremely highly ionised. Extreme states of ionisation of this sort only occur in nature within the confines of extremely hot incandescent plasmas, with temperatures well beyond that required to turn elements such as rhenium into gases, and rhenium has the highest boiling point of all the elements, at a whopping 5,600°C (courtesy of Kaye & Laby again). Even this temperature isn’t enough to strip 75 electrons off a rhenium atom, and the temperature required is more likely to be of the order of 65,000 Kelvins or beyond, which means that trying to invoke this mechanism as a means of “accelerating” nuclear decay involves heating the Earth’s crust to temperatures more normally associated with the chromospheres of supergiant O class stars.

Apart from the fact that this mechanism requires ludicrously absurd conditions to have occurred in the Earth’s crust in order for it to happen, and apart from the fact that this mechanism is useless for U²³⁸ and several other important isotopes used in radionuclide dating, because they do not decay via electron capture, it’s actually a waste of time trying to argue against the constancy of the decay law under normal conditions, because this one exception requires such extreme conditions that the Earth would not have remained a solid planet if they had been present.

Plus, the mere fact that the half-life for Re¹⁸⁷ is of the order of 4.16 × 10¹⁰ years, means that the presence of any detectable Re¹⁸⁷ in the Earth’s crust means that the planet is ancient. If by some bizarre freak of physics, Re¹⁸⁷ decay had been sped up by this mechanism without vaporising the Earth’s crust, the anomaly would have been so immediately detectable alongside the other radionuclides, that scientists would have noticed something was amiss long ago. That word ‘consilience’ is back to haunt those who prefer mythology to reality.

I think this covers relevant bases for now.

[TheMagus]

Very interesting! I find it funny when I ask my Christian friends how do they measure that the earth is 6000 years old and they can’t tell me. They’ve thrown out every dating tool we have because it disagrees with their static model. Never mind that there is a 9000 year old tree.

[dogalmighty]

Hey Cal…once you are down to quantifying decay rates, the average religitard is going to be quoting their respected comic books verses.

[Whitefire13]

@Calilasseia

fact that this mechanism is useless for U238. I

@Calilasseia

the mere fact that the half-life for Re187 is of the order of 4.16 × 1010 years, means that the presence of any detectable Re187 in the Earth’s crust means that the planet is ancient.

What is this strange language being spoken here??? This is above my brain-grade. Grabs black book

Flips page.

HEY butterfly - “In the beginning God created the heavens and the earth…”

Now this- THIS I understand.

walking away with belief and faith still in tact

[algebe]

@Calilasseia

Nice summary of a complex but very useful science.

If creationitwits can’t understand radioactive decay, maybe they can wrap their limited intellects around trees and snow. Both dendrochronology and ice cores go back substantially further than the creationist fantasy of a 6,000yo Earth.