wimpydd is hosted by Hepforge, IPPP Durham


WimPyDD provides a list of predefined elements that can be listed typing:

 >>> WD.list_elements() 
 >>> WD.Al WD.Ar WD.C WD.Ca WD.F WD.Fe WD.Ge WD.H WD.He WD.I  
 >>> WD.Mg WD.N WD.Na WD.Ne WD.Ni WD.O WD.S WD.Si WD.W WD.Xe 

Element targets are instantiated by the element class, that takes as an argument a string with the element symbol. By convention element symbols must begin with a capital letter. The string must correspond to a file in the directory WimPyDD/Targets containing the required information. For instance the xenon element, with symbol Xe, is instantiated with:

 >>> Xenon=WD.element('Xe') 

and the class requires the text file Xe.tab in WimPyDD/Targets. The content of such plain text file is the following:

The first three rows contain the element full name, symbol, and atomic number \(Z\). For each of the element isotopes each of the following lines initializes the mass atomic mass number \(A\) in the first column, the fractional abundance in the second column and the spin \(j_T\) in the third column. Finally for pre-defined elements the last column contains and internal integer code itar is used to initialize the nuclear response functions \(W_{l,T}^{\tau\tau^{\prime}}(q)\) provided in Anand et al. and R. Catena, B Schwabe for the 30 isotopes \(T\)=\(^{12}C\), \(^{19}F\), \(^{23}Na\), \(^{28}Si\), \(^{70}Ge\), \(^{72}Ge\), \(^{73}Ge\), \(^{74}Ge\), \(^{76}Ge\), \(^{127}I\), \(^{128}Xe\), \(^{129}Xe\), \(^{130}Xe\), \(^{131}Xe\), \(^{132}Xe\) \(^{134}Xe\), \(^{136}Xe\), \(^{16}O\), \(^{40}Ar\), \(^{40}Ca\), \(H\), \(^{3}He\), \(^{4}He\), \(^{14}N\), \(^{20}Ne\), \(^{24}Mg\), \(^{27}Al\), \(^{32}S\), \(^{14}N\), \(^{20}Ne\), \(^{24}Mg\), \(^{27}Al\), \(^{32}S\), \(^{56}Fe\), \(^{59}Ni\).

Adding new nuclear targets

Besides adding a file in WimPyDD/Targets, for new elements or isotopes for which numerical calculations of the \(W_{l,T}^{\tau,\tau^\prime}(q)\) functions are not available the internal code itar must be set to zero. It is possible to provide the corresponding nuclear response function (or to override the default definitions) adding the file name_func_w.py in the directory WimPyDD/Targets/Nuclear_response_functions directory WimPyDD/Targets/Nuclear_response_functions with name one of the entries of the isotopes attribute of the element class. The file must contain a function func_w_custom(q) with the same contain a function func_w_custom(q) with the same behaviour of each of the functions contained in the func_w method. i.e. it must depend on the single argument q and return an array with shape (8,2,2)). For instance, in order to change \(W_{l,T}\) for \(T\)=\(^{129}Xe\) the custom definition must be contained in the file 129Xe_func_w.py, since 129Xe is one of the entries of the WD.Xe.isotopes attribute. A few examples of custom definition files are provided with the WimPyDD distribution. If for a nuclear isotope itar=0 and no file is present in WimPyDD/Targets/Nuclear_response_functions the corresponding nuclear response functions are set to zero and its contribution to the output of signal routines set to zero. contribution to the output of signal routines set to zero.

Main attributes and methods of the element class

 >>> dir(WD.Ge) 
     ['a', 'abundance', 'average_a', 'average_mass', 'func_w', 
     'isotopes', 'itar', 'mass', 'n_isotopes', 'name', 
     'nt_kg', 'spin', 'symbol', 'z']
 >>> WD.Ge.a
     array([70, 72, 73, 74, 76])

Returns an array with the mass numbers of each isotope.

 >>> WD.Ge.abundance
     array([0.2052, 0.274 , 0.077 , 0.3652, 0.0775])

Returns an array with the fractional abundances of each isotope.

 >>> WD.Ge.average_a

Returns the average over isotopes of the mass number.

 >>> WD.Ge.average_mass

Returns the average over isotopes of the mass in GeV.

 >>> WD.Ge.func_w[2](q=0.1)[3]
     array([[6.02002596, 3.13170007],
     3.13170007, 1.66589114]])

The attribute func_w is an array of functions of the tranferred momentum \(q\) in GeV with the nuclear response function \(W_{l,T}^{\tau\tau^{\prime}}(q) \) for each isotope (the index of the array corresponds to each entry of the isotopes attribute array). It returns an array with shape (8,2,2)=(l,tau,tau_prime) corresponding to one 2\(\times\)2 matrix for each of the 8 nuclear currents \(l\) with \(l=(0,1,2,3,4,5,6,7) \rightarrow (M,\Sigma^{\prime\prime},\Sigma^{\prime},\Phi^{\prime\prime}, \tilde{\Phi}^{\prime}, \Delta, \Phi^{\prime\prime}M, \Delta\Sigma^{\prime})\)
The example above returns the 2\(\times\)2 matrix for the nuclear response function \(W^{\tau\tau^{\prime}}_{l,T}\) = \(W^{\tau\tau^{\prime}}_{\Phi^{\prime\prime},^{73}Ge}\).

 >>> WD.Ge.isotopes
     array(['70Ge', '72Ge', '73Ge', '74Ge', '76Ge'])

Contains an array with strings with the isotopes names.

 >>> WD.itar
     array([5, 6, 7, 8, 9])

Contains an array with the internal codes itar for each isotope of the element. itar=1,...30 is a progressive number corresponding to each of the 30 isotopes for which the nuclear response functions are available in the literature. Used to load the default nuclear response functions. itar must be set to 0 for isotopes for which default nuclear response functions are not available.

 >>> WD.Ge.mass
     array([65.17 , 67.032, 67.963, 68.894, 70.756])

Returns an array with the masses of the isotopes in GeV.

 >>> WD.n_isotopes

Integer with the number of isotopes.

 >>> WD.Ge.name

String with the element name.

 >>> WD.Ge.nt_kg
     array([1.69899841e+24, 2.26864310e+24, 6.37538389e+23, 3.02375350e+24, 6.41678249e+23])

Array with the number of targets per kg for each isotope.

 >>> WD.Ge.spin
     array([0. , 0. , 4.5, 0. , 0. ])

Array with the spin \(j_T\) of each isotope.

 >>> WD.Ge.symbol

String with the element symbol.

 >>> WD.Ge.z

Integer with the element atomic number.

print() can be used to get a summary of the element properties:

 >>> print(WD.Ge)  
 >>> germanium, symbol Ge, atomic number 32, average mass 67.691, 5 isotopes.


A target is a combination of elements with stoichiometric coefficients. In WimPyDD this can be achieved by simply taking a linear combination of elements with integer coefficients. For example, to define the Octafluoropropane target \(C_3 F_8\):

 >>> carbon=WD.element('C') 
 >>> fluorine=WD.element('F') 
 >>> octafluoropropane=3*sodium+8*fluorine 

Alternatively, the user can also define a target in the following way:

 >>> c3f8=WD.target('C3F8') 

The element symbols must begin with a capital letter in order for the input string passed to the target class to be parsed correctly (so the string 'NaI' initializes sodium iodine, while 'nai' returns an error) . A short summary of the targer can be obtained by using print():

 >>> print(c3f8) 
     C3F8 contains: 
     carbon, symbol C, atomic number 6, average mass 11.182, 2 isotopes. 
     fluorine, symbol F, atomic number 9, average mass 17.689, 1 isotopes. 
     Isotope-averaged mass: 175.058723 

Main attributes of the target class

 >>> dir(c3f8) 
     ['element', 'formula', 'mass', 'n', 'n_targets', 'nt_kg']
 >>> [element.name for element in c3f8.element]
     ['carbon', 'fluorine'] 

The element attribute is an array with the element objects contained in the target.

 >>> c3f8.formula

A string with the target formula.

 >>> c3f8.mass

A float with the molecular total mass in GeV.

 >>> c3f8.n
     array([3, 8])

An array with the stoichiometric coefficients.

 >>> c3f8.n_targets

An integer with the number of elements in the target.

 >>> c3f8.nt_kg
     array([9.60469705e+24, 2.56125255e+25])

An array with the number of targets per kg of molecular composite for each element.