# Is there a simple way to compute some physical constant from Feynman diagram statistics?

I've been playing around writing some software to generate Feynman diagrams for QED, respecting the vertex "rules" described here, and avoiding creating isomorphic duplicates.

So from a starter fermion arc

the only possibility - at least the way I create these things, by adding pairs of vertices by splitting fermion arcs and linking the fresh vertices with boson arcs - is:

(boson arcs are in red) and from that there are 3 possibilities

and then 15

and then 105, and then 945 (for 10 interior vertices) and so on (I just killed it after a few hours runtime and it was working it's way through the many thousands of 12 interior vertex possibilities).

My question is: given all these diagrams, is there some simple procedure I can apply and have some fundamental dimensionless constant of the universe (or at least of QED) pop out the other end ?

Background: my interest in this stuff got started with the iconography of these diagrams and their visual aspects (e.g how to layout such things nicely, and how to plot them with some wiggly lines as per the usual convention for the photons will get some attention next), and the original driving idea was to do some nice big poster which gives some visual impression of how the number of possibilities explodes. But it'd be the icing on the cake for the poster to also be able to show that some physical property of the universe follows as a direct consequence of some statistics or whatever of the diagrams (capped with the "Nobody knows how it can be like that" Feynman quote no doubt).

Disclaimer: I know mathematics and computing but most of my knowledge of this area of physics is from reading "pop science" books and I'm well aware they grossly simplify things, so maybe I'm hoping for too much from the diagrams alone.

Update: following helpful input from Vibert below, I'm now generating the below for 2 and 4 vertex diagrams (a grey background indicates "one particle reducible" diagrams); these now include the expected virtual particle loops (missed in my first attempt above). Diagrams which are isomorphic by simple reversal of all the arrows are omitted in the interests of saving space.

Outcome: Well after digging into this a bit more it seems there's no getting away from Doing The Math, integrals and all. I was originally hoping there might be something as simple and understandable as "well there are N diagrams in this category, M diagrams in that category and so the correction to the 1st order approximation is (N-M)/(N+M), which is more accurate". But material on Lamb Shift corrections gives an idea of what's actually involved (e.g the contribution from the simplest diagrams is ~1058 MHz, most of it from one case), while the next order's contributions sum to on the order of 100kHz.

Update 2022: Almost a decade on and I recently found myself revisiting this stuff... and creating a big A0 poster (tikz-feynman rendered, but I wasn't too happy with its layout algorithm and found something in the python-world's graph-tool which worked better):

https://pdfhost.io/v/TKiPIg~hO_Feynman_Diagrams_for_Quantum_Electrodynamics

(2.2MByte PDF; may take a while to display via the apparently notoriously slow pdfhost site). The number of diagrams is - as should be expected - 2, 10, 74, 706 (next would be 8162, which has no chance of fitting!) which is OEIS A000698... maybe landing on such a well studied series is as near as I'll get to the original idea of some "constant" emerging!

Here's an overview of the PDF

(or in higher resolution - although still nowhere near as good as the PDF - here).

What you're showing are diagrams that contribute to the two-point function/propagator of the electron. Essentially these are so-called loop corrections to the electron self-energy, which plays an important role in the renormalisation of QED. If you organise the calculation properly, then these diagrams encode how the physical mass of the electron depends on its interaction with the photon.

There are some subtleties: first of all, you need only a subclass of these diagrams (called 1PI - look the exact meaning up in some textbook, such as Peskin-Schroeder). Next, you really need to calculate the corresponding integrals; the combinatorics alone doesn't tell you much. By hand, you can more or less calculate the diagram with one red line; for the rest, people use computer programs.

If you consider the result $f(\alpha)$ as a function of the coupling constant $\alpha$, it should mathematically be possible to show that $f$ has (formally) zero radius of convergence mainly by counting the number of diagrams, but that's not a very interesting result, physics-wise.

• Thanks; indeed, a quick look into "one particle (ir)reducibility" has made me realize that all my diagrams are all one particle reducible, because I never glue any fermion-antifermion pairs into my boson arcs (this will be rectified promptly!). Alas, it seems I was a bit optimistic hoping to get something quantitative from diagrams alone. Peskin-Schroeder sounds rather intimidating, based on some amazon (uk) reviews flagging it as definitely not a primer (although it certainly looks interesting via "look inside"); any pointers to something more basic welcome. Dec 17, 2012 at 14:11
• Less than half of the ones you drew are reducible! The test is if you can cut the diagram in two pieces by making a single cut (not on the 'legs' which run in and out of the diagram). I understand that the reference is quite difficult, but you might get a flavor of what these diagrams really mean. Another relatively 'easy' book is A. Zee - QFT in a Nutshell. Good luck! Dec 17, 2012 at 14:37
• Ooops yes I see; no idea how I arrived at that conclusion now - somehow I was fixated on the fermion path alone. Amazingly, you've just reminded me I do actually have an unread copy of Zee in a box somewhere, bought on the back of reading - if skipping all the math still counts as reading - Griffith's "Introduction to Elementary Particles" about a decade ago... I must dig that out too and see what it has on this stuff. Dec 17, 2012 at 18:06

There is a nice set of lecture notes, by Cvitanovic http://chaosbook.org/FieldTheory/pdf.html, from which you could extract a bit more about the combinatorics of Feynman diagrams, it also has some pretty pictures. He was also among the people that computed higher order corrections (i.e. calculated Feynman diagrams with up to three loops as you showed them above).

The thing that you are talking about reminds me of something that I've heard from condensed-matter people. You can have a look at that:

There is a peculiar method called "Diagrammatic Monte-Carlo". The authors of the method claim that the technique allows one to avoid the numerical sign problem. Though keep in mind that, as I understand, the attitude to that in the cond-mat community is rather skeptical.

• I tried to parse that DiagMC paper (the one that relates to QED) but there was too much shop-talk... did you ever look into that more closely or heard more about it? Jan 16, 2023 at 1:16