\documentclass[twoside,12pt]{article}
\usepackage{amsmath,amssymb, amsthm}
\usepackage{mathtools}
%\usepackage{calligra}
%\usepackage[T1]{fontenc}%\calligra
\usepackage[mathscr]{euscript}
\usepackage{mathrsfs}
%\usepackage{calrsfs}
\usepackage{dsfont}
\usepackage[colorlinks,urlcolor=red]{hyperref}
\usepackage{graphicx}
\usepackage{orcidlink}
\renewcommand{\baselinestretch}{1.2}
\setlength{\oddsidemargin}{-0.5in}
\setlength{\evensidemargin}{-0.5in} \setlength{\textwidth}{6.65in}
\setlength{\topmargin}{0.15in} \setlength{\parindent}{0.3in}
\setlength{\parskip}{0.3truecm} \setlength{\hoffset}{0.40in}
\newenvironment{narrownum}%


% Set the height of the header
\setlength{\headheight}{0.25in}


% Set vertical distance between the header and the text
\setlength{\headsep}{0.25in}


% Set the height of the text
\setlength{\textheight}{8.8in}

% Set vertical distance between the text and the bottom of the footer
\setlength{\footskip}{0.25in}

% THEOREMS -------------------------------------------------------
\newtheorem{thm}{Theorem}[section]
\newtheorem{cor}[thm]{Corollary}
\newtheorem{lem}[thm]{Lemma}
\newtheorem{prop}[thm]{Proposition}
\theoremstyle{definition}
\newtheorem{definition}{Definition}[section]
\newtheorem{theorem}[definition]{Theorem}
\newtheorem{algorithm}[definition]{Algorithm}
\newtheorem{example}[definition]{Example}
\newtheorem{problem}[definition]{Problem}
\newtheorem{proposition}[definition]{Proposition}
\newtheorem{corollary}[definition]{Corollary}
\newtheorem{remark}[definition]{Remark}
\newtheorem{lemma}[definition]{Lemma}
\newtheorem{defn}[thm]{Definition}
\newtheorem{algo}[thm]{Algorithm}
\newtheorem{Ex}[thm]{Example}
\theoremstyle{remark}
\newtheorem{rem}[thm]{Remark}
\numberwithin{equation}{section}
% MATH -----------------------------------------------------------
\newcommand{\Z}{\Bbb Z}
\newcommand{\lm}{$\lim_{n \to \infty}$}
\newcommand{\snx}{\frac{1}{n}\sum_{j=1}^{n}(PT)^jx_n}
\newcommand{\sn}{\frac{1}{n}\sum_{j=1}^{n}}
\newcommand{\snn}{\frac{1}{n+1}\sum_{j=0}^{n}}
\newcommand{\lsn}{$\limsup_{n\to\infty}$}
\newcommand{\epg}{\left\{%
\begin{array}{ll}
    F(u_n,y)=\frac{1}{r_n}\langle y-u_n,u_n-x_n\rangle, & \hbox{$\forall y\in C$;} \\
    x_{n+1}=\alpha_n\gamma f(x_n) + (I-\alpha_nA)Tu_n, & . \\
\end{array}%
\right.}
\newcommand{\tn}{\alpha_{n-1}\gamma f(x_{n-1}) + (I-\alpha_{n-1}A)Su_{n-1}}
\newcommand{\tnn}{\alpha_n\gamma f(x_n) + (I-\alpha_nA)Su_n}
\newcommand{\inx}{\frac{1}{t_n}\int_0^{t_n}T(s)x_nds}
\newcommand{\iny}{\frac{1}{t_{n-1}}\int_0^{t_{n-1}}T(s)x_{n-1}ds}
\newcommand{\inz}{\frac{1}{t_n}\int_0^{t_n}T(s)x_{n-1}ds}
\newcommand{\ep}{\epsilon}
\renewcommand{\to}{\longrightarrow}
\newcommand{\dprime}{^{\prime\prime}}
\DeclareMathAlphabet{\pazocal}{OMS}{zplm}{m}{n}
\newenvironment{pf}
    {{\noindent  \textrm{\textbf{Proof }}}~~}


\begin{document}
\pagestyle{myheadings} \markboth{\rm\centerline{%Shih-sen Chang and
H. Gissy, A. A. H. Ahmadini, Salahuddin
}}{\rm\centerline{Solutions for the system of nonlinear mixed variational inequality problems}}

\begin{titlepage}
\title{\bf Solutions for the system of nonlinear mixed variational inequality problems}
\author{%Shih-sen Chang$^\dag$ and
H. Gissy, A. A. H. Ahmadini, Salahuddin\orcidlink{0000-0002-0496-3379}\thanks{Corresponding author.\textit{Email address}:
%changss2013@163.com(Shih-sen Chang),
dr$\_$salah12@yahoo.com(Salahuddin); smohammad@jazanu.edu.sa (Salahuddin), hgissy@jazanu.edu.sa (Gissy), aahmadini@jazanu.edu.sa(Ahmadini)
}
 \\
%\textit{\small{$^\dag$Center for General Education, China Medical University}}\\
%\textit{\small{Taichung-40402, Taiwan, China}}\\
\textit{\small{Department of Mathematics, College of Science, Jazan University, P.O. Box 114}}\\
\textit{\small{Jazan-45142, Kingdom of Saudi Arabia.}}}
\end{titlepage}
\date{}
\maketitle
%==================================================================================================================================
\vspace{-1cm} \noindent\hrulefill \vspace{-.2cm}
\begin{abstract}
In this paper, we propose a system of nonlinear mixed variational inequality problems, which consists of two elliptic mixed variational inequality problems on Banach spaces. Under suitable assumptions, using the Kakutani-Ky Fan fixed point theorem and Minty techniques, we prove the solution set to the system of nonlinear mixed variational inequality problem is nonempty, weakly compact and unique. Additionally, we suggest a stability result for the system of nonlinear mixed variational inequality problem by perturbing the duality mappings. Furthermore, we present an optimal control problem governed by the system of nonlinear mixed variational inequality problems and establish a solvability result.\par
\vspace{.2cm}
\noindent
\emph{Keywords}: System of nonlinear mixed variational inequality problem, Inverse relaxed monotonicity, Existence, Uniqueness, Stability, Optimal control problem.\par
\noindent
{\bf AMS Mathematics Subject Classification:} 47J20, 49J30, 35A15, 49N45, 47H20, 49H52.
\end{abstract}
\vskip-.6cm \hrulefill
%=================================================================================================
\section{Introduction}
Lin \cite{Li} introduced the system of generalized quasi-variational inclusion problems. This system consists of a set of problems defined on a product set. It includes several well-known problems such as variational inequalities, equilibrium problems, vector equilibrium problems, vector quasi-equilibrium problems, variational inclusions problems and variational disclusion problems.\\
Undoubtedly, in the realm of engineering, sciences, technology, chemical processes, and economics, several challenging and complex problems frequently result in inequalities instead of straightforward equations. In this scenario, variational inequalities have emerged as a powerful mathematical resource. Variational inequalities essentially arise from applied models with an underlying convex foundation and have been the subject of extensive research since the 1960s, encompassing mathematical theories, numerical techniques, and practical applications, among other significant sources, \textit{see,} \cite{ChSaAhWaWa, ChSaWaMa, HaSo, HlHaNeLo, KiSaAh, KiSt}. \\
It should be noted that the results mentioned earlier cannot be applied to coupled systems that consist of two elliptical mixed variational inequalities. However, coupled variational inequalities can be a useful mathematical tool for examining various coupled mixed boundary value problems, control problems, and similar problems. More details can be found in \cite{AlAhSa, DeMiPa03, GaPa, LiXuLi, Ze}.\par
In this paper, we introduce a system of nonlinear mixed variational inequality problems and using the Kakutani-Ky Fan fixed point theorem, the Minty techniques and inverse relaxed monotonicity  to establish the existence, convergence, uniqueness, stability, and optimal control of the problems.\\
Before we proceed, let us first define the problem that will play a crucial role in this paper. Consider two reﬂexive Banach spaces $(X, \| \cdot \|_X)$ and $(Y, \|\cdot\|_Y)$, each with its dual spaces $(X^\ast,\|\cdot\|_{X^\ast} )$ and $(Y^\ast,\|\cdot\|_{ Y^\ast} )$, \textit{respectively}. We use $\langle\cdot,\cdot\rangle_X$ (resp., $\langle\cdot,\cdot\rangle_Y$ ) to denote the duality pairing between $X^\ast$  and $X$ (resp., the duality pairing between $Y^\ast$  and $Y$). We use the symbols $\xlongrightarrow[]{w}$ and $\longrightarrow$ to denote the weak and the strong convergence in the space $X$ and employ the notation $X_w$ for the space $X$ equipped with the weak topology. The limits, lower limits and upper limits are considered as $n\to \infty,$ even if we do not mention it explicitly. 
%Recall that a function $\gamma\colon X\to\bar{\mathbb{R}}=\mathbb{R}\cup\{+\infty\}$ is called proper, convex, and lower semicontinuous if it fulﬁlls the conditions
%\begin{align*}
%D(\gamma)=&\{u\in X \mid \gamma(u) < +\infty\}\neq\emptyset,\\
%\gamma(\lambda u+(1-\lambda)v)&\leq\lambda \gamma(u)+(1-\lambda)\gamma(v),~\forall \lambda\in [0,1]~\text{and}~u,v\in X,\\
%\gamma(u)\leq&\liminf_{n\to\infty}\gamma(u_n),~\text{for all sequences} ~~\{u_n\}\subset X~\text{with}~~ u_n\to u,
%\end{align*}
%\textit{respectively}.\par
We can formulate a system of nonlinear mixed variational inequality problems on Banach spaces:
\begin{problem} \label{P1.1}
Find $(x, y)\in \Omega\times \mho$ such that
\begin{equation}\label{eq:1.1}
\langle \mathcal{A}(y, x), v-x\rangle_ X+\varphi(x,v)-\varphi(x,x)\geq \langle \gamma, v- x\rangle_ X,~ \forall v\in\Omega 
\end{equation} 
and
\begin{equation}\label{eq:1.2}
\langle \mathcal{B}(x, y),w-y\rangle_Y+\phi(y,w)-\phi(y,y)\geq\langle \zeta,w-y\rangle_Y, ~\forall w\in \mho.
\end{equation}
\end{problem}
We note that if $\varphi(x,x)=\varphi(x) $ and $\phi(x,x)=\phi(x)$, then Problem \ref{P1.1} reduces to the coupled system of variational inequality problems for finding $(x, y)\in \Omega\times \mho$ such that
\begin{equation}\label{eq:1.3}
\begin{cases}
\langle \mathcal{A}(y, x), v-x\rangle_ X+\varphi(v)-\varphi(x)\geq \langle \gamma, v- x\rangle_ X,~ \forall v\in\Omega \\ 
\langle \mathcal{B}(x, y),w-y\rangle_Y+\phi(w)-\phi(y)\geq\langle \zeta,w-y\rangle_Y, ~\forall w\in \mho
\end{cases}
\end{equation} 
studied in \cite{LiYaZeZh}.
\begin{definition}\label{d1.1}
Let $\Omega$ be a nonempty subset of a Banach space $X$. Let $\varphi\colon\Omega\to\overline{\mathbb{R}}$ be a proper convex and lower semicontinuous function, and $\mathcal{Q}\colon\Omega\to X^\ast$. Then $\mathcal{Q}$ is called
\begin{itemize}
\item[(i)]monotone, if it holds 
$$\langle \mathcal{Q}u-\mathcal{Q}v, u-v\rangle_X\geq0,~ \forall u,v\in\Omega;$$
\item[(ii)] strictly monotone, if it holds 
$$\langle \mathcal{Q}u-\mathcal{Q}v, u-v\rangle_X>0,~ \forall u,v \in\Omega\text{~~and~} u\neq v;$$
\item[(iii)]strongly monotone with constant	$\alpha_\mathcal{Q}>0,$ if it holds
$$\langle \mathcal{Q}u-\mathcal{Q}v, u-v\rangle_X\geq\alpha_\mathcal{Q}\|u-v\|^2_X,~\forall u,v\in\Omega;$$
\item[(iv)]relaxed strongly monotone with constant $\alpha_\mathcal{Q}>0,$ if it holds
$$\langle \mathcal{Q}u-\mathcal{Q}v, u-v\rangle_X\geq-\alpha_\mathcal{Q}\|u-v\|^2_X,~\forall u,v\in\Omega;$$
\item[(v)]inverse relaxed monotone with constant $\alpha_\mathcal{Q}>0,$ if it holds
$$\langle \mathcal{Q}u-\mathcal{Q}v, u-v\rangle_X\geq-\alpha_\mathcal{Q}\|Qu-Qv\|^2_X,~\forall u,v\in\Omega;$$
\item[(vi)]Lipschitz continuous with constant $\beta_\mathcal{Q}>0,$ if it holds
$$\| \mathcal{Q}u-\mathcal{Q}v\|_X\leq\beta_\mathcal{Q}\|u-v\|^2_X,~\forall u,v\in\Omega;$$
\item[(vii)]pseudomonotone, if for any $u,v\in\Omega$ we have 
$$\langle \mathcal{Q}u, v-u\rangle_X	\geq 0$$ 
then it entails that 
$$\langle \mathcal{Q}v, v-u\rangle_X\geq 0;$$
\item[(viii)]stable pseudo monotone with respect to the set $\mathcal{W}\subset X^\ast$, if $\mathcal{Q}$ and $u\rightarrowtail \mathcal{Q}u-w$ are pseudo monotone for all $w\in\mathcal{W};$
\item[(ix)]	$\varphi$-pseudomonotone, if for any $u,v\in\Omega$ we have 
$$\langle \mathcal{Q}u, v-u\rangle_X +\varphi(v)-\varphi(u)\geq 0$$
then it entails that 
$$\langle \mathcal{Q}v, v-u\rangle_X+\varphi(v)-\varphi(u)\geq 0;$$
\item[(x)] stable $\varphi$-pseudomonotone with respect to the set $\mathcal{W}\subset X^\ast$, if $\mathcal{Q}$ and $u\rightarrowtail \mathcal{Q}u-w$ are $\varphi$-pseudomonotone for each $w \in\mathcal{W}.$
\end{itemize}
\end{definition}
\begin{definition}\label{d1.1} \cite{Cl}
Let  $\varphi\colon X\to\bar{\mathbb{R}}=\mathbb{R}\cup\{+\infty\}$ be a proper, convex and lower semicontinuous function. The mapping $\partial_c\varphi\colon X\to 2^{X^\ast}$ defined by
$$\partial_c\varphi(u)=\{u^\ast\in X^\ast \vert \langle u^\ast, v-u\rangle_X\leq \varphi(v)-\varphi(u),~ \forall v \in	X\}$$
for $u \in X$, is called the subdifferential of $\varphi$. An element $x^\ast\in \partial_c\varphi(x)$ is called a subgradient of $\varphi$ in $u$.
\end{definition}
\begin{definition}\label{d1.2}\cite{Cl} 
Given a locally Lipschitz function $\varphi\colon X\to \mathbb{R}$, we denote by $\varphi^0(u;v)$ the Clarke generalized directional derivative of $\varphi$ at the point $u\in X$ in the direction $v\in X$ defined by
$$\varphi^0(u;v)=\limsup_{\lambda\to0^{+}, \eta\to u}\dfrac{\varphi(\eta+\lambda v)-\varphi(\eta)}{\lambda}.$$
The Clarke subdifferential or the generalized gradient of $\varphi$ at $u\in X$, denoted by $\partial \varphi(u)$, is a subset of $X^\ast$ given by
$$\partial \varphi(u)=\{u^\ast\in X^\ast\vert \varphi^0(u;v)\geq\langle u^\ast, v\rangle_X, \forall v\in X\}.$$
\end{definition}
Let $Z$ and $Y$ be topological spaces and $V$ be a nonempty subset of $Z$. We denote by $2^V$ the collection of subsets of $V$. Given a set-valued mapping $\mathcal{B}\colon Z\to 2^Y$, we use $Gr(\mathcal{B})$ to represent the graph of $\mathcal{B}$. In other word, $Gr(\mathcal{B})$ is defined as follows:
$$Gr(\mathcal{B})=\{(x,y) \in Z \times Y\vert y\in \mathcal{B}(x)\}\subset Z\times Y.$$
We say that the graph of $\mathcal{B}$ is sequentially closed in $Z\times Y$, or $\mathcal{B}$ is sequentially closed, if for any sequence $\{(x_n,y_n)\}\subset  Gr(\mathcal{B})$. If this sequence converges  to
$ (x,y)\in Z\times Y$ as $n\to\infty,$ 
then we have
$$(x,y)\in Gr (\mathcal{B})~~(\textit{i.e.,}~~ y\in \mathcal{B}(x)).$$
\begin{theorem}\label{t1.1}\cite{PaKy},
Let $Y$ be a reﬂexive Banach space and	$D$ be a nonempty, bounded, closed and convex set subset of  $Y$. Let $\Lambda\colon D\to 2^D$ be a set-valued map with nonempty, closed and convex values such that its graph is sequentially closed in $Y_w \times  Y_w$ topology. Then, $\Lambda$ has a ﬁxed point.
\end{theorem}
%------------------------------------------------------------------------------------------------------------------------------------------------------------------
\section{Main results}

In this section, we focus on the uniqueness of solutions and their existence to Problem \ref{P1.1}. We use the Minty methodology, Theorem \ref{t1.1}, and the Kakutani-Ky Fan fixed point theorem to establish the existence theorem for the solutions for the system of nonlinear mixed variational inequality problems under given modest assumptions. Additionally, we deploy the inverse relaxed monotonicity and Lipschitz continuity argument to present two uniqueness results for Problem \ref{P1.1}.\par
Furthermore, we propose the set-valued mappings $\mathcal{S}\colon\mho\to 2^\Omega$ and $\mathcal{T}:\Omega\to 2^\mho$, described by
$$\mathcal{S}(y)=\{x\in\Omega\vert\; x\text{~is a solution of \eqref{eq:1.1} corresponding to~} y\}, \forall y \in \mho,$$
and
$$\mathcal{T}(x)= \{y\in \mho\vert\; y \text{~is a solution of  \eqref{eq:1.2} corresponding to~} x\}, \forall x \in \Omega,$$
\textit{respectively}.\par
To investigate the existence of solutions to Problem \ref{P1.1}, we now suggest the following appropriate assumptions.
\begin{itemize}
\item[\textbf{(A):}] $\Omega\subset X$ and $\mho\subset Y$ are both nonempty, closed and convex.
\item[\textbf{(B):}] $\gamma\in X^\ast$ and $\zeta\in Y^\ast$.
\item[\textbf{(C):}] $\varphi\colon X\times X \to \overline{\mathbb{R}}$ is such that
\begin{itemize}
\item[(i)]$\varphi(\eta,\cdot)\colon X \to \overline{\mathbb{R}}$ is a proper, convex and lower semicontinuous function,
\item[(ii)] there exists $\varrho_\varphi\geq0$ such that
$$\varphi(\eta_1,y_2)-\varphi(\eta_1,y_1)+\varphi(\eta_2,y_1)-\varphi(\eta_2,y_2)\leq \varrho_\varphi\|\eta_1-\eta_2\|_X \|y_1-y_2\|_X, \forall \eta_1, \eta_2, y_1,y_2\in X,$$
\item[(iii)] for each $\eta\in X$, there exists $\varrho_\varphi(\eta)>0$ such that
$$\varphi(\eta,y_1)-\varphi(\eta,y_2)\leq\varrho_\varphi(\eta)\|y_1-y_2\|_X, \forall y_1,y_2\in X, \text{{\it see} \cite{CaSoXi}}.$$
\end{itemize}
\item[\textbf{(D):}] $\mathcal{A}:Y\times X\to X^\ast$	is such that
\begin{itemize}
\item[(i)] for each $y\in Y, x \rightarrowtail \mathcal{A}(y,x)$ is stable $\varphi$-pseudomonotone with  $\{\gamma\}$ and fulfills
$$\lim_{\lambda\to 0} \sup\langle \mathcal{A}(y, \lambda v +(1-\lambda)x), v-x\rangle_X\leq\langle \mathcal{A}(y, x), v-x\rangle_X, \forall y\in Y \text{~and~} x,v \in X;$$ 
\item[(ii)] it holds
$$\lim_{n\to\infty} \sup\langle \mathcal{A}(y_n,v), v-x_n\rangle_X\leq\langle \mathcal{A}(y,v), v-x\rangle_X,$$
whenever $v\in X, (x,y)\in X\times Y, \{y_n\}\subset Y	\text{~and~} \{x_n\}\subset X$ are such that
$$y_n\xlongrightarrow[]{w}y \in	Y \text{~~and~~} x_n\xlongrightarrow[]{w} x \in X \text{~as~} n\to\infty;$$
\item[(iii)] there exists a function $r\colon\mathbb{R}^{+}\times \mathbb{R}^{+}\to \mathbb{R}$ such that
$$\langle \mathcal{A}(y,x), x\rangle_ X\geq r (\|x\|_X, \|y\|_Y)\|x\|_X, \forall x\in X \text{~and~} y\in Y,$$
and\par
$\spadesuit$ every nonempty and bounded set $D\subset \mathbb{R}^{+},$ we have 
$$r(t,s)\to +\infty \text{~as~}  t\to +\infty, \forall s\in  D,$$
$\spadesuit$ for any constants $\varrho_1, \varrho_2 \geq 0,$ it holds $r(t, \varrho_1t + \varrho_2)\to +\infty$ as $t\to +\infty.$ 
\item[(iv)] there exists a constant $\varrho_\mathcal{A} > 0$ such that
$$ \|\mathcal{A}(y, x)\|_{X^\ast}\leq \varrho_\mathcal{A} (1+\|x\|_X+\|y\|_Y), \forall (x,y)\in X\times Y.$$
\end{itemize}
\item[\textbf{(E):}] $\phi\colon Y\times Y \to \overline{\mathbb{R}}$ is such that
\begin{itemize}
\item[(i)]$\phi(\eta,\cdot)\colon Y \to \overline{\mathbb{R}}$ is a proper, convex and lower semicontinuous function,
\item[(ii)] there exists $\varrho_\phi\geq0$ such that
$$\phi(\eta_1,y_2)-\phi(\eta_1,y_1)+\phi(\eta_2,y_1)-\phi(\eta_2,y_2)\leq \varrho_\phi\|\eta_1-\eta_2\|_Y \|y_1-y_2\|_Y, \forall \eta_1, \eta_2, y_1, y_2 \in Y,$$
\item[(iii)] for each $\eta\in X$, there exists $\varrho_\phi(\eta)>0$ such that
$$\phi(\eta,y_1)-\phi(\eta,y_2)\leq\varrho_\phi(\eta)\|y_1-y_2\|_X, \forall y_1,y_2\in X, \text{{\it see} \cite{CaSoXi}}.$$
\end{itemize}
\item[\textbf{(F):}] $\mathcal{B}\colon X\times Y\to Y^\ast$ is such that
\begin{itemize}
\item[(i)] for each $x\in X, y\rightarrowtail \mathcal{B}(x,y)$ is stable $\phi$-pseudomonotone with respect to $\{\zeta\}$ and satisﬁes
$$\lim_{\lambda\to 0} \sup\langle  \mathcal{B}(x, \lambda w +(1-\lambda)y), w-y\rangle_Y\leq\langle  \mathcal{B}(x, y), w-y\rangle_Y, \forall w, y\in Y \text{~and~} x \in X;$$ 
\item[(ii)] it holds
$$\lim_{n\to\infty} \sup \langle \mathcal{B}(x_n,w), w-y_n\rangle_Y\leq \langle \mathcal{B}(x,w),w -y\rangle_Y,$$
whenever $w\in Y, (x,y)\in X\times Y, \{y_n\}\subset Y$ and $\{x_n\}\subset X$	are such that
$$y_n\xlongrightarrow[]{w} y \in Y \text{~~and~} x_n\xlongrightarrow[]{w}x\in X\text{~~as~} n\to\infty;$$
\item[ (iii)]	there exists a function $\ell:\mathbb{R}^{+}\times \mathbb{R}^{+}\to\mathbb{R}$ such that
$$\langle \mathcal{B}(x,y), y\rangle_Y\geq \ell(\|y\|_Y,\|x\|_X)\|y\|_Y, \forall x\in X \text{~and~} y \in Y,$$
and\par
$\spadesuit$ every nonempty and bounded set $D\subset \mathbb{R}^{+}$, we	 have	
$$\ell (t,s)\to +\infty	\text{~as~} t\to +\infty, \text{~for all~} s\in D,$$
$\spadesuit$ for any constants $\varrho_1, \varrho_2 \geq 0$, it holds $\ell(t,\varrho_1t + \varrho_2)\to+\infty$ as $t\to +\infty.$ 
\item[(iv)] there exists a constant $\varrho_\mathcal{B} > 0$ such that
$$\|\mathcal{B}(x,y)\|_{Y^\ast}\leq\varrho_\mathcal{B}(1 +	\|x\|_X+\| y\|_Y), \forall (x,y)\in X\times Y.$$
\end{itemize}
\end{itemize}
\begin{remark}\label{r2.1} 
If function $r$ given in \textbf{(D)}(iii) (resp. $\ell$ given in \textbf{(F)}(iii)) is independent of its second variable, then condition \textbf{(D)}(iii) (resp. \textbf{(F)}(iii)) reduces to the following uniformly coercive condition:\par
\noindent
 \textbf{(D)}(iii)$^\prime$: there exists a function $r\colon\mathbb{R}^{+}\to \mathbb{R}$ with $r(s) \to + \infty$ as $s\to +\infty$ such that
$$\langle \mathcal{A}(y,x), x\rangle_X \geq r(\|x\|_X)\|x\|_X, \forall x\in X \text{~and~} y\in Y$$
(resp. \textbf{(F)}(iii)$^\prime$: there exists a function $\ell \colon \mathbb{R}^{+}\to \mathbb{R}$ with $\ell(s) \to +\infty $ as $ s\to +\infty$ such that
$$\langle \mathcal{B}(x,y), y\rangle_Y\geq\ell(\|y\|_Y)\|y\|_Y, \forall x\in X, \mbox{~and~} y\in Y).$$
\end{remark}
The main result of this paper, related to the existence of solutions to Problem \ref{P1.1}, is as follows.
\begin{theorem}\label{t2.1} 
 Suppose that \textbf{(A), (B), (C), (D), (E)} and \textbf{(F)} are held. Then, the solution set denoted by $\Gamma (\gamma,\zeta)$ of Problem \ref{P1.1} corresponding to $(\gamma,\zeta)\in X^\ast\times Y^\ast$ is nonempty and weakly compact in $X\times Y$.
 \end{theorem}
We require the following lemmas to prove this theorem.
\begin{lemma}\label{l2.1}
Assume that \textbf{(A), (B), (C)} and \textbf{(D)} are satisfied. Then, the following statements hold
\begin{itemize}
\item[(i)]for each fix $y\in Y$, $x\in\Omega$ is a solution of \eqref{eq:1.1}, if and only if, $x$ solves
the following Minty inequality for ﬁnding $x\in \Omega$ such that
\begin{equation}\label{eq:2.1}
\langle \mathcal{A}(y, v), v-x\rangle_X+\varphi(x,v)-\varphi(x,x)\geq\langle \gamma, v-x\rangle_X, \forall	v\in \Omega;
\end{equation}
\item[(ii)]for each fix $y\in Y$, the solution set denoted by $\mathcal{S}(y)$ of \eqref{eq:1.1} is nonempty, bounded, closed and convex;
\item[(iii)] the graph of the set-valued mapping $\mathcal{S}\colon \mho\to 2^\Omega$ is sequentially closed in $Y_w \times X_w,$ \textit{i.e.}, $\mathcal{S}$ is sequentially closed from $Y$ endowed with the weak topology into the subsets of $X$ with the weak topology;
\item[(iv)]	for each fix $y\in Y$, if the mapping $x \rightarrowtail \mathcal{A}(y,x)$ is strictly monotone, then $\mathcal{S}$ is a single-valued mapping and weakly continuous.
\end{itemize}
\end{lemma}
\begin{proof} The assumptions (i) and (ii) are the straightforward consequences of [\cite{LiMoZe}, Theorem 3.3] and [\cite{LiMoZe18}, Lemma 3.3].
Next, we present the conclusion (iii).\par
Let $\{(y_n,x_n)\} \subset Gr (\mathcal{S})$ be such that
\begin{equation}\label{eq:2.2}
y_n\xlongrightarrow []{w} y \in Y	\text{~~and~} x_n\xlongrightarrow[]{w} x\in X \text{~~as~} n\to\infty \text{~for~} (x,y)\in X\times Y.
\end{equation} 
Then, for each $n\in  \mathbb{N}$, we have $x_n\in \mathcal{S}(y_n),$ \textit{i.e.},
$$\langle \mathcal{A}(y_n,x_n), v-x_n\rangle_X + \varphi(x_n,v)-\varphi(x_n,x_n) \geq\langle \gamma, v-x_n\rangle_X,  \forall v \in \Omega.$$ 
The assertion (i) asserts that
\begin{equation}\label{eq:2.3}
\langle \mathcal{A}(y_n, v), v-x_n\rangle_X+\varphi(x_n,v)-\varphi(x_n,x_n) \geq \langle \gamma, v-x_n\rangle_X, \forall v\in \Omega.
\end{equation} 
To establish the upper limit as $n\to\infty$, use assumption \textbf{(D)}(ii) and weak lower semicontinuity of $\varphi$ (due to the convexity and lower semicontinuity of $\varphi$) to determine
\begin{align*}
\langle \mathcal{A}(y,v), v-x\rangle_X&+\varphi(x,v)-\varphi(x,x)
\geq \lim_{n\to\infty} \sup\langle \mathcal{A}(y_n,v), v -x_n\rangle_X\\
&\quad+\lim_{n\to\infty} \inf[\varphi(x_n,v)] -\lim_{n\to\infty} \inf [\varphi(x_n,x_n)]\\
&\geq\lim_{n\to\infty} \sup [\langle \mathcal{A}(y_n,v), v-x_n\rangle_X+\varphi(x_n,v)-\varphi(x_n,x_n)]\\
&\geq\lim_{n\to \infty} \sup\langle \gamma, v-x_n\rangle_X\\
&=\langle \gamma, v-x\rangle_X, \forall v\in \Omega.
\end{align*}
Using the assumption (i) again, we obtain 
$$x\in \mathcal{S}(y).$$ 
Consequently, $(y,x)\in Gr (\mathcal{S})$, \textit{i.e.,} the graph of the set-valued mapping $\mathcal{S}:\mho\in 2^\Omega$ is sequentially closed in $Y_w \times  X_w.$\par
Additionally, assume that $x\rightarrowtail \mathcal{A}(y,x)$ is strictly monotone. Let us consider $x_1, x_2 \in \Omega$ be two solutions to \eqref{eq:1.1}. Then, we have
\begin{equation}\label{eq:2.4}
\langle \mathcal{A}(y,x_1 ), v-x_1\rangle_X+\varphi(x_1,v)-\varphi(x_1,x_1)\geq\langle \gamma, v-x_1\rangle_X, \forall v\in \Omega
\end{equation} 
and 
\begin{equation}\label{eq:2.5}
\langle \mathcal{A}(y,x_2 ), v-x_2\rangle_X+\varphi(x_2,v)-\varphi(x_2,x_2)\geq\langle \gamma, v-x_2\rangle_X, \forall v\in \Omega.
\end{equation} 
Putting $v=x_2$ in \eqref{eq:2.4} and $v=x_1$ in \eqref{eq:2.5}, we have 
\begin{equation}\label{eq:2.6}
\langle \mathcal{A}(y,x_1 ), x_2-x_1\rangle_X+\varphi(x_1,x_2)-\varphi(x_1,x_1)\geq\langle \gamma, x_2-x_1\rangle_X,
\end{equation} 
and 
\begin{equation}\label{eq:2.7}
\langle \mathcal{A}(y,x_2 ), x_1-x_2\rangle_X+\varphi(x_2,x_1)-\varphi(x_2,x_2)\geq\langle \gamma, x_1-x_2\rangle_X.
\end{equation} 
Adding \eqref{eq:2.6} and \eqref{eq:2.7}, we have
\begin{equation}\label{eq:2.8}
\langle \mathcal{A}(y, x_1)-\mathcal{A}(y, x_2), x_1-x_2\rangle_X-\varphi(x_1,x_2)+\varphi(x_1,x_1)-\varphi(x_2,x_1)+\varphi(x_2,x_2)\leq0.
\end{equation}
Hence, from the assumption \textbf{(C)} and the strict monotonicity of $x \rightarrowtail \mathcal{A}(y,x)$ guarantees that $x_1 = x_2$. Therefore, $\mathcal{S}$ is a single-valued mapping. But, by virtue of assumption (iii), we can observe that $\mathcal{S}$ is weakly continuous.
\end{proof}
Similarly, for problem \eqref{eq:1.2}, we have the following lemma.
\begin{lemma}\label{l2.2}
Assume that \textbf{(A), (B), (E)} and \textbf{(F)} are satisfied. Then, the following statements hold
\begin{itemize}
\item[(i)] for each fix $x\in X$, $y\in \mho$ is a solution of \eqref{eq:1.2}, if and only if $y$ solves the following Minty inequality for ﬁnding $y\in	\mho$ such that
\begin{equation}\label{eq:2.9}
\langle \mathcal{B}(x,w), w-y\rangle_Y+\phi(y,w)-\phi(y,y) \geq\langle \zeta, w-y\rangle_Y, \forall w \in \mho;
\end{equation}
\item[(ii)] for each fix $x\in X$, the solution set $\mathcal{T}(x)$ namely, of \eqref{eq:1.2} is nonempty, bounded, closed and convex;
\item[(iii)] the graph of the set-valued mapping $\mathcal{T}\colon \Omega\to 2^\mho$ is sequentially closed in $X_w\times Y_w$;
\item[(iv)]	for each fix $x\in X$, if the mapping $y\rightarrowtail \mathcal{B}(x,y)$ is strictly monotone, then $\mathcal{T}$ is a single-valued mapping and weakly continuous.
\end{itemize}
\end{lemma}
Furthermore, we focus on providing a priori estimates for the solutions of Problem \ref{P1.1}.
\begin{lemma} \label{l2.3} 
Assume that \textbf{(A), (B), (C), (D), (E)} and \textbf{(F)} are satisfied. If the solution set $\Gamma(\gamma,\zeta)$ namely, of Problem \ref{P1.1} is nonempty, then there exists a constant $M > 0$ such that
\begin{equation}\label{eq:2.10} 
\|x\|_X \leq M \text{~~and~}  \|y\|_Y \leq M,~ \forall (x,y) \in \Gamma(\gamma,\zeta).
\end{equation}
\end{lemma}
\begin{proof} 
Suppose that $\Gamma(\gamma,\zeta)\neq \emptyset$. Let $(x,y)\in \Gamma(\gamma,\zeta)$ be arbitrary and $(x_0,y_0)\in (D(\varphi) \cap \Omega)\times(D(\phi)\cap \mho).$ By swapping $v = x_0$ and $w=y_0$ into \eqref{eq:1.1} and \eqref{eq:1.2}, \textit{respectively}, we obtain
\begin{equation}\label{eq:2.11}
\langle \mathcal{A}(y,x), x\rangle_X\leq\langle \mathcal{A}(y,x), x_0\rangle_X+\varphi(x,x_0)-\varphi(x,x) +\langle \gamma, x_0-x\rangle_X
\end{equation}
and
\begin{equation}\label{eq:2.12}
\langle \mathcal{B}(x,y),y\rangle_Y\leq \langle \mathcal{B}(x,y), y_0\rangle_Y+\phi(y,y_0)-\phi(y,y) +\langle \zeta, y_0-y\rangle_Y.
\end{equation}
%Recall that $\varphi$ and $\phi$ are both proper, convex and l.s.c., it follows from Proposition 2 that there are constants $\alpha_\phi, \alpha_\varphi, \beta_\phi, \beta_varphi\geq 0$ satisfying
%\begin{equation}\label{eq:2.8}
%\varphi(v) \geq-\alpha_\varphi\|v\|_X-\beta_\varphi \text{~and~}\phi(w)\geq-\alpha_\phi\|w\|_Y-\beta_\phi\forall (v,w) \in X \times Y
%\end{equation}
Taking account of \eqref{eq:2.11}, we use hypotheses \textbf{(C)}(i),(iii) and \textbf{(D)}(iii)–(iv) to obtain
\begin{align*}
r (\|x\|_X,\|y\|_Y )\|x\|_X &\leq\langle \mathcal{A}(y, x), x\rangle_X\\
&\leq  \langle \mathcal{A}(y, x), x_0\rangle_ X+\varphi(x,x_0)-\varphi(x,x)+\langle \gamma, x_0-x \rangle_X\\
&\leq\|\mathcal{A}(y, x)\|_{X^\ast}\| x_0\|_X+\varphi(x,x_0)-\varphi(x,x)+\|\gamma\|_{X^\ast} (\|x_0\|_X+\|x\|_X)\\
&\leq \varrho_\mathcal{A} (1+\|x\|_X+\|y\|_Y)\|x_0\|_X+\varrho_\varphi(x)(\|x_0\|_X+\|x\|_X)+\|\gamma\|_{X^\ast} (\|x_0\|_X+\|x\|_X ).
\end{align*}
This implies that
\begin{equation}\label{eq:2.13}
r (\|x\|_X ,\|y\|_Y)\leq\dfrac{\varrho_\mathcal{A} (1+\|x\|_X+\|y\|_Y)\|x_0\|_X}{\|x\|_X }+\dfrac{(\varrho_\varphi(x)+\|\gamma\|_{X^\ast})\|x_0\|_X}{\|x\|_X}+ \varrho_\varphi(x)+\|\gamma\|_{X^\ast}.
\end{equation}
Similarly, taking account of \eqref{eq:2.12}, we use hypotheses  \textbf{(E)}(i), (iii) and  \textbf{(F)}(iii)–(iv) to obtain
\begin{align}\label{eq:2.14}
\ell(\|y\|_Y,\| x\|_X)\|y\|_Y
&\leq \varrho_\mathcal{B} (1+ \|x\|_X + \|y\|_Y)\|y_0\|_Y +\varrho_\phi(y)(\|y_0\|_Y+\|y\|_Y)+\|\zeta\|_{Y^\ast} (\|y_0\|_Y+\| y\|_Y)\notag\\
\Longrightarrow&\notag\\
\ell(\|y\|_Y ,\| x\|_X)&\leq\dfrac{\varrho_\mathcal{B} (1+\|x\|_X+\|y\|_Y)\|y_0\|_Y}{\|y\|_Y }+\dfrac{(\varrho_\phi(y)+\|\zeta\|_{Y^\ast})\|y_0\|_Y}{\|y\|_Y}\notag\\
&\quad\quad+ \varrho_\phi(y)+\|\zeta\|_{Y^\ast}.
\end{align}
Contrary, suppose $\Gamma(\gamma,\zeta)$ is unbounded. Then, taking a subsequence if necessary, it is possible to ﬁnd a sequence $\{(x_n, y_n)\}\subset\Omega \times \mho$ such that it holds
\begin{equation}\label{eq:2.15}
\|x_n\|_X	\uparrow +\infty	\text{~as~} n \to\infty,
\end{equation}
or
\begin{equation}\label{eq:2.16}
\|y_n\|_Y	\uparrow +\infty	\text{~as~} n \to\infty.
\end{equation}
Let us segregate the subsequent cases:\par
$\textcircled{a}$ Assuming that \eqref{eq:2.15} is fulfilled and the sequence $\{y_n\}$ is bounded in the space $Y$,\par
$\textcircled{b}$ Assuming that \eqref{eq:2.16} is fulfilled and the sequence $\{x_n\}$ is bounded in the space $X$,\par
$\textcircled{c}$ Assuming that both \eqref{eq:2.15} and \eqref{eq:2.16} are fulfilled.\par
Assume that $\textcircled{a}$ is valid, then we put $x=x_n$  and $y=y_n$	in \eqref{eq:2.13} to get
\begin{align}\label{eq:2.17}
r (\|x_n\|_X ,\|y_n\|_Y)&\leq\dfrac{\varrho_\mathcal{A} (1+\|x_n\|_X+\|y_n\|_Y)\|x_0\|_X}{\|x_n\|_X }+\dfrac{(\varrho_\varphi(x)+\|\gamma\|_{X^\ast})\|x_0\|_X}{\|x_n\|_X}\notag\\
&\quad\quad+ \varrho_\varphi(x)+\|\gamma\|_{X^\ast}.
\end{align}\
When we let $n$ approach infinity in the inequality \eqref{eq:2.17} and make use of \eqref{eq:2.15} along with property \textbf{(D)}(iii), we obtain the following:
\begin{align}\label{eq:2.18}
+\infty&=\lim_{n\to\infty}r (\|x_n\|_X ,\|y_n\|_Y)\notag\\
&\leq \lim_{n\to\infty}\left[\dfrac{\varrho_\mathcal{A} (1+\|x_n\|_X+\| y_n\|_Y)\|x_0\|_X}{\|x_n\|_X}
+\dfrac{(\varrho_\varphi(x)+\|\gamma\|_{X^\ast})\| x_0\|_X}{\|x_n\|_X}+\varrho_\varphi(x)+\|\gamma\|_{X^\ast}\right]\notag\\
&=\varrho_\mathcal{A}\|x_0\|_X+\varrho_\varphi(x)+\|\gamma\|_{X^\ast}.
\end{align}
Hence, \eqref{eq:2.18} generates a contradiction. Similarly, for the case $\textcircled{b}$ we could also use \eqref{eq:2.14} to get a contradiction. However, we assume that $\textcircled{c}$ hold, we will proceed to discuss two additional situations:
\begin{itemize}
\item[(1)]	$\dfrac{\|y_n\|_Y}{\|x_n\|_X}\to+\infty	\text{~~as~~} n\to\infty;$
\item[(2)]	there exist $n_0\in \mathds{N}$ and $\widehat{\varrho_0}>0$ such that
 $$\dfrac{\|y_n\|_Y}{\|x_n\|_X}\leq\widehat{\varrho_0},~\forall n\geq n_0.$$
\end{itemize}
If the item (1) is true, then we put $x=x_n$ and $y=y_n$ into \eqref{eq:2.14} to yield
\begin{equation*}
\ell (\|y_n\|_Y, \|x_n\|_X)\leq \dfrac{\varrho_\mathcal{B}(1+\|x_n\|_X+\|y_n\|_Y)\|y_0\|_Y}{\|y_n\|_Y}
+\dfrac{(\varrho_\phi(y)+\|\zeta\|_{Y^\ast})\|y_0\|_Y}{\|y_n\|_Y}
+\varrho_\phi(y)+\|\zeta\|_{Y^\ast}.
\end{equation*}
Taking the limit as $n$ approaches infinity for the inequality mentioned above yields:
\begin{align}\label{eq:2.19}
+\infty&=\lim_{n\to\infty}\ell (\|y_n\|_Y, \|x_n\|_X)\notag\\
&\leq \lim_{n\to\infty}\left[\dfrac{\varrho_\mathcal{B}(1+\|x_n\|_X+\|y_n\|_Y)\|y_0\|_Y}{\|y_n\|_Y}
+\dfrac{(\varrho_\phi(y)+\|\zeta\|_{Y^\ast})\|y_0\|_Y}{\|y_n\|_Y}+\varrho_\phi(y)+\|\zeta\|_{Y^\ast}\right]\notag\\
&=\varrho_\mathcal{B}\|y_0\|_Y+\varrho_\phi(y)+\|\zeta\|_{Y^\ast}.
\end{align}
It is obviously impossible, however for a situation (2), we can deduce from \eqref{eq:2.13} that
\begin{align}\label{eq:2.20}
+\infty&\leftarrow r(\|x_n\|_X, \|y_n\|_Y )~~~(\text{as~} n\to\infty)\notag\\
&\leq \dfrac{\varrho_\mathcal{A} (1+\|x_n\|_X+\|y_n\|_Y)\|x_0\|_X}{\|x_n\|_X}+\dfrac{(\varrho_\varphi(x) +\|\gamma\|_{X^\ast})\|x_0\|_X}{\|x_n\|_X}+\varrho_\varphi(x)+\|\gamma\|_{X^\ast}\notag\\
&=\varrho_\mathcal{A}(2+\widehat{\varrho}_0)\|x_0\|_X+\varrho_\varphi(x)+\|\gamma\|_{X^\ast}\|x_0\|_X+\varrho_\varphi(x)+\|\gamma\|_{X^\ast}, \text{~for}~ n\geq n_1,
\end{align}
where $n_1\geq n_0$ is such that 
$$\|x_{n_1}\|_ X > 1.$$ 
This leads to a contradiction. Thus, we conclude that $\Gamma(\gamma,\zeta)$ is bounded in $X \times Y$, allowing us to find a constant $M > 0$ satisfying \eqref{eq:2.10}.
\end{proof}
Consider the set-valued mapping $\Lambda:\Omega\times \mho\to2^{\Omega\times \mho}$ given by
\begin{equation}\label{eq:2.21}
\Lambda(x,y)=(\mathcal{S}(y),\mathcal{T}(x)),  \forall (x,y)\in \Omega\times \mho.
\end{equation}
By invoking Lemma \ref{l2.1} and Lemma \ref{l2.2}, it can be seen that $\Lambda$ is well-defined. Additionally, there exists a bounded, closed, and convex set $\mathfrak{D}$ in $\Omega \times \mho$ such that $\Lambda$ maps $\mathfrak{D}$ into itself.
\begin{lemma}\label{l2.4}
Assume that conditions \textbf{(A), (B), (C), (D), (E)} and \textbf{(F)} are satisfied. Then, there exists	a constant $\widehat{M}>0$ satisfying	$\Lambda\overline{(\mathscr{B}(0,\widehat{M}))}\subset \overline{\mathscr{B}(0,\widehat{M})}$, where $\overline{\mathscr{B}(0, \widehat{M})}\subset X\times Y$ is deﬁned by
$$\overline{\mathscr{B}(0, \widehat{M})}=	\{(x,y)\in\Omega\times \mho\, \vert\;\, \|x\|_X\leq \widehat{M} \text{~~and~~} \|y\|_Y\leq \widehat{M}\}.$$
\end{lemma}
\begin{proof} Our proof will be based on contradiction. Suppose for each $n\in\mathds{N}$, it holds $\Gamma(\overline{\mathscr{B}(0,n)})\not\subset \overline{\mathscr{B}(0,n)}$. Then, for every $n \in \mathds{N}$, we can find $(x_n, y_n) \in \overline{\mathscr{B}(0, n)}$ and $(z_n, w_n) \in \Gamma(x_n, y_n)$  ( \text{i.e.}, $z_n\in\mathcal{S}(y_n$)  and $w_n\in \mathcal{T}(x_n))$ such that
\begin{equation}\label{eq:2.22}
\|z_n\|_X > n \text{~~~~or~~~~} \|w_n\|_Y > n.
\end{equation}
Thus, passing to a relabeled subsequence if needed, we can assume that $\|z_n\|_X > n$ for each $n\in \mathds{N}$ (since the proof for the case that $\|w_n\|_Y > n$ for each $n\in \mathds{N}$ is similar). Employing \eqref{eq:2.13}, it determines
\begin{equation*}
r(\|z_n\|_X, \|y_n\|_Y)\leq\dfrac{\varrho_\mathcal{A}(1+\|z_n\|_X+\|y_n\|_Y)\|x_0\|_X}{\|z_n\|_X}
+\dfrac{(\varrho_\varphi(x) +\|\gamma\|_{X^\ast})\|x_0\|_X}{\|z_n\|_X}+\varrho_\varphi(x)+\|\gamma\|_{X^\ast}.
\end{equation*}
Since 
$$\|y_n\|_Y\leq n < \|z_n\|_X.$$ 
Therefore, passing to the limit as  $n\to\infty$ for the inequality above, it gives
\begin{align*}
+\infty&=\lim_{n\to\infty} r(\|z_n\|_X, \|y_n\|_Y)\\
&\leq\lim_{n\to\infty}\left[\dfrac{\varrho_\mathcal{A}(1+\|z_n\|_X+\|y_n\|_Y)\|x_0\|_X}{\|z_n\|_X}
+\dfrac{(\varrho_\varphi(x) +\|\gamma\|_{X^\ast})\|x_0\|_X}{\|z_n\|_X}+\varrho_\varphi(x)+\|\gamma\|_{X^\ast}\right]\\
&\leq 2\varrho_\mathcal{A} \|x_0\|_X+\varrho_\varphi(x)+\|\gamma\|_{X^\ast}.
\end{align*}
This leads to a contradiction. Thus, there exists a constant $\widehat{M}>0$ satisfying 
$$\Lambda(\overline{\mathscr{B}(0, \widehat{M})})\subset \overline{\mathscr{B}(0,\widehat{M})}.$$
\end{proof}
\begin{proof} \textbf{(Proof of Theorem \ref{t2.1})}
Let us see that if $(x^\ast,y^\ast)$ is a ﬁxed point of $\Lambda$, then we have $x^\ast \in \mathcal{S}(y^\ast)$ and $y^\ast\in \mathcal{T}(x^\ast)$. Using the deﬁnitions of $\mathcal{S}$ and $\mathcal{T}$, it provides
$$\langle \mathcal{A}(y^\ast, x^\ast), v-x^\ast\rangle_X+\varphi(x^\ast,v)-\varphi(x^\ast,x^\ast) \geq \langle \gamma, v-x^\ast\rangle_X, \forall v\in \Omega,$$
and
$$\langle \mathcal{B}(x^\ast, y^\ast), w-y^\ast\rangle_Y+\phi(y^\ast,w)-\phi(y^\ast,y^\ast)\geq \langle \zeta, w-y^\ast\rangle_Y, \forall w\in \mho.$$
Consequently, it is clear that $(x^\ast, y^\ast)$ is also a solution to Problem \ref{P1.1}. We will apply Theorem \ref{t1.1}, the Kakutani-Ky Fan fixed point theorem, to determine the existence of a fixed point for $\Lambda$.\par
Moreover, Lemma \ref{l2.1}, Lemma \ref{l2.2}, and Lemma \ref{l2.4}, in fact, infer that $\Lambda\colon\overline{\mathscr{B}(0,\widehat{M})}\to2^{\overline{\mathscr{B}(0,\widehat{M})}}$ has nonempty, closed, and convex values, and that the graph of $\Lambda$ is sequentially closed in $(X\times Y)_w \times (X\times Y)_w.$  The conditions stated in Theorem \ref{t1.1} have been verified. By using this theorem, there exists a solution $(x^\ast, y^\ast) \in \Omega \times \mho$ to Problem \ref{P1.1} such that 
$$(x^\ast, y^\ast) \in \Lambda(x^\ast, y^\ast)$$ 
and hence,
$$\Gamma(\gamma,\zeta) \neq\emptyset.$$
It is evident from Lemma \ref{l2.3} that $\Gamma(\gamma,\zeta)$ is bounded in $X\times Y$. Next, we will demonstrate that $\Gamma(\gamma,\zeta)$ is weakly closed. Let $\{(x_n, y_n)\}\subset\Gamma(\gamma,\zeta)$ be such that
\begin{equation}\label{eq:2.23}
(x_n,y_n)\xlongrightarrow[]{w}(x,y)\in X\times Y\text{~~as~~}n\to\infty,\text{~~for some~} (x,y)\in \Omega \times \mho.
\end{equation}
It is easy to see that for each natural number $n$, the pair $(x_n,y_n)\in \Lambda(x_n,y_n).$ Since $\Lambda$ is sequentially closed from $(X \times Y)_w$ to $(X \times Y)_w$ (\textit{see} Lemma \ref{l2.1} and Lemma \ref{l2.2}), we can conclude that $$(x,y)\in \Lambda(x,y).$$ 
This implies that 
$$(x,y) \in \Gamma(\gamma,\zeta).$$ 
Thus, due to the boundedness of $\Gamma(\gamma,\zeta)$, we can conclude that $\Gamma(\gamma,\zeta)$ is weakly compact.
\end{proof}
Theorem \ref{t2.1} shows that the solution set of Problem \ref{P1.1} is both non-empty and weakly compact. However, it raises the question of whether it is possible to prove the uniqueness of the solution under certain assumptions. Fortunately, the following theorems provide a positive answer to this issue.
\begin{theorem}\label{t2.2} 
Assume that \textbf{(A), (B),  (C), (D), (E)} and \textbf{(F)} are satisﬁed. If, in addition, the following inequality holds,
\begin{align}\label{eq:2.24}
&\langle \mathcal{A}(y_1, x_1)-\mathcal{A}(y_2, x_2), x_1-x_2\rangle_X+\langle \mathcal{B}(x_1, y_1)-\mathcal{B}(x_2, y_2), y_1-y_2\rangle_Y+\varrho_\varphi\|x_1-x_2\|^2_X\notag\\
&~~~~~~~~~+\varrho_\phi\|y_1-y_2\|^2_Y>0, \forall (x_1, y_1), (x_2, y_2) \in X\times Y \text{~with~} (x_1, y_1) \neq (x_2, y_2).
\end{align}Then Problem \ref{P1.1} has a unique solution.
\end{theorem}
\begin{proof} 
Theorem \ref{t2.1} assures that $\Gamma(\gamma,\zeta)\neq\emptyset$. We now prove the uniqueness of Problem \ref{P1.1}. Let $(x_1, y_1), (x_2, y_2) \in \Gamma(\gamma,\zeta)$. Then, we have
\begin{equation}\label{eq:2.25}
\langle \mathcal{A}(y_i, x_i), v-x_i\rangle_X +\varphi(x_i,v) -\varphi(x_i,x_i) \geq \langle \gamma, v - x_i\rangle_X,  \forall v\in\Omega,
\end{equation} 
and
\begin{equation}\label{eq:2.26}
\langle \mathcal{B}(x_i, y_i), w-y_i\rangle_Y+\phi(y_i,w)-\phi(y_i,y_i) \geq \langle \zeta, w-y_i\rangle_Y,  \forall w\in \mho.
\end{equation} 
After setting $i=1$ to correspond to $v=x_2$  and $i=2$ to correspond to $v=x_1$ in equation \eqref{eq:2.25}, we add the two equations to obtain
\begin{equation}\label{eq:2.27}
\langle \mathcal{A}(y_1, x_1)-\mathcal{A}(y_2, x_2), x_1-x_2\rangle_X-\varphi(x_1,x_2) +\varphi(x_1,x_1)-\varphi(x_2,x_1) +\varphi(x_2,x_2)\leq0.
\end{equation}
Similarly, assigning $i=1$ to correspond to $w=y_2$  and $i=2$ to correspond to $w=y_1$ in equation \eqref{eq:2.26}, we add the two equations to get
\begin{equation}\label{eq:2.28}
\langle \mathcal{B}(x_1, y_1)-\mathcal{B}(x_2, y_2), y_1-y_2\rangle_Y-\phi(y_1,y_2) +\phi(y_1,y_1)-\phi(y_2,y_1) +\phi(y_2,y_2)\leq0.
\end{equation}
By using the \eqref{eq:2.27}, \eqref{eq:2.28}, and assertions \textbf{(C)} and \textbf{(E)}, we have
$$\langle \mathcal{A}(y_1, x_1)-\mathcal{A}(y_2, x_2), x_1-x_2\rangle_X+\langle \mathcal{B}(x_1, y_1)-\mathcal{B}(x_2, y_2), y_1-y_2\rangle_Y+\varrho_\varphi\|x_1-x_2\|^2_X+\varrho_\phi\|y_1-y_2\|^2_Y\leq0.$$
This combined with the condition \eqref{eq:2.24} implies that $x_1=x_2$ and $y_1=y_2$. Thus, Problem \ref{P1.1} has a unique solution.
\end{proof}
By adding an additional condition to \eqref{eq:2.24}, the resulting theorem establishes a unique solution for Problem \ref{P1.1}.
%---------------------------theorem 2.3---------------------------
\begin{theorem} \label{t2.3}
Assume that \textbf{(A), (B), (C), (D), (E)} and \textbf{(F)} are satisﬁed. If the following conditions also hold:
\begin{itemize}
\item[\textcircled{1}] for each $y\in Y$, the function $x\rightarrowtail \mathcal{A}(y,x)$ is inversely relaxed monotone and Lipschitz continuous with constants $\alpha_\mathcal{A} > 0$ and $\beta_\mathcal{A}>0$, \textit{respectively}, and for each $x\in X$ the function $y \rightarrowtail \mathcal{A}(y,x)$ is Lipschitz continuous with constant $\mathscr{L}_\mathcal{A} > 0$,
\item[\textcircled{2} ] for each $x\in X$, the function $y\rightarrowtail \mathcal{B}(x,y)$ is inversely relaxed  monotone and Lipschitz continuous with constants $\alpha_\mathcal{B} > 0$ and $\beta_\mathcal{B}>0$, \textit{respectively}, and for each $y\in Y$ the function $x\rightarrowtail \mathcal{B}(x,y)$ is Lipschitz continuous with constant $\mathscr{L}_\mathcal{B} > 0$,
\item[\textcircled{3}] $\dfrac{\mathscr{L}_\mathcal{A} \mathscr{L}_\mathcal{B}}{(\alpha_\mathcal{A}\beta_\mathcal{A}+\varrho_\varphi)(\alpha_\mathcal{B}\beta_\mathcal{B}+\varrho_\phi)} < 1.$
\end{itemize}
Then Problem \ref{P1.1} has a unique solution.
\end{theorem}
\begin{proof}	
Let $(x_1,y_1)$ and $(x_2, y_2)$ be two solutions to Problem \ref{P1.1}. Then, it has
\begin{equation}\label{eq:2.29}
\langle \mathcal{A}(y_1, x_1)-\mathcal{A}(y_2, x_2), x_1-x_2\rangle_X-\varrho_\varphi\|x_1-x_2\|^2_X\leq 0,
\end{equation}
and
\begin{equation}\label{eq:2.30}
\langle \mathcal{B}(x_1, y_1)-\mathcal{A}(x_2, y_2), y_1-y_2\rangle_Y-\varrho_\phi\|y_1-y_2\|^2_Y\leq 0.
\end{equation}
Again, from the inverse relaxed monotonicity and Lipschitz continuity of $\mathcal{A}$, we have
\begin{align}\label{eq:2.31}
-\alpha_\mathcal{A}\beta_\mathcal{A}\|x_1-x_2\|^2_X&\leq \langle \mathcal{A}(y_1, x_1-\mathcal{A}(y_1, x_2), x_1-x_2\rangle_X\notag\\
&\leq \langle \mathcal{A}(y_2, x_2)-\mathcal{A}(y_1, x_2), x_1-x_2\rangle_X\notag\\
&\leq \mathscr{L}_\mathcal{A}\|y_1-y_2\|_Y\|x_1-x_2\|_X.
\end{align}
Thus, from \eqref{eq:2.29} and \eqref{eq:2.31}, we have
\begin{equation}\label{eq:2.32}
(-\alpha_\mathcal{A}\beta_\mathcal{A}-\varrho_\varphi)\|x_1-x_2\|^2_X\leq \mathscr{L}_\mathcal{A}\|y_1-y_2\|_Y\|x_1-x_2\|_X.
\end{equation}
Similarly, we obtain
\begin{align}\label{eq:2.33}
-\alpha_\mathcal{B}\beta_\mathcal{B}\|y_1 - y_2\|^2_Y	&\leq \langle \mathcal{B}(x_1, y_1)-\mathcal{B}(x_1, y_2), y_1-y_2\rangle_Y\notag\\
&\leq \langle \mathcal{B}(x_2, y_2) -\mathcal{B}(x_1, y_2), y_1-y_2\rangle_Y\notag\\
&\leq \mathscr{L}_\mathcal{B}\|y_1-y_2\|_Y\|x_1-x_2\|_X.
\end{align}
Again, from \eqref{eq:2.30} and \eqref{eq:2.33}, we have
\begin{equation}\label{eq:2.34}
(-\alpha_\mathcal{B}\beta_\mathcal{B}-\varrho_\phi)\|y_1-y_2\|^2_Y\leq \mathscr{L}_\mathcal{B}\|y_1-y_2\|_Y\|x_1-x_2\|_X.
\end{equation}
Combining equations \eqref{eq:2.32} and \eqref{eq:2.34} yields 
\begin{equation}\label{eq:2.35}
\|x_1-x_2\|_X\leq \dfrac{\mathscr{L}_\mathcal{A} \mathscr{L}_\mathcal{B}}{(\alpha_\mathcal{A}\beta_\mathcal{A}+\varrho_\varphi)(\alpha_\mathcal{B}\beta_\mathcal{B}+\varrho_\phi)}\|x_1-x_2\|_X.
\end{equation}
However, the inequality $\frac{\mathscr{L}_\mathcal{A}\mathscr{L}_\mathcal{B}}{(\alpha_\mathcal{A}\beta_\mathcal{A}+\varrho_\varphi)(\alpha_\mathcal{B}\beta_\mathcal{B}+\varrho_\phi)} < 1$ implies that $x_1=x_2$	and	$y_1=y_2.$ Therefore, Problem \ref{P1.1} has a unique solution.
\end{proof}
%========================================================
\section{ Stability results}%for Coupled Mixed Variational Inequality Problems}
In this section, we delve into examining the stability of the system of nonlinear mixed variational inequality problems. Firstly, we present a set of regularized problems perturbed by duality mappings that correspond to Problem \ref{P1.1}. Secondly, we produce a stability result that shows that any solution sequence of a regularized problem has at least one subsequence that leads to a solution of the original problem, which is Problem \ref{P1.1}.\par
Recall that $X$ and $Y$ are two reflexive Banach spaces that can be renormed to become strictly convex. Without loss of generality, we may assume that	$X$	and $Y$ are strictly convex. Let $J_X\colon X\to X^\ast$ and $J_Y\colon Y\to Y^\ast$ be the duality mappings of the spaces $X$ and $Y$, \textit{respectively}, namely:
$$J_X(x)=\{x^\ast\in X^\ast\vert\langle x^\ast, x\rangle_X=\|x\|^2_X=\|x^\ast\|^2_{X^\ast}\},$$
$$J_Y(y)=\{y^\ast\in Y^\ast\vert\langle y^\ast , y\rangle_Y	=\|y\|^2_Y=\|y^\ast\|^2_{Y^\ast}\}.$$
Let real sequences $\{\varepsilon_n\}$ and $\{\delta_n\}$ be such that
\begin{equation}\label{eq:3.1}
\varepsilon_n > 0,~~ \delta_n>0,~~~\varepsilon_n\to 0 \text{~and~} \delta_n\to 0.
\end{equation}
For each $n\in \mathds{N},$ consider the following perturbated problem corresponding to Problem \ref{P1.1}.
\begin{problem}\label{P3.1}	
Find $(x_n, y_n) \in \Omega\times \mho$ such that
\begin{equation}\label{eq:3.2}
\langle \mathcal{A}(y_n, x_n)+\varepsilon_n J_X (x_n), v-x_n\rangle_X+\varphi(x_n,v)-\varphi(x_n,x_n)\geq\langle \gamma, v-x_n\rangle_X, \forall v\in\Omega,
\end{equation}
and
\begin{equation}\label{eq:3.3}
\langle \mathcal{B}(x_n, y_n) +\delta_n J_Y (y_n), w-y_n\rangle_Y+\phi(y_n,w)-\phi(y_n,y_n)\geq\langle \zeta, w-y_n\rangle_Y, \forall w \in \mho.
\end{equation}
\end{problem}
We make the following assumptions.
\begin{itemize}
\item[\textbf{(G):}]  $x\rightarrowtail \mathcal{A}(y,x)$ and $y\rightarrowtail \mathcal{B}(x,y)$ are monotone, and satisfy
$$\limsup_{\lambda\to 0} \langle \mathcal{A}(y, \lambda v +(1-\lambda)x), v-x\rangle_X\leq \langle \mathcal{A}(y, x), v-x\rangle_X,$$
$$\limsup_{\lambda\to 0} \langle \mathcal{B}(x,\lambda w +(1-\lambda) y), w-y\rangle_Y\leq\langle \mathcal{B}(x, y), w-y\rangle_Y, \forall w, y \in Y\text{~~and~} v, x\in X.$$
\item[\textbf{(H):}] $x \rightarrowtail \mathcal{A}(y,x)$ is inverse relaxed  monotone with constant $\alpha_\mathcal{A} > 0$  and Lipschitz continuous with constant $\beta_\mathcal{A} > 0$; similarly, $y\rightarrowtail \mathcal{B}(x, y)$ is inverse relaxed monotone with constant $\alpha_\mathcal{B} > 0$ and Lipschitz continuous with constant $\beta_\mathcal{B} > 0$, and and satisfy
$$\limsup_{\lambda\to 0} \langle \mathcal{A}(y,\lambda v +(1-\lambda)x), v-x\rangle_X\leq \langle \mathcal{A}(y, x), v-x\rangle_X,$$
$$\limsup_{\lambda\to 0} \langle \mathcal{B}(x, \lambda w +(1-\lambda) y), w-y\rangle_Y\leq\langle \mathcal{B}(x, y), w-y\rangle_Y, \forall w, y\in Y\text{~and~} v, x\in X.$$
\end{itemize}
The theorem below guarantees the existence and convergence of solutions to Problem \ref{P3.1}.
%-------------------------theorem 3.1----------------------
\begin{theorem}\label{t 3.1} 
Assume that \textbf{(A), (B), (C), (D)}(ii)–(iv), \textbf{(E)} and \textbf{(F)}(ii)–(iv) are satisfied. Then the following assertions hold,
\begin{itemize}
\item[(i)] if,	in addition, \textbf{(G)} holds, then for each $n\in\mathds{N}$, Problem \ref{P3.1} has at least one solution $(x_n, y_n)\in\Omega\times \mho;$
\item[(ii)] furthermore, if \textbf{(G)} holds, then for any sequence of solutions $\{(x_n, y_n)\}$ of Problem \ref{P3.1}, there exists a subsequence $\{(x_n, y_n)\}$, such that
\begin{equation}\label{eq:3.4}
(x_n, y_n)\xlongrightarrow[]{w}	(x, y) \in X\times Y \text{~~as~~} n\to\infty,
\end{equation}
where $(x,y)\in	\Omega\times \mho$ is a solution of Problem \ref{P3.1};
\item[(iii)] if \textbf{(H)} holds, then for any sequence of solutions $\{(x_n, y_n)\}$ of Problem \ref{P3.1}, there exists a subsequence $\{(x_n, y_n)\}$ such that
\begin{equation}\label{eq:3.5}
(x_n, y_n)\to(x,y) \in X\times Y \text{~~as~~} n\to\infty,
\end{equation}
where $(x, y)\in\Omega\times \mho$ is a solution of Problem \ref{P1.1}.
\end{itemize}
\end{theorem}
\begin{proof} 
\begin{itemize}
\item[(i)] Set 
$$\mathcal{A}_n(y,x) = \mathcal{A}(y,x) + \varepsilon_n J_X (x)$$ 
and 
$$\mathcal{B}_n(x, y)= \mathcal{B}(x,y)+\delta_n J_Y(y),~ \forall~(x, y) \in X \times Y.$$ 
We will verify that $\mathcal{A}_n$ and $\mathcal{B}_n$ satisfy  \textbf{(D)} and  \textbf{(F)}, \textit{respectively}. Note that $J_X$ is demicontinuous and
\begin{equation}\label{eq:3.6}
0\leq(\|x\|_X-\|v\|_X)^2\leq\langle  J_X(x)-J_X(v), x-v\rangle_X, \forall x,v\in X.
\end{equation}
Using hypotheses \textbf{(G)}, we determine that \textbf{(D)}(i) is satisfied for each $y\in Y, x\rightarrowtail \mathcal{A}_n(y,x)$. By using the facts, $\|J_X(x)\|_X=\|x\|_X$ and 
$$\langle J_X(x),x\rangle_X=\|x\|^2_X, \forall x\in X.$$ 
It is easy to show that $\mathcal{A}_n$ satisfies \textbf{(D)}(ii)–(iv). Similarly, $\mathcal{B}_n$ satisfies \textbf{(F)}. Consequently, by using Theorem \ref{t2.1}, we can argue that Problem \ref{P3.1} has a solution. 
\item[(ii)]  Let $\{(x_n,y_n)\}$ be an arbitrary sequence of solutions of Problem \ref{P3.1}. Next, a meticulous calculation yields
\begin{align}\label{eq:3.7}
r(\|x_n\|_X,\|y_n\|_Y) &\leq r(\|x_n\|_X, \|y_n\|_Y)+\dfrac{\varepsilon_n\langle J_X(x_n), x_n\rangle_X}{\|x_n\|_X}\notag\\
&\leq \dfrac{\varrho_\mathcal{A}(1+\|x_n\|_X+\|y_n\|_Y)\|x_0\|_X}{\|x_n\|_X}+\dfrac{(\varepsilon_n\|J_X(x_n)\|_{X^\ast}+\|\gamma\|_{X^\ast})\|x_0\|_X+\varrho_\varphi(x)}{\|x_n\|_X}\notag\\
&\quad\quad+ \varrho_\varphi(x) +\|\gamma\|_{X^\ast}\notag\\
&=\dfrac{\varrho_\mathcal{A}(1+\|x_n\|_X+\|y_n\|_Y)\|x_0\|_X}{\|x_n\|_X}+\dfrac{(\varepsilon_n\|x_n\|_X+\|\gamma\|_{X^\ast})\|x_0\|_X+\varrho_\varphi(x)}{\|x_n\|_X}\notag\\
&\quad\quad+\varrho_\varphi(x)+\|\gamma\|_{X^\ast},
\end{align}
and
\begin{align}\label{eq:3.8}
\ell(\|y_n\|_Y, \|x_n\|_X) &\leq \dfrac{\varrho_\mathcal{B}(1+\|x_n\|_X+\|y_n\|_Y)\|y_0\|_Y}{\|y_n\|_Y}
+\dfrac{(\delta_n\|y_n\|_Y+\|\zeta\|_{Y^\ast})\|y_0\|_Y+\varrho_\phi(y)}{\|y_n\|_Y}\notag\\
&\quad\quad +\varrho_\phi(y)+\|\zeta\|_{Y^\ast}.
\end{align}
The same arguments used in the proof of Lemma \ref{l2.3} that $\{(x_n,y_n)\}$ is bounded in $X\times Y$.\\
If necessitated we can proceed to a relabeled subsequence and presume that
\begin{equation}\label{eq:3.9}
(x_n,y_n)\xlongrightarrow[]{w}(x,y)\in X\times Y \text{~~as~~} n\to\infty, \text{~~for some~}(x,y)\in \Omega\times \mho.
\end{equation}
By using the monotonicity of $x\rightarrowtail \mathcal{A}(y,x)$ and $y\rightarrowtail \mathcal{B}(x,y)$, we can make the following deduction
\begin{equation}\label{eq:3.10}
\langle \mathcal{A}(y_n,v)+\varepsilon_n J_X(x_n),v-x_n\rangle_X+\varphi(x_n,v)-\varphi(x_n,x_n)\geq\langle \gamma, v-x_n\rangle_X, \forall v\in\Omega,
\end{equation} 
and
\begin{equation}\label{eq:3.11}
\langle \mathcal{B}(x_n,w)+\delta_n J_Y(y_n),w-y_n\rangle_Y+\phi(y_n,w)-\phi(y_n,y_n)\geq\langle \zeta, w-y_n\rangle_Y, \forall w\in \mho.
\end{equation}
Taking the upper limit as $n\to\infty$ and using hypotheses \textbf{(D)}(ii) and \textbf{(F)}(ii), we infer that 
$$\langle \mathcal{A}(y,v),v-x\rangle_X+\varphi(x,v)-\varphi(x,x)\geq\langle \gamma,v-x\rangle_X, \forall v\in \Omega,$$ 
and
$$\langle \mathcal{B}(x,w), w-y\rangle_Y+\phi(y,w)-\phi(y,y)\geq\langle \zeta,w-y\rangle_Y, \forall w\in \mho,$$
here we used the boundedness of $\{(x_n,y_n)\}\in X\times Y$. Using Minty approach, we find $(x,y)\in\Omega\times \mho$ solve Problem \ref{P1.1}, \textit{i.e.,} 
$$(x,y)\in\Gamma(\gamma,\zeta).$$
\item[(iii)] It can be deduced from assertion (ii) that if we have a sequence of solutions denoted by $\{(x_n,y_n)\}$ for Problem \ref{P3.1}, there will always exist a subsequence of $\{(x_n,y_n)\}$ which satisfies \eqref{eq:3.4}. We claim that the sequence $\{(x_n,y_n)\}$ converges strongly to $(x,y)$. It is simple to demonstrate that:
\begin{align}\label{eq:3.12}
-\alpha_\mathcal{A}\beta_\mathcal{A}\|x_n-x\|^2_X&\leq\langle \mathcal{A}(y_n,x_n)-\mathcal{A}(y_n,x),x_n-x\rangle_X\notag\\
&\leq\langle \mathcal{A}(y,x)-\mathcal{A}(y_n,x), x_n-x\rangle_X+\varepsilon_n\langle J_X(x_n), x-x_n\rangle_X\notag\\
\Longrightarrow&\notag\\
\alpha_\mathcal{A}\beta_\mathcal{A}\|x_n-x\|^2_X&\leq\langle \mathcal{A}(y,x)-\mathcal{A}(y_n,x), x-x_n\rangle_X+\varepsilon_n\langle J_X(x_n), x_n-x\rangle_X.
\end{align}
By using hypothesis \textbf{(D)}(ii) and taking the upper limit as $n$ approaches infinity on the above inequality, we obtain
\begin{align*}
0 &\leq\lim_{n\to\infty}\inf\alpha_\mathcal{A}\beta_\mathcal{A}\|x_n-x\|^2_X\\
&\leq\lim_{n\to\infty} \sup\|x_n-x\|^2_X\\
&\leq\langle \mathcal{A}(y, x)-\mathcal{A}(y_n,x), x-x_n\rangle_X+\lim_{n\to\infty} \sup \varepsilon_n\|x_n\|_X\|x_n-x\|_X\\
&\leq0.
\end{align*}
 This implies that 
 $$x_n\to x \in X \text{~~as~~} n\to\infty.$$ 
 On the other hand, it has 
 $$y_n\to y \in Y \text{~~as~~} n\to\infty.$$ 
\end{itemize}
\end{proof}
%++++++++++++++++++++++++++++++++++++++++++++++++
\section{Optimal control}
In this section, we explore optimal control for the system of nonlinear mixed variational inequality problems. Additionally, we consider an optimal control problem driven by the system of nonlinear mixed variational inequality problems and prove its solvability.\par
Consider two Banach spaces $Z_1$ and $Z_2$ with continuous embeddings from $X$ to $Z_1$ and from $Y$ to $Z_2$. Let $x_0\in Z_1$ and $y_0\in Z_2$ be two target profiles. We define subspaces $U\subset X^\ast$ and $V\subset Y^\ast$ such that the embeddings from $U$ to $X^\ast$ and $V$ to $Y^\ast$ are compact. Now, we investigate the following optimal control problem:\par
\begin{problem}\label{P4.1} Find $(\gamma^\ast,\zeta^\ast)\in U \times V$	such that
\begin{equation}\label{eq:4.1}
\daleth(\gamma^\ast,\zeta^\ast)=\inf_{(\gamma,\zeta)\in U\times V} \daleth(\gamma,\zeta),
\end{equation}
where the cost function $\daleth\colon U\times V\to \mathbb{R}$ is deﬁned by
\begin{equation}\label{eq:4.2}
\daleth(\gamma,\zeta)=\inf_{(x,y)\in\Gamma(\gamma,\zeta)}\left(\dfrac{\rho}{2}\|x-x_0\|^2_{Z_1}+\dfrac{\theta}{2}\|y-y_0\|^2_{Z_2}\right)+\varUpsilon(\gamma,\zeta).
\end{equation}
Here, $\Gamma(\gamma,\zeta)$ represents the solution set of Problem \ref{P1.1} for $(\gamma,\zeta)\in X^\ast \times Y^\ast$, with $\rho > 0$ and $\theta > 0$ as regularized parameters.
\end{problem}
We assume that the function $\varUpsilon$ satisfies the following conditions:\par
\textbf{(K):} $\varUpsilon:U\times V\to\mathbb{R}$ is such that
\begin{itemize}
\item[(i)] $\varUpsilon$ is bounded from below;
\item[(ii)] $\varUpsilon$ is coercive on $U\times V$, namely it holds
$$\lim\limits_{\substack{(\gamma,\zeta)\in U\times V\\
\|\gamma\|_U+\|\zeta\|_V\to\infty}} \varUpsilon(\gamma,\zeta)\to+\infty;$$
\item[(iii)] $\varUpsilon$ is weakly lower semicontinuous on $U\times V$, \textit{i.e.},
$$\lim_{n\to\infty} \inf\varUpsilon(\gamma_n,\zeta_n)\geq \varUpsilon(\gamma,\zeta),$$	
whenever $\{(\gamma_n,\zeta_n)\}\subset U\times V$ and $(\gamma,\zeta)\in	U\times V$ are such that
$$(\gamma_n,\zeta_n)\xlongrightarrow[]{w}(\gamma,\zeta) \in U \times V\text{~~as~~} n\to\infty.$$
\end{itemize}
In this context, we are exploring the existence result for Problem \ref{P4.1}.
\begin{theorem}\label{t4.1}  
Assume that \textbf{(A), (B), (C), (D)}(ii)–(iv), \textbf{(E)} and \textbf{(F)}(ii)–(iv) hold. If, in addition,  \textbf{(K)} and \textbf{(G)} are fulﬁlled. Then Problem \ref{P4.1} has an optimal control pair.
\end{theorem}
\begin{proof}	
For each fix $(\gamma,\zeta)\in U\times V$, the closedness of $\Gamma(\gamma,\zeta)$ (\textit{see} Theorem \ref{t2.1}) guarantees that there exists $(\hat{x}, \hat{y}) \in \Gamma(\gamma,\zeta)$ such that
\begin{equation}\label{eq:4.3}
\dfrac{\rho}{2}\|\hat{x}-x_0\|^2_{Z_1}+\dfrac{\theta}{2}\|\hat{y}-y_0\|^2_{Z_2}=\inf_{(x,y)\in \Gamma(\gamma,\zeta)}\left(\dfrac{\rho}{2}\|x -x_0\|^2_{Z_1}+\dfrac{\theta}{2}\|y-y_0\|^2_{Z_2}\right),
\end{equation}
%\textit{i.e.,}
%$$\inf_{(x,y)\in \Gamma(\gamma,\zeta)}\left(\dfrac{\rho}{2}\|x-x_0\|^2_{Z_1}+\dfrac{\theta}{2}\|y-y_0\|^2_{Z_2}\right)$$	
is attainable.\par
From the definition of $\daleth$ and hypothesis \textbf{(K)}(i), there exists a minimizing sequence $\{(\gamma_n, \zeta_n)\}\subset U\times V$ such that
\begin{equation}\label{eq:4.4}
\lim_{n\to\infty}\daleth(\gamma_n,\zeta_n)=\inf_{(\gamma,\zeta)\in U\times V} \daleth(\gamma,\zeta).
\end{equation}
We assert that the sequence $\{(\gamma_n,\zeta_n)\}$ is bounded in $U \times V$. Arguing by contradiction, we suppose that
$$\|\gamma_n\|_U + \|\zeta_n\|_V \to +\infty \text{~~as~~} n\to\infty.$$ 
The latter together with hypothesis \textbf{(K)}(ii), leads to the conclusion:
\begin{align}\label{eq:4.5}
\inf_{(\gamma,\zeta)\in U\times V}\daleth(\gamma,\zeta)&=\lim_{n\to\infty}\daleth(\gamma_n,\zeta_n) \notag\\
&\geq \lim_{n\to\infty}\varUpsilon(\gamma_n,\zeta_n)\notag\\
&=+\infty.
\end{align}
This leads to a contradiction, so, $\{(\gamma_n,\zeta_n)\}$ is bounded in $U\to V$. Passing to a relabeled subsequence if necessary, we may assume that
\begin{equation}\label{eq:4.6}
(\gamma_n,\zeta_n)\xlongrightarrow[]{w}(\gamma^\ast, \zeta^\ast) \in U\times V \text{~~as~~} n\to\infty,\text{~~for some~~} (\gamma^\ast, \zeta^\ast)\in U\times V.
\end{equation}
Let $\{(x_n, y_n)\}\subset \Omega\times \mho$ satisfy \eqref{eq:4.3} by taking $\hat{x}=x_n$, $\hat{y}=y_n$, and $(\gamma,\zeta)=(\gamma_n,\zeta_n)$. We will now prove that $\{(x_n, y_n)\}\subset \Omega\times \mho$ is uniformly bounded in $X\times Y$. A direct computation shows that
\begin{align}\label{eq:4.7}
r(\|x_n\|_X, \|y_n\|_Y)&\leq\dfrac{\varrho_\mathcal{A}(1+\|x_n\|_X+\|y_n\|_Y)\|x_0\|_X}{\|x_n\|_X}
+\dfrac{(\varrho_\varphi(x)+\|\gamma_n\|_{X^\ast})\|x_0\|_X}{\|x_n\|_X}\notag\\
&\quad\quad+\varrho_\varphi(x)+\|\gamma_n\|_{X^\ast},
\end{align}
and
\begin{align}\label{eq:4.8}
\ell(\|y_n\|_Y, \|x_n\|_X)&\leq\dfrac{\varrho_\mathcal{B}(1+\|x_n\|_X+\|y_n\|_Y)\|y_0\|_Y}{\|y_n\|_Y}
+\dfrac{(\varrho_\phi(y)+\|\zeta_n\|_{Y^\ast}\|y_0\|_Y}{\|y_n\|_Y}\notag\\
&\quad\quad +\varrho_\phi(y)+\|\zeta_n\|_{Y^\ast}.
\end{align}
Since the embeddings from $U$ to $X^\ast$ and from $V$ to $Y^\ast$ are both continuous. We can apply the same arguments as in the proof of Lemma \ref{l2.3} to infer that $\{(x_n,y_n)\}\subset\Omega \times \mho$ is uniformly bounded in $X\times Y$. Without loss of generality, we assume that
\begin{equation}\label{eq:4.9}
(x_n,y_n)\xlongrightarrow[]{w}(x^\ast, y^\ast) \in X \times Y, \text{~~and~~} Z_1\times Z_2\text{~~as~~} n\to\infty,\text{~~for some~~} (x^\ast, y^\ast)\in	\Omega\times \mho.
\end{equation} 
Using the Minty approach yields
\begin{equation}\label{eq:4.10}
\langle \mathcal{A}(y_n,v),v-x_n\rangle_X+\varphi(x_n,v)-\varphi(x_n,x_n)\geq\langle \gamma_n, v-x_n\rangle_X,  \forall v\in \Omega,
\end{equation} 
and
\begin{equation}\label{eq:4.11}
\langle \mathcal{B}(x_n,w), w-y_n\rangle_Y+\phi(y_n,w) -\phi(y_n,y_n)\geq \langle \zeta_n, w-y_n\rangle_Y, \forall w\in \mho.
\end{equation}
The compactness of the embedding from $(U, V)$ into $(X^\ast,Y^\ast)$ and \eqref{eq:4.6}, implies that
$$(\gamma_n,\zeta_n)\to(\gamma^\ast,\zeta^\ast) \in X^\ast\times Y^\ast \text{~~as~~} n\to\infty.$$ 
Taking the upper limit as $n\to\infty$ for inequalities \eqref{eq:4.10}–\eqref{eq:4.11}, we gain
$$\langle \mathcal{A}(y^\ast,v),v-x^\ast\rangle_X+\varphi(x^\ast,v)-\varphi(x^\ast,x^\ast)\geq\langle \gamma^\ast, v-x^\ast\rangle_X,  \forall v\in \Omega,$$
and
$$\langle \mathcal{B}(x^\ast, w), w-y^\ast\rangle_Y+\phi(y^\ast,w)-\phi(y^\ast,y^\ast)\geq\langle \zeta^\ast, w-y^\ast\rangle_Y, \forall w\in \mho,$$
where we used the conditions \textbf{(F)}(ii) and \textbf{(D)}(ii). Once more using the Minty trick, we get 
$$(x^\ast, y^\ast) \in\Gamma(\gamma^\ast,\zeta^\ast).$$
However, the weak lower semicontinuity of $\|\cdot\|_{Z_1}$ and $\|\cdot\|_{Z_2}$ suggests
\begin{equation}\label{eq:4.12}
\dfrac{\rho}{2}\|x^\ast-x_0\|^2_{Z_1}+\dfrac{\theta}{2}\|y^\ast-y_0\|^2_{Z_2}\leq\liminf_{n\to\infty}  \left[\dfrac{\rho}{2}\|x_n-x_0\|^2_{Z_1}+\dfrac{\theta}{2}\|y_n	-y_0\|^2_{Z_2}\right].
\end{equation}
Note that $\varUpsilon$ is weakly lower semicontinuous on $U\times V$, it implies
\begin{equation}\label{eq:4.13}
\varUpsilon(\gamma^\ast, \zeta^\ast) \leq	\liminf_{n\to\infty} \varUpsilon(\gamma_n,\zeta_n).
\end{equation}
Referring to equations \eqref{eq:4.12} and \eqref{eq:4.13}, we have
\begin{align}\label{eq:4.14}
\liminf_{n\to\infty} \daleth(\gamma_n,\zeta_n)&\geq\liminf_{n\to\infty}\inf_{(x,y)\in \Gamma(\gamma_n,\zeta_n)}\left( \dfrac{\rho}{2}\|x-x_0\|^2_{Z_1}+\dfrac{\theta}{2}\|y-y_0\|^2_{Z_2}\right)+\liminf_{n\to\infty} \varUpsilon(\gamma_n,\zeta_n)\notag\\
&=\liminf_{n\to\infty} \left(\dfrac{\rho}{2}\|x_n-x_0\|^2_X+\dfrac{\theta}{2}\|y_n-y_0\|^2_Y\right)+\liminf_{n\to\infty} \varUpsilon(\gamma_n,\zeta_n)\notag\\
&\geq\dfrac{\rho}{2}\|x^\ast-x_0\|^2_{Z_1}+\dfrac{\theta}{2}\|y^\ast-y_0\|^2_{Z_2}+ \varUpsilon(\gamma^\ast,\zeta^\ast), ~~(\text{where~}(x^\ast,y^\ast)\in\Gamma(\gamma^\ast,\zeta^\ast))\notag\\
&\geq\inf_{(x,y)\in\Gamma(\gamma^\ast,\zeta^\ast)}\left(\dfrac{\rho}{2}\|x-x_0\|^2_{Z_1}+\dfrac{\theta}{2}\|y-y_0\|^2_{Z_2}\right)+ \varUpsilon(\gamma^\ast,\zeta^\ast)\notag\\
&=\daleth(\gamma^\ast,\zeta^\ast).
\end{align}
We can use equation \eqref{eq:4.14} along with \eqref{eq:4.4} to arrive at the following conclusion:
$$\daleth(\gamma^\ast, \zeta^\ast) \leq\inf_{(\gamma,\zeta)\in U\times V} \daleth(\gamma,\zeta),$$
namely $(\gamma^\ast,\zeta^\ast)$ is an optimal control pair of Problem \ref{P4.1}.
\end{proof}

%============================================
%\section{Some special type of set-valued variational-hemivariational inequality problems}
%============
%\section{Conclusions}
%The paper introduces and examines a new type of coupled variational inequality problem in Banach spaces. The results concerning the existence and uniqueness of the solution to coupled variational inequality problems are delivered using the Kakutani-Ky Fan fixed point theorem, combined with the Minty method and the arguments of monotonicity. Additionally, a stability result for the coupled variational inequality problem is established, and an optimal control problem driven by the coupled variational inequality problem is considered. Problems of this nature are encountered in transport optimization, Nash equilibrium problems of multiple players, contact mechanics problems, and related fields. In the future, the theoretical results established in the current paper will be applied to an optimization problem, and coupled quasi-variational inequality problems will be investigated.
%=========================================================================================
\section{Declaration statements}
\noindent
\textbf{Availability of data and material:}
The data sets used and/or analysed during the current study are available from the corresponding author upon reasonable request.\\
\noindent
\textbf{Competing interests:}
No potential conflict of interest was reported by the author(s).\\
\noindent
\textbf{Funding:}
Deputyship for Research $\&$ Innovation, Ministry of Education in Saudi Arabia for funding this  research  work through Project Number: ISP-2024.\\
\noindent
\textbf{Author's contributions:}
%Conceptualization,  \textit{Ahmadini}; Writing of original draft, \textit{Salahuddin}; Methodology, \textit{Gissy}; Writing of review and editing, \textit{Salahuddin};  Validation, \textit{Gissy}.\par
\noindent
All authors read and approved the final manuscript.\\
\noindent
\textbf{Acknowledgements:}
The authors extend their appreciation to the Deputyship for Research $\&$ Innovation, Ministry of Education in Saudi Arabia for funding this  research  work through Project Number: ISP-2024.\\%The authors are grateful to the anonymous referee for their useful comments and suggestions, which have improved the quality of the paper.\\
\noindent
\textbf{Author's information:}
Not applicable.\\
%===========================================================
%\textbf{Acknowledgement.}

%The authors would like to express their heartfelt gratitude to the referees for their insightful criticism, which helped to enhance the paper's presentation.

%\noindent
%\textbf{Disclaimer/Publisher’s Note:} Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
%The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
%=====================================================================================================
\begin{thebibliography}{1}
\small
\bibitem{Li} Lin, L.J.: Systems of generalized quasi-variational inclusions problems with applications with applications to variational analysis and optimization problems. J. Glob. Optim., 38, 21-39 (2007).
\bibitem{ChSaAhWaWa} Chang, S.S., Salahuddin, Ahmadini, A.A.H., Wang, L., Wang, G.: The penalty method for generalized mixed variational-hemivariational inequality problems. Carpathian J. Math. \textbf{38}(2),  357-381 (2022).
\bibitem{ChSaWaMa} Chang, S. S., Salahuddin, Wang, L., Ma, Z.: The convergence results of differential variational inequality problems. Symmetry, MDPI. 14, 760 (2022), doi.org/10.3390/sym 14040760.
\bibitem{HaSo} Han, W., Sofonea, M.: Quasistatic contact problems in viscoelasticity and viscoplasticity. American Mathematical Society, International Press 2002).
\bibitem{HlHaNeLo}	Hlaváˇcek, I., Haslinger, J., Neˇcas, J., Lovíšek, J.: Solution of variational inequalities in mechanics, Applied Mathematical Sciences, vol. 66. Springer, New York 1988.
\bibitem{KiSaAh} Kim, J. K., Salahuddin, Ahmadini, A. A. H.: The study of error bounds for generalized vector inverse mixed quasi-variational inequalities. Filomat, 37 (2), 627-642 (2023).
\bibitem{KiSt}Kinderlehrer, D., Stampacchia, G.: An Introduction to Variational Inequalities and Their Application. Academic Press, New York 1980.
\bibitem{AlAhSa} Abd-Allah, A.A., Ahmadini, A.A.H., Salahuddin: Hybrid iterative algorithms for finding common solutions of a system of generalized mixed quasi-equilibrium problems and fixed point problems of nonexpansive semigroups. Int. J. Nonlinear Anal. Appl. \textbf{14}(1) 2271-2788 (2023).
\bibitem{DeMiPa03} Denkowski, Z.,	Migórski, S., Papageorgiou, N.S.: An	Introduction to Nonlinear Analysis: Theory. Kluwer Academic/Plenum Publishers, Boston, Dordrecht, London, New York  2003.
\bibitem{GaPa} Gasi´nski, L., Papageorgiou, N.S.: Nonlinear analysis. Series in mathematical analysis and applications. Chapman $\&$ Hall/CRC, Boca Raton, FL  2006. 
\bibitem{LiXuLi} Liu, Y., Xu, H.,	Lin,	G.H.:	Stability analysis of two-stage stochastic mathematical programs with complementarity constraints via NLP regularization. SIAM J. Optim. 21(3), 669–705 (2011).
\bibitem{Ze} Zeidler, E.: Nonlinear functional analysis and its applications:II/B: Nonlinear monotone operators. Springer, Berlin, Germany 2007.
\bibitem{LiYaZeZh} Liu, J., Yang, X, Zeng, S., Zhao, Y.: Coupled variational inequalities: Existence, Stability and Optimal Control. J. Optim. Theo. Appl., 193, 877-909 (2022).
\bibitem{Cl}Clarke, F.H.: Optimization and nonsmooth analysis. Wiley, New York 1983.
\bibitem{PaKy} Papageorgiou, N.S., Kyritsi, S.: Handbook of applied analysis. Springer, New York 2009.
\bibitem{CaSoXi} Cai, D. L., Sofonea, M, Xiao, Y. B.: Convergence results for elliptic variational-hemivariational inequalities. Adv. Nonlinear Anal. 10, 2-23 (2021).
\bibitem{LiMoZe} Liu, Z.H., Motreanu, D., and Zeng, S.D.: Nonlinear evolutionary systems are driven by mixed variational inequalities and their applications. Nonlinear Analysis, 42, 409-421 (2018).
\bibitem{LiMoZe18}Liu, Z.H., Motreanu, D., Zeng, S.D.: Nonlinear evolutionary systems driven by quasi-hemivariational inequalities. Math. Meth. Appl. Sci. 41, 1214–1229 (2018).
%\bibitem{SoMa}Sofonea, M., Matei, A.: Mathematical	models in	contact mechanics. Cambridge	University	Press, UK 2012.

\end{thebibliography}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\end{document}

