diff --git a/kwant/builder.py b/kwant/builder.py
index abf1cc69541ffc5964771085d6b5d6d3d13b90bf..9a0bfcffbe6a03e01af674e3f29a0173b35d89b1 100644
--- a/kwant/builder.py
+++ b/kwant/builder.py
@@ -1875,6 +1875,10 @@ def add_peierls_phase(syst, peierls_parameter='phi', fix_gauge=True):
 
         return f
 
+    if syst.vectorize:
+        raise TypeError("'add_peierls_phase' does not work with "
+                        "vectorized Builders")
+
     ret = _add_peierls_phase(syst, peierls_parameter).finalized()
 
     if fix_gauge:
diff --git a/kwant/wraparound.py b/kwant/wraparound.py
index 08efec8c668567acb374b12bfbed7777f0571bb5..e5688bab42919e664124fd73d8da0bdf7bb52ae6 100644
--- a/kwant/wraparound.py
+++ b/kwant/wraparound.py
@@ -183,6 +183,9 @@ def wraparound(builder, keep=None, *, coordinate_names='xyz'):
         f.__signature__ = inspect.Signature(params.values())
         return f
 
+    if builder.vectorize:
+        raise TypeError("'wraparound' does not work with vectorized Builders.")
+
     try:
         momenta = ['k_{}'.format(coordinate_names[i])
                    for i in range(len(builder.symmetry.periods))]