## ON AN INVERSE TRIDIAGONAL EIGENVALUE PROBLEM AND ITS APPLICATION TO SYNCHRONIZATION OF NETWORK MOTION

LUCA DIECI, CINZIA ELIA, AND ALESSANDRO PUGLIESE

Abstract. In this work, motivated by the study of stability of the synchronous orbit of a network with tridiagonal Laplacian matrix, we first solve an inverse eigenvalue problem which builds a tridiagonal Laplacian matrix with eigenvalues λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N and null-vector e = [ 1 . . . 1 ] . Then, we show how this result can be used to guarantee -if possible- that a synchronous orbit of a connected tridiagonal network associated to the matrix L above is asymptotically stable, in the sense of having an associated negative Master Stability Function (MSF). We further show that there are limitations when we also impose symmetry for L .

Notation . We let e k be the k -th column of the identity matrix, and e be the vector of all 1's. Boldface will indicate vectors, whose number of elements will be clear from the context. As needed, we will use Matlab notation for entries of vectors and matrices, e.g. A i, j ( ) will be the ( i, j )-th entry of a matrix A .

## 1. A linear algebra and a dynamical systems problem

The purpose of this work is two-fold. First, we will solve the following inverse eigenvalue problem: ' Given values λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N , find an unreduced tridiagonal matrix L with these as eigenvalues, with positive diagonal and negative off diagonal entries, and whose entries on each row add up to 0 (that is, L e = 0 '. ) Then, we will use this result to: ' Show how and when we can guarantee that a synchronous orbit of a connected tridiagonal network associated to the matrix L above is asymptotically stable, in the sense that the Master Stability Function (MSF) is negative '.

In the process of solving the above problems, we will also realize that it is generally not possible to also require that L is symmetric. Further, we will show that if the synchronous orbit of the network associated to the standard tridiagonal matrix T in (3) does not lead to a negative value of the MSF, then there cannot be any connected network with tridiagonal symmetric Laplacian matrix giving a negative MSF.

From the theoretical point of view, our results are an important addition to the vast literature on inverse eigenvalue problems for tridiagonal matrices, by adding the nontrivial constraint of having a specific null vector. From the application point of view, our results provide a rigorous and clear indication of how to weigh the arcs of a network with the nearest neighbor topology in order to obtain a stable synchronous orbit, thereby overcoming the limitations imposed by requiring symmetry.

A plan of the paper is as follows. In the remaining part of this section, we review some relevant linear algebra and network (dynamical systems) results, in particular about the concept of MSF. Then, in Section 2, we give our main theoretical results about the inverse eigenvalue problem, by giving and rigorously justifying novel algorithms to obtain a tridiagonal matrix with given spectrum

2020 Mathematics Subject Classification. 65F18, 15A18.

Key words and phrases. Tridiagonal network, inverse eigenvalue problem, synchronization.

The authors gratefully acknowledge the support of the School of Mathematics of Georgia Tech. The work has been partially supported by the GNCS-Indam group and the PRIN2022PNR P2022M7JZW ' SAFER MESH- Sustainable mAnagement oF watEr Resources modEls and numerical MetHods ' research grant.

and null vector. In Section 3, we look at some theoretical results about symmetric connected tridiagonal network structure, and in Section 4 we show performance of our techniques. Finally, conclusions are in Section 5.

## 1.1. Linear algebra background.

/negationslash

Definition 1.1. Given a matrix A ∈ R N × N , we say that A is reducible if there exists a permutation P and integers N 1 = 0 = N 2 , N 1 + N 2 = n , such that P T AP = [ A 11 0 0 A 22 ] where A 11 ∈ R N 1 × N 1 , A 22 ∈ R N 2 × N 2 . If a matrix A is not reducible, then it is called irreducible .

/negationslash

Definition 1.2. A matrix L ∈ R N × N is tridiagonal and unreduced if it has the form

$$( 1 ) & L = \left [ \begin{matrix} a _ { 1 } & b _ { 2 } & 0 & 0 & \cdots & 0 \\ c _ { 2 } & a _ { 2 } & b _ { 3 } & 0 & \cdots & 0 \\ 0 & c _ { 3 } & a _ { 3 } & b _ { 4 } & \ddots & 0 \\ \vdots & \ddots & \ddots & \ddots & \ddots & \vdots \\ 0 & \cdots & 0 & c _ { N - 1 } & a _ { N - 1 } & b _ { N } \\ 0 & \cdots & \cdots & 0 & c _ { N } & a _ { N } \end{matrix} \right ] \,, \\ \intertext { w i t h $ k $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ = m $ }$$

/negationslash

$$\begin{bmatrix} \colon & \ddots & \ddots & \ddots & \colon \\ 0 & \cdots & 0 & c _ { N - 1 } & a _ { N - 1 } & b _ { N } \\ 0 & \cdots & \cdots & 0 & c _ { N } & a _ { N } \end{bmatrix} \\ 2, \dots, N. & M o r e o v e r, \ w e \ c a l l \ L \text{ network-tridiag}$$

with b j = 0 , c j = 0 , j = 2 , . . . , N . Moreover, we call L network-tridiagonal if b j &lt; 0 , c j &lt; 0 , j = 2 , . . . , N , and a j = -b j +1 -c j , j = 1 , . . . , N ( c 1 = 0 = b N +1 ), in which case L is singular with null vector e .

/negationslash

Remark 1.3. We note that a network tridiagonal matrix L is effectively the Laplacian matrix of a given connected network. Moreover, we also remark that an unreduced tridiagonal matrix L is also automatically irreducible.

The following well known result (e.g., see [10, 7]) will be used below.

Lemma 1.4. Given a general tridiagonal, unreduced, matrix L as in (1) , then its eigenvalues are real and distinct, hence L is diagonalizable by a real matrix of eigenvectors V : V -1 LV = diag( λ , i i = 1 , . . . , N ) . Moreover, the eigenvalues do not change as long as the products b k c k , k = 2 , . . . , N , do not change either.

Proof. The proof follows from the fact that the characteristic polynomial of (1) can be recursively defined as follows (Sturm sequence):

$$p _ { 0 } ( \lambda ) & = 1 \, \quad p _ { 1 } ( \lambda ) = \lambda - a _ { 1 } \, \\ p _ { j } ( \lambda ) & = ( \lambda - a _ { j } ) p _ { j - 1 } ( \lambda ) - ( b _ { j } c _ { j } ) p _ { j - 2 } ( \lambda ) \, \quad j = 2, \dots, N \.$$

Note that we must have b c j j = 0, for all j = 2 , . . . , N , since L is unreduced.

$$& p _ { 0 } ( \lambda ) = 1 \, \quad p _ { 1 } ( \lambda ) = \lambda - a _ { 1 } \, \\ & p _ { j } ( \lambda ) = ( \lambda - a _ { j } ) p _ { j - 1 } ( \lambda ) - ( b _ { j } c _ { j } ) p _ { j - 2 } ( \lambda ) \, \quad j = 2, \dots, N \. \\ & \text{Note that we must have } b _ { j } c _ { j } \neq 0, \text{ for all } j = 2, \dots, N, \, \text{since } L \text{ is unreduced.}$$

/negationslash

Example 1.5. The most commonly studied instance of network tridiagonal matrix L of interest to us is that associated to symmetric nearest-neighbor coupling, called diffusive coupling in [6] . That is, one has

$$L & = \sigma T, \ \ T = \begin{bmatrix} 1 & - 1 & 0 & 0 & \cdots & 0 \\ - 1 & 2 & - 1 & 0 & \cdots & 0 \\ \vdots & \ddots & \ddots & \ddots & \vdots \\ 0 & \cdots & 0 & - 1 & 2 & - 1 \\ 0 & \cdots & \cdots & 0 & - 1 & 1 \end{bmatrix} \,,$$

where σ &gt; 0 is called the coupling strength. The eigenvalues µ j 's and (orthonormal) eigenvectors v j 's of T are well known:

$$& F o r \ \ j = 1, 2, \dots, N \colon \quad \mu _ { j } = 4 \sin ^ { 2 } \left ( \frac { ( j - 1 ) \pi } { 2 N } \right ) \, \quad a n d \\ & v _ { i 1 } = \sqrt { \frac { 1 } { N } } \, \ v _ { i j } = \sqrt { \frac { 2 } { N } } \cos \left ( \frac { ( i - 0. 5 ) ( j - 1 ) \pi } { N } \right ) \, \ \ i = 1, \dots, N \. \\ 1. 2. \text{ Networks and the MSF. Consider the following system of counted identical different} \quad \text{different} \quad \text{.}$$

1.2. Networks and the MSF. Consider the following system of coupled identical differential equations

$$\dot { x } _ { i } & = f ( x _ { i } ) + \sum _ { j = 1 } ^ { N } a _ { i j } E ( x _ { j } - x _ { i } ), \ i = 1, \dots, N, \\. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.$$

/negationslash where A = ( a ij ) i,j =1: N , is the matrix describing the interaction of N different agents x i ∈ R n , a ii = 0 , i = 1 , . . . , N , a ij ≥ 0 for j = , and i E ∈ R n × n is the matrix describing which components of each agent interact with one another. In general, A represents the structure of a directed graph with weighted edges, and we will henceforth assume that the graph is connected, hence one can get to any node starting from any other node, moving along edges (in particular, no row of A can be 0). We recall that a graph being connected is equivalent to the matrix A being irreducible.

In (5), separately each agent satisfies the same differential equation

$$\dot { x } = f ( x )$$

and by the structure of the system obviously taking N copies of the same solution of (6), x 1 = x 2 = · · · = x N = x , one obtains a solution of (5). This is called synchronous solution and we denote it as x s . Further, with

$$\Xi = \{ x \in \mathbb { R } ^ { n N } \, \colon \, x _ { 1 } = x _ { 2 } = \dots = x _ { N } \, \, x _ { j } \in \mathbb { R } ^ { n }, \, j = 1, \dots, N \}$$

we indicate the synchronous subspace .

/negationslash

It is convenient to rewrite (5) by defining the new matrix L : L ij = -a ij , for i = j and L ii = ∑ j a ij , i, j = 1 , . . . , N . Obviously, 0 is an eigenvalue of L , and since the graph is connected, 0 is a simple eigenvalue of L . In general, L is not symmetric, and it corresponds to what is known as out-degree Laplacian (see [14]). In this work, we will restrict to L satisfying the structure of network tridiagonal matrix given in Definition 1.2.

$$\text{network tridiagonal matrix given in Definition 1.2.} \\ \text{Let } x = \begin{bmatrix} x _ { 1 } \\ \vdots \\ x _ { N } \end{bmatrix} \in \mathbb { R } ^ { n N }, \, F ( x ) = \begin{bmatrix} f ( x _ { 1 } ) \\ \vdots \\ f ( x _ { N } ) \end{bmatrix}, \, \text{and rewrite } ( 5 ) \text{ as} \\ \dot { x } = F ( x ) + M x, \quad \text{where} \quad M = - L \otimes E \.$$

$$\dot { x } = F ( x ) + M x, \quad \text{where} \ \ M = - L \otimes E \.$$

In order to ascertain the stability of a synchronous orbit x s one needs to study the behavior of solutions of (8) transversal to Ξ and the Master Stability Function (MSF) does precisely that. In fact, the MSF tool (originally devised in [11]) is a widely adopted indicator of linearized stability of the synchronous orbit for the system (8). The power of the technique consists in the replacement of the large nN -dimensional linear system arising from linearizing (8), with a single n -dimensional parametrized linear system. Indeed, linearization of (8) about the synchronous solution x S gives

the linear system

$$\dot { y } & = \begin{bmatrix} D f ( x ) & & \\ & D f ( x ) & \\ & & \ddots & \\ & & D f ( x ) \end{bmatrix} y - ( L \otimes E ) \dot { y }, \\ \int _ { t } \, \text{Next, let } V \, \text{ be the matrix of eigenvectors of } L, \, \text{and perform the cha}$$

$$\dot { y } _ { i } = D f ( x ) y _ { i } - \lambda _ { i } E y _ { i } \, \\ \vdots \theta \, \left. \infty \, \left. \infty \, \right. \infty \, \left. \infty \, \left. \, \left. \infty \, \right.$$

where y = [ y 1 . . . y N ] . Next, let V be the matrix of eigenvectors of L , and perform the change of variable ( V -1 ⊗ I n ) y → y , to obtain the N linear systems of dimension n where λ 1 = 0 &lt; λ 2 &lt; . . . λ N are the eigenvalues of L . As a consequence, one considers the single parametrized linear system

$$& ( 9 ) & & \dot { z } = ( A ( t ) - \eta E ) z \, \quad \text{where} \quad \eta \geq 0 \quad \text{and} \quad A ( t ) = D f ( x ( t ) ) \. \\ & \varpi \ e. \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 \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 \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \\ & \varpi \ e. \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 \cdots \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colon \cdots \colots \colon \cdots \colots \colon \cdots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots \colots$$

Definition 1.6. The MSF is defined as the largest Lyapunov exponent (Floquet exponent, in case the synchronous solution is a periodic orbit) of (9) as η ranges over the eigenvalues of L . A negative value of the MSF implies stability of the synchronous orbit.

In this work, we will study directly the parametrized linear system (9) and a-priori decide what range of values of η (if any) will give a negative MSF, then ask whether or not it is possible to find a network tridiagonal structure (that is, a Laplacian matrix of the form in (1)) whose eigenvalues fit the stability region inferred by the MSF. More in details, if we find that there is a negative value of the MSF for η ∈ [ η , η 1 2 ], then the synchronous solution is asymptotically stable for (8) if

$$\frac { \lambda _ { 2 } } { \lambda _ { N } } > \frac { \eta _ { 1 } } { \eta _ { 2 } }.$$

Remark 1.7. Obviously, the MSF depends on E and it is easy to give examples where the MSF is negative for some E , but positive for some other coupling matrices E ; for example, see [8, Figure 1] . Indeed, there is a very vast literature on network synchronization, the MSF, and the interplay between network topology and negative MSF. This is fairly evident already at a graph theoretic level: e.g., by choosing the network structure of a complete graph, and the matrix E to be the identity, surely increases synchronizability. But this is not a desirable way to proceed, since the type of connections between agents is not just a mathematical artifact. For this reason, we will assume that the nearest neighbor interaction (that is, the network tridiagonal form of L ), and the matrix E , are specified by the underlying application and that there is no freedom in changing them, in particular the structure of L . Our goal will be to choose, if possible, the entries of L so that the MSF is negative.

## 2. Our inverse network tridiagonal eigenvalue problem

We will solve our inverse eigenvalue problem in two steps. First, in Section 2.1, we give an algorithm that builds a symmetric, unreduced, tridiagonal matrix with a given set of distinct eigenvalues. Then, in Section 2.2, we propose an algorithm that takes a given symmetric, singular, unreduced, tridiagonal matrix, and it produces an unreduced, generally non-symmetric, tridiagonal matrix with a specified null vector.

There exists an extensive literature in the general area of inverse eigenvalue problems, (e.g., see [2]). In fact, our first algorithm (to build a symmetric tridiagonal unreduced matrix with a preassigned spectrum) is effectively a known result. However, to the best of our knowledge, our result on specifying spectrum and null vector appears to be new, and it is what we need for obtaining a network leading to a negative MSF.

- 2.1. Symmetric unreduced tridiagonal with given spectrum. We are interested in solving the following inverse problem.

Given N real values λ 1 &lt; λ 2 &lt; · · · &lt; λ N , find an unreduced, symmetric, tridiagonal matrix S , with negative off diagonal, having these λ i 's as eigenvalues . To clarify, we are seeking S of the form

$$S = | \begin{array} {$$

with all b i &lt; 0, i = 2 , . . . , N , and eigenvalues λ 1 &lt; λ 2 &lt; · · · &lt; λ N . To solve this problem we adopted Algorithm diag2trid below.

$$S = \begin{bmatrix} a _ { i } & b _ { 2 } & 0 & 0 &$$

## Algorithm 1 Algorithm diag2trid

- 1: Set D = diag( λ , . . . , λ 1 N ).
- 3: Let Q be a Householder reflection such that Q e 1 = q .
- 2: Let q ∈ R N be the vector q = e / √ N .
- 4: Set A = Q DQ T .
- 5: Perform the Householder tridiagonalization algorithm on A , so that

$$H ^ { T } A H = S$$

where S is tridiagonal and H is orthogonal.

- 6: By possibly changing the signs of the b i 's, we further enforce that b i &lt; 0, i = 2 , . . . , N , to obtain the sought form of S .

## Remarks 2.1.

- (i) With the trivial exception of the last step, Algorithm diag2trid is effectively the same as the one we described, and justified, in [3] , where it is also shown to be equivalent to a technique of Schmeisser, see [12] , whose interest was in building an unreduced symmetric tridiagonal matrix whose characteristic polynomial is given. A related construction is also summarized in [1, Theorem 4.7] . We note that performing a possible change of sign of the b i 's is legitimate in light of Lemma 1.4.
- (ii) Although the choice q = e / √ N is our default choice, and the one we adopted in the experiments in Section 4, there is freedom in choosing the unit vector q in Algorithm diag2trid . This is natural, since the transformation H is such that H , (: 1) = H (1 :) , T = e 1 , which means that S is diagonalized by ( QH ) T , whose first row is given by q . Appealing to the Implicit Q Theorem, see [5, Theorem 8.3.2] , we know that a real symmetric tridiagonal matrix is completely characterized by its (real) eigenvalues and by the first row of the (orthogonal) matrix of its eigenvectors, in the sense that any two real symmetric, tridiagonal, and unreduced matrices, S and T , that are diagonalized by two orthogonal matrices having the same first row, must be equal up to the sign of their off diagonal entries.

The next Lemma will be needed below and it clarifies the structure of S produced by Algorithm 1 when the eigenvalues are λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N .

Lemma 2.2. Let S as in (11) be produced by Algorithm 1 relative to eigenvalues λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N . Then, the values a i in S are all positive: a i &gt; 0 , i = 1 , . . . , N .

Proof. By construction, we have S = H AH T = ( QH ) T D QH ( ). Let U = QH , so that S = U T DU . Therefore, for any j = 1 , . . . , N :

$$& \cdot \\ & a _ { j } = e _ { j } ^ { T } S e _ {$$

and so -since λ k &gt; 0 for k = 2 , . . . , N - we have a j &gt; 0 unless u kj = 0 for k = 2 , . . . , N . By contradiction, suppose that u kj = 0 for k = 2 , . . . , N , and some j = 1 , . . . , N . But then U e j = ± e 1 , and since U is orthogonal this means that e T 1 U = ± e T j . But, see Remark 4-(ii), S is diagonalized into D by an orthogonal matrix whose first row is q T , and so this would contradict the choice of q in Step 2 of Algorithm 1, and the thesis follows. /square

- 2.2. Singular tridiagonal with specific null-vector. Next, we consider the following modification of the previous problem.

Given N real values λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N , find an unreduced, network tridiagonal matrix L , whose spectrum is given by these values, and such that the eigenvector associated to the 0 -eigenvalue is aligned with [ 1 1 . . . 1 ] .

Our construction will produce a generally nonsymmetric tridiagonal matrix L as in (1) and we will see that -in general- one cannot require that there is a symmetric tridiagonal matrix satisfying our requests. The technique we used to resolve this problem is encoded in the following Algorithm, which will be justified in Theorem 2.3 below.

## Algorithm 2 TridZeroRowSum

- 1: Using Algorithm 1, diag2trid , generate an unreduced symmetric tridiagonal matrix S as in (11) with eigenvalues λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N .
- (a) Let α 2 such that a 1 + α b 2 2 = 0.
- 2: Modify the first N -1 rows of S as follows.
- (b) For k = 3 , . . . , N , let α k : α b k k + a k -1 + b k -1 α k -1 = 0.
- 3: The desired L is then:

$$& 3 \colon The desired L is then \colon \\ & L = \begin{bmatrix} a$$

Theorem 2.3. The above Algorithm 2, TridZeroRowSum , is well defined and terminates with an unreduced tridiagonal matrix with spectrum given by { 0 , λ 2 , · · · , λ n } , and eigenvector associated to

1

the

0

-eigenvalue given (up to normalization) by

e

=

 

.

.

.

.

1

 

/negationslash

Proof. The algorithm is well defined as long as α i = 0 for all i = 2 , . . . , N . Next, using Lemma 2.5 below, we show that the algorithm must complete the first N -1 steps, that is that, with b 1 = 0,

 

for k = 2 , . . . , N -1 we have:

/negationslash

$$a _ { k - 1 } + \frac { b _ { k - 1 } } { \alpha _ { k - 1 } } \neq 0. \\ \intertext { a, \, } \neq 0 \ B u t { \, } \text{thist fail}.$$

$$\text{there exists, unique, $\alpha_{k}\neq 0$ such that $\alpha_{k}b_{k}+a_{k-1}+\frac{b_{k-1}}{\alpha_{k-1}}=0$,} \\........................................................................................................................................................................................................$$

or, equivalently,

/negationslash

/negationslash

Now, for k = 2, we need to show that a 1 = 0. But this follows from Lemma 2.2. For k &gt; 2, we need to show that

/negationslash

By contradiction, suppose that a k -1 + b k -1 α k -1 = 0 and consider the principal minor of order k -1 of the matrix B k -1 obtained during the reduction process, that is

$$a _ { k - 1 } + \frac { b _ { k - 1 } } { \alpha _ { k - 1 } } \neq 0. \\. \ b _ { \nu - 1 } \cdots \in \cdots$$

$$B _ { k - 1 } = \begin{bmatrix} a _ { 1 } & \alpha _ { 2 } b _ { 2 } & 0 & 0 & \cdots & 0 \\ b _ { 2 } / \alpha _ { 2 } & a _ { 2 } & \alpha _ { 3 } b _ { 3 } & 0 & \cdots & 0 \\ 0 & b _ { 3 } / \alpha _ { 3 } & a _ { 3 } & \alpha _ { 4 } b _ { 4 } & \ddots & 0 \\ \vdots & & \ddots & \ddots & \vdots \\ 0 & \cdots & 0 & b _ { k - 1 } / \alpha _ { k - 1 } & a _ { k - 1 } \end{bmatrix} \,.$$

.

But then this matrix is singular, since we would have that

$$B _ { k - 1 } \begin{bmatrix} 1 \\ 1 \\ \vdots \\ 1 \end{bmatrix} = 0 \\ \intertext { \text{line of $B_{k-1}$ $ But, $1\sim\Gamma$ $C$} } \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \alpha \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Delta \sim \Lambda$$

and in particular 0 would be an eigenvalue of B k -1 . But, the eigenvalues of B k -1 are the same as the eigenvalues of the principal minor of order k -1 of S , that is of

$$S _ { k - 1 } = \begin{bmatrix} a _ { 1 } & b _ { 2 } & 0 & 0 & \cdots & 0 \\ b _ { 2 } & a _ { 2 } & b _ { 3 } & 0 & \cdots & 0 \\ 0 & b _ { 3 } & a _ { 3 } & b _ { 4 } & \ddots & 0 \\ \vdots & & \ddots & \ddots & \ddots & \vdots \\ 0 & \cdots & 0 & b _ { k - 1 } & a _ { k - 1 } \end{bmatrix} \,, \\ \text{e found that a principal minor of $S$ is singular, contradict} \, \text{and} \, \text{and} \, \text{$\frac{1}{2}$ was $0$} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and} \, \text{and}$$

and so we would have found that a principal minor of S is singular, contradicting Lemma 2.5.

Finally, since by construction the first N -1 rows of L have 0-sum and are linearly independent, then the last row must also have 0-sum (since it is linearly dependent on the previous N -1, given that 0 is an eigenvalue).

Therefore, the algorithm produces the matrix L in (12), which clearly has the desired spectrum and rows that sum to zero. This concludes the proof. /square

Lemma 2.4. Let S ∈ R N × N be a symmetric matrix with a 0 eigenvalue. Assume that 0 is a simple eigenvalue and let v be an associated eigenvector of length 1 . If no component of v is 0 , then S is irreducible. Conversely, if S is unreduced and tridiagonal, hence irreducible, with eigenvalues λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N , then v i = 0 , for all i = 1 , . . . , N .

/negationslash

Proof. We show that, if no component of v is 0, then S is irreducible, by showing that if S is reducible, then there exists some i : v = 0. Indeed, if S is reducible, then for some permutation n 1 + n 2 = N . So, we have PSP T P v = 0, or [ S 1 0 0 S 2 ][ w 1 w 2 ] = 0 with [ w 1 w 2 ] = P v . Then, since the

i P we have PSP T = [ S 1 0 0 S 2 ] , with S 1 = S T 1 ∈ R n 1 × n 1 and S 2 = S T 2 ∈ R n 2 × n 2 , n , n 1 2 ≥ 1 and kernel of S is 1-dimensional, we must have either w 1 = 0 or w 2 = 0, giving the claim.

/negationslash

To show the converse statement, for S unreduced and tridiagonal, suppose that v k = 0 for some 1 ≤ k ≤ N . First, note that if k = 1, then - writing S = [ a 1 b T b B ] - since S [ 0 v 2 ] = 0, v 2 ∈ R N -1 = 0, then B v 2 = 0 and this means that B has a 0 eigenvalue, which contradicts Lemma 2.5 below. Next, suppose v k = 0 and k &gt; 1. Partition S as follows:

$$- \circ \max \nu \succ \tilde { \Lambda } \, \cos \max \cos \max \cos. \\ S = \left [ \begin{matrix} T _ { 1 } & \left [ b _ { k } e _ { k - 1 } & 0 \right ] \\ 0 \end{matrix} \right ] & T _ { 2 } \right ] \\ \text{duced, of size } \tilde { \Gamma } _ { \L \L } \, \sin \tilde { \Gamma } _ { \L \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Gamma } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tilde { \Psi } _ { \L } \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \sin \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau } \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \sin \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau \, \cos \tau$$

where T 1 is tridiagonal, unreduced, of size ( k -1 , k -1), and T 2 is tridiagonal, unreduced, of size ( N -k +1 , N -k +1). Writing v =   v 1 0 v 2   , v 1 ∈ R k -1 , v 2 ∈ R N -k , then we must have T 1 v 1 = 0,

/negationslash contradicting Lemma 2.5.

and -because of Lemma 2.5T 1 is invertbile and this implies that v 1 = 0. Thus, we also have T 2 [ 0 v 2 ] = 0 and thus, either v 2 = 0 if T 2 is invertible, contradicting that v = 0, or T 2 is singular, /square

Lemma 2.5. Given a symmetric, unreduced, tridiagonal matrix S as in (11) , with N &gt; 1 , and with eigenvalues λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N . Then, any leading (respectively, trailing) principal submatrix of S of size ( p, p ) , 1 ≤ p ≤ N -1 , is positive definite. 1

Proof. The proof follows from a refinement of classic results on interlacing of eigenvalues for symmetric tridiagonal matrices. In particular, the following result holds (see [7, Problem 4.3.P17]):

M = C c c T d , symmetric, tridiagonal and unreduced, with C R N -1 ,N -1 . Let µ 1 &lt; µ 2 &lt;

· · · &lt; µ N be the eigenvalues of M , and let ν 1 &lt; ν 2 &lt; · · · &lt; ν N -1 be the eigenvalues of C . Then, the ν 's interlace properly the µ 's. That is, we have

'Given [ ] ∈ j i

$$\mu _ { 1 } < \nu _ { 1 } < \mu _ { 2 } < \nu _ { 2 } < \dots \mu _ { N - 1 } < \nu _ { N - 1 } < \mu _ { N }.$$

We show the result for the leading principal submatrices. Let S k , k = 1 2 , . . . , N -1, be the principal submatrices of order k of S , and let λ ( k ) 1 &lt; · · · &lt; λ ( k ) k be their eigenvalues. Using the proper interlacing result quoted above, in particular we must have:

$$\begin{array} { c } \text{The result also holds if we partition $M = \begin{bmatrix} a & b ^ { T } \\ b & C \end{bmatrix}. \end{array} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{                                                                                                                                                                                                        
\end{array}$$

$$\lambda _ { 1 } ^ { ( N ) } = \lambda _ { 1 } = 0 < \lambda _ { 1 } ^ { ( N - 1 ) } < \lambda _ { 1 } ^ { ( N - 2 ) } < \dots < \lambda _ { 1 } ^ { ( 1 ) } \\ \intertext { w o n d o f t u l i n e m i n o 1 c h m o t u l i o o i o d i o i o d i o 1 }$$

and the result follows. The case of trailing principal submatrices is identical. /square

We will also need the following result that refines Lemma 2.4 relative to the eigenvector of S associated to the 0 eigenvalue.

1 A leading (respectively, trailing) principal submatrix of size ( N -p, N -p ), p = 1 , . . . , N -1, is the matrix obtained by deleting the bottom (respectively, top) p rows and columns of S .

Lemma 2.6. Let S be a symmetric, unreduced, tridiagonal matrix as in (11) , with N &gt; 1 , and with eigenvalues λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N . Let v be a unit eigenvector of S associated to the 0 -eigenvalue. Then, the entries of v all have the same sign.

Proof. We are going to use a beautiful relation between the entries of the eigenvector and the Sturm sequence (2). Using [9, Formula (15)], it holds that

$$v _ { k } = c \, \frac { p _ { k - 1 } ( 0 ) } { b _ { 2 } \dots b _ { k } }, \ k = 2, \dots, N,$$

where c is a nonzero constant fixing v 1 . Because of Lemma 2.4, we know that all entries of v are not 0 and thus we can write

$$\frac { v _ { k + 1 } } { v _ { k } } = \frac { 1 } { b _ { k + 1 } } \frac { p _ { k } ( 0 ) } { p _ { k - 1 } ( 0 ) } \. \\ \text{ns in the right-hand-side are} \ n$$

In this last expression, both fractions in the right-hand-side are negative values. In fact, for the first fraction this is obvious, since b k +1 &lt; 0; the second fraction is the ratio between the characteristic polynomials of the leading principal minors of order k and k -1, evaluated at 0. But, because of the proper interlacing result on the eigenvalues of the principal minors (see the proof of Lemma 2.5), the polynomials p k and p k -1 assume opposite values at the origin, and so the ratio p k (0) /p k -1 (0) &lt; 0 and the result follows. /square

Finally, we conclude this section with the following result that summarizes the fact that with our construction we obtain a network tridiagonal matrix satisfying the structural form of Definition 1.2, which is what we wanted to achieve.

Theorem 2.7. The matrix L in (12) satisfies the structural assumptions of Definition 1.2; in particular, all values a i , i = 1 , . . . , N , in (11) are strictly positive and the off diagonal entries α b i i (and of course b /α i i ) , i = 2 , . . . , N , are strictly negative.

Proof. By looking at L in (12), we observe that the a i 's are the a i 's of S produced by Algorithm 1, hence they are strictly positive because of Lemma 2.2. Also, the values of b i , i = 2 , . . . , N , are negative because they come from S . So, we now show that the α i 's are positive and the result will follow, since (by construction) the sum of the entries in each row is 0.

Let v be the eigenvector of S associated to the 0-eigenvalue. As in the proof of Theorem 2.3, and because the α j 's are not zero in light of Lemma 2.4, we have the following relation between the entries of v and the α j 's, for some nonzero value of c :

$$c = v _ { 1 }, \ c \alpha _ { 2 } = v _ { 2 }, \ c \alpha _ { 2 } \alpha _ { 3 } = v _ { 3 }, \dots, \ c \alpha _ { 2 } \cdots \alpha _ { N } = v _ { N } \. \\ \cdot \quad. \.$$

Therefore, we have that and using Lemma 2.6 the result follows.

$$\alpha _ { k } = \frac { v _ { k } } { v _ { k - 1 } } \, \ k = 2, \dots, N, \\ \text{t follows}.$$

/square

2.3. Impact of symmetry. In Example 2.8, we point out that, in general, one cannot also require that the sought tridiagonal matrix be symmetric.

Example 2.8. Suppose that, given λ 1 = 0 &lt; λ 2 &lt; λ 3 , there exists a 3 × 3 real symmetric unreduced tridiagonal matrix T such that (i) T has eigenvalues { 0 , λ 2 , λ 3 } , and (ii) the kernel of T is spanned

$$\ t r d a g o n a l \ m a t r x { \, I } \ s u c h \ t h a t \left ( \imath \, I \ h a s \ e g e n v a l e u s \left \{ 0, \lambda _ { 2 }, \lambda _ { 3 } \right \}, \, a n d \left ( \imath \right ) \ t h e \ k r e n e l \ o f \, I \ s \ s p a n n e d \\ \ b y \begin{bmatrix} 1 \\ 1 \end{bmatrix}. \, T o \ s a t i s f y \, \condition { ( i ) }, T \ m u s t \ h a v e \ t h e \, \text{form} \, T = \begin{bmatrix} x & - x & 0 \\ - x & x + y & - y \\ 0 & - y & y \end{bmatrix}, \, a n d \, \text{to\, satisfy\, also}$$

condition (i), x and y must satisfy { 2( x + ) = y λ 2 + λ 3 3 xy = λ λ 2 3 . (uniquely) two pairs of solutions:

Solving with respect to x and y yields

$$\text{$v$ panes $oy$ sovousons} \colon \\ \begin{cases} x = \frac { 1 } { 1 2 } \left ( \pm \sqrt { 3 } \sqrt { 3 \lambda _ { 2 } ^ { 2 } - 1 0 \lambda _ { 2 } \lambda _ { 3 } + 3 \lambda _ { 3 } ^ { 2 } } + 3 \lambda _ { 2 } + 3 \lambda _ { 3 } \right ) \\ y = \frac { 1 } { 1 2 } \left ( \mp \sqrt { 3 } \sqrt { 3 \lambda _ { 2 } ^ { 2 } - 1 0 \lambda _ { 2 } \lambda _ { 3 } + 3 \lambda _ { 3 } ^ { 2 } } + 3 \lambda _ { 2 } + 3 \lambda _ { 3 } \right ). \\ \i d y \text{ to be real valued, we must have} \end{cases}$$

 But, for x and y to be real valued, we must have

$$1 0 \lambda _ { 2 } \lambda _ { 3 } \leq 3 \lambda _ { 2 } ^ { 2 } + 3 \lambda _ { 3 } ^ { 2 } \\ \rho \quad$$

In conclusion, there exists a real matrix T =   x -x 0 -x x + y -y 0 -y y   having eigenvalues 0 &lt; λ 2 &lt; λ 3 if and only if 10 λ λ 2 3 ≤ 3 λ 2 2 +3 λ 2 3 , which is not necessarily satisfied.

## 3. Symmetric Tridiagonal Networks and the MSF

In this section we give some new results on the ratio between the first positive and the largest eigenvalue of symmetric positive semi-definite network tridiagonal matrices and further highlight the impact that these results have on synchronizability of symmetric tridiagonal networks.

For every N ∈ N , with N ≥ 2, we let S denote the set of real N × N matrices as in (11) that further are network tridiagonal matrices (that is, for which the entries on each row add up to zero) and that are positive semi-definite. Further, we indicate with T the N × N matrix in (3) (obviously, an element of S ).

As usual, we let the eigenvalues be ordered in increasing fashion. Namely, for given S ∈ S , we denote its eigenvalues by

$$\lambda _ { 1 } = 0 < \lambda _ { 2 } < \dots < \lambda _ { N } \text{ \ and for } \ T \text{ \ we call them } \mu _ { 1 } = 0 < \mu _ { 2 } < \dots < \mu _ { N } \text{.}$$

Remark 3.1. We observe that any G ∈ S is a linear combination of elementary matrices E , . . . , E 1 N -1 :

$$& ( 1 3 ) & & G = \sum _ { k = 1 } ^ { N - 1 } x _ { k } E _ { k }, \, \text{for some } x _ { 1 }, \dots, x _ { N - 1 } \in \mathbb { R }, \\ & u h e r e \, \text{for amu } k = 1 \quad & N = 1 \quad E _ { i } \ \text{is the } N \ \times \ N \ \text{matrix} \, \text{different defined} \, \text{hu}$$

where for any k = 1 , . . . , N -1 , E k is the N × N matrix defined by

$$\cdots _ { \sigma } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda } \cdots _ { \lambda \lambda } \cdots _ { \lambda \lambda } \cdots _ { \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \xi } } \\ \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e,w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext { c e, w h e n } \intertext }$$

 For instance, when N = 4 we have

$$E _ { 1 } = \begin{bmatrix} 1 & - 1 & 0 & 0 \\ - 1 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{bmatrix}, \ E _ { 2 } = \begin{bmatrix} 0 & 0 & 0 \\ 0 & 1 & - 1 & 0 \\ 0 & - 1 & 1 & 0 \\ 0 & 0 & 0 \end{bmatrix}, \ E _ { 3 } = \begin{bmatrix} 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & - 1 \\ 0 & 0 & - 1 & 1 \end{bmatrix}.$$

The following result for parameter dependent symmetric matrices with distinct eigenvalues is well known (e.g., see [13])), and it will be used in the proof of Theorem 3.3 below.

Lemma 3.2 (First and second Hadamard variation formulae) . Let A t ( ) be a real symmetric matrix valued function that is twice continuously differentiable in t , and such that A t ( ) has simple eigenvalues for all t . Then, also its eigenvalues λ k ( ) t = λ k ( A t ( )) and corresponding orthonormal eigenvectors u k ( ) = t u k ( A t ( )) depend smoothly on t and we have:

$$\dot { \lambda } _ { k } ( t ) = u _ { k } ^ { T } ( t ) \dot { A } ( t ) u _ { k } ( t ),$$

/negationslash

$$\ddot { \lambda } _ { k } ( t ) & = u _ { k } ^ { T } ( t ) \ddot { A } ( t ) u _ { k } ( t ) + 2 \sum _ { j \neq k } \frac { | u _ { j } ^ { T } ( t ) \dot { A } ( t ) u _ { k } ( t ) | ^ { 2 } } { \lambda _ { k } ( t ) - \lambda _ { j } ( t ) }. \\ \intertext { O o r m o r o n } \, \frac { 1 } { 2 } \, \text{$\text{or all in one}$} \ N > 2 \text{ and with $the$ naturalian above $one$} \ n o n e.$$

Theorem 3.3. For all integers N ≥ 3 , and with the notation above, we have:

and the maximum is attained only at scalar multiples of T .

$$\max _ { S \in \mathcal { S } } \frac { \lambda _ { 2 } ( S ) } { \lambda _ { N } ( S ) } = \frac { \mu _ { 2 } } { \mu _ { N } } = \frac { \sin \left ( \frac { \pi } { 2 N } \right ) ^ { 2 } } { \sin \left ( \frac { ( N - 1 ) \pi } { 2 N } \right ) ^ { 2 } }, \\ \text{and the maximum is attained only at scalar multiplies of $T$}.$$

Proof. With notation from (3), we have

$$T v _ { k } = \mu _ { k } v _ { k }, \text{ for all } k = 1, \dots, N \.$$

Given any S ∈ S that is not a scalar multiple of T , define the following matrix valued function and the following real valued function, which is well defined and strictly positive for all t ∈ [0 , 1]:

$$A ( t ) = \underset { \begin{array} { c } \dots \\ \dots \end{array} } { \dots } \underset { \begin{array} { c } \dots \\ \dots \end{array} } { \dots } \dots$$

$$r ( t ) & = \frac { \lambda _ { 2 } ( A ( t ) ) } { \lambda _ { N } ( A ( t ) ) }.$$

We are going to show that:

- (i) ˙(1) = 0; r
- (ii) r is non-decreasing for all t ∈ [0 , 1];

Then the conclusion of the theorem will follow from (i)-(iii). In particular, (iii) shows that the maximum in (16) is attained only at scalar multiples of T .

- (iii) r is strictly increasing for all t ∈ [1 -δ, 1], for some 0 &lt; δ ≤ 1.

Proof of point (i). Consider the map

↦

$$& \quad \ \ \ \ \ \ \ f \colon x = \left [ \begin{smallmatrix} \begin{smallmatrix} x _ { 1 } \\ x _ { 2 } \\ \vdots \\ x _ { N - 1 } \end{smallmatrix} \right ] \in \mathbb { R } ^ { N - 1 } \mapsto \frac { \lambda _ { 2 } ( T + G ( x ) ) } { \lambda _ { N } ( T + G ( x ) ) }, \\ \text{where} \, \mathcal { C } - \mathcal { C } ( \mathcal { C } ) \text{ to defined as in } T \times \text{ } W a \text{ want to show that } \nabla f - \mathcal { D } \text{ at } \mathcal { C } - \mathcal { D } \end{smallmatrix}$$

where G = G ( x ) is defined as in (13). We want to show that ∇ f = 0 at x = 0. Using (14), we obtain

(19) ∂f ∂x j ∣ ∣ ∣ x =0 = µ N v T 2 E j v 2 -µ 2 v T N E j v N µ 2 N . Hence, we have

(20) ∂f ∂x j ∣ ∣ ∣ x =0 = 0 if and only if v T 2 E j v 2 v T N E j v N = µ 2 µ N . Direct computation yields

$$\frac { \mu _ { 2 } ^ { T } E _ { j } v _ { 2 } } { v _ { N } ^ { T } E _ { j } v _ { N } } = \frac { \left ( v _ { 2 } ( j ) - v _ { 2 } ( j + 1 ) \right ) ^ { 2 } } { \left ( v _ { N } ( j ) - v _ { N } ( j + 1 ) \right ) ^ { 2 } }.$$

Substituting (4) into the right-hand side of (21) and applying trigonometric identities, we obtain

$$2$$

This proves the statement on the right of (20), and hence concludes the proof of point (i).

$$\frac { \left ( v _ { 2 } ( j ) - v _ { 2 } ( j + 1 ) \right ) ^ { 2 } } { \left ( v _ { N } ( j ) - v _ { N } ( j + 1 ) \right ) ^ { 2 } } = \frac { \left ( \cos \frac { ( 2 j - 1 ) \pi } { 2 N } - \cos \frac { ( 2 j + 1 ) \pi } { 2 N } \right ) ^ { 2 } } { \left ( \cos \frac { ( N - 1 ) ( 2 j - 1 ) \pi } { 2 N } - \cos \frac { ( N - 1 ) ( 2 j + 1 ) \pi } { 2 N } \right ) ^ { 2 } } = \\ = \frac { \left ( \sin \frac { j \pi } { N } \sin \frac { \pi } { 2 N } \right ) ^ { 2 } } { \left ( \sin \frac { ( N - 1 ) j \pi } { N } \sin \frac { ( N - 1 ) \pi } { 2 N } \right ) ^ { 2 } } = \frac { \left ( \sin \frac { \pi } { 2 N } \right ) ^ { 2 } } { \left ( \sin \frac { ( N - 1 ) \pi } { 2 N } \right ) ^ { 2 } } = \frac { \mu _ { 2 } } { \mu _ { N } }. \\ \text{This proves the statement on the right of } ( 2 0 ), and hence concludes the proof of point \left ( i i. } \text{Recall that } A ( t ) \text{ is linear in } t. \text{ and that } 0 = \lambda _ { 1 } ( t ) < \lambda _ { 2 } ( t ) < \dots < \lambda _ { N } ( t ) \text{ for}$$

Proof of point (ii). Recall that A t ( ) is linear in t , and that 0 = λ 1 ( ) t &lt; λ 2 ( ) t &lt; . . . &lt; λ N ( ) for t all t ∈ (0 1]. , Then, it follows from (15) that λ N ( ) and t λ 2 ( ) are, respectively, convex and concave t functions of t in [0 , 1]. Differentiating r , we obtain

$$\dot { r } = \frac { \dot { \lambda } _ { 2 } \lambda _ { N } - \dot { \lambda } _ { N } \lambda _ { 2 } } { \lambda _ { N } ^ { 2 } }.$$

We have

$$\div \ddot { \lambda } _ { 2 } \lambda _ { N } - \ddot { \lambda } _ { N } \lambda _ { 2 } \leq 0. \\ \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \div \colon$$

$$& \text{Define } \varphi ( t ) \coloneqq \dot { \lambda } _ { 2 } ( t ) \lambda _ { N } ( t ) - \dot { \lambda } _ { N } ( t ) \lambda _ { 2 } ( t ) \text{ for all } t \in [ 0, 1 ]. \text{ We have} \\ & \dot { \varphi } = \ddot { \lambda } _ { 2 } \lambda _ { N } - \ddot { \lambda } _ { N } \lambda _ { 2 } \leq 0.$$

Since ˙(1) = 0, we have r ϕ (1) = 0, and therefore ϕ t ( ) ≥ 0 for all t ∈ [0 , 1]. This proves ˙( r t ) ≥ 0 for all t ∈ [0 , 1], i.e. r is non-decreasing on [0 , 1].

Proof of point (iii). We are going to show that ¨ λ 2 ( ) t &lt; 0 in (1 -δ, 1) for some δ &gt; 0. This implies that r is strictly decreasing in (1 -δ, 1), and hence r (0) &gt; r (1).

/negationslash

$$s _ { 1 }, \dots, s _ { N - 1 }, \alpha & \text{ as follows:} \\ ( 2 4 ) & \begin{bmatrix} w _ { 1 } - w _ { 2 } & & & - w _ { 1 } \\ - w _ { 1 } + w _ { 2 } & w _ { 2 } - w _ { 3 } & & & - w _ { 2 } \\ - w _ { 2 } + w _ { 3 } & \ddots & & \vdots \\ & & \ddots & w _ { N - 1 } - w _ { N } & - w _ { N - 1 } \\ & & - w _ { N - 1 } + w _ { N } & - w _ { N } \end{bmatrix} \begin{bmatrix} s _ { 1 } \\ s _ { 2 } \\ \vdots \\ s _ { N - 1 } \\ \alpha \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \\ \vdots \\ 0 \\ 0 \end{bmatrix}. \\ \text{Performing row-reduction on (24), we obtain the equivalent system:}$$

We know that ¨ λ 2 (1) ≤ 0. Suppose ¨ λ 2 (1) = 0. Because of (15), we have that v T j ( T -S ) v 2 = 0 for all j = 2 (recall that v T 1 ( S -T ) is the zero vector). Then, we have that ( T -S ) v 2 must be a scalar multiple of v 2 , possibly zero. In any case, this means that v 2 is an eigenvector of S . A non-zero vector w = [ w . 1 . . w N ] ∈ R N is an eigenvector of S if and only if there exists α ∈ R such that ( S -αI ) w = 0. Setting s j = -S j,j +1 for all j = 1 , . . . , N -1, we recast this as a linear system in s , . . . , s 1 N -1 , α as follows:

Performing row-reduction on (24), we obtain the equivalent system:

$$( 2 5 ) & \left [ \begin{matrix} w _ { 1 } - w _ { 2 } & & - w _ { 1 } \\ & w _ { 2 } - w _ { 3 } & & - w _ { 1 } - w _ { 2 } \\ & \ddots & \vdots \\ & & - w _ { 1 } - w _ { 2 } - \dots - w _ { N } \end{matrix} \right ] \left [ \begin{matrix} s _ { 1 } \\ s _ { 2 } \\ \vdots \\ s _ { N - 1 } \end{matrix} \right ] = \left [ \begin{matrix} 0 \\ 0 \\ 0 \\ 0 \end{matrix} \right ]. \\ \text{Let us denote the matrix in (25) by W and set $w=v_{2}$. Then, since $w_{j} \neq w_{j+1}$ for all $j=1,\dots,N-1$} \\ ( \text{see (4)}). and $v_{1}^{T}v_{2}=0$. it follows that rank(W) = N-1$. Recall that a non-trivial solution to}$$

/negationslash

Let us denote the matrix in (25) by W and set w = v 2 . Then, since w j = w j +1 for all j = 1 , . . . , N -1 (see (4)), and v T 1 v 2 = 0, it follows that rank( W ) = N -1. Recall that a non-trivial solution to

(25) is known and is given by s 1 = . . . = s N -1 = 1 and α = µ 2 . It follows that S must be a scalar multiple of T , which is something we have excluded. Therefore, we must have ¨ λ 2 (1) &lt; 0. Being ¨ λ 2 continuous, then it must remain strictly negative in a (left) neighborhood of 1. This concludes the proof. /square

Corollary 3.4. Suppose we are given 0 = ν 1 &lt; ν 2 &lt; . . . &lt; ν N , with

$$& \supset \dots \sum _ { 1 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N }, \, \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N }, \, \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N }, \, \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu ^ { 2 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu ^ { 2 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu ^ { 2 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu ^ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nmid \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu { \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 n } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu \dots \nu _ { 1 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 N } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 4 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \narrow \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 1 N } \nu _ { 2 N } \nu _ { 3 } \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ { 1 } \nu _ { 2 } \nu _ { 3 } \dots \nu _ { 5 } \nu _ { 6 } \nu _ { 7 } \nu _ { 8 } \nu _ { 9 } \nu _ {$$

Then, there is no matrix S ∈ S such that λ j ( S ) = ν j , j = 1 , . . . , N .

Remark 3.5. For our purposes, the key consequence of Corollary 3.4 is the following. Suppose that there is a finite interval of values η for which the MSF is negative. Then, because of (10) , if

$$\ e v e \ e v e \ e _ { J } \ e v e \ e v e \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ { J } \ e _ } \\ \ e t o \ e b t a i n \ a \ n e g a t i v e \ v a l u e \ e f o t \ e \ M S F \ w i t h \ a n y \ s y m m e t r i c \ t r i.$$

it is not possible to obtain a negative value of the MSF with any symmetric tridiagonal network Laplacian.

## 4. Numerical Results

Here we show how our technique works in practice. We give two examples, one is a network of van der Pol oscillators with periodic synchronous orbit, the other is a network of R¨ssler oscillators with o chaotic synchronous orbit. Among our goals is to show that, in general, the symmetric tridiagonal structure (3) may fail to give stability of the synchronous orbit (that is, it won't give a negative value of the MSF, no matter how large is σ in (3)), but our technique is always able to give a negative value of the MSF, if there is an η -interval where the MSF is negative.

For both our examples, we proceed as follows.

- 1. We consider the parametrized linear system (9) and compute the MSF in function of η .
- 2. If there exist any, we select an interval [ η , η 1 2 ] so that for η ∈ [ η , η 1 2 ], the MSF is negative.
- 4.1. Van der Pol. We consider a network of N identical Van der Pol oscillators. The single agent satisfies the following equation
- 3. We use Algorithm TridZeroRowSum to build the Laplacian L as in (12) with nonzero eigenvalues in [ η , η 1 2 ], and null vector 1 √ N [ 1 . . . 1 ] .

$$\lambda \cdots \hat { \Phi } _ { 1 } \Phi _ { 2 } \Phi _ { 3 } \Phi _ { 4 } \Phi _ { 5 } \Phi _ { 6 } \Phi _ { 7 } \Phi _ { 8 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phicon _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9} \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi_ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 } \Phi _ { 9 }$$

$$\begin{matrix} \end{matrix}$$

and we choose the following coupling matrix E = ( 0 0 1 0 ) . In Figure 1 on the left we plot the MSF

for η ∈ [0 , 0 5]. . The MSF is negative for η &gt; 0 39 and remains negative. . at least for η ≤ 50. If we couple the agents via symmetric diffusive coupling as in Example 1.5, then in order to synchronize the network we must impose σλ 2 &gt; 0 39 with . σ constant coupling strength. We then need large values of the coupling strength, namely σ &gt; 40 5 for . N = 32. We instead employ our technique and consider an asymmetric tridiagonal coupling.

We have some freedom on what values we select for the N -1 eigenvalues of L . We experimented with many different choices for these values and below we report on two different experiments: i) select the eigenvalues linearly spaced in [1 , 10], and ii) take the eigenvalues to be Chebyshev points

Figure 1. Van der Pol, N=32. Left: MSF. Right: 2-norm of the difference between agents.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01066v2-figure-8d8162ac-2ad5-4e9c-8feb-7bcc035ecd3d.jpg)

of the first kind (rescaled to [1 , 10] or to larger intervals, as needed). For N = 32, all the elements of the coupling matrix are less than 6 2. . We take an initial condition on the attractor and perturb it with a normally distributed perturbation vector. Then we integrate the network with a 4-th order Runge Kutta method and fixed stepsize h . For these methods, theoretical results guarantee that for h sufficiently small the numerical method has a closed invariant curve and that the distance of this curve from the synchronous periodic orbit is O h ( 4 ) (see [4]).

We synchronize 32 agents with both choices of eigenvalues. To witness, in Figure 1 on the right we plot max i =1 ,..., 31 ‖ x i ( ) t -x i +1 ( ) t ‖ 2 at the grid points. The greatest value of the MSF is obtained for η = 1, and it is /similarequal -0 13. . We plot e -0 13 . t as well in order to appreciate the convergence speed to the synchronous solution. The plots are obtained for one initial condition, but the behavior we observe is consistent for every normally distributed random perturbation we considered. The transient in this case is relatively short.

4.2. R¨ ossler. Next, we consider a network of N identical R¨ssler oscillators. o Each agent satisfies the system

$$( 2 8 ) & & \begin{cases} \dot { y } _ { 1 } = - y _ { 2 } - y _ { 3 } \\ \dot { y } _ { 2 } = y _ { 1 } + 0. 2 y _ { 2 } \\ \dot { y } _ { 3 } = 0. 2 + ( y _ { 1 } - 9 ) y _ { 3 } \end{cases} \\ \end{cases} \\ \text{and we choose the running matrix } E \equiv \mathbb { R } ^ { 3 \times 3 } \text{ in } ( 8 \rangle \text{ given } h v \text{ } E = \begin{cases} 1 \\ 0 \end{cases}$$

left, we plot the MSF in function of η . The MSF is negative in the interval η ∈ [0 19 . 4 61]. .

and we choose the coupling matrix E ∈ R 3 × 3 in (8) given by E =   1 0 0 0 0 0 0 0 0   . In Figure 2, on the

If we use diffusive coupling with a constant coupling strength σ , like in Example 1.5, then (see the explicit values of the eigenvalues given in Example 1.5), in order to synchronize N agents, we need σλ 2 &gt; 0 19 and . σλ N &lt; 4 61. . A necessary condition for N agents to synchronize is then, see (10), λ N λ 2 &lt; 4 61 . 0 19 . and as soon as N &gt; 7 this condition is not satisfied and hence we cannot hope to synchronize more than 7 agents with diffusive coupling and constant coupling strength. In [11] the

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01066v2-figure-86585ff4-170b-4849-9850-0e527c337032.jpg)

Figure 2. Left: MSF for R¨ ossler. Right: 2-Norm of the difference between agents.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01066v2-figure-12f37b29-48d8-466f-adb3-e31a06ec92a2.jpg)

authors use diffusive coupling in a circular array and can synchronize up to 10 agents, but already with 16 agents the necessary conditions for synchronization are not met anymore. 2

In what follows we couple and synchronize a network of 64 R¨ssler agents with our technique, o about the chaotic orbit of a single R¨ssler oscillator. o Of course, as we observed in our numerical experiments, the lack of symmetry causes a large transient and the basin of attraction of the synchronous periodic orbit is in general affected by it.

We report on experiments with two different sets of eigenvalues: linearly spaced in [0 5 . , 3], and Chebyshev points of the first kind (rescaled to [0 5 . , 3]). We consider initial conditions obtained by adding a normally distributed perturbation with variance 1 of a synchronous initial condition and integrate (8) to verify that indeed we obtain synchronization. In Figure 2 we plot max i =1 ,..., 63 ‖ x i ( ) t -x i +1 ( ) t ‖ 2 for the two choices of eigenvalues. For η = 0 5, the corresponding . value of the MSF is /similarequal -0 15 and the dashed line in the right plot is the graph of . e -0 15 . t . It is clear that, after a transient, the convergence speed of the two perturbed solutions is also e -0 15 . t .

## 5. Conclusions

In this work, we gave two contributions. First, we built a tridiagonal matrix L with eigenvalues λ 1 = 0 &lt; λ 2 &lt; · · · &lt; λ N and null-vector e = [ 1 . . . 1 ] . Then, we used this result to achieve -if possiblethat a synchronous orbit of a tridiagonal network associated to the matrix L above is asymptotically stable, in the sense of having an associated negative Master Stability Function (MSF). Our matrix L , in general, is not symmetric, and we further gave some results highlighting the limitations occurring if we also require symmetry. From the practical point of view, the main consequence of our results is that, if at all possible to have a negative MSF with tridiagonal network structure for the Laplacian, then we can have a stable synchronous motion, and we gave numerical examples to show this. At the same time, we also observed that, when we increase the number of agents, it is not always possible to achieve synchronization with arbitrary precision. We leave investigation of this fact to a future numerical study.

## References

- [1] M. Chu and G.H. Golub . Structured inverse eigenvalue problems. Acta Numerica , pp. 1-71, 2002.

2 The authors of [11] point our that, with symmetric tridiagonal Laplacians there will always be an upper limit in the size of a the network in order to obtain a stable synchronous chaotic orbit. With our technique, in principle there is not such limitation.

- [2] M. Chu and G.H. Golub . Inverse Eigenvalue Problems. Oxford Science Publications. Oxford University Press, New York, 2013.
- [3] L. Dieci and A. Pugliese . Forming a symmetric, unreduced, tridiagonal matrix with a given spectrum. arXiv preprint 2311.02677, 2023.
- [4] T. Eirola . Invariant curves of one-step methods. BIT , 28-1, pp. 113-122 (1988).
- [5] L.H. Golub and C.F. Van Loan . Matrix Computations. Johns Hopkins Studies in the Mathematical Sciences. Johns Hopkins University Press, 2013.
- [6] J. Hale. Diffusive Coupling, Dissipation, and Synchronization. Journal of Dynamics and Differential Equations , 9-1, pp. 1-52 (1997).
- [7] R. Horn and C. Johnson . Matrix Analysis, 2nd Edition. Cambridge Univ Press 2013.
- [8] L. Huang, Q. Chen, Y-C Lai and L.M. Pecora . Generic behavior of master-stability functions in coupled nonlinear dynamical systems. Physical Review E , 80, 036204-1-11 (2009).
- [9] L.G. Molinari . Lesson 6: Tridiagonal Matrices. https://api.semanticscholar.org/CorpusID:232239257 , 2019.
- [10] T. Muir . Treatise on the Theory of Determinants. Revised and Enlarged by William H. Metzler. New York: Dover Publications 1960.
- [11] L. Pecora and T. Carroll. Master Stability Functions for Synchronized Coupled Systems. Physical Review Letters , vol. 80, n. 10, (1997), pp. 2109-1112.
- [12] G. Schmeisser . A real symmetric tridiagonal matrix with a given characteristic polynomial. Linear Algebra Appl.s , 193, pp. 11-18 (1993).
- [13] T. Tao . Topics on Random Matrix Theory. Graduate Studies in Mathematics, Vol. 132, American Mathematical Society, Providence, RI, 2012.
- [14] J.J.P. Veerman and R. Lyons . A primer on Laplacian dynamics in directed graphs, Nonlinear Phenomena in Complex Systems , 2020, url=https://api.semanticscholar.org/CorpusID:211066395.

School of Mathematics, Georgia Tech, Atlanta, GA 30332 U.S.A. Email address : dieci@math.gatech.edu

Dipartimento di Matematica, Univ. of Bari, I-70100, Bari, Italy Email address : cinzia.elia@uniba.it

Dipartimento di Matematica, Univ. of Bari, I-70100, Bari, Italy Email address : alessandro.pugliese@uniba.it
Luca Dieci Cinzia Elia Alessandro Pugliese
On an inverse tridiagonal eigenvalue problem and its application to synchronization of network motion
In this work, motivated by the study of stability of the synchronous orbit of
a network with tridiagonal Laplacian matrix, we first solve an inverse
eigenvalue problem which builds a tridiagonal Laplacian matrix with eigenvalues
$\lambda_1=0<\lambda_2<\cdots <\lambda_N$ and null-vector $\boldsymbol{e} =
\begin{bmatrix} 1 \\ \vdots \\ 1 \end{bmatrix}$. Then, we show how this result
can be used to guarantee -- if possible -- that a synchronous orbit of a
connected tridiagonal network associated to the matrix $L$ above is
asymptotically stable, in the sense of having an associated negative Master
Stability Function (MSF). We further show that there are limitations when we
also impose symmetry for $L$.