diff --git a/src/thermophysicalModels/specie/reaction/Reactions/Reaction/Reaction.H b/src/thermophysicalModels/specie/reaction/Reactions/Reaction/Reaction.H
index a2ea0e3e..9489ac31 100644
--- a/src/thermophysicalModels/specie/reaction/Reactions/Reaction/Reaction.H
+++ b/src/thermophysicalModels/specie/reaction/Reactions/Reaction/Reaction.H
@@ -141,10 +141,10 @@ private:
         List<specieCoeffs> rhs_;
 
         //- Electron temperature parameter
-        scalar Te_;
+        mutable scalar Te_;
 
         //- Reduced electric field parameter
-        scalar En_;
+        mutable scalar En_;
 
 
     // Private Member Functions
@@ -310,9 +310,11 @@ public:
             inline word& name();
             inline const word& name() const;
 
+            inline scalar Te(const scalar newTe) const;
             inline scalar& Te();
             inline scalar Te() const;
 
+            inline scalar En(const scalar newEn) const;
             inline scalar& En();
             inline scalar En() const;
 
diff --git a/src/thermophysicalModels/specie/reaction/Reactions/Reaction/ReactionI.H b/src/thermophysicalModels/specie/reaction/Reactions/Reaction/ReactionI.H
index 1448a68f..0d23cf51 100644
--- a/src/thermophysicalModels/specie/reaction/Reactions/Reaction/ReactionI.H
+++ b/src/thermophysicalModels/specie/reaction/Reactions/Reaction/ReactionI.H
@@ -46,6 +46,14 @@ inline const word& Reaction<ReactionThermo>::name() const
 }
 
 
+template<class ReactionThermo>
+inline scalar Reaction<ReactionThermo>::Te(const scalar newTe) const
+{
+    Te_ = newTe;
+    return Te_;
+}
+
+
 template<class ReactionThermo>
 inline scalar& Reaction<ReactionThermo>::Te()
 {
@@ -60,6 +68,14 @@ inline scalar Reaction<ReactionThermo>::Te() const
 }
 
 
+template<class ReactionThermo>
+inline scalar Reaction<ReactionThermo>::En(const scalar newEn) const
+{
+    En_ = newEn;
+    return En_;
+}
+
+
 template<class ReactionThermo>
 inline scalar& Reaction<ReactionThermo>::En()
 {