This paper characterizes the zeros of slice functions over dual quaternions, providing insights into motion polynomial factorization and advancing understanding of algebraic structures relevant to mechanism science.
Contribution
It offers a full characterization of zero sets of slice functions over dual quaternions, a specific algebra, which was previously unknown for general alternative $*$-algebras.
Findings
01
Full zero set characterization for dual quaternion slice functions
02
Enhanced understanding of motion polynomial factorization
03
Implications for mechanism science applications
Abstract
This work studies the zeros of slice functions over the algebra of dual quaternions and it comprises applications to the problem of factorizing motion polynomials. The class of slice functions over an alternative ∗-algebra A was defined by Ghiloni and Perotti in 2011, extending the class of slice regular functions introduced by Gentili and Struppa in 2006. Both classes strictly include the polynomials over A. We focus on the case when A is the algebra of dual quaternions DH. The specific properties of this algebra allow a full characterization of the zero sets, which is not available over general alternative ∗-algebras. This characterization sheds some light on the study of motion polynomials over DH, introduced by Heged\"us, Schicho, and Schr\"ocker in 2013 for their relevance in mechanism science.
Tables1
Table 1. Table 1: Scheme from Theorem 6.4 . The first column indicates the nature of the zeros of f 𝑓 f in 𝕊 x subscript 𝕊 𝑥 {\mathbb{S}}_{x} , the first row the nature of the zeros of g 𝑔 g in 𝕊 x subscript 𝕊 𝑥 {\mathbb{S}}_{x} . Each crossing lists the possible intersections between the zero set of f ⋅ g ⋅ 𝑓 𝑔 f\cdot g and 𝕊 x subscript 𝕊 𝑥 {\mathbb{S}}_{x} , without mentioning the conditions that distinguish the various possibilities.
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Full text
Zeros of slice functions and polynomials over dual quaternions
This work studies the zeros of slice functions over the algebra of dual quaternions and it comprises applications to the problem of factorizing motion polynomials. The class of slice functions over a real alternative *-algebra A was defined by Ghiloni and Perotti in 2011, extending the class of slice regular functions introduced by Gentili and Struppa in 2006. Both classes strictly include the polynomials over A. We focus on the case when A is the algebra of dual quaternions DH. The specific properties of this algebra allow a full characterization of the zero sets, which is not available over general real alternative *-algebras. This characterization sheds some light on the study of motion polynomials over DH, introduced by Hegedüs, Schicho, and Schröcker in 2013 for their relevance in mechanism science.
Acknowledgements. This work was partly supported by INdAM, through: GNSAGA; INdAM project “Hypercomplex function theory and applications”. It was also partly supported by MIUR, through the projects: Finanziamento Premiale FOE 2014 “Splines for accUrate NumeRics: adaptIve models for Simulation Environments”; PRIN 2017 “Real and complex manifolds: topology, geometry and holomorphic dynamics”. The authors warmly thank the anonymous reviewers for their precious suggestions.
1 Introduction
The concept of slice regular function is one of the possible generalizations of the notion of complex analytic function to higher dimensional real algebras. It was introduced over quaternions in [4, 5] and over other algebras in [6, 7] (see also the related work [1]). Then the work [8] extended the concept of slice regular function to all real alternative *-algebras and introduced the broader class of slice functions. The theory of slice regular functions includes polynomials and it generalizes many classical results of complex analysis. In some cases the generalization is straightforward; in many others, it is more intricate and it reveals a richer environment. This is already evident when studying the zeros of slice regular functions: as proven in [9], some form of discreteness of the zero set is always present, but the exact characterization of the zero set is algebra-specific because it depends on the nature of the zero divisors within the algebra.
Deepening the study of the zero sets of slice functions and slice regular functions by focusing on a specific algebra is exactly the purpose of this work. We focus on the algebra DH of dual quaternions, whose set of zero divisors is well understood, not only for its intrinsic interest but also for its applications to mechanism science, see [11, 12, 13, 14]. These articles introduced motion polynomials, which correspond to rational rigid body motions in the Euclidean 3-space, and studied their factorizations, which correspond to linkages producing the same motions.
The work is organized as follows. Section 2 describes the construction of the algebra DH, its properties, and its use to represent the group of proper rigid body transformations. The definitions and basic properties of slice functions and slice regular functions are recalled in Section 3, along with the definition of motion polynomial. Section 4 constructs the primal part of a slice function over DH, extending a known construction over polynomials and providing a tool to fully exploit the peculiarities of the algebra of dual quaternions. The study of the zeros of slice functions over DH is first addressed in Section 5. It continues in Section 6, which focuses on the zeros of products of slice functions, and in Section 7, which is devoted to slice regular functions. Section 8 presents successful applications to the problem of factorizing motion polynomials.
2 The algebras of quaternions and dual quaternions
This section presents the algebra of quaternions H and the algebra of dual quaternions DH, along with relevant actions of their multiplicative groups.
2.1 Quaternions
Let H denote the real algebra of quaternions. Recall that it is obtained by endowing R4 with the multiplication operation defined on the standard basis {1,i,j,k} by
[TABLE]
and extended by distributivity to all quaternions q=x0+x1i+x2j+x3k. We set
[TABLE]
Re(q),Im(q) and ∣q∣ are called the real part, the imaginary part and the modulus of q, respectively. The quaternion
[TABLE]
is called the conjugate of q and it satisfies
[TABLE]
The inverse of any element q∈H∗:=H∖{0} is given by
[TABLE]
Two quaternions p,q commute if, and only if, Im(p)=cIm(q) for some c∈R. The commutative center
[TABLE]
of H coincides with R. Furthermore, for each pair of quaternions p,q, the standard scalar product between p and q equals 21(pqc+qpc); when p,q are orthogonal, i.e., pqc+qpc=0, we write p⊥q. For a more detailed discussion we refer the reader to [3].
2.2 Dual quaternions
Let DH denote the algebra of dual quaternions, which is an associative algebra over the real field R defined as
[TABLE]
with the following definitions for all h=h1+ϵh2,h′=h1′+ϵh2′:
Addition:
h+h′=(h1+h1′)+ϵ(h2+h2′);
Multiplication:
hh′=h1h1′+ϵ(h1h2′+h2h1′).
For each h=h1+ϵh2∈DH, we will refer to h1 as the primal part and h2 as the dual part of h.
Remark 2.1**.**
The center Z(DH) of DH coincides with the real (commutative) subalgebra DR:=R+ϵR, whose elements are called dual numbers.
By construction, ϵ2=0. As a consequence, we can make the following remark.
Remark 2.2**.**
A dual quaternion h∈DH∗=DH∖{0} is a zero divisor if, and only if, h=ϵh2, i.e., its primal part vanishes. In other words, the set of zero divisors in DH is ϵH∗ and the set of zero divisors in DR is ϵR∗. In particular, the product of any two zero divisors in DH equals zero.
Since DH is associative, an element admits a multiplicative inverse if, and only if, it is neither zero nor a zero divisor.
Remark 2.3**.**
For every h∈DH∖ϵH, the inverse of h is
[TABLE]
The algebra DH is a *-algebra with a *-involution called conjugation. For each h∈DH, it is defined as follows.
Conjugation:
hc=h1c+ϵh2c,
where h1c,h2c are the standard conjugates of the quaternions h1 and h2.
2.3 Imaginary units
Definition 2.4**.**
For every h∈DH⊇H, the trace of h and the (squared) norm of h are defined as
[TABLE]
Remark 2.5**.**
Trace and norm of any h∈DH are dual numbers, namely
[TABLE]
As a consequence, for all h∈H the trace t(h)=2Re(h) and the norm n(h)=∣h∣2 are real numbers.
Both for A=DH and for A=H, we call the elements of
[TABLE]
the imaginary units of A. The set
[TABLE]
is the unit 2-sphere in the 3-space of purely imaginary quaternions, while
[TABLE]
can be seen as the total space of the tangent bundle over SH. Indeed, any element J∈SDH splits into a primal part J1∈SH and a dual part J2 which equals mI for some m∈R and some I∈SH orthogonal to J1 and we may think of the set
[TABLE]
as the tangent plane to SH at J1.
Remark 2.6**.**
It is well-known that SH={J∈H:J2=−1}. As a consequence of this fact and of the equality h2=h12+ϵ(h1h2+h2h1), we also have that
[TABLE]
Both for A=DH and for A=H, the set of elements x that can be uniquely expressed as x=α+βJ for some α,β∈R, β>0 and J∈SA is the quadratic cone
where CJ=R+RJ is the unitary *-subalgebra generated by 1 any J. Clearly, CJ is isomorphic to the complex field and CI∩CJ=R for every I,J∈SA with I=±J. Moreover, QA splits as the disjoint union
[TABLE]
In other words, each x=α+βJ∈QA lies in a unique
[TABLE]
which is obtained by real translation and dilation from SA when β>0 and is a singleton {α} when β=0.
The sets QDH and SDH have interesting expressions in Cartesian coordinates. Let us consider the standard basis {1,i,j,k} of H. Any dual quaternion can be expressed by
[TABLE]
with rs∈R for all s∈[0,…,7]. The condition h1⊥h2 is equivalent to the equation of the Study quadric
[TABLE]
which will play an important role throughout the paper. The condition h2∈Im(H) is equivalent to the equation r4=0. Finally, h1∈H∖R if, and only if, r12+r22+r32=0. To sum up, QDH is determined by the conditions
[TABLE]
and it has dimension 6. In particular, the elements of SDH are determined by the equations
[TABLE]
For future reference, we point out that the intersection between S7 and the subspace r4=0 is the union
[TABLE]
where the first operand is a 4-space that intersects QDH along the axis R.
2.4 Conjugacy
Let us consider the multiplicative subgroup
[TABLE]
of the algebra DH. This group acts on DH as follows.
Definition 2.7**.**
We define
[TABLE]
We will also use the notation Ch(l):=C(h,l).
Thanks to the associativity of DH, it follows immediately that C is an action. Le us denote by [a,b]:=ab−ba the commutator of any a,b∈DH. As a consequence of Remark 2.3,
[TABLE]
Proposition 2.8**.**
For each h∈DH×, the map Ch maps SDH bijectively into itself. As a consequence, for all x∈QDH, the map Ch maps Sx bijectively into itself.
Proof.
Let us consider an arbitrary element h∈DH× and verify that Ch(J)∈SDH for every J∈SDH. Recalling the definition of SDH, it is sufficient to evaluate the norm and the trace of the element h−1Jh:
[TABLE]
[TABLE]
Moreover, we observe that (Ch)∣SDH:SDH→SDH is a bijection with inverse (Ch−1)∣SDH.
Finally, for all x=α+βJ∈QDH we have that (Ch)∣Sx is a transformation of Sx=α+βSDH as a consequence of the equality
[TABLE]
valid for all I∈SDH.
∎
The previous proposition implies that the action C on DH is not transitive. However, we can prove transitivity on a single Sx.
Proposition 2.9**.**
Fix x∈QDH. The action of DH× on Sx is transitive. If we quotient by DR×=DR∖ϵR then the action of DH×/DR× on Sx is faithful, but not free.
Proof.
Thanks to equality (4), it suffices to prove the theses for SDH (i.e., for x∈SDH).
Let us fix J∈SDH and prove that
[TABLE]
is a surjective map DH×→SDH. For any K1∈SH, the equality
[TABLE]
is fulfilled for h1=J1+K1 when K1=−J1; it is fulfilled for any h1∈SH with h1⊥J1 when K1=−J1. Now we want to prove that, for K1∈SH fixed and for the aforementioned h1∈H, the map
[TABLE]
from H to the tangent plane at K1 to SH in H is surjective. To this end, it suffices to prove that the real linear map
[TABLE]
has rank 2. Let us choose L1∈SH with L1⊥K1 and set M1:=K1L1=21[K1,L1]. Then the last displayed linear map transforms the vectors h1L1,h1M1 into the linearly independent vectors [K1,L1]=2M1,[K1,M1]=−2L1, as desired.
We have thus proven that the action of DH× on SDH is transitive.
This action is clearly not faithful, but it reduces to a faithful action if we quotient DH× with its center, which is the maximal normal subgroup of DH× included in all stabilizers of the action. As a consequence of the fact that the center of DH is DR, we conclude that the center of DH× is DR×=DR∖ϵR.
Finally, we can see that the quotient action is not free by considering the stabilizer of i. Indeed, since ih=hi is equivalent to h∈DC:=C+ϵC, we conclude that the stabilizer of i equals DC×/DR×, where DC×=DC∖ϵC.
∎
The next definition will be useful in the sequel. It extends C to DH×DH, although the first factor in the Cartesian product is no longer a group.
Definition 2.10**.**
We define
[TABLE]
and Ch(l):=C(h,l).
This peculiar extension is justified by the next remarks.
Remark 2.11**.**
For all x∈QDH, the extended C maps DH×Sx surjectively into Sx. Moreover, for all h∈DH, the equalities Re(Ch(x))=Re(x),Im(Ch(x))=Ch(Im(x)) hold.
Remark 2.12**.**
For all h,l∈DH, it holds
[TABLE]
This equality is obvious if h is invertible or h=0. If h∈ϵH∗, it follows by direct computation: hCh(l)=ϵh2(h2−1lh2)=l(ϵh2)=lh. In the last case, we actually have hl~=lh for all l~∈TCh(l1).
Remark 2.13**.**
Let us fix l∈DH. Then h2−1lh2 is a limit point of h−1lh as h1→0. Indeed,
[TABLE]
Moreover, l is a limit point of h−1lh as h→0 because limR∋t→0C(t,l)=limR∋t→0l=l.
2.5 Euclidean displacement
It is well-known, see [2, 15], that dual quaternions can be used to represent the group of proper rigid body transformations SE(3). Let us outline this characterization, starting with a few necessary tools.
First of all, R3 can be identified with 1+ϵIm(H) by seeing the vector (x1,x2,x3) as the dual quaternion 1+ϵ(x1i+x2j+x3k). Proper rigid body transformations can be obtained by means of an appropriate action on 1+ϵIm(H) of the following subgroup of the multiplicative group DH×:
[TABLE]
In order to define this action, let us introduce a new *-involution on DH (different from h↦hc):
*Alternate -involution:
h~=h1c−ϵh2c.
Now, to each h∈G, we can associate the transformation
[TABLE]
If we restrict to the case when t(h)∈R, i.e., h2∈Im(H), we find two special cases.
Translations.
If h∈R+ϵIm(H), then the previous transformation is a pure translation of R3 with translation vector 2h2h1−1. In particular, if h2=0 then the transformation is the identity.
Rotations.
If h∈QDH∖R, then the previous transformation is a pure rotation of R3. The rotation angle θ is determined by the equality cos(2θ)=∣h1∣Re(h1). The rotation axis has Plücker coordinates \big{(}\frac{\operatorname{Im}(h_{1})}{|\operatorname{Im}(h_{1})|},\frac{h_{2}}{|\operatorname{Im}(h_{1})|}\big{)}. In particular, if h=h1∈H∖R, then the rotation axis is the line through the origin parallel to Im(h1). This corresponds to the double covering of SO(3) by means of group of unitary quaternions described in [10, Theorem 3.17].
holds and the group of dual quaternions h with n(h)=1 is a double covering of SE(3).
3 The algebra of slice functions
In this section, we let the symbol A refer to H or DH indistinctly. We overview some material from [8, 9] concerning the theory of slice functions over the *-algebra A.
The *-algebra A can be seen as 2n-dimensional vector space over R for n=2 or n=4. The left multiplication by a element J∈SA induces a complex structure on A, thus there exist vectors e1J,…,en−1J∈A such that the set {1,J,e1J,Je1J,…,en−1J,Jen−1J} is a real vector basis, called a splitting basis of A associated to J. We consider on A the natural Euclidean topology and differential structure. The relative topology on each CJ with J∈SA clearly agrees with the topology determined by the natural identification between CJ and C, through the *-isomorphism
[TABLE]
Given a subset D of C, its circularization ΩD is defined as the following subset of QA:
[TABLE]
A subset of QA is termed circular if it equals ΩD for some D⊂C. For instance, given x=α+βJ∈QA the smallest circular set including x is
[TABLE]
From now on, we assume D to be invariant under complex conjugation z=α+βi↦α−βi. As a consequence, for each J∈SA the “slice” ΩJ:=Ω∩CJ is equivalent to D under the natural identification between CJ and C.
The class of A-valued functions we consider was defined in [8] by means of the complexified algebra AC=A⊗RC={x+ιy:x,y∈A,ι2=−1} of A, endowed with the following product:
[TABLE]
In this setting, A can be found as the subalgebra A+ι0 of AC. An isomorphic copy of C can be found as the subalgebra RC=R⊗RC=R+ιR of AC. Henceforth, we identify C (which includes D) with RC.
Remark 3.1**.**
We notice that if A=H, then Z(HC)=RC. Similarly, if A=DH then Z(DHC)=DRC.
In addition to the complex conjugation
[TABLE]
the *-involution on A induces a *-involution on AC, namely
[TABLE]
For all J∈SA, we can extend the previously defined map ϕJ:C→CJ to
[TABLE]
Let D be a subset of C and consider a function
[TABLE]
with A-valued components F1 and F2. The function F is called a stem function on D if F(z)=F(z) for every z∈D or, equivalently, if F1(z)=F1(z) and F2(z)=−F2(z) for every z∈D.
Definition 3.2**.**
A function f:ΩD⟶A is called a (left) slice function if there exists a stem function F:D⟶AC such that the diagram
[TABLE]
commutes for each J∈SA. In this situation, we say that f is induced by F and we write f=I(F). If F is RC-valued, then we say that the slice function f is slice preserving.
The algebraic structure of slice functions is described by the following proposition. A detailed proof can be found in [8].
Proposition 3.3**.**
*The stem functions D⟶AC form a *-algebra over R with pointwise addition (F+G)(z)=F(z)+G(z), multiplication (FG)(z)=F(z)G(z) and conjugation Fc(z)=F(z)c. This -algebra is associative and its center includes the subset of RC-valued stem functions. Let Ω=ΩD and consider the mapping
[TABLE]
*Besides the pointwise addition (f,g)↦f+g, there exist unique operations of multiplication (f,g)↦f⋅g and conjugation f↦fc on S(Ω) such that the mapping I is a *-algebra isomorphism. The *-algebra S(Ω) is associative and its center Z(S(Ω)) includes the -subalgebra SR(Ω) of slice preserving functions.
The next result describes the centers of the *-algebras of quaternionic and dual quaternionic slice functions.
Proposition 3.4**.**
*If A=H and Ω=ΩD⊆H, then Z(S(Ω))=SR(Ω). If A=DH and Ω=ΩD⊆QDH, then Z(S(Ω)) is the *-subalgebra SDR(Ω) of slice functions induced by DRC-valued stem functions, which properly includes the -subalgebra SR(Ω) of slice preserving functions.
Proof.
Let us prove the statement for A=DH: the case A=H is well-known and can be proven using the same technique.
Thanks to the previous proposition, it suffices to prove that the center of the *-algebra of stem functions D→DHC is the *-subalgebra of stem functions D→DRC.
If F:D→DRC is a stem function, then it commutes with any stem function G:D→DHC by Remark 2.1.
Conversely, suppose F:D→DHC to commute with all stem functions G:D→DHC. In particular, F commutes with all G=G1+ιG2 with G1≡a for a∈DH and G2≡0. If we fix z∈D, it follows that
The product f⋅g of two functions f,g∈S(Ω) is called slice product of f and g. If f belongs to the center of S(Ω) then, by direct inspection (f⋅g)(x)=f(x)g(x) for all x∈Ω. In general to compute f⋅g, one needs instead to compute FG and then f⋅g=I(FG). Similarly, for the slice conjugatefc of f=I(F) we compute fc=I(Fc). The normal function of f in S(Ω) is defined as
[TABLE]
Formulae to express the operations on slice functions without computing the corresponding stem functions can be given by means of two further operations. To each f∈S(Ω), we associate a function fs∘:Ω⟶A, called spherical value of f, and a function fs′:Ω∖R⟶A, called spherical derivative of f, by setting
[TABLE]
Spherical value and spherical derivative are slice functions, too: if f=I(F1+ιF2) then fs∘=I(F1) and fs′=I(F2) with F2(α+ιβ):=β−1F2(α+ιβ). Clearly,
[TABLE]
What is less obvious, but a consequence of the definition of slice functions, is the fact that fs∘,fs′ are constant on each Sx⊆Ω. As a consequence, f∈S(Ω) is uniquely determined by its restriction f∣ΩJ to any slice ΩJ (with J∈SA) of its domain.
We are now ready to state the aforementioned formulae for the operations on slice functions (see [9]): for all x∈Ω∖R,
[TABLE]
Remark 3.5**.**
Consider h,h′∈A then, by direct computation, n(h)=n(hc) and t(hh′)=t(h′h). As a consequence of formula (12), for all f∈S(Ω) we have that N(f)=N(fc). Moreover, since n and t take values in the center of A, it follows that N(f)∈Z(S(Ω)).
Definition 3.6**.**
A function f∈S(Ω) is tame if N(f)=N(fc) is slice preserving.
For A=H, all slice functions are tame. For A=DH, the tame elements of S(Ω) form a proper subset of S(Ω), which is closed under multiplication by [9, Remark 2.7].
Within the class of slice functions, we consider a special subclass having nice properties that recall those of holomorphic functions of a complex variable. Suppose Ω=ΩD to be open in QA, then for any J∈SA, Ω∩CJ is open in the relative topology of CJ; therefore, D itself is open. We let S0(Ω) and S1(Ω) denote the real vector spaces of slice functions on Ω induced by continuous stem functions and by stem functions of class C1, respectively. Now consider a function f=I(F)∈S1(Ω); for z=α+ιβ, set
[TABLE]
[TABLE]
Both ∂z∂F and ∂z∂F are still stem functions. They induce the slice functions
[TABLE]
Definition 3.7**.**
Let Ω be open in QA. A slice function f∈S1(Ω) is called slice regular if ∂xc∂f=0 in Ω. We denote by SR(Ω) the real vector space of slice regular functions on Ω. For each f∈SR(Ω), the slice regular function ∂x∂f is called the slice derivative (or complex derivative) of f.
The following lemma explains the connection between slice regularity and complex holomorphy. The proof of this result can be found in [8].
Lemma 3.8**.**
Suppose Ω=ΩD to be open in QA. Let J∈SA and let {1,J,e1J,Je1J,…,en−1J,Jen−1J} be an associated splitting basis of A. For f∈S1(Ω), let f0,…,fn−1:ΩJ⟶CJ be the C1 functions such that f∣ΩJ=∑l=0n−1flelJ, where e0J:=1. Then f is slice regular if, and only if, for each l∈{0,…,n−1}, fl is holomorphic from ΩJ to CJ, both equipped with the complex structure associated to left multiplication by J.
It has been proven in [8] that slice regularity is closed under addition, slice multiplication and slice conjugation. Thus, SR(Ω) is a *-subalgebra of S(Ω).
The next result, also from [8], provides a relevant class of examples of slice regular functions over DH, that will be particularly useful throughout the paper.
Proposition 3.9**.**
*Let DH[t] denote the -algebra of polynomials ∑n=0dtnan over dual quaternions, with the standard operations
[TABLE]
*Let denote H[t] the *-subalgebra of quaternionic polynomials. Mapping each polynomial P(t) into the function P∣QDH:QDH→DH defines an injective *-algebra homomorphism DH[t]→SR(QDH). Moreover, the inclusion H[t]⟶SR(H) is an injective -algebra homomorphism.
If we take into account that DR is both the center of DH and the subspace of the points of DH preserved by conjugation, a direct inspection in the previous definition allows the following remark.
Remark 3.10**.**
Let P(t),Q(t)∈DH[t], let R(t):=P(t)⋅Q(t) and let us evaluate these three polynomials at a point h∈DH. While (P+Q)(h)=P(h)+Q(h), the equalities
[TABLE]
are only guaranteed if h belongs to DR or if P(t) belongs to DR[t].
As explained in [11], polynomials over dual quaternions are particularly relevant for applications in kinematics. Let us sketch the relevant construction. We point out out that we are considering polynomials with coefficients on the right-hand side, as opposed to the convention adopted in [11]. The polynomials ∑n=0dtnan and ∑n=0dantn coincide when evaluated at real points, but not when evaluated at other points of DH.
Definition 3.11**.**
Consider a polynomial P(t)=∑n=0dtnan∈DH[t] having degree d. The polynomial N(P)(t):=P(t)⋅Pc(t)∈DR[t] is called the norm of P(t). If N(P)(t) belongs to R[t] and the leading coefficient ad of P(t) belongs to DH×, then P(t) is called a motion polynomial.
If P(t) is a motion polynomial, then we can make the following observations:
•
for each t0∈R, the value h=P(t0) has the property n(h)=N(P)(t0)∈R, whence it belongs to the Study quadric S7;
•
the leading coefficient of N(P)(t) is the 2d-th coefficient, namely n(ad)∈R∗; for any t0∈R that is not a root of N(P)(t), the value h=P(t0) is an element of G=DH×∩S7.
As a consequence, for any t0∈R that is not a root of N(P)(t) it is possible to consider the proper rigid body transformation (5) with h=P(t0):
[TABLE]
If we fix a point, say 1+ϵx0, then its trajectory
[TABLE]
will be a rational curve.
The following remark relates motion polynomials to tame functions and their properties.
Remark 3.12**.**
Let us consider a slice regular function f defined as f:=P∣QDH for some P(t)=∑n=0dtnan∈DH[t] having degree d. Then its normal function N(f) coincides with N(P)∣QDH. We can draw the following consequences.
•
P(t)* is a motion polynomial if, and only if, f is tame and ad∈DH×.*
•
If P(t) is a motion polynomial then f cannot be a zero divisor in S0(QDH) by **[9, Proposition 5.18]** and P(t) cannot be a zero divisor in DH[t].
4 Primal part function
In this section we associate to each DH-valued slice regular function f an H-valued slice regular function, called the primal part of f. This notion extends the analogous notion defined for polynomials in [11]. We begin with some preliminary definitions and results.
Definition 4.1**.**
Let π:DH⟶H be the function that maps a dual quaternion into its primal part, i.e., π(h)=π(h1+ϵh2)=h1.
Lemma 4.2**.**
*The function π:DH⟶H is a surjective real *-algebra homomorphism. Its extension πC:DHC⟶HC, defined as πC(x+ιy)=π(x)+ιπ(y), is a surjective complex -algebra homomorphism.
Proof.
By construction, π is a surjective R-linear map. Moreover, for all h,l∈DH, the following equalities hold:
[TABLE]
Thus, π is a real *-algebra homomorphism. Moreover,
[TABLE]
Thus, πC is a surjective complex *-algebra homomorphism.
∎
We can easily study the effect of π on the the map C of Definition 2.10.
Remark 4.3**.**
Let h,l∈DH. If h is invertible, equality (3) implies that
[TABLE]
If h∈ϵH∗, then
[TABLE]
If h=0 then both of the previous formulae are true.
Proposition 4.4**.**
For each stem function F:D⟶DHC, the function πF:=πC∘F is a stem function
[TABLE]
*Moreover, if F is holomorphic then πF is holomorphic too. Finally, the map F↦πF is a -algebra homomorphism.
Proof.
By definition, for every z∈D
[TABLE]
Thus, πF is a stem function. The fact that F↦πF is a *-algebra homomorphism follows at once from the previous lemma.
By applying again the previous lemma, we observe that πC is a C-linear map, whence it coincides with its differential. Thus, if F is holomorphic then
[TABLE]
whence πF is a holomorphic function, too.
∎
The previous result allows us to give the next definition and to derive the subsequent corollary.
Definition 4.5**.**
Let Ω=ΩD⊆QDH and let f=I(F) be a function in S(Ω). We define the primal part of f as the quaternionic slice function πf:=I(πF)∈S(Ω∩H).
Corollary 4.6**.**
Let Ω=ΩD⊆QDH. The map
[TABLE]
*is a -algebra homomorphism. If Ω is open in QDH, the homomorphism maps slice regular functions into slice regular functions.
The entire construction is designed to satisfy the following property:
Remark 4.7**.**
Let Ω=ΩD⊆QDH and consider a function f∈S(Ω), induced by the stem function F. For all x=α+βJ (whence x1=α+βJ1) and for z=α+ιβ, it holds
[TABLE]
As a consequence,
[TABLE]
The previous remark will be useful in Section 5 to relate the zeros of a slice function over dual quaternions with the zeros of its primal part.
Remark 4.8**.**
*We can repeat the construction by using any -algebra homomorphism ψ instead of π, indeed
[TABLE]
where ψ(J)∈SDH because n(ψ(J))=ψ(n(J))=ψ(1)=1 and t(ψ(J))=ψ(t(J))=ψ(0)=0. In our case ψ=π, the map J↦ψ(J) for J∈SDH is the natural projection of SDH onto SH.
Let us relate our construction with the concept of primal part of a polynomial considered in [11] (again, with a different convention about the side of the coefficients).
Definition 4.9**.**
The map primal:DH[t]→H[t] is defined as
[TABLE]
Remark 4.10**.**
Let us consider P(t)∈DH[t] and evaluate it at h=h1+ϵh2∈DH. Then
[TABLE]
Proposition 4.11**.**
*If f∈S(QDH) is defined as f:=P∣QDH for some P(t)∈DH[t], then πf=primal(P). As a consequence, primal:DH[t]→H[t] is a -algebra homomorphism.
Proof.
If P(t)=∑n=0dtnan then f(x)=∑n=0dxnan. By Lemma 4.2 and Remark 4.7, for all x∈DH it holds
[TABLE]
Since π:DH→H is surjective, we conclude that
[TABLE]
for all w∈H, which is our first statement.
The second statement now follows from Proposition 3.9 and Corollary 4.6.
∎
Remark 4.12**.**
If P(t)∈DH[t] is a motion polynomial of degree d, then primal(P)(t)∈H[t] has degree d. The norm of primal(P)(t) is a 2d-degree real polynomial and it coincides with N(P)(t).
5 Zeros of slice functions
In this section, we describe some algebraic and geometric properties of the zero set
[TABLE]
of a slice function f∈S(Ω) with Ω=ΩD. Moreover, we study how this set is related to the zero sets of fc,πf and N(f).
Before proceeding towards the main results, let us establish two useful equalities.
Lemma 5.1**.**
Let f∈S(Ω). For all x∈Ω∩R, it holds fc(x)=f(x)c. For all x∈Ω∖R and all y∈Sx, it holds fc(C(fs′(x),yc))=f(y)c.
Proof.
If x∈Ω∩R, then fc(x)=(fc)s∘(x)=fs∘(x)c=f(x)c.
Suppose, instead, x∈Ω∖R and y∈Sx. We claim that C(h,l)hc=hcl for all h∈DH and we compute
[TABLE]
which proves our thesis.
Our claim can be easily derived from the equality hC(h,l)=lh proven in Remark 2.12. If h is invertible, it suffices to multiply each hand of the equality by hc both on the left and on the right and to divide it by the dual number n(h). If, instead, h=ϵh2, it suffices to multiply each hand of the equality by h2c both on the left and on the right and to divide it by the real number n(h2).
∎
We are now ready to study the zero set V(f) and its relation to V(fc). For each x∈Ω∖R, let us set the notation
[TABLE]
for the tangent plane to the 2-sphere Sx∩H at x1. This is a consistent extension of the notation TJ1 we have already set for J1∈SH.
Theorem 5.2**.**
Let f∈S(Ω). If x∈Ω∩R, then Sx={x} is included either in both V(f) and in V(fc) or in none of the two. If instead x∈Ω∖R, then one of the following properties holds:
V(f)* does not intersect Sx;*
2. 2.
V(f)∩Sx={y}, fs′(x) is invertible and y=Re(x)−fs∘(x)fs′(x)−1;
3. 3.
V(f)∩Sx=Ty1* for some y1∈Sx∩H and fs′(x),fs∘(x)∈ϵH∗;*
4. 4.
V(f)* includes Sx and fs′(x)=fs∘(x)=0.*
In each of the aforementioned cases, respectively:
V(fc)* does not intersect Sx;*
2. 2.
V(fc)∩Sx={C(fs′(x),yc)};
3. 3.
V(fc)∩Sx=TC(fs′(x),y1c);
4. 4.
V(fc)* includes Sx.*
Proof.
Our statement for x∈Ω∩R follows from the equality fc(x)=f(x)c. Now let us suppose x∈Ω∖R. For all y∈Sx, the following decomposition holds:
[TABLE]
If V(f)∩Sx=∅, then consider y∈V(f)∩Sx.
•
If fs′(x) is invertible then, by equality (15), we have
[TABLE]
whence y=Re(x)−fs∘(x)fs′(x)−1.
•
If instead fs′(x)∈ϵH∗ then, starting again from equality (15), we have
[TABLE]
In this case, fs∘(x)∈ϵH∗. Moreover, for every z∈Sx, we decompose f(z) as follows:
[TABLE]
By Remark 2.2, we conclude that f(z)=0 if and only if, z−y∈ϵH. This is a same as z1=y1, i.e., z∈Ty1.
•
Finally, if fs′(x)=0, then f≡fs∘(x) in Sx. Since f(y)=0, we immediately conclude that fs∘(x)=0 and V(f)⊇Sx.
If x∈Ω∩R then Sx={x}. It is straightforward that either V(f)⊇Sx or V(f)∩Sx=∅, depending on whether f(x) vanishes or not.
We now prove the statement concerning fc. If x∈Ω∩R, whence Sx={x}, then either V(fc)⊇Sx or V(fc)∩Sx=∅, depending on whether fc(x)=f(x)c vanishes or not. Now let x∈Ω∖R. First suppose V(f)∩Sx includes a point y. Then, by Lemma 5.1, V(fc)∩Sx includes the point C(fs′(x),yc). According to what we have proven so far, there are three possibilities.
•
If V(f)∩Sx={y}, then (fc)s′(x)=fs′(x)c is invertible and V(fc)∩Sx={C(fs′(x),yc)}.
•
If V(f)∩Sx=Ty1, then (fc)s′(x)=fs′(x)c∈ϵH∗ and V(fc)∩Sx=TC(fs′(x),y1c).
•
If V(f)⊇Sx then (fc)s′(x)=fs′(x)c=0 and V(fc)⊇Sx.
The only remaining case is V(f)∩Sx=∅. In this case, fc cannot have any zero z∈Sx: otherwise, by Lemma 5.1, f=(fc)c would vanish at C(fs′(x)c,zc)∈Sx and we would obtain a contradiction.
∎
The next result connects the zero set of a function f∈S(Ω) with the zero set of its primal part.
Proposition 5.3**.**
Let f∈S(Ω). If x∈Ω∩R, then f(x)∈ϵH, if and only if, πf(x)=0. If x∈Ω∖R, then there are three possibilities:
all values of f in Sx are invertible and V(πf) does not intersect Sx∩H;
2. 2.
f* maps exactly one tangent plane Ty1 into ϵH and V(πf)∩Sx={y1};*
3. 3.
f* maps Sx into ϵH and V(πf) includes Sx∩H.*
Moreover:
if V(f)∩Sx={y}, then V(πf)∩Sx={y1};
2. 2.
if V(f)∩Sx=Ty1 or V(f)⊇Sx, then V(πf) includes Sx∩H.
Proof.
By Remark 4.7, for all w∈Ω it holds f(w)∈ϵH if, and only if, πf(w1)=0.
If w=x∈Ω∩R, it holds w1=x and the first statement follows.
We can prove the second statement as follows. If x∈Ω∖R, then πf may vanish:
at no point of Sx∩H;
2. 2.
at exactly one point of Sx∩H, say y1;
3. 3.
at all points of Sx∩H.
In case 1. it holds f(w)∈ϵH for all w∈Sx. In case 2. it holds f(y1+ϵγ)∈ϵH for all γ∈Im(H),γ⊥Im(y1), while f(w)∈ϵH for all w∈Sx with w1=y1. In case 3. it holds f(w)∈ϵH for all w∈Sx.
Let us now prove the third statement. If V(f)∩Sx={y} with y∈Sx then πf(y1)=0. Moreover, if πf had another zero z1∈Sx∩H then (πf)s′(x1)=π(fs′(x)) would vanish and fs′(x) would be a zero divisor, contradicting Theorem 5.2. If, instead, V(f)∩Sx=Ty1 or V(f)⊇Sx, then by the same theorem, fs′(x) and fs∘(x) belong to ϵH. As a consequence, f maps Sx into ϵH and V(πf)⊇Sx∩H.
∎
We now study the zero set of the normal function N(f), taking full advantage of the properties of DH. We begin by establishing that it is circular.
Proposition 5.4**.**
Let f∈S(Ω). Then V(N(f)) is circular, i.e., V(N(f))∩Sx=∅ implies V(N(f))⊇Sx.
Proof.
Suppose Sx=α+βSDH. For the sake of simplicity, we denote fs∘(x) and βfs′(x) by a and b, respectively. For all J∈SDH, by formula (12),
[TABLE]
If N(f)(α+βI)=0 then n(a)−n(b),I1t(a1b1c) and ϵI1t(a1b2c+a2b1c)+ϵI2t(a1b1c) vanish, separately. This is, in turn, equivalent to n(a)−n(b)=t(a1b1c)=t(a1b2c+a2b1c)=0. If this is the case, then N(f)(α+βJ)=0 independently of J and V(N(f))⊇Sx.
∎
Our next aim is studying the relation between V(N(f)) and V(f). The next lemma will be useful to this end, because it connects the spherical derivative and the spherical value of N(f) to the values of f.
Lemma 5.5**.**
Let f∈S(Ω) and let x∈Ω∖R. For all y∈Sx, the following equalities hold.
[TABLE]
Proof.
For all y=α+βJ∈Sx, we have
[TABLE]
where as usual we denote fs∘(x) and βfs′(x) by a and b, respectively. By Equation (12), N(f)s′(x)=t(abc) and N(f)s∘(x)=n(a)−n(b). Recalling that DR is the center of DH and noticing that it is invariant under *-involution, the following equalities hold:
[TABLE]
Moreover,
[TABLE]
We are now able to state and prove the following theorem.
Theorem 5.6**.**
Let f∈S(Ω) with Ω=ΩD. Then
[TABLE]
Moreover, for each x∈Ω∖R the normal function N(f) vanishes (identically) in Sx, if and only if, either f has a unique zero in Sx or πf vanishes identically in Sx∩H.
Proof.
As a first step, let us prove that V(N(f)) includes ⋃V(f)∩Sx=∅Sx. If y∈V(f)∩Sx then, by Lemma 5.5, N(f)s′(x)=N(f)s∘(x)=0, whence V(N(f))⊇Sx.
As a second step, we prove that V(N(f)) includes ⋃Sx∩H⊆V(πf)Sx. If πf≡0 in Sx∩H then a=fs∘(x) and b=βfs′(x) have a1=0=b1. In other words, a and b belong to ϵH, whence n(a)=n(b)=abc=0. As a consequence, for all y=α+βJ∈Sx, the expression
[TABLE]
vanishes.
As a third step, let us take any Sx contained in
[TABLE]
and prove that Sx includes exactly one zero of f. We observe that N(πf)=πN(f) by Corollary 4.6, whence V(N(πf))⊇Sx∩H. By [3, Proposition 3.9], the function πf has a zero y1=α+βJ1 in Sx∩H, a zero which is unique because we have assumed πf not to vanish identically in Sx∩H. We complete our proof by finding a unique zero of f in Ty1. If a=fs∘(x) and b=βfs′(x) (whence β(πf)s′(x)=b1 and (πf)s∘(x)=a1=−J1b1), we have to prove that there exists a unique γ∈Im(H) with γ⊥J1 such that
[TABLE]
This happens if, and only if, a2b1c+J1b2b1c is an element of Im(H), orthogonal to J1, i.e.,
[TABLE]
After recalling that t(pq)=t(qp) and t(pqc)=t(qpc) for all p,q∈H, we can observe that
[TABLE]
is the dual part of t(abc) and that
[TABLE]
is the dual part of n(a)−n(b). Since N(f) vanishes identically in Sx, we know that n(a)−n(b)=0=t(bac). Thus, system (17) is fulfilled and our proof is complete.
∎
We conclude the section with some examples that illustrate the previous results.
Example 5.7**.**
For all x∈QDH, let
[TABLE]
Then fc=f and N(f)(x)=(x2+1)2. Moreover, πf(x1)=x12+1. Thus,
[TABLE]
Example 5.8**.**
For all x∈QDH, let
[TABLE]
It holds πf(x1)=x12+1 and N(f)=(x2+1)(x2+1−2xϵ), whence
[TABLE]
It is easy to observe that i is a zero of f. Moreover, since f(x) coincides with −xϵ+ϵi for all x∈SDH, we conclude that fs′(i)=−ϵ. Thus,
[TABLE]
where we took into account the fact that C(−ϵ,ic)=−i.
Example 5.9**.**
For all x∈QDH, let
[TABLE]
By direct computation, fc(x)=x−1+i+ϵj, N(f)(x)=x2−2x+2 and πf(x1)=x1−1−i. As a consequence,
[TABLE]
Example 5.10**.**
On QDH, let
[TABLE]
By direct computation, fc≡−ϵi, N(f)≡0 and πf≡0. Thus,
[TABLE]
6 Zeros of slice products
This section describes in great detail the zero set of the slice product of two slice functions over DH. We begin with a result that expresses the values of a slice product as products of values of its two factors.
Theorem 6.1**.**
Let f,g∈S(Ω) and fix x∈Ω. If x∈Ω∩R, then (f⋅g)(x)=f(x)g(x).
Suppose instead x∈Ω∖R. If y,z∈Sx fulfill one of the following (mutually equivalent) conditions:
yf(y)−f(y)z=0;
2. 2.
zfc(z)−fc(z)y=0;
then
[TABLE]
Condition 1. is equivalent to z=C(f(y),y), when f(y) is invertible; it is equivalent to z∈TC(f(y),y1), when f(y) is a zero divisor; it is automatically fulfilled when f(y)=0. Similarly, condition 2. is equivalent to y=C(fc(z),z) when fc(z) is invertible; it is equivalent to y∈TC(fc(z),z1) when fc(z) is a zero divisor; it is automatically fulfilled when fc(z)=0.
whenever Im(y)f(y)=f(y)Im(z), which is equivalent to condition 1.
•
We can prove the equivalence between conditions 1. and 2., as follows. For the sake of simplicity, we denote fs∘(x) and βfs′(x) by a and b, respectively. Supposing y=α+βJ,z=α+βK for some J,K∈SDH, it holds f(y)=a+Jb and, by formula (10), fc(z)=ac+Kbc. Condition 1. is equivalent to
[TABLE]
The last equality is equivalent to
[TABLE]
which is, in turn, equivalent to condition 2.
•
The characterization of conditions 1. and 2. follows directly from Remark 2.12.
∎
Corollary 6.2**.**
Let f,g∈S(Ω). The formulae
[TABLE]
hold for all y,z∈Ω. As a consequence, V(f⋅g) includes both the zero set V(f) of f and the set {C(fc(z),z):z∈V(g)}.
To deepen the study of V(f⋅g), we recall Theorem 5.2: the different types of zeros of f⋅g correspond to different properties of its spherical derivative (f⋅g)s′. Therefore, it is useful to establish the next result.
Lemma 6.3**.**
Let f,g∈S(Ω) and let x∈Ω∖R.
If y∈V(f)∩Sx, then (f⋅g)s′(x)=fs′(x)g(C(fs′(x),yc)).
2. 2.
If z∈V(g)∩Sx, then (f⋅g)s′(x)=(fc(z))cgs′(x).
3. 3.
If y∈V(f)∩Sx and z∈V(g)∩Sx, then
(f\cdot g)^{\prime}_{s}(x)=f^{\prime}_{s}(x)\big{(}\mathscr{C}(f^{\prime}_{s}(x),y^{c})-z\big{)}g^{\prime}_{s}(x).
Proof.
According to formula (11), it holds (f⋅g)s′=fs∘gs′+fs′gs∘.
If f has a zero y in Sx, then fs∘(x)=−Im(y)fs′(x) and
If g has a zero z in Sx, then gs∘(x)=−Im(z)gs′(x) and
[TABLE]
Finally, if both f and g have zeros in Sx, namely y and z, then
[TABLE]
where we have used again Remark 2.12, along with the equalities −Im(y)=Im(yc) and Re(yc)=Re(y)=Re(z).
∎
Theorem 6.4**.**
Let f,g∈S(Ω). If x∈Ω∩R then x∈V(f⋅g) is equivalent to x∈V(f)∪V(g)∪(V(πf)∩V(πg)). If, instead, x∈Ω∖R, then the following statements hold.
If V(f)⊇Sx or V(g)⊇Sx, then V(f⋅g)⊇Sx.
2. 2.
If V(f)∩Sx=Ty1 and V(g)∩Sx=Tz1, then V(f⋅g)⊇Sx.
3. 3.
Suppose V(f)∩Sx=Ty1 and V(g)∩Sx={z}.
•
If z1=C(fs′(x),y1c), then V(f⋅g)⊇Sx.
•
Otherwise, V(f⋅g)∩Sx=Ty1.
4. 4.
Suppose V(f)∩Sx=Ty1 and V(g)∩Sx=∅.
•
If C(fs′(x),y1c)∈V(πg), then V(f⋅g)⊇Sx.
•
Otherwise, V(f⋅g)∩Sx=Ty1.
5. 5.
Suppose V(f)∩Sx={y} and V(g)∩Sx=Tz1.
•
If z1=C(πfs′(x1),y1c), then V(f⋅g)⊇Sx.
•
Otherwise, V(f⋅g)∩Sx=Ty1.
6. 6.
Suppose V(f)∩Sx={y} and V(g)∩Sx={z}.
•
If z=C(fs′(x),yc), then V(f⋅g)⊇Sx.
•
If z is a point of TC(πfs′(x1),y1c) other than C(fs′(x),yc), then V(f⋅g)∩Sx=Ty1.
•
Otherwise, V(f⋅g)∩Sx={y}.
7. 7.
Suppose V(f)∩Sx={y} and V(g)∩Sx=∅.
•
If C(πfs′(x1),y1c)∈V(πg), then V(f⋅g)∩Sx=Ty1.
•
Otherwise, V(f⋅g)∩Sx={y}.
8. 8.
Suppose V(f)∩Sx=∅ and V(g)∩Sx=Tz1.
•
If z1∈V(πfc), then V(f⋅g)⊇Sx.
•
Otherwise, V(f⋅g)∩Sx=TC(πfc(z1),z1).
9. 9.
Suppose V(f)∩Sx=∅ and V(g)∩Sx={z}.
•
If z1∈V(πfc), then V(f⋅g)∩Sx=TC(fc(z),z1).
•
Otherwise, V(f⋅g)∩Sx={C((fc(z),z)}.
10. 10.
If V(f)∩Sx=∅ and V(g)∩Sx=∅, then
[TABLE]
Proof.
The first statement follows follows from the fact that
[TABLE]
for all x∈Ω∩R. To prove the second statement, we proceed step by step by step and we repeatedly apply Theorem 5.2, Corollary 6.2 and Lemma 6.3.
If V(f)⊇Sx then V(f⋅g)⊇Sx. If V(g)⊇Sx, take any z∈Sx. Then V(f⋅g)∩Sx includes the point C(fc(z),z) and gs′(x)=0. Since
[TABLE]
it follows that V(f⋅g)⊇Sx.
2. 2.
If V(f)∩Sx=Ty1 and V(g)∩Sx=Tz1, then V(f⋅g)⊇Ty1 and fs′(x),gs′(x)∈ϵH∗. Since
[TABLE]
we immediately conclude that V(f⋅g)⊇Sx.
3. 3.
If V(f)∩Sx=Ty1 and V(g)∩Sx={z}, then V(f⋅g)⊇Ty1, fs′(x)∈ϵH∗ and gs′(x) is invertible. Either V(f⋅g)⊇Sx or V(f⋅g)∩Sx=Ty1, depending on whether (f⋅g)s′(x) vanishes or not. But
[TABLE]
vanishes if, and only if, the second factor belongs to ϵH. This happens if, and only if, z1=C(fs′(x),y1c).
4. 4.
If V(f)∩Sx=Ty1 and V(g)∩Sx=∅, then V(f⋅g)⊇Ty1 and fs′(x)∈ϵH∗. Either V(f⋅g)⊇Sx or V(f⋅g)∩Sx=Ty1, depending on whether (f⋅g)s′(x) vanishes or not. The expression
[TABLE]
vanishes if, and only if, the second factor belongs to ϵH. This happens if, and only if, C(fs′(x),y1c)∈V(πg).
5. 5.
If V(f)∩Sx={y} and V(g)∩Sx=Tz1, then V(f⋅g)⊇{y}, fs′(x) is invertible and gs′(x)∈ϵH∗. The expression
[TABLE]
vanishes if, and only if, the second factor belongs to ϵH; this is, in turn, equivalent to z1=C(πfs′(x1),y1c). If this is the case, then V(f⋅g)⊇Sx. Otherwise, (f⋅g)s′(x)∈ϵH∗ and V(f⋅g)∩Sx=Ty1.
6. 6.
If V(f)∩Sx={y} and V(g)∩Sx={z}, then V(f⋅g)⊇{y} and fs′(x),gs′(x) are invertible. By the equality
[TABLE]
there are three possibilities. If z=C(fs′(x),yc) then (f⋅g)s′(x)=0 and V(f⋅g)⊇Sx. If not, but if we still have z1=C(πfs′(x1),y1c), then (f⋅g)s′(x)∈ϵH∗ and V(f⋅g)∩Sx=Ty1. Otherwise, (f⋅g)s′(x) is invertible and V(f⋅g)∩Sx={y}.
7. 7.
If V(f)∩Sx={y} and V(g)∩Sx=∅, then V(f⋅g)⊇{y} and fs′(x) is invertible. Either V(f⋅g)∩Sx=Ty1 or V(f⋅g)∩Sx={y}, depending on whether or not
[TABLE]
belongs to ϵH∗. This happens if, and only if, g(C(fs′(x),yc))∈ϵH∗, which is, in turn, equivalent to C(πfs′(x1),y1c)∈V(πg).
8. 8.
If V(f)∩Sx=∅ and V(g)∩Sx=Tz1, then V(f⋅g)⊇{C(fc(z1),z1)} and gs′(x)∈ϵH∗. The expression
[TABLE]
vanishes if, and only if, fc(z1)∈ϵH; this is, in turn, equivalent to z1∈V(πfc). If this is the case, then V(f⋅g)⊇Sx. Otherwise, (f⋅g)s′(x)∈ϵH∗ and V(f⋅g)∩Sx=TC(πfc(z1),z1).
9. 9.
If V(f)∩Sx=∅ and V(g)∩Sx={z}, then V(f⋅g)∩Sx⊇{C(fc(z),z)} and gs′(x) is invertible. The expression
[TABLE]
belongs to ϵH∗ if, and only if, the first factor does. This is, in turn, equivalent to z1∈V(πfc). If this is the case, then V(f⋅g)∩Sx=TC(fc(z),z1). Otherwise, V(f⋅g)∩Sx={C(fc(z),z)}.
10. 10.
If V(f)∩Sx=∅ and V(g)∩Sx=∅, then for all y∈Sx it holds
[TABLE]
This product vanishes if, and only if, both factors belong to ϵH∗. This happens if, and only if, y1∈V(πf) and C(f(y),y1)∈V(πg).
∎
Examples 6.5** (Case 6.).**
Fix J∈SDH. For all x∈QDH, let
[TABLE]
Then V(f)={i}, V(g)={J} and fs′(i)=1=πfs′(i). Let us determine the zeros of the product
[TABLE]
Since V(πf)={i} and V(πg)={J1} are both included in SDH, the zero set V(f⋅g) must also be included in SDH. There are three possibilities.
•
If J=C(1,−i)=−i then V(f⋅g)=SDH.
•
If J is a point of T−i other than −i, then V(f⋅g)=Ti.
•
Otherwise, V(f⋅g)={i}.
Example 6.6** (Case 9.).**
For all x∈QDH, let
[TABLE]
The function f has no zeros in QDH, while V(g)={−i−ϵj}⊂SDH. Let us determine the zeros of the product
[TABLE]
For SDH, we compute fc(x)=2x+2i−ϵi, whence V(πfc)={−i}, and
[TABLE]
Thus, V(f⋅g)∩SDH=T53i−54j. Since V(πf)={i} and V(πg)={−i} are both included in SDH, there are no other possible zeros of f⋅g and
[TABLE]
Example 6.7** (Cases 4. and 10.).**
For all x∈QDH, let
[TABLE]
Then V(f)=T2k and g never vanishes in QDH. Let us determine the zero set of the product
[TABLE]
We first focus on 2SDH∋2k: we compute fs′(2k)=ϵ and C(ϵ,2k)=2k, which does not belong to V(πg)={i}; we conclude that
[TABLE]
Since V(πf)=H and V(πg)={i}, the other possible zeros y∈QDH of f⋅g are determined by the equation C(f(y),y1)=i, which is equivalent to
[TABLE]
We conclude that
[TABLE]
7 Zeros of slice regular functions
In this section, we study the discreteness of the zeros of slice regular functions. We also show that they can be factored out, a fact which will be particularly useful for subsequent applications to motion polynomials.
Our discreteness results require two preliminary definitions.
Definition 7.1**.**
If Ω=ΩD, where D is an open connected subset of C that intersects the real line R and is preserved by complex conjugation, then Ω is called a slice domain. If Ω=ΩD, where D is an open subset of C that does not intersect R and has two connected components swapped by complex conjugation, then Ω is called a product domain.
It is not restrictive to study slice regular functions on either slice domains or product domains, because all circular open subsets of DH are unions of slice domains and product domains. We now come to the announced discreteness results, where CJ+:={α+Jβ:α,β∈R,β>0} and CJ−:={α+Jβ:α,β∈R,β<0}.
Theorem 7.2**.**
Let Ω=ΩD, where D is an open subset of C, and let f∈SR(Ω). If Ω is a slice domain then either
(1)
for each J∈SDH, the intersection V(f)∩CJ is closed and discrete in ΩJ; or
(2)
f* vanishes identically.*
If Ω is a product domain then, in addition to cases (1) and (2), there is one further possibility:
(3)
There exists J∈SDH such that
[TABLE]
If πfs′≡0 then J is unique; if instead πfs′≡0, then the imaginary units having the same property are exactly the elements of TJ1. For all other K∈SDH, the intersection V(f)∩CK is closed and discrete in ΩK.
Proof.
After applying [9, Theorem 4.11], it only remains to study the uniqueness of J in case (3) by finding the solutions J′ of the equation (J′−J)fs′≡0. If at least one value of fs′ is invertible (that is, πfs′≡0), then the only solution is J′=J. Otherwise, the image of fs′ is included in the set ϵH and πfs′≡0. In this situation, the solutions are the J′∈SDH such that J′−J∈ϵH or, equivalently, the elements of TJ1.
∎
Theorem 7.3**.**
Let f∈SR(Ω), where either
•
Ω⊆QDH* is a slice domain and πf≡0; or*
•
Ω⊆QDH* is a product domain and N(πf)≡0*
Then the zero set V(f) is a union of singletons {y}, tangent planes Ty1 or “spheres” Sy, each isolated from the rest of V(f). If, moreover, f=P∣QDH for some polynomial P(t)∈DH[t], then the union is finite.
Proof.
By Theorem 5.2, the zero set V(f) is a union of singletons {y}, tangent planes Ty1 or “spheres” Sy. Moreover, each singleton {y} included in V(f) corresponds to a singleton {y1} included in V(πf); each tangent plane Ty1 or “sphere” Sy corresponds to a 2-sphere Sy1 included in V(πf).
We claim that, under either of our two hypotheses, V(πf) consists of isolated points or isolated 2-spheres of the form Sy1, whence the first statement immediately follows.
As for the second statement, the equality f=P∣QDH implies that πf coincides with the quaternionic polynomial primal(P)(t), which has finitely many isolated zeros or 2-spheres of zeros (see, e.g., [3, Proposition 3.30]).
We prove our previous claim as follows.
•
If Ω is a slice domain in QDH then Ω∩H is a slice domain in H. If πf≡0, then V(πf) consists of isolated points y1 or isolated 2-spheres of the form Sy1 by [3, Theorem 3.12].
•
If Ω is a product domain in QDH then Ω∩H is a product domain in H. If N(πf)≡0, then, according to [8, page 1681], V(πf) consists of isolated points y1 or isolated 2-spheres of the form Sy1.
This completes the proof.
∎
All examples provided in Sections 5 and 6 are slice regular functions on slice domains and they have zero sets of the types described in the last two theorems. We now provide a few pathological examples that do not fulfill the hypotheses of these theorems.
Example 7.4**.**
Consider the slice regular functions
[TABLE]
on the product domain QDH∖R=ΩC+: for each J∈SDH, it holds f∣CJ+≡1+Ji and g∣CJ+≡(1+Ji)ϵ. As a consequence,
[TABLE]
We point out that πfs′≡0≡πf and N(πf)≡0, while πgs′,πg,N(πg) all vanish identically.
We now proceed forward, aiming at factoring out the zeros of slice regular functions.
Definition 7.5**.**
Let Ω=ΩD⊆QDH. For f,h∈S(Ω), we say that h is a left factor of f if there exists g∈S(Ω) such that
[TABLE]
in Ω. If this is the case, we also say that h(x)dividesf(x)on the left and write h(x)∣f(x). If h(x) does not divide f(x) on the left, we write h(x)∣f(x).
We adopt analogous terminologies and notations for other algebras, such as the algebra of quaternionic slice functions S(Ω∩H) and the algebra of polynomials over dual quaternions DH[t].
Theorem 7.6**.**
Let f∈SR(Ω) with Ω=ΩD and let y∈Ω.
•
The zero set V(f) includes y if, and only if, x−y∣f(x). If this is the case, then we have x2−xt(y)+n(y)∣N(f)(x).
If y∈Ω∖R then
•
V(f)⊇Ty1* if, and only if, x−y∣f(x) and x12−x1t(y)+n(y)∣πf(x1) ;*
•
V(f)⊇Sy* if, and only if, x2−xt(y)+n(y)∣f(x).*
Proof.
We begin with the first statement.
•
If x−y∣f(x) then f(y)=0 by Corollary 6.2. Now let us prove the converse implication. Let J∈SDH be an imaginary unit such that y∈ΩJ and let {1,J,e1J,Je1J,e2J,Je2J,e3J,Je3J} be an associated splitting basis of DH.
By Lemma 3.8, there exist holomorphic functions f0,f1,f2,f3:ΩJ⟶CJ such that
[TABLE]
where e0J:=1. If f(y)=0 then, for each l∈{0,1,2,3}, it holds fl(y)=0 and there exists a holomorphic function gl:ΩJ⟶CJ such that
[TABLE]
for all z∈ΩJ.
Claim 1.There exists g∈SR(Ω) with g∣ΩJ=∑l=04glelJ.
Claim 2.f(x)=(x−y)⋅g(x)* in Ω.
*We postpone the proofs of our claims until later and we complete the proof of the first statement by computing
[TABLE]
We now suppose y∈Ω∖R and prove the third and second statement, reversing the order of presentation for the sake of clarity.
•
V(f)⊇Sy is equivalent to f(y)=0=fs′(y) by Theorem 5.2. By what we have already proven, f(y)=0 is equivalent to the existence of g∈SR(Ω) with f(x)=(x−y)⋅g(x). By formula (11), the last equality implies
[TABLE]
Thus, f(y)=0=fs′(y) is equivalent to the existence of g,h∈SR(Ω) with f(x)=(x−y)⋅g(x) and with g(x)=(x−yc)⋅h(x). This is, in turn, equivalent to asking for
[TABLE]
to divide f(x).
•
V(f)⊇Ty1 is equivalent to f(y)=0,fs′(y)∈ϵH∗ by Theorem 5.2. This is equivalent to the existence of g∈SR(Ω) with f(x)=(x−y)⋅g(x) and with πg(y1c)=0 (by Remark 4.7). This is, in turn, equivalent to asking for x−y to divide f(x) and for
[TABLE]
to divide πf(x1).
Proof of claim 1. Setting
[TABLE]
for all α+iβ∈D defines a stem function G=G1+ιG2:D→DH. G is holomorphic by direct inspection, whence it induces a slice regular function g=I(G)∈SR(Ω). Now, for all z=α+Jβ∈ΩJ,
[TABLE]
as desired.
Proof of claim 2. We wish to prove that h(x):=(x−y)⋅g(x) coincides with f(x) in Ω. As explained in Section 3, it suffices to prove that they coincide in ΩJ. Now, for all z∈ΩJ, Theorem 6.1 guarantees that h(z)=(z−y)g(z). Thus,
[TABLE]
as desired.
∎
Corollary 7.7**.**
Let f∈SR(Ω) with Ω=ΩD and let y∈Ω∖R. If
[TABLE]
then f has a unique zero w∈Sy and x−w∣f(x). If
[TABLE]
then either x−w∣f(x) for all w∈Sy or the subset of those w∈Sy such that x−w∣f(x) is a tangent plane Th1⊆Sy. Finally, if
[TABLE]
then x−w∣f(x) for all w∈Sy.
Proof.
Let us prove our first statement. According to Theorem 7.6, V(N(f)) includes Sy but V(πf) does not include Sy∩H. By Theorem 5.6, f has a unique zero w∈Sy. We easily conclude that x−w∣f(x) by a further application of Theorem 7.6.
The second statement follows directly from Theorem 7.6.
As for the third statement, it follows from the fact that x−w∣x2−xt(y)+n(y) for all w∈Sy.
∎
The previous corollary is false when y∈Ω∩R, as proven by the next examples.
Example 7.8**.**
Let us fix y=y1∈R and v∈Im(H). Define f∈SR(QDH) by the formula
[TABLE]
Then πf(x1)=x1−y1 and N(f)(x)=x2−2xy1+y12=x2−xt(y)+n(y). Nevertheless, in Sy={y} it holds f(y)=0 and x−y∣f(x).
Example 7.9**.**
Let us fix y=y1∈R and define f∈SR(QDH) by the formula
[TABLE]
Then πf(x1)=x12−x1t(y)+n(y) and x2−xt(y)+n(y)∣f(x). Nevertheless, in Sy={y} it holds f(y)=0 and x−y∣f(x).
The difference between the non real case and the real case is explained by the next remark, which involves the concept characteristic polynomial of a dual quaternion h∈DH: a polynomial having minimal degree among all monic real polynomial vanishing at h.
Remark 7.10**.**
Let us fix h∈DH. The norm of t−h, namely
[TABLE]
is a real polynomial if, and only if, h∈(R+ϵIm(H))∪QDH. Moreover:
•
If h∈QDH∖R, then Δh(t) is the characteristic polynomial of each element of its zero set Sh.
•
If h∈R+ϵIm(H), then Δh(t)=(t−h1)2 is the characteristic polynomial of each element of its zero set h1+ϵH except h1, whose characteristic polynomial is t−h1.
Finally, for each monic quadratic real polynomial M(t): if M(t) is irreducible in R[t], then it equals Δh(t) for some h∈QDH∖R; if M(t)=(t−h1)2 with h1∈R, then M(t)=Δh1(t); if M(t) has two real roots, then it has no other roots in DH and it does not coincide with Δh(t) for any h∈DH.
8 Applications to the study of motion polynomials
This section explains the meaning of possible factorizations of the motion polynomial P(t) in formula (14) and it studies existence and uniqueness of such factorizations. This is done combining material from [11] with applications of our new results about the zeros of slice functions over dual quaternions.
Throughout this section, we only consider monic polynomials. There is no loss of generality in doing so, because multiplying the motion polynomial P(t) in formula (14) by the inverse of its leading coefficient will only result in a change in the coordinate frame.
8.1 Meaning of factorization of motion polynomials
Let us consider a monic linear polynomial
[TABLE]
with h∈DH. P(t) is a motion polynomial if, and only if, N(P)(t)=Δh(t) belongs to R[t]. This is, in turn, equivalent to h∈(R+ϵIm(H))∪QDH. We distinguish two cases for the transformations
[TABLE]
Translations.
If h∈R+ϵIm(H), then the previous transformation is a translation with translation vector 2h2(t−h1)−1. The direction of the vector does not depend on t, although its length does: the trajectories are straight lines parallel to h2.
Rotations.
If h∈QDH∖R, then the previous transformation is a rotation. The rotation axis, which has Plücker coordinates (∣Im(h1)∣Im(h1),∣Im(h1)∣h2), does not depend on t. The rotation angle θ is determined by the equality cos(2θ)=∣t−h1∣t−Re(h1). Thus, the trajectories are circles whose centers lie on the fixed axis.
Motions of former type have been identified in [11] as those of linkages consisting of a single prismatic joint. Similarly, motions of the latter type are associated to linkages consisting of a single revolute joint. Let us give an explicit example of the latter type.
Example 8.1**.**
For all t∈R, consider the monic linear polynomial
[TABLE]
It is a motion polynomial because N(P)(t)=t2+1∈R[t]. For each t0∈R, we consider the proper rigid transformation
[TABLE]
which is a rotation around the axis with Plücker coordinates (i,−j) (i.e., the axis with direction i through the point −i∧j=−k). For instance, the trajectory of 1+ϵk, corresponding to the point (0,0,1) in R3, is the rational curve
[TABLE]
It is a parametrization of the circle of radius 2 centered at (0,0,−1) in the plane of points (x1,x2,x3) with x1=0, except for the point (0,0,1), which is the limit of the curve as t→±∞.
Now let us consider a motion polynomial P(t)∈DH[t] of degree n>0. Suppose P(t) admits a factorization
[TABLE]
Then the transformation in formula (13) is the composition of a number s∈{1,…,n} of pure rotations and of n−s pure translations, completely determined by the dual quaternions h(1),…,h(n) and by the parameter t0. The resulting trajectories have been identified in [11] as the motions of linkages consisting of s revolute joints and n−s prismatic joints. This motivated the search for sufficient conditions on a motion polynomial P(t)∈DH[t] that guarantee the existence of a factorization of the form (19).
8.2 Sufficient conditions for the existence of factorizations
The work [11] proved that the following property is a sufficient condition for the existence of a factorization of the form (19), with h(1),…,h(n)∈QDH∖R (corresponding to a linkage with n revolute joints).
Definition 8.2**.**
A motion polynomial P(t)∈DH[t] is called generic if every real polynomial that divides primal(P)(t) is constant.
The key ingredient in the proof of this sufficient condition was [11, Lemma 3], which we can restate and prove as follows (again, with a different convention about the side of the coefficients of polynomials in DH[t]).
Lemma 8.3**.**
Let P(t)∈DH[t] and let y∈(R+ϵIm(H))∪QDH. If
[TABLE]
then P(t) and Δy(t) have a unique common root h and t−h∣P(t). Moreover, h∈(R+ϵIm(H))∪QDH.
Proof.
If y∈QDH∖R, then the thesis immediately follows from Corollary 7.7. If, instead, y∈R+ϵIm(H), then the results of Section 7 do not generally apply. Therefore, we provide a direct proof of our thesis.
Our first step is proving that P(t) and Δy(t) have a unique common zero h. Since Δy(t)=Δy1(t)=(t−y1)2 divides N(P)(t), Corollary 6.2 tells us that N(P)(t) vanishes at the real point y1. By Theorem 5.6, so does primal(P)(t). Thus,
[TABLE]
We have P(t)=(t−y1)⋅Q(t)+r for some other polynomial Q(t)∈DH[t] and some r∈DH. In this situation, primal(P)(t)=(t−y1)⋅primal(Q)(t)+r1 and conditions (20) imply that r1=0, while primal(Q)(y1)=q∈H∗. As a consequence, P(t)=(t−y1)⋅Q(t)+ϵr2 and, for all v∈H,
[TABLE]
where the second equality follows by applying Remark 3.10 to t−y1∈R[t]⊂DR[t] and the third equality follows from Remark 4.10. As a consequence, the unique zero of P(t) in the zero set y1+ϵH of (t−y1)2=Δy(t) is the point
[TABLE]
As a second step, we prove that h belongs to y1+ϵIm(H), i.e., that t(r2q−1)=0. By construction,
[TABLE]
The hypothesis (t−y1)2∣N(P)(t) implies that t−y1∣R(t). We conclude that
[TABLE]
where the second and third equalities follow by applying Remark 3.10 at y1∈R⊂DR. As a consequence, t(r2q−1)=n(q)t(r2qc)=0, as desired.
Our third step is proving that t−h∣P(t). To do so, let us start again from the equality P(t)=(t−y1)⋅Q(t)+ϵr2. By Remark 3.10,
[TABLE]
whence there exists Q(t)∈DH[t] such that Q(t)=(t−y1)⋅Q(t)+q. Taking into account that (t−y1)2=Δh(t), we conclude that
[TABLE]
Since t−h∣Δh(t) by definition, our thesis t−h∣P(t) immediately follows.
∎
The same work [11] hinted it was possible to generalize the sufficient condition for the existence of a factorization of the form (19) to allow both revolute and prismatic joints. This is done in the next definition and in the subsequent result, which uses both Lemma 8.3 and the Quaternionic Fundamental Theorem of Algebra (see, e.g., [3, Theorem 3.18]).
Definition 8.4**.**
A motion polynomial P(t)∈DH[t] is called Δ-free if
[TABLE]
for all y∈(R+ϵIm(H))∪QDH; or, equivalently, for all y=y1∈H.
Theorem 8.5**.**
Let P(t)∈DH[t] be a Δ-free monic motion polynomial of degree n>0. Then P(t) admits a factorization P(t)=(t−h(1))⋅…⋅(t−h(n)) with h(1),…,h(n)∈(R+ϵIm(H))∪QDH.
Proof.
We proceed by induction on the degree n. For the case n=1, we already observed by direct computation that P(t)=t−h(1) with h(1)∈(R+ϵIm(H))∪QDH. Now let us suppose our thesis proven for all degrees 1,…,n−1 and let us establish it for degree n.
According to the Quaternionic Fundamental Theorem of Algebra, the quaternionic polynomial primal(P) admits a zero y1∈H. Thus,
[TABLE]
As a consequence, Δy1(t) divides the norm of primal(P)(t), which coincides with N(P)(t) according to Remark 4.12. By the definition of Δ-free motion polynomial, Δy1(t)∣primal(P)(t). According to Lemma 8.3,
[TABLE]
for some polynomial Q(t)∈DH[t] of degree n−1 and some point h(1)∈(R+ϵIm(H))∪QDH that is a zero of Δy1(t). By direct computation,
[TABLE]
whence Q(t) is still a Δ-free motion polynomial. Our inductive hypothesis guarantees the existence of a factorization Q(t)=(t−h(2))⋅…⋅(t−h(n)) with h(2),…,h(n)∈(R+ϵIm(H))∪QDH. It immediately follows that P(t) admits a factorization of the form (19), as desired.
∎
Remark 8.6**.**
A Δ-free motion polynomial P(t) of degree n has as many different factorizations as the quaternionic polynomial primal(P)(t) does. This follows by direct inspection in the previous proof. In particular, P(t) has at most n! distinct factorizations.
Techniques to classify all possible factorizations of the quaternionic polynomial primal(P)(t) are described in [3, §3.5]. We present here two significant examples.
Example 8.7**.**
Let P(t) be the following Δ-free motion polynomial:
[TABLE]
Its primal part
[TABLE]
has exactly two roots, namely −2j and 53i−54j. As a consequence of [3, Theorem 3.24], primal(P) admits exactly two factorizations, the second one being
[TABLE]
If we repeat the proof of Theorem 8.5 starting with this second factorization of primal(P)(t), we find a second factorization of P(t), namely
[TABLE]
Example 8.8**.**
The Δ-free motion polynomial
[TABLE]
admits a unique factorization. Indeed, its primal part
[TABLE]
only vanishes at {i} and has a unique factorization by [3, Proposition 3.23].
8.3 Factorization in general
In this final subsection, we treat the problem of factorization of motion polynomials in general. In addition to the results of the previous subsection, we will make use of the next lemma.
Lemma 8.9**.**
Let P(t)∈DH[t]. Suppose y∈(R+ϵIm(H))∪QDH is such that Δy(t)∣primal(P)(t). The set of common zeros of Δy(t) and P(t), which is the set of zeros h of Δy(t) such that t−h∣P(t), may be:
the whole zero set of Δy(t);
2. 2.
a tangent plane Tw1⊂Sy, provided y∈QDH∖R;
3. 3.
the empty set.
Proof.
If y∈QDH∖R, then the thesis immediately follows from Corollary 7.7. Let us therefore suppose that y∈R+ϵIm(H), whence the zero set of Δy(t)=(t−y1)2 is y1+ϵH.
Let us divide P(t) by Δy(t): P(t)=Δy(t)Q(t)+R(t) for some Q(t),R(t)∈DH[t] with degR(t)≤1. If R(t) is constant then we are either in case 1. or in case 3. Suppose, instead, R(t)=ta+b with a=0.
The hypothesis Δy(t)∣primal(P)(t) implies that primal(R)(t)=0. Thus, R(t)=tϵa2+ϵb2 with a2∈H∗ and
[TABLE]
has a unique quaternionic zero, namely −b2a2−1. If this zero coincides with y1, then the zero set of R(t)=(t−y1)ϵa2 is y1+ϵH. Moreover, the equality R(t)=(t−h)ϵa2 holds for all h∈y1+ϵH, whence t−h∣P(t), and we are in case 1. If, instead, −b2a2−1 does not coincide with y1 then, for every h∈y1+ϵH, the polynomial R(t) does not vanish at h and t−h∣R(t). Thus, we are in case 3.
∎
Example 8.10**.**
The motion polynomials
[TABLE]
have
[TABLE]
for all ℓ∈{1,2,3}. The zero sets of P1(t),P2(t),P3(t) are, respectively, S,Ti,∅.
We are now ready for a general description of motion polynomials. To this end, the next definition will be useful
Definition 8.11**.**
A motion polynomial P(t)∈DH[t] is Δ-covered if, for each zero y1∈H of primal(P)(t), it holds
[TABLE]
Proposition 8.12**.**
Let P(t)∈DH[t] be a monic motion polynomial of degree n>0. Then
[TABLE]
where Q(t) is a Δ-free motion polynomial of degree q≤n and R(t) is a Δ-covered motion polynomial. The polynomial Q(t) admits at least 1 and at most q! factorizations of the form
[TABLE]
each corresponding to exactly one factorization primal(Q)(t)=(t−h1(1))⋅(t−h1(2))⋅…⋅(t−h1(q)) of its primal part.
The polynomial R(t) admits either infinitely many or no factorization of the form
[TABLE]
If it does, then Δh1(ℓ)∣primal(R)(t) for all ℓ∈{q+1,q+2,…,n}.
Proof.
Consider the quaternionic polynomial primal(P)(t): thanks to the techniques described in [3, §3.5], it is possible to find a factorization
for some h(1)∈(R+ϵIm(H))∪QDH with h1(1)=p1(1) and for some P(1)(t)∈DH[t] with
[TABLE]
After q−1 further applications of Lemma 8.3, we find that
[TABLE]
with h(1),…,h(q)∈(R+ϵIm(H))∪QDH and with
[TABLE]
If we set R(t):=P(q)(t), then the first and second statements are proven.
To prove the third statement, we observe that if R(t) admits a factorization of the form
[TABLE]
then it admits infinitely many because of Lemma 8.9. Furthermore, the previous equality implies that
[TABLE]
As a consequence, the zero set of primal(R)(t) intersects Sh1(ℓ)∩H for all ℓ∈{q+1,q+2,…,n}. Since R(t) is Δ-covered by construction, it follows that Δh1(ℓ)∣primal(R)(t) for all ℓ∈{q+1,q+2,…,n}.
∎
Example 8.13**.**
The motion polynomial
[TABLE]
has primal(P)(t)=t3+t2(i+j)=(t+i+j)⋅t2 and N(P)(t)=(t2+2)⋅t4. It factorizes as
[TABLE]
Example 8.14**.**
The motion polynomial
[TABLE]
has primal(P)(t)=(t2+1)2 and N(P)(t)=(t2+1)4. It holds
[TABLE]
and R(t) never vanishes in (R+ϵIm(H))∪QDH.
Bibliography15
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] F. Colombo, I. Sabadini, and D. C. Struppa. Slice monogenic functions. Israel J. Math. , 171:385–403, 2009.
2[2] K. Daniilidis. Hand-eye calibration using dual quaternions. Internat. J. Robotics Res. , 18(3):286–298, March 1999.
3[3] G. Gentili, C. Stoppato, and D. C. Struppa. Regular functions of a quaternionic variable . Springer Monographs in Mathematics. Springer, Heidelberg, 2013.
4[4] G. Gentili and D. C. Struppa. A new approach to Cullen-regular functions of a quaternionic variable. C. R. Math. Acad. Sci. Paris , 342(10):741–744, 2006.
5[5] G. Gentili and D. C. Struppa. A new theory of regular functions of a quaternionic variable. Adv. Math. , 216(1):279–301, 2007.
6[6] G. Gentili and D. C. Struppa. Regular functions on a Clifford algebra. Complex Var. Elliptic Equ. , 53(5):475–483, 2008.
7[7] G. Gentili and D. C. Struppa. Regular functions on the space of Cayley numbers. Rocky Mountain J. Math. , 40(1):225–241, 2010.
8[8] R. Ghiloni and A. Perotti. Slice regular functions on real alternative algebras. Adv. Math. , 226(2):1662–1691, 2011.