How do calculate the energy release of an A-bomb from first principles (eg, Trinity)

I was thinking about the infinitismal, but real energy created/released (I don’t think I really understand the difference) when, on the rare occasion, a neutrino zaps into my retina. (Neutrinos just came up in a thread; neutrinos zapping retinas is discussed at length in an old thread.)

So then: you can add up those events, and after a while they add up to real money, as they say. Is that how the physicists at Los Alamos did it, but with the what-all atomic events going on in that reaction?

At Trinity, they notoriously weren’t sure - there was a betting pool. At one point, they considered the possibility that it might ignite the atmosphere, although I think they were at least fairly sure this couldn’t happen. Fermi had side bets going on the destruction of the entire state.

http://www.techrepublic.com/blog/geekend/geek-trivia-which-manhattan-project-scientist-won-the-groups-infamous-doomsday-betting-pool/

I think after a few basic empirical factors were sorted out, simple A-bombs could then be simulated accurately.

But they then messed up on the predicted yield of one of the early H-bomb tests, Castle Bravo, missing one reaction path completely.

How did he plan to collect his winnings?

I seem to recall that one of them scattered some papers about, so that he could measure how far they blew and calculate the energy release from that.

Also Fermi. Hence “Fermi estimation”.

A bit of a joker and a smart fellow, Enrico.

The amount of energy created by fission can be both calculated and measured. Calculated by the change in mass of the products vs the original atom’s mass, and measured by boiling water in a reactor. The problem during Trinity was the unknown efficiency of the bomb - how much of the mass of the bomb actually underwent fission. So during the Trinity test, they set up many different instruments to measure the force of the blast. During later tests (“Mike,” for example), a rough estimate of the magnitude of the explosion could be determined by the size of the fireball, which one scientist did by using a quarter held at arm’s length.

Before anyone actually made any atomic explosions, they had only very roughly an idea how long criticality would be maintained such that the chain reaction would continue, and how much would be consumed during that time period. They knew that the energy from the explosion would cause itself to blow apart eventually such that it was extremely unlikely that even most of the material that was fissile actually would undergo fission (counting by atom; most likely it involved nuclei in every region of the critical mass, just not all of them), but it’s telling that the yield was “predicted” the best by the last person who entered the betting pool and had to take the last choice (so really he wasn’t right, everyone else was wrong).

Thermonuclear weapons do not produce their additional energy through fusion as much as through additional fissioning of U-238 that requires the fast neutrons from the fusion reactions in order to fission. Still in these designs, they will blow themselves up before consuming a majority of the fuel. It would require building the entire device atom-by-atom to be able to ensure fissioning of almost all the fuel, and the required amount of work/expense to do that is prohibitive compared to just mashing it all together with conventional explosives and letting most of it go to waste. Maybe someday we’ll have nanotech advanced far enough to make a significant explosion that uses all of its fuel, but it won’t be any time soon. (Even then, it probably wouldn’t be good enough because in room temperature stable conditions you’d have to have those pesky electrons to keep everything organized, and their spatial probability distribution suggests that they might be able to mess things up beyond predictability thresholds - but I don’t know enough about the process to make solid claims about this sort of thing)

G. I. Taylor famously published rather accurate estimates of the Trinity test’s yield using only the series of time-stamped still images published by Life magazine (example image at that link). The yield was still classified at the time.

Taylor’s calculation was rather involved, but you can get a pretty good answer using the single image in the Life magazine link above and some dimensional analysis. To wit…

Things we know: At time t=0.025 seconds after detonation the shockwave reaches a radius r=130 meters, give or take. En route it was pushing air out of the way, and air has a density d=1.2 kg/m[sup]3[/sup]. (Instead of the density of air you might find it more natural to think in terms of the total mass of the air displaced, but this mass is simply the air’s density times the volume under that hemispherical shockwave, both of which we already know. So, we can work with density or total mass displaced.)

We want to calculate the energy released by the bomb. Energy has SI units of kg/m[sup]2[/sup]/s[sup]2[/sup], so we need an expression involving t, r, and d that gives us the right units(*). There is only one way to do so:

E = d r[sup]5[/sup] / t[sup]2[/sup] .

Assuming any missing numerical (unitless) factors are of order unity (i.e., not too much bigger or smaller than 1, a common – and commonly good – assumption), we can just run the numbers. Google can do the unit conversions for you to get to “kilotons TNT”:

Google: (1.2 kg/m^3)*(130 m)^5/(0.025 s)^2 in kilotons TNT = 17 ktons TNT.

The actual yield is estimated at 21 kton TNT.
(*) self-nitpick: we need an expression that has the correct “dimensions”. Units are a matter of convention.

Neat. (To say the least.)

But back to the “first principles” type of answer: is that (crudely spelled out in OP) remotely how the thing was estimated beforehand?

I haven’t seen my copy of Rhodes in a dog’s age.

I wasn’t sure what you meant in the OP. Estimating the yield before actually detonating the bomb is primarily an engineering question. The physics of fission reactions, neutron moderation, etc., were known to a decent enough level by the end of the Manhattan project. The driving question was how much of the reactants would actually react before the whole thing was blown to smithereens. That’s a tough engineering question, but it’s philosophically no different from other engineering questions.

Or are you after how the energy release of a single particular nuclear fission reaction is known (never minding how many of those you get when the bomb goes off)? Some of this has been said upthread, but more can be said. I’m just not sure which question you’re asking.

Of course, it’s always ambiguous just what “first principles” actually means, since there are always other principles before the ones you consider “first”. As Sagan once said, to bake an apple pie from scratch, you must first create the Universe.

They almost certainly knew what the absolute upper-bound was- with a known quantity of fissile material, if you assume every atom fissions, you’ll be able to calculate the absolute maximum amount of possible energy released.

But like others have said, knowing just how much of the pit would actually fission is something they couldn’t predict. They could (and probably did) do calculations like “If 1% of the U235 inthe bomb fissions, yield will be X kilotons. IF 2%, then Y kilotons” and so forth. That was relatively easy- they knew how fission worked, as they already could see it in action in various experiments. The only trick was getting it to actually go prompt critical, and there was a whole lot of variables and uncertainty surrounding the bomb itself. For example, if the pit wasn’t perfectly spherical or was not uniformly dense, that would have an effect on the yield. So would imperfect explosive lenses or detonation timing.

For example, Little Boy was only 2% efficient- 2% of the uranium actually fissioned (an amount the size of a peppercorn), while Trinity and Little Boy were more efficient-

After a while, they got a fairly good handle on just how much might actually fission, but at first, they had no idea.

As an example, the “2nd generation” Fat Man-type bombs used a “levitated” core, where there was a gap between the core and the “pusher” (tamper). This resulted in significantly greater compression of the core, and more than doubled the efficiency. (The analogy is pushing a nail into a board with a hammer, verses striking the nail with the hammer).

Right, and they couldn’t really know ahead of time how well that might work, other than knowing that it would be more efficient.

Estimate of the yield of the Nagasaki bomb.
From the Oak Ridge National Laboratory ORNL data, critical mass versus natural uranium reflector thickness ( dr ) for 93% U 235 sphere.
The following is an empirical formula that gives approximately the effective critical mass versus a natural uranium reflector/tamper thickness dr.
Rc(dr) = [ 0.92 +{ 0.48 / (1+ dr/4.25) } ] Rc(dr => 21.6 cm.) (a)
Fat Boy ( Nagasaki )
Uranium tamper 2.5” thick, 6.35 cm, estimate 115 Kg
Aluminium Pusher 4.5”, 11.43 cm, “ 118 Kg
Plutonium 6.2 Kg
Mt approx… 240 Kg
ORNL Plutonium 16.5 gm/cm3 critical bare mass 16.25 Kg, infinite reflector 5.8 Kg.
Rc(bare dr = 0) = 6.17 cm, Rc(dr = infin.) = 4.37 cm. R at 6.2 Kg = 4.47 cm.
From (a) Rc (dr = 6.35 cm) = 4.87 cm and Mc = 8.0 Kg.
ORNL compression factor Kc quoted as “twice ”and “2.5 possibly somewhat less”. The numerical mean is 2.25 and the phrase “possibly somewhat less” would put this figure in the region of 2.2 to 2.1.
Pu fission MFP = 14.2 cm giving t = 8.5 x 10^-9 sec. Neutron yield n = 3.
E = Mt/2 { [R Kc^2/3 – Rc][1- ( Rc/(R Kcom^2/3))^2 ][n – 1 ]/4t }^2 ergs.**
Energy equivalent 4.2×10^16 ergs/ton TNT.
It is estimated that 10% of the energy released was due to the fission of the tamper.
Kc = 2.1 yield 18.58 x 10/9 = 20,64 tons TNT
2.15 21,39 23,76
2.2 24,40 27,11
Dr Penney observed the Nagasaki explosion from an escort plane and estimated from an on the ground survey a yield of 22± 2 Kilo tons TNT.
** Based on a lecture given by Werner Heisenberg 1945.

Correction, that should have been a 93% U238 sphere.

Wouldn’t they have studied the paramaters of neutron bombardment ?
take a sample, measure its radioactive elements, temperature, bombard it, measure it again.

Change bombardment intensity, change thickness, change temperature… change all the variables.
So you can find Kc, to calibrate Genro’s equation … Genro lists the equation to calculate an effective reflector… I think that was a far different question, what size of reflector to use for the core of about 6.2kg, and what yield… is there an optimum relfector size ? The size, Rc appears in the Energy release equation and they went for a sweet spot of 6.2kg of Plutinium ,with a thickness of reflector of 4.47cm, with the equation being based on geometry (^2 and ^2/3 ) , with the rate of reaction having to be calibrated by also measuring Plutonium neutron MFP, mean free path , the mean distance a neutron can travel through Plutonium before a collision , and left there to point out that the neutrons from different sources have different speeds. But it was convenient to leave that as time given a limitted selection of fissionables to use.

Good question. Unlike the Hiroshima bomb which used Tungsten Carbide as a reflector, the Nagasaki bomb used natural Uranium. Tungsten Carbide has a large neutron scatter cross section and little neutron absorption. The optimum thickness can be mathematically defined using hyperbolic functions. In the case of natural uranium, it not only has a scatter cross section but also a small fission and in-elastic scatter cross section which makes it extremely difficult the optimise mathematically. For this reason a great many f experimental measurements were undertaken. These results are tabulated in an ORNL documents and for convenience I have produced a simple equation.

Rc(dr) =[ 0.92 + { 0.48/(1+dr/4.25) } ] Rc(dr=>21.6 cm )

RC (dr=>21.6 cm) is the critical radius for Plutonium with a natural Uranium of reflector of 21.6 cm.

Criticality is the product of density and critical radius. Doubling the density reduces the radius to 0.8 thus making criticality greater than 1. Achievable compression factor were established with the RaLa experiments and exactly what they are is either a top secret, to which I am not privy or nobody really knows. Hence the use of values from 2.1 to 2.15.

Assuming the reflector has the same neutron diffusion coefficient as Plutonium and zero absorption the critical radius using published nuclear parameters yields 4.62 cm (6.6 Kg) and from ad ORNL report, with natural Uranium 4.42 cm (5.8 Kg).

With a natural Uranium reflector of 6.35 cm the critical mass radius is 4.87 cm (8Kg). The actual radius used was 4.47 cm (6.2 Kg) which is about 0.9 of criticality. Compressing this by a factor of 2.15 gives a criticality of 1.54.

Only one Hiroshima bomb was made and unlike the Nagasaki bomb it was simpler and more information is available.

Estimate of the yield of Little Boy the Hiroshima bomb.

U 235 enrichment 80 %, q = 0.8 Mf = 64 Kg.
Tungsten carbide tamper Mr = 311 Kg. The thickness of this gives a critical mass of 17 Kg i.e. 100% U235 .Compression factor K is of course =1.

M 2335 + Mtamper = Mt = 375 Kg

Over the range of q = 0.8 to 1, the critical radius Rc to a first approximation is given by :-

Rc = ( pi Aw) / (2 p Av) { 1/[ 3 (n – 1) Gf* Gs* ]}^1/2.

Aw = atomic mass, Av = Avogadro constant, p = density, n = neutron yield, and Gf*, Gs* are the effective neutron cross-section for fission and scatter respectively.

Gf* = Gf235q + Gf238( 1– q ) and Gs* = Gs235q + Gs238( 1 –q ).
________ Gf_____ Gs.
U235____ 1.24____ 4.56 x10^-24 cm2.
U238____ 0.3_____ 4.8.

Gf* = 1.052×10^-24cm2, Gs* = 4.6 x 10^-24cm2, n = 3, p =18.95 gm/cm3.

Rc =6.63cm , Mc = 23 Kg.

M =64 Kg, R = 9.3 cm, n = 2.64, t = 1 x10-8 sec.

Et = Mt/2 { ( R – Rc )( 1 – (Rc/R)^2)( n – 1)/4t }^2 **.

Et = 5.39 x10^20 ergs and at 4.2 x10^16 ergs/ton TNT :- yield = 12.800 tons TNT.

Dr. Penney estimated from a ground survey a yield of 12±1 Kilo ton TNT.
** Based on a lecture given by Werner Heisenberg 1945.

Thanks to Genro for such a detailed answer. Always good to have good answers here on the Dope.
:slight_smile: