1. Introduction
Given a complex and monic polynomial
(see Comment A1), it is always possible to define a matrix with a specified arrangement of the polynomial coefficients as its entries, such that
coincides with the characteristic polynomial of the matrix. The set
of all these
matrices with complex entries, sharing the same characteristic polynomial
, is infinite and can include both derogatory and non-derogatory matrices (see Comment A2). In fact, we observe that, by definition, the
Frobenius Matrix [
1] of the shared characteristic polynomial always belongs to
. A remarkable property of this matrix, which stems directly from its construction, is the coincidence between its characteristic and minimal polynomials, whatever
. Therefore this matrix is classified as non-derogatory and, following Horn and Johnson [
2], it is known as the
Companion Matrix (CM) of its characteristic or minimal polynomial (see Comment A3).
Henceforth, we refer to the Frobenius Matrix as the Frobenius Companion matrix (
) of
. When the algebraic multiplicity of each of the
n eigenvalues of the
is 1, any matrix in
is non-derogatory, since these
n distinct eigenvalues are all necessarily roots of its minimal polynomial, which therefore has degree
n [
2]. We also remark that this condition guarantees that each matrix in the
set is diagonalizable [
2] and that, consequently, all matrices
can be generated from the
by means of a similarity transformation. In this way, by definition and under the conditions established for the spectrum
of
,
includes all and only the companion matrices of
.
When, instead, the distinct roots of the common characteristic polynomial are
,
includes infinite derogatory matrices, which cannot be structurally similar to the Frobenius matrix (see Comment A4) and to any non-derogatory matrix belonging to
. For example, consider that the set
contains the diagonal matrices whose
n entries are nothing but the
p distinct roots of the characteristic polynomial repeated as many times as their multiplicities, whose sum is
n. The degree of their characteristic polynomial is
n, while the degree of their minimal polynomial is
[
2]. Therefore, these matrices are derogatory, as are the infinite matrices generated from them by similarity.
The Frobenius matrix dates back to 1879 and was given in the form [
1]:
by the German mathematician Ferdinand Georg Frobenius [
3]. Sometimes, in the literature it is presented in three other unitarily transformed forms, all parameterized in terms of the
n coefficients of the polynomial
and with the same number of entries equal to 0 or 1 as in (
1) [
4]. When the
n eigenvalues of the Frobenius matrix are distinct, the unitary matrix that yields the diagonal form of (
1) is the Vandermont matrix of its
n eigenvalues. This property, as well as other properties and some applications of the Frobenius matrix may be found in [
4].
Despite the fact that the first CM was proposed more than 140 years ago, the generation of different CMs (with further properties) has attracted a great deal of applied research activity. CMs emerge naturally in mathematical methods for finding and characterizing the roots of polynomials [
5,
6,
7,
8,
9,
10] and can be applied as well in the determination of solutions of high-order scalar linear differential and difference equations [
11]. CMs can give matrix representations of some fields [
12] and are widely used in control theory, for example in writing the controllable canonical form associated with the transfer function of a system [
13]. The product of CMs is also used in the study of random walks and Markov chains [
14].
The original structure of the FCM shows a modest level of flexibility and, in fact, has stimulated the search and the emergence of generalizations leading to new proposals of CMs that pave the way to applications beyond the FCM. To this end, a successful strategy is based on a radical change of the basis in which the characteristic polynomial
of the FCM is represented. By definition,
is monic and is written as the sum of
and a linear combination with coefficients in
of
. All these
powers of
z constitute the monomial basis, which can be replaced by other polynomial bases. In this way, one can introduce new (still non-derogatory) CMs with nonvanishing elements on the main, sub, and super diagonals. For example, the Chebyshev basis, independently adopted by Specht [
15,
16] and Good [
17], led to a new CM called the
colleague matrix. This approach has been further generalized by Barnett [
4], who considered a basis of orthogonal polynomials and named the newly emerged CMs
comrade matrices. Subsequently, Barnett proposed to call
confederate matrices the CMs arising from the use of a general polynomial basis [
18].
The applications dedicated to the classic problem of finding the real zeros of a real coefficient polynomial of arbitrary degree deserve a special mention, because in this context the CMs have inspired a different approach, alongside exquisitely mathematical and computational investigations [
19,
20,
21,
22,
23,
24]. In the last decade, new quantum theory-based root-finding algorithms exploiting the construction of Hermitian companion matrices [
25,
26,
27,
28,
29,
30] have also been proposed, thus increasing the interest in the study of CMs and related matrices in the rapidly growing research field of quantum computing. Furthermore, recent studies [
31] have shown the opportunity of using CMs in mathematical constructions useful to the investigation of quantum entanglement, quantum state tomography, and quantum information in general.
The above is the general context for the main question addressed by this paper: is it possible to find a CM of a real or complex polynomial that is also hermitian/unitary, for example, or possesses some other prescribed special matrix structure? This question has so far been answered only partially.
In [
32] the CM for a
with real coefficients and real zeros is constructed as a real symmetric tridiagonal Hermitian matrix. This provides a complete solution to a problem raised and partly solved by M. Fiedler in [
33], which has rekindled the interest in the general structure of CMs. Note that the Frobenius matrix is itself a Fiedler matrix after a reverse permutation matrix similarity. In [
34] all CMs are characterized in terms of combinatorial structure to generate new CMs. It is interesting to note that both the Frobenius and Fiedler CMs are sparse matrices, as they have
nonzero elements [
35]. A new class of sparse CMs, also known as intercyclic CMs, was introduced in [
34] and includes the Fiedler matrices as a special case. In [
35] the non-sparse CMs are introduced noting that they are not connected with the sparse ones by a reverse permutation matrix similarity.
To the best of our knowledge, the question of whether the CM of a complex polynomial can be sought as unitary or Hermitian is still open.
1.1. Purpose and Contribution of this Study
CMs generally show relatively limited versatility due to their combinatorial structure. For example, the FCM is never Hermitian or unitary. Searching for CMs that satisfy additional constraints of this kind is important when a given class of parametric polynomials is designed to reach a reliable theoretical control of a quantum physical scenario. Finding exact flexible solutions to well-defined inverse problems of this kind is a target of the present study. We stress that, for application purposes, we often do not need to associate non-derogatory matrices to a given polynomial. To emphasize this particular aspect of our matrix construction, we introduce the term almost-companion matrix (ACM) to refer to matrices that have a given polynomial as their characteristic polynomial but can be derogatory. Clearly, every CM is also an ACM, i.e., the set of all ACMs is a superset of the set of all CMs of a given polynomial.
In short, in this study, we address the following inverse problem: given a real or complex polynomial of any kind (for example, it may belong to a special class of polynomials, which is reflected in some special condition satisfied by its coefficients), we find a parametric ACM (see Comment A5) that satisfies preassigned conditions (e.g., those required to be Hermitian or unitary) and whose characteristic polynomial coincides with the given one. For definiteness, our investigation is here limited to complex polynomials of the third degree (). The extension of the analysis to higher-degree polynomials is discussed.
The solution of the inverse problem outlined above for polynomials of order is accompanied by some useful applications. The relaxed constraints that characterize an ACM, as compared with a CM, allow the freedom to search, from the outset, for a trial ACM of a given that is a Hermitian, unitary, or positive matrix, for example. If our inverse problem can be solved systematically through an approach that finds such ACMs whatever the given , then we readily have at our disposal a good platform for successful applications to problems such as those mentioned below.
1.2. Physical Applications
Constructing an ACM of a given generally complex polynomial, in addition to being interesting in itself, also has considerable applications. In elementary algebra, for example, it could be a solution tool for counting the number of the real roots (and consequently that of the complex roots in the case of a complex polynomial) of a real or complex polynomial. In addition, it may help determine a (or the only) real root of a real polynomial of odd degree. In quantum mechanics or quantum information, it could provide new parametric representations of the density operator or the evolution operator of a physical system living in a finite-dimensional Hilbert space.
The results found for a generic complex polynomial can be applied to the important particular case of a real polynomial. Investigating such a link is certainly of interest in physics. For example, in classical physics, cubic real polynomials appear when looking for the principal axes of symmetric Cartesian tensors of rank two, such as inertia or magnetic/electric dipolar tensors [
36]. In quantum mechanics, they enter the scene as characteristic polynomials of any observable of a physical system that lives in a three-dimensional Hilbert space, such as a three-level atom or a qutrit. Recipes for constructing an ACM of a cubic polynomial possessing real roots after the appropriate assignment of parametric real coefficients could provide an easy way to build, e.g., Hamiltonian qutrit models on demand for control purposes, or even the density operator describing a mixed state of a three-level atom.
The paper is organized as follows. In
Section 2, we formulate the inverse problem consisting of the search of the ACM for a generic complex polynomial. In
Section 3, we construct the ACMs for a generic cubic complex polynomial. Through these ACMs, we introduce a way to find the roots of the given polynomial without using the Cardano-Dal Ferro formulas. The case of the polynomial with real coefficients is also discussed in detail. In
Section 4 we present an application in quantum mechanics, constructing on demand the density matrix of a qutrit system as an ACM.
Section 5 shows the construction on demand of the unitary ACM of a qutrit. Possible extensions to higher-degree polynomials, as well as further possible applications, are discussed in
Section 6.
2. Formulation of the Inverse Problem
Consider a matrix
, where
is the set of all
matrices over the complex field
. Denoting
the identity matrix, the monic polynomial in the complex variable
z
is, by definition, the characteristic polynomial of
A and belongs to the set
of all complex monic polynomials of degree
n.
It is well known that its
n coefficients
,
contain information about the elements of
A that is invariant under arbitrary similarity transformations (see Comment A6). In fact,
is the sum of all the principal minors of order
k of
A. In particular,
and
. The profound interrelationship between a matrix and its characteristic polynomial becomes even more surprising considering the Cayley–Hamilton theorem (see Comment A7) [
2] and/or Newton’s identities [
37], which reveal the existence of finite algebraic expressions for the coefficients of the characteristic polynomial of a matrix in terms of traces of powers (up to
n) of the matrix. [
38,
39].
The function is surjective but not injective and, hence, it cannot be inverted. In fact, it is easy to convince oneself that, for any element , is indeed an infinite subset of , since by definition it consists of all and only the matrices belonging to whose characteristic polynomial is .
Thus, while the direct or forward problem of finding the characteristic polynomial of a given
matrix is certainly well-posed according to Hadamard [
40], conversely, the problem of finding a matrix
generating a given complex polynomial
is an ill-posed inverse problem [
41,
42], as it manifestly violates Hadamard’s uniqueness requirement, considering that every element
is a solution to the problem.
It is possible to overcome such an ill-posedness by introducing a restriction
of the function
C to a subset
of
, which is injectively and surjectively valued on
. To this end, let us first observe that the function
C is surjective and, by definition,
is the characteristic polynomial of all and only the matrices belonging to
. Moreover,
when
. Therefore, the infinite subsets
of
corresponding to the infinite
n-degree polynomials
represent a partition of
. We can say equivalently that we are introducing in
the equivalence relation
consisting in the condition that
and
share the characteristic polynomial and thus belong to a given equivalence class
. At this point, we define the subset
by choosing one element from each equivalence class. According to Zermelo’s postulate,
can always be constructed (in infinitely many ways in the present case, since each equivalence class is infinite), and the cardinality of its intersection with
is precisely one for any
by construction. Therefore, every
one-to-one function
obtained by applying the axiom of choice to the quotient set
to generate
is invertible. Hence, function
defines a
-dependent Hadamard well-posed inverse problem, whose solution, by construction, can be given in terms of
in the form
We remark that different legitimate choices of the subset
lead to different inverse problems, all well-posed in the fixed domain
, and the corresponding solutions (
3) differ in the generally non-similar images of one or more polynomials
.
We also point out that any derogatory matrix
D cannot be classified as a companion matrix of its characteristic polynomial
, since
D annihilates a polynomial having a degree lower than that of
[
2].
In this paper, a matrix whose characteristic polynomial coincides with a given polynomial
is called an
almost-companion matrix of
. Clearly, any CM of
is an ACM too. A derogatory matrix
D such that
is an ACM. In addition, a matrix similar to an ACM is still an ACM. The converse of this statement is generally false: two ACMs of the same given polynomial are not necessarily similar [
43]. The set of all the ACMs of
cannot be generated by similarity transformations starting from an assigned ACM, since this set always includes both derogatory and non-derogatory matrices.
3. Almost-Companion Matrices of a Cubic Complex Polynomial
In this section, we focus on the search for an ACM of the third-degree polynomial
which is the canonical form of
obtained by the translation
The generally complex numbers
p and
q in (
4) are related to the coefficients of
as follows:
We denote
the ACM of
defined by
which means that
is the simple recipe to obtain the corresponding ACM
of
from
. This analysis sheds light on the advantage of first deriving
for the simpler canonical form of a given polynomial and then finding
from the straightforward relation (
9).
Next, we formulate a trial ACM
of (
4). To this end, we observe that, in accordance with the Vieta-Girard formula for the sum of the roots of (
4) [
44], the absence of the quadratic term in
implies that
. Moreover, every matrix with elements in
is unitarily equivalent to a matrix with equal main diagonal elements [
2]. Thus, it is legitimate to set the diagonal elements of our trial
equal to zero. In constructing an ACM of (
4), we aim to write its non-diagonal elements in such a way that, in the particular case of a real cubic
and, hence,
, the trial matrix
becomes structurally Hermitian provided that
p and
q in (
4) satisfy specific conditions, which will also be derived from our approach. The feasibility of this approach will highlight the greater flexibility of the ACMs compared to that of the CMs.
Following this strategy, we propose the following trial
:
where the minus sign was introduced for convenience, considering the form of the characteristic polynomial. In Equation (
10),
is real and positive,
is real, whereas
is, in general, a complex number. It is readily seen that, when
or
and
is real, the matrix
is Hermitian, consistent with our search strategy. It is useful to note that the complex conjugate of
is
, where
denotes the conjugate of
.
The characteristic polynomial of
is
Then, identifying polynomial (
4) with (
11) yields:
Given
p, the first Equation (
12) allows us to fix
and select
(in an infinite set) as follows:
Defining, for
, the complex parameter
the second Equation (
12) becomes an elementary trigonometric equation in
:
which admits infinitely many solutions for any
; in fact, similarly to the cosine function of a real variable, the complex cosine function is even and periodic with period
.
Using Euler’s formula, Equation (
15) is easily transformed into a quadratic equation in the variable
, whose solution leads to
Due to the presence of the multi-valued complex function
, expression (
16) represents the set of infinite images of
generated by the inverse of the non-injective cosine function over
. Therefore, strictly speaking, the expression found for
cannot be introduced as it is in the matrix
. In fact, the three parameters appearing in the trial ACM (
11) of (
4) must be single-valued functions of the complex coefficients
p and
q. For our purposes, therefore, we now need to extract a specific single-valued complex function from the multi-valued function
.
The single-valued complex function that we use here is the principal value
of
, which is obtained from (
16) by substituting the multi-valued functions arg and ln with their principal values, denoted Arg and Ln, respectively. This choice amounts, by definition, to constructing the principal value
of function
, which is mostly used in the literature [
45,
46]. It is worth noting that equation (
16) can also be written in terms of
, but then the use of the principal value in the resulting expression for
would have some drawbacks, as is discussed in detail in [
47].
The procedure described above gives
where
ln denotes the ordinary real logarithm of its positive argument and the first term in (
17) is the principal value of
, which, by definition, generates real images in
. Equation (
17) provides the algebraic representation of the complex single-valued function
whatever
. Note that, for any real
such that
, the imaginary (real) component of
identically vanishes.
We determined all of the ingredients for constructing the trial ACM of (
4) when
. In accordance with (
17), this ACM is a generally non-Hermitian matrix that can be written as follows:
where
Thus, based on (
9), the ACM of (
5) has the form
We can find an ACM of (
4) for
through the same kind of approach, beginning with a matrix different from (
10). Our solution, denoted by
, can be cast as follows:
For completeness, we also write
:
since
from (
17). It is not difficult to see that the eigenvalues of (
22) are 0 and
.
Above, we formulated and solved the inverse problem of finding ACMs of generic cubic complex polynomials. Next, by exploiting these ACMs, we will present a way to find the roots of the given polynomial without resorting to the Cardano-Dal Ferro formulas, and then we will delve into the implications of having a polynomial with real coefficients on the form of the ACM and its roots.
3.1. Roots Characterization
Next, we assume that
. Then, the eigenvalues of
are the roots of (
4) and, hence, the roots of the characteristic polynomial of the matrix appearing in the right-hand side of (
18), each multiplied by the pre-factor
. This characteristic polynomial
in the unknown
has the form
The cosine representation of the free term in polynomial (
23) is remarkable because it allows one to guess, at first glance, one of its three roots and then to exactly construct the other two by simply reducing the cubic polynomial
to a quadratic polynomial. In fact, without resorting to the well-known Cardano-Dal Ferro formulas [
48], and using instead the elementary triplication formula for the cosine function
, which also holds in the complex field, it is immediate to see that the generally complex expression (see Comment A8)
is a root of (
23), whatever the complex coefficients
and
q. The algebraic representation (
17) of
is the key to explicitly write the
p- and
q-dependencies of the real and imaginary components of the algebraic expression for
, which are obtained using Euler’s formula as
where
and
are defined in (
17). The other two roots are easily found to be [
49]
Equations (
24) and (
27) express the roots of
as functions of parameter
, which appears in the first and last anti-diagonal terms of matrices (
18) and (
20). We, thus, conclude that our procedure to construct the ACM also yields the roots of the (generally complex) characteristic polynomial.
It is interesting to highlight the conditions for a complex polynomial to admit real roots (note that the roots can never be all real, however, if the imaginary part of at least one of the polynomial coefficients is nonzero). To this end, it is convenient to start from the polynomial form (
5). We write the three coefficients of (
5) as
, with
. If a real root
r of (
5) exists, it must satisfy the equation
, which results from equating to zero the imaginary part of the polynomial.
is clearly a necessary condition for the existence of the root
r, and the two only possible expressions of
r are
. Then, any of these two expressions is indeed a root of (
5) only if it satisfies the additional condition
, which results from the real part of the polynomial. We can thus state that the cubic polynomial (
5) has at least one real root if and only if the inequality
and the last condition are both met. In particular, when
,
r is a double root. Finally, we examine the case
(and
, since otherwise no real root exists for a complex polynomial). Applying the same procedure, one finds that a real root
of (
5) exists if and only if the condition
holds. It is easy to convince oneself that this real root has multiplicity two, being a real root of the first derivative of (
5).
3.2. Real Polynomial Case
We now investigate the special form of
in the case in which the three coefficients of (
5) are real, with the aim of establishing properties of the polynomial roots based on its almost-companion representation built above.
In this case, (
7) implies that
p and
q are real, and thus (
4) is also a real polynomial over
. Moreover, since
can only be 0 (
) or
(
), the parameter
defined in (
14) is purely imaginary or real, respectively, and, therefore, its square
is real for any
p. As a consequence, using (
17) it is not difficult to prove that, if and only if (see also Comment A9)
(which implies
, i.e.,
, and
), the imaginary part of
given in Equation (
17) vanishes, while its real part assumes the simple expression
where we used the identity
valid for any real
x such that
. It is worth noting that in the present case (
27) takes the simpler form
since (
30) shows that
.
Under condition (
29) and considering that
, the specialization of (
18) to the case under scrutiny produces the Hermitian ACM of (
4) as
where the real angle
is given by Equation (
30). The Hermitian nature of the ACM built assures that (
4), as well as (
5), has three real roots. These roots are distinct when
, while two of them are coincident if
[
50].
We can write the real roots
of (
5) as follows:
where
are the three (real) roots of (
23). Equations (
24) and (
27), together with (
30), yield
It is possible to check that this formula gives the well-known trigonometric and translated forms of the three roots of (
5) when they are real (see [
51]).
When
, only one of the three roots of Equation (
4) or (
5) (with real polynomial coefficients) is real, while the other two roots are complex conjugates. In particular, the real root corresponds to
for
and to
for
. In more detail, the three roots are
where
and
with
for
, and
where
with
for
(see derivation in
Appendix A.1). A more elaborate derivation of the roots leading to their formulation in terms of trigonometric functions can be found in [
50].
4. Almost-Companion Density Matrices of a Qutrit on Demand
A quantum system living in the Hilbert space
spanned by three orthonormal states
,
, and
is called a qutrit [
52]. A pure state of the qutrit can always be represented as a normalized linear combination of these three states. To describe an arbitrary pure or mixed state of the qutrit with the same formalism, one uses instead a linear operator
called the density operator [
53]. It acts on
and, by definition, is positive semi-definite with trace 1:
. It is well known that any positive semi-definite operator is Hermitian since its skew-Hermitian part vanishes [
2,
54]. As a consequence, any positive semi-definite operator is diagonalizable, and it is possible to show that its eigenvalues are real non-negative numbers. In particular, any density operator
is Hermitian. The three eigenvalues of the operator
describing the state of a qutrit are the populations of the three eigenstates of
. The
basis-dependent matrix representation of
, called density matrix and denoted by
, is also positive semi-definite and, hence, Hermitian. We observe incidentally that, conversely, any Hermitian matrix with non-negative eigenvalues is positive semi-definite and, if its trace is 1, it is a density matrix.
The purpose of this section is to demonstrate that our recipe for constructing the Hermitian ACM of a generic third-degree polynomial admitting three real roots provides an effective tool for writing density matrices of a qutrit on demand. It is worth noting that our approach itself does not require the support of a vector space, while the relationship to a basis of physical states appears in this application to a qutrit.
The aforementioned definition of density matrix results in unambiguous properties of the real coefficients of its characteristic polynomial. First of all, writing this polynomial in the canonical form
the condition
stemming from Equation (
29) ensures that
is real, so that the ACM (
31) of (
44) is Hermitian and, hence, has real eigenvalues. Then, turning to form (
5) of the monic polynomial through the translation (
6), the Vieta–Girard formula for the sum of the roots [
55] implies that the coefficient of the quadratic term is
. Moreover, Descartes’s sign rule [
56] requires that the four coefficients of polynomial (
5) have alternate signs in order to have three positive roots. Therefore, the characteristic polynomial of an arbitrary density matrix of a qutrit is necessarily a third-degree real and monic polynomial of the form
where
a and
b are real numbers that satisfy condition (
45) after translation (
6). One or two roots are zero if
or
, respectively, while the inequality (
45) is never satisfied for
and
.
In conclusion, under conditions (
45) and (
46),
, with
given by Equation (
31), is a density matrix. Incidentally, we point out that our inverse problem admits infinitely many non-Hermitian solutions, i.e., non-Hermitian ACMs of (
44) or (
46), such as, for example, the corresponding Frobenius companion matrix. Therefore, the explicit construction of the Hermitian ACM of (
46) and, more generally, of any real third-degree polynomial with only real roots is a successful outcome of our search strategy (
10). This recipe, in turn, forms the basis of the application presented below.
Let us introduce the set
of all density matrices of a qutrit in a given basis
of
.
be the binary relation in
defined as follows:
and
are in the relation
if they are ACMs of the same polynomial
defined in (
46). This relation, expressed by writing
, is an equivalence relation as it is manifestly reflexive, symmetric, and transitive.
is thus partitioned by
. The quotient set
consists of all equivalence classes of
with respect to
. Each equivalence class, which comprises all density matrices with the same characteristic polynomial, is represented by one (arbitrarily chosen) of its elements,
, and is commonly denoted by
. This is where our result (
31) enters the scene, providing an easy way to parameterize the quotient set of
.
It is always possible to use the matrix
as the representative element of the equivalence class consisting of all elements of
sharing the characteristic polynomial (
46), which, in turn, is uniquely associated with its canonical form (
44). In this way, we establish a one-to-one correspondence between
and the set
of polynomials (
44). This correspondence amounts to parameterizing the quotient set of
in terms of
p and
q. The most ambitious target of parameterizing set
is discussed in a recent topical issue [
57]. It is worth emphasizing that a density matrix
belongs to the class of equivalence
if and only if it can be unitarily generated from
, since its characteristic polynomial, trace, and Hermiticity are unitarily invariant, thus implying the invariance of the positive semi-definiteness. Therefore, while two similar matrices are ACMs of the same polynomial, a similarity transformation of a density matrix does not generate, in general, a density matrix [
2].
Our parameterization of
in terms of the coefficients of its characteristic polynomial written in canonical form provides the theoretical basis for constructing, on demand and in a prefixed basis of
, almost-companion density matrices of any assigned polynomial
fulfilling the condition
. We illustrate the concrete applicability of our recipe by constructing an almost-companion density matrix starting from the polynomial
The translation:
yields
so that, in this case,
while
q vanishes. Exploiting Equation (
30), we easily have:
We have, thus, obtained the few ingredients necessary to build an almost-companion density matrix of the given polynomial (
47) as the sum
of the representative element of the corresponding equivalence class and the matrix
, in accordance with the realization (
20) of (
9). That is, the density matrix has the form
Note that, while in this case it is easy to find the roots of (
48) directly, and then those of (
47), in general the roots of the polynomial
corresponding to a given
can be found using (
33).
We stress that any matrix equivalent to (
50) through a unitary transformation
is an almost-companion density matrix of (
47) and vice versa. For a given basis, each density matrix thus obtained describes a different (generally mixed) state of the qutrit. If, instead, the unitary transformation is interpreted as the generator of a change in the basis
, the matrix obtained represents the same density operator in the new basis
.
6. Discussion and Conclusions
Finding a matrix with an assigned characteristic polynomial is a classic inverse problem solved by Frobenius a long time ago. In
Section 2 we observed that there are infinitely many other solutions, generally not equivalent to the one found by Frobenius. The exhaustive description of the set
of all complex matrices sharing the same characteristic polynomial is still an open problem. Among the reasons for the missing solution to this problem, it is useful to consider that, even if
is invariant under similarity transformations, it includes non-similar (and in general not even equivalent) matrices. Another problem is that the structure of the non-empty subset of non-sparse CMs belonging to
, unlike the set of sparse CMs [
35], has not yet been fully characterized.
A related inverse problem, stimulated by applications of current interest to both physicists and mathematicians, is the search for ACMs constrained to possess prescribed structural properties, such as, for example, unitarity, positive semi-definiteness, or Hermiticity.
The first focus of this study is the construction of a new ACM of a generic monic and complex third-degree polynomial , characterized by versatility for applications. This objective is pursued by parameterizing the elements of the ACM in such a way that they lend themselves to additional constraints dictated by specific problems (of which only the structural properties are exploited). To the best of our knowledge, our investigation of inverse problems of this kind opens a new fruitful chapter in this research area whose central goal is the proposal of new ACMs which, in particular, can be CMs. We address the three above-mentioned constrained inverse problems, providing methodology and results that aim for broad applicability and are potentially transferable to solving analogous problems involving higher-degree polynomials.
The adopted step-by-step approach builds new specific classes of unconstrained or constrained matrices as ACMs of suitably given polynomials, relaxing from the outset any condition on the degree of the minimal polynomial. In particular, the strategy implemented, as well as the mathematical tools used, is not influenced by any FCM-based or FCM-inspired technique.
The elements of the ACM that we construct as a solution to the general inverse problem are single-valued complex functions of the coefficients of the given generic polynomial. Exploiting the structural properties of this ACM, we find the algebraic expressions for the three, generally complex, roots of its complex characteristic polynomial.
It is remarkable that, when the polynomial becomes real, the associated general ACM smoothly becomes Hermitian under easy-to-find necessary and sufficient conditions on the coefficients of the polynomial described by (
45). Using simply the fact that the ACM becomes Hermitian if and only if
, we are able to extend the trigonometric representation of the three roots of the real polynomial to all possible cases, that is even when (
29) does not hold. This representation is obtained without resorting to the well-known Cardano-Del Ferro formulas.
We emphasize that the FCM of a characteristic polynomial does not undergo any structural change when the coefficients of the complex polynomial become real, thus providing no additional information on the polynomial roots. For this reason, we claim that the FCM has a lower flexibility than our ACM. We show how to use this flexibility to obtain an ACM of a prescribed characteristic polynomial on demand, applying our procedure to the important problem of finding a density matrix, particularly that of a qutrit.
A second, novel constrained inverse problem addressed here consists in finding a unitary ACM of a generic polynomial with three roots of modulus one. Excluding the trivial case in which the unitary ACM can be directly given in diagonal form, we first reach the intermediate goal (interesting in itself) of parameterizing
,
, and
in the polynomial (
5) so as to set the necessary conditions on the structure of a polynomial to admit an ACM. By setting appropriate relations between the parameters through coupled geometric and analytical considerations, we further constrain the polynomial structure in such a way as to identify the set
of all and only the third-degree polynomials that admit a unitary ACM. Then, we conclude our analysis by constructing the associated ACM.
The results of this study can be further explored and usefully applied to physical and mathematical contexts, including, at their intersection, the research area of quantum computing. For example, for time-dependent parameters, the prescription of a time-dependent characteristic polynomial of or leads to a unitary time-dependent matrix, hence to the pertinent time-dependent Hamiltonian that generates the time evolution of a qutrit, whose properties can thus be traced back to the prescription of . In mathematical contexts, the analysis developed in this study can be extended to polynomials of a higher degree. For example, in the case of a fifth-degree polynomial, our protocol may provide conditions on the coefficients of the polynomial such that its roots are real. The rich analysis enabled by the use of third-degree polynomials in this study sets a clearer basis to conceive the extension of the analysis to higher-degree polynomials.