diff --git a/src/REPLICA/fix_pimd_langevin.cpp b/src/REPLICA/fix_pimd_langevin.cpp
index 2e4b1c17b8..17ea07c845 100644
--- a/src/REPLICA/fix_pimd_langevin.cpp
+++ b/src/REPLICA/fix_pimd_langevin.cpp
@@ -451,7 +451,7 @@ void FixPIMDLangevin::init()
 
   nmpimd_init();
 
-  Langevin_init();
+  langevin_init();
 
   c_pe = modify->get_compute_by_id(id_pe);
   c_press = modify->get_compute_by_id(id_press);
@@ -888,7 +888,7 @@ void FixPIMDLangevin::press_o_step()
 
 /* ---------------------------------------------------------------------- */
 
-void FixPIMDLangevin::Langevin_init()
+void FixPIMDLangevin::langevin_init()
 {
   double beta = 1.0 / kBT;
   _omega_np = np / beta / hbar;
diff --git a/src/REPLICA/fix_pimd_langevin.h b/src/REPLICA/fix_pimd_langevin.h
index 79a2ee7477..1ef1dfae52 100644
--- a/src/REPLICA/fix_pimd_langevin.h
+++ b/src/REPLICA/fix_pimd_langevin.h
@@ -115,7 +115,7 @@ class FixPIMDLangevin : public Fix {
   class RanMars *random;
 
   int tstat_flag;    // tstat_flat = 1 if thermostat if used
-  void Langevin_init();
+  void langevin_init();
   void b_step();    // integrate for dt/2 according to B part (v <- v + f * dt/2)
   void a_step();    // integrate for dt/2 according to A part (non-centroid mode, harmonic force between replicas)
   void qc_step();    // integrate for dt/2 for the centroid mode (x <- x + v * dt/2)