diff --git a/tests/HBO.xyz b/tests/HBO.xyz
new file mode 100644
index 00000000..796a1e19
--- /dev/null
+++ b/tests/HBO.xyz
@@ -0,0 +1,5 @@
+3
+HBO Geo: Experiment Mult: 1 symmetry: 14
+B    0.0    0.0    1.166
+H    0.0    0.0    0.0
+O    0.0    0.0    2.366
diff --git a/tests/unit_test/unit_test.py b/tests/unit_test/unit_test.py
index 19d129d1..7ce1c0f6 100755
--- a/tests/unit_test/unit_test.py
+++ b/tests/unit_test/unit_test.py
@@ -169,6 +169,7 @@ def run_hf(geo, basis, mult=1, pseudo=False, remove_after_sucess=True):
 
     ref_energy["sto-3g"]["methane"] = Energy(-39.7267433402, None)
     ref_energy["vdz"]["SO2"] = Energy(None, -41.48912297776174)
+    ref_energy["vdz"]["HBO"] = Energy(None, -19.11982537379352)
 
     # ~#~#~#~#~#~#~#~#~#~#~#~#~#~#~ #
     # G l o b a l _ v a r i a b l e #
@@ -378,7 +379,7 @@ class ValueTest(unittest.TestCase):
                                          pseudo=False))
 
     def test_hf(self):
-        self.assertTrue(run_hf(geo="SO2",
+        self.assertTrue(run_hf(geo="HBO",
                                basis="vdz",
                                mult=1,
                                pseudo=True))