Abstract: The definition of Gielis curves allows for the representation of self intersecting curves. The analysis and the understanding of these representations is of major interest for the analytical representation of sectors bounded by multiple subsets of curves (or surfaces), as this occurs for instance in many natural objects. We present a construction scheme based on R-functions to build signed potential fields with guaranteed differential properties, such that their zero-set corresponds to the outer, the inner envelop, or combined subparts of the curve. Our framework is designed to allow for the definition of composed domains built upon Boolean operations between several distinct objects or some subpart of self-intersecting curves, but also provides a representation for soft blending techniques in which the traditional Boolean union and intersection become special cases of linear combinations between the objects' potential fields. Finally, by establishing a connection between R-functions and Lame curves, we can extend the domain of the p parameter within the R-p-function from the set of the even positive numbers to the real numbers strictly greater than 1, i.e. p is an element of]1, +infinity[.

Abstract: Our scientific and technological worldviews are largely dominated by the concepts of entropy and complexity. Originating in 19th-century thermodynamics, the concept of entropy merged with information in the last century, leading to definitions of entropy and complexity by Kolmogorov, Shannon and others. In its simplest form, this worldview is an application of the normal rules of arithmetic. In this worldview, when tossing a coin, a million heads or tails in a row is theoretically possible, but impossible in practice and in real life. On this basis, the impossible (in the binary case, the outermost entries of Pascal's triangle xn and yn for large values of n) can be safely neglected, and one can concentrate fully on what is common and what conforms to the law of large numbers, in fields ranging from physics to sociology and everything in between. However, in recent decades it has been shown that what is most improbable tends to be the rule in nature. Indeed, if one combines the outermost entries xn and yn with the normal rules of arithmetic, either addition or multiplication, one obtains Lamé curves and power laws respectively. In this article, some of these correspondences are highlighted, leading to a double conclusion. First, Gabriel Lamé's geometric footprint in mathematics and the sciences is enormous. Second, conic sections are at the core once more. Whereas mathematics so far has been exclusively the language of patterns in the sciences, the door is opened for mathematics to also become the language of the individual. The probabilistic worldview and Lamé's footprint can be seen as dual methods. In this context, it is to be expected that the notions of information, complexity, simplicity and redundancy benefit from this different viewpoint.

Abstract: Electroporation is a non-thermal electromagnetic phenomenon widely used in medical diseases treatment. Different mathematical models of electroporation have been proposed in literature to study pore evolution in biological membranes. This paper presents a nonlinear dispersive multiphysic model of electroporation in irregular shaped biological cells in which the spatial and temporal evolution of the pores size is taken into account. The model solves Maxwell and asymptotic Smoluchowski equations and it describes the dielectric dispersion of cell media using a Debye-based relationship. Furthermore, the irregular cell shape has been modeled using the Gielis superformula. Taking into account the cell in mitosis phase, the electroporation process has been studied comparing the numerical results pertaining the model with variable pore radius with those in which the pore radius is supposed constant. The numerical analysis has been performed exposing the biological cell to a rectangular electric pulse having duration of 10 μs. The obtained numerical results highlight considerable differences between the two different models underling the need to include into the numerical algorithm the differential equation modeling the spatial and time evolution of the pores size.

Abstract: A new class of irregularly shaped dielectric lens antennas with a supershaped microstrip antenna feeder is presented and detailed in this work. The surface of the lens antenna and the feeder shape have been modelled by using the three and two-dimensional Gielis formula, respectively. The antenna design has been carried out by integrating an home-made software tool with the CST Microwave Studio®. The radiation properties of the whole antenna system have been evaluated using a dedicated high-frequency technique based on the tube tracing approximation. Moreover, the effects due to the multiple internal reflections have been properly modeled. The proposed model was applied to study unusual and complex lens antenna systems with the aim to design special radiation characteristics.

Abstract: The original motivation to study Generalized Möbius-Listing GML surfaces and bodies was the observation that the solution of boundary value problems greatly depends on the domains. Since around 2010 GML’s were merged with (continuous) Gielis Transformations, which provide a unifying description of geometrical shapes, as a generalization of the Pythagorean Theorem. The resulting geometrical objects can be used for modeling a wide range of natural shapes and phenomena. The cutting of GML bodies and surfaces, with the Möbius strip as one special case, is related to the field of knots and links, and classifications were obtained for GML with cross sectional symmetry of 2, 3, 4, 5 and 6. The general case of cutting GML bodies and surfaces, in particular the number of ways of cutting, could be solved by reducing the 3D problem to planar geometry. This also unveiled a range of connections with topology, combinatorics, elasticity theory and theoretical physics.

Abstract: Stomata are essential for the exchange of water vapour and atmospheric gases between vascular plants and their external environments. The stomatal geometries of many plants appear to be elliptical. However, prior studies have not tested whether this is a mathematical reality, particularly since many natural shapes that appear to be ellipses are superellipses with greater or smaller edge curvature than predicted for an ellipse. Compared with the ellipse equation, the superellipse equation includes an additional parameter that allows generation of a larger range of shapes. We randomly selected 240 stomata from each of four Magnoliaceae species to test whether the stomatal geometries are superellipses or ellipses. The stomatal geometries for most stomata (943/960) were found to be described better using the superellipse equation. The traditional “elliptical stomata hypothesis” resulted in an underestimation of the area of stomata, whereas the superellipse equation accurately predicted stomatal area. This finding has important implications for the estimation of stomatal area in studies looking at stomatal shape, geometry, and function.

Abstract: There is convincing evidence for a scaling relationship between leaf dry weight (DW) and leaf surface area (A) for broad-leaved plants, and most estimates of the scaling exponent of DW vs. A are greater than unity. However, the scaling relationship of leaf fresh weight (FW) vs. A has been largely neglected. In the present study, we examined whether there is a statistically strong scaling relationship between FW and A and compared the goodness of fit to that of DW vs. A. Between 250 and 520 leaves from each of 12 bamboo species within 2 genera (Phyllostachys and Pleioblastus) were investigated. The reduced major axis regression protocols were used to determine scaling relationships. The fit for the linearized scaling relationship of FW vs. A was compared with that of DW vs. A using the coefficient of determination (i.e., r2). A stronger scaling relationship between FW and A than that between DW and A was observed for each of the 12 bamboo species investigated. Among the 12 species examined, five had significantly smaller scaling exponents of FW vs. A compared to those of DW vs. A; only one species had a scaling exponent of FW vs. A greater than that of DW vs. A. No significant difference between the two scaling exponents was observed for the remaining 6 species. Researchers conducting future studies might be well advised to consider the influence of leaf fresh weight when exploring the scaling relationships of foliar biomass allocation patterns.

Abstract: Many geometries of plant organs can be described by the Gielis equation, a polar coordinate equation extended from the superellipse equation, . Here, r is the polar radius corresponding to the polar angle φ; m is a positive integer that determines the number of angles of the Gielis curve when φ ∈ [0 to 2π); and the rest of the symbols are parameters to be estimated. The pentagonal radial symmetry of calyxes and corolla tubes in top view is a common feature in the flowers of many eudicots. However, prior studies have not tested whether the Gielis equation can depict the shapes of corolla tubes. We sampled randomly 366 flowers of Vinca major L., among which 360 had five petals and pentagonal corolla tubes, and six had four petals and quadrangular corolla tubes. We extracted the planar coordinates of the outer rims of corolla tubes (in top view) (ORCTs), and then fitted the data with two simplified versions of the Gielis equation with k = 1 and m = 5: (Model 1), and (Model 2). The adjusted root mean square error (RMSEadj) was used to evaluate the goodness of fit of each model. In addition, to test whether ORCTs are radially symmetrical, we correlated the estimates of n2 and n3 in Model 1 on a log-log scale. The results validated the two simplified Gielis equations. The RMSEadj values for all corolla tubes were smaller than 0.05 for both models. The numerical values of n2 and n3 were demonstrated to be statistically equal based on the regression analysis, which suggested that the ORCTs of V. major are radially symmetrical. It suggests that Model 1 can be replaced by the simpler Model 2 for fitting the ORCT in this species. This work indicates that the pentagonal or quadrangular corolla tubes (in top view) can both be modeled by the Gielis equation and demonstrates that the pentagonal or quadrangular corolla tubes of plants tend to form radial symmetrical geometries during their development and growth.

Abstract: Measuring the inequality of leaf area distribution per plant (ILAD) can provide a useful tool for quantifying the influences of intra- and interspecific competition, foraging behavior of herbivores, and environmental stress on plants’ above-ground architectural structures and survival strategies. Despite its importance, there has been limited research on this issue. This paper aims to fill this gap by comparing four inequality indices to measure ILAD, using indices for quantifying household income that are commonly used in economics, including the Gini index (which is based on the Lorenz curve), the coefficient of variation, the Theil index, and the mean log deviation index. We measured the area of all leaves for 240 individual plants of the species Shibataea chinensis Nakai, a drought-tolerant landscape plant found in southern China. A three-parameter performance equation was fitted to observations of the cumulative proportion of leaf area vs. the cumulative proportion of leaves per plant to calculate the Gini index for each individual specimen of S. chinensis. The performance equation was demonstrated to be valid in describing the rotated and right shifted Lorenz curve, given that >96% of root-mean-square error values were smaller than 0.004 for 240 individual plants. By examining the correlation between any of the six possible pairs of indices among the Gini index, the coefficient of variation, the Theil index, and the mean log deviation index, the data show that these indices are closely related and can be used interchangeably to quantify ILAD.

Abstract: The Robin problem for the Helmholtz equation in normal-polar shells is addressed by using a suitable spherical harmonic expansion technique. Attention is in particular focused on the wide class of domains whose boundaries are defined by a generalized version of the so-called superformula introduced by Gielis. A dedicated numerical procedure based on the computer algebra system Mathematica? is developed in order to validate the proposed methodology. In this way, highly accurate approximations of the solution, featuring properties similar to the classical ones, are obtained.

Abstract: Many natural shapes exhibit surprising symmetry and can be described by the Gielis equation, which has several classical geometric equations (for example, the circle, ellipse and superellipse) as special cases. However, the original Gielis equation cannot reflect some diverse shapes due to limitations of its power-law hypothesis. In the present study, we propose a generalized version by introducing a link function. Thus, the original Gielis equation can be deemed to be a special case of the generalized Gielis equation (GGE) with a power-law link function. The link function can be based on the morphological features of different objects so that the GGE is more flexible in fitting the data of the shape than its original version. The GGE is shown to be valid in depicting the shapes of some starfish and plant leaves.

Abstract: Preston's equation is a general model describing the egg shape of birds. The parameters of Preston's equation are usually estimated after re-expressing it as the Todd-Smart equation and scaling the egg's actual length to two. This method assumes that the straight line through the two points on an egg's profile separated by the maximum distance (i.e., the longest axis of an egg's profile) is the mid-line. It hypothesizes that the photographed egg's profile is perfectly bilaterally symmetrical, which seldom holds true because of photographic errors and placement errors. The existing parameter estimation method for Preston's equation considers an angle of deviation for the longest axis of an egg's profile from the mid-line, which decreases prediction errors to a certain degree. Nevertheless, this method cannot provide an accurate estimate of the coordinates of the egg's center, and it leads to sub-optimal parameter estimation. Thus, it is better to account for the possible asymmetry between the two sides of an egg's profile along its mid-line when fitting egg-shape data. In this paper, we propose a method based on the optimization algorithm (optimPE) to fit egg-shape data and better estimate the parameters of Preston's equation by automatically searching for the optimal mid-line of an egg's profile and testing its validity using profiles of 59 bird eggs spanning a wide range of existing egg shapes. We further compared this method with the existing one based on multiple linear regression (lmPE). This study demonstrated the ability of the optimPE method to estimate numerical values of the parameters of Preston's equation and provide the theoretical egg length (i.e., the distance between two ends of the mid-line of an egg's profile) and the egg's maximum breadth. This provides a valuable approach for comparing egg shapes among conspecifics or across different species, or even different classes (e.g., birds and reptiles), in future investigations.

Abstract: Many cross-sectional shapes of plants have been found to approximate a superellipse rather than an ellipse. Square bamboos, belonging to the genus Chimonobambusa (Poaceae), are a group of plants with round-edged square-like culm cross sections. The initial application of superellipses to model these culm cross sections has focused on Chimonobambusa quadrangularis (Franceschi) Makino. However, there is a need for large scale empirical data to confirm this hypothesis. In this study, approximately 750 cross sections from 30 culms of C. utilis were scanned to obtain cross-sectional boundary coordinates. A superellipse exhibits a centrosymmetry, but in nature the cross sections of culms usually deviate from a standard circle, ellipse, or superellipse because of the influences of the environment and terrain, resulting in different bending and torsion forces during growth. Thus, more natural cross-sectional shapes appear to have the form of a deformed superellipse. The superellipse equation with a deformation parameter (SEDP) was used to fit boundary data. We find that the cross-sectional shapes (including outer and inner rings) of C. utilis can be well described by SEDP. The adjusted root-mean-square error of SEDP is smaller than that of the superellipse equation without a deformation parameter. A major finding is that the cross-sectional shapes can be divided into two types of superellipse curves: hyperellipses and hypoellipses, even for cross sections from the same culm. There are two proportional relationships between ring area and the product of ring length and width for both the outer and inner rings. The proportionality coefficients are significantly different, as a consequence of the two different superellipse types (i.e., hyperellipses and hypoellipses). The difference in the proportionality coefficients between hyperellipses and hypoellipses for outer rings is greater than that for inner rings. This work informs our understanding and quantifying of the longitudinal deformation of plant stems for future studies to assess the influences of the environment on stem development. This work is also informative for understanding the deviation of natural shapes from a strict rotational symmetry.

Abstract: The size and shape of plant leaves change with growth, and an accurate description of leaf shape is crucial for describing plant morphogenesis and development. Bilateral symmetry, which has been widely observed but poorly examined, occurs in both dicot and monocot leaves, including all nominated bamboo species (approximately 1,300 species), of which at least 500 are found in China. Although there are apparent differences in leaf size among bamboo species due to genetic and environmental profiles, bamboo leaves have bilateral symmetry with parallel venation and appear similar across species. Here, we investigate whether the shape of bamboo leaves can be accurately described by a simplified Gielis equation, which consists of only two parameters (leaf length and shape) and produces a perfect bilateral shape. To test the applicability of this equation and the occurrence of bilateral symmetry, we first measured the leaf length of 42 bamboo species, examining >500 leaves per species. We then scanned 30 leaves per species that had approximately the same length as the median leaf length for that species. The leaf-shape data from scanned profiles were fitted to the simplified Gielis equation. Results confirmed that the equation fits the leaf-shape data extremely well, with the coefficients of determination being 0.995 on average. We further demonstrated the bilateral symmetry of bamboo leaves, with a clearly defined leaf-shape parameter of all 42 bamboo species investigated ranging from 0.02 to 0.1. This results in a simple and reliable tool for precise determination of bamboo species, with applications in forestry, ecology, and taxonomy.

Abstract: Plant leaves exhibit diverse shapes that enable them to utilize a light resource maximally. If there were a general parametric model that could be used to calculate leaf area for different leaf shapes, it would help to elucidate the adaptive evolutional link among plants with the same or similar leaf shapes. We propose a simplified version of the original Gielis equation (SGE), which was developed to describe a variety of object shapes ranging from a droplet to an arbitrary polygon. We used this equation to fit the leaf profiles of 53 species (among which, 48 bamboo plants, 5 woody plants, and 10 geographical populations of a woody plant), totaling 3310 leaves. A third parameter (namely, the floating ratio c in leaf length) was introduced to account for the case when the theoretical leaf length deviates from the observed leaf length. For most datasets, the estimates of c were greater than zero but less than 10%, indicating that the leaf length predicted by the SGE was usually smaller than the actual length. However, the predicted leaf areas approximated their actual values after considering the floating ratios in leaf length. For most datasets, the mean percent errors of leaf areas were lower than 6%, except for a pooled dataset with 42 bamboo species. For the elliptical, lanceolate, linear, obovate, and ovate shapes, although the SGE did not fit the leaf edge perfectly, after adjusting the parameter c, there were small deviations of the predicted leaf areas from the actual values. This illustrates that leaves with different shapes might have similar functional features for photosynthesis, since the leaf areas can be described by the same equation. The anisotropy expressed as a difference in leaf shape for some plants might be an adaptive response to enable them to adapt to different habitats.