diff --git a/src/determinants/density_matrix_complex.irp.f b/src/determinants/density_matrix_complex.irp.f
index e6c81d33..18bae7db 100644
--- a/src/determinants/density_matrix_complex.irp.f
+++ b/src/determinants/density_matrix_complex.irp.f
@@ -106,7 +106,7 @@ END_PROVIDER
   complex*16, allocatable  :: tmp_a(:,:,:), tmp_b(:,:,:)
   integer                        :: krow, kcol, lrow, lcol
 
-  PROVIDE psi_det
+  PROVIDE psi_det psi_coef_complex
 
   one_e_dm_mo_alpha_complex = (0.d0,0.d0)
   one_e_dm_mo_beta_complex  = (0.d0,0.d0)
diff --git a/src/determinants/determinants_complex.irp.f b/src/determinants/determinants_complex.irp.f
index 692bd253..38e4126a 100644
--- a/src/determinants/determinants_complex.irp.f
+++ b/src/determinants/determinants_complex.irp.f
@@ -15,9 +15,9 @@ BEGIN_PROVIDER [ complex*16, psi_coef_complex, (psi_det_size,N_states) ]
   character*(64)                 :: label
 
   PROVIDE read_wf N_det mo_label ezfio_filename
-  psi_coef = (0.d0,0.d0)
+  psi_coef_complex = (0.d0,0.d0)
   do i=1,min(N_states,psi_det_size)
-    psi_coef(i,i) = (1.d0,0.d0)
+    psi_coef_complex(i,i) = (1.d0,0.d0)
   enddo
 
   if (mpi_master) then