diff --git a/integrators/__init__.py b/integrators/__init__.py
index f8d59cd4e7129d90690c5235df49aea1501d6a69..2a28c5e3784850fd8d0ffb9852983200f4d0fa03 100644
--- a/integrators/__init__.py
+++ b/integrators/__init__.py
@@ -27,6 +27,8 @@ from commics.integrators.kw1_mps import KW1_MPS
 from commics.integrators.tps1_kw2 import TPS1_KW2
 from commics.integrators.tps1_mps import TPS1_MPS
 
+from commics.integrators.tps2ab_initial_tps2ee import TPS2AB_INITIAL_TPS2EE
+
 __all__ = \
   [ "Integrator" \
   , "DUMMY" \
@@ -36,4 +38,5 @@ __all__ = \
   , "TPS1HEE" \
   , "TPS1HEE_MPS" \
   , "KW1_MPS", "TPS1_KW2", "TPS1_MPS" \
+  , "TPS2AB_INITIAL_TPS2EE"
   ]
diff --git a/integrators/tps2ab_initial_tps2ee.py b/integrators/tps2ab_initial_tps2ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..2cb36bab1846262a8f69186e81a61974016d1411
--- /dev/null
+++ b/integrators/tps2ab_initial_tps2ee.py
@@ -0,0 +1,25 @@
+###  TPS2+AB + TPS2EE for first time step ###
+from . import tps2ab
+
+
+################################################################################
+
+
+class TPS2AB_INITIAL_TPS2EE(tps2ab.TPS2AB):
+  '''
+  implements TPS2AB from second time-step on,
+  BUT: At first time step use TPS2EE = EE approximation of pi, i.e., pi(v) := 0
+  '''
+
+
+#------------------------------------------------------------------------------#
+
+
+  def _DerivLinearTerms_pi_Implicit_RHS(self):
+    return self._ZeroCF()
+
+
+#------------------------------------------------------------------------------#
+
+
+################################################################################