+int harmonic_oscillator_check_2b_bond(t_moldyn *moldyn,
+ t_atom *ai,t_atom *aj,u8 bc);
+
+/* harmonic oscillator potential parameter defines */
+
+// silicon
+#define HO_SC_SI 1
+#define HO_ED_SI (0.25*sqrt(3.0)*LC_SI)
+
+// carbon
+#define HO_SC_C 1
+#define HO_ED_C (0.25*sqrt(3.0)*LC_C)