Preprint
Article

This version is not peer-reviewed.

Descartes Curvature Geometry, Compact-Cycle Amplitudes, and the GST Cabibbo Seed

Submitted:

20 August 2026

Posted:

21 August 2026

You are already at the latest version

Abstract
Several apparently distinct formulations of the charged-lepton mass relation \( p_2/e_1^2=2/3 \) are shown to describe the same algebraic structure. The compact-cycle Fourier-amplitude condition, the Descartes curvature condition, and the smaller Soddy completion are related exactly. After normalization by the total mass scale, the smaller completion is universal on the positive exact \( 2/3 \) locus and selects \( \alpha=\sqrt{3/2}-1 \). Its square \( r=\alpha^2 \) therefore defines the associated homogeneous map \( \mathcal M(x)=rx \), whose iterates form an exact geometric-mean orbit. The lepton pole masses and this universal ratio fix \( \mu_\star \), \( x_1 \), and \( x_2 \) without quark-sector input. The first iterate is empirically identified with the running strange-quark mass; once that correspondence is made, the next iterate is fixed. Thus the exact orbit ratio \( x_2/x_1=\alpha^2 \) becomes the conditional quark statement \( m_d/m_s=\alpha^2 \) and the leading GST/Fritzsch factor \( \sqrt{m_d/m_s}=\alpha \). A separate terminal boundary hypothesis gives a target for \( m_u \). Explicit PDG inputs and four-loop running make the numerical comparison auditable. The construction fixes a leading Cabibbo seed, not a complete CKM prediction.
Keywords: 
;  ;  ;  ;  

1. Introduction

The Cabibbo angle has been measured to high precision through the CKM element | V u s | , with the current world average [1] | V u s | = 0.22503 ± 0.00068 . The Gatto–Sartori–Tonin (GST) relation [2,3]
sin θ C m d m s ,
connects this angle to the ratio of down- and strange-quark masses. Gatto, Sartori, and Tonin originally derived the Cabibbo angle from the cancellation of quadratically divergent weak self-masses; [2] the quark-mass-ratio form (1) is due to the Fritzsch texture-zero ansatz. [3] Equation (1) is well established; what remains unexplained, within this phenomenological relation, is why the ratio m d / m s has its observed value.
Separately, the empirical relation
m m 2 = 2 3
has held for the charged-lepton pole masses to one part in 10 5 since it was first noted by Koide. [4] Several authors have recast this relation in different mathematical languages, [6,7,8] raising the question of whether these are independent observations or coordinate descriptions of a single underlying structure.
A useful feature of the construction is that the quantities later compared with quark masses are fixed without quark-sector input. The scale μ = m pole is constructed solely from measured charged-lepton pole masses, while the dimensionless factor r = α 2 follows from the normalized charged-lepton geometry. Consequently the orbit targets x 1 = r μ and x 2 = r 2 μ are determined before any quark mass, quark renormalization convention, or lattice result enters. The quark scheme and scale dependence therefore resides on the comparison side of the test rather than in the construction of the targets.
In this paper we pursue two aims. First, we show that the Fourier-mode constraint of Shulga’s compact-cycle model, [7] the Descartes curvature reading of Kocik, [6] and the Soddy completion of the resulting circle configuration are related by exact algebra, with no approximations (Sec. Section 3). We then isolate the normalized content of that completion and show that it supplies a universal ratio and an associated homogeneous map whose iterates form an exact geometric-mean orbit (Sec. Section 4). Second, beginning from the empirical strange-quark correspondence, we compare the fixed next iterate with the down-quark mass and examine the resulting conditional GST/Fritzsch Cabibbo seed (Secs. Section 5Section 6).

2. The Algebraic Object

Let e 1 = i m i , p 2 = i m i for the three charged leptons, and define the lepton-sum scale μ p 2 = 1883.1 MeV. The condition (2) is p 2 = 2 3 e 1 2 , equivalently e 1 = 3 / 2 p 2 .
Define the companion amplitude
F e 1 p 2 = α μ ,
where
α 3 2 1 0.2247 .
The squared companion is
F 2 = α 2 μ , α 2 = 5 2 6 0.0505 .
This is the central algebraic object. The three-input symmetric-polynomial identity (2) collapses to one dimensionless constant α times one mass scale μ , with α fixed by the condition (2).

3. Equivalence of Descriptions

We now show that three apparently distinct formulations of the condition (2) describe the same constraint.

3.1. Compact-Cycle Amplitudes (Shulga)

Shulga [7] constructs a real amplitude on an internal circle,
Z ( ϕ ) = c 0 + 2 r c cos ( ϕ + δ ) ,
built from the two lowest antiperiodic modes, with the three charged-lepton families sampling at equally spaced points: m a = Z ( 2 π a / 3 ) for a = 0 , 1 , 2 . The parameter c 0 is the flavor-symmetric (democratic) component and r c the flavor-breaking amplitude; the cosine coefficient is 2 r c , so Shulga’s notation Z = s [ 1 + 2 cos ( ϕ + θ ) ] corresponds to c 0 = s and 2 r c / c 0 = 2 , i.e. r c / c 0 = 1 / 2 .
By the orthogonality of equally spaced phases on the circle, a cos ( 2 π a / 3 + δ ) = 0 for any δ , and similarly for the double-angle sum. Therefore
e 1 = a Z ( 2 π a / 3 ) = 3 c 0 ,
p 2 = a Z ( 2 π a / 3 ) 2 = 3 ( c 0 2 + 2 r c 2 ) .
The ratio p 2 / e 1 2 becomes
p 2 e 1 2 = c 0 2 + 2 r c 2 3 c 0 2 = 1 3 + 2 r c 2 3 c 0 2 .
Setting this equal to 2 / 3 gives
r c c 0 = 1 2 .
The condition (2) is exactly the statement that the flavor-breaking and flavor-symmetric amplitudes on the compact cycle stand in the ratio 1 : 2 . Shulga derives this ratio from the spinor-squaring of antiperiodic modes and a reality condition on the amplitude. [7]

3.2. Descartes Curvatures (Kocik)

The Descartes circle theorem [5] for four mutually tangent circles reads ( i k i ) 2 = 2 i k i 2 . Kocik [6] observed that identifying curvatures k i = m i places (2) within the generalized Descartes family. The present work uses the standard (mutually tangent) theorem and computes the fourth curvature.
The outer Soddy curvature of the three-circle configuration is k 4 = e 1 2 e 2 , where e 2 = i < j m i m j . Under (2), one has 4 e 2 = p 2 , so k 4 = e 1 p 2 = F : the Soddy curvature is the companion amplitude (3).

3.3. The Equivalence

Expressing the Soddy curvature in Shulga’s variables via (7)–(10):
F = e 1 p 2 = 3 c 0 6 c 0 = c 0 ( 3 6 ) .
Then α = F / μ = c 0 ( 3 6 ) / 6 c 0 2 = ( 3 6 ) / 6 = 3 / 2 1 , confirming (4).
The three descriptions are related by exact algebra:
Framework Constraint Language
Shulga r c / c 0 = 1 / 2 Fourier modes
Condition (2) p 2 / e 1 2 = 2 / 3 Sym. polynomials
Soddy F = α μ Inversive geometry
No approximations are involved in any direction. The constant α = 3 / 2 1 appears in each language as the unique dimensionless number implied by the shared constraint. The plus branch e 1 + p 2 gives 3 / 2 + 1 2.22 , which has no interpretation in the quark sector; the minus branch is the only phenomenologically relevant root.

4. Universal Normalized Completion and the Companion Map

The value of the smaller completion is not tied to the particular orientation of the physical charged-lepton triple. For any positive triple satisfying the exact condition (2), define
y i m i μ , μ i m i .
Then i y i 2 = 1 and i y i = 3 / 2 . Normalizing the fourth Descartes curvature by the same scale,
f k 4 μ ,
the Descartes equation becomes
3 2 + f 2 = 2 ( 1 + f 2 ) ,
with the two roots
f ± = 3 2 ± 1 .
Thus the normalized smaller completion is
f = α = 3 2 1 ,
independent of the orientation of the normalized square-root-mass vector on the positive exact- 2 / 3 locus. This strengthens Eq. (5): the quantity
r α 2 = 5 2 6
is a universal squared-completion ratio rather than a number peculiar to one charged-lepton configuration.
We therefore define the associated homogeneous companion map
M : R > 0 R > 0 , M ( x ) = r x .
This is a mathematical definition, not a dynamical evolution law, and by itself contains no quark-flavor assignment. Its iterates satisfy
x n = M n ( x 0 ) = r n x 0 , x n 2 = x n 1 x n + 1 .
Hence every interior orbit point is exactly the geometric mean of its two neighbors. The geometric progression used below is therefore the orbit of one fixed map; no additional step ratio is introduced when the map is iterated.

5. The Mass Cascade

For numerical work we use the charged-lepton pole masses quoted by PDG, [1]
m e = 0.5109989507 MeV , m μ = 105.6583755 MeV , m τ = 1776.93 ( 9 ) MeV ,
which give
μ = 1883.0994 ± 0.090 MeV .
Here μ is a physical energy constructed from measured lepton pole masses. The distinction is useful rather than problematic: the PDG quark-mass review emphasizes that observable leptons admit pole-mass definitions, whereas confined quark masses require an explicitly specified renormalization scheme and scale. [1] In the present construction no quark quantity enters the definition of μ , r, x 1 , or x 2 . The quark masses are only then evaluated, independently, in the MS ¯ scheme at the numerical energy μ R = μ . Thus no conversion of a lepton pole mass into a quark running mass is implied, and the fixed targets cannot be an artifact of the quark renormalization convention used to test them. The phenomenological hypothesis is instead that this externally specified lepton-defined physical energy is relevant to the light-quark spectrum.
Taking x 0 = μ , the first two nontrivial orbit points of Eq. (19) are
x 1 = M ( μ ) = α 2 μ = 95.11583 MeV ,
x 2 = M 2 ( μ ) = α 4 μ = 4.804325 MeV .
The empirical strange-quark correspondence predates the present manuscript and is reported in Ref. [9]. The fixed light-quark targets used here were subsequently deposited in a dated pre-registration on 18 May 2026, Ref. [10], before the 2026 PDG edition dated 1 June 2026. That record fixes the s, d, and u targets at 95.12, 4.804, and 2.216 MeV, respectively, with their stated mass conventions and scales. The exact-locus value x 1 = 95.11583 MeV used here differs from the archived physical-point value only by the deterministic exact- 2 / 3 idealization quantified below, not by refitting. When confronted with the 2026 PDG input used below, the strange central value lies farther from the fixed target but remains compatible at current precision; no parameter, scale, or target was adjusted. Because the earlier comparison used 2024 FLAG/PDG inputs whereas Table 1 uses the 2026 PDG final lattice estimate, the change in normalized residual is not interpreted as a time-series trend. The prospective point is that the numerical targets were frozen before this subsequent data cycle, although the original strange-quark identification was retrospective.
We retain that empirical identification,
x 1 m s MS ¯ ( μ ) ,
and then apply the already defined map once more. The down-quark target x 2 is therefore fixed by the same universal ratio; no second step size is chosen. The resulting phenomenological comparison is
x 2 m d MS ¯ ( μ ) .
The charged-lepton geometry does not independently derive the flavor labels in Eqs. (24)–(25); the claim is that the fixed continuation from the observed strange correspondence lands on the running down-quark mass.
The u-quark closure has a different logical status. Retaining the physical pair threshold T = 2 m e as a lower boundary, we impose the separate midpoint hypothesis
m u MS ¯ ( μ ) 2 = hyp T m d model ( μ ) , T = 2 m e pole ,
which gives
m u model ( μ ) = 2 m e x 2 = 2.21585 MeV .
This terminal condition is not a third iterate of M and is not derived from the normalized-completion result.

5.1. Sensitivity to the Exact 2 / 3 Idealization

The algebra above is exact on the 2 / 3 locus, while the measured charged-lepton masses are slightly off it. For general Q = p 2 / e 1 2 , Eq. (3) gives the normalized smaller completion
f ( Q ) = Q 1 / 2 2 Q 1 1 .
Using the central values in Eq. (20) gives
Q , phys = 0.66666446 , f ( Q , phys ) = 0.22474194 , r phys = 0.05050894 .
Replacing the exact r by this physical-point value shifts the first two orbit targets by only
Δ x 1 = 0.00248 MeV , Δ x 2 = 0.000251 MeV ,
and changes the terminal u target by about 5.8 × 10 5 MeV. Propagating the present m τ uncertainty through the off-locus completion gives about 0.0103 MeV and 0.00081 MeV on the corresponding x 1 and x 2 values. These effects are small compared with the current quark-mass uncertainties, so the exact- 2 / 3 idealization is kept in the closed-form relations and its measured-point shift is treated as a separate sensitivity check.

5.2. Auditable Quark-Mass Comparison

For a single reproducible data prescription, we use the PDG final lattice-QCD estimates in the four-flavor MS ¯ theory at 2 GeV, [1]
m s = 92.74 ± 0.54 MeV , m d = 4.69 ± 0.05 MeV , m u = 2.20 ± 0.07 MeV .
All three are evolved to μ R = μ = 1.883099 GeV with the same four-loop pure-QCD running at fixed n f = 4 , using RunDec/CRunDec. [11] No heavy-flavor threshold is crossed in this short interval. With α s ( M Z ) = 0.1180 , the common mass-evolution factor is
R m m q ( μ ) m q ( 2 GeV ) = 1.01744 .
The model-target uncertainties from Eq. (21) are approximately 0.00455 MeV for x 1 , 0.000230 MeV for x 2 , and 5.3 × 10 5 MeV for the terminal u target.
Using the published individual mass uncertainties, the three absolute central-value differences correspond to approximately 1.38 , 0.64 , and 0.32 input standard deviations for s, d, and u, respectively. These are one-observable residuals; they are not treated as statistically independent pieces of evidence for a common model.
The short running is not numerically delicate. At the same local coupling, the three-, four-, and five-loop factors are 1.017156 , 1.017440 , and 1.017530 . A deliberately broad variation α s ( 4 ) ( 2 GeV ) = 0.30082 ± 0.010 changes the evolved m s by at most 0.081 MeV and m d by at most 0.0041 MeV; adjacent loop-order variation adds at most 0.027 MeV and 0.0014 MeV, respectively. Leading QED running over the same short interval is about 2.3 × 10 5 fractionally for d and s and 9.3 × 10 5 for u. [12] Thus the absolute s , d comparison is presently dominated by the quoted quark inputs rather than the short RG/QED translation. For the terminal u relation, the different electric charge means that QED/isospin separation does not cancel as cleanly as it does for the d / s ratio, so the u comparison is kept at the descriptive percent level.

5.3. Scale Sensitivity of the Orbit Comparison

The comparison point μ is fixed by the charged-lepton pole-mass sum; it is not obtained by fitting the quark masses. To quantify how strongly the numerical agreement depends on that specified scale, we solve for the scales at which the same four-loop running masses cross the fixed orbit targets:
m s MS ¯ ( μ s ) = x 1 μ s = 1832.8 MeV , m d MS ¯ ( μ d ) = x 2 μ d = 1840.3 MeV .
These central crossings lie 2.67 % and 2.27 % below μ = 1883.099 MeV. Propagating the individual PDG mass uncertainties through the same running gives crossing-scale uncertainties of approximately ± 36 MeV for s and ± 66 MeV for d. The broader RG/QED prescription stress envelope quoted above corresponds to only about ± 7 MeV in either crossing location and is therefore subdominant to the present quark-input uncertainties. Equation (33) is a scale-sensitivity diagnostic, not an additional fit or independent validation: present precision does not sharply localize a preferred renormalization scale, while the independently specified μ lies within the few-percent region over which the correspondence remains numerically close.

6. The GST Cabibbo Seed

From Eqs. (22)–(), the exact mathematical result is
x 2 x 1 = r = α 2 = 5 2 6 .
Under the empirical ordered correspondence in Eqs. (24)–(25), this fixed orbit ratio gives the conditional quark-sector statement
m d m s = α 2 .
The leading GST relation (1) then yields
m d m s = α = 3 2 1 = 0.22474
This is the down-sector Cabibbo seed in the sense of the GST/Fritzsch leading-order relation. It is not, by itself, a precision prediction of | V u s | . The full GST expression contains the up-sector contribution and a relative phase, [3]
| V u s | m d m s e i ϕ m u m c ,
and neither the corresponding charm-sector input nor the phase ϕ is fixed by the present construction. The result in Eq. (36) is therefore the exact leading down-sector factor conditional on the s , d orbit correspondence; a complete CKM test remains open.

7. The Fritzsch Texture

The GST relation arises from the four-zero texture [3]
M d = 0 A 0 A * 0 B 0 B * C .
In the present construction the Fritzsch texture is interpretive rather than part of the Descartes proof. Conditional on the orbit assignments, its leading down-sector mass-ratio factor is fixed by the same α = 3 / 2 1 selected by the normalized completion.
The two map-generated quark targets and the separate lower boundary can be summarized as
μ M m s M m d , m u = m d · 2 m e .
The fixed multiplicative step of the orbit is r = α 2 —the square of the leading GST Cabibbo seed.

8. Weak-Basis Considerations

The condition (2) holds for charged-lepton pole masses, where the mass and weak bases are approximately aligned: in the Standard Model, the charged-lepton Yukawa matrix can be taken diagonal without loss of generality when neutrino masses are neglected. The amplitude structure described in Sec. Section 3 is therefore visible directly in the physical spectrum.
For quarks the situation is different. The CKM matrix rotates between up-type and down-type mass bases, so the condition (2) should not be expected to hold for quark pole masses directly. The orbit comparisons (22)–() and the terminal condition (26) use MS ¯ running masses at a specific scale, and the GST relation (1) itself arises from the Fritzsch texture in the weak basis. Whether the amplitude structure of Sec. Section 3 is native to the weak-basis Yukawa matrix, or requires a more involved embedding, is an open question.
The cross-sector comparison is therefore scale-specified rather than a claim that pole and running masses are the same kind of parameter. The lepton pole masses define the physical comparison energy μ and the dimensionless constant α ; the quark quantities are all evaluated in one MS ¯ prescription at that numerical scale. This use of a lepton-defined energy is not excluded by the sector labels themselves: in the Standard Model the charged-lepton and quark Yukawa matrices are all couplings to the same Higgs field. [1] That common electroweak origin does not derive the present relation, but it provides a physical setting in which an additional cross-sector flavor structure can be meaningfully asked about. A dynamical mechanism explaining why the same constant and characteristic energy should recur across the two sectors remains open.

9. Tests and Failure Modes

The exact mathematics in Secs. Section 3Section 4 is distinct from the quark interpretation. The latter can fail in concrete ways. Future quark-mass determinations may move m s MS ¯ ( μ ) or m d MS ¯ ( μ ) away from the fixed orbit targets in Eqs. (22)–(); the scale, ratio, ordering, and flavor assignments stated here are not to be retuned when those inputs are updated. The separate terminal condition (26) can likewise fail as m u determinations and QCD–QED separation improve.
A completed GST/Fritzsch treatment can also fail if an up-sector texture and relative phase fixed without fitting | V u s | do not reproduce CKM mixing. Finally, the cross-sector interpretation would be weakened if a consistent higher-precision QCD+QED treatment moved the quark quantities away from the present correspondence, or if no dynamical embedding of the compact-cycle structure into the quark Yukawa sector can be found. The strange target had already been publicly fixed before the 2026 PDG update, which moved its central value away from the target without prompting any retuning; future updates continue that prospective comparison. Publication of the present manuscript additionally freezes the d-orbit assignment and the separate terminal u boundary for subsequent data cycles.

11. Summary

The charged-lepton mass condition p 2 / e 1 2 = 2 / 3 , the compact-cycle amplitude constraint r c / c 0 = 1 / 2 , and the smaller Descartes–Soddy completion are exact algebraic descriptions of the same positive charged-lepton constraint. Normalizing the completion reveals the additional result that f = α = 3 / 2 1 is universal over the positive exact- 2 / 3 locus. Its square therefore defines one scale-independent ratio r = α 2 and the homogeneous map M ( x ) = r x , whose iterates obey the exact geometric-mean identity x n 2 = x n 1 x n + 1 .
The physical comparison begins after that mathematical result. Importantly, μ , r, x 1 , and x 2 are fixed entirely from charged-lepton observables and the normalized completion, with no quark input; the renormalization convention enters only when the independently defined quark masses are evaluated against those targets. The empirically identified first lepton-anchored orbit point lies near the running strange-quark mass; once the same fixed map is applied again, the next target is fixed and lies near the running down-quark mass. With the 2026 PDG inputs evolved to the common scale μ , the central differences are 0.80 % for s and 0.68 % for d. These are two components of one ordered orbit correspondence rather than independent evidential claims. The strange target had been publicly posted before the 2026 PDG update; that update moved the central value away from the target without any retuning. The exact central crossing scales are 1832.8 and 1840.3 MeV for s and d, respectively, showing that the fixed μ comparison is stable at the present few-percent level in renormalization scale rather than sharply localized there. A separate lower-boundary hypothesis gives an m u target within about 1.0 % of the evolved PDG central value; it is an auxiliary terminal conjecture rather than evidence for the s , d orbit. The measured departure of the charged-lepton masses from the exact 2 / 3 locus and the short RG/QED translation are both small relative to the present quark uncertainties.
The exact orbit statement is x 2 / x 1 = α 2 . Conditional on the observed s , d correspondence, this becomes m d / m s = α 2 , and the leading GST/Fritzsch down-sector factor is exactly m d / m s = α = 0.224744871 . This is a leading texture-level Cabibbo seed, not a complete CKM prediction. What remains open is the physical origin of the cross-sector correspondence and a dynamical embedding that fixes the up-sector structure and phase without additional fitted freedom.
AI Usage Disclosure: Generative AI tools were used solely for programming, writing assistance, and adversarial review. All scientific ideas, claims, and conclusions presented in this paper originate with and are approved by the author.

Author Contributions

Conceptualization, formal analysis, software, investigation, visualization, writing—original draft preparation, and writing—review and editing: A.M.B.

Funding

This research received no external funding.

Data Availability Statement

All numerical results reported here are derived from published sources cited in the text and are reproducible from the inputs and running prescriptions stated in Section 5 and Section 5.3. The pre-registered cascade targets are archived at Zenodo [10].

Acknowledgments

The author thanks A. Rivero for extensive correspondence on the algebraic connections between lepton and quark sectors, the lattice QCD community for successive improvements to light-quark mass determinations, and K. Shulga for the compact-cycle construction that clarifies the amplitude-level content of the lepton mass relation.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Takahashi, F.; et al. [Particle Data Group], Review of Particle Physics. Int. J. Mod. Phys. A 2026, 41, 2630011. [Google Scholar]
  2. Gatto, R.; Sartori, G.; Tonin, M. Phys. Lett. B 1968, 28, 128.
  3. Fritzsch, H. Nucl. Phys. B 1979, 155, 189. [CrossRef]
  4. Koide, Y.; Lett. Nuovo Cim. 1982, 34, 201. [CrossRef]
  5. Coxeter, H. S. M. Amer. Math. Mon. 1968, 75, 5.
  6. Kocik, J. Koide Lepton Mass Formula Geom. Circ. arXiv:1201.2067. [CrossRef]
  7. Shulga, K. Charg.-Lepton Koide Geom. From A Green-Dressed Compact Fam. Cycle, arXiv:2605.10245. [CrossRef]
  8. Żenczykowski, P. Phys. Rev. D. 2012, 86, 117303. [CrossRef]
  9. Brilliant, A. M. A Monte Carlo Null-Model Test of an Outer-Soddy Completion of the Koide Lepton Triple . SSRN. [CrossRef]
  10. Brilliant, A. M. Pre-Registered Fermion Mass Predictions for Comparison with PDG 2026 . Zenodo, 18 May 2026. [CrossRef]
  11. Herren, F.; Steinhauser, M. Comput. Phys. Commun. 2018, 224, 333. [CrossRef]
  12. Tarasov, O. V. Phys. Part. Nucl. Lett. 2020, arXiv:1910.1223117, 109.
  13. Rivero, A. A New Koide Tuple Strange-Charm-Bottom arXiv:1111.7232. [CrossRef]
Table 1. Light-quark comparison using one MS ¯ prescription. All masses are in MeV.
Table 1. Light-quark comparison using one MS ¯ prescription. All masses are in MeV.
Model at μ PDG at 2 GeV PDG at μ Central diff.
s 95.11583 92.74 ± 0.54 94.35739 ± 0.54942 + 0.804 %
d 4.804325 4.69 ± 0.05 4.771794 ± 0.050872 + 0.682 %
u 2.21585 2.20 ± 0.07 2.238368 ± 0.071221 1.006 %
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.