Function calculate_bond_gradient 

fn calculate_bond_gradient(
    geometry: &Geometry,
    a1: usize,
    a2: usize,
) -> [f64; 3]