# How did Enrico Fermi compute when the Chicago Pile-1 nuclear reactor would become critical?

I'm trying to understand the first nuclear reactor, the Chicago Pile-1, specifically the math Fermi did to figure out when the reactor would go critical. There's a nice report available from Fermi, where he tracks the value of $$R_{eff}^2/A$$, where $$R_{eff}$$ is the effective radius of the pile, and A is the measured neutron intensity at the center of the pile, some screenshots below.

Fermi then claims that when $$R_{eff}^2/A$$ reaches 0, the pile will become critical. This is where I get lost - I read through the report and some other sources, but I don's see where this math is coming from -> why should the pile become critical when $$R_{eff}^2/A=0$$?

• Clearly, criticality is associated to $A$ becoming divergent, which kind of makes sense. I guess the question really is why he prefers to plot $R_{eff}^2/A$ instead of just $A$, right ?
Feb 14 at 13:14
• At some volume of the pile ($R_{eff}$ estimated from an ellipse as they built up towards the spherical pile end state) (see page 4), neutron losses to the outside world become irrelevant, and you go critical. When you go critical, the neutron activity (A) at the center of the pile goes to infinity. Feb 14 at 13:25
• physics.stackexchange.com/questions/696091/… may be of some interest. Feb 14 at 15:37

This kind of plot is typically called a "1/M Plot", where M is the reactor multiplication. To measure the multiplication, you need a neutron source and a detector. As you add more and more material to the reactor, the detector signal increases. When the reactor goes critical (or more precisely supercritical), the detector signal (M) will grow to infinity. If you plot the inverse of the signal (1/M), it will become zero when the reactor goes critical.

If you plot 1/M as a function of the reactor size (or control rod position, or number of rods, etc.), you can do a curve fit to estimate where the point of criticality is. In Fermi's plots, they were adding material (layers) and the plot shows the reactor when critical at layer 57.

I cannot be sure I understand it correctly, but this is how it looks to me at the moment. The report that you quote says:

In a spherical structure having the reproduction factor 1 for infinite dimensions the activation of a detector placed at the center due to the natural [my emphasis] neutrons is proportional to the square of the radius.

I guess the question really is why he prefers to plot $$R^2_{eff}/A$$ instead of just $$A$$

According to the quote from the report, if you just have $$A\to\infty$$ when the radius increases, it is possible that you just have natural neutrons. If, however, $$R^2_{eff}/A\to 0$$, it means that the number of natural neutrons becomes negligible compared to the number of neutrons from chain reaction, therefore, the reactor is getting critical.

I have made a video response here: https://youtu.be/m-LwNJmicmU

EDIT: in response to comments I have included a typed explanation below, exclusively for stack exchange...

There are 2 key quantities here:

• $$R_{eff}^{2}$$ = square of the effective radius of the reactor.
• A = induced activity in the indium foil from the neutron irradiation (counts/minute), measured the next day.

Unfortunately it's a bit of conjecture from here on, but my understanding is that Fermi used $$R^{2}$$ as a surrogate for the cross sectional area of his reactor, and A as a surrogate for the 'neutron intensity' (presumably, neutron intensity = number of neutrons per unit area).

Fermi was likely most interested in how the value of A changes over time.

• Increasing A = chain reaction (i.e. nuclear bomb)
• Constant A = sustainable reaction that continues until you run out of fissile material.
• Decreasing A = reaction that will create a small bit of heat but die out.

Fermi could have just plotted A against the number of layers and found where his trendline shoots off to infinity, and determined that to be the point of criticality. But A isn't the only influence factor, there's an interplay effect between A and the cross sectional area. Since he's using effective radius ($$R_{eff}$$) it mightn't be obvious how that behaves, and the result will also be very specific to his setup.

Either way we can see that $$R^{2}$$ will increase with the square of the number of layers, and we know A will increase exponentially. Hence A grows faster than $$R^{2}$$, so when we take a ratio of $$R^{2}/A$$ it trends to $$0$$ for criticality. Importantly since both cross sectional area and neutrons/area are strictly positive and non-zero, the ratio will never truly reach zero.

There are a few aspects of the experiment that my (rather naive) video solution cannot speak to. For example, I assumed a perfectly spherical reactor with radius '$$R$$' when in reality it wasn't spherical and Fermi used the effective radius '$$R_{eff}$$' instead. There is also uncertainty surrounding the indium foil detector, I have assumed that the quantity Fermi was actually interested in was something similar to neutron fluence and his measurements were aimed at evaluating it through a proxy-quantity (induced activity in the foil) since he was unable to measure it directly.