diff --git a/plugins/local/tc_bi_ortho/tc_utils.irp.f b/plugins/local/tc_bi_ortho/tc_utils.irp.f
index 4dfd4316..067e8719 100644
--- a/plugins/local/tc_bi_ortho/tc_utils.irp.f
+++ b/plugins/local/tc_bi_ortho/tc_utils.irp.f
@@ -13,6 +13,8 @@ subroutine write_tc_energy()
   ! GS 
   ! ---
 
+  call htilde_mu_mat_opt_bi_ortho(psi_det(1,1,1), psi_det(1,1,1), N_int, hmono, htwoe, hthree, htot)
+
   allocate(E_TC_tmp(N_det), E_1e_tmp(N_det), E_2e_tmp(N_det), E_3e_tmp(N_det))
 
   !$OMP PARALLEL                                                                &
diff --git a/plugins/local/tc_int/compute_tc_int.irp.f b/plugins/local/tc_int/compute_tc_int.irp.f
index e6881f34..60c0d53b 100644
--- a/plugins/local/tc_int/compute_tc_int.irp.f
+++ b/plugins/local/tc_int/compute_tc_int.irp.f
@@ -236,9 +236,13 @@ subroutine provide_int2_grad1_u12_ao()
 
   ! ---
 
+  double precision :: tmp_omp
+
   call wall_time(time1)
 
   PROVIDE ao_integrals_map
+  tmp_omp = get_ao_two_e_integral(1, 1, 1, 1, ao_integrals_map)
+
   !$OMP PARALLEL DEFAULT(NONE)                         &
   !$OMP SHARED(ao_num, tc_int_2e_ao, ao_integrals_map) &
   !$OMP PRIVATE(i, j, k, l)
diff --git a/plugins/local/tc_int/write_tc_int_cuda.irp.f b/plugins/local/tc_int/write_tc_int_cuda.irp.f
index 7d0a0385..212518ee 100644
--- a/plugins/local/tc_int/write_tc_int_cuda.irp.f
+++ b/plugins/local/tc_int/write_tc_int_cuda.irp.f
@@ -136,11 +136,14 @@ subroutine do_work_on_gpu()
 
   integer :: i, j, l
   double precision :: t1, t2
+  double precision :: tmp
   double precision, external :: get_ao_two_e_integral
 
   call wall_time(t1)
 
   PROVIDE ao_integrals_map
+  tmp = get_ao_two_e_integral(1, 1, 1, 1, ao_integrals_map)
+
   !$OMP PARALLEL DEFAULT(NONE)                      &
   !$OMP SHARED(ao_num, int_2e_ao, ao_integrals_map) &
   !$OMP PRIVATE(i, j, k, l)