Abstract
Depending on whether a mathematical program with equilibrium constraints (MPEC) is considered in its original or its enhanced (via KKT conditions) form, the assumed qualification conditions as well as the derived necessary optimality conditions may differ significantly. In this paper, we study this issue when imposing one of the weakest possible qualification conditions, namely the calmness of the perturbation mapping associated with the respective generalized equations in both forms of the MPEC. It is well known that the calmness property allows one to derive the so-called M-stationarity conditions. The restrictiveness of assumptions and the strength of conclusions in the two forms of the MPEC is also strongly related to the qualification conditions on the “lower level”. For instance, even under the linear independence constraint qualification (LICQ) for a lower level feasible set described by \(\mathscr {C}^1\) functions, the calmness properties of the original and the enhanced perturbation mapping are drastically different. When passing to \(\mathscr {C}^{1,1}\) data, this difference still remains true under the weaker Mangasarian–Fromovitz constraint qualification, whereas under LICQ both the calmness assumption and the derived optimality conditions are fully equivalent for the original and the enhanced form of the MPEC. After clarifying these relations, we provide a compilation of practically relevant consequences of our analysis in the derivation of necessary optimality conditions. The obtained results are finally applied to MPECs with structured equilibria.
Similar content being viewed by others
References
Aussel, D., Correa, R., Marechal, M.: Electricity spot market with transmission losses. J. Ind. Manag. Optim. 9, 275–290 (2013)
Dempe, S., Dutta, J.: Is bilevel programming a special case of a mathematical program with complementarity constraints? Math. Program. A 131(1–2), 37–48 (2012)
Dontchev, A., Rockafellar, R.T.: Implicit Functions and Solution Mappings. Springer, Berlin, Heidelberg, New York (2009)
Gfrerer, H.: First order and second order characterizations of metric subregularity and calmness of constraint set mappings. SIAM J. Optim. 21(4), 1439–1474 (2011)
Gfrerer, H., Outrata, J.: On computation of limiting coderivatives of the normal-cone mapping to inequality systems and their applications. Optimization. doi:10.1080/02331934.2015.1066372
Henrion, R., Jourani, A., Outrata, J.: On the calmness of a class of multifunctions. SIAM J. Optim. 13(2), 603–618 (2002)
Henrion, R., Outrata, J., Surowiec, T.: On the co-derivative of normal cone mappings to inequality systems. Nonlinear Anal. Theory. Methods Appl. 71(3–4), 1213–1226 (2009)
Henrion, R., Outrata, J., Surowiec, T.: Analysis of M-stationary points to an EPEC modeling oligopolistic competition in an electricity spot market. ESAIM Control Optim. Calc. Var. 18(2), 295–317 (2012)
Henrion, R., Römisch, W.: On M-stationary points for a stochastic equilibrium problem under equilibrium constraints in electricity spot market modeling. Appl. Math. 52(6), 473–494 (2007)
Ioffe, A.D., Outrata, J.: On metric and calmness qualification conditions in subdifferential calculus. Set Valued Anal. 16(2–3), 199–227 (2008)
Janin, R.: Directional derivative of the marginal function in nonlinear programming. Math. Program. Study 21, 110–126 (1984)
Klatte, D., Kummer, B.: Constrained minima and Lipschitzian penalties in metric spaces. SIAM J. Optim. 13(2), 619–633 (2002)
Mordukhovich, B.S.: Variational Analysis and Generalized Differentiation I. Springer, Berlin, Heidelberg, New York (2006)
Mordukhovich, B.S., Outrata, J.: Coderivative analysis of quasi-variational inequalities with applications to stability and optimization. SIAM J. Optim. 18(2), 389–412 (2007)
Outrata, J.: A generalized mathematical program with equilibrium constraints. SIAM J. Control Optim. 38(5), 1623–1638 (2000)
Outrata, J.V.: Optimality conditions for a class of mathematical programs with equilibrium constraints. Math. Oper. Res. 24(3), 627–644 (1999)
Robinson, S.M.: Some continuity properties of polyhedral multifunctions. In: König, H., Korte, B., Ritter, K. (eds.) Mathematical Programming at Oberwolfach, Math. Programming Studies, vol. 14, pp. 206–214. Springer, Berlin (1981)
Rockafellar, R.T., Wets, R.J.B.: Variational Analysis. Springer, Berlin, Heidelberg, New York (1998)
Scholtes, S., Stöhr, M.: How stringent is the linear independence assumption for mathematical programs with complementarity constraints? Math. Oper. Res. 26(4), 851–863 (2001)
Surowiec, T.: Explicit stationarity conditions and solution characterization for equilibrium problems with equilibrium constraints. Ph.D. thesis, Humboldt University Berlin (2010)
Ye, J.J.: Constraint qualifications and necessary optimality conditions for optimization problems with variational inequality constraints. SIAM J. Optim. 10(4), 943–962 (2000)
Ye, J.J., Ye, X.Y.: Necessary optimality conditions for optimization problems with variational inequality constraints. Math. Oper. Res. 22(4), 977–997 (1997)
Acknowledgements
The authors would like to thank two anonymous referees for their critical comments which led to a substantial improvement of the paper. We are particularly indebted to Prof. Helmut Gfrerer for pointing us to an incorrect result in the first version of the paper.
Author information
Authors and Affiliations
Corresponding author
Additional information
This paper is dedicated to Prof. R. T. Rockafellar on the occasion of his 80th birthday.
This work was partially supported by the Grant Agency of the Czech Republic under grant 15-00735S (first and third author), by the Australian Research Council under grant DP-160100854 (second and third author) and by the DFG Sonderforschungsbereich/Transregio 154 (project B04, second author).
Appendix: A strong counterexample to the reversion of Proposition 2 under MFCQ and \(\mathscr {C}^2\) data for \(\varGamma \)
Appendix: A strong counterexample to the reversion of Proposition 2 under MFCQ and \(\mathscr {C}^2\) data for \(\varGamma \)
In Example 2 we have shown that under MFCQ and smooth inequalities describing the set \(\varGamma \), the mapping M may be calm, whereas the enhanced mapping \(\tilde{M}\) fails to be calm for some multiplier. In the following stronger counterexample we construct a set \(\varGamma \) described by \(\mathscr {C}^2\) inequalities satisfying MFCQ at given \(\bar{y}\) and a function F such that M is calm at \((0,\bar{x},\bar{y})\) while \(\tilde{M}\) is not calm at \((0,0,\bar{x},\bar{y},\lambda )\) for any \(\lambda \in \varLambda (\bar{x},\bar{y})\).
Define first \(\varphi _1,\varphi _2:[-1,1]\rightarrow {\mathbb {R}}\) and \(q_1,q_2:[-1,1]\times {\mathbb {R}}\rightarrow {\mathbb {R}}\) as
put \(\omega ={\mathbb {R}}\) and as the reference point take \((\bar{x},\bar{y}_1,\bar{y}_2)=(0,0,0)\). These functions are depicted in Fig. 1. Note first that MFCQ is indeed satisfied for \(\varGamma \) and that \(\varphi _1\) and \(\varphi _2\) are twice continuously differentiable. Define further
Note that for any given k only one of functions \(\varphi _1\) and \(\varphi _2\) will be active in the definition of \(\phi \) on interval \((\frac{1}{k+1},\frac{1}{k})\). Because \(\phi '(\frac{1}{k})=\phi ''(\frac{1}{k})=0\) for all \(k\in {\mathbb {N}}\), it remains to verify the twice continuous differentiability of \(\phi \) at 0. But we have
which implies that \(|\phi '(0)|=0\). Similarly we obtain \(\phi ''(0)=0\) and that \(\phi \) is twice continuously differentiable. Finally, we define \(F(x,y):=(-\phi '(y_1),\ 1)\). By construction of \(\phi \), we obtain that F is continuously differentiable. Since \(\varGamma ={\text {epi}}\phi \) we have that
As \(M(p)\subset M(0)\) for all p small enough, we obtain that M is calm at \((0,\bar{x},\bar{y})\).
It is easy to see that \(\varLambda (\bar{x},\bar{y})=\{\lambda \ge 0|\,\lambda _1+\lambda _2=1\}\). We will show now that \(\tilde{M}\) is not calm at \((0,0,\bar{x},\bar{y},\lambda )\) for any \(\lambda \in \varLambda (\bar{x},\bar{y})\). Define
and note that for all \(t\in \varOmega _2\cup \varOmega _3\) small enough it holds that \(|\varphi _2(t)|<|\varphi _1(t)|\) and for all \(t\in \varOmega _3\) small enough we have \(|\varphi _2'(t)|<|\varphi _1'(t)|\).
We will show first that \(\hat{T}_{\{1\}}\) defined in (20) is not calm at \((0,\bar{y})\). From the definition we see that
and thus
Now pick any sequence \(y_{k1}> 0\), \(y_{k1}\rightarrow 0\) such that \(y_{k1}\in \varOmega _2\) and \(\varphi _1(y_{k1})<0\) and define \(p_{k1}:=0\), \(y_{k2}:=\varphi _1(y_{k1})\) and \(p_{k2}:=\varphi _2(y_{k1})-y_{k2}\). Then \(y_k\in \hat{T}_{\{1\}}(p_k)\). Moreover, as \(\varphi _1\) and \(\varphi _2\) have the same signs
Consider now a point \(\tilde{y}_{k1}\in \varOmega _1\) at which \(d(y_{k1},\varOmega _1)\) is realized. Since \(\varOmega _1\subset \hat{T}_{\{1\}}(0)\) and \(\varphi _1\) is zero on \(\varOmega _1\), we obtain
where in the last equality we have used the mean value theorem to find some \(\xi _k\) which lies in the line segment connecting \(y_{k1}\) and \(\tilde{y}_{k1}\). Since \(\varphi _1\) is twice continuously differentiable with \(\varphi _1'(0)=0\), we have proved that \(\hat{T}_{\{1\}}\) is not calm at \((0,\bar{y})\). For \(\hat{T}_{\{2\}}\) we proceed with a similar construction. In this case we have
and for the contradicting sequence we choose some \(y_{k1}>0\), \(y_{k1}\rightarrow 0\) such that \(y_{k1}\in \varOmega _2\) and \(\varphi _1(y_{k1})>0\) and define again \(p_{k1}:=0\), \(y_{k2}:=\varphi _1(y_{k1})\) and \(p_{k2}:=\varphi _2(y_{k1})-y_{k2}\) and perform the estimates as in the previous case. Since for \(\hat{T}_{\{1,2\}}\) we have
either of the previous contradicting sequences can be chosen.
Fix now any \(\bar{\lambda }\in \varLambda (\bar{x},\bar{y})\) and consider the corresponding index set \(I=\{i|\ \bar{\lambda }_i>0\}\). In the previous several paragraphs we have shown that \(\hat{T}_I\) is not calm at \((0,\bar{y})\) and found a sequence \((\tilde{p}_k,\tilde{y}_k)\) violating the calmness property. By virtue of Lemma 4 we obtain that T is not calm at \((0,\bar{y},\bar{\lambda })\). Moreover, from the proof of this lemma we see that the sequence \((p_k,y_k,\lambda _k)\), which violates the calmness of T at \((0,\bar{y},\bar{\lambda })\), can be taken in such a way that \(p_k=\tilde{p}_k\), \(y_k=\tilde{y}_k\) and \(\lambda _k=\bar{\lambda }\) with \((\tilde{y}_k,\bar{\lambda })\in T(\tilde{p}_k)\) and
Furthermore, in all the previous cases we have chosen \(\tilde{y}_k\) in such a way that \(\tilde{y}_{k1}\in \varOmega _2\).
We will show that \(\tilde{M}\) is not calm at \((0,0,\bar{x},\bar{y},\bar{\lambda })\). Consider sequence
and show first that \((\bar{x},\tilde{y}_{k1},\tilde{y}_{k2},\bar{\lambda }_1,\bar{\lambda }_2)\in \tilde{M}(0,0,\tilde{p}_{k1},\tilde{p}_{k2})\), which amounts to showing
We know that \((\tilde{y}_k,\bar{\lambda })\in T(\tilde{p}_k)\) and hence the inclusion is satisfied. Moreover, as \(\tilde{y}_{k1}\in \varOmega _2\) by construction of this sequence and as \(\bar{\lambda }_1+\bar{\lambda }_2=1\), we indeed obtain
From the respective definitions of \(\tilde{M}\) and T, we infer that \(\tilde{M}(0,0,0,0)\subset {\mathbb {R}}^n\times T(0,0)\) and consequently due to (55) we obtain
This together with (56) and (57) implies that \(\tilde{M}\) is indeed not calm at \((0,0,\bar{x},\bar{y},\bar{\lambda })\). Since \(\bar{\lambda }\) was chosen arbitrarily from \(\varLambda (\bar{x},\bar{y})\), the construction has been completed.
Rights and permissions
About this article
Cite this article
Adam, L., Henrion, R. & Outrata, J. On M-stationarity conditions in MPECs and the associated qualification conditions. Math. Program. 168, 229–259 (2018). https://doi.org/10.1007/s10107-017-1146-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10107-017-1146-3
Keywords
- Mathematical programs with equilibrium constraints
- Optimality conditions
- Constraint qualification
- Calmness
- Perturbation mapping