Preprint
Article

This version is not peer-reviewed.

Differential Geometric Analysis of Curves and Surfaces Generated by the Gielis Superformula

Submitted:

14 April 2026

Posted:

14 April 2026

You are already at the latest version

Abstract
The Gielis superformula is a powerful parametric tool that generates an infinite variety of natural and organic curves and surfaces through a compact set of parameters. However, classical differential geometry has lacked a unified framework for analyzing their curvature, torsion, and intrinsic geometric properties. This study addresses this gap by developing a novel superelliptic geometric framework that integrates the superformula 6with the differential geometry of curves and surfaces. We define the superelliptic inner and cross products, the star derivative, and the superelliptic Frenet frame to extend Euclidean and Riemannian interpretations of curvature and torsion to a more flexible parametric structure. The framework provides a uniform geometric characterization of all Gielis curves and surfaces, independent of their classical parametric expressions; even singular cases are regularized so that their curvature and torsion reduce exactly to those of a circle. This unifies the entire family under a common, robust foundation while preserving orthonormality and differentiability. This superelliptic approach offers a consistent and computationally 14tractable model that bridges mathematical abstraction with real-world morphology, with 15the superformula serving as a representative example of the framework’s broad generality for diverse geometric structures.
Keywords: 
;  ;  ;  

1. Introduction

One of the most traditional and basic areas of differential geometry is the theory of curves, which focuses on the local and global characteristics of smooth curves in Euclidean spaces [1,9]. The formal definition of a curve is a smooth mapping
γ : I R R 3 , γ ( t ) = ( x ( t ) , y ( t ) , z ( t ) ) ,
where t is either a differentiable reparameterization of the arc-length or a real parameter that represents it. The derivatives of such a curve, which provide details about its velocity, curvature, and torsion, are used to examine its geometric behavior. The Frenet-Serret frame is an orthonormal moving trihedron made up of the tangent vector T ( t ) , the normal vector N ( t ) , and the binormal vector B ( t ) . The definition of this frame is as follows:
T = d γ / d t d γ / d t , N = d T / d t d T / d t , B = T × N .
The Frenet-Serret equations control how this frame changes along the curve:
d d t T N B = 0 κ 0 κ 0 τ 0 τ 0 T N B ,
where κ ( t ) and τ ( t ) represent the curve’s curvature and torsion, respectively. Up to rigid motion, these numbers fully characterize the curve [7].
Simple yet fundamental examples of curves, like circles and circular helices, illustrate the geometric meaning of these invariants. A circle has zero torsion τ = 0 and constant curvature κ = 1 / r , but a helix has both constant curvature and constant torsion and maintains a constant ratio τ / κ . These canonical forms serve as reference models for understanding more complex curves in physics, biology, and computer-aided design [5].
Curve theory’s current advancement uses parametric and computational modeling to expand on traditional concepts. Specifically, the introduction of the Superformula, which is a generalization of supercircle and superellipses [3,14], offers a potent parametric generalization that can describe a variety of closed and open curves seen in nature. The names superformula and supershapes originate from the names superellipses and superquadrics. The name superformula was changed by mathematicians to the Gielis Formula [16,18,19], and supershapes to Gielis’ curves and surfaces [15,17]. It has the following definition in polar coordinates:
r ( t ) = [ | cos m t 4 a | n 2 + | sin m t 4 b | n 3 ] 1 n 1
The superformula introduces the parameters m, a, b, n 1 , n 2 , and n 3 , which govern the symmetry, curvature variation, and geometric complexity of the curve, where a, b, and n 1 are assumed to be nonzero. It generalizes Lamé’s superellipse [8] and Piet Hein’s superellipse [6]. This formula can produce smooth transitions between circular, polygonal, star-like, and biomorphic shapes when the right parameters are used [3,13,14,15]. Gielis transformations denote the superformula applied as a transformation on other functions. The name "Universal natural shapes" was introduced by Leopold Verstraelen to denote that the formula can be used to describe/model shapes and forms at any level [15,20]. The way back to Cartesian via Chebyshev polynomials is first used in [34]. This includes a wide range of possibilities, including trigonometric and exponential functions, as well as polynomials. The cosine and sine functions in the superformula can also be other functions, like cosh/sinh or even many other functions. One could set math requirements for this (closed curves, going through the origin,….), see [21]. It has been widely used in computer-aided geometric design, biological modeling, and form optimization [22,23,24,25,26,27,28,29,30,31,32,33].
Recent advances by Özdemir and Parlak [10,11] extend the Superformula by integrating superelliptic and quaternionic algebraic structures, thereby enabling the unified treatment of superelliptic rotational forms and spatial transformations within a single analytical framework. These methods enable the study of curvature-driven morphologies and symmetric deformations in higher dimensions by fusing traditional Frenet-Serret theory with contemporary parametric geometry.
In this study, the differential geometric properties of the curves generated by the superformula are investigated. Using the Frenet-Serret formalism, we aim to determine the curvatures, torsions, and frame evolutions of these parametrically defined curves both analytically and numerically. Furthermore, the theoretical framework is supported by illustrative examples and visualizations. This approach enables a clearer, visually intuitive understanding of the geometric behavior and kinematic evolution of the curves. In this study, the intrinsic geometry of superformula-based curves is elucidated, and the interaction between geometric invariants and parameter management is emphasized. As a result, integrating parametric models such as the superformula with classical curve theory creates new opportunities for the precise differential geometric definition and construction of complex, nature-inspired structures.
Although Gielis curves are capable of modeling a wide variety of natural and synthetic shapes within a large parametric family, they exhibit singularities and sharp corner behavior, particularly for certain specific parameter choices, such as asteroids. This complicates the definition of differential geometric quantities such as curvature, torsion, and differentiability at every point.
In this study, we demonstrate that all Gielis shapes can be reduced to a common geometric class by using a new Frenet-frame structure defined for Gielis curves. This structure provides a transformation that allows each Gielis curve to be reinterpreted within a framework equivalent to a circle with constant curvature and zero torsion. Thus, the Gielis family is given a globally regular structure despite its local singularities, providing a holistic basis for differential-geometric investigations.

2. Basic Concepts and Notions

The most essential part of this paper is the definition and some properties of the Gielis formula. Gabriel Lamé expanded the formulas for circles and ellipses to encompass squares, rectangles, parallelograms, astroids, and all conic sections. Subsequently, he formulated a singular equation expressed as follows: where a , b , n are positive constants.
x a n + y b n = 1 .
As a result, the notion of Lamé curves is incorporated into geometry [36,37]. This study [3] demonstrates the influence of parameter values in the Lamé curve equation on the geometric forms of ellipses and rectangles. An expression for Lamé’s extended equation in polar coordinates is described below:
r ( θ ) = | cos θ a | + | sin θ b | n 1 n .
In the 1990s, superellipses were generalized for any symmetry [3].
r ( t ) = [ | cos m t 4 a | n 2 + | sin m t 4 b | n 3 ] 1 n 1 .
For the Lamé equation, when t [ π , π ] , this formulation is referred to as the superformula. The parameters m, a, and b are real-valued, while n 1 , n 2 , and n 3 are typically taken as positive real constants, although negative values have also been considered in the literature [38]. By varying the parameters in a novel, straightforward, and unifying equation that transcends the limitations of Euclidean and Pythagorean measurements. In the literature, "superformula" and "supershapes" are also referred to as "Superellipse," "Gielis formula," and "Gielis curve".
We introduced the following three-dimensional form, which is symmetric and positive definite bilinear, which was utilized in our research: for any two vectors u = ( u 1 , u 2 , u 3 ) , v = ( v 1 , v 2 , v 3 ) R 3
, S E : R 3 × R 3 R ; u , v S E = 1 r ( t ) 2 σ 1 u 1 v 1 + 1 r ( t ) 2 σ 2 u 2 v 2 + 1 r ( t ) 2 σ 3 u 3 v 3 .
where r ( t ) = [ | cos m t 4 a | n 2 + | sin m t 4 b | n 3 ] 1 n 1 and σ i R + , i = { 1 , 2 , 3 } . The vectors u , v will be defined as superelliptic vectors in this context. A superelliptic inner product describes this type of multiplication. The space R 3 equipped with the superelliptic inner product is referred to as the superelliptic 3-space, denoted as R S E 3 = ( R 3 , , S E ) . The superelliptic inner product , S E can also be written as u , v S E = u t G v , where the associated matrix is
G = 1 r ( t ) 2 σ 1 0 0 0 1 r ( t ) 2 σ 2 0 0 0 1 r ( t ) 2 σ 3 .
The superelliptic inner product’s length is supplied by the superelliptic vector u R S E 3
u S E = u , u S E .
A pair of superelliptic vectors u and v are considered superelliptic orthogonal in R S E 3 if and only if u , v S E = 0 . If the superelliptic norm of a superelliptic vector equals 1, the vector is termed a superelliptic orthonormal vector. We will refer to the conventional basis of Euclidean 3-space as { e 1 , e 2 , e 3 } . Based on this, we can represent the base vectors of R S E 3 as { e 1 ´ , e 1 ´ , e 1 ´ } , and they meet the condition
e 1 ´ = r ( t ) σ 1 e 1 = ( r ( t ) σ 1 , 0 , 0 ) , e 2 ´ = r ( t ) σ 2 e 2 = ( 0 , r ( t ) σ 2 , 0 ) , e 3 ´ = r ( t ) σ 3 e 3 = ( 0 , 0 , r ( t ) σ 3 ) ,
and this form can be used to express any vector u = ( u 1 , u 2 , u 3 ) R S E 3 .
u = u 1 e 1 ´ + u 2 e 2 ´ + u 3 e 3 ´ = u 1 r ( t ) σ 1 e 1 + u 2 r ( t ) σ 2 e 2 + u 3 r ( t ) σ 3 e 3 .
Given two non-zero superelliptic vectors u and v , find the cosine function of the angle between them.
cos ϕ = u , v S E u S E v S E .
Here ϕ agrees with the angular parametric equations of a superellipse or a superellipsoid. In R S E 3 , the superelliptic cross product of two vectors { u , v } is computed as
u × S E v = Λ * r ( t ) 2 e 1 σ 1 r ( t ) 2 e 2 σ 2 r ( t ) 2 e 3 σ 3 u 1 u 2 u 3 v 1 v 2 v 3 ,
where Λ * = 1 r ( t ) 3 σ 1 σ 2 σ 3 , σ 1 , σ 2 , σ 3 R + (See [11], for details).
Before developing the superelliptic geometric framework, let us define the following function that will be used in establishing this structure.
Definition 1.
Consider a vector function F = ( f 1 , f 2 , f 3 , . . . , f n ) such that the function F : I R R n is a smooth function. Then, we can define the function ϝ as
ϝ : R n R S E n F ( t ) ϝ ( F ( t ) ) = ( f 1 ( t ) r ( t ) σ 1 , f 2 ( t ) r ( t ) σ 2 , f 3 ( t ) , . . . , f n ( t ) r ( t ) σ n ) ,
where r ( t ) is the Gielis’ superformula and σ i , i { 1 , 2 , . . . , n } , positive real number. This function is referred to as the Gielis radial scaling transformation.
The Gielis radial scaling transformation satisfies the following conditions:
1.
Since each component function f i ( t ) r ( t ) σ i , i { 1 , 2 , . . . , n } , is of class C , the mapping ϝ is also a smooth.
2.
Since r ( t ) 0 , the function ϝ is one-to-one and its inverse can be expressed as:
f i ( t ) = ϝ i ( t ) · σ i r ( t ) , i = 1 , , n .
3.
The Gielis transformation can be written in matrix form:
ϝ ( F ( t ) ) = D ( t ) · F ( t ) ,
where D ( t ) is a time-dependent diagonal matrix:
D ( t ) = diag r ( t ) σ 1 , r ( t ) σ 2 , , r ( t ) σ n .
Under this transformation, the inner product changes to:
ϝ ( F ( t ) ) , ϝ ( H ( t ) ) S E = i = 1 n f i ( t ) h i ( t ) = F ( t ) , H ( t ) .
which shows that ϝ preserves the inner product, i.e., it preserves both angles and lengths. That is, in R n S E , the Gielis radial scaling transformation ϝ is an isometry.
Therefore, the Gielis radial scaling transformation ϝ defines an isometric diffeomorphism on R n S E .

3. Superelliptic Curve

In this section, we examine the properties of superelliptic curves from the perspective of differential geometry.
Definition 2.
Let I be an interval in space R and α : I R R n ; t α ( t ) = ( α 1 ( t ) , α 2 ( t ) , , α n ( t ) ) be a curve. In R n S E , the superelliptic curve defined via the following Gielis radial scaling transformation ϝ as follows:
ϝ : R n R S E n α ϝ ( α ( t ) ) = α S E = ( α 1 ( t ) r ( t ) σ 1 , α 2 ( t ) r ( t ) σ 2 , α n ( t ) r ( t ) σ n )
where r ( t ) is the Gielis’ super formula and σ i , i { 1 , 2 , . . . , n } positive real number.
Example 1.
In R 2 , the unit circle can be expressed as α ( t ) = ( cos t , sin t ) , then the corresponding superelliptic circle (superellipse) can be obtained as ϝ ( α ( t ) ) = α S E = ( cos t r ( t ) σ 1 , sin t r ( t ) σ 2 ) . Similarly, a right helix can be given as α ( t ) = ( cos t 2 , sin t 2 , t 2 ) , then the corresponding superelliptic helix can be obtained as ϝ ( α ( t ) ) = α S E = ( cos t 2 r ( t ) σ 1 , sin t 2 r ( t ) σ 2 , t r ( t ) 2 σ 3 ) . Figure 1 and Figure 2 illustrate the superellipses and superelliptic helices for some values of the Gielis formula r ( t ) .
In order to remain compatible with the metric · , · S E , the direct use of the classical derivative is no longer sufficient, since it does not incorporate the scaling effects induced by the superelliptic structure. For this reason, we introduce a modified differentiation operator that reflects the intrinsic geometry of the space R S E n .
Via the transformation ϝ , the classical derivative of a curve is mapped into the superelliptic space, giving rise to a new derivative operator, which we call the star derivative (or *-derivative). This operator simultaneously encodes the directional behavior of the curve and the metric-dependent scaling factor.
Definition 3.
The star derivative operator, denoted by D * , is defined as the composition of the radial scaling ϝ with the classical differential operator D = d d t :
D * : = ϝ D .
Definition 4.
The tangent vector at point t on the curve α : I R R n is denoted by α ( t ) = ( α 1 ( t ) , α 2 ( t ) , , α n ( t ) ) . In this case, under ϝ transformation, the tangent vector T S E of the superelliptic curve α S E is determined as follows:
α S E * = ϝ ( α ) = T S E = ( d α 1 ( t ) d t r ( t ) σ 1 , d α 2 ( t ) d t r ( t ) σ 2 , d α n ( t ) d t r ( t ) σ n ) .
Thus, the k-th order star derivative of the curve α is given by the pushforward of the k-th order classical derivative under the map ϝ:
α S E ( k ) * ( t ) : = ϝ α ( k ) ( t ) = r ( t ) α 1 ( k ) ( t ) σ 1 , , r ( t ) α n ( k ) ( t ) σ n
The star derivative is a direct transfer of the classical derivative to superelliptic space. This allows quantities such as the velocity and acceleration of the curve to be expressed under the new metric. Since the transformation ϝ preserves angles, the vectors obtained with the star derivative also preserve angular relationships. Since derivatives are taken component-wise, all the properties of the classical derivative and the chain rule apply to the star derivative.
The star derivative represents the tangent vector field measured with respect to the Gielis-type metric rather than the standard Euclidean one. It accounts for the local stretching and compression of the ambient space.
Since ϝ is a radial transformation, the star derivative preserves the angular orientation of the classical derivative while rescaling its magnitude. It describes the curve’s evolution within a deformed geometry where distances are governed by the superelliptic scaling factor.
This construction ensures that local geometric properties, such as curvature and torsion, are consistently defined within the superelliptic framework. The resulting Frenet-Serret frame in R S E n corresponds to a scaled version of the Euclidean frame, maintaining compatibility with the induced inner product of the space.
In Figure 3, The transformation of a fundamental Euclidean unit circle into a superelliptic geometry, and the subsequent modification of its associated tangent vector fields, is depicted. This visual representation elucidates the operational impact of the Gielis radial scaling on the curve’s differential characteristics. Consequently, Figure 3 demonstrates that the star derivative is the natural velocity vector of the curve when the ambient space itself is deformed. This approach provides a more robust framework for analyzing organic morphologies where growth rates and surface tangents are dictated by non-Euclidean radial constraints.
Definition 5.
Let α S E : I R R S E n be a superelliptic curve at each point t I . Then α S E is called a unit-speed superelliptic curve if and only if it satisfies α S E * ( t ) S E = 1 .
Example 2.
Consider that the superelliptic curve α S E = ( cos t r ( t ) σ 1 , sin t r ( t ) σ 2 ) in R S E 2 then we compute that it is a superelliptic tangent vector as the form α S E * = ( sin t r ( t ) σ 1 , cos t r ( t ) σ 2 ) and since α S E * ( t ) S E = 1 . Therefore, the curve is of unit-speed in the superelliptic sense.
Proposition 1.
Let α : I R R n be a unit-speed curve, that is, α ( t ) = 1 . Then, at each point t I , the related superelliptic curve ϝ ( α ( t ) ) = α S E ( t ) also satisfies α S E * ( t ) S E = 1 , which confirms that the corresponding superelliptic curve is of unit-speed.
Definition 6.
A superelliptic curve α S E : I R R S E n is called a regular superelliptic curve if and only if, at each point t I , it satisfies α S E * ( t ) S E 0 .
Proposition 2.
Let α : I R R n be a regular curve such that α ( t ) 0 . Then at each point t I , the related superelliptic curve ϝ ( α ( t ) ) = α S E ( t ) is also a superelliptic regular curve.

4. An Orthonormal Frame of a Superelliptic Curve in R S E 2

In this section, we explain how to construct an orthonormal frame at any point on the regular superelliptic curve α S E . For a unit-speed curve α S E = ( α 1 ( t ) r ( t ) σ 1 , α 2 ( t ) r ( t ) σ 2 ) , we know that α * ( t ) = T S E ( t ) such that T S E , T S E S E = 1 . This gives T S E * , T S E S E = 0 . Thus, we can write T S E * = T S E * , N S E S E N S E = κ S E N S E . Finally, we compute N S E * = T S E * , N S E S E T S E = κ S E T S E . Accordingly, we can write the following matrix equations;
T S E * N S E * = 0 κ S E κ S E 0 T S E N S E .
Theorem 1.
Assume α is a unit-speed curve with the Frenet frame { T , N , κ } . Then, for the unit-speed superelliptic curve ϝ ( α ) = α S E , we have α S E * S E = 1 . Thus, the superelliptic Frenet frame { T S E , N S E , κ S E } along the superelliptic curve ϝ ( α ) = α S E is obtained as { ϝ ( T ) = T S E , ϝ ( N ) = N S E , κ S E } .
Proof. 
Let α be a unit-speed curve with the Frenet frame apparatus { T , N , κ } . Then, we compute that ϝ ( α ) = α S E is a unit superelliptic curve, and we obtain:
ϝ ( T ( t ) ) = ϝ ( α ( t ) ) = α S E * ( t ) = T S E ( t ) ,
and since κ = α = ϝ ( α ) S E = κ S E , we reach
ϝ ( N ) = ϝ ( T κ ) = T S E * κ S E = N S E ( t ) .
Theorem 2.
Let { T , N } be the orthonormal frame of any regular curve α in R 2 such that ν = α and κ , τ are the curvature function and the torsion function, respectively. Then the related superelliptic regular curve α S E with the superelliptic orthonormal frame { ϝ ( T ) = T S E , ϝ ( N ) = N S E } such that ν S E = ϝ ( α ) S E = α S E * S E has the following superelliptic frame equation
T S E * N S E * = 0 ν S E κ S E ν S E κ S E 0 T S E N S E .
Proof. 
First, since we have T T , N N , and N T , it is clear that under the ϝ transformation, it will be ϝ ( T ) ϝ ( T ) , ϝ ( N ) ϝ ( N ) , and ϝ ( T ) ϝ ( N ) . Since ϝ ( T ) ϝ ( N ) , we have ϝ ( T ) = a ϝ ( N ) and ϝ ( N ) = b ϝ ( T ) . Then the derivative of the equation ϝ ( N ) , ϝ ( T ) = 0 is calculated as ϝ ( N ) , ϝ ( T ) = ϝ ( N ) , ϝ ( T ) b = a . If we take ν S E = ϝ ( α ) S E = α S E * S E for the curve α S E = ( α 1 ( t ) r ( t ) σ 1 α 2 ( t ) r ( t ) σ 2 ) , then we can take the superelliptic unit tangent vector as ϝ ( T ) = ϝ ( α ( t ) ) ν S E . Then the vector ϝ ( N ) is obtained by rotating counterclockwise to the vector ϝ ( T ) as follows;
ϝ ( N ) = 1 ν S E ( α 2 ( t ) r ( t ) σ 2 , α 1 ( t ) r ( t ) σ 1 )
Accordingly, if the derivative of the equation ϝ ( α ( t ) ) = ν S E ϝ ( T ) is taken, we reach
ϝ ( α ( t ) ) = ν S E ϝ ( T ) + ν S E ϝ ( T ) = ν S E ϝ ( T ) + ν S E a ϝ ( N ) .
If the inner product of this equation with ϝ ( N ) is taken and using the equality κ S E , it becomes ϝ ( α ( t ) ) , ϝ ( N ) = ν S E and thus we obtain T S E * = ϝ ( T ) = ν S E κ S E N S E , N S E * = ϝ ( N ) = ν S E κ S E T S E , where curvature satisfies
κ S E ( t ) = Λ * det ( α S E * ( t ) , α S E * * ( t ) ) α S E * ( t ) S E 3 .

5. An Orthonormal Frame of a Superelliptic Curve in R S E 3

Consider a unit-speed superelliptic curve α S E with tangent vector ϝ ( α ) = α S E * = T S E . Then we have
T S E , T S E S E = 1 .
If we take the derivative of the last equation, we reach
T S E * , T S E S E + T S E , T S E * S E = 0 T S E * , T S E S E = 0 .
Therefore, if we normalize the vector T S E * and represent it by N S E , we get
N S E = T S E * T S E * S E = α S E * * α S E * * S E .
Hence, two mutually orthogonal unit vectors are obtained along the related superelliptic curve. By easily obtaining the third vector, which is perpendicular to these two vectors, through the vector cross product, we have three mutually perpendicular unit vectors on the curve α S E . We denote the third vector obtained by the cross product of the two vectors T S E and N S E , as B S E . Accordingly, we have
B S E = T S E × S E N S E .
Consequently, the superelliptic Frenet frame fields are given by { T S E , N S E , B S E } . For a unit-speed curve α and a related superelliptic curve ϝ ( α ) = α S E , we can express the following relationship under the transformation ϝ ,
T ϝ ( T ) = T S E N ϝ ( N ) = N S E B ϝ ( B ) = B S E
where { T S E , N S E , B S E } is an orthonormal frame for a unit-speed superelliptic curve and satisfies
T S E × S E N S E = B S E , N S E × S E B S E = T S E , B S E × T S E = N S E .
Example 3.
In R S E 3 , α S E ( t ) = 1 2 cos t r ( t ) σ 1 , 1 2 sin t r ( t ) σ 2 , 3 2 t r ( t ) σ 3 is a superelliptic curve, the orthonormal frame { T S E , N S E , B S E } satisfies the equations
T S E ( t ) = α S E * ( t ) = 1 2 sin t r ( t ) σ 1 , 1 2 cos t r ( t ) σ 2 , 3 2 r ( t ) σ 3 ,
and
N S E ( t ) = T S E * T S E * S E = cos t r ( t ) σ 1 , sin t r ( t ) σ 2 , 0 ,
then using the superelliptic vector product with T S E ( t ) and N S E ( t ) the superelliptic binormal vector B S E ( t ) is found as
B S E ( t ) = T S E ( t ) × S E N S E ( t ) = 3 2 sin t r ( t ) σ 1 , 3 2 cos t r ( t ) σ 2 , 1 2 r ( t ) σ 3 .
Similarly,
N S E ( t ) = B S E ( t ) × S E T S E ( t ) = cos t r ( t ) σ 1 , sin t r ( t ) σ 2 , 0
and
T S E ( t ) = N S E ( t ) × S E B S E ( t ) = 1 2 sin t r ( t ) σ 1 , 1 2 cos t r ( t ) σ 2 , 3 2 r ( t ) σ 3
can be calculated.
Figure 4 shows the motion of the Frenet frame vectors along a helical curve and the motion of superelliptic Frenet vector fields along the superelliptic curve associated with this curve. The first figure (Figure 4(a)) shows the motion of the standard Frenet frame along the helical curve, and as can be seen, the frame is an orthonormal frame in its motion along the curve. The second figure (Figure 4(b)) shows the motion of the superelliptic Frenet frame for different values of the Gielis formula. It can be observed that the length of the frame vectors along the superelliptic curve is scaled by the Gielis radial scaling function. However, according to the superelliptic inner product defined in the article, the superelliptic frame is still an orthonormal frame along the superelliptic curve.
Definition ;7.
Let α S E be a unit-speed curve with the Frenet frame apparatus { T E , N E , B E , κ E , τ E } . Then the superelliptic curvature κ E and torsion τ E of superelliptic curve can be calculated as
κ S E = T E * , N S E S E
and
τ S E = N E * , B S E S E .
Additionally, since
N S E = T S E * T S E * S E T S E * = T S E * S E N S E
it is also defined as
κ S E ( t ) = T S E * ( t ) , N S E ( t ) S E = T S E * ( t ) S E N S E ( t ) , N S E ( t ) S E = T S E * ( t ) S E = α S E * * ( t ) S E .
Theorem 3.
Let α S E be a unit-speed superelliptic curve with the superelliptic frame apparatus { T S E , N S E , B S E , κ S E , τ S E } . Then the superelliptic frame derivative formulas are computed by
T S E * N S E * B S E * = 0 κ S E 0 κ S E 0 τ S E 0 τ S E 0 T S E N S E B S E .
Proof. 
According to the definition of the curvature function, we have T S E * = κ S E N S E , T S E * = T S E * N S E , and κ S E ( t ) = T S E * ( t ) S E = α S E * * ( t ) S E . Let us write the vector N S E * in terms of the vectors { T S E = ϝ ( T ) , N S E = ϝ ( N ) , B S E = ϝ ( B ) } as N S E * = u ϝ ( T ) + v ϝ ( N ) + w ϝ ( B ) . If the derivative of the expression ϝ ( N ) , ϝ ( N ) S E = 1 is taken, it is clear that v = 0 which means that it is in the form of N S E * = u ϝ ( T ) + w ϝ ( B ) . If we take the inner product of this equation with ϝ ( T ) and since ϝ ( T ) ϝ ( B ) and ϝ ( T ) , ϝ ( T ) S E = 1 , we get
ϝ ( N ) , ϝ ( T ) S E = u .
Conversely
ϝ ( N ) , ϝ ( T ) S E = 0 ϝ ( N ) , ϝ ( T ) S E = ϝ ( N ) , ϝ ( T ) S E
and if we use the equation ϝ ( T ) = κ S E ϝ ( N ) = T S E * = κ S E N S E * , we obtain
u = ϝ ( N ) , ϝ ( T ) S E ϝ ( T ) , ϝ ( N ) S E = κ S E .
If similar operations are performed to find w, the result will be w = ϝ ( N ) , ϝ ( B ) = τ S E . In this case, we compute
N S E * = κ S E T S E + τ S E B S E .
Finally, using similar calculations, and using the definition of τ S E , we reach B S E * = τ S E N S E . □
Theorem 4.
Let us take the superelliptic curve α S E with the superelliptic frame apparatus { T S E , N S E , B S E , κ S E , τ S E } such that ν S E = ϝ ( α ) S E = α S E * S E . Then the superelliptic frame fields derivative formulas are given by
T S E * N S E * B S E * = 0 ν S E κ S E 0 ν S E κ S E 0 ν S E τ S E 0 ν S E τ S E 0 T S E N S E B S E .
Theorem 5.
Consider the superelliptic curve α S E with curvature κ S E and torsion τ S E . Then the following equations are satisfied
T S E = α S E * α S E * S E , B S E = α S E * × S E α S E * * α S E * × S E α S E * * S E , N S E = B S E × S E T S E κ S E = α S E * × S E α S E * * S E α S E * S E 3 , τ S E = α S E * × S E α S E * * , α S E * * * S E α S E * × S E α S E * * S E 2 .
Proof. 
Since ν S E = α S E * and α S E * = ν S E ϝ ( T ) , we compute:
ϝ ( T ) = α S E * ( t ) α S E * ( t ) S E
Furthermore, considering the following relations:
ϝ ( α ) = ν ϝ ( T ) + ν 2 κ S E ϝ ( N ) , ϝ ( T ) × S E ϝ ( T ) = 0 , ϝ ( T ) × S E ϝ ( N ) = ϝ ( B )
it follows that:
ϝ ( α ) × S E ϝ ( α ) = ν S E 3 κ S E B S E
By taking the norm of both sides and observing that B S E S E = 1 , we obtain:
ϝ ( α ) × S E α S E = ν 3 κ S E κ S E = ϝ ( α ) × S E α S E ϝ ( α ) 3 S E
On the other hand, given that ν S E 3 κ S E = ϝ ( α ) × S E α S E , the binormal vector is defined as:
B S E = α S E * × S E α S E * * α S E * × S E α S E * * S E
The principal normal vector N S E is then derived via the vector product:
N S E = B S E × S E T S E
If the third-order star derivative α S E * * * is expanded with respect to the superelliptic Frenet basis, we have:
α S E * * * = α S E * * * , T S E S E T S E + α S E * * * , N S E S E N S E + α S E * * * , B S E S E B S E
According to the Frenet formulas, we obtain:
α S E * * * , T S E S E = ν S E 3 κ S E 2 , α S E * * * , N S E S E = ν S E 3 κ S E , α S E * * * , B S E S E = ν S E 3 κ S E τ S E
Thus, the third-order derivative is computed as:
α S E * * * = ( ν 3 κ 2 ) T S E + ( ν 3 κ ) N S E + ( ν 3 κ τ ) B S E
From this result, the superelliptic torsion τ S E is derived as:
τ S E = α S E * × S E α S E * * , α S E * * * S E α S E * × S E α S E * * S E 2

5.1. Some Geometric Meaning of the Superelliptic Curvatures

Theorem 6.
A superelliptic curve α S E is a superelliptic line; the necessary and sufficient condition is that it is κ S E = 0 .
Proof. ( ) If the superelliptic curve is a superelliptic line, it can be written in the form α S E ( t ) = p + t t , where t is a constant superelliptic unit vector. Consequently, Consequently, since the vector t is constant, we compute T S E * = 0 and thus we get κ S E = T S E * = 0 .
( ) Let the curvature κ S E of the unit-speed superelliptic curve α S E be zero, ie., κ S E = 0 . According to the orthonormal frame equations, we obtain
T S E * = κ S E N S E = 0 .
Therefore, the superelliptic tangent vector T S E is a constant superelliptic unit vector. Since α S E * = T S E is constant, we reach
α S E ( t ) = t T S E + a
and this is the equation of the superelliptic line whose direction is T S E and which passes through the point a. So, α S E is a superelliptic line. □
Next, we introduce the following definition, which will be used to prove Theorem 7.
Definition 8.
Assume that α S E is a superelliptic curve in R S E 3 such that all the points on this superelliptic curve have a constant distance r S E from the constant point m S E . Then, this curve is called a superelliptic circle with the center m S E and the radius r S E , satisfying as follows:
m S E ( t ) = α S E ( t ) + 1 κ S E ( t ) N S E ( t ) r S E ( t ) = 1 κ S E .
Theorem 7.
A superelliptic curve α S E is a superelliptic circle if and only if κ S E is constant and τ S E = 0 .
Proof. (⇒)Let α S E be a unit-speed superelliptic circle with center m S E and radius r S E . Then, we know that α S E is a superelliptic planar curve and so τ S E = 0 . Therefore, we must show that κ S E is a positive constant. A circle with radius r S E and center point m S E can be written as
α S E ( t ) m S E S E = r S E
and from here we have
α S E ( t ) m S E , α S E ( t ) m S E S E = r S E 2
If we take the star derivative of both sides of the equation, we get
2 α S E * ( t ) , α S E ( t ) m S E S E = 0 T S E , α S E ( t ) m S E S E = 0 .
If we take the star derivative again, we find
T S E * , α S E ( t ) m S E S E + T S E , T S E S E = κ S E N S E , α S E ( t ) m S E S E + 1 = 0
According to the last equation, we obtain
κ S E > 0 and N S E , α S E ( t ) m S E S E 0
Taking the star derivative of the last equation, we compute
κ S E * N S E , α S E ( t ) m S E S E + κ S E N S E * , α S E ( t ) m S E S E + κ S E N S E , T S E S E = 0
Then, since N S E , T S E S E = 0 and N S E * = κ S E T S E + τ S E B S E , we reach
κ S E * N S E , α S E ( t ) m S E S E κ S E 2 T S E , α S E ( t ) m S E S E + κ S E τ S E B S E , α S E ( t ) m S E S E = 0
Since T S E , α S E ( t ) m S E S E = 0 and α S E ( t ) is planar, then we have τ S E = 0 . Thus, we obtain
κ S E * N S E , α S E ( t ) m S E S E = 0 .
Since
N S E , α S E ( t ) m S E S E 0
it becomes κ S E * = 0 . So, we calculate that κ S E is constant.
(⇐) Let κ S E be constant and τ S E =0. To show that the superelliptic curve α S E is a unit-speed circle, we must show that the distance from the point α S E ( t ) to a fixed point m S E is equal to a constant number r S E . First, let us see that the center m S E is fixed. In equation
m S E ( t ) = α S E ( t ) + 1 κ S E ( t ) N S E ( t ) r S E ( t ) = 1 κ S E ,
since τ S E = 0 , κ S E is a constant and N S E * = κ S E T S E + τ S E B S E , we find that m S E * = 0 . Therefore, m S E is a constant. On the other hand, let us show that the distance from every point on the curve to the point m S E is r S E ( t ) = 1 κ S E . Indeed,
α S E m S E = 1 κ S E N S E
is found to be α S E m S E S E = 1 κ S E N S E = r S E . Therefore, α S E is a circle. □
Example 4.
Let us see α S E = ( ( 1 cos t ) r ( t ) σ 1 , ( 2 + cos t ) r ( t ) σ 2 , 2 sin t r ( t ) σ 3 ) ) is a superelliptic circle by finding κ S E and τ S E . If we take the following star derivatives
α S E * = ( sin t r ( t ) σ 1 , sin t r ( t ) σ 2 , 2 cos t r ( t ) σ 3 ) ) , α S E * * = ( cos t r ( t ) σ 1 , cos t r ( t ) σ 2 , 2 sin t r ( t ) σ 3 ) ) .
then we compute the cross product as
α S E * × S E α S E * * = Λ * r ( t ) 2 e 1 σ 1 r ( t ) 2 e 2 σ 2 r ( t ) 2 e 3 σ 3 sin t r ( t ) σ 1 sin t r ( t ) σ 2 2 cos t r ( t ) σ 3 cos t r ( t ) σ 1 cos t r ( t ) σ 2 2 sin t r ( t ) σ 3 = ( 2 , 2 , 0 )
Thus we calculate α S E * S E = 2 and α S E * × S E α S E * * S E = 2 and we reach
κ S E = α S E * × S E α S E * * S E α S E * S E 3 = 1 2 .
If we take the third star derivative
α S E * * * = ( sin t r ( t ) σ 1 , sin t r ( t ) σ 2 , 2 cos t r ( t ) σ 3 )
we calculate
τ S E = α S E * × S E α S E * * , α S E * * * S E α S E * × S E α S E * * S E 2 = 0 .
Definition 9.
In R S E 3 , a superelliptic curve α S E is a helix such that κ S E 0 if and only if the superelliptic tangent vector T S E makes a constant angle with a constant vector U S E . Accordingly, a superelliptic helix is defined by
  • U S E * = 0
  • T S E , U S E S E = cos ϕ .
Theorem 8.
Consider the superelliptic curve α S E with curvature κ S E and torsion τ S E . Then α S E is a superelliptic helix if and only if the ratio τ S E κ S E is constant.
Proof. (⇒)According to the definition above, if α is a helix, there exists a U S E vector that satisfies the conditions U S E * = 0 and we know that T S E ( t ) , U S E S E = c o s ϕ and T * ( t ) , U S E = 0 . Then we get κ S E N S E , U S E S E = 0 , which gives N S E U S E . In this case, since the vector U S E is in the plane formed by T S E and B S E , it is written as
U S E = cos ϕ T S E + sin ϕ B S E .
If the star derivative of both sides of the equation is taken and since U S E * = 0 , we compute
cos ϕ κ S E N S E sin ϕ τ S E N S E = 0 .
Thus, it is found that the ratio τ S E κ S E is constant.
(⇐) Now assume that the ratio τ S E κ S E is constant. We claim that U S E is a constant vector and the superelliptic curve α S E makes a constant angle with the superelliptic tangent vector T S E . Then we have
U S E * = cos ϕ κ S E N S E sin ϕ τ S E N S E .
If τ S E = κ S E cot ϕ is written, then U S E * = 0 is obtained. Therefore, the vector U S E is constant. □
Theorem 9.
Let α be a helix in R 3 and α S E be a related superelliptic curve. Then the superelliptic curve α S E is a superelliptic helix. The converse of this theorem is also true.
Proof. (⇒)According to the definition above, if α is a helix, there exists a U vector that satisfies condition U = 0 and we know that T ( t ) , U = c o s ϕ and T ( t ) , U = 0 . If U = 0 , then we have ϝ ( U ) = U S E * = 0 . On the other hand, since ϝ is an isometry, we have
T S E * ( t ) , U S E S E = ϝ ( T ( t ) ) , U S E S E = T ( t ) , U = 0 .
In this case, we obtain T S E * ( t ) , U S E S E = c o n s t a n t . Therefore, α S E is a superelliptic helix. Similarly, we can prove that the converse holds as well. □
Example 5.
Let us construct the superelliptic frame apparatus { T S E , N S E , B S E , κ S E , τ S E } structure for the unit-speed superelliptic helix α S E = ( cos t 2 r ( t ) σ 1 , sin t 2 r ( t ) σ 2 , t r ( t ) 2 σ 3 ) . If we take the star derivative, we get
α S E * = ( 1 2 sin t 2 r ( t ) σ 1 , 1 2 cos t 2 r ( t ) σ 2 , r ( t ) 2 σ 3 )
such that α S E * S E = 1 . Thus, we compute
T S E = α S E * = ( 1 2 sin t 2 r ( t ) σ 1 , 1 2 cos t 2 r ( t ) σ 2 , r ( t ) 2 σ 3 ) , N S E = α S E * * α S E * * S E = ( cos t 2 r ( t ) σ 1 , sin t 2 r ( t ) σ 2 , 0 ) ,
and
B S E = T S E × S E N S E = Λ * r ( t ) 2 e 1 σ 1 r ( t ) 2 e 2 σ 2 r ( t ) 2 e 3 σ 3 1 2 sin t 2 r ( t ) σ 1 1 2 cos t 2 r ( t ) σ 2 r ( t ) 2 σ 3 cos t 2 r ( t ) σ 1 sin t 2 r ( t ) σ 2 0 = ( 1 2 sin t 2 r ( t ) σ 1 , 1 2 cos t 2 r ( t ) σ 2 , r ( t ) 2 σ 3 ) .
Next, we calculate κ S E and τ S E for the superelliptic helix α S E . For these, we compute
T S E * = ( 1 2 cos t 2 r ( t ) σ 1 , 1 2 sin t 2 r ( t ) σ 2 , 0 ) and N S E = ( cos t 2 r ( t ) σ 1 , sin t 2 r ( t ) σ 2 , 0 ) . Then the superelliptic curvature is calculated as
κ S E = T S E * , N S E S E = 1 2 .
If we use
N S E * = ( 1 2 sin t 2 r ( t ) σ 1 , 1 2 cos t 2 r ( t ) σ 2 , 0 )
and
B S E = ( 1 2 sin t 2 r ( t ) σ 1 , 1 2 cos t 2 r ( t ) σ 2 , r ( t ) 2 σ 3 )
the superelliptic torsion is obtained as
τ S E = N S E * , B S E S E = 1 2 .

5.2. The Superelliptic Darboux Vector of the Superelliptic Orthonormal Frame

Definition 10.
The rate at which a point rotates about a center, that is, the rate of change of its angular position, is briefly called angular speed, and during this rotation, there is also a rotation axis. This axis of rotation is called the angular velocity vector. This vector is denoted by ω . The angular velocity vector of the superelliptic orthonormal frame of a space curve is called the superelliptic Darboux vector. Accordingly, the vector ω that satisfies the equation
T S E * = ω × S E T S E N S E * = ω × S E N S E B S E * = ω × S E B S E
for the T S E , N S E , B S E orthogonal vector fields of an α S E curve is called the superelliptic Darboux vector of the superelliptic orthonormal frame. This vector is usually denoted by the letter D S E .
Theorem 10.
Let D S E be the superelliptic Darboux vector of the T S E , N S E , B S E orthonormal frame of any unit-speed superelliptic curve α S E . It is determined by
D S E = τ S E T S E + κ S E B S E ,
Proof. 
Assume that D S E = u T S E + v N S E + w B S E .
T S E * = D S E × S E T S E κ S E N S E = ( u T S E + v N S E + w B S E ) × S E T S E
κ S E N S E = v B S E + w N S E
from the equality, v = 0 and w = κ S E are found. Therefore, it has to be structured as D S E = u S E T S E + κ S E B S E . On the other hand
B S E * = D S E × S E B S E τ N S E = ( u S E T S E + κ S E B S E ) × S E B S E
τ N S E = u N S E
in the equation, u = τ S E is obtained. Consequently, D S E = τ S E T S E + κ S E B S E is acquired. □

6. Superelliptic Surface

As is well known, in this study a conformal elliptic metric defined on the space R n in the following form is used:
u , v S E = 1 r ( t ) 2 i = 1 n σ i u i v i ,
where σ i > 0 are positive constant coefficients and r ( t ) is a direction-dependent radial function. This metric can be interpreted as a conformal scaling of the elliptic inner product.
In the previous section, a curve on the elliptic sphere, α ( t ) = ( α 1 ( t ) , , α n ( t ) ) , was defined with respect to the superelliptic metric in the following form:
α S E ( t ) = α 1 ( t ) r ( t ) σ 1 , α 2 ( t ) r ( t ) σ 2 , , α n ( t ) r ( t ) σ n .
This definition expresses the normalization of the curve under the superelliptic metric and its mapping onto the superelliptic sphere. Using the superelliptic metric and the mapping;
ϝ : R n R S E n ,
by
ϝ ( α ( t ) ) = α S E = α 1 ( t ) r ( t ) σ 1 , α 2 ( t ) r ( t ) σ 2 , , α n ( t ) r ( t ) σ n .
We can give the following definition.
Definition 11.
Let S : U R 2 R 3 ; ( t , v ) S ( t , v ) = ( f 1 ( t , v ) , f 2 ( t , v ) , f ( t , v ) ) be a smooth surface. In R 3 S E , the superelliptic surface is defined via the following Gielis radial scaling transformation ϝ as follows:
ϝ : R 3 R S E 3 S ( t , v ) ϝ ( S ( t , v ) ) = X ( t , v ) = ( f 1 ( t , v ) r ( t ) σ 1 , f 2 ( t , v ) r ( t ) σ 2 , f 3 ( t , v ) r ( t ) σ 3 )
where r ( t ) is the Gielis’ super formula and σ i , i { 1 , 2 , 3 } are positive real numbers.
Next, we introduce distinguished surfaces within superelliptic three-dimensional space, characterized by their specific geometric and algebraic properties. These surfaces play a significant role in understanding the structural behavior and symmetry conditions inherent to superelliptic geometries

6.1. Superelliptic Sphere

In space, under the superelliptic metric, the geometric locus consisting of all points whose distance to a fixed point is constant is called a superelliptic sphere. This fixed point is referred to as the center of the superelliptic sphere, and the constant distance is called the superelliptic radius. Since the superelliptic metric defines a direction-dependent and scaled notion of distance, the resulting superelliptic sphere exhibits anisotropic geometric properties when compared to the classical sphere. For appropriate choices of the metric and parameters, this definition reduces to the classical sphere, whereas in the general case, the superelliptic sphere represents a broader class of spherical surfaces characterized by superelliptic structures.
The superelliptic sphere (superellipsoid) is obtained as
X ( t , v ) = r ( t ) r ( v ) cos t sin v σ 1 , r ( t ) r ( v ) sin t sin v σ 2 , r ( t ) r ( v ) cos v σ 3 .
The superellipsoids obtained for various values of the Gielis formula with r ( t ) , r ( v ) are illustrated in Figure 5.

6.2. Superelliptic Cylindrical Surface

In superelliptic space, the surface generated by straight lines that move along a fixed superelliptic curve and remain parallel to a fixed direction is called a superelliptic cylindrical surface. In this definition, the fixed superelliptic curve is called the directrix, and the fixed direction is referred to as the director of the superelliptic cylinder.
According to this definition, let the superelliptic directrix be
α S E = α 1 ( t ) r ( t ) σ 1 , α 2 ( t ) r ( t ) σ 2 , α 3 ( t ) r ( t ) σ 3
and the director be
w = a r ( t ) σ 1 , b r ( t ) σ 2 , c r ( t ) σ 3 .
Then, the parametric equation of the corresponding superelliptic cylindrical is obtained as
X ( t , v ) = ( α 1 ( t ) + v a ) r ( t ) σ 1 , ( α 2 ( t ) + v b ) r ( t ) σ 2 , ( α 3 ( t ) + v c ) r ( t ) σ 3 .
Here, a , b , c are real constants. Examples of superelliptic cylinders for some values of r are given in Figure 6.
Example 6.
Let us determine the parametric equation of the superelliptic cylinder with the base curve α S E ( t ) = ( cos t r ( t ) σ 1 , sin t r ( t ) σ 2 , 0 ) and the direction vector w = ( 0 , 0 , r ( t ) σ 3 ) . The parametric equation for v r ( v ) = ( | cos ( 6 v / 4 ) | 0.28 + | sin ( 6 v / 4 ) | 0.28 ) 1 / 0.28 is,
X ( t , v ) = α S E ( t ) + r ( v ) w = cos t r ( t ) σ 1 , sin t r ( t ) σ 2 , r ( v ) r ( t ) σ 3 , ( t , v ) [ 0 , 2 π ] × [ 0 , 2 π ] .
and it is given in Figure 7

6.3. Superelliptic Conical Surface

In space, the surface generated by a superelliptic line passing through a fixed point and moving along a superelliptic curve is called a superelliptic conical surface. The curve lying on the superelliptic curve and forming the base of the surface is called the base curve, and the fixed point of the cone is called the vertex. The moving line is referred to as the generatrix of the superelliptic cone.
Let the base curve be
α S E = α 1 ( t ) r ( t ) σ 1 , α 2 ( t ) r ( t ) σ 2 , α 3 ( t ) r ( t ) σ 3
and let the vertex point on the curve α S E be
P = a r ( t ) σ 1 , b r ( t ) σ 2 , c r ( t ) σ 3 .
Then, the equation of the superelliptic conical surface is given by
X ( u , t ) = a r ( t ) σ 1 , b r ( t ) σ 2 , c r ( t ) σ 3 + v ( α 1 ( t ) a ) r ( t ) σ 1 , ( α 2 ( t ) b ) r ( t ) σ 2 , ( α 3 ( t ) c ) r ( t ) σ 3 .
Example 7.
Let us take the base curve
α S E = r ( t ) σ 1 cos t , r ( t ) σ 2 sin t , 0
and the vertex point
P = 0 , 0 , r ( t ) σ 3 .
Then, we obtain the superelliptic conical surface as
α S E = r ( t ) σ 1 v cos t , r ( t ) σ 2 v sin t , v r ( t ) σ 3 .
Some superelliptic conical surfaces corresponding to different values of the superformula are presented in Figure 8.

6.4. Superelliptic Surfaces of Revolution

Surfaces obtained by rotating a superelliptic curve about a specified axis are called superelliptic surfaces of revolution. For example, a superellipsoid can be regarded as a surface of revolution obtained by rotating a superellipse. A curve given in the x z -plane by
α ( v ) = ( r ( v ) f ( v ) , 0 , r ( v ) g ( v ) )
when rotated superelliptically (see [11]) about the z-axis, yields a surface of revolution whose parametric representation is
X ( t , v ) = r ( v ) f ( v ) cos t r ( t ) σ 1 , r ( v ) f ( v ) sin t r ( t ) σ 2 , r ( v ) g ( v ) r ( t ) σ 3 , v [ 0 , 2 π ) .
Example 8.
Consider a superelliptic circle of radius a whose center is located at ( R , 0 ) in the x z -plane.
x ( v ) = ( R + a cos v ) r ( v ) , z ( v ) = a sin v r ( v ) ,
Under the superelliptic rotation, a point ( x ( v ) , 0 , z ( v ) ) about the angle t [ 0 , 2 π ) is mapped to
X ( t , v ) = x ( v ) cos t r ( t ) σ 1 , y ( t , v ) = x ( v ) sin t r ( t ) σ 2 , z ( t , v ) = z ( v ) r ( t ) σ 3 .
Substituting the expressions of x ( v ) and z ( v ) yields
X ( t , v ) = r ( v ) ( R + a cos v ) cos t r ( t ) σ 1 , y ( t , v ) = r ( v ) ( R + a cos v ) sin t r ( t ) σ 2 , z ( t , v ) = a sin v r ( t ) r ( v ) σ 3 .
X ( t , v ) = r ( v ) ( R + a cos v ) cos t r ( t ) σ 1 , r ( v ) ( R + a cos v ) sin t r ( t ) σ 2 , a sin v r ( t ) r ( v ) σ 3 .
For different choices of the parameters in r, the corresponding superelliptic tori are illustrated in Figure 9.
Example 9.
Let us consider the curve in the x z -plane given by
α S E ( v ) = a cosh v a , 0 , v
and rotate this curve superelliptically about the z-axis by an angle v. As a result of this rotation, the parametric representation of the resulting superelliptic surface of revolution is given by
X ( t , v ) = a r ( t ) σ 1 cosh v a cos t , a r ( t ) σ 2 cosh v a sin t , r ( t ) σ 3 v .
This surface is referred to as thesuperelliptic catenoidand is illustrated in Figure 10. Superelliptic catenoids are also Constant Anisotropic Mean Curvature Surfaces [16].

6.5. Superelliptic Ruled Surface

These surfaces are called superelliptic ruled surfaces. In other words, the surfaces obtained by sweeping superelliptic lines that move along a superelliptic curve and whose directions may vary depending on the parameter are superelliptic ruled surfaces.
For any differentiable superelliptic curve, let α S E be the base curve and γ S E be the director curve. Then, the ruled surface is given by
X ( t , v ) = α S E ( t ) + v γ S E ( t ) .
Example 10.
The superelliptic helicoid surface given by the parametrization
X ( t , v ) = v cos t r ( t ) σ 1 , v sin t r ( t ) σ 2 , t r ( t ) σ 3
can be written as
X ( t , v ) = 0 , 0 , t r ( t ) σ 3 + v cos t r ( t ) σ 1 , sin t r ( t ) σ 2 , 0 .
This surface is a ruled surface and is illustrated in Figure 11.
Example 11.
Let the superelliptic Möbius strip be defined by the parametrization
X ( t , v ) = ( cos t + v cos ( t 2 ) cos t ) r ( t ) σ 1 , ( sin t + v cos ( t 2 ) sin t ) r ( t ) σ 2 , v sin ( t 2 ) r ( t ) σ 3
This expression can be written in the affine form
X ( t , v ) = cos t r ( t ) σ 1 , sin t r ( t ) σ 2 , 0 + v cos ( t 2 ) cos t r ( t ) σ 1 , cos ( t 2 ) sin t r ( t ) σ 2 , sin ( t 2 ) r ( t ) σ 3 .
Therefore, X is aruled surface: for each fixed t, the map is linear in v and parametrizes a straight ruling in R 3 . The corresponding superelliptic Möbius strip is depicted in Figure 12.

6.6. Superelliptic Gauss and Mean Curvatures

The surface normal is a vector perpendicular to the surface’s tangent plane spanned by the vectors { X t * , X v * } at a given point, in a superelliptic sense. For a parametrized surface X ( t , v ) , the normal vector is given by
N = X t * × S E X v * ,
and the unit normal is
n = X t * × S E X v * X t * × S E X v * S E .
For a parametrized surface X ( t , v ) , the first fundamental form describes the metric properties (lengths and angles) on the surface:
I S E = E S E d t 2 + 2 F S E d t d v + G S E d v 2 ,
where
E S E = X t * , X t * S E , F S E = X t * , X v * S E , G S E = X v * , X v * S E .
The second fundamental form
I I S E = L S E d t 2 + 2 M S E d t d v + N S E d v 2 .
describes how the surface bends in space, where
L S E = X t t * , n S E , M S E = X t v * , n S E , N S E = X v v * , n S E .
Gaussian curvature can be expressed using the fundamental form coefficients as
K S E = L S E N S E M S E 2 E S E G S E F S E 2 .
It indicates whether the superelliptic surface is locally shaped like a superelliptic sphere ( K S E > 0 ), a superelliptic saddle ( K S E < 0 ), or a superelliptic plane ( K S E = 0 ). The mean curvature describes how the surface bends on average and appears in many geometric and physical problems such as minimal surfaces, which is given by
H S E = E S E N S E 2 F S E M S E + G S E L S E 2 ( E S E G S E F S E 2 ) .
Here, the coefficients E S E , F S E , G S E represent the intrinsic metric properties of the superelliptic surface S E , characterizing the first fundamental form. In contrast, the coefficients L S E , M S E , N S E correspond to the second fundamental form and describe how the surface bends within the ambient space R S E 3 .
A surface is said to be minimal if its mean curvature satisfies H S E = 0 at every point. As a canonical example, the superelliptic catenoid constitutes a minimal surface, since its mean curvature vanishes identically across the entire surface.

6.7. Superelliptic Darboux Frame Apparatus in R S E 3

A natural moving frame constructed on a surface is known as a Darboux frame. Let X be a superelliptic surface in R S E 3 and α S E : I R R S E 3 be a superelliptic curve on the superelliptic surface. We denote the superelliptic normal vector field with n S E and the superelliptic tangent vector with T S E , then since they are perpendicular to each other, the n S E × S E T S E vector is perpendicular to both the vectors n S E and T S E . If we say n S E × S E T S E = g S E , we obtain a superelliptic orthonormal basis on the surface consisting of the vectors { T S E , n S E , g S E } along the curve α S E . Moreover, for a unit-speed curve α and a related superelliptic curve ϝ ( α ) = α S E , we can express the following relationship under the transformation ϝ ,
T ϝ ( T ) = T S E n ϝ ( n ) = n S E g ϝ ( g ) = g S E
where { T S E , n S E , g S E } is a Darboux frame for a unit-speed superelliptic curve and satisfies
T S E × S E n S E = g S E , n S E × S E T S E = g S E , g S E × T S E = n S E , n S E × S E b S E = T S E .
Additionally, the function that measures how much a superelliptic curve bends toward its surface normal is called the superelliptic normal curvature κ n S E and can be calculated as follows:
κ n S E = α S E * * , n S E α S E * 2
The function that measures the deviation of a superelliptic curve from being a straight line on the surface is called the superelliptic geodesic curvature κ g S E and can be calculated as follows:
κ g S E = α S E * × S E α S E * * , n S E S E α S E * S E 3 .
The function showing how the direction of a superelliptic curve on a surface rotates with respect to the curvature direction of the surface is called the superelliptic geodesic torsion τ g S E and can be calculated as follows:
τ g S E = n S E * × S E n S E , α S E * S E α S E * S E 2
On the other hand, if α S E : I R R S E n is a unit-speed superelliptic curve, then we compute the superelliptic normal curvature, superelliptic geodesic curvature, and superelliptic geodesic torsion, respectively, as follows:
κ n S E = T S E * , n S E S E κ g S E = T S E * , g S E S E τ g S E = n S E * , g S E S E .
Next, we give the formulas for the star derivatives of the moving frame { T S E , n S E , g S E } .
Theorem 11.
Let α S E be a unit-speed superelliptic curve on the superelliptic surface X, and { T S E , n S E , g S E , κ n S E , κ g S E , τ g S E } be the Darboux frame apparatus. Then, the superelliptic Darboux frame derivative formulas are computed by
T S E * n S E * g S E * = 0 κ n S E κ g S E κ n S E 0 τ g S E κ g S E τ g S E 0 T S E n S E g S E .
Proof. 
We can write T S E * = x 1 T S E + x 2 n S E + x 3 g S E . If the superelliptic inner product of this equation with T S E , n S E , g S E is taken, respectively, it becomes
x 1 = T S E * , T S E S E = 0 , x 2 = T S E * , n S E S E = κ n S E , x 3 = T S E * , g S E S E = κ g S E .
Thus we can write
T S E * = κ n S E n S E + κ g S E g S E .
Similarly, if we write n S E * = x 1 T S E + x 2 n S E + x 3 g S E and if the superelliptic inner product of this equation with T S E , n S E , g S E is taken, respectively, it becomes
y 2 = n S E * , n S E S E = 0 , y 1 = T S E , n S E * S E = T S E * , n S E S E = κ n S E , x 3 = g S E , n S E * S E = τ g S E ,
then we compute
n S E * = κ n S E T S E + τ g S E g S E .
Finally, the calculation is done in a similar way for g S E , we obtain
g S E * = κ g S E T S E τ g S E n S E .
These completes the proof. □
Example 12.
Let us take the superelliptic cyclinder is
X ( t , v ) = ( r ( t ) σ 1 cos t , r ( t ) σ 2 sin t , v r ( t ) σ 3 ) .
In this case, the superelliptic moving frame along the superelliptic cyclinder are computed as follows: For X t * = ( r ( t ) σ 1 sin t , r ( t ) σ 2 cos t , 0 ) and X v * = ( 0 , 0 , r ( t ) σ 3 ) we compute the superelliptic unit normal vector
n S E = X t * × S E X v * = ( r ( t ) σ 1 cos t , r ( t ) σ 2 sin t , 0 ) ,
superelliptic tangent vector
T S E = α S E * α S E * S E = ( r ( t ) 2 σ 1 sin t , r ( t ) 2 σ 2 cos t , r ( t ) 2 σ 3 ) ,
and superelliptic vector
g S E = T S E × S E n S E = ( r ( t ) 2 σ 1 sin t , r ( t ) 2 σ 2 cos t , 1 2 ) .
Then, the superelliptic normal curvature, superelliptic geodesic curvature, and superelliptic geodesic torsion are calculated as follows:
κ n S E = α S E * * , n S E S E α S E * S E 2 = 1 2 ,
κ g S E = α S E * × S E α S E * * , n S E S E α S E * S E 3 = 1 2 2 ,
τ g S E = n S E * × S E n S E , α S E * S E α S E * S E 2 = 1 2 .

7. Conclusion

In this study, we introduced a novel geometric framework that integrates Gielis’ superformula with the classical differential geometry of curves and surfaces. By defining the superelliptic inner and cross products, as well as the star derivative and the superelliptic Frenet frame, we extended the Euclidean and Riemannian interpretations of curvature and torsion to a more parametric and flexible geometric structure.
The proposed superelliptic geometry provides a new perspective in modeling natural and organic forms, where symmetry, scaling, and deformation can be controlled through a small set of parameters. This structure preserves essential geometric properties such as orthonormality and differentiability while allowing an infinite family of curves and surfaces to be generated from the same analytical foundation.
The framework definition presented in this study demonstrates that Gielis curves can be treated within the same formal and geometric characterization, independent of their classical parametric expressions. The resulting framework reduces even singular Gielis curves to the curvature and torsion properties of a circle, thus providing a "regulatory" layer that eliminates the geometric effects of singularities. This approach unifies all members of the Gielis family, including generalized Möbius-Listing surfaces [35], within a common framework, providing a more robust foundation for both computational geometry and natural shape modeling.
Finally, it is worth emphasizing that the proposed framework is highly general. In particular, when the superformula is considered as an illustrative example, it can be shown that, through the use of the star derivative, a wide variety of geometric structures and families of curves can be represented within the same formalism. This observation indicates that the proposed approach is not restricted to a specific special case; rather, it possesses the flexibility to encompass a broad class of functions and shapes. Hence, the superformula, analyzed via the star derivative, provides a representative example that highlights the generality and potential applicability of the developed framework. The results show that the superelliptic framework offers a consistent and computationally tractable model for analyzing complex geometries, bridging the gap between mathematical abstraction and real-world morphology.
Future work may include the extension of this framework to non-Euclidean manifolds and applications in physics, architecture, and biological modeling. Furthermore, this framework paves the way for the development of new differential-geometric methods that can be applied to Gielis surfaces, growth models, and minimal surfaces in the future.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Do Carmo, M.P. Differential geometry of curves and surfaces; Prentice-Hall: New Jersey, USA, 1979. [Google Scholar]
  2. Gielis, J. Wiskundige supervormen bij bamboes. Newsl. Belg. Bamboo Soc. 1996, 13, 20–26. [Google Scholar]
  3. Gielis, J. A generic geometric transformation that unifies a wide range of natural and abstract shapes. Am. J. Bot. 2003, 90(3), 333–338. [Google Scholar] [CrossRef]
  4. Gielis, J.; Beirinckx, B.; Bastiaens, E. Superquadrics with rational and irrational symmetry. In Proceedings of the eighth ACM symposium on Solid modeling and applications, Seattle, WA, USA, 16–20 June 2003; pp. 262–265. [Google Scholar]
  5. Gray, A. Modern differential geometry of curves and surfaces with Mathematica, 2nd ed.; CRC Press: Boca Raton, FL, USA, 1998. [Google Scholar]
  6. Gardner, M. Piet Hein’s Superellipse, Mathematical Carnival. A New Round-Up of Tantalizers and Puzzles from Scientific American; Vintage Press: New York, NY, USA, 1977. [Google Scholar]
  7. Kühnel, W. Differential geometry, 3rd ed.; American Mathematical Soc.: Providence, RI, USA, 2015; Vol. 77. [Google Scholar]
  8. Lamé, G.; Courcier, Mme. Ve. Examen des différentes méthodes employées pour résoudre les problémes de géométrie; Imprimeur-Libraire: Paris, France, 1818. [Google Scholar]
  9. O’Neill, B. Elementary differential geometry, 2nd ed.; Academic Press: San Diego, CA, USA, 1997. [Google Scholar]
  10. Özdemir, Z.; Parlak, E. Superquadric Motion and Superquadric Hyperbolic Split Quaternion Algebra Via Gielis Formula. Math. Meth. Appl. Sci. 2025, 48(16), 15035–15048. [Google Scholar] [CrossRef]
  11. Parlak, E.; Özdemir, Z. Superelliptic quaternions and superelliptic rotations with applications. Phys. Scr. 2025, 100(9), 095212. [Google Scholar] [CrossRef]
  12. Brandi, P.; Salvadori, A. A magic formula of nature. In Proceedings of the APLIMAT 17th Conference on Applied Mathematics, Bratislava, Slovakia, 6–8 February 2018; pp. 110–126. [Google Scholar]
  13. Gielis, J. The Geometrical Beauty of Plants; Atlantis/Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  14. Gielis, J. De Uitvinding van de Cirkel (Dutch, 2001); Inventing the Circle (English, 2003); Geniaal 557 Publishing: Antwerp, Belgium, 2001. [Google Scholar]
  15. Gielis, J.; Haesen, S.; Verstraelen, L. Universal natural shapes: from the supereggs of Piet Hein to the cosmic egg of Georges Lemaître. Kragujevac J. Math. 2005, 28, 57–68. [Google Scholar]
  16. Koiso, M.; Palmer, B. Equilibria for anisotropic surface energies and the Gielis formula. Forma 2008, 23, 1–8. [Google Scholar]
  17. Verstraelen, L. A concise mini history of geometry. Kragujevac J. Math. 2014, 38, 5–21. [Google Scholar] [CrossRef]
  18. Natalini, P.; Patrizi, R.; Ricci, P.E. The Dirichlet problem for the Laplace equation in a starlike domain of a Riemann surface. Numer. Algorithms 2008, 49(1-4), 299–313. [Google Scholar] [CrossRef]
  19. Matsuura, M. Gielis’ superformula and regular polygons. J. Geom. 2015, 106(2), 383–403. [Google Scholar] [CrossRef]
  20. Gielis, J.; Caratelli, D.; Fougerolle, Y.; Ricci, P.E.; Tavkelidze, I.; Gerats, T. Universal natural shapes: From unifying shape description to simple methods for shape analysis and boundary value problems. PloSOne 2012, 7(9), e29324. [Google Scholar] [CrossRef]
  21. Gielis, J.; Natalini, P.; Ricci, P.E. A note about generalized forms of the Gielis formula. In Modeling in Mathematics: Proceedings of the Second Tbilisi-Salerno Workshop on Modeling in Mathematics; Atlantis Press: Dordrecht, The Netherlands, 2017; pp. 107–116. [Google Scholar]
  22. Haesen, S.; Nistor, A.; Verstraelen, L. On growth and form and geometry I. Kragujevac J. Math. 2012, 36, 5–25. [Google Scholar]
  23. Shi, P.J.; Huang, J.G.; Hui, C.; Grissino-Mayer, H.D.; Tardif, J.C.; Zhai, L.H.; Wang, F.S.; Li, B.L. Capturing spiral radial growth of conifers using the superellipse to model tree-ring geometric shape. Front. Plant Sci. 2015, 6, 856. [Google Scholar] [CrossRef]
  24. Faisal, T.R.; Hristozov, N.; Western, T.L.; Rey, A.; Pasini, D. The twistto-bend compliance of the Rheum rhabarbarum petiole: integrated computations and experiments. Comput. Methods Biomech. Biomed. Eng. 2016, 19(15), 1–12. [Google Scholar]
  25. Shi, P.; Ratkowsky, D.A.; Li, Y.; Zhang, L.; Lin, S.; Gielis, J. A general leaf area geometric formula exists for plants—Evidence from the simplified Gielis equation. Forests 2018, 9, 714. [Google Scholar] [CrossRef]
  26. Huang, W.; Li, Y.; Niklas, K.J.; Gielis, J.; Ding, Y.; Cao, L.; Shi, P. A superellipse with deformation and its application in describing the cross-sectional shapes of a square bamboo. Symmetry 2020, 12, 2073. [Google Scholar] [CrossRef]
  27. Gielis, J.; Shi, P.; Beirinckx, B.; Caratelli, D.; Ricci, P.E. Lamé-Gielis curves in biology and geometry. In Proceedings of the International Conference Riemannian Geometry and Applications (RIGA 2021), Bucharest, Romania, 15–17 January 2021. [Google Scholar]
  28. Mescia, L.; Bia, P.; Caratelli, D.; Chiapperino, M.A.; Stukach, O.; Gielis, J. Electromagnetic mathematical modeling of 3D supershaped dielectric lens antennas. Math. Probl. Eng. 2016, 2016, 8130160. [Google Scholar] [CrossRef]
  29. Geerarts, J.M.; Theis, G.; Smolders, A.B.; Caratelli, D. Constrained deterministic synthesis of conformal supershaped antenna arrays. IEEE Trans. Antennas Propag. 2023, 71, 8726–8738. [Google Scholar] [CrossRef]
  30. Leung, K.W.; Huang, C.T.; Fang, Z.; Lu, K.; Yang, N. Wideband Circularly Polarized Horn Antenna with Twisted Ridges of Gielis SuperFormula Contour. IEEE Trans. Antennas Propag. 2024, 73, 1197–1202. [Google Scholar] [CrossRef]
  31. Rodrigo, J.A.; Alieva, T. Polymorphic beams and Nature inspired circuits for optical current. Sci. Rep. 2016, 6, 35341. [Google Scholar] [CrossRef] [PubMed]
  32. Xue, S.; Guan, J.; Yue, H.; Sun, Y.; Gao, Y.; Qi, X. Multidimensional controllable sunflower beams: turbulence-resistant encoding in free-space optical communication. J. Opt. Soc. Am. A 2026, 43, B14–B22. [Google Scholar] [CrossRef]
  33. Fougerolle, Y.D.; Gribok, A.; Foufou, S.; Truchetet, F.; Abidi, M.A. Boolean operations with implicit and parametric representation of primitives using R-functions. IEEE Trans. Vis. Comput. Graph. 2005, 11, 529–539. [Google Scholar] [CrossRef] [PubMed]
  34. Gielis, J.; Caratelli, D.; Moreno de Jong van Coevorden, C.; Ricci, P.E. The common descent of biological shape description and special functions. In Proceedings of the International Conference on Differential and Difference Equations and Applications; Springer International Publishing: Cham, Switzerland, 2017; pp. 119–131. [Google Scholar]
  35. Tavkhelidze, I.; Caratelli, D.; Gielis, J.; Ricci, P.E.; Rogava, M.; Transirico, M. On a geometric model of bodies with “Complex” configuration and some movements. In Modeling in Mathematics: Proceedings of the Second Tbilisi-Salerno Workshop on Modeling in Mathematics; Atlantis Press: Paris, France, 2017; pp. 129–158. [Google Scholar]
  36. Barr, A. H. Superquadrics and angle-preserving transformations. IEEE Computer graphics and Applications. 1981, 1(1), 11–23. [Google Scholar] [CrossRef]
  37. Jaklic, A.; Leonardis, A.; Solina, Franc. Segmentation and Recovery of Superquadrics of Computational Imaging and Vision; Springer, 2000; volume 20. [Google Scholar]
  38. Wang, L.; Ratkowsky, D. A.; Gielis, J.; Ricci, P. E.; Shi, P. Effects of the numerical values of the parameters in the Gielis equation on its geometries. Symmetry 2022, 14(12), 2475. [Google Scholar] [CrossRef]
Figure 1. (a)-(j) Superelliptic circle models for some values of r.
Figure 1. (a)-(j) Superelliptic circle models for some values of r.
Preprints 208301 g001
Figure 2. (a)-(j) Superelliptic helix models for some values of r.
Figure 2. (a)-(j) Superelliptic helix models for some values of r.
Preprints 208301 g002
Figure 3. Superelliptic curve (red) with Star derivative tangent vectors ( m = 8 , n 1 = n 2 = n 3 = 4 , a = b = 1 ) vs. standard derivative
Figure 3. Superelliptic curve (red) with Star derivative tangent vectors ( m = 8 , n 1 = n 2 = n 3 = 4 , a = b = 1 ) vs. standard derivative
Preprints 208301 g003
Figure 4. Moving frames along a helix α and a related superelliptic curve α S E , for some values of m = 12 , n 1 = 2 , n 2 = 7 , n 3 = 1 , a = b = 1 .
Figure 4. Moving frames along a helix α and a related superelliptic curve α S E , for some values of m = 12 , n 1 = 2 , n 2 = 7 , n 3 = 1 , a = b = 1 .
Preprints 208301 g004
Figure 5. Superellipsoid models for some values of superformula r.
Figure 5. Superellipsoid models for some values of superformula r.
Preprints 208301 g005
Figure 6. superelliptic cylinder surface for some values of superformula r.
Figure 6. superelliptic cylinder surface for some values of superformula r.
Preprints 208301 g006
Figure 7. Twisted superelliptic cylinder surface .
Figure 7. Twisted superelliptic cylinder surface .
Preprints 208301 g007
Figure 8. superelliptic cone surface for some values of superformula r.
Figure 8. superelliptic cone surface for some values of superformula r.
Preprints 208301 g008
Figure 9. superelliptic tori surface for some values of superformula r.
Figure 9. superelliptic tori surface for some values of superformula r.
Preprints 208301 g009
Figure 10. superelliptic catenoid surface for some values of superformula r.
Figure 10. superelliptic catenoid surface for some values of superformula r.
Preprints 208301 g010
Figure 11. superelliptic helicoid surfaces for some values of superformula r.
Figure 11. superelliptic helicoid surfaces for some values of superformula r.
Preprints 208301 g011
Figure 12. superelliptic Möbius strip for some values of superformula r.
Figure 12. superelliptic Möbius strip for some values of superformula r.
Preprints 208301 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated