diff --git a/source/source_relax/relax_driver.cpp b/source/source_relax/relax_driver.cpp
index ef8c9b62e6f..ef97baf4dd7 100644
--- a/source/source_relax/relax_driver.cpp
+++ b/source/source_relax/relax_driver.cpp
@@ -172,6 +172,14 @@ void Relax_Driver::relax_driver(
         ++istep;
     } // end while (istep <= inp.relax_nmax && !stop)
 
+    // output the final relaxed structure in CIF format
+    if (inp.calculation == "relax" || inp.calculation == "cell-relax")
+    {
+        ModuleIO::CifParser::write(PARAM.globalv.global_out_dir + "STRU_FINAL.cif",
+                                   ucell,
+                                   "# Generated by ABACUS ModuleIO::CifParser",
+                                   "data_?");
+    }
 
 	if (inp.calculation == "relax" || inp.calculation == "cell-relax")
 	{