Global convergence of a modified BFGS-type method based on function information for nonconvex multiobjective optimization problems

Yingxue Yang ∗

Xin Deng †

Liping Tang ‡

October 14, 2024

Abstract: In this paper, based on function information, we propose a modified BFGS-type method for nonconvex multiobjective optimization problems (MFQNMO). In the multiobjective quasi-Newton method (QNMO), each iteration involves separately approximating the Hessian matrix for each component objective function, which results in significant storage and computational burdens. MFQNMO employs a common BFGS-type matrix to approximate the Hessian matrix of all objective functions in each iteration. This matrix is updated using function information from the previous step. This approach strikes a balance between efficacy and computational cost. We confirm the convergence of the method without relying on convexity assumptions, under mild conditions, we establish a local superlinear convergence rate for MFQNMO. Furthermore, we validate its effectiveness through experiments on both nonconvex and convex test problems.

Keywords: Nonconvex multiobjective optimization, Quasi-Newton methods, Wolfe line search, Global convergence, Superlinear convergence.

AMS subject classifications:

## 1 Introduction

In this paper, we consider the following unconstrained multiobjective optimization problem (MOP):

$$\min _ { x \in \mathbb { R } ^ { n } } f ( x ) = \left ( f _ { 1 } \left ( x \right ), f _ { 2 } \left ( x \right ), \dots, f _ { m } \left ( x \right ) \right ) ^ { \top },$$

where f i ( x ) : R n → R is continuously differentiable function, for any i ∈ [ m ]. Multiobjective optimization is widely used in many research areas, which can be found in engineering [18, 23], environmental analysis [1,13], machine learning [28,32], etc.

In 2000, Fliege and Svaiter [7] proposed the steepest descent method for multiobjective optimization (SDMO). Building on their pioneering work, researchers have conducted extensive

∗ National Center for Applied Mathematics in Chongqing, Chongqing Normal University, Chongqing, 401331, China, E-mails: 2022110518022@stu.cqnu.edu.cn

‡ National Center for Applied Mathematics in Chongqing, Chongqing Normal University, Chongqing, 401331, China, E-mails: tanglipings@163.com

† National Center for Applied Mathematics in Chongqing, Chongqing Normal University, Chongqing, 401331, China, E-mails: 2022110518004@stu.cqnu.edu.cn

studies on the first-order descent algorithm for solving (MOP), (see e.g. [2,15,19,31]). The firstorder algorithm is often favored for its lower cost than the second-order method, as it avoids the computational burden of determining the Hessian matrix. However, this advantage comes with a trade-off: first-order methods typically exhibit slower convergence rates, especially ill-conditioned problems.

Ill-conditioned problems exhibit large curvature variations at certain points, resulting in highly unstable curvature values. The Newton method [8] is effective in such cases as it captures the local curvature of functions and adjusts the direction more efficiently, leading to a faster rate compared to first-order methods. Fliege et al. [8] proved the superlinear convergence rate of the Newton method. However, despite its efficiency, the computational cost per iteration is significantly higher, particularly for large-scale problems.

To reduce the computational costs incurred by the direct calculation of the Hessian matrix, the quasi-Newton method is considered. Unlike single objective optimization, for individual objective functions, the quasi-Newton methods [24, 27] employ different matrices to approximate their Hessian matrices, capturing the local curvature. Numerous variants of the quasi-Newton method have been proposed for multi-objective problems (MOPs). In [20], the global convergence of the quasi-Newton method is established under average nonmonotone line search. Morovati et al. [21] extended both the self-scaling BFGS (SS-BFGS) method and the Huang BFGS (H-BGFS) method for multiobjective optimization problems. Prudente and Souza [26] proposed a modified BFGS update formula under nonconvex assumptions, ensuring that the algorithm is well-defined. Furthermore, under the assumption of strong convexity, they establish the global convergence and local superlinear convergence rate of the method with Wolfe line search.

In large-scale multiobjective optimization problems, an increase in the number of objective functions results in an exponential growth in the scale of the matrices that need to be stored and computed, particularly in cases with high-dimensional variable spaces. Consequently, it is impractical to calculate or approximate the Hessian matrix for each objective function separately in large-scale scenarios. To address this, Ansary and Panda [3] utilized a common matrix to approximate all Hessian matrices in each iteration and establish the global convergence under an aggregated Wolf line search. Lapucci and Mansueto [14] extended the limited memory QuasiNewton approach (L-BFGS) for large-scale MOPs, and the global convergence with R-linear convergence rate is obtained under the strong convexity assumption. However, all the improved algorithms mentioned above require convexity conditions to ensure global convergence. Chen et al. [5] proposed a variable metric method for unconstrained MOPs. Based on an assumption weaker than the convexity condition, this method establishes strong convergence. Additionally, they mentioned that solving the quadratic constrained subproblem of the second-order method is more difficult than solving the direction-finding dual problem of the first-order method. In [25], the global convergence can be guaranteed with Wolf line search in the nonconvex case.

For MOPs, there are relatively few studies on quasi-Newton methods in nonconvex cases. This raises the question of whether alternative quasi-Newton method exists that can strike a balance between global convergence and per-iteration cost for nonconvex MOPs. In this paper, inspired by [3,25], by utilizing function information, we construct a common modified BFGS-type matrix to approximate the Hessian matrices of all objective functions in each iteration. We proposed a modified BFGS-type method based on function information for nonconvex MOPs (MFQNMO). We verify that the modified BGFS-type matrix is positive definite for each iteration and establish a global convergence for MFQNMO without convexity assumption. Under some mild conditions, we state the superlinear convergence property of MFQNMO. It also shows superior performance in numerical experiments. Our work has the following advantages:

- · Rather than approximating the Hessian matrix for each objective function individually, our method utilizes a single common matrix to capture the curvature information. This approach strategically balances the effectiveness of the optimization algorithm with its computational cost.
- · For nonconvex functions, most versions of quasi-Newton methods encounter difficulties in achieving convergence. Therefore, the MFQNMO method proposed in this paper for nonconvex unconstrained problems is particularly significant and valuable.
- · our method's update process for the BFGS-type formula does not rely on any specific parameter selections. Instead, it only utilizes information from prior iterations.

The structure of this paper is as follows: In Sect. 2, we give some notations and definitions for our later use. Sect. 3 recalls several multiobjective descent methods, including SDMO, QNMO and MQNMO. In Sect. 4, we introduce the proposed MFQNMO method and establish its convergence analysis. Numerical results are presented in Sect. 5. At the end of the paper, we draw some conclusions.

## 2 Preliminaries

In this section, we recall some notations, definitions and results used in the sequel. Let R n be the n -dimensional Euclidean space, and R n + denote nonnegative orthant of R n . N represents the set of natural numbers. For any x := ( x , . . . , x 1 n ) and y := ( y , . . . , y 1 n ) ∈ R n , we define:

$$n$$

̸

$$x & \geq y \quad \Leftrightarrow \quad x - y \in \mathbb { R } _ { + } ^ { n } \, \wedge \, \{ 0 \} ; \\ x & \gtrsim y \quad \Leftrightarrow \quad x - y \notin \mathbb { R } _ { + } ^ { n } \, \wedge \, \{ 0 \} ; \\ x & > y \quad \Leftrightarrow \quad x - y \in \text{int} \mathbb { R } _ { + } ^ { n } ; \\ x & \gtrsim y \quad \Leftrightarrow \quad x - y \notin \text{int} \mathbb { R } _ { + } ^ { n }.$$

$$+$$

The following are some notations that will be used.

- · [ m ] = 1 2 { , , · · · , m } .

$$\bullet \Delta _ { m } = \left \{ \lambda \, \colon \sum _ { i \in [ m ] } \lambda _ { i } = 1, \lambda _ { i } \geqslant 0, \ i \in [ m ] \right \}.$$

- · ∥·∥ is the Euclidean norm in R n .
- · Jf ( x ) ∈ R m n × , ∇ f i ( x ) ∈ R n , ∇ 2 f i ( x ) the Jacobian matrix, the gradient of f i and the Hessian matrix at x , respectively.
- · D ( x, d ) := max i ∈ [ m ] ∇ f i ( x ) ⊤ d , d ∈ R n .
- · For all B ∈ R n × n , B ≻ 0 means that B is positive definite, B ⪰ 0 means that B is positive semi-definite.
- · For x, y ∈ R n , ⟨ x, y ⟩ B := x ⊤ By and ∥ x ∥ B := √ ⟨ x, x ⟩ B .
- · For x ∈ R , ⌈ x ⌉ is the least integer greater than or equal to x.
- · A ⪯ B ⇔ B -A ⪰ 0 , A,B ∈ R n × n

̸

Definition 2.1. [7] A point x ∗ ∈ R n is called the Pareto (weakly Pareto ) solution for MOP, if there exists no point x ∈ R n such that f ( x ) ⩽ f ( x ∗ ) and f ( x ) = f ( x ∗ ) ( f ( x ) &lt; f ( x ∗ ) ) .

Definition 2.2. [7] A point x ∗ ∈ R n is called a Pareto critical point for MOP, if exists i ∈ [ m ] such that 〈 ∇ f i ( x ∗ ) , d 〉 ⩾ 0 , for all d ∈ R n .

Definition 2.3. [7] A vector d ∈ R n is called descent direction for f at x ∗ , if 〈 ∇ f i ( x ∗ ) , d 〉 &lt; , 0 for all i ∈ [ m ].

Lemma 2.1. [8, Theorem 3.1] The following statements hold:

- (i) if x ∗ is local weak Pareto optimal, then x ∗ is a Pareto critical point for f ;
- (ii) if f is convex and x ∗ is Pareto critical for f , then x ∗ is weak Pareto optimal;
- (iii) if f is strictly convex and x ∗ is Pareto critical for f , then x ∗ is Pareto optimal.

## 3 Descent methods for MOPs.

In this section, we revisit the steepest descent method [7] and the quasi-Newton method [24] for MOPs.

## 3.1 The steepest descent method for MOPs.

For a given x ∈ R n , consider the following subproblem (3.1):

$$\min _ { d \in \mathbb { R } ^ { n } } \max _ { i \in [ m ] } \left \langle \nabla f _ { i } \left ( x \right ), d \right \rangle + \frac { 1 } { 2 } \| d \| ^ { 2 }.$$

Since 〈 ∇ f i ( x , d ) 〉 + 1 2 ∥ d ∥ 2 is strongly convex for i ∈ [ m ], then subproblem (3.1) has a unique minimizer. We denote that θ ( x ) and d ( x ) are the optimal value and the optimal solution in problem (3.1), respectively. Therefore,

$$d \left ( x \right ) = \underset { d \in \mathbb { R } ^ { n } } { \arg \min } \max _ { i \in [ m ] } \left \langle \nabla f _ { i } \left ( x \right ), d \right \rangle + \frac { 1 } { 2 } \| d \| ^ { 2 },$$

and

$$\theta \left ( x \right ) = \min _ { d \in \mathbb { R } ^ { n } } \max _ { i \in [ m ] } \left \langle \nabla f _ { i } \left ( x \right ), d \right \rangle + \frac { 1 } { 2 } \| d \| ^ { 2 }.$$

The problem can be transformed into the following equivalent optimization problem (3.2), which is easier to solve.

$$\min _ { ( t, d ) \in \mathbb { R } \times \mathbb { R } ^ { n } } \quad t + \frac { 1 } { 2 } \| d \| ^ { 2 } \\ \text{s.t.} \quad \langle \nabla f _ { i } ( x ), d \rangle \leqslant t, \quad i \in [ m ].$$

We consider that the dual problem of problem (3.2):

$$\max _ { \lambda \in \Delta _ { m } } \min _ { ( t, d ) \in \mathbb { R } \times \mathbb { R } ^ { n } } L \left ( \lambda, ( t, d ) \right ),$$

where L λ, ( ( t, d ) ) is a Lagrange function,

$$L \left ( \lambda, ( t, d ) \right ) = t + \frac { 1 } { 2 } \| d \| ^ { 2 } + \sum _ { i \in [ m ] } \lambda _ { i } \left ( \left \langle \nabla f _ { i } \left ( x \right ), d \right \rangle - t \right ).$$

According to the KKT condition, we can obtain the steepest descent direction:

$$d _ { S D } \left ( x \right ) = - \sum _ { i \in \left [ m \right ] } \lambda _ { i } ^ { S D } \left ( x \right ) \nabla f _ { i } \left ( x \right ),$$

where λ SD ( x ) = ( λ SD 1 ( x , λ ) SD 2 ( x , ) · · · , λ SD m ( x ) ) ⊤ is the solution of the following dual problem

$$- \min _ { \lambda \in \Delta _ { m } } \frac { 1 } { 2 } \left \| \sum _ { i \in [ m ] } \lambda _ { i } \nabla f _ { i } \left ( x \right ) \right \| ^ { 2 },$$

owing to strong duality, we have:

$$\theta _ { S D } \left ( x \right ) = - \frac { 1 } { 2 } \left \| \sum _ { i \in \left [ m \right ] } \lambda _ { i } ^ { S D } \left ( x \right ) \nabla f _ { i } \left ( x \right ) \right \| ^ { 2 } = - \frac { 1 } { 2 } \left \| d _ { S D } \left ( x \right ) \right \| ^ { 2 } \leqslant 0,$$

$$\left \langle \nabla f _ { i } \left ( x \right ), d _ { S D } \left ( x \right ) \right \rangle \leqslant - \left \| d _ { S D } \left ( x \right ) \right \| ^ { 2 } \leqslant 0.$$

The following lemma describes the proposition of θ SD ( x ) and d SD ( x . )

Lemma 3.1. [7, Proposition 3.1] For problem (3.1) , θ SD ( x ) and d SD ( x ) are the optimal value and solution, respectively. Then,

- ( a ) It can be conclude that the equivalent relation is as follows.
- ( ) i The point x ∈ R n is not Pareto critical;

̸

- ( ii ) d SD ( x ) = 0 ;
- ( iii ) θ SD ( x ) &lt; 0 .
- ( b ) The function d SD ( x ) and θ SD ( x ) are continuous, for all x ∈ R n .

## 3.2 The quasi-Newton methods for MOPs (QNMO).

Quasi-Newton direction d q ( x ) in [24] is the optimal solution to the following subproblem:

$$\min _ { d \in \mathbb { R } ^ { n } } \max _ { i \in [ m ] } \nabla f _ { i } ( x ) ^ { \top } d + \frac { 1 } { 2 } d ^ { \top } B _ { i } ( x ) d,$$

where B x i ( ) ∈ R n × n , approximates the Hessian matrix ∇ 2 f i ( x . ) For i ∈ [ m B x ], i ( ) ≻ 0 such that ∇ f i ( x ) ⊤ d + 1 2 d ⊤ B x d i ( ) is strongly convex. Then, the subproblem 3.4 has a unique minimizer. We denote that θ q ( x ) and d q ( x ) be the optimal value and solution in problem (3.4), respectively,

$$d _ { q } ( x ) = \underset { d \in \mathbb { R } ^ { n } } { \arg \min } \max _ { i \in [ m ] } \nabla f _ { i } ( x ) ^ { \top } d + \frac { 1 } { 2 } d ^ { \top } B _ { i } ( x ) d,$$

and

$$\theta _ { q } ( x ) = \min _ { d \in \mathbb { R } ^ { n } } \max _ { i \in [ m ] } \nabla f _ { i } ( x ) ^ { \top } d + \frac { 1 } { 2 } d ^ { \top } B _ { i } ( x ) d.$$

- (3.4) is equivalent to the convex quadratic optimization problem:

$$\min _ { \substack { ( t, d ) \in \mathbb { R } \times \mathbb { R } ^ { n } \\ s. t. } } \quad t \\ \nabla f _ { i } ( x ) ^ { \top } d + \frac { 1 } { 2 } d ^ { \top } B _ { i } ( x ) d \leqslant t, \quad \forall i \in [ m ].$$

According to the KKT condition, we can obtain the quasi-Newton direction:

$$d _ { q } ( x ) = - \left [ \sum _ { i = 1 } ^ { m } \lambda _ { i } ^ { q } ( x ) B _ { i } ( x ) \right ] ^ { - 1 } \left ( \sum _ { i = 1 } ^ { m } \lambda _ { i } ^ { q } ( x ) \nabla f _ { i } ( x ) \right ),$$

where λ q ( x ) = ( λ q 1 ( x , λ ) q 2 ( x , ) · · · , λ q m ( x ) ) ⊤ is the solution of the following dual problem

$$- \min _ { \lambda \in \Delta _ { m } } \frac { 1 } { 2 } \left \| \sum _ { i \in [ m ] } \lambda _ { i } \nabla f _ { i } \left ( x \right ) \right \| _ { \left [ \sum _ { i \in [ m ] } \lambda _ { i } B _ { i } ( x ) \right ] ^ { - 1 } } ^ { 2 }.$$

owing to strong duality, we have

$$\theta _ { q } ( x ) = - \frac { 1 } { 2 } d _ { q } ( x ) ^ { \top } \left [ \sum _ { i = 1 } ^ { m } \lambda _ { i } ^ { q } ( x ) B _ { i } ( x ) \right ] d _ { q } ( x ).$$

## 3.3 The modified quasi-Newton methods for MOPs (MQNMO).

Even if quasi-Newton avoids calculating the second-order Hessian matrix, the per-iteration cost is still expensive. Recently, Ansary and Panda [3] proposed a modified quasi-Newton method (MQNMO) that utilized a single positive matrix to approximate all the Hessian matrices. The Modified quasi-Newton direction d x ( ) is the optimal solution to the following subproblem:

$$\min _ { d \in \mathbb { R } ^ { n } } \max _ { i \in [ m ] } \nabla f _ { i } ( x ) ^ { \top } d + \frac { 1 } { 2 } d ^ { \top } B ( x ) d,$$

where B x ( ) ∈ R n × n and B x ( ) ≻ 0. Problem (3.6) is equivalent to the convex quadratic optimization problem as follows

$$\min _ { \substack { ( t, d ) \in \mathbb { R } \times \mathbb { R } ^ { n } \\ s. t. } } \ t \\ \nabla f _ { i } ( x ) ^ { \top } d + \frac { 1 } { 2 } d ^ { \top } B ( x ) d \leqslant t, \quad \forall i \in [ m ]. \\ \nabla f _ { i } ( x ) ^ { \top } d + \frac { 1 } { 2 } d ^ { \top } B ( x ) d \leqslant t, \quad \forall i \in [ m ].$$

Similar to the analysis in the previous section, we can obtain the modified quasi-Newton direction using the KKT condition

$$d ( x ) = - B ( x ) ^ { - 1 } \left ( \sum _ { i = 1 } ^ { m } \lambda _ { i } ( x ) \nabla f _ { i } ( x ) \right ),$$

and

$$\theta ( x ) = - \frac { 1 } { 2 } d ( x ) ^ { \top } B ( x ) d ( x ),$$

where λ x ( ) = ( λ 1 ( x , λ ) 2 ( x , ) · · · , λ m ( x ) ) ⊤ is the solution of the following dual problem (3.9)

$$- \min _ { \lambda \in \Delta _ { m } } \frac { 1 } { 2 } \left \| \sum _ { i \in [ m ] } \lambda _ { i } \nabla f _ { i } \left ( x \right ) \right \| _ { B \left ( x \right ) ^ { - 1 } } ^ { 2 }.$$

According to the above analysis, we can get

$$\mathcal { D } ( x _ { k }, d _ { k } ) \leqslant - \| d _ { k } \| _ { B _ { k } } ^ { 2 }$$

The following rank one DFP formula of generating B k +1 is defined in [3] by

$$B _ { k + 1 } = ( I - \frac { y _ { k } s _ { k } ^ { \top } } { s _ { k } ^ { \top } y _ { k } } ) B _ { k } ( I - \frac { s _ { k } y _ { k } ^ { \top } } { s _ { k } ^ { \top } y _ { k } } ) + \frac { y _ { k } y _ { k } ^ { \top } } { s _ { k } ^ { \top } y _ { k } },$$

where s k = x k +1 -x k and y k = m ∑ i =1 λ k i ( ∇ f i ( x k +1 ) -∇ f i ( x k )) .

## 4 MFQNMO:Modified BFGS-type method based on function information for MOPs.

In this section, we propose a modified BFGS-type method based on function information for MOPs and analyze its properties. The formal description of the algorithm is as follows.

## Algorithm 1 MFQNMO

Require: Choose starting point x 0 ∈ R n . Let 0 &lt; σ 1 &lt; σ 2 &lt; 1, k = 0, ε ⩾ 0 and B 0 ≻ 0 .

Step 1. Compute the solution λ k = ( λ , λ k 1 k 2 , · · · , λ k m ) ⊤ of problem (3.9) and d k , θ k as (3.7), (3.8).

Step 2. If | θ k | &lt; ε , return x . k

Step 3. Compute step size α k &gt; 0 satisfies

$$f _ { i } ( x _ { k } + \alpha _ { k } d _ { k } ) \leqslant f _ { i } ( x _ { k } ) + \sigma _ { 1 } \alpha _ { k } \mathcal { D } ( x _ { k }, d _ { k } ), \ \forall i \in [ m ],$$

$$\mathcal { D } ( x ^ { k } + \alpha _ { k } d _ { k }, d _ { k } ) \geqslant \sigma _ { 2 } \mathcal { D } ( x _ { k }, d _ { k } ),$$

Step 4. Set x k +1 = x k + α d k k .

Step 5. Compute

$$\eta _ { k } = \frac { y _ { k } ^ { \top } s _ { k } } { \| s _ { k } \| ^ { 2 } },$$

$$m _ { k } = \max \{ - \eta _ { k }, 0 \} + \sum _ { i = 1 } ^ { m } \lambda _ { i } ^ { k } ( f _ { i } ( x _ { k } ) - f _ { i } ( x _ { k + 1 } ) ).$$

$$\gamma _ { k } = y _ { k } + m _ { k } s _ { k }.$$

k

i

m

+1

)

-∇

f

i

(

x

k

) and

y

k

=

∑

k

k

i

=1

i

i

Step 6. Update the BFGS-type matrices,

$$B _ { k + 1 } = B _ { k } - \frac { B _ { k } s _ { k } s _ { k } ^ { \top } B _ { k } } { s _ { k } ^ { \top } B _ { k } s _ { k } } + \frac { \gamma _ { k } \gamma _ { k } ^ { \top } } { \gamma _ { k } ^ { \top } s _ { k } },$$

$$B _ { k + 1 } ^ { - 1 } = ( I - \frac { s _ { k } \gamma _ { k } ^ { \top } } { s _ { k } ^ { \top } \gamma _ { k } } ) B _ { k } ^ { - 1 } ( I - \frac { \gamma _ { k } s _ { k } ^ { \top } } { s _ { k } ^ { \top } \gamma _ { k } } ) + \frac { s _ { k } s _ { k } ^ { \top } } { s _ { k } ^ { \top } \gamma _ { k } }.$$

Step 7. Set k ← k +1 and go to Step 1.

Remark 4.1. In step 1, solving the descent direction through the dual problem of subproblem 3.9, it is noted that the initial BFGS matrix is selected as a positive definite matrix, usually chosen as the unit matrix. Step 2 is the termination criterion of the algorithm. In step 3, we search the step size through the Wolfe line search. In steps 5-6, we notice that the BFGS-type update formula (4.6) and (4.7) is different from [25]. We use the previous iterative function information to update the m k adaptively. It can note that

$$B _ { k + 1 } s _ { k } & = y _ { k } + m _ { k } s _ { k }.$$

The following result shows that Algorithm 1 is well-defined without convexity assumption.

Theorem 4.1. x k is not the Pareto critical point, then B k updated by Algorithm 1 is positive definite, for all k ≥ 0 .

λ µ

.

where

s

k

=

x

k

+1

-

x

k

,

µ

=

∇

f

i

(

x

k

Proof. By the definition of γ k , η k and m k , we have

$$\text{of. By the definition of } \gamma _ { k }, \, \eta _ { k } \text{ and } m _ { k }, \, \text{we have} \\ \gamma _ { k } ^ { \top } \, s _ { k } & = y _ { k } ^ { \top } s _ { k } + m _ { k } \| s _ { k } \| ^ { 2 } \\ & = \left ( \frac { y _ { k } ^ { \top } } { \| s _ { k } \| ^ { 2 } } + m _ { k } \right ) \| s _ { k } \| ^ { 2 } \\ & = ( \eta _ { k } + m _ { k } ) \| s _ { k } \| ^ { 2 } \\ & = \left [ \eta _ { k } + \max \{ - \eta _ { k }, 0 \} + \sum _ { i = 1 } ^ { m } \lambda _ { i } ^ { k } \left ( f _ { i } ( x _ { k } ) - f _ { i } ( x _ { k + 1 } ) \right ) \right ] \| s _ { k } \| ^ { 2 } \\ & \geqslant \sum _ { i = 1 } ^ { m } \lambda _ { i } ^ { k } \left ( f _ { i } ( x _ { k } ) - f _ { i } ( x _ { k + 1 } ) \right ) \| s _ { k } \| ^ { 2 } \\ & > 0. \\ \text{en for all vector } x \in \mathbb { R } ^ { n },$$

Then for all vector x ∈ R n ,

$$x ^ { \top } B _ { k + 1 } x = \| x \| _ { B _ { k } } ^ { 2 } - \frac { \langle x, s _ { k } \rangle _ { B _ { k } } ^ { 2 } } { \| s _ { k } \| _ { B _ { k } } ^ { 2 } } + \frac { ( x ^ { \top } \gamma _ { k } ) ^ { 2 } } { \gamma _ { k } ^ { \top } s _ { k } },$$

it follows by the Cauchy-Schwarz inequality and (4.10) that

$$x ^ { \top } B _ { k + 1 } x \geqslant \frac { ( x ^ { \top } \gamma _ { k } ) ^ { 2 } } { \gamma _ { k } ^ { \top } s _ { k } } \geqslant 0.$$

Since x ⊤ B k +1 x = 0 if and only if x ⊤ γ k = 0 and x = cs k , where c is a constant. It follows that cs ⊤ k γ k = 0. By (4.10), we can obtain that c = 0 and x = 0 . Then for all nonzero vector x ∈ R n ,

$$x ^ { \top } B _ { k + 1 } x > 0,$$

we conclude the proof.

Assumption 4.1. For any z ∈ R n , the level set { x : f ( x ) ⩽ f ( z ) } is a bounded set.

Next, we show that the Wolf line search in Algorithm 1 can terminate at a finite step.

Proposition 4.1. Suppose Assumption 4.1 holds and d k is the descent direction. Then there exists an interval of values [ α , α 1 2 ] , with 0 &lt; α 1 &lt; α 2 , such that for all α ∈ [ α , α 1 2 ] , inequalities (4.1) and (4.2) hold.

Proof. The proof is similar to that in [14, Proposition 2], we omit it here.

Lemma 4.1. Let θ x ( ) and d x ( ) be the optimal value and solution in problem (3.6) . Then we have:

̸

- ( a ) Vector d x ( ) = 0 is a descent direction.
- ( b ) The following conditions are equivalent:

̸

- ( ) i The point x ∈ R n is not Pareto critical; ( ii ) d x ( ) = 0 ; ( iii ) θ x ( ) &lt; 0 .

Proof. These assertions can be derived by using the same proof demonstrated in the [25, Lemma 2.2].

## 4.1 Global convergence

In this section, we present the global convergence analysis of Algorithm 1 without assuming convexity. First, we require the following assumptions.

Assumption 4.2. (i) For i ∈ [ m ], ∇ f i ( x ) is Lipschitz continuous with constant L . i (ii)For any x ∈ R n , there exists constant a and b such that

$$a I _ { n } \preceq B ( x ) \preceq b I _ { n }.$$

In the subsection, the following properties imply that under the Wolfe line search, the objective function will be significantly decreased along the direction generated by Algorithm 1.

Lemma 4.2. Let B k be generated by the BFGS update formula (4.6) in Algorithm 1. If there exists constants c 1 and c 2 such that γ ⊤ k s k ∥ s k ∥ 2 ⩾ c 1 &gt; 0 , ∥ γ k ∥ 2 γ ⊤ k s k ⩽ c 2 . Then, for any p ∈ (0 1) , , there exists a constant δ &gt; 0 such that, for any k ∈ N ,

$$\cos \beta _ { j } = \frac { s _ { j } ^ { \top } B _ { j } s _ { j } } { \| s _ { j } \| \| B _ { j } s _ { j } \| } \geqslant \delta,$$

holds for at least ⌈ p k ( +1) ⌉ values of j ∈ [0 , k . ]

Proof.

The proof is similar to that in [4, Theorem 2.1], we omit it here.

Proposition 4.2. Let { x k } be the sequence generated by Algorithm 1. Suppose that Assumption 4.2 holds. Then, there exists constants c 1 and c 2 such that

$$\frac { \gamma _ { k } ^ { \top } s _ { k } } { \| s _ { k } \| ^ { 2 } } \geqslant c _ { 1 } \| d _ { k } \| _ { B _ { k } } ^ { 2 },$$

$$\frac { \| \gamma _ { k } \| ^ { 2 } } { \gamma _ { k } ^ { \top } s _ { k } } \leqslant \frac { c _ { 2 } } { \| d _ { k } \| _ { B _ { k } } ^ { 2 } }.$$

Proof. From Wolfe line search, (3.10) and (4.9), we can get

$$f _ { i } ( x _ { k } ) - f _ { i } ( x _ {$$

$$\gamma _ { k } ^ { \top } s _ { k } \geqslant$$

Let c 1 = σ α 1 1 &gt; 0, it follows by Proposition 4.1 that

γ

⊤

k

s

k

⩾

2

∥

s

c

1

∥

d

k

∥

.

k

∥

2

B

k

Next, we consider that inequality (4.15). By the definition of the 4.5,

$$\begin{array} {$$

Since ∇ f i ( x ) is Lipschitz continuous with constant L i , for i ∈ [ m ], we can get

$$\frac { \| y _ { k } \| } { \| s _ { k } \$$

By the Cauchy-Schwarz inequality,

$$\frac { | \langle y _ { k }, s _ { k } \rangle |$$

Therefore,

$$\frac { \| y _ { k } \| } { \| s _ { k } \| } + \max \{ - \eta _ { k }, 0 \} \leqslant \frac { 2 \| y _ { k } \| } { \| s _ { k } \| } \leqslant 2 \max _ { i \in [ m ] } L _ { i }.$$

Assumption 4.1 and Assumption 4.2 (i) imply the existence of f 1 ∈ R such that f i ( x k ) ⩾ f 1 for all k ⩾ 0 and i ∈ [ m ]. It follows by (4.16) and (4.17) that

$$\| \gamma _ { k } \| \leqslant ( 2 \max _ { i \in [ m ] } L _ { i } + \max _ { i \in [ m ] } ( f _ { i } ( x _ { 0 } ) - f _ { 1 } ) ) \| s _ { k } \|,$$

hence,

$$\frac { \| \gamma _ { k } \| ^ { 2 } } { \gamma _ { k } ^ { \top } s _ { k } } \leqslant \frac { c _ { 2 } } { \| d _ { k } \| _ { B _ { k } } ^ { 2 } },$$

$$\text{where } c _ { 2 } = \frac { ( 2 \max _ { i \in [ m ] } L _ { i } + \max _ { i \in [ m ] } ( f _ { i } ( x _ { 0 } ) - f _ { 1 } ) ) ^ { 2 } } { c _ { 1 } }.$$

$$\frac { i \in [ m ] } { c _ { 1 } }.$$

Remark 4.2. β k is the angle between s k and B s k k . By Lemma 4.2 and Proposition 4.2, we can obtain that the angles β k remain away from 0, for an arbitrary fraction p of the iterates.

Lemma 4.3. Let { x k } be the sequence generated by Algorithm 1. Then,

$$\mathcal { D } ( x _ { k }, d _ { k } ) \leqslant - \cos \beta _ { k } \| d _ { k } \| \| d _ { S D } ^ { k } \|.$$

Proof. Since s k = x k +1 -x k = α d k k , then

$$\cos \beta _ { k } = \frac { s _ { k } ^ { \top } B _ { k } s _ { k } } { \| s _ { k } \| \| B _ { k } s _ { k } \| } = \frac { d _ { k } ^ { \top } B _ { k } d _ { k } } { \| d _ { k } \| \| B _ { k } d _ { k } \| }.$$

It follows that

This completes the proof.

$$\mathcal { D } ( x _ { k }, d _ { k } ) & \leqslant - \| d _ { k } \| _ { B _ { k } } ^ { 2 } \\ & = - \cos \beta _ { k } \| d _ { k } \| \| B _ { k } d _ { k } \| \\ & = - \cos \beta _ { k } \| d _ { k } \| \left \| \sum _ { i \in [ m ] } \lambda _ { i } ^ { k } \nabla f _ { i } ( x _ { k } ) \right \| \\ & \leqslant - \cos \beta _ { k } \| d _ { k } \| \| d _ { S D } ^ { k } \|. \\ \dots \sim \epsilon$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-055b28cd-b67d-4464-a474-c466feebfd48.jpg)

The following result establishes that the sequence { x , d k k } satisfies the Zoutendijk condition, which is the basis for establishing the global convergence of Algorithm 1.

Lemma 4.4. Let { x k } be the sequence generated by Algorithm 1. Suppose that Assumption 4.1 and Assumption 4.2 hold. Then

$$\i n c { n } \\ \sum _ { k \geqslant 0 } \frac { \mathcal { D } ( x _ { k }, d _ { k } ) ^ { 2 } } { \| d _ { k } \| ^ { 2 } } < \infty.$$

Proof.

The proof is similar to that in [15, Proposition 3.3], we omit it here.

The following theorem establishes the global convergence of Algorithm 1.

Theorem 4.2. Let { x k } be the sequence generated by Algorithm 1. Suppose that Assumption 4.1 and Assumption 4.2 hold. Then, the following statements hold.

- (i) lim inf k →∞ ∥ d k SD ∥ = 0 ,
- (ii) Every accumulation point of the sequence { x k } is Pareto critical.

Proof. Assume that there exists a constant ε &gt; 0 such that

$$| |.$$

$$\| d _ { S D } ^ { k } \| \geqslant \varepsilon, \quad \forall k \geqslant 0.$$

By the definition of d k and Assumption 4.2 (ii), we get

$$\i \text{ of } d _ { k } \text{ and Assumption } 4. 2 \text{ (ii), we get } \\ - \frac { 1 } { 2 } \| d _ { k } \| ^ { 2 } & = \min _ { d \in \mathbb { R } ^ { n } } \max _ { i \in [ m ] } \nabla f _ { i } ( x _ { k } ) ^ { \top } d + \frac { 1 } { 2 } \| d \| ^ { 2 } _ { B _ { k } } \\ & \leqslant \min _ { d \in \mathbb { R } ^ { n } } \max _ { i \in [ m ] } \nabla f _ { i } ( x _ { k } ) ^ { \top } d + \frac { b } { 2 } \| d \| ^ { 2 } \\ & \leqslant \frac { 1 } { b } \min _ { d \in \mathbb { R } ^ { n } } \max _ { i \in [ m ] } \nabla f _ { i } ( x _ { k } ) ^ { \top } ( b d ) + \frac { 1 } { 2 } \| b d \| ^ { 2 } \\ & = - \frac { 1 } { 2 b } \| d _ { S D } ^ { k } \| ^ { 2 }. \\ d C _ { 2 } = \underline { C } _ { 2 } ^ { b } \text{ if follows by Proposition } 4. 2 \text{ that }$$

Set C 1 = c 1 ε 2 b and C 2 = c 2 b ε 2 , it follows by Proposition 4.2 that

$$\frac { \gamma _ { k } ^ { \top } s _ { k } } { \| s _ { k } \| ^ { 2 } } \geqslant C _ { 1 } > 0,$$

$$\frac { \dots } { \| s _ { k } \| ^ { 2 } } \geq \cup _ { 1 } > 0, \\ \frac { \| \gamma _ { k } \| ^ { 2 } } { \gamma _ { k } ^ { \top } s _ { k } } \leqslant C _ { 2 }. \\ \dots \dots \dots + \dots + \dots \dots.$$

Then, by Lemma 4.3 and Lemma 4.2, there exist a constant δ &gt; 0 and a infinite subsequence of indices K ⊂ N such that

$$- \frac { \mathcal { D } ( x _ { k }, d _ { k } ) } { \| d _ { k } \| } \geqslant \delta \| d _ { S D } ^ { k } \| \geqslant \delta \varepsilon, \quad \forall k \in \mathcal { K },$$

it follows that

$$\sum _ { k \in \mathbb { N } } \frac { \mathcal { D } ( x _ { k }, d _ { k } ) ^ { 2 } } { \| d _ { k } \| ^ { 2 } } \geqslant \sum _ { k \in \mathcal { K } } \frac { \mathcal { D } ( x _ { k }, d _ { k } ) ^ { 2 } } { \| d _ { k } \| ^ { 2 } } \geqslant \sum _ { k \in \mathcal { K } } \delta ^ { 2 } \varepsilon ^ { 2 } = \infty,$$

this contradicts the result of Lemma 4.4. Then,

$$\liminf _ { k \to \infty } \| d _ { S D } ^ { k } \| = 0.$$

Assumption 4.1 implies there exists a infinite subsequence of indices K 1 ⊂ N such that x k K 1 -→ x ∗ and ∥ d k SD ∥ -→ K 1 0 . By the Lemma 3.1, d SD ( x ) is continuous, we can get

$$d _ { S D } ( x ^ { * } ) = 0.$$

Furthermore, by Lemma 4.1, x ∗ is Pareto critical.

We have obtained the weak convergence of the sequence { x k } generated by Algorithm 1. Under strict convexity, the following theorem demonstrates that the sequence { x k } generated by Algorithm 1 strongly converges to the Pareto optimal point x . ∗

Theorem 4.3. Let { x k } be the sequence generated by Algorithm 1. Suppose that f ( x ) is strictly convex, Assumption 4.1 and Assumption 4.2 hold. Then, { x k } converges to a Pareto optimal point x ∗ of (MOP) .

Proof.

$$P r o o f. \text{ The proof is similar to that in [25, Theorem 3.6], we omit it here.}$$

## 4.2 Superlinear convergence analysis

In this section, we consider the local superlinear convergence of Algorithm 1 under Assumption 4.3.

Assumption 4.3. (i) F is twice continuously differentiable.

- (ii) There exists constants U, L &gt; 0 such that UI n ⪯ ∇ 2 f i ( x ) ⪯ LI , n for i ∈ [ m. ]
- (iii) ∇ 2 f i ( x ) is Lipschitz continuous with constant M i , ∈ [ m. ]
- (iv) For any x ∈ R n . there exists constant a and b such that

$$a I _ { n } \preceq B ( x ) \preceq b I _ { n }.$$

The Assumption 4.3 ( ii ) implies that f i ( x ) is strongly convex with modulus U and ∇ 2 f i ( x ) is Lipschitz continuous with constant L , for i ∈ [ m ].

Next, we investigate the R-linear convergence of Algorithm 1.

Lemma 4.5. Suppose that Assumption 4.1 and Assumption 4.3 hold. Let { x k } be the sequence generated by Algorithm 1. Then, for any k ⩾ 0 ,

- (i) ∥ x k -x ∗ ∥ ⩽ 2 U ∥ d k SD ∥ ,
- (ii) ∥ s k ∥ ⩾ 1 -σ 2 L cos β k ∥ d k SD ∥ .

Proof. The two assertions can be obtained by using the same arguments as in the proof of [26, Lemma 4.4], we omit them here.

Proposition 4.3. Let { x k } be the sequence generated by Algorithm 1. Suppose that Assumption 4.1 and Assumption 4.3 hold. Then, there exists a constant r such that

$$\frac { \gamma _ { k } ^ { \top } s _ { k } } { \| s _ { k } \| ^ { 2 } } \geqslant U,$$

$$\frac { \| \gamma _ { k } \| ^ { 2 } } { \gamma _ { k } ^ { \top } s _ { k } } \leqslant \frac { r ^ { 2 } } { U }.$$

Proof. By the definition of γ k , s k and m &gt; k 0, we have

$$\frac { \gamma _ { k } ^ { \top } s _ { k } } { \| s _ { k } \| ^ { 2 } } = \frac { ( y _ { k } + m _ { k } s _ { k } ) ^ { \top } s _ { k } } { \| s _ { k } \| ^ { 2 } } = \frac { y _ { k } ^ { \top } s _ { k } } { \| s _ { k } \| ^ { 2 } } + m _ { k } \geqslant U.$$

Since each objective function is U-strongly convex, for i ∈ [ m ] and any x ∈ R n , f i ( x ) has a lower bound f min i .

It follows from the proof of Proposition 4.2 that

$$\| \gamma _ { k } \| \leqslant r \| s _ { k } \|,$$

where r = 2 max i ∈ [ m ] L i +max( i ∈ [ m ] f i ( x 0 ) -f min i ), the proof is complete.

Theorem 4.4. Suppose that Assumption 4.1 and Assumption 4.3 hold. Let { x k } be the sequence generated by Algorithm 1 and x ∗ is defined as Theorem 4.3. Then, (i) { x k } converges R-linearly to x , ∗

$$\stackrel { \nu } { ( i i ) } \sum _ { k \geqslant 0 } ^ { \nu } \| x _ { k } - x ^ { * } \| < \infty.$$

Proof. The proof of the theorem is similar to the corresponding discussion in [14, Proposition 4.1], so we omit the proof here.

Remark 4.3. Using the convexity of f ( x ), Proposition 4.1 and Assumption 4.3( iv ), we can get

$$0 & \leqslant \sum _ { i \in [ m ] } \lambda _ { i } ^ { k } ( f _ { i } ( x _ { k } ) - f _ { i } ( x _ { k + 1 } ) ) & \leqslant ( \sum _ { i \in [ m ] } \lambda _ { i } ^ { k } \nabla f _ { i } ( x _ { k } ) ) ^ { \top } \left ( x _ { k } - x _ { k + 1 } \right ) \\ & = d _ { k } ^ { \top } B _ { k } s _ { k } \\ & \leqslant \frac { b } { \alpha _ { 1 } } \| s _ { k } \| ^ { 2 } \quad \forall k \geqslant 0, \quad \text{(4.24)} \\, & \quad \, \stackrel { \sim } { \top }$$

where λ k = ( λ , λ k 1 k 2 , · · · , λ k m ) ⊤ is the solution of (3.9). In the convex case, we can observe that η k &gt; 0 and by Theorem 4.3 get

$$\| s _ { k } \| = \| x _ { k + 1 } - x _ { k } \| \rightarrow 0, \, k \rightarrow \infty.$$

Then, i.e.,

$$m _ { k } & = \max \{ - \eta _ { k }, 0 \} + \sum _ { i \in [ m ] } \lambda _ { i } ^ { k } ( f _ { i } ( x _ { k } ) - f _ { i } ( x _ { k + 1 } ) ) \\ & = \sum _ { i \in [ m ] } \lambda _ { i } ^ { k } ( f _ { i } ( x _ { k } ) - f _ { i } ( x _ { k + 1 } ) ) \to 0, \ k \to \infty.$$

in (4.5). Based on the analysis results above, for a sufficiently large k , B k +1 s k ≈ y k . Then (4.8) can be seen as an approximation of the secant equation in [3]. Remarkably, for some convex

$$\sum _ { i \in [ m ] } \lambda _ { i } ^ { k } ( f _ { i } ( x _ { k } ) - f _ { i } ( x _ { k + 1 } ) ) \to 0, \ k \to \infty,$$

functions, both MFQNMO and MQNMO in [3] may exhibit comparable convergence rates. This observation can be supported by the numerical experiments presented in the following section. It can be seen that MFQNMO not only ensure the effectiveness of the algorithm for convex problems, but also guarantee the global convergence for nonconvex MOPs.

Now, we show some useful tools that will be used in superlinear convergence analysis. These tools help to characterize the Dennis-Mor´ condition and demonstrate the convergence rate. e Denote the average matrix ¯ G k as follows

i

$$\bar { G } _ { i } ^ { k } = \int _ { 0 } ^ { 1 } \nabla ^ { 2 } f _ { i } ( x _ { k } + \tau s _ { k } ) d \tau, \ \ i \in [ m ].$$

Then,

$$\bar { G } _ { i } ^ { k } s _ { k } & = \mu _ { i } ^ { k }.$$

Next, we define another two auxiliary functions as follows

$$f _ { \lambda } ^ { k } ( x ) = \sum _ { i \in [ m ] } \lambda _ { i } ^ { k } f _ { i } ( x ),$$

$$G _ { \lambda } ^ { k } = \int _ { 0 } ^ { 1 } \sum _ { i \in [ m ] } \lambda _ { i } ^ { k } \nabla ^ { 2 } f _ { i } ( x _ { k } + \tau s _ { k } ) d \tau.$$

where λ k = ( λ , λ k 1 k 2 , · · · , λ k m ) ⊤ is the solution of (3.9). Let us define

$$\bar { s } _ { k } = \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { \frac { 1 } { 2 } } s _ { k }, \ \bar { y } _ { k } = \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } } y _ { k }, \ \bar { \gamma } _ { k } = \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } } \gamma _ { k },$$

$$\bar { B } _ { k } = \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } } B _ { k } \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } },$$

$$\cos \bar { \beta } _ { k } = \frac { \bar { s } _ { k } ^ { \top } \bar { B } _ { k } \bar { s } _ { k } } { \| \bar { s } _ { k } \| \| \bar { B } _ { k } \bar { s } _ { k } \| }, \quad \bar { q } _ { k } = \frac { \bar { s } _ { k } ^ { \top } \bar { B } _ { k } \bar { s } _ { k } } { \bar { s } _ { k } ^ { \top } \bar { s } _ { k } }, \quad \bar { a } _ { k } = \frac { \bar { \gamma } _ { k } ^ { \top } \bar { s } _ { k } } { \| \bar { s } _ { k } \| ^ { 2 } }, \quad \bar { b } _ { k } = \frac { \| \bar { \gamma } _ { k } \| ^ { 2 } } { \bar { \gamma } _ { k } ^ { \top } \bar { s } _ { k } }.$$

By the BFGS-tpye matrix updated fomula (4.6), it can be noted that (4.26) leads to the formula

$$\bar { B } _ { k + 1 } = \bar { B } _ { k } - \frac { \bar { B } _ { k } \bar { s } _ { k } \bar { s } _ { k } ^ { \top } \bar { B } _ { k } } { \bar { s } _ { k } ^ { \top } \bar { B } _ { k } \bar { s } _ { k } } + \frac { \bar { \gamma } _ { k } \bar { \gamma } _ { k } ^ { \top } } { \bar { \gamma } _ { k } ^ { \top } \bar { s } _ { k } }.$$

Another useful tool used in [4], given that the following function

$$\psi ( B ) = \text{trace} ( B ) - \ln ( \det ( B ) ),$$

$$\psi ( \bar { B } _ { k + 1 } ) = \psi ( \bar { B } _ { k } ) + ( \bar { b } _ { k } - \ln ( \bar { a } _ { k } ) - 1 ) + \ln ( \cos ^ { 2 } \bar { \beta } _ { k } ) + \left [ 1 - \frac { \bar { q } _ { k } } { \cos ^ { 2 } \bar { \beta } _ { k } } + \ln ( \frac { \bar { q } _ { k } } { \cos ^ { 2 } \bar { \beta } _ { k } } ) \right ] \quad ( 4. 2 7 )$$

In the next theorem, we prove that the Dennis-Mor´ condition of superlinear convergence e holds, a similar result is discussed in [25] and [4].

Theorem 4.5. Suppose that Assumption 4.1 and Assumption 4.3 hold. Let { x k } be the sequence generated by Algorithm 1. Then

$$\lim _ { k \to \infty } \frac { \| ( B _ { k } - \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ) s _ { k } \| } { \| s _ { k } \| } = 0.$$

Proof. By (4.25), we can get

$$\bar { G } _ { \lambda } ^ { k } s _ { k } = y _ { k },$$

hence,

$$( \bar { G } _ { \lambda } ^ { k } - \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ) s _ { k } = y _ { k } - \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) s _ { k }, \ \ k \geqslant 0.$$

It follows by Assumption 4.3 (iii) , (4.28) and the definition of ¯ , ¯ y k s k that

$$\dots \cdots \cdots \cdots _ { r } \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon \\ \| \bar { y } _ { k } - \bar { s } _ { k } \| & = \| \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } } ( \bar { G } _ { \lambda } ^ { k } - \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ) \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } } \bar { s } _ { k } \| \| \bar { G } _ { \lambda } ^ { k } - \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) \| \\ & \leqslant \| \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } } \| ^ { 2 } \| \bar { s } _ { k } \| \int _ { 0 } ^ { 1 } \| x _ { k } + \tau s _ { k } - x ^ { * } \| d \tau \\ & \leqslant \bar { c } \| \bar { s } _ { k } \| \varepsilon _ { k }, \\ _ { \, - \,. \, \quad, \, \quad -$$

where ¯ = c M ∥∇ 2 f k λ ( x ∗ ) -1 2 ∥ 2 , ε k = max {∥ x k +1 -x ∗ ∥ ∥ , x k -x ∗ ∥} . Since

$$| \| \bar { y } _ { k } \| - \| \bar { s } _ { k } \| | \leqslant \| \bar { y } _ { k } - \bar { s } _ { k } \|,$$

by (4.29), we get

$$( 1 - \bar { c } \varepsilon _ { k } ) \| \bar { s } _ { k } \| \leqslant \| \bar { y } _ { k } \| \leqslant ( 1 + \bar { c } \varepsilon _ { k } ) \| \bar { s } _ { k } \|,$$

$$( 1 - \bar { c } \varepsilon _ { k } ) ^ { 2 } \| \bar { s } _ { k } \| ^ { 2 } - 2 \bar { y } _ { k } ^ { \top } \bar { s } _ { k } + \| \bar { s } _ { k } \| ^ { 2 } \leqslant \| \bar { y } _ { k } \| ^ { 2 } - 2 \bar { y } _ { k } ^ { \top } \bar { s } _ { k } + \| \bar { s } _ { k } \| ^ { 2 } \leqslant \bar { c } ^ { 2 } \| \bar { s } _ { k } \| ^ { 2 } \varepsilon _ { k } ^ { 2 }.$$

It follows that

$$\frac { \bar { y } _ { k } ^ { \top } \bar { s } _ { k } } { \| \bar { s } _ { k } \| ^ { 2 } } \geqslant 1 - \bar { c } \varepsilon _ { k }.$$

Using the Assumption 4.3 (ii) , and (4.31) we obtain

$$\bar { a } _ { k } = \frac { \bar { \gamma } _ { k } ^ { \top } \bar { s } _ { k } } { \| \bar { s } _ { k } \| ^ { 2 } } = \frac { ( \bar { y } _ { k } + m _ { k } \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { - 1 } \bar { s } _ { k } ) ^ { \top } \bar { s } _ { k } } { \| \bar { s } _ { k } \| ^ { 2 } } & = \frac { \bar { y } _ { k } ^ { \top } \bar { s } _ { k } } { \| \bar { s } _ { k } \| ^ { 2 } } + m _ { k } \frac { \bar { s } _ { k } ^ { \top } \nabla ^ { 2 } f _ { \lambda } ^ { k } ( x ^ { * } ) ^ { - 1 } \bar { s } _ { k } } { \| \bar { s } _ { k } \| ^ { 2 } } \\ & \geqslant 1 - \bar { c } \varepsilon _ { k } + \frac { m _ { k } } { L } \\ & \geqslant 1 - \bar { c } \varepsilon _ { k }.$$

By the L-Lipschitz continuity of f i , we have

$$\frac { \| y _ { k } \| } { \| s _ { k } \| } \leqslant L.$$

From (4.21), (4.33) and the right side of (4.30),

$$& \text{From } ( 4. 2 1 ), ( 4. 3 3 ) \text{ and the right side of } ( 4. 3 0 ), \\ & \bar { b } _ { k } = \frac { \| \bar { \gamma } _ { k } \| ^ { 2 } } { \bar { \gamma } _ { k } \bar { s } _ { k } } \\ & \quad = \frac { \| \bar { y } _ { k } \| ^ { 2 } } { \bar { \gamma } _ { k } \bar { s } _ { k } } + 2 m _ { k } \frac { \bar { y } _ { k } ^ { T } \nabla ^ { 2 } f _ { X } ^ { k } ( x ^ { * } ) ^ { - 1 } \bar { s } _ { k } } { \bar { \gamma } _ { k } ^ { T } \bar { s } _ { k } } + ( m _ { k } ) ^ { 2 } \frac { \bar { s } _ { k } ^ { T } \nabla ^ { 2 } f _ { X } ^ { k } ( x ^ { * } ) ^ { - 2 } \bar { s } _ { k } } { \bar { \gamma } _ { k } ^ { T } \bar { s } _ { k } } \\ & \quad \leq \frac { ( 1 + \bar { c } _ { c } _ { k } ) ^ { 2 } } { 1 - \bar { c } _ { c } _ { k } } + 2 m _ { k } \frac { \| \nabla ^ { 2 } f _ { X } ^ { k } ( x ^ { * } ) ^ { - 1 } \frac { \| \| y _ { k } \| \nabla ^ { 2 } f _ { X } ^ { k } ( x ^ { * } ) ^ { - 1 } \| \| \nabla ^ { 2 } f _ { X } ^ { k } ( x ^ { * } ) ^ { \frac { 1 } { 2 } \| \| \nabla \| _ { s } \| } } { U \| _ { s } k \| ^ { 2 } } + ( m _ { k } ) ^ { 2 } \frac { \| \nabla ^ { 2 } f _ { X } ^ { k } ( x ^ { * } ) ^ { - 2 } \| \| \bar { s } _ { k } \| ^ { 2 } } { U \| _ { s } k \| ^ { 2 } } \\ & \quad \leq 1 + \frac { 3 \bar { c } + \bar { c } ^ { 2 } \varepsilon _ { k } } { 1 - \bar { c } _ { c } _ { k } } \varepsilon _ { k } + \eta _ { 1 } m _ { k } + \eta _ { 2 } m _ { k } ^ { 2 },$$

where η 1 = 2 ∥ ∇ 2 f k λ ( x ∗ ) -1 2 ∥∥∇ 2 f k λ ( x ∗ -) 1 ∥∥∇ 2 f λ ( x ∗ ) 1 2 ∥ L U , η 2 = ∥∇ 2 f k λ ( x ∗ -) 2 ∥∥∇ 2 f k λ ( x ∗ ) 1 2 ∥ 2 U . Using (4.27) and (4.32),

$$\psi ( \bar { B } _ { k + 1 } ) \leqslant \psi ( \bar { B } _ { k } ) & + \left ( \frac { 3 \bar { c } + \bar { c } ^ { 2 } \varepsilon _ { k } } { 1 - \bar { c } \varepsilon _ { k } } \varepsilon _ { k } + \eta _ { 1 } m _ { k } + \eta _ { 2 } m _ { k } ^ { 2 } - \ln ( 1 - \bar { c } \varepsilon _ { k } ) \right ) \\ & + \ln ( \cos ^ { 2 } \bar { \beta } _ { k } ) + \left [ 1 - \frac { \bar { q } _ { k } } { \cos ^ { 2 } \bar { \beta } _ { k } } + \ln \left ( \frac { \bar { q } _ { k } } { \cos ^ { 2 } \bar { \beta } _ { k } } \right ) \right ],$$

By adding the above inequality from 1 to k , then we obtain

$$0 < \psi ( \bar { B } _ { k + 1 } ) \leqslant \psi ( \bar { B } _ { 1 } ) & + \sum _ { j = 1 } ^ { k } \left ( \frac { 3 \bar { c } + \bar { c } ^ { 2 } \varepsilon _ { j } } { 1 - \bar { c } \varepsilon _ { j } } \varepsilon _ { j } + \eta _ { 1 } m _ { j } + \eta _ { 2 } m _ { j } ^ { 2 } - \ln ( 1 - \bar { c } \varepsilon _ { j } ) \right ) \\ & + \left [ 1 - \frac { \bar { q } _ { j } } { \cos ^ { 2 } \bar { \beta } _ { j } } + \ln \left ( \frac { \bar { q } _ { j } } { \cos ^ { 2 } \bar { \beta } _ { j } } \right ) \right ] + \ln \left ( \cos ^ { 2 } \bar { \beta } _ { j } \right ),$$

From Theorem 4.4, there exists k 0 ∈ N such that

$$\bar { c } \varepsilon _ { k } < \frac { 1 } { 2 }, \ \forall k > k _ { 0 },$$

then,

We can get

$$\sum _ { k \geqslant 1 } m _ { k } < \infty,$$

from (4.24) and Theorem 4.4. It follows by (4.36) and (4.37) that

$$\sum _ { j > k _ { 0 } } \left ( \frac { 3 \bar { c } + \bar { c } ^ { 2 } \varepsilon _ { j } } { 1 - \bar { c } \varepsilon _ { j } } \varepsilon _ { j } + \eta _ { 1 } m _ { j } + \eta _ { 2 } m _ { j } ^ { 2 } - \ln ( 1 - \bar { c } \varepsilon _ { j } ) \right ) & < \sum _ { j > k _ { 0 } } 9 \bar { c } \varepsilon _ { j } + \sum _ { j > k _ { 0 } } ( \eta _ { 1 } m _ { j } + \eta _ { 2 } m _ { j } ^ { 2 } ) \\ & < \sum _ { j = 1 } ^ { \infty } 9 \bar { c } \varepsilon _ { j } + \sum _ { j = 1 } ^ { \infty } ( \eta _ { 1 } m _ { j } + \eta _ { 2 } m _ { j } ^ { 2 } ) \\ & < \infty. \\ \sum _ { j = 1 } ^ { \prime \prime } \infty \quad \cdots \quad \cdots.$$

By (4.38) and (4.35), we get

$$\sum _ { j = 1 } ^ { \infty } \left ( - \left [ 1 - \frac { \bar { q } _ { j } } { \cos ^ { 2 } \bar { \beta } _ { j } } + \ln \left ( \frac { \bar { q } _ { j } } { \cos ^ { 2 } \bar { \beta } _ { j } } \right ) \right ] - \ln \left ( \cos ^ { 2 } \bar { \beta } _ { j } \right ) \right ) < \infty.$$

Since -ln(cos 2 ¯ ) β j &gt; 0 ,for all j &gt; k 0 , and h t ( ) = 1 -t + lnt ⩽ 0 when t &gt; 0, hence, we get

$$- \left [ 1 - \frac { \bar { q } _ { j } } { \cos ^ { 2 } \bar { \beta } _ { j } } + \ln \left ( \frac { \bar { q } _ { j } } { \cos ^ { 2 } \bar { \beta } _ { j } } \right ) \right ] \geqslant 0,$$

$$\ln ( 1 - \bar { c } \varepsilon _ { k } ) \geqslant - 2 \bar { c } \varepsilon _ { k }.$$

$$\lim _ { k \to \infty } \ln \left ( \frac { 1 } { \cos ^ { 2 } \bar { \beta } _ { k } } \right ) = 0, \quad \lim _ { k \to \infty } \left [ 1 - \frac { \bar { q } _ { k } } { \cos ^ { 2 } \bar { \beta } _ { k } } + \ln \left ( \frac { \bar { q } _ { k } } { \cos ^ { 2 } \bar { \beta } _ { k } } \right ) \right ] = 0.$$

The above relations imply that

$$\lim _ { k \to \infty } \ln \left ( \cos ^ { 2 } \bar { \beta } _ { j } \right ) = 1, \quad \lim _ { k \to \infty } \bar { q } _ { k } = 1.$$

Now,

$$& \text{Now,} \\ & \quad \text{In} \frac { \| ( B _ { k } - \nabla ^ { 2 } f _ { k } ^ { k } ( x ^ { * } ) ) s _ { k } \| ^ { 2 } } { \| s _ { k } \| ^ { 2 } } & = \lim _ { k \to \infty } \frac { \| \nabla ^ { 2 } f _ { k } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } } ( B _ { k } - \nabla ^ { 2 } f _ { k } ^ { k } ( x ^ { * } ) ) \nabla ^ { 2 } f _ { k } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } } \bar { s } _ { k } \| ^ { 2 } } { \| \nabla ^ { 2 } f _ { k } ^ { k } ( x ^ { * } ) ^ { - \frac { 1 } { 2 } } s _ { k } \| ^ { 2 } } \\ & \quad = \lim _ { k \to \infty } \frac { \| ( \bar { B } _ { k } - I _ { n } ) \bar { s } _ { k } \| ^ { 2 } } { \| \bar { s } _ { k } \| ^ { 2 } } \\ & \quad = \lim _ { k \to \infty } \frac { \| \bar { B } _ { k } \bar { s } _ { k } \| ^ { 2 } - 2 \bar { s } _ { k } ^ { \top } \bar { B } _ { k } \bar { s } _ { k } + \| \bar { s } _ { k } \| ^ { 2 } } { \| \bar { s } _ { k } \| ^ { 2 } } \\ & \quad = \lim _ { k \to \infty } \left [ \frac { \bar { q } _ { k } ^ { 2 } } { \cos ^ { 2 } \bar { \beta } _ { k } } - 2 \bar { q } _ { k } + 1 \right ] \\ & \quad = 0 \\ & \quad \text{This completes the proof.}$$

This completes the proof.

Theorem 4.6. Suppose that Assumption 4.1 and Assumption 4.3 hold. Let { x k } be the sequence generated by Algorithm 1. Then,

- (i) (4.1) and (4.2) hold with α k = 1 for sufficiently large k .
- (ii) { x k } converges to x ∗ at a superlinear rate.

Proof.

The proof is similar to that in [25, Theorem 4.3], we omit it here.

## 5 Numerical experiments

In this section, we present some numerical experiments to verify the efficiency of MFQNMO. We compare MFQNMO with QNMO [24] with Wolfe line search and MQNMO [3] with Wolfe line search. The selected test problems include convex and nonconvex problems. Our evaluation criteria include the method's effectiveness in finding optimal solutions and its ability to approximate real Pareto surfaces accurately. These experiments were implemented using Python 3.9.

For MQNMO and MFQNMO methods, we used the Frank-Wolf gradient method to solve the dual problem (3.9). For QNMO, (3.4) is solved using scipy.optimize, which is an optimization Solver. All test problems and main characteristics are listed in Table 1. Columns 'n' and 'm' indicate the number of variables and objectives, respectively. The 'Convex' column indicates whether the problem is convex or not. 'Y' stands for convex MOP. 'N' stands for nonconvex MOP. Problems are solved under box constraints x L ⩽ x ⩽ x U . The initial points are randomly selected within these box constraints x L ⩽ x ⩽ x U . We used the termination criterion ∣ ∣ θ ( x ) ∣ ∣ ⩽ 10 -8 for all test problems. The maximum number of iterations is set to 500. To obtain the stepsize that satisfies the Wolfe conditions (4.1) and (4.2), we used the Wolf line search algorithm proposed in [14]. We used σ 1 = 10 -4 , σ 2 = 0 1, and set . B 0 = I n . We perform 200 computations using the same initial points for all methods.

Table 1: The description of all problems used in numerical experiments

| problem   |   n |   m | x L                     | x U               | convex   | Reference   |
|-----------|-----|-----|-------------------------|-------------------|----------|-------------|
| SD        |   4 |   2 | (1 , - √ 2 , - √ 2 , 1) | (3 , 3 , 3 , 3)   | Y        | [30]        |
| PNR       |   2 |   2 | ( - 2 , - 2)            | (2 , 2)           | Y        | [22]        |
| JOS1a     |  50 |   2 | ( - 2 , · · · , - 2)    | (2 , · · · , 2)   | Y        | [12]        |
| JOS1b     | 100 |   2 | ( - 2 , · · · , - 2)    | (2 , · · · , 2)   | Y        | [12]        |
| DGO1      |   1 |   2 | - 10                    | 13                | N        | [9]         |
| DGO2      |   1 |   2 | - 9                     | 9                 | Y        | [9]         |
| Lov1      |   2 |   2 | ( - 10 , - 10)          | (10 , 10)         | Y        | [17]        |
| Lov2      |   2 |   2 | ( - 0 . 75 , - 0 . 75)  | (0 . 75 , 0 . 75) | N        | [17]        |
| Lov3      |   2 |   2 | ( - 20 , - 20)          | (20 , 20)         | N        | [17]        |
| Lov4      |   2 |   2 | ( - 20 , - 20)          | (20 , 20)         | N        | [17]        |
| SK1       |   1 |   2 | - 100                   | 100               | N        | [9]         |
| BK1       |   2 |   2 | ( - 5 , - 5)            | (10 , 10)         | Y        | [9]         |
| SLCDT1    |   2 |   2 | ( - 1 . 5 , - 1 . 5)    | (1 . 5 , 1 . 5)   | N        | [29]        |
| MOP1      |   1 |   2 | - 10 5                  | 10 5              | Y        | [9]         |
| MOP2      |   2 |   2 | ( - 4 , - 4)            | (4 , 4)           | N        | [9]         |
| LDTZ      |   3 |   3 | (0 , 0 , 0)             | (1 , 1 , 1)       | N        | [16]        |
| Hil1      |   2 |   2 | (0 , 0)                 | (1 , 1)           | N        | [10]        |
| AP2       |   1 |   2 | - 100                   | 100               | Y        | [3]         |
| AP3       |   2 |   2 | ( - 100 , - 100))       | (100 , 100)       | N        | [3]         |
| FF1       |   2 |   2 | ( - 1 , - 1)            | (1 , 1)           | N        | [9]         |
| KW2       |   2 |   2 | ( - 3 , - 3)            | (3 , 3)           | N        | [11]        |
| MHHM1     |   1 |   3 | 0                       | 1                 | Y        | [9]         |
| MHHM2     |   2 |   3 | (0 , 0)                 | (1 , 1)           | Y        | [9]         |

Table 2: Number of average iterations (iter), average CPU time (time(ms)), number of average function evaluations (feval), and number of failure points (NF) of MFQNMO, QNMO, and MQNMO.

| Problem   | QNMO   | QNMO     | QNMO     | QNMO   | MFQNMO   | MFQNMO   | MFQNMO   | MFQNMO   | MQNMO   | MQNMO   | MQNMO   | MQNMO   |
|-----------|--------|----------|----------|--------|----------|----------|----------|----------|---------|---------|---------|---------|
| Problem   | iter   | time     | feval    | NF     | iter     | time     | feval    | NF       | iter    | time    | feval   | NF      |
| SD        | 13.49  | 486.47   | 15.54    | 1      | 7.35     | 2.55     | 18.51    | 1        | 11.17   | 5.53    | 43.54   | 1       |
| PNR       | 2.56   | 118.66   | 5.48     | 0      | 7.87     | 2.81     | 28.31    | 0        | 8.12    | 3.84    | 54.67   | 0       |
| JOS1a     | 3.01   | 270.14   | 7.39     | 0      | 5.96     | 3.84     | 32.89    | 0        | 3.05    | 5.38    | 24.51   | 0       |
| JOS1b     | 2.89   | 12725.46 | 7.39     | 0      | 7.65     | 19.38    | 10.05    | 0        | 11.23   | 43.337  | 103.01  | 0       |
| DGO1      | 175.51 | 1906.25  | 5058.35  | 70     | 1.39     | 0.64     | 4.14     | 0        | 1.44    | 0.72    | 4.37    | 0       |
| DGO2      | 3.59   | 52.34    | 9.56     | 0      | 5.40     | 1.64     | 17.41    | 0        | 70.35   | 31.75   | 409.96  | 0       |
| Lov1      | 268.91 | 3969.94  | 9213.41  | 107    | 4.27     | 1.45     | 12.72    | 0        | 4.30    | 1.72    | 15.04   | 0       |
| Lov2      | 30.43  | 448.92   | 77.89    | 146    | 12.36    | 8.41     | 31.49    | 124      | 29.70   | 19.44   | 107.27  | 124     |
| Lov3      | 78.77  | 1462.87  | 3400.81  | 30     | 11.86    | 6.27     | 29.42    | 3        | 24.13   | 9.81    | 71.69   | 3       |
| Lov4      | 248.53 | 4516.95  | 12193.75 | 101    | 1.30     | 0.57     | 4.04     | 0        | 1.52    | 0.71    | 4.85    | 0       |
| SK1       | 119.17 | 488.25   | 2747.57  | 69     | 2.36     | 1.65     | 26.54    | 0        | 2.66    | 1.93    | 46      | 0       |
| BK1       | 98.41  | 1085.38  | 100.41   | 39     | 1        | 0.13     | 3        | 0        | 1       | 0.16    | 3       | 0       |
| SLCDT1    | F      | F        | F        | 200    | 3.35     | 0.64     | 8.51     | 0        | 4.18    | 1.38    | 11.28   | 0       |
| MOP1      | 4.66   | 237.83   | 146.13   | 47     | 1.09     | 0.21     | 4.89     | 18       | 1       | 0.15    | 3       | 0       |
| MOP2      | 89.24  | 1489.39  | 34344.11 | 35     | 3.67     | 1.53     | 11.41    | 0        | 7.18    | 2.28    | 121.99  | 1       |
| LDTZ      | F      | F        | F        | 200    | 18.02    | 87.05    | 119.06   | 18       | F       | F       | F       | 200     |
| Hil1      | 258.88 | 7409.41  | 6274.12  | 110    | 11.71    | 6.41     | 62.52    | 16       | 98.39   | 186.89  | 529.98  | 53      |
| AP2       | 1.04   | 35.95    | 3.17     | 1      | 0.985    | 0.41     | 2.95     | 0        | 0.985   | 0.46    | 2.95    | 0       |
| AP3       | 147.91 | 3363.81  | 5353.84  | 115    | 22.36    | 10.41    | 109.45   | 0        | 248.21  | 204.68  | 3858.64 | 77      |
| FF1       | 245.17 | 3394.94  | 10089.34 | 98     | 32.05    | 15.97    | 190.12   | 0        | 16.12   | 8.94    | 81.77   | 0       |
| KW2       | 92.97  | 1697.45  | 298.53   | 103    | 13.81    | 7.98     | 51.01    | 24       | 99.81   | 69.37   | 314.52  | 32      |
| MHHM1     | 0.88   | 26.87    | 2.82     | 0      | 0.88     | 0.40     | 2.65     | 0        | 0.88    | 0.42    | 2.65    | 0       |
| MHHM2     | 245.54 | 4752.67  | 9049.92  | 99     | 6.21     | 43.96    | 206.47   | 31       | 78.34   | 1356.45 | 3508.36 | 31      |

In Table 2, the average values obtained from these runs are recorded to determine the number of iterations (iter), function evaluations (feval), and CPU time (time(ms)). Among the 200 initial points, the number of points that did not converge in 500 iterations is recorded in column 'NF' for MFQNMO, QNMO, and MQNMO. We denote 'Failed' as 'F,' which indicates that the corresponding method fails to solve the test problem under the stopping criterion. For large-scale problems, such as JOS, the QNMO method's CPU time consumption is substantial although it requires the fewest iterations. For all nonconvex test problems in Table 1, the QNMO method exhibited poor performance and the MQNMO method surpasses the MFQNMO method only for the FF1 test problem. For certain problems in Table 1, the experimental results of the MFQNMO method and the MQNMO method are equal approximately, such as BK1, MOP1, MHHM1 and AP2. We mentioned the possible reasons for this phenomenon in Remark 4.3. As can be seen from Table 2, the MFQNMO method is more efficient and competitive.

Figure 1-3 shows the Pareto frontier obtained using the MFQNMO method for some convex problems listed in Table 1. Figure 4-10 shows the Pareto frontier obtained using the MFQNMO method for certain nonconvex problems listed in Table 1. For the chosen test problems, it can be observed that the MFQNMO method was able to satisfactorily estimate the Pareto frontiers when considering evenly distributed starting points.

Figure 1: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the PNR problem.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-aa7b84da-5cbe-486f-9151-8fd7c6b97e96.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-fb6b2155-1a65-4264-b257-66d0579ef7db.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-8d897c54-9476-4ec7-ad69-ba50aaa0340a.jpg)

Figure 2: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the DGO2 problem.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-fdd3b598-8ad2-4342-a784-d650e53275f3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-6ecde381-905a-4460-a617-1082c3789e96.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-dcd6bc46-83a6-48fb-add9-8f4ad106a6c2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-833a2a9d-dca1-47ac-975c-1295d50b4dcb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-3c74914d-381a-442e-86ac-fc71a0dfd7b6.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-8a3267b7-9f2b-4c9b-88ae-d2826f2bd81f.jpg)

Figure 3: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the SD problem.

Figure 4: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the FF1 problem.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-a64267f1-a1f4-4048-b65b-4458e983b7f4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-e76fda43-c266-49da-a38d-ab40facfdeb9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-cb38b582-ceba-4d06-862a-e6fb2fcbd091.jpg)

Figure 5: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the Hil problem.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-e3c3b6a7-8fb5-4661-889c-caa50bc13055.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-dae214a1-c7bf-43b2-9514-a67fbd1e9fba.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-07a41b8f-05a2-40db-972d-c1f9e644dc68.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-913ad10c-9217-4252-8c9d-057862f30744.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-f2906faa-60d1-49b7-b702-2e0eb10c5b24.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-34f3e6b4-aa80-4f7b-8954-5f9e97bdc725.jpg)

Figure 6: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the KW2 problem.

Figure 7: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the SLCDT1 problem.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-420e5ee9-69c3-46ae-8434-e39cb1eccc45.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-0e1eb0f7-8e76-4668-a94e-6ba6c3292bf8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-ede928bf-ae85-41f8-b172-22b698236d26.jpg)

Figure 8: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the MOP2 problem.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-c46f22e5-b512-433e-8e1f-298cc2588a17.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-fb0d3525-0315-4c5d-9316-d55a34c9e415.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-f832572d-34f7-4906-a0d9-7e90221c2f3c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-adf65b33-bd3a-4cea-9836-c9784e2e7d2e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-1de851bb-32c0-4236-93b3-2fade61d4865.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-bbd75b95-d45e-46cd-9bcb-92450a1e3732.jpg)

Figure 9: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the SK1 problem.

Figure 10: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the LDTZ problem.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-68523731-06b7-4d9d-86e6-c55a31851c09.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-abdef0f4-4ed2-4d84-855b-20669c2ff290.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-2575159f-6ae6-40f1-9e72-c9d98609ef0b.jpg)

Figure 11: The approximated nondominated frontiers generated by QNMO, MQNMO, and MFQNMO for the AP3 problem.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-00daff90-b259-49de-b241-f280d80c4e8c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-7ba70177-0790-49de-9767-cecd1cb62f06.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00543v2-figure-58a25cfe-96db-46b9-af08-eb0f7f74e177.jpg)

## 6 Conclusion

We propose a modified BFGS-type method based on function information for nonconvex multiobjective optimization problems (MFQNMO), in which the BFGS-type formula is updated by using

the function information. To obtain the global convergence property, we only need the smoothness of the objective function. Compared with [25], our work has the following advantages.

- (1) Instead of approximating the Hessian matrix of all objective functions, we use a common matrix to capture the curvature information of the objective function, which balances the effectiveness of the algorithm and the computational cost.
- (2) For nonconvex functions, most versions of quasi-Newton methods may not converge, so the MFQNMO method proposed in this paper for nonconvex unconstrained problems is meaningful.
- (3) Compared with [25], the update of the BFGS-type formula does not depend on the selection of parameters, but only needs to use the information of past iterations.

In our future work, we will consider applying the method to nonconvex constrained optimization problems.

## References

- [1] Asha, L.N., Dey, A., Yodo, N., Aragon, L.G.: Optimization approaches for multiple conflicting objectives in sustainable green supply chain management. Sustainability-Base. 14(19), 12790 (2022)
- [2] Assun¸ c˜o, P.B., Ferreira, O.P., Prudente, L.F.: a Conditional gradient method for multiobjective optimization. Comput. Optim. Appl. 78, 741-768 (2021)
- [3] Ansary, M.A., Panda, G.: A modified quasi-Newton method for vector optimization problem. Optimization. 64(11), 2289-2306 (2015)
- [4] Byrd, R.H., Nocedal, J.: A tool for the analysis of quasi-Newton methods with application to unconstrained minimization. SIAM J. Numer. Anal. 26(3), 727-739 (1989)
- [5] Chen, J., Li, G.X., Yang, X.M.: Variable Metric Method for Unconstrained Multiobjective Optimization Problems. J. Oper. Res. Soc. China. 11(3), 409-438 (2023)
- [6] Chen, J., Tang L.p., Yang X.M.: Preconditioned Barzilai-Borwein Methods for Multiobjective Optimization Problems.
- [7] Fliege J., Svaiter B.F.: Steepest descent methods for multicriteria optimization. Math. Meth. Oper. Res. 51, 479-494 (2000)
- [8] Fliege, J., Drummond, L. G., Svaiter, B.F.: Newton's method for multiobjective optimization. SIAM J. Optim. 20, 602-626 (2009)
- [9] Huband, S., Hingston, P., Barone, L., While, L.: A review of multiobjective test problems and a scalable test problem toolkit. IEEE T. Evolut. Comput. 10(5), 477-506 (2006)
- [10] Hillermeier, C.: Generalized homotopy approach to multiobjective optimization. J. Optimiz. Theroy. App. 110(3), 557-583 (2001)
- [11] I. Kim and O. de Weck. Adaptive weighted-sum method for bi-objective optimization: Pareto front generation. Struct. Multidiscip. Optim. 29(2), 149-158 (2005)
- [12] Jin, Y., Olhofer, M., Sendhoff, B.: Dynamic weighted aggregation for evolutionary multiobjective optimization: Why does it work and how. In: L.A. Spector, (ed.) Proceedings of the genetic and evolutionary computation conference, pp. 1042-1049. San Francisco, CA, USA (2001)

| [13]   | Leschine, T.M., Wallenius, H., Verdini, W.A.: Interactive multiobjective analysis and assim- ilative capacity-based ocean disposal decisions. Eur. J. Oper. Res. 56(2), 278-289 (1992)                                                                                                            |
|--------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [14]   | Lapucci, M., Mansueto, P.: A limited memory Quasi-Newton approach for multi-objective optimization. Comput. Optim. Appl. 85(1), 33-73 (2023)                                                                                                                                                      |
| [15]   | Lucambio P´rez, e L.R., Prudente, L.F.: Nonlinear conjugate gradient methods for vector optimization. SIAM J. Optim. 28(3), 2690-2720 (2018)                                                                                                                                                      |
| [16]   | Laumanns, M., Thiele, L., Deb, K., Zitzler, E.: Combining convergence and diversity in evolutionary multiobjective optimization. Evol. Comput. 10(3), 263-282 (2002)                                                                                                                              |
| [17]   | Lovison, A.: Singular continuation: Generating piecewise linear approximations to pareto sets via global analysis. SIAM J. Optim. 21(2), 463-490 (2011)                                                                                                                                           |
| [18]   | Marler, R.T., Arora, J.S.: Survey of multi-objective optimization methods for engineering. Struct. Multidiscip. O. 26, 369-395 (2004)                                                                                                                                                             |
| [19]   | Morovati V., Pourkarimi L., Basirzadeh H.: Barzilai and Borwein's method for multiobjective optimization problems. Numer. Algorithms. 72, 539-604 (2016)                                                                                                                                          |
| [20]   | Mahdavi-Amiri, N., Salehi Sadaghiani, F.: A superlinearly convergent nonmonotone quasi- Newton method for unconstrained multiobjective optimization. Optim. Method. Softw. 35(6), 1223-1247 (2020)                                                                                                |
| [21]   | Morovati, V., Basirzadeh, H., Pourkarimi, L.: Quasi-Newton methods for multiobjective optimization problems. 4or-q. J. Oper. Res. 16, 261-294 (2018)                                                                                                                                              |
| [22]   | Preuss, M., Naujoks, B., Rudolph, G.: Pareto set and EMOA behavior for simple multi- modal multiobjective functions. In: Runarsson, T.P., et al. (ed.) Proceedings of the Ninth International Conference on Parallel Problem Solving from Nature (PPSN IX), pp. 513-522. Springer, Berlin. (2006) |
| [23]   | Pereira, J.L.J., Oliver, G.A., Francisco, M.B., Cunha Jr, S.S., Gomes, G.F.: A review of multi-objective optimization: methods and algorithms in mechanical engineering problems. Arch. Comput. Method E. 29(4), 2285-2308 (2022)                                                                 |
| [24]   | Povalej ˇ Z.: Quasi-Newton's method for multiobjective optimization. J. Comput. Appl. Math. 255, 765-777 (2014)                                                                                                                                                                                   |
| [25]   | Prudente, L.F., Souza, D.R.: Global convergence of a BFGS-type algorithm for nonconvex multiobjective optimization problems. Comput. Optim. Appl. 88(3), 719-757 (2024)                                                                                                                           |
| [26]   | Prudente, L.F., Souza, D.R.: A quasi-Newton method with Wolfe line searches for multiob- jective optimization. J. Optim. Theory Appl. 194(3), 1107-1140 (2022)                                                                                                                                    |
| [27]   | Qu, S., Mark G., Felix T.C.: Quasi-Newton methods for solving multiobjective optimization. Oper. Res. Lett. 39(5), 397-399 (2011)                                                                                                                                                                 |
| [28]   | Sener, O., Koltun, V.: Multi-task learning as multi-objective optimization. Adv. Neural Inf. Process. 31, (2018)                                                                                                                                                                                  |

- [29] Sch¨tze, O., Laumanns, M., Coello Coello, C. A., Dellnitz, M., Talbi, E.G.: a Convergence of stochastic search algorithms to finite size Pareto set approximations. J. Global. Optim. 41(4), 559-577 (2008)
- [30] Stadler W., Dauer J.: Multicriteria optimization in engineering: a tutorial and survey. AIAA. 150, 209-249 (1993)
- [31] Tanabe H., Fukuda E.H., Yamashita N.: Proximal gradient methods for multiobjective optimization and their applications. Comput. Optim. Appl. 72, 339-361 (2019)
- [32] Ye, F., Lin, B., Yue, Z., Guo, P., Xiao, Q., Zhang, Y.: Multi-objective meta learning. Adv. Neural Inf. Process. 34, 21338-21351 (2021)
Yingxue Yang
Global convergence of a modified BFGS-type method based on function information for nonconvex multiobjective optimization problems
In this paper, based on function information, we propose a modified BFGS-type
method for nonconvex multiobjective optimization problems (MFQNMO). In the
multiobjective quasi-Newton method (QNMO), each iteration involves separately
approximating the Hessian matrix for each component objective function, which
results in significant storage and computational burdens. MFQNMO employs a
common BFGS-type matrix to approximate the Hessian matrix of all objective
functions in each iteration. This matrix is updated using function information
from the previous step. This approach strikes a balance between efficacy and
computational cost. We confirm the convergence of the method without relying on
convexity assumptions, under mild conditions, we establish a local superlinear
convergence rate for MFQNMO. Furthermore, we validate its effectiveness through
experiments on both nonconvex and convex test problems.