diff --git a/src/operators/mpohamiltonian.jl b/src/operators/mpohamiltonian.jl
index 91bdd03d4..a7d876678 100644
--- a/src/operators/mpohamiltonian.jl
+++ b/src/operators/mpohamiltonian.jl
@@ -432,12 +432,12 @@ function Base.:+(H₁::MPOH, H₂::MPOH) where {MPOH<:MPOHamiltonian}
             A₁ = H₁.A[i]
             A₁_inds = CartesianIndices((2:(size(H₁[i], 1) - 1), 1:1, 1:1,
                                         2:(size(H₁[i], 4) - 1)))
-            iszero(A₁) || copyto!(W, A₁_inds, A₁, CartesianIndices(A₁))
+            copyto!(W, A₁_inds, A₁, CartesianIndices(A₁))
 
             A₂ = H₂.A[i]
             A₂_inds = CartesianIndices((size(H₁[i], 1):(size(W, 1) - 1), 1:1, 1:1,
                                         size(H₁[i], 4):(size(W, 4) - 1)))
-            iszero(A₂) || copyto!(W, A₂_inds, A₂, CartesianIndices(A₂))
+            copyto!(W, A₂_inds, A₂, CartesianIndices(A₂))
         end
 
         if H₁ isa InfiniteMPOHamiltonian || i != 1