1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92

use crate::particle::Particle;
use crate::decay::decayimpl::DecayMode;

/// Shared trait for atoms 
pub trait Atom: Clone {
    /// Atomic number
    fn atomic_num(&self) -> u64;
    
    /// Atomic mass in Daltons
    fn am(&self) -> f64;
    /// Atomic mass in kilograms
    fn am_kg(&self) -> f64;

    //fn am_ev(&self) -> f64;
    /// Mass defect or the difference between the empirical mass and the mass of the constituents, in Daltons
    fn mass_deficit(&self) -> f64;
    // Mass defect in Kilograms
    fn mass_deficit_kg(&self) -> f64;
    /// Mass defect in Joules
    fn mass_deficit_j(&self) -> f64;
    /// Mass defect in MeV  fix these, same as binding energy
    fn mass_deficit_ev(&self) -> f64;
     
    fn binding_energy(&self) -> f64;

    fn binding_energy_j(&self) -> f64;
    /// Spin as a i8 pair
    fn spin_parity(&self) -> (i8, i8);
    /// Electron affinity in kj/mol
    fn electron_affinity(&self) -> f64;
    /// Electron affinity in MeV
    fn electron_affinity_ev(&self) -> f64;
    ///Returns the ionization energies for all known levels. Values are in kj/mol
    fn ionization_energies(&self, level: usize) -> Option<f64>;
    ///Returns the ionization energies for all known levels. Values are in MeV
    fn ionization_energies_ev(&self, level: usize) -> Option<f64>;
    /// Returns Oganov-Tantardini values, the current best evaluation
    fn electronegativity(&self) -> f64;
    /// Mullikan electronegativity
    fn mullikan_en(&self) -> f64;
    /// Allen electronegativity
    fn allen_en(&self) -> f64;
    /// Pauling electronegativity
    fn pauling_en(&self) -> f64;
    /// Covalent radii of the first three bonds
    fn covalent_radii(&self, bond: usize) -> Option<f64>;
    /// Ionic radii
    fn ionic_radii(&self) -> f64;
    /// Van der Waal radius in crystalline structure
    fn vdr_crystal(&self) -> f64;
    /// Van der Waal radius in isolated atoms
    fn vdr_isolated(&self) -> f64;
    
    /// Half-life of nuclide/isomer. TotalDecay mode returns the half-life, all other modes return  the partial half-life 
    /// # Nan
    /// Particle is unstable but the selected decay mode is not supported
    /// # Inf
    /// Particle is stable
    fn half_life<T: DecayMode>(&self) -> f64;
    /// The mean lifetime of nuclide/isomer
    fn mean_lifetime<T: DecayMode>(&self) -> f64;
    /// Returns the probable decay modes as a string
    fn decay_mode(&self) -> String;
    
    fn daughter_energetic<T: DecayMode>(&mut self) -> (f64, Vec<Particle>);
    
    /// Returns the daughter nuclide 
    /// # None
    /// If nuclide has not been observed to decay by the mode, returns None
    fn daughter<T: DecayMode>(&self) -> Option<Self>;
    
    /// Decay constant in seconds of nuclide/isomer. TotalDecay mode returns the decay constant, all other modes return  the partial decay constant 
    /// # Nan
    /// Particle is unstable but the selected decay mode is not supported
    /// # Inf
    /// Particle is stable
    fn decay_constant<T: DecayMode>(&self) -> f64;
    /// Returns the probability of the nuclide to decay after the time in seconds provided
    fn decay_probability<T: DecayMode>(&self, time: f64) -> f64;
    
    /// Checks if nuclide/isomer would decay in the selected time
    fn decay_time<T: DecayMode>(&self, time: f64) -> bool;

    /// Continously performs decay throughout the time selected, collecting all particles into a vector with decay energies.
    fn decay<T: DecayMode>(&mut self, time: f64) -> (f64, Vec<Particle>);
    
    /// Q-value (total energy) of a nuclear decay, regardless of whether it is observed
    /// # NAN
    /// Returns NAN if this decay mode results in a nonexistent nuclide
    fn decay_q<T: DecayMode>(&self) -> f64;
        
}