diff --git a/atomistics/workflows/quasiharmonic.py b/atomistics/workflows/quasiharmonic.py
index 856bd39a..24c4b1e5 100644
--- a/atomistics/workflows/quasiharmonic.py
+++ b/atomistics/workflows/quasiharmonic.py
@@ -67,6 +67,9 @@ def get_thermal_properties(
     volume_lst = np.array(get_volume_lst(structure_dict=structure_dict)) / np.prod(
         repeat_vector
     )
+    eng_internal_dict = {
+        key: value / np.prod(repeat_vector) for key, value in eng_internal_dict.items()
+    }
     if quantum_mechanical:
         tp_collect_dict = _get_thermal_properties_quantum_mechanical(
             phonopy_dict=phonopy_dict,
diff --git a/tests/test_quasiharmonic_lammps.py b/tests/test_quasiharmonic_lammps.py
index 2c214777..1d1380b2 100644
--- a/tests/test_quasiharmonic_lammps.py
+++ b/tests/test_quasiharmonic_lammps.py
@@ -59,16 +59,16 @@ def test_calc_phonons(self):
         self.assertTrue(tp_collect_dict["free_energy"][-1] > -2.7)
         self.assertTrue(tp_collect_dict["entropy"][0] < 0.1)
         self.assertTrue(tp_collect_dict["entropy"][0] > 0.0)
-        self.assertTrue(tp_collect_dict["entropy"][-1] < 270)
-        self.assertTrue(tp_collect_dict["entropy"][-1] > 269)
+        self.assertTrue(tp_collect_dict["entropy"][-1] < 273)
+        self.assertTrue(tp_collect_dict["entropy"][-1] > 272)
         self.assertTrue(tp_collect_dict["heat_capacity"][0] < 0.1)
         self.assertTrue(tp_collect_dict["heat_capacity"][0] > 0.0)
         self.assertTrue(tp_collect_dict["heat_capacity"][-1] < 100)
         self.assertTrue(tp_collect_dict["heat_capacity"][-1] > 99)
-        self.assertTrue(tp_collect_dict["volumes"][-1] < 66.6)
-        self.assertTrue(tp_collect_dict["volumes"][-1] > 66.5)
-        self.assertTrue(tp_collect_dict["volumes"][0] < 66.5)
-        self.assertTrue(tp_collect_dict["volumes"][0] > 66.4)
+        self.assertTrue(tp_collect_dict["volumes"][-1] < 68.6)
+        self.assertTrue(tp_collect_dict["volumes"][-1] > 68.5)
+        self.assertTrue(tp_collect_dict["volumes"][0] < 66.8)
+        self.assertTrue(tp_collect_dict["volumes"][0] > 66.7)
         thermal_properties_dict = workflow.get_thermal_properties(
             temperatures=[100, 1000],
             output_keys=["temperatures", "volumes"],
diff --git a/tests/test_quasiharmonic_lammps_functional.py b/tests/test_quasiharmonic_lammps_functional.py
index 6220fa76..db3b2e36 100644
--- a/tests/test_quasiharmonic_lammps_functional.py
+++ b/tests/test_quasiharmonic_lammps_functional.py
@@ -77,16 +77,16 @@ def test_calc_phonons(self):
         self.assertTrue(tp_collect_dict["free_energy"][-1] > -2.7)
         self.assertTrue(tp_collect_dict["entropy"][0] < 0.1)
         self.assertTrue(tp_collect_dict["entropy"][0] > 0.0)
-        self.assertTrue(tp_collect_dict["entropy"][-1] < 270)
-        self.assertTrue(tp_collect_dict["entropy"][-1] > 269)
+        self.assertTrue(tp_collect_dict["entropy"][-1] < 273)
+        self.assertTrue(tp_collect_dict["entropy"][-1] > 272)
         self.assertTrue(tp_collect_dict["heat_capacity"][0] < 0.1)
         self.assertTrue(tp_collect_dict["heat_capacity"][0] > 0.0)
         self.assertTrue(tp_collect_dict["heat_capacity"][-1] < 100)
         self.assertTrue(tp_collect_dict["heat_capacity"][-1] > 99)
-        self.assertTrue(tp_collect_dict["volumes"][-1] < 66.6)
-        self.assertTrue(tp_collect_dict["volumes"][-1] > 66.5)
-        self.assertTrue(tp_collect_dict["volumes"][0] < 66.5)
-        self.assertTrue(tp_collect_dict["volumes"][0] > 66.4)
+        self.assertTrue(tp_collect_dict["volumes"][-1] < 68.6)
+        self.assertTrue(tp_collect_dict["volumes"][-1] > 68.5)
+        self.assertTrue(tp_collect_dict["volumes"][0] < 66.8)
+        self.assertTrue(tp_collect_dict["volumes"][0] > 66.7)
         thermal_properties_dict = get_thermal_properties(
             eng_internal_dict=eng_internal_dict,
             phonopy_dict=phonopy_dict,