diff --git a/src/apps/molresponse/FrequencyResponse.cpp b/src/apps/molresponse/FrequencyResponse.cpp
index c7d59b640de..624ce28f92f 100644
--- a/src/apps/molresponse/FrequencyResponse.cpp
+++ b/src/apps/molresponse/FrequencyResponse.cpp
@@ -539,8 +539,6 @@ auto QuadraticResponse::setup_XBC(World &world) -> std::pair<X_space, X_space> {
     return {r_XB, r_XC};
 }
 
-X_space QuadraticResponse::compute_gamma_bc_virt(World &world, const X_space &XB, const X_space &XC) {}
-
 Tensor<double> QuadraticResponse::compute_beta_unrelaxed(World &world, const X_space &AB_left, const X_space &AB_right,
                                                          X_space &BA_left, X_space &BA_right) {
 
diff --git a/src/apps/molresponse/FrequencyResponse.hpp b/src/apps/molresponse/FrequencyResponse.hpp
index c2ca3bbd520..c8356678c6e 100644
--- a/src/apps/molresponse/FrequencyResponse.hpp
+++ b/src/apps/molresponse/FrequencyResponse.hpp
@@ -206,7 +206,6 @@ class QuadraticResponse : public ResponseBase {
     std::pair<X_space, X_space> dipole_perturbation(World &world, const X_space &left, const X_space &right) const;
     Tensor<double> compute_beta_1(const X_space &upsilon);
     Tensor<double> compute_beta_1(World &world, const X_space &upsilon);
-    X_space compute_gamma_bc_virt(World &world, const X_space &VB, const X_space &XC);
     X_space compute_gamma_occ(World &world, const X_space &XB, const X_space &XC);
     X_space compute_g1_term(World &world, const X_space &left, const X_space &right, const X_space &apply) const;
     X_space compute_coulomb_term(World &world, const X_space &B, const X_space &C, const X_space &D) const;