diff --git a/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h b/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h
index 1882e1e31903f..3bd92b4e11664 100644
--- a/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h
+++ b/DataFormats/simulation/include/SimulationDataFormat/O2DatabasePDG.h
@@ -296,6 +296,19 @@ inline void O2DatabasePDG::addALICEParticles(TDatabasePDG* db)
                     2.5e-15, 6, "Ion", ionCode);
   }
 
+  // Lithium 4 ground state
+  ionCode = 1000030040;
+  if (!db->GetParticle(ionCode)) {
+    db->AddParticle("Lithium4", "Lithium4", 3.74976, kFALSE,
+                    0.005, 9, "Ion", ionCode);
+  }
+  // anti Lithium 4 ground state
+  ionCode = -1000030040;
+  if (!db->GetParticle(ionCode)) {
+    db->AddParticle("AntiLithium4", "AntiLithium4", 3.74976, kFALSE,
+                    0.005, 9, "Ion", ionCode);
+  }
+
   ionCode = 1010020050;
   if (!db->GetParticle(ionCode)) {
     db->AddParticle("Hyperhelium5", "Hyperhelium5", 4.841, kFALSE,
diff --git a/Steer/src/O2MCApplication.cxx b/Steer/src/O2MCApplication.cxx
index 64f79d4c46ac9..118374751ee69 100644
--- a/Steer/src/O2MCApplication.cxx
+++ b/Steer/src/O2MCApplication.cxx
@@ -264,6 +264,11 @@ void addSpecialParticles()
   //Anti-Hyper helium 4 excited state
   TVirtualMC::GetMC()->DefineParticle(-1010020041, "AntiHyperhelium4*", kPTHadron, 3.9231, 2.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
 
+  // Lithium 4 ground state
+  TVirtualMC::GetMC()->DefineParticle(1000030040, "Lithium4", kPTHadron, 3.74976, 3.0, 9.1e-23, "Ion", 0.005, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
+  // Anti-Lithium 4 ground state
+  TVirtualMC::GetMC()->DefineParticle(-1000030040, "AntiLithium4", kPTHadron, 3.74976, 3.0, 9.1e-23, "Ion", 0.005, 0, 1, 0, 0, 0, 0, 0, 4, kFALSE);
+
   //Hyper helium 5
   TVirtualMC::GetMC()->DefineParticle(1010020050, "Hyperhelium5", kPTHadron, 4.841, 2.0, 2.632e-10, "Ion", 0.0, 0, 1, 0, 0, 0, 0, 0, 5, kFALSE);
   //Anti-Hyper helium 5
@@ -597,6 +602,38 @@ void addSpecialParticles()
   //Decay for the excited state (after em transition)
   TVirtualMC::GetMC()->SetDecayMode(-1010020041, abratio4, amode4);
 
+  // Define the 2-body phase space decay for the Lithium 4
+  Int_t model4[6][3];
+  Float_t bratiol4[6];
+
+  for (Int_t kz = 0; kz < 6; kz++) {
+    bratiol4[kz] = 0.;
+    model4[kz][0] = 0;
+    model4[kz][1] = 0;
+    model4[kz][2] = 0;
+  }
+  bratiol4[0] = 100.;
+  model4[0][0] = 1000020030; // Helium3
+  model4[0][1] = 2212;       // proton
+
+  TVirtualMC::GetMC()->SetDecayMode(1000030040, bratiol4, model4);
+
+  // Define the 2-body phase space decay for the Anti-Lithium 4
+  Int_t amodel4[6][3];
+  Float_t abratiol4[6];
+
+  for (Int_t kz = 0; kz < 6; kz++) {
+    abratiol4[kz] = 0.;
+    amodel4[kz][0] = 0;
+    amodel4[kz][1] = 0;
+    amodel4[kz][2] = 0;
+  }
+  abratiol4[0] = 100.;
+  amodel4[0][0] = -1000020030; // Anti-Helium3
+  amodel4[0][1] = -2212;       // Anti-proton
+
+  TVirtualMC::GetMC()->SetDecayMode(-1000030040, abratiol4, amodel4);
+
   // Define the 3-body phase space decay for the Hyper Helium 5
   Int_t mode41[6][3];
   Float_t bratio41[6];