WOLFRAM NOTEBOOK

  • Note: I wrote this Wolfram notebook after reading this interesting blog post by Jeremy E. Monat: https://bertiewooster.github.io/2023/12/26/Molecular-isotopes-1-permutations.html
  • In[]:=
    ChlorineAbundances=
    chlorine
    ELEMENT
    ["IsotopeAbundances"]
    Out[]=
    chlorine-35
    75.78
    %
    ,
    chlorine-37
    24.22
    %
    In[]:=
    SulfurAbundances=
    sulfur
    ELEMENT
    ["IsotopeAbundances"]
    Out[]=
    sulfur-32
    94.93
    %
    ,
    sulfur-33
    0.76
    %
    ,
    sulfur-34
    4.29
    %
    ,
    sulfur-36
    0.02
    %
    In[]:=
    MoleculePlot[Interpreter["Chemical"]["Cl2S"]]
    Out[]=
    In[]:=
    ds=Dataset[Flatten[Table[<|"Chlorine1"->c1,"Chlorine2"->c2,"Sulfur1"->s1,"MolarMass"->c1["MolarMass"]+c2["MolarMass"]+s1["MolarMass"],"RelativeAbundance"->(QuantityMagnitude[ChlorineAbundances[c1]]QuantityMagnitude[ChlorineAbundances[c2]]QuantityMagnitude[SulfurAbundances[s1]])/10^6|>,{c1,Keys[ChlorineAbundances]},{c2,Keys[ChlorineAbundances]},{s1,Keys[SulfurAbundances]}],2]]
    Out[]=
    Chlorine1
    Chlorine2
    Sulfur1
    MolarMass
    RelativeAbundance
    chlorine-35
    chlorine-35
    sulfur-32
    101.909777g/mol
    0.5451
    chlorine-35
    chlorine-35
    sulfur-33
    102.909164g/mol
    0.0044
    chlorine-35
    chlorine-35
    sulfur-34
    103.905572g/mol
    0.0246
    chlorine-35
    chlorine-35
    sulfur-36
    105.904786g/mol
    0.0001
    chlorine-35
    chlorine-37
    sulfur-32
    103.906826g/mol
    0.1742
    chlorine-35
    chlorine-37
    sulfur-33
    104.906214g/mol
    0.0014
    chlorine-35
    chlorine-37
    sulfur-34
    105.902622g/mol
    0.00787
    chlorine-35
    chlorine-37
    sulfur-36
    107.901836g/mol
    0.00004
    chlorine-37
    chlorine-35
    sulfur-32
    103.906826g/mol
    0.1742
    chlorine-37
    chlorine-35
    sulfur-33
    104.906214g/mol
    0.0014
    chlorine-37
    chlorine-35
    sulfur-34
    105.902622g/mol
    0.00787
    chlorine-37
    chlorine-35
    sulfur-36
    107.901836g/mol
    0.00004
    chlorine-37
    chlorine-37
    sulfur-32
    105.903876g/mol
    0.05569
    chlorine-37
    chlorine-37
    sulfur-33
    106.903264g/mol
    0.00045
    chlorine-37
    chlorine-37
    sulfur-34
    107.899672g/mol
    0.00252
    chlorine-37
    chlorine-37
    sulfur-36
    109.898886g/mol
    0.00001
    In[]:=
    ListPlot[ds[All,{QuantityMagnitude[#["MolarMass"]],#["RelativeAbundance"]}&],Frame->True,FrameLabel->{"Molar mass","Relative abundance"},PlotLabel->"Sulfur dichloride isotope abundance",Filling->Axis,ImageSize->Large]
    Out[]=
    Wolfram Cloud

    You are using a browser not supported by the Wolfram Cloud

    Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.


    I understand and wish to continue anyway »

    You are using a browser not supported by the Wolfram Cloud. Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.