CFEL-CMI/cmifly
Trajectory-based simulations of the deflection of molecular beams in inhomogeneous electric fields
CMI fly
A simple Python program for trajectory-based simulations of the deflection of molecular beams in
inhomogeneous electric fields.
The initial version of this code was published in the following paper and any use of the code should
reference that publication:
Yuan-Pin Chang, Daniel Horke, Sebastian Trippel, Jochen Küpper: Spatially-controlled complex
molecules and their applications, International Reviews in Physical Chemistry 34, 557–590
(2015); arXiv:1505.05632
[physics].
See the accompanying license and the documentation for further
details.
Getting Started with CMIfly
For now, get the latest version from github and run
CMIfly.py.
This requires a Stark-energy file from CMIstark for the
molecule/state of interest.
Documentation
For now, this code is undocumented; this needs to be changed;-)
Some notes on the usage for different setups
For the a-type deflector, use the following lines in CMIfly.py:
deflector_fieldnorm_filename = 'deflector_field_norm'
deflector_fieldgradient_filename = 'deflector_field_gradient'
# deflector_voltage for which the fields were calculated (kV)
deflector_field_voltage = 5.
# geometric boundary of the deflector
rod_center = [0, 0.0036] # rod center (x,y)
rod_radius = 0.003 # rod radius
trough_center = [0, 0.002408] # trough center (x,y)
trough_radius = 0.0032 # trough radiusFor the b-type deflector, use the following lines in CMIfly.py:
deflector_fieldnorm_filename = 'Ec.norm.txt'
deflector_fieldgradient_filename = 'gt.grad.txt'
# deflector_voltage for which the fields were calculated (kV)
deflector_field_voltage = 60.
# geometric boundary of the deflector
rod_center = [0, 0.0036] # rod center (x,y)
rod_radius = 0 #0.003 # rod radius
trough_center = [0, 0.002408] # trough center (x,y)
trough_radius = 0 #0.0032 # trough radius
Since the rod_radius = 0, it will now use the fields to determine the edge of the deflector