From 9916d70d84bc3e36a669cddce31fc9d6b759068c Mon Sep 17 00:00:00 2001 From: "Ryan M. L. McFadden" Date: Tue, 3 Jan 2023 12:35:26 -0400 Subject: [PATCH] use updated stopping coefficients for hydrogen-like projectiles - The stopping coefficients from ICRU Report 49 (published in 1993) are now fairly old, with many originating from the even older compilation by Anderson & Ziegler (published in 1977). For some target elements, little experimental data was available at the time of publication and the tabulated values are not necessarily reliable. - This patch provides revised values for stopping coefficients in select elemental targets derived from Varelas-Biersack fits to stopping data compiled in the IAEA's database (https://www-nds.iaea.org/stopping/). For further details on their determination, see e.g.: https://doi.org/10.48550/arXiv.2212.11081 - When creating input files for trimspNL using the GUI, the default behaviour is changed to prefer using the revised coefficients. In instances where updated values are unavailable (e.g., due to lack of experimental data), fall back to using the coefficients from ICRU Report 49. --- TrimSPlib.js | 70 +++++++++++++++++++++++++++++++++++++++++++++++++++- 1 file changed, 69 insertions(+), 1 deletion(-) diff --git a/TrimSPlib.js b/TrimSPlib.js index 4b1679d..d98d44a 100644 --- a/TrimSPlib.js +++ b/TrimSPlib.js @@ -448,7 +448,75 @@ function StoppingCoefficients(Element) "U":"7.29000,8.20400,19180.00000,586.29999,0.00267", "":"0.00000,0.00000,0.00000,0.00000,0.00000" }; - return StoppingCoefficientsICRU[Element]; + + // Updated Anderson-Ziegler table values for hydrogen projectiles using + // up-to-date stopping data from the IAEA database (https://www-nds.iaea.org/stopping/). + // For further details on their determination, see e.g.: + // R. M. L. McFadden at al., arXiv:2212.11081 [cond-mat.supr-con] + // https://doi.org/10.48550/arXiv.2212.11081 + var StoppingCoefficientsIAEA = { + "H": "1.21941306,1.39705253,316.72610292,3622.08960714,0.03609996", + "Li": "2.55867734,2.98398784,672.22623329,223.63937274,0.05756479", + "Be": "2.13846776,2.43024550,1353.80005691,500.78353377,0.01143076", + "B": "2.92598306,3.33271910,1523.83401267,574.31712932,0.01671298", + "C": "2.51057212,2.84874359,1325.63086028,2111.57332043,0.03815981", + "N": "3.03231463,3.44621744,1287.27189959,3445.25706665,0.06909550", + "O": "2.58218404,2.91913321,2066.69556610,1937.28839450,0.01708201", + "Mg": "4.32219620,4.95332934,2711.05757349,99.32422512,0.01692829", + "Al": "4.16266645,4.73942640,3007.21967155,192.15430476,0.01478852", + "Si": "4.56277193,5.17964542,3812.93366288,249.45105835,0.00951865", + "Cl": "5.18954121,5.87010049,5726.20198003,338.90592713,0.00527894", + "Ca": "7.38419020,8.41800474,4257.94547049,311.84732092,0.01482602", + "Sc": "6.27366940,7.11486087,5202.32117316,349.64931132,0.00808596", + "Ti": "4.89215322,5.52996090,4733.05002706,821.54929980,0.01141794", + "V": "4.39418898,4.95984106,5216.28902299,873.00768319,0.00982301", + "Cr": "4.22355616,4.76724205,4805.40719904,1062.25237873,0.01508696", + "Mn": "3.45870888,3.89580224,5452.18588822,1778.16054656,0.01057380", + "Fe": "3.69846643,4.16916938,5392.90710239,1040.38447053,0.01152298", + "Co": "3.55890909,4.00775461,6796.04145497,933.86447935,0.00558368", + "Ni": "3.42262695,3.86234880,4759.65515188,520.18934465,0.02026669", + "Cu": "3.54912421,4.02142201,4344.22605325,128.32905834,0.02886218", + "Zn": "4.31895518,4.91763977,5346.03925560,46.16548699,0.01518164", + "Ga": "4.11015447,4.63998531,7681.41419449,149.38563152,0.00628772", + "Ge": "5.80984815,6.60667983,6039.04808738,87.19744473,0.01158192", + "Se": "5.58734301,6.31119973,8423.01872583,221.61559131,0.00530434", + "Br": "6.28480648,7.09489448,11732.76411988,150.96205589,0.00281298", + "Kr": "6.47378266,7.31581518,9199.06183849,216.98680923,0.00437111", + "Rb": "8.48759187,9.62951209,7271.11124671,272.67996745,0.00818815", + "Sr": "8.66447838,9.84912713,5979.31960977,334.15705556,0.01184544", + "Y": "7.54854264,8.54789357,6885.34110669,429.53252110,0.00930644", + "Zr": "6.24168401,7.04236496,10674.06252597,269.66181445,0.00339805", + "Nb": "5.96194954,6.72632300,10250.57513827,275.25886405,0.00386599", + "Mo": "5.58867862,6.29983727,9573.16605660,523.57577519,0.00473082", + "Pd": "4.74762031,5.34871025,8086.17802763,942.05721488,0.00842671", + "Ag": "5.42052240,6.17824432,4642.18214950,90.75350922,0.06061521", + "Cd": "4.89128643,5.50754704,11831.77188508,396.06905848,0.00412049", + "In": "5.98978289,6.77050040,6343.54720062,564.65410792,0.02270520", + "Sn": "6.73191335,7.61854414,6946.18830598,347.19391310,0.01677246", + "Sb": "7.92400389,8.99628409,5736.29881090,403.11722206,0.02973016", + "Xe": "8.96826400,10.16224223,8543.27304860,283.90067433,0.01019984", + "La": "8.16949326,9.21144395,13933.96869212,429.99442382,0.00347440", + "Nd": "8.60801371,10.26060598,6798.65217629,7.06731816,0.02713704", + "Gd": "8.38379061,9.50012183,10804.38767518,110.71066098,0.00761735", + "Tb": "8.69907359,10.36476708,5985.09307671,7.76175759,0.04237031", + "Dy": "6.69350026,7.81543484,5289.65107666,12.11560717,0.05746489", + "Er": "4.31497049,4.85346547,12995.54504003,691.00380321,0.00499022", + "Yb": "6.38304258,7.21491284,12327.74719685,82.30064162,0.00516845", + "Lu": "7.00205766,8.09268611,4589.67801229,31.06321362,0.11142124", + "Hf": "5.32763157,6.00140558,13853.49562768,206.93150053,0.00411522", + "Ta": "5.15280708,5.80545310,12417.25115167,231.48549646,0.00530723", + "W": "4.81261962,5.41579909,15245.84079976,286.09077182,0.00344862", + "Re": "4.26021967,4.80273372,3191.05524918,76051.48730893,0.00003161", + "Ir": "4.16543126,4.68849191,6887.94139926,4263.90641919,0.00000004", + "Pt": "4.62228847,5.20033612,16308.49878855,266.02263275,0.00309560", + "Au": "4.84253074,5.44769072,15602.21662947,412.25485588,0.00331383", + "Pb": "9.42473604,17.46235676,11946.86905435,0.00000001,0.00652912", + "Bi": "10.25823796,12.26670424,7765.05387481,7.37274504,0.02310050" + }; + + // prefer using stopping coefficients derived from data in the IAEA database; + // however, if they are unavailable, return those from ICRU Report 49 + return (Element in StoppingCoefficientsIAEA) ? StoppingCoefficientsIAEA[Element] : StoppingCoefficientsICRU[Element]; }