uff_relax/params/
mod.rs

1pub mod data;
2
3use crate::atom::UffAtomType;
4use data::UFF_DATA;
5
6/// UFF Parameters for a specific atom type.
7#[derive(Debug, Clone, Copy)]
8pub struct UffParams {
9    pub r1: f64,      // Valence bond radius (Angstroms)
10    pub theta0: f64,  // Valence angle (Degrees)
11    pub x1: f64,      // vdW distance (Angstroms)
12    pub d1: f64,      // vdW energy (kcal/mol)
13    pub zeta: f64,    // vdW scale term
14    pub z_star: f64,  // Effective charge
15    pub chi: f64,     // Electronegativity (Paulings)
16    pub u_i: f64,     // Torsional barrier parameter (kcal/mol)
17}
18
19/// Returns UFF parameters for a given atom type label.
20pub fn get_uff_params(atom_type: &UffAtomType) -> Option<UffParams> {
21    let label = atom_type.as_str();
22    let label_stem = label.trim_end_matches('_');
23    
24    UFF_DATA.iter()
25        .find(|&&(id, ..)| {
26            let id_stem = id.trim_end_matches('_');
27            id == label || 
28            id.starts_with(label) || 
29            label.starts_with(id) ||
30            id.starts_with(label_stem) ||
31            label_stem.starts_with(id_stem)
32        })
33        .map(|&(_, r1, theta0, x1, d1, zeta, z_star, chi, u_i)| UffParams {
34            r1, theta0, x1, d1, zeta, z_star, chi, u_i,
35        })
36}
37
38/// Helper to get the element symbol from atomic number (Z).
39pub fn element_symbol(z: usize) -> &'static str {
40    let symbols = [
41        "H", "He", "Li", "Be", "B", "C", "N", "O", "F", "Ne",
42        "Na", "Mg", "Al", "Si", "P", "S", "Cl", "Ar", "K", "Ca",
43        "Sc", "Ti", "V", "Cr", "Mn", "Fe", "Co", "Ni", "Cu", "Zn",
44        "Ga", "Ge", "As", "Se", "Br", "Kr", "Rb", "Sr", "Y", "Zr",
45        "Nb", "Mo", "Tc", "Ru", "Rh", "Pd", "Ag", "Cd", "In", "Sn",
46        "Sb", "Te", "I", "Xe", "Cs", "Ba", "La", "Ce", "Pr", "Nd",
47        "Pm", "Sm", "Eu", "Gd", "Tb", "Dy", "Ho", "Er", "Tm", "Yb",
48        "Lu", "Hf", "Ta", "W", "Re", "Os", "Ir", "Pt", "Au", "Hg",
49        "Tl", "Pb", "Bi", "Po", "At", "Rn", "Fr", "Ra", "Ac", "Th",
50        "Pa", "U", "Np", "Pu", "Am", "Cm", "Bk", "Cf", "Es", "Fm",
51        "Md", "No", "Lr",
52    ];
53    if z > 0 && z <= symbols.len() {
54        symbols[z - 1]
55    } else {
56        "X"
57    }
58}
uff_relax/params/mod.rs

uff_relax/params/
mod.rs
