anchor.diff:   double precision :: delta_e_ja(N_states)
    double precision :: hja
    double precision :: contrib_hij 
  + double precision :: fock_operator_local(n_act_orb,n_act_orb,2)
    accu_contrib = 0.d0
   !matrix_1h2p = 0.d0
    
  
positive.code:                  ! | det_tmp_bis > = a^{\dagger}_{borb,kspin} a_{aorb,kspin} | det_tmp > 
                   call do_mono_excitation(det_tmp_bis,list_act(aorb),list_act(borb),kspin,i_ok)
                   if(i_ok .ne. 1)cycle
  
                   ! < det_tmp     | H | det_tmp_bis >  = F_{aorb,borb}
                   call i_H_j(det_tmp_bis,det_tmp,N_int,hab)
                   ! < jdet | H | det_tmp_bis > = phase * (ir|cv)
  !                call i_H_j(det_tmp_bis,psi_det(1,1,idx(jdet)),N_int,hja)
                   call get_double_excitation(det_tmp_bis,psi_det(1,1,idx(jdet)),exc,phase,N_int)
  
negative.code: !!!!!!!!!!!!!!!!!!!!!!!!!!! determination of the connections between I and the other J determinants mono excited in the CAS
  !!!!!!!!!!!!!!!!!!!!!!!!!!!! the determinants I and J must be connected by the following operator 
  !!!!!!!!!!!!!!!!!!!!!!!!!!!! <Jdet | a^{\dagger}_b a_{a}  | Idet> 
               integer :: i_hole,i_part
         do jdet = 1, idx(0)
           if(idx(jdet).ne.idet)then
  !         print*, degree(jdet)
  