Info: Particle Initial Conditions

From SIMION

Table of contents

1 Overview 2 User Programming Examples 2.1 Uniformly randomized cone angle 3 Procedure for creating a .ION file in Excel 4 Flying ions in the reverse direction or other ion transformations

Overview

Initial conditions on particles can be defined in various ways in SIMION. You can use an ION File, a FLY File, a FLY2 File, or a user program.

The ION File provides a very general method and just consists of a list of particle parameters in a straightforward "comma-separated-value" (CSV) text format as described in Section F.6 "Ion Definition File" (SIMION 8.0 manual), p. 8-6 (SIMION 7.0 manual), or here. There are a number of ways to create this file from external programs or from SIMION itself.

The FLY File allows a linear sequence of particles to be defined rather easily via the SIMION GUI, but it can only be used if your particles can be defined as a linear sequence, and it is not easy created from an external program. Therefore, FLY files are only useful in limited circumstances. FLY files are built using "Define Ions by Groups" in the SIMION ion definition screen.

The FLY2 file (http://www.simion.com/beta/fly2.html) (new in SIMION 8.0) allows particles to be defined in a very flexible but simple way. FLY2 is a replacement for the FLY format. A FLY2 file can also be converted into an ION file.

A workbench user program (initialize segment) can define particle parameters at the start of ion flight. User programs are just as flexible as ION files (e.g. it is even possible for a user program to essentially load and use an ION file), but they are even more flexible since ion parameters can be calculated in the middle of the simulation or parameterized by adjustable variables. They can require a bit more experience to create though and are not always needed.

An ION itself can be generated in many ways:

• Create with the SIMION GUI ("Define Ions Individually")
• Create from a FLY/FLY2 file
• Create from SIMION Data Recording output (click "Ion" button for .ION file format)
• Create from Excel or any other external program (e.g. C, Perl, TRIM, etc.)
• Virtual Device (http://www.simion.com/virtual_device/) includes some routines for generating .ION files for common beam types.

Creating the ION file from the SIMION GUI is done with "How are particles defined? Individually (.ION)" (or "Define Ions Individually" in version 7) on the partial definition screen. It would be tedious to use if you have many particles. Still you can use any other tool. You might be more familiar with Excel or you might generate the initial ion parameters from another simulation program (e.g. TRIM).

User Programming Examples

Uniformly randomized cone angle

The follow PRG program uniformly randomizes the directions of ions within a certain cone angle. A cone vertex angle of 180 degrees achieves a completely random direction. This code is actually an improvement on SIMION's random.prg/random.sl (http://www.simion.com/sl/examples/random.html) examples, which generate random direction that are not entirely uniform. A uniformly random direction is achieved by applying transformation derived from the fundamental transformation law of probabilities (http://www.taygeta.com/random/weibull.xml).

; PRG program to randomize ion direction within a cone.

; angle in degrees of cone vertex (use 180 for all directions)
DEFA vertex_angle 180

SEG initialize
   ; get ion's specified velocity components in polar form
   RCL ion_vz_mm
   RCL ion_vy_mm
   RCL ion_vx_mm
   >P3D
   STO speed RLUP
   STO az RLUP
   STO el RLUP
   
   ; compute uniformly randomized direction from +Y axis.
   ; c = 1 - cos(vertex_angle * PI / 180)
   1 RCL vertex_angle 3.141592 * 180 / COS - STO c
   ; el2 = arccos(1 - rand() * c) * 180 / PI + 90
   1 RAND RCL c * - ACOS 180 * 3.141592 / 90 + STO el2
   ; az2 = rand() * 360
   RAND 360 * STO az2
   
   ; convert to rectangular coordinates
   RCL el2
   RCL az2
   RCL speed
   >R3D
   
   ; rotate center back to +X axis.
   ; then rotate center back to original direction
   -90 RCL el + >ELR
   RCL az >AZR
   
   ; save new velocity vector
   STO ion_vx_mm RLUP
   STO ion_vy_mm RLUP
   STO ion_vz_mm RLUP

Here is an equivalent FLY2 file (which avoids user programming):

 standard_beam_define {
   n = 100,
   x = 0, y = 0, z = 0,
   direction = cone_distribution {
     axis_direction = {1, 0, 0},
     vertex_angle = 180
   }
 }

Here is an equivalent way to do that in Excel: simion_cone.xls

Procedure for creating a .ION file in Excel

The follow procedure is recommended for building an ION file in Excel:

• Create and save a simple .ION file using SIMION (View | Def | Define Ions Individually | Save). This will be used as a template so that you don't have to create the file entirely from scratch.
• Load the .ION file into Excel. When prompted during import, choose "delimited" format and "comma" delimiters (i.e. comma-delimited).
• Edit the data in Excel. You can even define ions using formulas.
• Save the file in Excel. Do this with "Save As" and choose "CSV (comma delimited) (*.csv)". Name it with a .ion extension (e.g. "hello.ion"), adding double quotes around the filename so that Excel doesn't try to add a .csv extension to the file name.

Note that a little care is required since Excel otherwise expects CSV file names to have a ".csv" extension rather than a SIMION ".ion" extension. If you are not careful, Excel will save/read your file in tab-delimited rather than comma-delimited format. Also, if you need to save any formulas, you must save a copy of your spreadsheet in a more intelligent format (e.g. native Excel .xls format).

For an example, see simion_cone.xls (updated on 2006-01-30 so that the angle is uniformly random).

Flying ions in the reverse direction or other ion transformations

The following procedure can be used to fly ions in the reverse direction of a previous simulation. The same general method can be used for other situations as well.

For example, consider SIMION's TOF example:

The ions start at the source at the bottom and hit the detector at the upper left. Given ion conditions at the source, SIMION calculates the ion conditions at the detector. But, what if we already have the conditions at the detector (e.g. from SIMION's data recording) and want SIMION to instead calculate the conditions at the source (or something similar)?

We start as usual by calculating ion conditions at the detector with SIMION's Data Recording feature. Enable recording to a file. Give it an extension ".ion" since we will tell SIMION to record in the ION File format:

In the Data Recording Screen, click "ION" so that SIMION records data in the ".ION" format. Note that data is recorded "when ions splat" (against the detector).

Fly the ions. SIMION writes the file "tof.ion" containing the ion conditions at the detector (and displays it's contents).

Now, we want to reverse this process. We first need to reverse the ion conditions (e.g. their velocity direction) just obtained. We can do this by editing the .ION file in Excel. When loading in Excel, be sure to specify "delimited" format and the "comma" delimiter.

Delete the first comment row (since SIMION will trip on it). Columns G and H are the azimuth and elevation angles (in degrees) respectively. To reverse the direction add (or subtract) 180 degrees to the azimuth angles and multiply the elevation angles by -1 (see figure of coordinate system on p. I-6 of the SIMION 7.0 manual).

Now save. Be sure to save it in the "CSV" (comma-separated-value) format. Include double quotes around the filename so that Excel uses the filename exactly as specified (not try to append a ".csv" to it).

Now, in SIMION, load the ION file. Ensure that the coordinates are defined relative to the workbench (in mm). These are the coordinates that SIMION Data Recordin uses, so you want the same coordinates. You should also clear the "to dev/file" box so as not to overwrite your newly modified ION file.

Refly:

Note the similarity in trajectory shape compared to before. One property of electrostatic fields is that they are reversible in this way--a property that we make use of here. How well this property holds provides measure some of simulation accuracy.

You might use similar methods to acquire data recording data on some plane and then restart the ions at the same position (maybe reversing, randomizing, or changing the ions in a sort-of like a secondary emissions effect). You might also use this type of method to transform ion data to and from another simulation program. One might alternatively use user programs to do the ion transformations.

This process can be made a bit simpler of course. A "reverse ions" button is a planned feature to be added to SIMION. [Virtual Device (http://www.simion.com/docs/virtual_device.html)] also has some functions for transforming SIMION .ION and Data Recording files.