diff --git a/python/sumk_dft.py b/python/sumk_dft.py
index 7a7db48a..b315f08c 100644
--- a/python/sumk_dft.py
+++ b/python/sumk_dft.py
@@ -1478,6 +1478,7 @@ class SumkDFT(object):
                                 for imu in range(self.n_orbitals[ik, isp]):
                                     fout.write("%.14f  %.14f " % (deltaN[sp][ik][
                                                inu, imu].real, deltaN[sp][ik][inu, imu].imag))
+                                fout.write("\n")
                             fout.write("\n")
                         fout.close()
         elif dm_type == 'vasp':