Optimal feedback control for a class of fed-batch fermentation processes using switched dynamical system approach

Xiang Wu; Yuzhou Hou; Kanjian Zhang; Xiang Wu; Yuzhou Hou; Kanjian Zhang

doi:10.3934/math.2022510

AIMS Mathematics

2022, Volume 7, Issue 5: 9206-9231. doi: 10.3934/math.2022510

Previous Article Next Article

Research article

Optimal feedback control for a class of fed-batch fermentation processes using switched dynamical system approach

1.
School of Mathematical Sciences, Guizhou Normal University, Guiyang 550001, China
2.
School of life sciences, Guizhou Normal University, Guiyang 550001, China
3.
School of Electrical Engineering, Southeast University, Nanjing 210096, China
4.
School of Automation, Southeast University, Nanjing 210096, China
5.
Key Laboratory of Measurement and Control of CSE, Ministry of Education, Southeast University, Nanjing 210096, China

Received: 29 October 2021 Revised: 22 February 2022 Accepted: 28 February 2022 Published: 09 March 2022

This paper considers an optimal feedback control problem for a class of fed-batch fermentation processes. Our main contributions are as follows. Firstly, a dynamic optimization problem for fed-batch fermentation processes is modeled as an optimal control problem of switched dynamical systems, and a general state-feedback controller is designed for this dynamic optimization problem. Unlike the existing switched dynamical system optimal control problem, the state-dependent switching method is applied to design the switching rule, and the structure of this state-feedback controller is not restricted to a particular form. Then, this problem is transformed into a mixed-integer optimal control problem by introducing a discrete-valued function. Furthermore, each of these discrete variables is represented by using a set of 0-1 variables. By using a quadratic constraint, these 0-1 variables are relaxed such that they are continuous on the closed interval $[0, 1]$ . Accordingly, the original mixed-integer optimal control problem is transformed intoa nonlinear parameter optimization problem. Unlike the existing works, the constraint introduced for these 0-1 variables are at most quadratic. Thus, it does not increase the number of locally optimal solutions of the original problem. Next, an improved gradient-based algorithm is developed based on a novel search approach, and a large number of numerical experiments show that this novel search approach can effectively improve the convergence speed of this algorithm, when an iteration is trapped to a curved narrow valley bottom of the objective function. Finally, numerical results illustrate the effectiveness of this method developed by this paper.

Keywords:

Citation: Xiang Wu, Yuzhou Hou, Kanjian Zhang. Optimal feedback control for a class of fed-batch fermentation processes using switched dynamical system approach[J]. AIMS Mathematics, 2022, 7(5): 9206-9231. doi: 10.3934/math.2022510

Related Papers:

[1]	Xiaoling Zhang, Cuiling Ma, Lili Zhao . On some $m$ -th root metrics. AIMS Mathematics, 2024, 9(9): 23971-23978. doi: 10.3934/math.20241165
[2]	Rajesh Kumar, Sameh Shenawy, Lalnunenga Colney, Nasser Bin Turki . Certain results on tangent bundle endowed with generalized Tanaka Webster connection (GTWC) on Kenmotsu manifolds. AIMS Mathematics, 2024, 9(11): 30364-30383. doi: 10.3934/math.20241465
[3]	Aliya Naaz Siddiqui, Mohammad Hasan Shahid, Jae Won Lee . On Ricci curvature of submanifolds in statistical manifolds of constant (quasi-constant) curvature. AIMS Mathematics, 2020, 5(4): 3495-3509. doi: 10.3934/math.2020227
[4]	Fahad Sikander, Tanveer Fatima, Sharief Deshmukh, Ayman Elsharkawy . Curvature analysis of concircular trajectories in doubly warped product manifolds. AIMS Mathematics, 2024, 9(8): 21940-21951. doi: 10.3934/math.20241066
[5]	Salah G. Elgendi . Parallel one forms on special Finsler manifolds. AIMS Mathematics, 2024, 9(12): 34356-34371. doi: 10.3934/math.20241636
[6]	Ibrahim Al-Dayel, Meraj Ali Khan . Ricci curvature of contact CR-warped product submanifolds in generalized Sasakian space forms admitting nearly Sasakian structure. AIMS Mathematics, 2021, 6(3): 2132-2151. doi: 10.3934/math.2021130
[7]	Tong Wu, Yong Wang . Super warped products with a semi-symmetric non-metric connection. AIMS Mathematics, 2022, 7(6): 10534-10553. doi: 10.3934/math.2022587
[8]	Fangmin Dong, Benling Li . On a class of weakly Landsberg metrics composed by a Riemannian metric and a conformal 1-form. AIMS Mathematics, 2023, 8(11): 27328-27346. doi: 10.3934/math.20231398
[9]	Ali H. Alkhaldi, Meraj Ali Khan, Shyamal Kumar Hui, Pradip Mandal . Ricci curvature of semi-slant warped product submanifolds in generalized complex space forms. AIMS Mathematics, 2022, 7(4): 7069-7092. doi: 10.3934/math.2022394
[10]	Songting Yin . Some rigidity theorems on Finsler manifolds. AIMS Mathematics, 2021, 6(3): 3025-3036. doi: 10.3934/math.2021184

Abstract

1. Introduction

Finsler geometry extends the classical Riemannian geometry by considering more general metric structures. A very important class of Finsler metrics is known as $(\alpha, \beta)$ -metrics, which were introduced by M. Matsumoto in $1972$ . An $(\alpha, \beta)$ -metric can be expressed as $F = \alpha \phi(s)$ , where $\alpha$ is a Riemannian metric and $s = \frac{\beta}{\alpha}$ , $\beta$ is a $1$ -form. Randers metric, Kropina metric, exponential metric, Matsumoto metric, and cubic metric are important classes of $(\alpha, \beta)$ -metric ^[9].

To study the curvature characteristics is a central problem in Finsler geometry. The Ricci curvature and S-curvature are very important non-Riemannian quantities in the Finslerian manifold ^[2]. The Ricci curvature in Finsler geometry is a natural extension of the Ricci curvature in Riemannian geometry and is defined as the trace of the Riemann curvature ^[5]. The S-curvature is a mathematical quantity and measures the rate of change of volume form of a Finsler space along the geodesics. Recent studies in differential geometry, such as those on Ricci solitons and conformal structures, have highlighted the importance of Ricci-type curvatures in understanding the geometric flow and structure of manifolds ^[6,7,8]. In Finsler geometry, the study of curvature involves understanding the deviation from flatness. The projective Ricci curvature is one aspect of this analysis. The concept of projective Ricci curvature in Finsler geometry is introduced by X. Cheng ^[1] in $2017$ . Projective geometry deals with the properties that are invariant under projective transformations. The projective Ricci curvature measures the deviation of the Finsler metric from being projectively flat. Projective Ricci curvature has applications in various areas of mathematics and physics. It plays a crucial role in understanding the geometry of Finsler manifolds and connects to the problems in the calculus of variations, differential equations, and geometric optics.

In $2020$ , H. Zhu ^[15] gave an expression of projective Ricci curvature for an $(\alpha, \beta)$ -metric. Later on, many geometers ^[4,12,13] have studied the geometric properties of projective Ricci curvature. In this article, we obtain the geometric properties and flatness condition of projective Ricci curvature for the cubic Finsler metric, which is defined as $F = \alpha \phi(s)$ with

$\begin{equation} \phi = (1+s)^3, \end{equation}$

(1.1)

i.e., $F = \frac{(\alpha+\beta)^3}{\alpha^2}$ . Cubic metric is Finsler metric for $b^2 < \frac{1}{4}$ ^[14].

The following notations will be used to state our main result:

$\begin{equation} \begin{split} &2s^{}_{jk} = b^{}_{j;k}-\! b^{}_{k;j}, \quad 2r^{}_{jk} = b^{}_{j;k}\!+\! b^{}_{k;j},\quad s^{j}_{k} = a^{jl}s^{}_{kl},\quad r^{j}_{k} = a^{jl}r^{}_{kl},\quad s^{}_{j} = b^{l}s^{}_{lj} = b^{}_{k}s^{k}_{j},\\& r^{}_{j} = b^{l}r^{}_{lj} = b^{}_{k}r^{k}_{j},\quad r^{}_{j0} = r^{}_{jk}y^{k},\quad r^{}_{00} = r^{}_{jk}y^{j}y^{k},\quad r = r^{}_{jk}b^{j}b^{k} = b^{j}r^{}_{j},\quad s^{}_{j0} = s^{}_{jk}y^{k}, \\& s^{}_{0} = s^{}_{j}y^{j}, \quad r^{}_{0} = r^{}_{j}y^{j},\quad b^{j} = a^{jk}b^{}_{k},\quad t^{}_{jk} = s^{}_{jm}s^{m}_{k}, \quad t^{}_{j} = b^{m}t^{}_{mj} = s^{i}_{}s^{i}_{j}, \end{split} \end{equation}$

(1.2)

where " $;$ " denotes the covariant derivative with respect to the Levi-Civita connection of the Riemannian metric $\alpha$ .

A $1$ -form $\beta$ is said to be a Killing form if $r^{}_{ij} = 0$ . The $1$ -form $\beta$ is said to be a constant Killing form if it is a Killing form and constant length concerning $\alpha$ , equivalently $r^{}_{ij} = 0$ and $s^{}_{i} = 0$ .

In this paper we will use the following lemma:

Lemma 1.1. If $\alpha^2 = 0 (mod \beta)$ , that is, $a_{ij}y^{i}y^{j}$ contains $b^{}_i(x)y^i$ as a factor, then the dimension is equal to two and $b^2$ vanishes. In this case, we have $\delta = d^{}_i(x)y^i$ satisfying $\alpha^2 = \beta\delta$ and $d^{}_ib^i = 2$ .

We first prove the following result:

Theorem 1.1. For the cubic Finsler metric $F = \frac{(\alpha+\beta)^3}{\alpha^2}$ on an $n$ -dimensional $(n > 2)$ Finsler manifold $M$ , the S-curvature vanishes if and only if $\beta$ is a constant Killing form.

Next, we obtain the flatness condition for the projective Ricci curvature as

Theorem 1.2. If the $n$ -dimensional $(n > 2)$ Finsler space with cubic metric $F = \frac{(\alpha+\beta)^3}{\alpha^2}$ is projective Ricci-flat $(PRic = 0),$ then $\beta$ is parallel with respect to the Riemannian metric $\alpha$ .

In view of the above result, we obtain

Corollary 1.1. If the $n$ -dimensional $(n > 2)$ Finsler space with cubic metric $F = \frac{(\alpha+\beta)^3}{\alpha^2}$ is projective Ricci flat then, it vanishes the S-curvature. Therefore, the Riemannian metric of $\alpha$ is Ricci flat $(Ric^{\alpha}$ $=$ $0)$ .

We also prove the following result:

Theorem 1.3. The $n$ -dimensional $(n > 2)$ Finsler space with cubic metric $F = \frac{(\alpha+\beta)^3}{\alpha^2}$ is weak PRic-curvature if and only if it is a PRic-flat metric.

2. Preliminaries

Let $F$ be an $n$ -dimensional Finsler manifold, and let $G^{j}$ be the geodesic coefficients of $F$ , which are defined as

$\begin{equation*} G^{j} = \frac{1}{4}g^{jl}\Big[\frac{\partial^2(F^2)}{\partial x^{k}\partial y^{l}}y^{k}-\frac{\partial(F^2)}{\partial x^{l}}\Big], \qquad y\epsilon T^{}_{x}M. \end{equation*}$

The geodesic coefficients of an $(\alpha, \beta)$ -metric are given as ^[3]

$\begin{equation} G^{j} = G^{j}_{\alpha}+\alpha Q s^{j}_{0}+(r^{}_{00}-2\alpha Q s^{}_{0})(\Psi b^{j}+\frac{\Theta y^{j}}{\alpha}), \end{equation}$

(2.1)

where $G^{i}_{\alpha}$ denotes the geodesic coefficients of the Riemannian metric $\alpha$ and

$\begin{equation} \begin{split} &Q = \frac{\phi^{'}}{\phi-s\phi^{'}},\quad \psi = \frac{\phi^{"}}{2\phi(\phi-s\phi^{'}+(B-s^{2})\phi^{"})},\quad \Theta = \frac{\phi\phi^{'}-s(\phi\phi^{"}+\phi^{'}\phi^{'})}{2\phi(\phi-s\phi^{'}+(B-s^{2})\phi^{"})}. \end{split} \end{equation}$

(2.2)

For any $x\epsilon M$ and $y\epsilon T^{}_{x}M\backslash\{0\}$ , the Riemann curvature $R^{}_{y}$ is defined as

$\begin{equation*} R^{}_{y}(v) = R^{j}_{k}(y)v^{k}\frac{\partial}{\partial x^{j}} , v = v^{j}\frac{\partial}{\partial x^{j}} , \end{equation*}$

where

$\begin{equation*} R^{j}_{k} = 2\frac{\partial G^{j}}{\partial x^{k}}+ 2G^{i}\frac{\partial^{2} G^{j}}{\partial y^{i}\partial y^{k}}- \frac{\partial^{2} G^{j}}{\partial x^{i}\partial y^{k}} y^{i}-\frac{\partial G^{j}}{\partial y^{i}}\frac{\partial G^{i}}{\partial y^{k}}. \end{equation*}$

The trace of Riemann curvature is called Ricci curvature $Ric = R^{m}_{m}$ , which is a mathematical object that regulates the rate at which a metric ball's volume in a manifold grows. A Finsler metric $F$ is called an Einstein metric if Ricci curvature satisfies the equation $Ric(x, y) = (n-1)\gamma F^{2}$ , where $\gamma = \gamma(x)$ is a scalar function.

In $1997$ , Z. Shen ^[11] discussed S-curvature, which measures the average rate of change of $(T_{x}M; F_{x})$ in the direction $y \epsilon T^{}_{x}M$ and is defined as

$\begin{equation*} S(x,y) = \frac{\partial G^{m}}{\partial y^{m}}-y^{m}\frac{\partial (log \sigma^{}_{F})}{\partial x^{m}} , \end{equation*}$

where $\sigma ^{}_{F}$ is defined as

$\begin{equation*} \sigma_{F} = \frac{Vol(B^{n})}{Vol\{y^{i} \epsilon \mathbb{R}^{n}| F( x,y) < 1\}} , \end{equation*}$

and Vol denotes the Euclidean volume, and $B^{n}(1)$ denotes the unit ball in $R^{n}$ .

The expression of S-curvature for an $(\alpha, \beta)$ -metric is given as ^[10]

$\begin{equation} S = (s_{0}+r_{0})\left(2\psi- \Pi\right)-\alpha^{-1}\frac{\Phi}{2\Delta^{2}}(r_{00}-2\alpha Q s_{0}), \end{equation}$

(2.3)

where

$\begin{equation} \begin{split} &\Pi = \frac{f^{'}(b)}{bf(b)} , \quad \Delta = 1+sQ+(B-s^2)Q^{}_{s} ,\quad B = b^{2} , \\& \Phi = -(Q-sQ^{}_{s})(n\Delta+1+sQ)-(B-s^2)(1+sQ) Q^{}_{ss} . \end{split} \end{equation}$

(2.4)

The projective Ricci curvature is first defined by X. Cheng ^[1] as

$\begin{equation} PRic = Ric+\frac{n-1}{n+1}S^{}_{|m}y^{m}+\frac{n-1}{(n+1)^2}S^{2}, \end{equation}$

(2.5)

where " $\vert$ " denotes the horizontal covariant derivative with respect to the Berwald connections of $F$ . A Finsler space $F$ is called weak projective Ricci curvature if

$\begin{equation} PRic = (n-1)\Big[\frac{3 \theta}{F}+\gamma\Big]F^{2}, \end{equation}$

(2.6)

where $\gamma = \gamma(x)$ is a scalar function and $\theta = \theta_{i}(x)y^{i}$ is a $1$ -form. If $\gamma =$ constant, then $F$ is called constant projective Ricci curvature. If $\theta = 0$ , then $F$ is called isotropic projective Ricci curvature PRic = $(n-1)\gamma F^{2}$ .

In $2020$ , H. Zhu ^[15] gave an expression of the projective Ricci curvature for the $(\alpha, \beta)$ -metrics as

$\begin{equation} \begin{split} PRic\! = \! Ric^{\alpha}&\!+\! \frac{1}{n\!+\!1}\Big[\frac{r^{2}_{00}}{\alpha^2}V_{1}-\frac{r^{}_{00}s^{}_{0}}{\alpha}V_{2}-\frac{r^{}_{00}r^{}_{0}}{\alpha}V_{3}+\frac{r^{}_{00|0}}{\alpha}V_{4}+s^{2}_{0}V_{5}+r^{}_{00}rV_{6}-4r^{2}_{0}V_{7}+2r^{}_{0}s^{}_{0}V_{8}\Big]\\&\!+\!(r^{}_{00}r^{i}_{i}\!+\!r^{}_{00|b}\!-\!b^{i}r^{}_{i0|0}\!-\!r^{}_{0i}r^{i}_{0})V_{9}\!+\!r^{}_{0i}s^{i}_{0}V_{10}\!+\!s^{}_{0|0}V_{11}\!+\!s^{}_{0i}s^{i}_{0}V_{12}\!+\!\alpha rs^{}_{0}V_{13}\!+\!\alpha s^{}_{j}s^{j}_{0}V_{14}\\& \!+\!\alpha\Big[\frac{4}{n\!+\!1}r^{}_{i}s^{i}_{0}\!-\!2s^{}_{0|b}\!-\!2r^{i}_{i}s^{}_{0}\!+\!3r^{}_{0i}s^{i}\!+\!b^{i}s^{i}_{|0}\Big]V_{15}\!+\!\alpha s^{i}_{0|i}V_{16}\!+\!\alpha^2 s^{}_{i}s^{i}V_{17}\!+\!\alpha^2 s^{i}_{j}s^{j}_{i}V_{18}\\& +r^{}_{0|0}V_{19}+\frac{2(n-1)}{n+1}\Big[\frac{\Psi^{}_{s}}{\alpha}(r^{}_{00}-2\alpha Q s^{}_{0})(B-s^2)+2 \Psi (r^{}_{0}+s^{}_{0})\Big]\rho^{}_{0}+(n-1)\Big[-2\Psi(r^{}_{00}\\& -2\alpha Q s^{}_{0})\rho^{}_{b}-2\alpha Q \rho^{}_{k}s^{k}_{0}+\rho^{2}_{0}+\rho^{}_{0|0}\Big], \end{split} \end{equation}$

(2.7)

where

$\begin{equation*} \begin{split} &V_{1} = 4s\Psi^{}_{s}+(4 \Psi \Psi^{}_{ss}-\Psi^{2}_{s})(B-s^2)^2-2(\Psi^{}_{ss}+6s\Psi\Psi^{}_{s})(B-s^2)-\frac{n^{2}+1}{n+1}\Psi^{2}_{s}(B-s^2)^2,\\ &V_{2} = 4\big[2\Psi(\Psi Q^{}_{ss}+2Q\Psi^{}_{ss}+Q^{}_{s}\Psi^{}_{s})-Q(\Psi^{}_{s})^2\big](B-s^2)^2+4\big[2(Q-sQ^{}_{s})(\Psi^{}_{s})^2-(1+10sQ)\Psi\Psi_{s}\\&\; \; \; \; \; \; \; -2\Psi Q^{}_{ss}-2Q\Psi^{}_{ss}-Q^{}_{s}\Psi_{s}+\Psi_{Bs}\big](B-s^2)+2Q^{}_{ss}+8\Psi_{s}-4Q\Psi+4s\Psi Q_{s}+20s Q \Psi_{s}\\&\; \; \; \; \; \; \; +(n-1)\big[4((\Psi)^2Q_{ss}-Q(\Psi_{s})^2)(B-s^2)^2+4((Q-sQ_{s})(\Psi)^2)+\Psi(\Psi_{s}-Q_{ss})(B-s^2)\\&\; \; \; \; \; \; \; +2s(Q_{s}\Psi+Q\Psi_{s})-2Q\Psi+Q_{ss}\big]+\frac{8}{n+1}\Psi_{s}\big[\Psi-Q\Psi_{s}(B-s^2)\big](B-s^2),\\& V_{3} = 2\Psi_{s}-2(3\Psi \Psi_{s}-\Psi_{Bs})(B-s^2)-(n-1)(1-2\Psi(B-s^2))\Psi_{s}+\frac{4}{n+1}\Psi \Psi_{s}(B-s^2),\\& V_{4} = -2\Psi_{s}(B-s^2),\\ &V_{5} = (n-1)\Big\{4\big[2Q\Psi^2Q_{ss}-\Psi^2(Q_{s})^2-Q^2(\Psi_{s})^2\big](B-s^2)^2+4[2Q\Psi(Q\Psi-2sQ_{s}\Psi-Q_{ss}+\Psi_{s})\\&\; \; \; \; \; \; \; +\Psi(Q_{s})^2+Q_{s}\Psi_{B}](B-s^2)-4Q\Psi(s^2Q\Psi-3sQ_{s}+2Q)+4sQ(Q\Psi_{s}+\Psi_{B})+8Q_{s}\Psi\\&\; \; \; \; \; \; \; +2QQ_{ss}-(Q_{s})^2+4\Psi_{B}\Big\}+4\big[4Q\Psi(\Psi Q_{ss}+Q\Psi_{ss}+Q_{s}\Psi_{s})-2(Q_{s})^2\Psi^2-Q^2(\Psi_{s})^2\big](B-s^2)^2\\&\; \; \; \; \; \; \; +8\big[Q\Psi(-4sQ_{s}\Psi-4sQ\Psi_{s}+2Q\Psi-2Q_{ss}-\Psi_{s})+(Q_{s})^2\Psi-Q^2\Psi_{ss}-QQ_{s}\Psi_{s}+Q\Psi_{Bs}\\&\; \; \; \; \; \; \; +Q_{s}\Psi_{B}\big](B-s^2)+24sQ^2\Psi_{s}- 8Q\Psi(s^2Q\Psi-3sQ_{s}+2Q)+8sQ\Psi_{B}+4\Psi(\Psi+4Q_{s})+4QQ_{ss}\\&\; \; \; \; \; \; \; +16Q\Psi_{s}-2(Q_{s})^2-\frac{8}{n+1}\big[\Psi-Q\Psi_{s}(B-s^2)\big]^2, \end{split} \end{equation*}$

$\begin{equation} \begin{split} & V_{6} = 4\big[(n+1)\Psi_{B}+2(\Psi)^2\big],\qquad V_{7} = \frac{n^2+n+2}{n+1}\Psi^2+2\Psi_{B},\\ &V_{8} = (n-1)\Big[2(2Q_{s}\Psi^2+2Q\Psi\Psi_{s}+Q_{s}\Psi_{B})(B-s^2)+2sQ(2\Psi^2+\Psi_{B})-2Q\Psi_{s}+2\Psi_{B}\Big]\\&\; \; \; \; \; \; \; +4\big[2Q_{s}\Psi^2-3Q\Psi\Psi_{s}+Q\Psi_{Bs}+Q_{s}\Psi_{B}\big](B-s^2)+4sQ(2\Psi^2+\Psi_{B})+4\Psi^2+4Q\Psi_{s}\\&\; \; \; \; \; \; \; -4\Psi_{B}-\frac{8}{n+1}\Psi\big[\Psi-Q\Psi_{s}(B-s^2)\big],\\& V_{9} = 2\Psi , \qquad V_{10} = 2\big[-2Q_{s}\Psi(B-s^2)-2sQ\Psi-\Psi+Q_{s}+\frac{4}{n+1}Q \Psi_{s}(B-s^2)\big],\\& V_{11} = 2Q_{s}\Psi(B-s^2)+2\Psi(1+2Q)-Q_{s}+\frac{4}{n+1}\big[Q\Psi_{s}(B-s^2)-\Psi\big],\\& V_{12} = 2Q_{s}-2Q(Q-sQ_{s}), \qquad V_{13} = -8Q[\Psi_{B}+\frac{2}{n+1}\Psi^2],\\& V_{14} = \frac{-2}{n+1}Q\big[(n-3)\Psi+4Q\Psi_{s}(B-s^2)\big], \qquad V_{15} = 2Q\Psi, \qquad V_{16} = 2Q,\\& V_{17} = -4Q^2\Psi, \qquad V_{18} = -Q^2, \qquad V_{19} = \frac{ 2(n-1)}{n+1}\Psi,\qquad \rho = \frac{ln \frac{\sigma^{}_{\alpha}}{\sigma}}{n+1}, \qquad \rho^{}_{0} = \rho^{}_{x^i}y^{i}. \end{split} \end{equation}$

(2.8)

3. S-curvature of cubic metrics

For Eq (1.1), we obtain the following values:

$\begin{equation} \begin{split} &Q = \frac{3}{1-2s}, \quad Q_{s} = \frac{6}{(1-2s)^2}, \quad Q_{ss} = \frac{24}{(1-2s)^{3}},\quad \psi = \frac{3}{1+6B-s-8s^2},\\& \psi_{s} = \frac{3+48s}{(1 + 6 B - s - 8s^2)^2},\quad \psi_{ss} = \frac{18 (3 + 16 B + 8 s + 64 s^2)}{(1 + 6 B - s - 8 s^2)^3}, \quad \psi_{B} = \frac{-18}{(1 + 6 B - s - 8 s^2)^2},\\& \psi_{Bs} = \frac{36 + 576 s}{(1 + 6 B - s - 8 s^2)^2}, \quad \Theta = \frac{3(1-4s))}{2(1 + 6 B - s - 8 s^2)}, \quad \Theta_{s} = -\frac{9 + 72 B - 48 s + 96 s^2}{2(1 + 6 B - s - 8 s^2)^2},\\ & \Theta_{B} = \frac{9 (-1 + 4 s)}{2(1 + 6 B - s - 8 s^2)^2}, \quad \Delta = \frac{1 + 6 B - s - 8 s^2}{(1 - 2 s)^2},\\& \Phi = \frac{-(3 (1 - 5 s - 6 s^2 + B (8 + 6 n + 8 s - 24 n s) + n (1 - 5 s - 4 s^2 + 32 s^3)}{(1 - 2 s)^4}. \end{split} \end{equation}$

(3.1)

By using Eqs (2.1) and (3.1), we obtain the spray coefficient $G^{j}$ for the cubic metric as

$\begin{equation} \begin{split} &G^{j} = G^{j}_{\alpha}+ \frac{1}{2\alpha(1-2s)(1 +6 B-s-8 s^2)}\Big[ (6 + 36 B-6 s - 48 s^2) \alpha^{2}s^{j}_{0} + [18 \alpha s^{}_{0}(4s-1)\\&\; \; \; \; \; \; \; \; + 3 r^{}_{00} (1-6s+8s^{2})]y^j - 6\alpha b^j [6 s^{}_{0}+(2s-1)r^{}_{00}]\Big]. \end{split} \end{equation}$

(3.2)

In view of Eqs (2.5) and (3.1) and using Mathematica program, we obtain the S-curvature for the cubic Finsler metric as

$\begin{equation} \begin{split} &S = \frac{1}{2(\alpha - 2 \beta)(\alpha^2+6B\alpha^2-\alpha \beta-8 \beta^2)^2}\Big[-2r^{}_{0} (\alpha - 2 \beta) ((1 + 6 B)\alpha^2 - \alpha \beta - 8 \beta^2) [-\alpha \beta \Pi - 8 \beta^2 \Pi \\&\; \; \; \; \; \; \; + \alpha^2 (-6 + \Pi + 6 B \Pi)] -2s^{}_{0} [3 \alpha^2 ((1 + 3 n + 6 B (2 + 3 n)) \alpha^3 - 3 (3 + 5 n + 8 B (-2 + 3 n)) \alpha^2 \beta \\&\; \; \; \; \; \; \; - 6 (1\! + \!2 n) \alpha \beta^2 \!+\! 32 (-1 \!+\! 3 n) \beta^3) \!+\! (\alpha \!-\! 2 \beta) (-(1 \!+\! 6 B) \alpha^2 \!+\! \alpha \beta\! +\! 8 \beta^2)^2 \Pi]\!+3r^{}_{00} (\alpha\! - \! 2 \beta) [(1 \!+\! n \\&\; \; \; \; \; \; \; + B (8 + 6 n)) \alpha^3- (5 - 8 B + 5 n + 24 B n) \alpha^2 \beta - 2 (3 + 2 n) \alpha \beta^2 + 32 n \beta^3]\Big]. \end{split} \end{equation}$

(3.3)

Now, we are in the position to prove Theorem 1.1.

Proof of Theorem 1.1. First we prove the converse part.

Let us assume that $\beta$ is a constant Killing form i.e., $s_{0} = 0$ and $r^{}_{00} = 0$ ; putting this in Eq (3.3) vanishes the S-curvature.

For the if part, let us take $S = 0$ ; then Eq (3.3) becomes

$\begin{equation} t^{}_{0}+t^{}_{1}\alpha+t^{}_{2}\alpha^{2}+t^{}_{3}\alpha^{3}+t^{}_{4}\alpha^{4}+t^{}_{5}\alpha^{5} = 0, \end{equation}$

(3.4)

where

$\begin{equation*} \begin{split} & t^{}_{0} = 64 \beta^4 (4 \beta \Pi r^{}_{0} + 4 \beta \Pi s^{}_{0} - 3 n r^{}_{00}), \qquad t^{}_{1} = 4 \beta^3 (-16 \beta \Pi r^{}_{0} - 16 \beta \Pi s^{}_{0} + 3 (3 + 10 n)r^{}_{00}),\\& t^{}_{2} = (192 \beta^3 - 92 \beta^3 \Pi - 384 B \beta^3 \Pi) r^{}_{0} + (192 \beta^3 - 576 n \beta^3 - 92 384 B \beta^3 \Pi) s^{}_{0}\\& \; \; \; \; \; \; \; + (12 \beta^2 - 48 B \beta^2 + 18 n \beta^2 + 144 B n \beta^2) r^{}_{00},\\& t^{}_{3} = (-72 \beta^2 + 22 \beta^2 \Pi + 144 B \beta^2 \Pi) r^{}_{0} + (36\beta^2 + 72 n \beta^2 + 22 \beta^2 \Pi + 144 B \beta^2 \Pi) s^{}_{0} \\&\; \; \; \; \; \; \; + (-21 \beta - 24 B \beta- 21 n \beta - 108 B n \beta) r^{}_{00},\\& t^{}_{4} = (-36 \beta - 144 B \beta + 8 \beta \Pi+ 72 B \beta \Pi + 144 B^2 \beta \Pi) r^{}_{0}+ (54 \beta - 288 B \beta \\&\; \; \; \; \; \; + 90 n\beta + 432 B n \beta + 8 \beta \Pi + 72 B \beta \Pi+ 144 B^2 \beta \Pi) s^{}_{0}+ (3 + 24 B + 3 n + 18 B n) r^{}_{00},\\& t^{}_{5} = (12 + 72 B - 2 \Pi - 24 B \Pi - 72 B^2 \Pi) r^{}_{0} + (-6 - 72 B - 18 n - 108 B n - 2 \Pi - 24 B \Pi - 72 B^2 \Pi) s^{}_{0}. \end{split} \end{equation*}$

Taking the rational and irrational parts of Eq (3.4), we obtain

$\begin{equation} t^{}_{0}+\alpha^{2}(t^{}_{2}+\alpha^{2}t^{}_{4}) = 0, \end{equation}$

(3.5)

$\begin{equation} t^{}_{1}+ \alpha^{2}(t^{}_{3}+\alpha^{2}t^{}_{5}) = 0. \end{equation}$

(3.6)

From Eqs (3.5) and (3.6), we can say that $\alpha^2$ will divide $t_{0}$ as well as $t_{1}$ . In view of Lemma 1.1 $, \alpha^2$ is coprime with $\beta$ for $n > 2$ . Solving Eqs (3.5) and (3.6), we get, respectively,

$\begin{equation*} 4\beta\Pi(r_{0}+s_{0})-3nr_{00} = \gamma_{1}\alpha^2,\quad \text{for} \quad \gamma_{1} = \gamma_{1}(x), \end{equation*}$

and

$\begin{equation*} 16\beta \Pi(r_{0}+s_{0})-3(10n+3)r_{00} = \gamma_{2}\alpha^2,\quad \text{for} \quad \gamma_{2} = \gamma_{2}(x). \end{equation*}$

From the above equations, we obtain

$\begin{equation} r^{}_{00} = c\alpha^{2}, \quad \text{and then} \quad r^{}_{0} = c\beta, \end{equation}$

(3.7)

for some scalar function $c = c(x)$ on $M$ .

Putting the above values in Eq (3.4) and simplifying, we get

$\begin{equation*} 256 \Pi \beta^{5}(c\beta+s^{}_{0}) = \alpha^{2}\; (....), \end{equation*}$

where $(...)$ denotes the polynomial term in $\alpha$ and $\beta$ . Here also $\alpha^{2}$ does not divide $\beta^{5}$ and $(c\beta+s^{}_{0})$ . Therefore, $c\beta+s^{}_{0} = 0$ . Differentiating it with respect to $y^{i}$ , we obtain $c b^{}_{i}+s^{}_{i} = 0$ , which, on contracting by $b^{i}$ , gives $c = 0$ , implying $s^{}_{0} = 0$ and $r^{}_{00} = 0$ . Which means $\beta$ is a constant Killing form.

This completes, the proof of Theorem 1.1. □

4. Ricci curvature of cubic metric

In this section we obtain the projective Ricci curvature for the aforesaid metric.

Proof of Theorem 1.2. For this, we first obtain all the values of Eq (2.8) by using Eq (3.1) and the Mathematica program as

$\begin{equation*} \begin{split} &V^{}_{1} = \frac{-1}{(1 + n) (-1 - 6 B + s + 8 s^2)^4}\big[3 (6 B (6 (1 + n) + 8 (1 + n) s + (36 - (-37 + n) n) s^2 \\&\; \; \; \; \; \; \; \; - 4 (45 + n (37 + 8 n)) s^3 - 256 (4 + n (3 + n)) s^4) + s (-4 (1 + n) - 92 (1 + n) s + 92 (1 + n) s^2 \\&\; \; \; \; \; \; \; \; + (-238 + n (-241 + 3 n)) s^3 + 32 (23 + n (20 + 3 n)) s^4 + 256 (14 + n (11 + 3 n)) s^5) \\&\; \; \; \; \; \; \; \; + 3 B^2 (66 + 24 s (1 + 32 s) + (n + 16 n s)^2 + n (65 + 8 s (-1 + 64 s))))\big],\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{2} = \frac{-1}{(1 + n) (-1 + 2 s) (-1 - 6 B + s + 8 s^2)^4}\big[6 (-5 + 192 B + 1224 B^2 - 6 n + 186 B n + 1206 B^2 n \\&\; \; \; \; \; \; \; \; - n^2 - 18 B n^2 - 54 B^2 n^2 - 6 (3 (5 + 6 n + n^2) + 12 B^2 (-10 + n + 9 n^2) + B (16 + 41 n + 51 n^2)) s \\&\; \; \; \; \; \; \; \; + 3 (-99 - 104 n - n^2 - 36 B (-12 - 5 n + n^2) + 384 B^2 (8 + n + 3 n^2)) s^2 + 2 (508 + 675 n + 245 n^2 \\&\; \; \; \; \; \; \; \; + 12 B (-244 - 125 n + 105 n^2)) s^3 - 6 (330 + 183 n - 45 n^2 + 64 B (70 + 23 n + 15 n^2)) s^4 \\&\; \; \; \; \; \; \; \; - 96 (-36 - 11 n + 23 n^2) s^5 + 2048 (8 + 3 n + n^2) s^6)\big],\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{3} = \frac{3 (1 + 16 s) (6 B (3 + 5 n) + n (-2 + 2 s - 26 s^2) + 3 (-1 + s - 4 s^2) - n^2 (-1 + s + 2 s^2))}{(1 + n) (-1 - 6 B + s + 8 s^2)^3},\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{4} = \frac{6 (1 + 16 s) (-B + s^2)}{(-1 - 6 B + s + 8 s^2)^2},\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{5} = \frac{1}{(1 + n) (-1 + 2 s)^3 (-1 - 6 B + s + 8 s^2)^4}\big[36 (-6 - 7 n + n^2 - (63 + n (82 + 35 n)) s\!+\! 15 (8 \\&\; \; \; \; \; \; \; \; + 3 n (3 + n)) s^2 + (849 + n (1354 + 785 n)) s^3 - (2562 + n (2351 + 1123 n)) s^4 - 6 (-398 + 495 n \\&\; \; \; \; \; \; \; \; + 771 n^2) s^5 + 64 (45 + n (-7 + 100 n)) s^6 + 2048 (-3 + n (7 + 2 n)) s^7 - 216 B^3 (1 + n)^2 (-1 + 8 s) \\&\; \; \; \; \; \; \; \; + 9 B^2 (74 + 9 n (9 + n) - 144 s - 6 n (47 + 37 n) s + 48 (1 + n (-11 + 8 n)) s^2 \\&\; \; \; \; \; \; \; \; + 160 (2 + 3 n (5 + n)) s^3) + 6 B (13 + 13 n + 2 n^2 - (70 + n (109 + 89 n)) s + (151 + 55 n (1 + 2 n)) s^2 \\&\; \; \; \; \; \; \; \; + 2 (-178 + n (317+ 505 n)) s^3 - 8 (47 + n (-136 + 199 n)) s^4 - 1024 (-1 + n (5 + n)) s^5))\big], \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{6} = \frac{-72 n}{(1 + 6 B - s - 8 s^2)^2},\qquad V^{}_{7} = \frac{9 (-2 - 3 n + n^2)}{(1 + n) (1 + 6 B - s - 8 s^2)^2},\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{8} = \frac{1}{(1 \!+ \! n) (-1 \!+ \!2 s) (1 \!+ 6 B - s - 8 s^2)^3}\big[18 (-7 \!+ \! n (-8 \!+ \! 3 n) \! - \!33 s + \! 3 n (-4 \!+ \! 3 n) s \\&\; \; \; \; \; \; \; \; + \! 6 (10 \!+ \! (7 \! - \! 5 n) n) s^2 \!- \! 256 (1 \!+ \! 2 n) s^3 - 6 B (1 \!+ n \!- 2 n^2 \! + \! 4 (-14 \!+ \!(-23 + n) n) s))\big],\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{9} = \frac{6}{1 + 6 B - s - 8 s^2},\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{10} = \frac{-6 (1 + n + 6 B (3 + n)) + 18 (1 + 4 B (-15 + n) + n) s + 36 (3 + n) s^2 - 96 (-11 + n) s^3}{(1 + n) (-1 + 2 s) (-1 - 6 B + s + 8 s^2)^2},\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{11} = \frac{12 (1 + 3 B - 3 s - 60 B s - 3 s^2 + 64 s^3)}{(1 + n) (-1 + 2 s) (-1 - 6 B + s + 8 s^2)^2},\qquad V^{}_{12} = \frac{6 (1 - 8 s)}{(-1 + 2 s)^3},\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation} V^{}_{13} = \frac{-432n}{(1 + n) (-1 + 2 s) (-1 - 6 B + s + 8 s^2)^2},\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{equation}$

(4.1)

$\begin{equation*} \begin{split} &V^{}_{14} = \frac{18 (3 + 6 B - n - 6 B n + 3 (-3 + 4 B (-19 + n) + n) s+6 (-1 + n) s^2 - 16 (-15 + n) s^3)}{(1 + n) (-1 + 2 s)^2 (-1 - 6 B + s + 8 s^2)^2},\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{15} = \frac{18}{(-1 + 2 s) (-1 - 6 B + s + 8 s^2)},\qquad V^{}_{16} = \frac{6}{1-2s},\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation*} \begin{split} &V^{}_{17} = \frac{108}{(-1 + 2 s)^2 (-1 - 6 B + s + 8 s^2)},\quad V^{}_{18} = \frac{-9}{(-1 + 2 s)^2}, \quad V^{}_{19} = \frac{6 (-1 + n)}{(1 + n) (1 + 6 B - s - 8 s^2)}\,.\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

Plugging all the values of the above Eq (4.1) into Eq (2.7) and simplifying by the using Mathematica program, we obtain the projective Ricci curvature for the aforesaid metric as

$\begin{equation*} PRic = \frac{1}{(1 + n)^2 (\alpha - 2 \beta)^3 (-(1 + 6 B) \alpha^2 + \alpha \beta + 8 \beta^2)^4}\sum\limits_{i = 0}^{i = 13}\alpha^{i}t^{'}_{i}, \end{equation*}$

where

$\begin{equation*} \begin{split} &t^{'}_{0} = -2048 \beta^9 (-3 r^{2}_{00} (14 + n (11 + 3 n)) + 8 (1 + n) \beta (3 r^{}_{00|0} + 2 (1 + n) Ric^\alpha \beta) \\&\; \; \; \; \; \; + 8 \beta (2 (-1 + n) (1 + n)^2 \beta (\rho^{2}_{0} + \rho_{0|0})) - 3 (-1 + n^2) \rho_{0}r^{}_{00}),\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{split} \end{equation*}$

$\begin{equation} \begin{split} &t^{'}_{1} = 256 \beta^8 (-3 r^{2}_{00} (145 + n (112 + 33 n)) + 4 (1 + n) \beta (57 r^{}_{00|0} + 32 (1 + n) Ric^\alpha \beta) \\&\; \; \; \; \; \; \; + 4 \beta (32 (-1 + n) (1 + n)^2 \beta (\rho^{2}_{0} + \rho_{0|0})) - 57 (-1 + n^2) \rho_{0}r^{}_{00}),\\& \vdots \\& t^{'}_{13} = -9 (1 + 6 B)^3 (12 s_{k}s^{k} + s^{i}_{k}s^{k}_{i} + 6 B s^{i}_{k}s^{k}_{i}) (1 + n)^2. \end{split} \end{equation}$

(4.2)

Next, we obtain the flatness condition under which the projective Ricci curvature vanishes.

Let the projective Ricci curvature $PRic = 0$ , which implies $U(\alpha, \beta) = 0$ , where

$\begin{equation} U(\alpha, \beta) = t^{'}_{0}+\alpha t^{'}_{1}+\alpha^{2}t^{'}_{2}+......+\alpha^{13}t^{'}_{13}. \end{equation}$

(4.3)

Using Mathematica, we can see that

$\begin{equation*} U(\alpha, \beta) = \frac{1}{4}(-6-7n+n^2)(\alpha-2\beta)^3 (\alpha+\beta)^2 (\alpha+16\beta)^2(r^{}_{00}(\alpha-2\beta)-6s^{}_{0}\alpha^{2})^2\; \; \; \text{mod}[(1 + 6 B) \alpha^2-\alpha \beta-8 \beta^2]. \end{equation*}$

Therefore

$\begin{equation*} (...)[(1 + 6 B) \alpha^2-\alpha \beta-8 \beta^2])-\frac{1}{4}(-6-7n+n^2)(\alpha-2\beta)^3 (\alpha+\beta)^2 (\alpha+16\beta)^2(r^{}_{00}(\alpha-2\beta)-6s^{}_{0}\alpha^{2})^2 = 0, \end{equation*}$

where $(....)$ are polynomial in $\alpha$ and $\beta$ . As $B < \frac{1}{4}$ , therefore $((1 + 6 B) \alpha^2-\alpha \beta-8 \beta^2)$ does not divide $(\alpha-2\beta)^3$ or $(\alpha+\beta)^2$ or $(\alpha+16\beta)^2$ . Therefore $((1 + 6 B) \alpha^2-\alpha \beta-8 \beta^2)$ will divide $(r^{}_{00}(\alpha-2\beta)-6s^{}_{0}\alpha^{2})^2$ ; then $((1 + 6 B) \alpha^2-\alpha \beta-8 \beta^2)$ will also divide $(r^{}_{00}(\alpha-2\beta)-6s^{}_{0}\alpha^{2})$ , i.e.,

$\begin{equation*} (r^{}_{00}(\alpha-2\beta)-6s^{}_{0}\alpha^{2}) = (c^{}_{1}+\alpha c^{}_{0})((1 + 6 B) \alpha^2-\alpha \beta-8 \beta^2), \end{equation*}$

where $c^{}_{1}$ is a $1$ -form and $c^{}_{0}$ is a scalar. Taking the rational and irrational parts of the above equation, we obtain

$\begin{equation} -2\beta r^{}_{00}- 6\alpha^{2}s^{}_{0} = c^{}_{1}\alpha^{2}(1+6B)-8\beta^{2}c^{}_{1}-c^{}_{0}\alpha^{2}\beta, \end{equation}$

(4.4)

and

$\begin{equation} r^{}_{00} = c^{}_{0}\alpha^{2}(1+6B)-\beta c^{}_{1}-8c^{}_{0}\beta^{2}. \end{equation}$

(4.5)

Solving the above equations, we get $c^{}_{1} = -\frac{8}{5}\beta c^{}_{0}$ , and then (4.5) gives

$\begin{equation} r^{}_{00} = c^{}_{0}[\alpha^{2}(1+6B)-\frac{32}{5}\beta^{2}]. \end{equation}$

(4.6)

Substituting the above values into Eq (4.4), we obtain

$\begin{equation} (4B-1)c^{}_{0}\beta+10 s^{}_{0} = 0. \end{equation}$

(4.7)

Differentiating the above equation with respect to $y^{i}$ gives $(4B-1)c^{}_{0}b^{}_{i}+s^{}_{i} = 0$ , which, on contracting by $b^{i}$ , we obtain $c^{}_{0} = 0$ . Then from Eqs (4.6) and (4.7), we obtain

$\begin{equation} r^{}_{00} = 0,\qquad s^{}_{0} = 0. \end{equation}$

(4.8)

In view of (4.8), Eq (4.3) becomes

$\begin{equation*} \begin{split} & 3 \alpha^2 (-2 s_{0k}s^{k}_{0} (\alpha - 8 \beta) + (-3 s^{i}_{k}s^{k}_{i} \alpha^2 + 2 s^{k}_{0;k} (\alpha - 2 \beta)) (\alpha - 2 \beta)) + Ric^\alpha (\alpha - 2 \beta)^3 \\&- (n-1) (\alpha - 2 \beta)^2 (-(\alpha - 2 \beta) \rho^{2}_{0} + 6 \alpha^2 s^{k}_{0}\rho_{k}- (\alpha - 2 \beta) \rho_{0|0}) = 0, \end{split} \end{equation*}$

which can be rewritten as

$\begin{equation*} \begin{split} &(\alpha - 2 \beta)\Bigg\{6 \alpha^2s_{0k}s^{k}_{0}+9\alpha^4s^{i}_{k}s^{k}_{i}-6\alpha^2s^{k}_{0;k} (\alpha - 2 \beta) -Ric^\alpha (\alpha - 2 \beta)^2\\&+ (n-1)(\alpha - 2 \beta)\big[6 \alpha^2 s^{k}_{0}\rho_{k}-(\alpha - 2 \beta) \rho^{2}_{0}- (\alpha - 2 \beta) \rho_{0|0}\big]\Bigg\} = 36s^{}_{0k}s^{k}_{0}\alpha^2\beta. \end{split} \end{equation*}$

Since $(\alpha-2\beta)$ does not divide $\alpha^2$ or $\beta$ , therefore $(\alpha-2\beta)$ will divide $s^{}_{0k}s^{k}_{0}$ . Thus

$\begin{equation*} s^{}_{0k}s^{k}_{0} = (d^{}_{1}+\alpha d^{}_{0})(\alpha-2\beta), \end{equation*}$

where $d^{}_{1}$ is a $1$ -form and $d^{}_{0}$ is a scalar. Taking the rational and irrational parts of the above equation and solving, we obtain

$\begin{equation} s^{}_{0k}s^{k}_{0} = d^{}_{0}(\alpha^2-4\beta^2). \end{equation}$

(4.9)

If $d^{}_{0} \neq 0$ then one can conclude by the above equation that $\alpha$ is not positive definite, which is not possible. Therefore, $d^{}_{0} = 0$ . This implies that

$\begin{equation} s^{}_{ik} = 0, \end{equation}$

(4.10)

i.e., $\beta$ is closed. In view of Eqs (4.8) and (4.10), we obtain $b^{}_{i; k} = 0$ , then $1$ -form $\beta$ is parallel with respect to $\alpha$ .

This completes the proof of Theorem 1.2. □

Now, we obtain the condition for the weak projective Ricci curvature of a cubic Finsler metric.

Proof of Theorem 1.3. Let $F$ be a cubic Finsler metric with weak projective Ricci curvature. Then from Eq (2.6) we obtain

$\begin{equation} (n\!-\!1)\big[3\theta(\alpha\!+\!\beta)^3 \alpha^2+\!\gamma(\alpha\!+\!\beta)^6\big] = \frac{\alpha^4}{(1 \!+ n\!)^2 (\alpha \!- \! 2 \beta)^3 ((1\! +\! 6 B) \alpha^2 \!-\! \alpha \beta \!- \! 8 \beta^2)^4}\sum\limits_{i = 0}^{i = 13}\alpha^{i}t^{'}_{i}. \end{equation}$

(4.11)

For the cubic metric, we have $B < \frac{1}{4}$ , which implies that $\alpha^4$ does not divide $(\alpha\!-\!2\beta)^3$ or $((1\! +\! 6 B) \alpha^2 \!-\! \alpha \beta$ $- 8 \beta^2)^4$ or $3\theta(\alpha\!+\!\beta)^3 \alpha^2$ . Consequently, it follows that $\alpha^2$ must divide $\gamma(\alpha\!+\!\beta)^6$ . However, such division is only possible if $\gamma = 0$ . Combining this result with Eq (4.11), then we deduce that $3\theta(\alpha\!+\!\beta)^3$ is divided by $\alpha^2$ . This is impossible unless $\theta = 0$ . Then $F$ reduces to a projective Ricci-flat metric.

The converse is obvious. This completes the proof. □

Example 4.1. The Finsler metric $\frac{ 1}{|y|^2}\big(|y|+ < a, y > \big)^3$ for $a =$ constant is projectively Ricci flat.

5. Conclusions

Projective Ricci curvature is a concept in differential geometry that generalizes the notion of Ricci curvature. It has various applications in the fields of general relativity, optimal transformation theory, complex geometry, Weyl geometry, Einstein metrics, and many more. In this article, we have proved that if the cubic metric $F = \frac{(\alpha+\beta)^3}{\alpha^2}$ is projective Ricci flat $(PRic = 0),$ then $\beta$ is parallel with respect to Riemannian metric $\alpha$ , and then from Eq (2.3), the S-curvature vanishes. Therefore, from Eq (2.5), we obtain that the Riemannian metric $\alpha$ is also Ricci-flat, which is Corollary 1.1.

Author contributions

M. K. Gupta and S. Sharma wrote the framework and the original draft of this manuscript. Y. Li and Y. Xie reviewed and validated the manuscript. All authors have read and agreed to the final version of the manuscript.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Conflict of interest

The authors declare no conflict of interest.

References

[1]	X. Lin, Z. Su, Y. Yang, S. Zhang, The potential of ionic liquids in biopharmaceutical engineering, Chinese J. Chem. Eng., 30 (2021), 236–243. https://doi.org/10.1016/j.cjche.2020.11.015 doi: 10.1016/j.cjche.2020.11.015
[2]	J. Yang, W. Li, Q. Liu, H. Liu, Dissolution of antibiotics mycelium in ionic liquids: Performance and mechanism, Chinese J. Chem. Eng., 26 (2017), 252–258. https://doi.org/10.1016/j.cjche.2017.04.003 doi: 10.1016/j.cjche.2017.04.003
[3]	S. Basharat, Z. Huang, M. Gong, X. Lv, A. Ahmed, I. Hussain, et al., A review on current conventional and biotechnical approaches to enhance biosynthesis of steviol glycosides in Stevia rebaudiana, Chinese J. Chem. Eng., 30 (2021), 92–104. https://doi.org/10.1016/j.cjche.2020.10.018 doi: 10.1016/j.cjche.2020.10.018
[4]	Y. Yu, G. Li, W. Han, L. Zhu, T. Si, H. Wang, et al., An efficient preparation of porous polymeric microspheres by solvent evaporation in foam phase, Chinese J. Chem. Eng., 29 (2021), 409–416. https://doi.org/10.1016/j.cjche.2020.09.002 doi: 10.1016/j.cjche.2020.09.002
[5]	M. N. Pantano, M. C. Fernandez, M. E. Serrano, O. A. Ortiz, G. J. E. Scaglia, Tracking control of optimal profiles in a nonlinear fed-batch bioprocess under parametric uncertainty and process disturbances, Ind. Eng. Chem. Res., 57 (2018), 11130–11140. https://doi.org/10.1021/acs.iecr.8b01791 doi: 10.1021/acs.iecr.8b01791
[6]	D. Niu, L. Zhang, F. Wang, Modeling and parameter updating for nosiheptide fed-batch fermentation process, Ind. Eng. Chem. Res., 55 (2016), 8395–8402. https://doi.org/10.1021/acs.iecr.6b01245 doi: 10.1021/acs.iecr.6b01245
[7]	F. A. Gorrini, S. I. Biagiola, J. L. Figueroa, A. V. Wouwer, Simple observer-based feedback strategy for controlling fed-batch hybridoma cultures, Ind. Eng. Chem. Res., 56 (2017), 15072–15082. https://doi.org/10.1021/acs.iecr.7b03395 doi: 10.1021/acs.iecr.7b03395
[8]	A. Simorgh, A. Razminia, J. A. T. Machado, Optimal control of nonlinear fed-batch process using direct transcription method, Comput. Chem. Eng., 130 (2019), 106561. https://doi.org/10.1016/j.compchemeng.2019.106561 doi: 10.1016/j.compchemeng.2019.106561
[9]	L. Veiter, J. Kager, C. Herwig, Optimal process design space to ensure maximum viability and productivity in Penicillium chrysogenum pellets during fed-batch cultivations through morphological and physiological control, Microb. Cell Fact., 19 (2020), 33. https://doi.org/10.1186/s12934-020-1288-5 doi: 10.1186/s12934-020-1288-5
[10]	J. Schorsch, C. C. Castro, L. D. Couto, C. Nobre, M. Kinnaert, Optimal control for fermentative production of fructo-oligosaccharides in fed-batch bioreactor, J. Process Contr., 78 (2019), 124–138. https://doi.org/10.1016/j.jprocont.2019.03.004 doi: 10.1016/j.jprocont.2019.03.004
[11]	X. Wu, K. Zhang, M. Cheng, Optimal control of bioprocess systems using hybrid numerical optimization algorithms, Optimization, 67 (2018), 1287–1306. https://doi.org/10.1080/02331934.2018.1466299 doi: 10.1080/02331934.2018.1466299
[12]	F. Rossi, S. Copelli, A. Colombo, C. Pirola, F. Manenti, Online model-based optimization and control for the combined optimal operation and runaway prediction and prevention in (fed-)batch systems, Chem. Eng. Sci., 138 (2015), 760–771. https://doi.org/10.1016/j.ces.2015.09.006 doi: 10.1016/j.ces.2015.09.006
[13]	H. S. Shin, H. C. Lim, Optimal fed-batch operation for recombinant cells with segregational plasmid instability, Chem. Eng. Commun., 195 (2008), 1122–1143. https://doi.org/10.1080/00986440801943396 doi: 10.1080/00986440801943396
[14]	R. F. Rao, X. D. Li, Input-to-state stability in the meaning of switching for delayed feedback switched stochastic financial system, AIMS Mathematics, 6 (2021), 1040–1064. https://doi.org/10.3934/math.2021062 doi: 10.3934/math.2021062
[15]	X. Wu, K. Zhang, M. Cheng, Adaptive numerical approach for optimal control of a single train, J. Syst. Sci. Complex, 32 (2019), 1053–1071. https://doi.org/10.1007/s11424-018-7277-7 doi: 10.1007/s11424-018-7277-7
[16]	X. Wu, K. Zhang, M. Cheng, Optimal control of constrained switched systems and application to electrical vehicle energy management, Nonlinear Anal. Hybrid Syst., 30 (2018), 171–188. https://doi.org/10.1016/j.nahs.2018.05.006 doi: 10.1016/j.nahs.2018.05.006
[17]	X. Wu, K. Zhang, M. Cheng, A switched dynamical system approach towards the optimal control of chemical processes based on a gradient-based parallel optimization algorithm, Comput. Chem. Eng., 118 (2018), 180–194. https://doi.org/10.1016/j.compchemeng.2018.08.007 doi: 10.1016/j.compchemeng.2018.08.007
[18]	D. Sahoo, G. Samanta, Impact of fear effect in a two prey-one predator system with switching behaviour in predation, Differ. Equ. Dyn. Syst., in press. https://doi.org/10.1007/s12591-021-00575-7
[19]	S. Saha, G. Samanta, Modelling of a two prey and one predator system with switching effect, Comput. Math. Biophys., 9 (2021), 90–113. https://doi.org/10.1515/cmb-2020-0120 doi: 10.1515/cmb-2020-0120
[20]	S. Saha, G. Samanta, Modeling of insect-pathogen dynamics with biological control, Math. Biol. Bioinform., 15 (2020), 268–294. https://doi.org/10.17537/2020.15.268 doi: 10.17537/2020.15.268
[21]	M. Huang, S. Liu, Y. Zhang, Mathematical modeling and analysis of biological control strategy of aphid population, AIMS Mathematics, 7 (2022), 6876–6897. https://doi.org/10.3934/math.2022382 doi: 10.3934/math.2022382
[22]	D. Liberzon, A. S. Morse, Basic problems in stability and design of switched systems, IEEE Control Syst. Mag., 19 (1999), 59–70. https://doi.org/10.1109/37.793443 doi: 10.1109/37.793443
[23]	H. Lin, P. J. Antsaklis, Stability and stabilizability of switched linear systems: A survey of recent results, IEEE Trans. Auto. Control, 54 (2009), 308–322. https://doi.org/10.1109/.2005.1466986 doi: 10.1109/.2005.1466986
[24]	D. Yang, C. X. Huang, G. D. Zong, Finite-time H-infinity bumpless transfer control for switched systems: A state-dependent switching approach, Int. J. Robust Nonlinear Control, 30 (2020), 1417–1430. https://doi.org/10.1002/rnc.4825 doi: 10.1002/rnc.4825
[25]	X. Huo, L. Ma, X. D. Zhao, G. D. Zong, Event-triggered adaptive fuzzy output feedback control of MIMO switched nonlinear systems with average dwell time, Appl. Math. Comput., 365 (2020), 124665. https://doi.org/10.1016/j.amc.2019.124665 doi: 10.1016/j.amc.2019.124665
[26]	A. Kundu, A new condition for stability of switched linear systems under restricted minimum dwell time switching, Syst. Control Lett., 135 (2020), 104597. https://doi.org/10.1016/j.sysconle.2019.104597 doi: 10.1016/j.sysconle.2019.104597
[27]	A. A. Kahloul, A. Sakly, Hybrid approach for constrained optimal control of nonlinear switched systems, J. Control Auto. Electr. Syst., 31 (2020), 865–873. https://doi.org/10.1007/s40313-020-00586-9 doi: 10.1007/s40313-020-00586-9
[28]	X. Xu, X. Mao, Y. Li, H. Zhang, New result on robust stability of switched systems with all subsystems unstable, IET Control Theory Appl., 13 (2019), 2138–2145. https://doi.org/10.1049/iet-cta.2019.0018 doi: 10.1049/iet-cta.2019.0018
[29]	S. Sui, S. Tong, C. L. Chen, K. Sun, Fuzzy adaptive optimal control for nonlinear switched systems with actuator hysteresis, Int. J. Adapt. Control Signal Process., 33 (2019), 609–625. https://doi.org/10.1002/acs.2975 doi: 10.1002/acs.2975
[30]	P. Wang, J. Zhao, Stability and guaranteed cost analysis of switched positive systems with mode-dependent dwell time and sampling, IET Control Theory Appl., 14 (2020), 378–385. https://doi.org/10.1049/iet-cta.2019.0466 doi: 10.1049/iet-cta.2019.0466
[31]	X. Xu, P. J. Antsaklis, Optimal control of switched systems based on parameterization of the switching instants, IEEE Trans. Auto. Control, 49 (2004), 2–16. https://doi.org/10.1109/TAC.2003.821417 doi: 10.1109/TAC.2003.821417
[32]	X. Wu, K. Zhang, M. Cheng, Computational method for optimal control of switched systems with input and state constraints, Nonlinear Anal. Hybrid Syst., 26 (2017), 1–18. https://doi.org/10.1016/j.nahs.2017.04.001 doi: 10.1016/j.nahs.2017.04.001
[33]	H. Chen, W. Zhang, On weak topology for optimal control of switched nonlinear systems, Automatica, 81 (2017), 409–415. https://doi.org/10.1016/j.automatica.2017.03.039 doi: 10.1016/j.automatica.2017.03.039
[34]	X. Liu, S. Li, K. Zhang, Optimal control of switching time in switched stochastic systems with multi-switching times and different costs, Int. J. Control, 90 (2017), 1604–1611. https://doi.org/10.1080/00207179.2016.1214879 doi: 10.1080/00207179.2016.1214879
[35]	X. Wu, K. Zhang, C. Sun, Numerical algorithm for a class of constrained optimal control problems of switched systems, Numer. Algor., 67 (2014), 771–792. https://doi.org/10.1007/s11075-013-9822-8 doi: 10.1007/s11075-013-9822-8
[36]	H. G. Bock, C. Kirches, A. Meyer, A. Potschka, Numerical solution of optimal control problems with explicit and implicit switches, Optim. Meth. Softw., 33 (2018), 450–474. https://doi.org/10.1080/10556788.2018.1449843 doi: 10.1080/10556788.2018.1449843
[37]	X. Xu, P. J. Antsaklis, Optimal control of switched systems via non-linear optimization based on direct differentiations of value functions, Int. J. Control, 75 (2002), 1406–1426. https://doi.org/10.1080/0020717021000023825 doi: 10.1080/0020717021000023825
[38]	X. Wu, K. Zhang, C. Sun, Constrained optimal control of switched systems based on modified BFGS algorithm and filled function method, Int. J. Comput. Math., 91 (2014), 1713–1729. https://doi.org/10.1080/00207160.2013.859678 doi: 10.1080/00207160.2013.859678
[39]	X. Wu, K. Zhang, Three-dimensional trajectory design for horizontal well based on optimal switching algorithms, ISA Trans., 58 (2015), 348–356. https://doi.org/10.1016/j.isatra.2015.04.002 doi: 10.1016/j.isatra.2015.04.002
[40]	R. Vasudevan, H. Gonzalez, R. Bajcsy, S. S. Sastry, Consistent approximations for the optimal control of constrained switched systems-Part 1: A conceptual algorithm, SIAM J. Control Optim., 51 (2013), 4463–4483. https://doi.org/10.1137/120901490 doi: 10.1137/120901490
[41]	R. Vasudevan, H. Gonzalez, R. Bajcsy, S. S. Sastry, Consistent approximations for the optimal control of constrained switched systems-Part 2: An implementable algorithm, SIAM J. Control Optim., 51 (2013), 4484–4503. https://doi.org/10.1137/120901507 doi: 10.1137/120901507
[42]	X. Wu, K. Zhang, M. Cheng, X. Xin, A switched dynamical system approach towards the economic dispatch of renewable hybrid power systems, Int. J. Electr. Power Energy Syst., 103 (2018), 440–457. https://doi.org/10.1016/j.ijepes.2018.06.016 doi: 10.1016/j.ijepes.2018.06.016
[43]	T. Sardarmehni, A. Heydari, Sub-optimal scheduling in switched systems with continuous-time dynamics: A gradient descent approach, Neurocomputing, 285 (2018), 10–22. https://doi.org/10.1016/j.neucom.2018.01.003 doi: 10.1016/j.neucom.2018.01.003
[44]	H. R. Tabrizidooz, M. Pourbabaee, M. Hedayati, Optimal control of switched systems by a modified pseudo spectral method, Iran J. Math. Chem., 8 (2017), 161–173. https://doi.org/10.22052/ijmc.2017.44718 doi: 10.22052/ijmc.2017.44718
[45]	G. Wu, J. Sun, J. Chen, Optimal linear quadratic regulator of switched systems, IEEE Trans. Auto. Control, 64 (2019), 2898–2904. https://doi.org/10.1109/TAC.2018.2872204 doi: 10.1109/TAC.2018.2872204
[46]	A. Heydari, Optimal switching with minimum dwell time constraint, J. Franklin Inst., 354 (2017), 4498–4518. https://doi.org/10.1016/j.jfranklin.2017.04.015 doi: 10.1016/j.jfranklin.2017.04.015
[47]	T. M. Caldwell, T. D. Murphey, Switching mode generation and optimal estimation with application to skid-steering, Automatica, 47 (2011), 50–64. https://doi.org/10.1016/j.automatica.2010.10.010 doi: 10.1016/j.automatica.2010.10.010
[48]	X. Wu, K. Zhang, M. Cheng, Sensitivity analysis for an optimal control problem of chemical processes based on a smoothing cost penalty function approach, Chem. Eng. Res. Design., 146 (2019), 221–238. https://doi.org/10.1016/j.cherd.2019.04.011 doi: 10.1016/j.cherd.2019.04.011
[49]	X. Wu, Q. Liu, K. Zhang, X. Xin, Optimal-tuning of proportional-integral-derivative-like controller for constrained nonlinear systems and application to ship steering control, J. Franklin Inst., 355 (2018), 5667–5689. https://doi.org/10.1016/j.jfranklin.2018.06.017 doi: 10.1016/j.jfranklin.2018.06.017
[50]	K. Dehingia, H. K. Sarmah, K. Hosseini, K. Sadri, S. Salahshour, C. Park, An optimal control problem of immuno-chemotherapy in presence of gene therapy, AIMS Mathematics, 6 (2021), 11530–11549. https://doi.org/10.3934/math.2021669 doi: 10.3934/math.2021669
[51]	A. Das, K. Dehingia, H. K. Sharmah, C. Park, J. R. Lee, K. Sadri, et al., Optimal control of effector-tumor-normal cells dynamics in presence of adoptive immunotherapy, AIMS Mathematics, 6 (2021), 9813–9834. https://doi.org/10.3934/math.2021570 doi: 10.3934/math.2021570
[52]	A. Das, H. K. Sarmah, D. Bhattacharya, K. Dehingia, K. Hosseini, Combination of virotherapy and chemotherapy with optimal control for combating cancer, Math. Comput. Simulat., 194 (2022), 460–488. https://doi.org/10.1016/j.matcom.2021.12.004 doi: 10.1016/j.matcom.2021.12.004
[53]	B. J. Nath, H. K. Sarmah, H. Maurer, An optimal control strategy for antiretroviral treatment of HIV infection in presence of immunotherapy, Qual. Theory Dyn. Syst., 21 (2022), 30. https://doi.org/10.1007/s12346-022-00564-3 doi: 10.1007/s12346-022-00564-3
[54]	N. U. Ahmed, Dynamic systems and control with applications, Singapore: World Scientific, 2006. https://doi.org/10.1142/6262
[55]	H. M. Saber, A. Ravindran, Nonlinear goal programming theory and practice: A survey, Comput. Oper. Res., 20 (1993), 275–291. https://doi.org/10.1016/0305-0548(93)90004-3 doi: 10.1016/0305-0548(93)90004-3
[56]	K. L. Teo, C. Goh, K. Wong, A unified computational approach to optimal control problems, Essex: Longman Scientific $&$ Technical, 1991. https://doi.org/10.1515/9783110883237.2763
[57]	S. S. Rao, Engineering optimization: theory and practice, Chichester: John Wiley & Sons, 2019. https://doi.org/10.1002/9781119454816
[58]	Y. Mu, H. Teng, D. J. Zhang, W. Wang, Z. L. Xiu, Microbial production of 1, 3-propanediol by Klebsiella pneumoniae using crude glycerol from biodiesel preparations, Biotechnol Lett., 28 (2006), 1755–1759. https://doi.org/10.1007/s10529-006-9154-z doi: 10.1007/s10529-006-9154-z
[59]	A. Nikoobin, M. Moradi, Indirect solution of optimal control problems with state variable inequality constraints: finite difference approximation, Robotica, 35 (2017), 50–72. https://doi.org/10.1017/S0263574715000521 doi: 10.1017/S0263574715000521

This article has been cited by:

Md Aquib, Slant Submersions in Generalized Sasakian Space Forms and Some Optimal Inequalities, 2025, 14, 2075-1680, 417, 10.3390/axioms14060417

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)