On Semi-Classical Orthogonal Polynomials Associated with a Modiﬁed Sextic Freud-Type Weight

: Polynomials that are orthogonal with respect to a perturbation of the Freud weight function by some parameter, known to be modiﬁed Freudian orthogonal polynomials, are considered. In this contribution, we investigate certain properties of semi-classical modiﬁed Freud-type polynomials in which their corresponding semi-classical weight function is a more general deformation of the classical scaled sextic Freud weight | x | α exp ( − cx 6 ) , c > 0, α > − 1. Certain characterizing properties of these polynomials such as moments, recurrence coefﬁcients, holonomic equations that they satisfy, and certain non-linear differential-recurrence equations satisﬁed by the recurrence coefﬁcients, using compatibility conditions for ladder operators for these orthogonal polynomials, are investigated. Differential-difference equations were also obtained via Shohat’s quasi-orthogonality approach and also second-order linear ODEs (with rational coefﬁcients) satisﬁed by these polynomials. Modiﬁed Freudian polynomials can also be obtained via Chihara’s symmetrization process from the generalized Airy-type polynomials. The obtained linear differential equation plays an essential role in the electrostatic interpretation for the distribution of zeros of the corresponding Freudian polynomials. the following result for recurrence coefﬁcients for the modiﬁed Freudian weight in (4).


Introduction
The theory of orthogonal polynomials on finite intervals, which was studied in the seminal works of G. Szegö [1], is essentially not the same as the theory of orthogonal polynomials on large intervals. In the late 1970s, Géza Freud has done prominent work on polynomials orthogonal with respect to exponential weights on the real line, which referred to as Freud-type orthogonal polynomials [2]. At that time, the aim of Freud's work was to extend the theory of best approximation and Jackson-Bernstein type estimates to the real line; by expecting that orthogonal expansion may serve as a near-best approximation (cf. [2][3][4][5]).
Classical orthogonal polynomials have their weight function that satisfies Pearson's equation where σ(x) is a monic polynomial of degree at most 2 and τ(x) is a polynomial with degree 1. However for semi-classical orthogonal polynomials, the weight function w(x) satisfies the Pearson Equation (1) with either deg(σ) > 2 or deg(τ) = 1 (cf. [6][7][8][9][10]). The information we presently know about classical orthogonal polynomials on infinite intervals is largely due to the presence of generating functions, differential equations and recursive formulas that can be exploited to discover certain properties of such polynomials. For instance, Hermite, Laguerre, Pollaczek, Sonine, Stieltjes-Wigert, Poisson-Charlier, Krawtchouk, Meixner, Lommel, Askey-Wilson, and Al-Salam-Carlitz polynomials are some of the classical examples of orthogonal polynomials with non-compact support. On the other hand, the study of semi-classical orthogonal polynomials associated with various deformed weights is still active research since some of the properties are not yet explicitly studied. For some work on semi-classical weights, we refer the reader to [11][12][13][14][15][16]. Géza Freud (1973) considered class of weight functions where K is a positive constant [5,10,16]. Certain properties of Freud-type polynomials can be found in [12,15,[17][18][19]. The authors in [12] considered monic orthogonal polynomials with respect to the quartic Freud-type weight function; that is, (m = 2 in Equation (2)) upon measure deformation by exp tx 2 . Some approximation-theoretic properties of orthonormal Freud-type polynomials is given in [19]. It is noted in [20] that slight generalizations of the weights of some of the classical orthogonal polynomials present great difficulties in exploring their properties, most of which have not been conquered yet.
Following the work in [21], a new semi-classical deformation for semi-classical Freudian weight function has been essential because of their practical applications and connections to random matrix theory, integrable systems and Painlevé equations. Motivated by the above facts, we investigate certain properties of monic orthogonal polynomials associated with the modified Freud-type inner product with real parameters λ > 0, c > 0, t ∈ [0, K], (K ∈ R + ). The modified Freudian weight in (3) emerges from deformation of the scaled Freudian weight exp(−cx 6 ), c > 0. (see also [3,5,12]). The motive for the choice of the orthogonality weight in Equation (3) is as follows. First, from some of the well-known classical orthogonal polynomials (for e.g., Laguerre polynomials), a new class of semi-classical orthogonal polynomials can be obtained by making slight modifications of the orthogonality measure; for example, the semi-classical Laguerre measure, which is well-studied in [11,[22][23][24]. Knowing the fact that such measure modification usually leads to difficulties as mentioned in the papers of P. Nevai [5,20]; and motivated by the work in [20], a slight modification of a new orthogonality measure on non-compact support always leads to a new class of orthogonal polynomials provided the corresponding moments are finite [1,10,25]. We noticed that the generic semi-classical weights in [8] are special instances of such measure modification given in [20]. Secondly, perturbation of an orthogonality measure allows one to investigate certain fresh properties such as new-type of Toda-like differential-recurrence relations, new non-linear higher-order recursions as well as differential equations, and some properties of the zeros. For more on this, one can refer to [26,27] and the references therein. The obtained non-linear differential/difference equations for such orthogonal polynomials have considerable applications in modeling non-linear phenomena, Soliton Theory, Random matrix theory, Quantum oscillators and in the crystal structure in solid-state physics (for e.g., see the works by Chen and Its [13], Clarkson and Jordaan [11], Van Assche [10], Marcellán [28] and others).
Our aim in this paper is to investigate certain results on the recurrence coefficients for the semi-classical modified Freudian polynomials and to explore some properties for these polynomials including differential, differential-recurrence and difference equations and the zeros. For the semi-classical weight under consideration, it is not easy to obtain a concise formulation of the recurrence coefficients straightforwardly (see [11,12,15,29] for some related works). Related to the weight in (3), we accomplish some new work, by investigating Toda-type relations, differential, and differential-recurrence relations for the recurrence coefficients and the polynomials themselves.
The paper is organized as follows. In Section 2, we give some background for semi-classical modified Freud-type orthogonal polynomials. Section 3 gives some results that the recurrence coefficients for the weight in (3) obey nonlinear recursion relations.
Besides, nonlinear differential-recurrence relations as well as differential equations satisfied by the recurrence coefficients as well as the polynomials are also obtained. Section 4 employs both the methods of ladder operators as well as Shohat's quasi-orthogonality approach to derive the ladder relations associated with the modified weight in (3) and note that the results obtained by both methods are the same. In Section 5, we show that the modified Freudian polynomials can be obtained from the semi-classical generalized Airy type polynomials using Chihara's symmetrization process [30]. Section 6 gives an electrostatic interpretation of the zeros of the modified Freudian polynomials from the obtained differential equation. Section 7 provides the discussion of the results and the final section ends with conclusions.

The Modified Sextic Freud-Type Weight
Consider {S n (x; t)} n≥0 be sequence of monic polynomials that are orthogonal with respect to the semi-classical weight function with real parameters λ > 0, c > 0, and the orthogonality condition is given by where δ mn denotes the Kronecker delta, andζ n is the square of L 2 -norm of the polynomial S n (x, t).
Then the three-term recurrence relation follows as subject to the initial conditions S 0 (x; t) := 1 and β 0 S −1 (x; t) := 0. Multiplying both sides of Equation (6) by S n−1 (x; t)W λ (x; t) and integrating this with respect to x on R, which, due to the orthogonality condition (5), gives us One can see that S n (x; t) contains only the terms x n−j , j ≤ n and is even, since the weight function W λ (x; t) is even on R. This implies that S n (−x; t) = (−1) n S n (x; t) and S n (0; t) S n−1 (0; t) = 0.
The following proposition is about the finiteness of moments for the semi-classical weight (4).
By using the theory of integration [31], the higher order moments η k , k ∈ N are also finite. The even moments η 2n (t; λ), n ∈ N for the weight (4) are while the odd moments η 2n+1 (t; λ) are For t = 0, the moments for the symmetric weight w(x) = |x| 2λ+1 exp −cx 6 , c > 0 are given by However, finding the moments explicitly has not been an easy task for some semi-classical weights. We refer to an interesting work in [15], which shows that the moments for the modified sextic and dodic weights are expressed in terms of generalized hypergoemetric functions and for quartic Freudian weights, see [12] (see also [18]).
The Hankel determinant of moments η k for the symmetric weight in (4) is defined by [11,32] The following factorization for the determinants holds for the given weight in (4): where H where The following lemma shows the recurrence coefficients in terms of moment determinants for symmetric weights in general. See the recent work by Clarkson and Jordaan [15] for its proof and the references therein. m be the moment determinants as given in Equation (13), then the recurrence coefficient β n corresponding to the semi-classical weight (4) is given by Since the modified Freud-type weight in (4) depends on the parameter t, which means that the recurrence coefficients, Hankel determinants as well as the polynomials themselves also rely on t. However, unless it is mandatory we do not always state the t-dependence. For the Hankel determinantal characterization of the recurrence coefficients in terms of tau functions, we refer to the recent work in [15] (see also [10] and the references therein).

Explicit Formulation of the First Few Polynomials
In view of Equation (6), the first few monic polynomials are given by The following proposition, with its proof given for completeness, provides a concise formulation for modified Freudian polynomials S n (x; t). For a similar result, see also [32][33][34].
Proof. From the fact that the Freudian polynomials S n (x; t) are symmetric and monic of degree n; i.e., we obtain, for a fixed t ∈ R, where e n−2k = Ψ k (n) with Ψ 0 (n) = 1 and Ψ k (n) = 0 for k > n 2 . Substituting Equation (14) into the recursion relation Equation (6) and comparing the coefficients yields By applying induction on k, we prove Equation (15). For k = 1, we see that By using a telescopic sum on Equation (18), we obtain β j 1 (t; λ), for every n ≥ 1.
Suppose that Equation (15) is true for values up to k − 1 for every n ∈ N, i.e., Now, iterating Equation (17), we obtain Substituting Equation (19) into Equation (20) yields Equation (15). Thus, the result is valid for k ∈ N.
An alternate representation for Proposition 2 is given in the result below.  The normalization constantζ n in Equation (5) for the semi-classical weight (4) is given bŷ where Ψ k (n) is given in Equation (15).

Recurrence Coefficients of the Sextic Freud-Type Polynomials
In this section, we derive a non-linear discrete equation satisfied by the recursion coefficient β n (t; λ) for the semi-classical weight (4), which is one of the main results of this paper.
Theorem 4. The recurrence coefficient β n (t; λ) in Equation (6) satisfies the non-linear difference equation where β 0 = 0 and β 1 (t; λ) is given by where the expressions Q n and Ω n are, respectively, given by and Proof.
(i) By employing Freud's method in [16,20], we consider, for t ∈ R + , where {S n (x; t)} n≥0 are sequences of monic modified Freudian polynomials andζ n is given in Equation (21). Then we have where R n−2 ∈ P n−2 . On the other hand, by employing technique of integration on Equation (25), we obtain only consists of even powers of x and hence will dominate the limit as x → ±∞.
Since the Freudian measure given in (4) is symmetric, the integral expression when n is even. For n is odd, the expression S n (x; t) x is a polynomial of degree n − 1 and hence where Ω n is given in Equation (24). Using the recurrence iterated relation from Equation (6), we do have where Q n is given in Equation (23). Applying the identities in Equation (29) and the Pearson equation for the weight (4) together with Equation (28) into Equation (27), we obtain which simplifies to from the fact thatζ n−1 = 0 and Ω n is given in Equation (24). Note that Equations (30) and (26) yield Equation (22). (31) can be given in [10,12,17,26].

Remark 5. Quite similar non-linear discrete equations like Equation
We obtain the following result for recurrence coefficients for the modified Freudian weight in (4).

Theorem 6.
The recurrence coefficients β n (t; λ) obey the differential-difference formula where the expression Q n is given by Equation (23).
Proof. By differentiating the normalization constantζ n (t) with respect to t, we have where the first integral vanishes by orthogonality since is a monic polynomial in x of degree n − 1. By using orthogonality and the recurrence relation in Equation (6), we rewrite Equation (33) On the other hand, applying differentiation on both sides of Equation (7) with respect to t gives and substituting Equation (34) into above equation, we obtain the desired result. (6) is also known to be non-linear discrete mKdV-like equation (see [26] for similar equations).

Remark 7. The discrete equation in Equation
Our next result provides a higher order non-linear differential-recurrence relation satisfied by the recursion coefficients β n (t; λ) in reference to the semi-classical weight (4). (4), the recursion coefficients β n (t; λ) satisfy the following differential-recurrence equation

Theorem 8. For the modified Freud-type weight in
where Ω n and Q n are given in Equations (24) and (23), respectively and the expression ϑ(t) is given by Proof. The above result is obtained by applying differentiation of Equation (32) with respect to t; i.e., together with employing the following identity (22) where γ := 2λ + 1 and Ω n is given in (24) and the expression ϑ(t) is given in (35b).

The Differential-Difference Equation Satisfied by Modified Sextic Freud-Type Polynomials
In view of an extension of the ladder operators technique given in [35], Chen and Feigin [13] examined the weight w(x)|x − t| K , x, t, K ∈ R, for any smooth reference weight w(x). They indicated that when w(x) is the Gaussian (Hermite) weight (e −x 2 , x ∈ R), the recurrence coefficients fulfill a particular two-parameter Painlevé IV condition. Filipuk et al. [23] found that the recurrence coefficients for the quartic Freudian weight |x| 2α+1 e −x 4 +tx 2 , x, t ∈ R, α > −1 are identified with solutions of the Painlevé IV and the first discrete Painlevé equation. Clarkson et al. [12,15] gave a methodical investigation on Freud weights and some generalized work for [13] (see also [18,36]).
The following result gives the differential-difference equation (i.e., lowing operator) satisfied by the modified Freud-type polynomials corresponding to the weight given in Equation (4). This result works also for semi-classical polynomials Φ n (z) that are orthogonal with reference to a general symmetric weight in the semi-classical class, and specifically, it holds for the modified Freudian polynomials S n (x, t).

Lemma 2.
For α ≥ 1, consider a sequence of monic polynomials Φ n (z), which are orthogonal with respect to the semi-classical weight of the form w(z) = |z| α w 0 (z), (36) (where the modified weight in (4) is a special case when w 0 (z) . Then the orthogonal polynomials Φ n (z) associated with the weight in (36) satisfy the differential-difference relation From orthogonality and integrating by parts, we have where v 0 (y) := cy 6 + t(y 4 − y 2 ) .

Lemma 3.
The ladder coefficients A n (z) and B n (z) defined by Lemma 2 obey the relation Proof. From the definition of A n (z), we have that which completes the proof. (42) is the well-known supplementary condition (S1), which given in [25,35]. The other supplementary condition [25] is given by

Remark 9. Equation
where v(z) = − ln w(z), since w(x) is non-negative in an interval [a, b] ⊆ R.

Lowering Operator for the Modified Freud-Type Weight in
and hence Theorem 10. For monic modified Freudian polynomials S n (x; t), we have the differential-difference relation where the expressions Q n and Ω n are given in Equations (23) and (24) respectively.

Proof. From Equations (37b) and (46), we have that
and the integrand in Equation (48) is odd and hence it vanishes.
Similarly, from Equation (37c), using orthogonality and three-term recurrence relation Equation (6), we obtain where the expressions for Q n and Ω n are given respectively in (23) and (24).

More Non-Linear Difference Equations for the Recursion Coefficients
This subsection provides several nonlinear difference equations fullfilled by the recurrence coefficients corresponding to the weight (4).

Remark 12.
It is important to acknowledge the work of Clarkson and Jordaan [15] on differential-difference and differential equations as we have similar results (see also [12] for the quartic case). In this contribution, our choice of the modified Freud-type weight (4) is a more general measure modification by dµ(x; t) = e t(x 4 −x 2 ) dµ(x; 0), (a similar one is given in [21]); while the one in [15] is measure deformation by dµ(x; t) = e tx 2 dµ(x; 0). This choice led us to a different version of results of new type of Toda-like relations, new non-linear recursion relations, differential equations and some properties of the zeros for the polynomials. For the semi-classical weight under consideration, we have obtained differential-recurrence relations using two different methods; the method ladder operators and Shohat's quasi-orthogonality method [37]. In the following subsection, we show that both methods provide us same results.

The Differential Equation Satisfied by Modified Sextic Freud-Type Polynomials
Theorem 13. For the modified Freudian weight (4), the monic orthogonal polynomials S n (x; t) satisfy the linear second-order ODE (with rational coefficients): and T n (x; t) = 18cβ n x 2 + 6cQ n − (2λ + 1)Ω n x 2 + 4tβ n where Q n and Ω n are given as in (23) and (24) respectively.
One can see that a more expanded expression of the coefficients of the differential equation in Equation (76) can be obtained using the symbolic package in Maple but will be very cumbersome. (4), the lowering operator (cf. Theorem 10) is rewritten as

Symmetrizing Semi-Classical Modified Airy-Type Weight
Orthogonal polynomials associated with a symmetric measure can be obtained from classical orthogonal polynomials by means of quadratic transformation (cf. [12,30]) and for general case, we refer to [39]. For instance, a class of generalized Hermite polynomials is generated from Laguerre polynomials while a class of generalized Ultraspherical polynomials is obtained from Jacobi polynomials [30].
We will, next, show that symmetrizing the semi-classical generalized Airy-type weight function (cf. [34]) gives rise to the modified sextic Freud-type weight function (4).
n (x; t)} ∞ n=0 be a sequence of monic semi-classical generalized Airy-type polynomials, orthogonal with reference to the semi-classical weight in (84).
Let's define the monic polynomials where the polynomial is of degree n and monic in nature. Since xP (λ+1) n ∈ P n+1 , we now write xP (λ+1) n in terms of the semi-classical generalized Airy-type basis {P where the expression a n+1,k (t), with fixed t ∈ R, is given by Using Equation (85), we have that We see that a n+1,n+1 (t) = 1 since P (λ+1) n (x; t) is monic, and Equation (86) becomes Computing Equation (87) at x = 0, we obtain a n+1,n (t) = − P and hence Hence, the polynomial sequence {S 2m (x; t)} ∞ m=0 is orthogonal with respect to the symmetric It is proved in [30] that since the integrand is odd. Thus, the polynomial sequence {S n (x; t)} ∞ n=0 is orthogonal with reference to the modified weight (4). We note that W λ (x, t) = |x| −1 W λ (x 2 ; t) = |x| 2λ−1 exp −[cx 6 + t(x 4 − x 2 )] is another symmetric dual weight for the semi-classical generalized Airy-type weight function (cf. [32]).

Remark 15.
By employing quadratic transformation to the non-symmetric measure (84), we can see that the symmetrization process due to Chihara [30] preserves orthogonality and hence the modified sextic Freud-type polynomials follows from generalized Airy-type polynomials [34].

Electrostatic Interpretation of the Zeros
In this section, let's first consider a sequence of orthonormal polynomials {Ψ n (x)} n≥0 with respect to a weight function w(x) = exp(−v(x)) satisfying a three-term recurrence relation [40] xΨ n (x) = a n+1 Ψ n+1 (x) + b n Ψ n (x) + a n Ψ n−1 (x), n ≥ 0, a n > 0, with Ψ −1 (x) ≡ 0 and Ψ 0 (x) ≡ 1. It is known in [25,41] that under certain assumptions on w the polynomials Ψ n satisfy differential-difference relation where the coefficients A n and B n are explicitly given in terms of w, the recurrence coefficients a n , and the values of Ψ n at the endpoints c and d. Consequently, Ψ n obeys a linear differential equation [42] ] + a n a n−1 A n (x)A n−1 (x) .
(see Lemma 5 for its equivalent monic form). Specifically, if X * ⊂ (α, β) n is the ordered set of zeros of Ψ n , then By considering an electrostatic model for X * with logarithmic interaction between particles and an external field (k n is any appropriate standardization constant, taken in [43] equal to a −1 n ), one can see that this external field has two components: the first term in the right hand side of Equation (90) has its origin in the orthogonality weight w(x) = exp(−v(x)), and Ismail [25] named it the long range potential. The second term, known to be the short range potential, is slightly troublesome, but informs many features of the classical models and at the same time, allows to provide a generalization of the electrostatic interpretation (see [42][43][44] for more details). M. E. H. Ismail in [43] states that, assuming w(x) > 0 on (a, b), both v and ln(A n ) in C 2 (a, b), and the external field (90) convex, the total energy has a unique point of global minimum, which is precisely X * (that is, the zeros of the orthonormal polynomial Ψ n ).
Clearly, what is used here is in fact the relation (89) equivalent to the fact that X * is a critical point for E(X), and convexity assures that this is the unique point of minimum (see also [25,43,45]).

Properties of Electrostatic Properties of the Zeros for the Modified Freud-Type Weight
For electrostatic interpretation of the zeros for classical orthogonal polynomials, we refer to [1]. In this discussion, an electrostatic representation of distribution of the zeros for modified Freudian polynomials is provided.
In [40], the authors considered a perturbation of the quartic Freud weight (w(x) = exp(−x 4 )) by the addition of a fixed charged point of mass λ at the origin. (See also [43] and the recent work in [28] for the Freudian-Sobolev case). It was shown in [40] that the semi-classical quartic Freud polynomials obey a differential equation of the form Equation (88), and the electrostatic model was studied. Recall from Equation (45) that the potential v for the modified Freudian weight (4) is given by The following result immediately follows from Remark 14. Proposition 16. For the semi-classical weight (4), we have the following equation that corresponds to Equation (80): where the ladder coefficientsĀ(x, n) andB(x, n) are given bȳ B(x, n) = 6cβ n x 5 + (6cQ n + 4tβ n ) where the coefficients A n (x) and B n (x) are given in Equation (47) and the expression for Q n and Ω n are also given in Equations (23) and (24) respectively.
The electrostatic interpretation of the zeros corresponding to the semi-classical Freudian weight (4) is attained by finding the zeros of the coefficientÂ(x, n), which provide us the location of some fixed charges. Corollary 1. For n ≥ 1, the coefficientĀ(x, n) in Equation (92) has the following six roots {r j } 6 j=1 are given by where the values of A, B and E are respectively where the parameter t is assumed to be a positive constant, c > 0 and Q n are given in (23) respectively.

Remark 17.
(i) Generally, one can see that the coefficientĀ(x, n) in Corollary 1 has one double root 0, and the other four roots in the complex plane. We also see from Equation (94) that A > 0, B > 0 and Q n > 0 since the recurrence coefficient β n is always positive. For the particular cases, we have the following: (ii) When B 2 − 4AE > 0, the coefficientĀ(x, n) in Corollary 1 has a root r 1,2 (n) = 0 of multiplicity 2, two real roots r 3,4 (n) and the remaining two roots r 5,6 which are complex for n) in Corollary 1 has a root r 1,2 (n) = 0 of multiplicity 2, the remaining other roots r 1,2,3,4 (n) are entirely complex.
We now propose an electrostatic model in the presence of a varying external potential from the second-order linear differential equation obtained in Section 4.4. In other words, the electrostatic interpretation of the zeros for polynomials orthogonal with respect to (4) using Theorem 13 and Remarks 14 is studied.
We begin by denoting {x n,k } n k=1 for the zeros of S n (x) associated with the semi-classical weight (4). Evaluating the second-order differential equation at x n,k , we have Then Applying the following property (cf. [43,46]) for 1 ≤ k ≤ n, and where t is assumed to be a fixed positive real constant. For the semi-classical weight in (4), the total external potential V(x) is the sum of an external field, v(x) = − ln W λ (x; t) = −(2λ + 1) ln |x| + cx 6 + t(x 4 − x 2 ), and a varying external potential 1 2 ln |A(x, n)| − ln |x|, called by Ismail [43] long range field and short range field respectively; i.e., which turns out to be We consider the potential energy at x of a point charge q located at s is −q ln |x − s|. We notice from Equation (98) that the zeros of the coefficientĀ(x, n) give us the position relying on n of six fixed charges.Then the external field is generated by a fixed charge +1 at the origin-due to a perturbation of the weight function-plus six fixed charges of magnitude −1/2; two of them located at the real positions r 3 (n), r 4 (n), and the remaining ones at the complex positions r 1,2 (n), r 5,6 (n) (see Remark 17).
Let us introduce the following electrostatic model corresponding to the weight (4): Consider the system of n movable unit positive charges in n distinct points {x n,i } n i=1 of the real line in the presence of the total external potential V(x).
By finding the zeros ofĀ(x, n), we obtain the position relying on n of six fixed charges distributed on the real line. By denoting the position vector as x = (x n,1 , x n,2 , . . . , x n,n ) where x n,j < x n,k if j < k and considering the total energy of the system as we attain Equation (97), which is the derivative of the energy function. This means the zeros of the modified Freud-type orthogonal polynomials are critical points of the energy function. From the physical point of view, these kind of possibilities may want to correspond to the steady or unstable equilibrium situations. The equilibrium is understood as the zero gradient of the complete energy of the system. Stable equilibrium means the existence of a global minimum of the total energy. However, studying the stability for the equilibrium configuration (global minimum) seeks a thorough discussion, which we will pass over here. Nevertheless, we partly resolve this quest by targeting in the study of the local minima of the energy function. For this purpose, the Hessian matrix is considered as we will infer when E(x) attains local minima at the zeros of the monic polynomials S n (x) with reference to the weight (4). Indeed, taking into account and the matrix H is a symmetric real matrix. It is noted in [47] that the Hessian matrix H is positive definite when the matrix H is strictly diagonally dominant and all its diagonal terms are positive. In this case, by setting conditions, which guarantee the equilibrium position of the proposed system will be reached at the zeros of the monic Freud-type polynomial [40] and coming along this way, one can find the electrostatic equilibrium position in the presence of the external field at the zeros {x n,k } n k=1 for the semi-classical modified Freud-type polynomials.

Discussion
In this work, we have investigated certain properties of monic orthogonal polynomials associated with a scaled Freud-type weight upon deformation of the Freudian measure |x| α exp −Cx 6 , C > 0 by exp t(x 4 − x 2 ) . This slightly different Freudian measure could be considered as a generalization in some sense for the weights in (2). For this case, we have investigated certain properties of orthogonal polynomials corresponding to the semiclassical weight given in (4). Of these properties, we have studied the finite moments, nonlinear recursion relations and differential-recurrence relations for the recurrence coefficients as well as a linear differential equation for the polynomials themselves. This work has certain similarity with the work in [12]. However, the motive of the current work is different as it deals with the modified Freud-type weight function, which already involves a higher-order of the polynomial in the exponential factor in the semi-classical weight (4). A similar work for semi-classical Laguerre weight has been given in [11] and for the Freud-type weight in [18,36].
Besides, this work mainly shows the relation between both the semi-classical polynomials that are orthogonal with respect to the modified weight in (4) as well as the expression of the recurrence coefficients corresponding to these polynomials. Since the recurrence coefficients for sextic Freudian measure (4) are not explicitly found, it is important to obtain several nonlinear recursion relations as well as differential-recurrence relations associated with the weight function (4). Following the method given by J. A. Shohat [12,37] using quasi-orthogonality, we applied this technique as a powerful tool to obtain the ladder relations associated with the weight under investigation. The differential-difference equations obtained are good references to characterize the corresponding semi-classical modified Freudian polynomials. The second-order linear differential equation also follows by combining the ladder relations with the three-term recurrence relation in a similar way to the work in [12]. As an application of the results, an electrostatic interpretation of the zeros from the obtained linear differential equation with rational coefficients is also explored.

Conclusions
In this paper, we have investigated certain properties of modified sextic Freudian polynomials. Knowing the fact that the recurrence coefficients associated with modified Freudian measure are not straightforwardly formulated; we have derived certain properties such as several non-linear recursion equations, Toda-like equation, differential and difference equations satisfied by the recurrence coefficients as well as the polynomials that are orthogonal with respect to the considered weight function.
Special attention, using the method of ladder operators, is given for the recurrence coefficients associated with the modified Freudian weight given in (4) as this weight function is derived from the generalized Airy-type weight (84) using Chihara's symmetrization process [30]. A differential-difference equation as well as linear second-order differential equation associated with the semi-classical weight (4) were obtained, both via methods of ladder operators as well as Shohat's method of quasi-orthogonality [37]. As an application of the properties, we have investigated an electrostatic interpretation of the zeros from the obtained linear second differential equation with rational coefficients. Following the work and methods used in [15], some asymptotic properties of the recurrence coefficients in terms of Hankel determinants, asymptotics of the differential equation satisfied by the studied orthogonal polynomials will be addressed in a future contribution. We also hope our results motivate further investigation of the recurrence coefficients via Hankel moments of semi-classical orthogonal polynomials in relation to certain (discrete) integrable systems and identifying these connections would also be an interesting continuation of this work.