On Row Sequences of Hermite–Padé Approximation and Its Generalizations

Hermite–Padé approximation has been a mainstay of approximation theory since the concept was introduced by Charles Hermite in his proof of the transcendence of e in 1873. This subject occupies a large place in the literature and it has applications in different subjects. Most of the studies of Hermite–Padé approximation have mainly concentrated on diagonal sequences. Recently, there were some significant contributions in the direction of row sequences of Type II Hermite–Padé approximation. Moreover, various generalizations of Type II Hermite–Padé approximation were introduced and studied on row sequences. The purpose of this paper is to reflect the current state of the study of Type II Hermite–Padé approximation and its generalizations on row sequences. In particular, we focus on the relationship between the convergence of zeros of the common denominators of such approximants and singularities of the vector of approximated functions. Some conjectures concerning these studies are posed.


Introduction
Charles Hermite [1] was the first who introduced the idea of Hermite-Padé approximation. Particularly, he used it for systems of exponential functions to prove that e is transcendental. A formal study of Hermite-Padé approximation for general systems of functions was initiated by Mahler [2] (see also the papers by his students, Coates and Jager [3,4], for important results in this regard). There are basically two different types of Hermite-Padé approximation, namely Types I and II. However, the combination of both (called mixed type) is possible. In this paper, we primarily consider Type II.
An optimal choice of the coefficients of the denominators and numerators of Hermite-Padé approximants makes it an important tool to study analytic continuation of functions and localization of their singularities. We can make this point more clearly by considering classical Padé approximation (scalar Hermite-Padé approximation) stated as follows.
Given a formal Taylor series at the origin for any integers n, m ≥ 0, we can find polynomials P ∈ P n and polynomials Q ∈ P m , Q ≡ 0, such that (QF − P)(z) = O(z n+m+1 ), as z → 0 (P n is the set of all polynomials of degree at most n). The rational function R n,m := P Q = P n,m Q n,m is uniquely defined and is called the (n, m) classical Padé approximant of F. The polynomials P n,m and Q n,m in Equation (2) are selected so that Q n,m is monic and gcd(P n,m , Q n,m ) = 1. For a function F as in Equation (1), we denote by R 0 (F) the radius of the largest open disk at the origin to which F can be extended analytically and by R m (F) the radius of the largest open disk at the origin to which F can be extended so that F has at most m poles counting multiplicities. Set B R := {z ∈ C : |z| < R}.
By Q m (F), we denote the monic polynomial whose zeros are the poles of F in B R m (F) counting multiplicities. The set of all distinct zeros of Q m (F) is denoted by P m (F). Analytic continuation and locations of poles including their multiplicities of F(z) = ∑ ∞ k=0 f k z k in B R m (F) when F has exactly m poles in B R m (F) can be completely described by the convergences of {R n,m } n≥0 and {Q n,m } n≥0 in the following theorem (see [29,30]). Theorem 1 (Montessus de Ballore-Gonchar's Theorem). Let F be defined as in Equation (1) and fix m ∈ N. Then, the following statements are equivalent: (a) R 0 (F) > 0 and F has exactly m poles in B R m (F) counting multiplicities.
There is a polynomial Q m of degree m, Q m (0) = 0, such that the sequence of {Q n,m } n≥0 satisfies lim sup n→∞ Q m − Q n,m 1/n = θ < 1, where · denotes the coefficient norm in the space of all polynomials.
Moreover, if either (a) or (b) holds, then Q m = Q m (F), and lim sup n→∞ F − R n,m where K is any compact subset of B R m (F) \ P m (F) and · K denotes the sup-norm on K.
In [29], Montessus de Ballore proved that (a) implies (b) with Q m = Q m (F) and the inequalities with the sign "≤" instead of "=" in Equations (3) and (4). This part of the above theorem is commonly known as Montessus de Ballore's theorem. The implication (b)⇒(a) with the inequalities with the sign "≥" instead of "=" in Equations (3) and (4) was proved by Gonchar in [30]. In this current paper, we refer to this part as the inverse statement of Montessus de Ballore-Gonchar's Theorem. We also note that the sequence {R n,m (z)} n≥0 diverges at every point z, |z| > R m (F) (see Section 7 of [31]).
Later, Gonchar [31] studied the rate of attraction of each individual pole of F in B R m (F) to zeros of Q n,m . He introduced several indicators describing the asymptotic behavior of the zeros of Q n,m to a point λ ∈ C \ {0}. Set |z − w| 1 := min{1, |z − w|}, z, w ∈ C.
The following theorem (see [31], Theorem 1) asserts that a pole of F of order ν in B R m (F) attracts with geometric rate exactly ν zeros of the polynomials Q n,m . Theorem 2 (Gonchar's theorem). Let F be defined as in Equation (1), m ∈ N, and let 0 = λ ∈ C. The following two assertions are equivalent: (a) λ ∈ B R m (F) and F has a pole at λ.
If either (a) or (b) takes place, then where ν is the order of the pole at λ. Moreover, if we assume further that lim inf n→∞ |λ − λ n,ν+1 | > 0, then Naively, Gonchar asked what we can say about λ if λ attracts a certain numbers of zeros of the polynomials Q n,m without such geometric rate in the above theorem. In [31], Gonchar also proposed the following conjecture.
Theorem 4 (Suetin's theorem). Assume that the formal power series F(z) = ∑ ∞ k=0 f k z k has coefficients such that for fixed m ∈ N and sufficiently large n ∈ N the approximants R n,m have precisely m finite poles λ n,1 , . . . , λ n,m , which are convergent: lim Then, (i) The power series defines a holomorphic function F in the disk B R min , where R min := min 1≤j≤m |λ j |.
(iii) All points λ 1 , . . . , λ m are singularities of F, the ones lying in the disk B R m−1 (F) are poles, and F has no other poles in this disk.
It is easy to check that if the Taylor coefficients of F(z) = ∑ ∞ k=0 f k z k satisfy f k = 0 and f k+1 = 0, then f k / f k+1 is a zero of Q n,1 . Therefore, when m = 1, Gonchar's conjecture and Suetin's theorem reduce to the following Fabry ratio theorem (see [36]).
Theorem 5 (Fabry ratio theorem). Suppose that the coefficients of a power series F(z) = ∑ ∞ k=0 f k z k are such that the limit exists. Then, the series converges uniformly on each compact subset of the disk B |λ| and λ is a singular point of the function F.
Recently, extensions of Montessus de Ballore-Gonchar's theorem to Type II Hermite-Padé approximation and its generalizations were proved. As time progresses, there were many results on Type II Hermite-Padé approximation in the direction of Gonchar's conjecture. The purpose of this survey paper is to review all those results and collect open problems in this respect.

Definition and Notation
Type II Hermite-Padé approximation involves the approximation of a vector of functions by a vector of rational functions with the same denominator. Let F = (F 1 , F 2 , . . . , F d ) be a system of d formal Taylor expansions at the origin; that is, for each = 1, 2, . . . , d, we have In what follows, N := {1, 2, 3, . . .}, P n is the set of all polynomials of degree at most n, and for a given multi-index m = (m 1 , m 2 , . . . , m d ) ∈ N d , we define |m| := m 1 + m 2 + · · · + m d .
. . , F d ) be a system of d formal Taylor expansions as in Equation (5) and m = (m 1 , m 2 , . . . , m d ) ∈ N d be a fixed multi-index. Then, for each n ≥ max{m 1 , m 2 , . . . , m d }, there exist Q ∈ P |m| and P ∈ P n−m for all = 1, 2, . . . , d such that Q ≡ 0 and The vector of rational functions R n,m = (P 1 /Q, P 2 /Q, . . . , P d /Q) is called an (n, m) Type II Hermite-Padé (HP) approximant of F.
Traditionally, the numbers of interpolation conditions of F , for = 1, 2, . . . , d, in Equation (6) at 0 are selected to be the same, which is n + 1. Since m = 0 does not provide any advantage in locating singularities of F , we can restrict our multi-indices m ∈ N d as stated in Definition 1. For any fixed (n, m) ∈ N × N d , in general, R n,m may not be unique and we assume that, given (n, m), one particular solution is assigned. We write R n,m = (R n,m,1 , R n,m,2 , . . . , R n,m,d ) = (P n,m,1 , P n,m,2 , . . . , P n,m,d )/Q n,m , where P n,m, , = 1, 2, . . . , d, and Q n,m are chosen so that Q n,m is a monic polynomial that has no common zero with all the P n,m, . When m remains fixed, we call the sequences {R n,m } n≥|m| m row sequences. When n and m have the relation, m 1 = · · · = m d = m, n = (d + 1)m, m ∈ N (or nearby configurations of multi-indices), the sequences {R n,m } are called diagonal sequences. Another construction called Type I HP approximants has very close algebraic relation to Type II HP approximants. For Type I HP approximation, one wants to approximate polynomial combination of the vector (F 1 , F 2 , . . . , F d ) by a polynomial. However, in Type II HP approximation, one wants to approximate F separately by rational functions with the same denominator. For interested readers, we refer to Chapter 4 of [37] for more details about Type I HP approximation. Since we only provide the survey of the studies of Type II HP approximants, we omit the word "Type II" when we refer to Type II HP approximants. Moreover, we would like to emphasize that, for d = 1, the HP approximant R n,m reduces to the (n − m, m) classical Padé approximant. However, because we are interested in the cases when m is fixed and n → ∞, all theorems in Section 1 hold for the scalar HP approximation.

Results and Conjectures
In the direction of row sequences, the paper [38] by Graves-Morris and Saff was a pioneering result in the sense that it initiated a convergence theory for HP approximants. In [38], they proved a Montessus de Ballore type theorem for HP approximants under the concept of polewise independence. Later, Cacoq, de la Calle, and López [39] improved the results [38] in several ways; namely, improving the estimates on the rates of convergences of {R n,m } n≥|m| and {Q n,m } n≥|m| and weakening the assumption of polewise independence. Note that in the Montessus de Ballore type theorem for HP approximants in [39], they found the exact rate of convergence of {R n,m } n≥|m| but not of {Q n,m } n≥|m| . A significant contribution in the direction of row sequences is due to Cacoq, de la Calle, and López [40] where they gave necessary and sufficient conditions for the convergence with geometric rate of {Q n,m } n≥|m| and calculated the exact rate of convergence of {Q n,m } n≥|m| . To explain the results in [40], we need to state some definitions.
Let R 0 (F) be the radius of the largest disk B R 0 (F) to which all the expansions F , = 1, 2, . . . , d can be extended analytically. If R 0 (F) = 0, we take B R m (F) = ∅, m ≥ 0; otherwise, R m (F) is the radius of the largest disk B R m (F) centered at the origin to which all the analytic elements (F , B R 0 (F ) ) can be extended so that F has at most m poles counting multiplicities.
To prove an analog of Montessus de Ballore-Gonchar's theorem for HP approximants, we have to decide what actually is the limit of {Q n,m } n≥|m| . We found that the convergence of {Q n,m } n≥|m| for the vector case is more complicated than the one for the scalar case, as the following example shows.
and fix m = (1, 1). DefineF := (F 1 ,F 2 ), wherê It is not difficult to check that HP approximants of (F, m) and (F, m) have the same Q n,m . Applying Theorem 4.4 in [39] forF, 2 and 4 which are poles of F attract zeros of Q n,m with geometric rate as n → ∞. To be precise, we obtain lim sup n→∞ Q n,m − (z − 2)(z − 4) 1/n < 1.
However, since both F 1 and F 2 are meromorphic up to B 3 , 4 is not a pole of F. Based on the above example, the idea of pole in Definition 2 is not suitable when we study a system of functions. The authors of [40] proposed a new definition of pole (called "system pole") below. Definition 3. Given F = (F 1 , F 2 , . . . , F d ) in as Equation (5) and m = (m 1 , m 2 , . . . , m d ) ∈ N d , we say that λ ∈ C \ {0} is a system pole of order τ of F with respect to m if τ is the largest positive integer such that for each t = 1, 2, . . . , τ, there exists at least one polynomial combination of the form which is holomorphic on a neighborhood of B |λ| except for a pole at z = λ of exact order t.
Notice that in Definition 3, instead of examining poles of each function F separately, we examine poles of the polynomial combinations of the vector (F 1 , F 2 , . . . , F d ). When d = 1, the statement F has a system pole of order τ with respect to m = m simply means that F has a pole of order τ in B R m (F) . From Example 1, 2 and 4 are system poles of order 1 of F with respect to m but 4 is not a pole of F. Conversely, the following example shows that a pole of F may not be a system pole.
and fix m = (1, 1). Clearly, 2 is a pole of F. Since 1 and 3 are only system poles of F, 2 is not a system pole of F.
In conclusion, a system pole may not be a pole of F or vice versa. Note that a system F cannot have more than |m| system poles with respect to m counting their order (see Lemma 3.5 of [40]).
To state the main result (see Theorem 6 below) in [40], we need a generalization of the notion R m (F) for a system of functions F. For each system pole λ of F with respect to m, we want to define a corresponding characteristic number R λ (F, m) as follows. Let τ be the order of λ as a system pole of F. For each t = 1, . . . , τ, denote by r λ,t (F, m) the largest of all the numbers R t (g) (the radius of the largest disk containing at most t poles of g), where g is a polynomial combination of type Equation (8) that is analytic on a neighborhood of B |λ| except for a pole at z = λ of order t. Then, we define As in the scalar case, we denote by Q m (F) the monic polynomial whose zeros are the system poles of F with respect to m taking account of their order. The set of distinct zeros of Q m (F) is denoted by P (F, m).
An analog of Montessus de Ballore-Gonchar's theorem for HP approximation (see Theorem 1.4 of [40]) is stated below. Theorem 6. Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m ∈ N d . Then, the following two assertions are equivalent: (a) R 0 (F) > 0 and F has exactly |m| system poles with respect to m counting multiplicities.
The denominators Q n,m , n ≥ |m|, of the HP approximants of F are uniquely determined for all sufficiently large n, and there exists a polynomial Q |m| of degree |m|, Q |m| (0) = 0, such that where · denotes the coefficient norm in the space of polynomials.
Moreover, if either (a) or (b) takes place, then Q |m| ≡ Q m (F) and An exact expression for the rate and region of convergence of {R n,m } n≥|m| to F is given in Theorem 3.7 of [40]. Theorem 6 is considered to be a cornerstone in the study of HP approximation on row sequences in the sense that it contains the first inverse type result.

Conjecture 1.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m ∈ N d . If the denominators Q n,m are uniquely determined for all sufficiently large n and σ(λ) ≥ 1, then λ is a system pole of F with respect to m of order τ = σ(λ).
The generalization of "singularity" when we consider a system of functions is stated below. (5) and m = (m 1 , m 2 , . . . , m d ) ∈ N d , we say that λ ∈ C \ {0} is a system singularity of F with respect to m if there exists at least one polynomial combination G of the form in Equation (8) analytic on B |λ| and λ is a singular point of G.

Conjecture 2.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m ∈ N d . If the denominators Q n,m are uniquely determined for all sufficiently large n and γ(λ) ≥ 1, then λ is a system singularity of F with respect to m.
In ( [41], p. 155), López and Zaldivar gave an example to support their conjectures. That example shows that a point in C \ {0} can be a system pole or a system singularity depending how it attracts zeros of Q n,m . For the self-contained purpose, we show their example here.
Conjecture 1 for the scalar case (d = 1) is the part "(b)⇒(a)" in Gonchar's theorem. However, Conjecture 2 for the scalar case is Gonchar's conjecture which remains open. The converse statement of Conjecture 1 proved in Theorem 2.1 of [41] is stated as follows.
Theorem 7. Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m ∈ N d . Assume that λ is a system pole of order τ of F with respect to m. Then, and According to Gonchar's theorem, a natural conjecture is Conjecture 3. Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m ∈ N d . Assume that λ is a system pole of order τ of F with respect to m. Then, the inequalities in Equation (10) are equalities. Moreover, if we assume further that lim inf n→∞ |λ − λ n,τ+1 | > 0, then In the past few years, some progress (see [41][42][43]) on Conjecture 2 was made. In those papers, the authors studied Conjecture 2 with an additional assumption, namely, Without loss of generality, let λ 1 , λ 2 , . . . , λ |m| be the zeros of Q |m| in such a way that 0 < |λ 1 | ≤ |λ 2 | ≤ · · · ≤ |λ |m| | (12) and λ n,1 , λ n,2 , . . . , λ n,|m| be the zeros of Q n,m indexed so that lim n→∞ λ n,k = λ k , k = 1, 2, . . . , |m|.
The latest progress on Conjecture 2 with the condition in Equation (11) is the following theorem (see Theorem 4.4 of [42]). Theorem 8. Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m ∈ N d . Assume that Equations (11) and (12) hold. Assume further that, if λ j , λ j+1 , . . . , λ j+N−1 , where N > 1, are the zeros of Q |m| lying on the same circle (namely, they have the same modulus), then Then, each λ k , k = 1, 2, . . . , |m|, is a system singularity of F with respect to m. Moreover, if λ k is a zero of multiplicity τ k of Q |m| which verifies Equation (13), then it is a system pole of F with respect to m of order τ k .
Considering the coefficients of Q n,m as the coefficients of a certain recurrence relation, López and Gerpe [42] employed Buslaev's results (extensions of Poincaré's theorem on recursion relations) in Theorems 1 and 2 of [44] to prove the above theorem.
Note that the following special case of Conjecture 2 (an analog of Suetin's theorem) remains unsolved.

Conjecture 4.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m ∈ N d . Assume that Q n,m is unique for all sufficiently large n, Equation (11) takes place, and Q |m| (λ) = 0. Then, λ is a system singularity of F with respect to m.

Some Remarks
• The statement (b)⇒(a) in Theorem 6 was the first inverse-type result on the study of HP approximation on row sequences. Its proof is very constructive. It suggests how to find a polynomial combination in Equation (8) verifying that all zeros of Q |m| in Equation (9) are the system poles of F with respect to m.

•
The study of zeros of P n,m, , = 1, 2, . . . , d, in Equation (7) is irrelevant to our interest in this paper. However, it is worth mentioning the paper [45] where la Calle Ysern and Mínguez Ceniceros studied the distribution of zeros of P n,m, , = 1, 2, . . . , d, as m is fixed and n → ∞.

Definitions and Notation
After an analog of Montessus de Ballore-Gonchar's theorem for HP approximation (Theorem 6) was proved, Theorem 6 was extended for various generalizations of HP approximation. These generalizations are formulated to approximate a vector of functions holomorphic on a neighborhood of the following sets E.
Let E be an infinite compact subset of the complex plane C such that C \ E is simply connected. Whenever we consider set E, E is described as above, unless we specifically say otherwise. Denote by H(E) the space of all functions holomorphic in some neighborhood of E. We define In Section 3, we are mainly interested in approximating a vector in the space H(E) d .

Orthogonal Hermite-Padé Approximations
The first two approximations are constructed from orthogonal polynomials on E. Let µ be a finite positive Borel measure with infinite support supp(µ) contained in E. We write µ ∈ M(E) and define the associated inner product, Using the Gram-Schmidt process, we can generate the sequence of the orthonormal polynomials with positive leading coefficients p n (z) := κ n z n + · · · , κ n > 0, n = 0, 1, 2, . . . , corresponding to the inner product in Equation (14). It is well-known that, for each n ≥ 0, such orthonormal polynomial p n is unique. Combining the concepts of HP approximation and orthogonal polynomials, we define two types of orthogonal Hermite-Padé approximations. The first one is a natural extension of HP approximation. When E = {z ∈ C : |z| ≤ 1} and µ is the normalized arc length on the unit circle, p n (z) = z n for all n ≥ 0 and the system of linear equations in Equation (15) reduces to . . , n and = 1, 2, . . . , d.
Then, the above polynomial Q µ,S n,m coincides with the polynomial Q in Equation (6), which further implies that, for each = 1, 2, . . . , d, the corresponding polynomial ∑ n−m j=0 Q µ,S n,m F , p j µ p j coincides with the polynomial P in Equation (6). Moreover, when d = 1, the rational function R µ,S n,m is the usual (n − m, m) orthogonal Padé approximant defined in [46]. Therefore, the approximation in Definition 5 is a natural generalization of orthogonal Padé approximation to the vector case.
The definition of the other orthogonal Hermite-Padé approximation was recently introduced by Bosuwan and López [47] in the problem concerning localization of system poles of F ∈ H(E) d around the set E.
The corresponding vector of rational functions The motivation of the above definition is in the paragraph that contains Equation (39).

Faber-Hermite-Padé Approximations
Now, we want to combine the ideas of HP approximation and Faber polynomials on E. Let Φ be the unique Riemann mapping function from C \ E to the exterior of the closed unit disk verifying Φ(∞) = ∞ and Φ (∞) > 0. For each ρ > 1, the level curve of index ρ and the canonical domain of index ρ are defined by respectively. Let F ∈ H(E) d . Denote by ρ m (F) the index ρ > 1 of the largest canonical domain D ρ to which F has at most m poles counting multiplicities.
The Faber polynomial of E of degree n is One can also define Φ n as the polynomial part of the Laurent expansion of Φ n at infinity. The nth Faber coefficient of F ∈ H(E) with respect to Φ n is defined by the formula where ρ ∈ (1, ρ 0 (F)).
As with the SOHP and MOHP approximations, we have two ways to define Faber-Hermite-Padé approximations.
The corresponding vector of rational functions R E,S n,m := (R E,S n,m,1 , R E,S n,m,2 , . . . , R E,S n,m,d ) are called an (n, m) modified Faber-Hermite-Padé (MFHP) approximant of F with respect to E.
The above MFHP approximation was recently introduced by Bosuwan and López [49]. In [49], MFHP approximation was proved to be an effective tool to locate system poles of F ∈ H(E) d around the set E.
This approximation was introduced in [2] by Mahler long before other generalizations. Note that, if all the interpolation points are 0, then the MHP approximants reduce to HP approximants.

Some Remarks
n,m , Q E,S n,m , orQ E,M n,m is equivalent to solving |m| + 1 unknowns from linear system of |m| equations in Equations (15)-(18), respectively. Moreover, finding P α n,m, , = 1, 2, . . . , d and Q α n,m is equivalent to solving (n + 1)d + 1 unknowns from linear system of (n + 1)d equations in Equation (19). Therefore, these polynomials, Q µ,S n,m , Q µ,M n,m , Q E,S n,m , Q E,M n,m , and Q α n,m , always exist but may not be unique. Since such polynomials are not the zero function, we normalize them to be "monic" polynomials. Moreover, Q α n,m in Equation (20) is chosen so that it does not have a common zero with all the P α n,m, . We would like to emphasize that for any (n, m) ∈ N × N d , R µ,S n,m , R µ,M n,m , R E,S n,m , R E,M n,m , and R α n,m may not be unique.

•
Extensions of generalized HP approximations in Definitions 5-9 to a compact set E whose complement is connected are possible. However, the results in this survey paper are restricted to the case when E is a compact subset of the complex plane with simply connected complement in the extended complex plane. This is because, for the sets E which are disconnected, the zeros of the corresponding orthonormal polynomials p n , second type functions s n (defined in Equation (21)), or Faber polynomials Φ n may lie in C \ E which may be the locations of system poles.

Classes of Measures in M(E)
For the studies of SOHP and MOHP approximations, we need to impose asymptotic properties of the sequences of the orthonormal polynomials {p n } n≥0 and the corresponding second type functions {s n } n≥0 defined below We keep classes of measures in M(E) here so that they do not disturb the flow of our paper. The readers who are not interested in SOHP or MOHP approximations may skip Section 3.1.5.
where the c n s are positive constants and S is a non-vanishing holomorphic function on C \ E. (22) and (23) are assumed to hold uniformly on compact subsets of C \ E and the ones in Equations (24)- (26) are assumed to hold uniformly on compact subsets of C \ E.

The limits in Equations
The classes Reg 1 (E) and Reg 2 (E) are more or less the same in some cases. In particular, if E is convex, then Reg 1 (E) = Reg 2 (E) and these two classes coincide with the regular class in the usual sense (see Definition 3.1.2 of [50] for the definition of the regular class in the usual sense). Clearly, S(E) ⊂ Rat 1 (E) ⊂ Reg 1 (E) and Rat 2 (E) ⊂ Reg 2 (E).
Denote by K 1 the collection of all compact sets E (stated at the beginning of Section 3) satisfying the condition that the inverse of the exterior conformal function Φ −1 can be extended continuously to C \ {w ∈ C : |w| < 1}. If we assume that E ∈ K 1 and µ ∈ Rat 1 (E), then uniformly on compact subsets of C \ E (see Lemma 3.1 of [46]). This statement implies that, if E ∈ K 1 , (see, e.g., Lemma 2 of [51]). Therefore, for each m ∈ N, Rat 1 (E) ⊂ Reg m 1 (E) and Rat 1,2 (E) ⊂ Reg m 1,2 (E). For a general compact set E, the class of measures Reg 1 (E) has been well studied and characterized in terms of the analytic properties of the measure or of the corresponding sequence of leading coefficients (see, e.g., [50], Theorem 3.1.1). The situation is not quite the same for other classes of measures. To discuss them would take us too far from our main direction; rather we refer the reader to Pages 532-533 in [52] for more details and references.

Classes of Tables of Interpolation Points
When we state results on MHP approximation, we refer to two classes of tables of interpolation points. (z − α n,k ).
(a) α ∈ Strong(E) if and only if the corresponding polynomials a n satisfy the following strong asymptotics: α ∈ Root(E) if and only if the corresponding polynomials a n satisfy the following nth root asymptotics:

The Scalar Case
We add a section dedicated to the scalar case of all generalizations (defined in Section 3.1) here. The purpose is not only to help the readers consolidate their understanding of the vector case but also to discuss some interesting theorems and conjectures corresponding to such scalar case. When d = 1, we have m = m and write where Q m (F) is the monic polynomial whose zeros are the poles of F in D ρ m (F) taking account of their order.
where P m (F) is the set of the distinct zeros of Q m (F).
lim sup lim sup lim sup under the convention that if K ⊂ E, then Φ K is replaced by 1. (31) is equality when α is a Newton type.

Remark 1. The inequality in Equation
where σ n,m := inf r F − r E , and the infimum is taken over the class of all rational functions of type (n, m) r(z) = a n z n + a n−1 z n−1 + · · · + a 0 b m z m + b m−1 z m−1 + · · · + b 0 .
We refer the reader to [56] for the proof of the fifth equality in Equation (32). For Gonchar's conjecture, we made good progress in the case when the denominators Q The proof of the above theorem relies on deep results on refinements of Poincaré's theorem on recurrence relations developed by Buslaev in Theorems 5 and 6 of [57]. These Buslaev's results connect Suetin's theorem (for classical Padé approximants) and analogous ones for generalized Padé approximants (see Theorem 3 of [57] and Theorems 2.7 and 2.8 of [58]).
An analog of Gonchar's theorem for Q α n.m when α is a Newton type was proved by himself in Theorem 2 of [31]. For other generalizations, we have proofs only for an analog of (a)⇒(b) in Gonchar's theorem (see Corollary 5.6 of [51], Corollary 1 of [59], Theorems 2.1 and 2.3 of [58]).
A direct analog of the Fabry ratio theorem for orthogonal and Faber polynomial expansions is Theorem 11. Let F ∈ H(E). If one of the following holds: then Φ −1 (ξ) is a singularity of F and ρ 0 (F) = |ξ|.
The part of the above theorem concerning the limit of the ratio of Faber coefficients follows from a straightforward change of variables: and the use of the Fabry ratio theorem for the analytic part of F • Φ −1 . The part of the above theorem concerning the limit of the ratio of Fourier coefficients is in Theorem 2.3 of [46]. We end this section by stating an application of Theorem 11 to orthogonal polynomial theory on the unit circle. Firstly, we recall some facts. If E = B 1 and the measure µ supported on the unit circle satisfies the Szegő condition, where dµ(θ) = w(θ)dθ/2π + dµ s (θ) is the Radon-Nikodym decomposition of µ, then it is well known that µ ∈ S(E), the leading coefficients of the orthonormal polynomials p n satisfy lim n→∞ κ n = κ := exp − 1 4π 2π 0 log w(θ)dθ , denotes the interior Szegő function (see ( [60], pp. [19][20] for the proof). Therefore, Theorem 11 enables us to locate the singularity nearest the origin of the reciprocal of the interior Szegő function 1/S int in terms of the Verblunsky coefficients α n := −p n (0)/κ n .

Theorem 12.
Let µ satisfy the Szegő condition in Equation (33) and assume that 1/S int ∈ H(B). Suppose that the Verblunsky coefficients α n corresponding to µ verify Then, λ is a singularity of 1/S int and 1/S int is holomorphic on B |λ| .
The above result was stated in Corollary 2.4 of [46].

The Vector Case
Graves-Morris and Saff were the first to prove an extension of Montessus de Ballore's theorem to MHP approximation in Theorem 3 of [38]. Their theorem relies on the concept of polewise independence, which influenced the current author to define the following adapted polewise independence.
Note that, if we replace Φ in (ii) of Definition 12 by the identity mapping, then the above definition reduces the definition of polewise independence in [38].
Putting together Theorem 2.3 in [51] and Theorem 1 in [61], we obtain the following analog of Montessus de Ballore's theorem for SOHP and SFHP approximations under the concept of polewise independence in Definition 12.
By Q F |m| , we denote the monic polynomial whose zeros are the poles of F in D ρ |m| (F) taking account of their order. The set of distinct zeros of Q F |m| is denoted by P |m| (F). Moreover, In fact, Cacoq and López [62] were the first who proved the above theorem but in the context of SOHP approximation for the case when E is the closed unit disk and µ is supported on the unit circle.
Making use of incomplete orthogonal Padé and incomplete Padé-Faber approximations (see Definition 5.1 of [51] and Definition 5 of [59]), Bosuwan proved another Montessus de Ballore type theorem for SOHP and SFHP approximations in Theorem 2.4 of [51] and Theorem 1 of [59]. The idea of incomplete approximants allows us to analyze each function F individually. Let us define some more notation about the region of convergence.
Given a system F = (F 1 , F 2 , . . . , F d ) and a multi-index m = (m 1 , m 2 , . . . , By Q F m , we denote the monic polynomial whose zeros are the poles of F in D m (F) counting multiplicities. This set of poles is denoted by P m (F). For = 1, 2, . . . , d, set P m, (F) := P m (F) ∩ D ρ m (F ) .
To each pole λ of F in this system of domains we associate an index (λ) ∈ {1, 2, . . . , d} as follows. The index (λ) verifies that λ ∈ D ρ m (λ) (F (λ) ) and λ is a pole of F (λ) of the same order as is a pole of F in D m (F). If there are several indices satisfying this condition, then we choose one among those with greatest ρ m (F ).

Theorem 14.
Let m ∈ N d be a fixed multi-index, F = (F 1 , F 2 , . . . , F d ) ∈ H(E) d , P m (F) = {λ 1 , λ 2 , . . . , λ q }, and µ ∈ Rat 1,2 (E). Suppose that F ∈ H(E) d has exactly |m| poles in D m (F). Then, R µ,S n,m and R E,S n,m is uniquely determined for all sufficiently large n and for each = 1, 2, . . . , d, R µ,S n,m, and R E,S n,m, converge uniformly to F on compact subsets of D ρ m (F ) \ P m, (F). Moreover, for all = 1, 2, . . . , d, n,m, where K is any compact subset of D ρ m (F ) \ P m, (F) and if K ⊂ E, then Φ K is replaced by 1. Additionally, we have Note that Theorems 13 and 14 have their own values. Section 3 in [51] and Section 2 in [59] gave examples showing when Theorem 13 is applicable but Theorem 14 is not applicable or vice versa.
During my visit to Universidad Carlos III de Madrid in the summer of 2016, López and I initiated a project emphasizing the study of analogs of the inverse statement of Montessus de Ballore-Gonchar's Theorem for SOHP and SFHP approximations. As guided by the principle of the definition of system pole in [40] (see Definition 3), we define Definition 13. Given F = (F 1 , F 2 , . . . , F d ) ∈ H(E) d and m = (m 1 , m 2 , . . . , m d ) ∈ N d , we say that λ ∈ C is a system pole of order τ of F with respect to m if τ is the largest positive integer such that, for each t = 1, 2, . . . , τ, there exists at least one polynomial combination of the form which is holomorphic in a neighborhood of D |Φ(λ)| except for a pole at z = λ of exact order t.
We had an application in my mind that we wanted to use the zeros of {Q µ,S n,m } n≥|m| and {Q E,S n,m } n≥|m| to detect |m| system poles of F ∈ H(E) d nearest the set E. Although we are not able to solve our initial aim (Conjecture 6 below), we observed that the system poles of F with respect to m and their orders are the same as the system poles of with respect to m = (1, 1, . . . , 1), with |m| = |m| and their orders. The creation of the vector in Equation (39) [47,49]).
To state main results in [47,49,55], we need a generalization of R λ (F, m). For each system pole λ of F with respect to m, we define a characteristic index ρ λ (F, m) as follows. Let τ be the order of λ as a system pole of F. For each t = 1, 2, . . . , τ, denote by ρ λ,t (F, m) the largest of all the numbers ρ t (G) (the index of the largest canonical domain containing at most t poles of G), where G is a polynomial combination of the type in Equation (38) that is holomorphic in a neighborhood of D |Φ(λ)| except for a pole at z = λ of order t. There is only a finite number of such possible values so the maximum is indeed attained. Then, we define ρ λ (F, m) := min t=1,2,...,τ ρ λ,t (F, m).

(d)
The polynomials Q α n,m are uniquely determined for all sufficiently large n, and there exists a polynomialQ |m| of degree |m| such that lim sup n→∞ Q α n,m −Q |m| 1/n =θ < 1.
Moreover, if one of Assertions (a)-(d) takes place, thenQ |m| =Q |m| =Q |m| = Q m (F) and   are much more complicated than HP approximants. There is no nice equality (similar to Equality (3.3) in [41]) which allows us to say something about the convergence and divergence of incomplete Padé approximation appear in SOHP, MOHP, SFHP, and MFHP approximations. Such equality is also a main ingredient for the proofs of Gonchar's theorem and Suetin's theorem for classical Padé approximation stated in Section 1 and an analog of Gonchar's theorem for the scalar MHP approximation when α is a Newton type in Theorem 2 of [31].

Conclusions
The study of classical Padé approximants on row sequences roots from the work of Montessus de Ballore [29] on the uniform convergence of row sequences of the approximants. The subject received renewed interest by Gonchar [30], who proved the converse statement of Montessus de Ballore's theorem. Soon after that, he studied an attraction of an individual pole of F in B R m (F) to the zeros of Q n,m when m is fixed as n → ∞ in [31]. In the same paper, Gonchar proposed his conjecture (which is commonly called Gonchar's conjecture and remains unsolved). Since its introduction, several positive answers have supported the conjecture (see, e.g., [32][33][34][35]).
Graves-Morris and Saff [38] extended Montessus de Ballore's theorem to the vector case (which we call HP approximation). Recently, Cacoq, de la Calle Ysern, and López [40] defined the notion of system pole and proved the inverse result of Graves-Morris and Saff's result. In the spirit of Gonchar's conjecture, López and Gerpe proposed the following conjectures.

Conjecture 1.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m ∈ N d . If the denominators Q n,m are uniquely determined for all sufficiently large n and σ(λ) ≥ 1, then λ is a system pole of F with respect to m of order τ = σ(λ).

Conjecture 2.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m ∈ N d . If the denominators Q n,m are uniquely determined for all sufficiently large n and γ(λ) ≥ 1, then λ is a system singularity of F with respect to m.
Some progress on Conjecture 2 was made in [41][42][43]. However, both conjectures are still open. Concurrently, generalizations of HP approximation were introduced and studied on row sequences [47,49,51,55,61]. Both conjectures can be asked for all generalizations of HP approximation. Although SOHP and SFHP approximations are natural, the proofs of analogs of Montessus de Ballore-Gonchar's theorem for such approximations are not available (see Conjecture 6 for our expectation).
Analogs of Montessus de Ballore-Gonchar's theorem and Suetin's theorem for the scalar case of all generalizations were completely proved. However, the exact rates of convergences of