This outcome is different from Eq.(65) in [1]} because of the fact that the time derivative scalar function associated with the vorticity tensor \(\Omega \)  is directly dependent on the \(f(R)\)  terms. Using the above equation along with Eq.(REF ), we acquire
\(\dot{G}=B^2r^2\tilde{g}(t,\theta ),\) 
