Plot Isotope Decay Networks

Construct isotope decay networks with the new RelationGraph function.

In[1]:=

DaughterNuclides[s_List] := DeleteCases[ Union[Apply[Join, Map[IsotopeData[#, "DaughterNuclides"] &, DeleteCases[s, _Missing]]]], _Missing]; ReachableNuclides[s_List] := FixedPoint[Union[Join[#, DaughterNuclides[#]]] &, s]; verts = ReachableNuclides[{"Uranium235"}]; vsym = IsotopeData[#, "Symbol"] & /@ verts; opts = {VertexLabels -> Thread[verts -> (Placed[#, Center] & /@ vsym)], PlotTheme -> "ClassicDiagram", VertexSize -> {.5, .3}}; DaughterNuclidesQ[s1_, s2_] := (s1 =!= s2 && MemberQ[DaughterNuclides[{s1}], s2]);

In[2]:=

RelationGraph[DaughterNuclidesQ, verts, opts]

Out[2]=