diff --git a/Resample/Flux/MatrixFlux.cpp b/Resample/Flux/MatrixFlux.cpp
index 3573c2ab73c6c2382e9c6e44bdd09aac59566f54..01b970b8b005d52a3dcc2d45adedc61a78e178ad 100644
--- a/Resample/Flux/MatrixFlux.cpp
+++ b/Resample/Flux/MatrixFlux.cpp
@@ -34,7 +34,7 @@ SpinMatrix PhaseRotationMatrix(const R3& B_direction, const Spinor& diagonal)
 
     if (std::abs(e.mag() - 1.) < eps) {
         const double factor1 = 2. * (1. + e.z());
-        SpinMatrix Q(1. + e.z(), -e.x() + I * e.y(), e.x() + I * e.y(), e.z() + 1.);
+        SpinMatrix Q(1. + e.z(), e.x() - I * e.y(), e.x() + I * e.y(), -e.z() - 1.);
         return Q * M * Q.adjoint() / factor1;
     }
     ASSERT(e.mag() < eps); // remaining case: no field