> TVS5@ w%bjbj22 (6XXhVVVVVVVj
8
,6j8VV:^".8080808080808$:R<fT8V{{{T8VVi8!!!{VV.8!{.8!!!5VVz7J .w
u r688086W=o!:W=0z7jjVVVVW=Vz7>D,!p$T8T8jjd!jjSpecial Applications of TRIM
Calculating Ion Energy Reducers (Energy Degraders)
TRIM may be used to calculate the energy loss through foils or thick blocks of materials used to lower the energy of a beam of particles. This is often done for light ions, such as protons, to obtain lower energies quickly. For example, if one starts with a beam of monoenergetic protons at 158.6 MeV and introduces a block of plexiglas (called perspex in Europe) 14.1 cm thick, the exiting beam of particles will be at 30 MeV. Degrader blocks made with light atoms can reduce particle energy with a minimum of straggle and little long term radioactivity (plexiglas is C6H4O2).
Two problems occur in planning "energy degraders". First is the calculation of the proper block thickness, and second is the calculation of the final beam energy spread. Below we illustrate both calculations.
The calculation of the required degrader thickness consists of four basic steps. We illustrate these steps below using the example of protons being slowed by a block of plexiglas from an initial energy of 158.6 MeV to a final energy of 30 MeV.
1 Make a table of proton stopping in plexiglas by starting SRIM and going to Stopping and Range Tables.
The default ion is H so this can remain.
Change the Highest Ion Energy to 160000 keV.
Press Compound Dictionary and open the file Common Target Materials. In this list is Plexiglas / Acrylic. Select this material and press Add to Target. A pop-up window will ask if you wish to keep the Compound Correction to the stopping of H ions in Plexiglas. Answer Yes.
At the bottom left of the setup window is Stopping Power Units. Open the drop-down menu as specify MeV/mm. This will tell you the energy loss of the beam in convenient units.
Press Calculate Table.
2 Using the table generated in (1), find the range of protons for the initial beam energy and estimate the depth at which the protons will be at 30 MeV by looking at the ranges and working backwards. For example, the table will show that 160 MeV protons have a range of 155 mm in plexiglas, and 30 MeV protons will go 7.7 mm. Hence your first guess for the thickness of plexiglas to reduce 160 MeV protons to 30 MeV will be 155 - 7.7 = 147.3 mm thick. Further, since the initial proton beam is not 160 MeV but 158.6 MeV, one interpolates the distance protons travel between 160 MeV to 158.6 MeV as 2.3 mm, so the initial block size for an incident 158.6 MeV beam will be 147.3 - 2.3 = 145 mm.
3 Close the Table window.
Press the Main Menu button.
Press TRIM to go to the TRIM setup window.
Specify ions of H (158.6 MeV) into a target of plexiglas (145 mm).
At the bottom middle of the Setup window is Output Disk Files. Check the box beside Transmitted Ions.
Double check your inputs, and then press Save Input and Run TRIM.
Calculate about 200 ions, then stop. Close TRIM (do not Save). Quit the input window.
3 You now need to look at the file: SRIM Outputs/TRANSMIT.txt. This will contain the transmitted protons through 145 mm of Plexiglas. You need to find the mean energy of these particles, and their average deviation from this mean (straggling). This can be done in Excel or other statistical packages. There is a BASIC program included with SRIM (/SRIM Tools/TRANSMIT.BAS) if you can run this type of program. It will read the file TRANSMIT.txt and give you a full statistical analysis.
4 Use the difference between your desired final energy (30 MeV) and mean final energy found in step 3 above, along with the stopping power tables, made in step 1 above, to estimate a new target thickness. For example, if the protons exit at 32 MeV and the stopping power for protons at 30 MeV is 2.16 MeV/mm, then you need (2 MeV)/(2.16 MeV/mm) = 0.92 mm more plexiglas in the block to bring the beam down to 30 MeV. You then use TRIM to make up a new plexiglas target with a width of 145 mm (initial estimate) + 0.92 mm (correction) = 145.92 mm. Loop through steps 2-4 until you come get as close as you want to the desired thickness.
The final energy straggle of the beam after it has passed through the energy reducing block is a mixture of two components: (1) the straggle introduced by random collisions in the degrader block and (2) the broadening of any incident beam energy straggle by the block (straggle is the spread of the ion beam energy spectrum). We calculate the straggle of the incident monoenergetic beam with the TRIM calculations above, steps 1-4, and we find, for example, a straggle of about 2.0 MeV when we degrade protons from 160 MeV to 30 MeV using plexiglas (this comes out of the analysis program TRANSMIT.BAS) To this broadening, we have to add the straggle due to expansion of the initial beam energy spread, which may be the dominant factor in the final energy width. Assuming that the incident beam is not monoenergetic, but has a spread of 1 MeV at 160 MeV, we will find that this spread expands to 3.6 MeV by the time the beam degrades to 30 MeV. This expansion occurs because of the non-linearity of stopping powers with particle energy as explained below.
Assume that the energy straggle of the incident beam is (. We show below that the final energy straggle of the exiting beam, (', can be estimated as :
(' = (S'/S) ( ,
where S and S' are the stopping powers, dE/dx, in the degrader material at the initial and final beam energy.
Consider two particles of the incident beam with energies E1 and E2, separated by the incident beam energy straggle, (. The final energy straggle of these two particles, after transiting the energy-lowering block of material, is defined as (' = E1' E2'. where E1' and E2' are the final energies of particles 1 and 2.
Considering the incident particles, assuming E1 > E2, , then particle #1 will travel a short distance into the degrader block before its energy will be reduced to that of particle #2. This distance, x, is about x = ( / S, where S is the particle stopping power (dE/dx) at energy E. After particle #1 reaches depth x, then both particles #1 and #2 will transit through the block similarly until particle #1 reaches the end of the block. Particle #2 will still have more of the block to traverse, an amount x' = (' / S' , where (' = E1' - E2', the final energy difference between the two particles, and S' is the stopping power at the final beam energy, E'. But the block length is the same for both particles, so x = x', and therefore ( / S = (' / S' . This can be solved for the energy straggle, (' = (S'/S) ( .
Note that several approximations are made in this argument, especially that S(E1) = S(E2) and S(E1') = S(E2'), and that the total energy loss is much larger than the straggling, (E - E') >> ( .
For an example, protons at 160 MeV in plexiglas have a stopping power of 0.593 MeV/mm and at 30 MeV have a stopping power of 2.157 MeV/mm. Assuming the initial beam straggle, (, is 1 MeV, then the final beam energy straggle, (' = (1 MeV) (2.157 MeV/mm) / (.593 MeV/mm) = 3.6 MeV.
The total beam spread through the energy-degrader block is the rms total value of the two separate quantities : (1) the beam straggle through the block, 2.0 MeV, and (2) the beam energy broadening due to non-linear stopping powers, 3.6 MeV, which, when added in quadrature, gives a total straggle of 4.1 MeV .when protons at 160 MeV are degraded to 30 MeV.
The above comments come from Dr. Bernard Gottschalk, Harvard Cyclotron Laboratory, Harvard University, USA.
SRIM Special Degrader Blocks.doc Page PAGE 1 of NUMPAGES 2 DATE \@ "M/d/yyyy" 3/13/2003
P
XY
2
E
\
'*,Vj24য় h3Rh3R5ehrh3Rh3R56h3Rh3R5h3RhqBhqB5 hqBhqB5ehrhqB56 hqB5hqBhqB56
hqBEHhqB
hqBCJh 5CJ(aJ(h h 5CJ(aJ(4P
dX
,PO
&F gdOBHh^h`gdOBH
&Fgd3RhP^h`gdqBPgdqBPgdqB $Pa$gdqB%v%*,:EMO8<OPQs456QR
OPQimnͼͼͷͭͣ۷۟۟۟۟۟۔۔ۅ{ jshqB5CJ
hqB5CJhOBHhqB5 jshOBHhqB5hOBHhOBHhOBH56h6h656 hqB5 h6h65ehrh6h6h3R56hqB h3Rh3R5ehrh3Rh3R h3R5h3R0O5i/^! ":#$
%%t%u%v%w%$
!($a$gdN
&dPgdqBP&dPgdqBPgdqBPgdOBHhP^h`gdOBHn{|}*+12de]^bcd- . / = > ? C D I J !!!!!L!M!Z![!^!!!!!!!!!併䪢䪢䪢䚕 hqB5 jshqB5hoXhqB5 jshoXhqB5hqBCJEHhoXhoXH*hOBHhOBHhqB5 jshOBHhqB5hOBHhOBHH*hqB
hqBCJ jshqB5CJ
hqB5CJ:!"""""###%%,%1%5%6%7%=%>%?%@%A%D%E%O%P%Q%R%S%T%h%i%r%s%t%u%ּ|vlvlal]hOBHhOBH0JmHnHujhOBH0JU
hOBH0JhoX0JmHnHujhNhOBH0JUhNhOBH0J hoX0J5CJ aJ mHnHuhNhOBH0J5CJ aJ $jhNhOBH0J5CJ UaJ hNhOBHhOBHCJaJhNhOBHCJaJhoXhqB5 jshoXhqB5hqB"u%v%w%hqBhoX&1h:pe$/ =!"#$%f@fNormal$<5$7$8$9DH$a$$B*CJ\_HaJmH phsH tH \@\ Heading 1$@&a$5B*CJ OJQJ\aJphZ@Z Heading 2
xx@&5B*CJOJQJ\aJph@i Heading 36$
p@!dh<5$7$8$9D@&H$56CJOJQJ\^JaJF@Fi Heading 4$
&F<@&56DA@DDefault Paragraph FontRiRTable Normal4
l4a(k@(No ListL@L
iTOC 3!
X !]^56B@BiTOC 2
X !]5R"@R Caption
xx5B*CJOJQJ\aJph4@"4NHeader
!4 @24NFooter
!.)@A.NPage Numberw6PdX,
P
O5
i/^ :
tux000p000 0 0 0 0 000 0 0 0 0p 0 0000p0p0p0p000000p0@0@000PdX5
i/^ :
x@06@0A{@0A{@0A{@0A{@0A{@0A{@0A{@0A{@0A{@0A{@0A{@0A{@0A{@0A{0wgggjn!u%w%Ow%v%(/16ACEZdj!x
_Toc318776473
_Toc318778560
_Toc318778918
_Toc318779514
_Toc318783041
_Toc318785022
_Toc318792983
_Toc318795914
_Toc334002300
_Toc334425237
_Toc336068119
_Toc336077306
_Toc336144424
_Toc336152753
_Toc336152926
_Toc336153123
_Toc336660235
_Toc336660524
_Toc336661800
_Toc336754075
_Toc336855488
_Toc337376343
_Toc337435022
_Toc340898170
_Toc344277587
_Toc345230552
_Toc345319577
_Toc347203137
_Toc318172250
_Toc318282438
_Toc318601812
_Toc318603896
_Toc318732873
_Toc318775188
_Toc290953917
_Toc307221020
_Toc307222540
_Toc307222916
_Toc307293625
_Toc307297653
_Toc307297829
_Toc307298023
_Toc307388601
_Toc307391658
_Toc307391840
_Toc307394103
_Toc307395680
_Toc307396068
_Toc307630643
_Toc311433945
_Toc312235320
_Toc314305082
_Toc317496405
_Toc317499082
_Toc317583992
_Toc317584337
_Toc317584487
_Toc317584777
_Toc317663012
_Toc317667889
_Toc317678558
_Toc317762835
_Toc317938427
_Toc317938609
_Toc317999172
_Toc317999386
_Toc317999551
_Toc317999820
_Toc318000315
_Toc318034479
_Toc318034679
_Toc318083272
_Toc318172259
_Toc318282447
_Toc318601821
_Toc318603905
_Toc318732882
_Toc318775197
_Toc318776483
_Toc318778570
_Toc318778928
_Toc318779524
_Toc318783051
_Toc318785032
_Toc318792993
_Toc318795924
_Toc334002312
_Toc334425249
_Toc336068132
_Toc336077320
_Toc336144438
_Toc336152767
_Toc336152940
_Toc336153137
_Toc336660244
_Toc336660532
_Toc336661808
_Toc336754083
_Toc336855496
_Toc337376351
_Toc337435030
_Toc340898178
_Toc344277595
_Toc345230560
_Toc345319585
_Toc347203145 DDE_LINK2 DDE_LINK1
DDE_LINK24
DDE_LINK23
DDE_LINK20
DDE_LINK19
DDE_LINK22
DDE_LINK21
DDE_LINK28
DDE_LINK27
DDE_LINK14
DDE_LINK13
DDE_LINK26
DDE_LINK25
mmLLx
rstuvw !"#$%&'()*+,-./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO||[[
xWܔW쌉WWܔxx8*urn:schemas-microsoft-com:office:smarttagsCityB*urn:schemas-microsoft-com:office:smarttagscountry-region9*urn:schemas-microsoft-com:office:smarttagsplaceZg|
}
dm;
D
!.5:FOy|AKux|
}
dm;
D
PRnrbd26MQux3333333333333333333 8Ouxux. ?@ABDEFGHIJLMNOPQRURoot Entry FWData
1Table$=WordDocument(6SummaryInformation(CDocumentSummaryInformation8KCompObjj
FMicrosoft Word Document
MSWordDocWord.Document.89q