OPEN ACCESS Modern Physics Letters A (2026) 2650073 (16 pages) . The Author(s) #c DOI: 10.1142/S0217732326500732 Uni¯ed fermion mass ratios from orbifold conformal ¯eld theory Mod. Phys. Lett. A Downloaded from www.worldscientific.com Samir Varma VS Asset Management, LLC 16 Windsor Lane, Cos Cob, CT 06807, USA samir@vsasset.com Received 4 November 2025 Revised 2 January 2026 Accepted 28 January 2026 Published 9 March 2026 Motivated by simple-current orbifolds, where twisted sectors come with characteristic chiral and left{right \degrees," we construct the minimal modular-covariant, permutation-symmetric, P 1=k k q ð mi Þ degree-zero functional of charged-fermion masses, R ¼ P i 1=ð2kÞ 2k q : This CFT-motivated ð i mi Þ ansatz reproducesR‘ ¼ 2=3 for charged leptons (the known Koide relation) and yields a by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. quark-sector ratio numerically consistent with Rq ¼ 1=2 when evaluated on current MS quark masses. The ratio is approximately invariant at one loop; full two-loop running causes only few- 103 -level variations in Rq over the range MZ {1014 GeV. Under SM running, this predicts Rq ð1014 GeVÞ ¼ 0:500  0:002; a precision target within reach of forthcoming lattice stepscaling calculations. The construction is numerically veri¯ed, RG-stable, and falsi¯able through improved mass determinations. This framework o®ers a CFT-motivated route to linking lowenergy °avor phenomenology and high-scale orbifold dynamics. Keywords: Orbifold conformal ¯eld theory; fermion mass relations; modular invariance; °avor physics; Koide relation; renormalization group invariance. PACS Nos.: 11.25.Mj, 12.15.Ff, 14.65.q We construct a minimal modular-covariant, permutation-symmetric mass ratio from orbifold conformal ¯eld theory that uni¯es the Koide relation for charged leptons (R‘ ¼ 2=3) with a parameter-free functional for quarks whose numerical value is consistent with Rq ¼ 1=2. The ratio is approximately renormalization-group invariant at one loop, with two-loop running shifting Rq by only  103 between MZ and 1014 GeV. This yields a precision target Rq ð1014 GeVÞ ¼ 0:500  0:002 testable by forthcoming lattice step-scaling calculations. This is an Open Access article published by World Scienti¯c Publishing Company. It is distributed under the terms of the Creative Commons Attribution 4.0 (CC BY) License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. 2650073-1 S. Varma The pattern of fermion masses remains one of the most persistent puzzles in particle physics. While the Standard Model successfully describes particle interactions, it provides no explanation for the observed mass hierarchies spanning ¯ve orders of magnitude. Among various empirical observations, the Koide relation1 for charged leptons stands out for its remarkable precision and stability: me þ m þ m 2 K  pffiffiffiffiffiffi pffiffiffiffiffiffi pffiffiffiffiffiffi 2  : ð1Þ ð me þ m þ m Þ 3 Mod. Phys. Lett. A Downloaded from www.worldscientific.com Despite signi¯cant e®orts including symmetry-based approaches,2,3 modular symmetries,4 and °avor models,5 a uni¯ed theoretical framework explaining both lepton and quark mass relations has remained elusive. Unlike earlier studies that searched data for appealing fractions and then built models around them, our ratio is motivated by orbifold-CFT consistency and contains no free parameters. Attempts to generalize Koide-like relations to quarks have faced challenges due to their wider mass hierarchies and °avor dynamics. Earlier attempts to extend Koidetype ratios to quarks    e.g. texture-zero ansätze and modular-symmetry ¯ts   have so far required adjustable phases or family-speci¯c parameters. The study of fermion mass relations has a rich history spanning several decades. E®orts to understand particle mass patterns date back to Gell-Mann's work on hadron mass relations.6 As the quark model developed, attention shifted to ¯nding patterns among fundamental fermion masses. The Georgi{Jarlskog relations2 proby 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. posed a framework connecting quark and lepton masses through GUT-scale physics, while the Fritzsch ansatz3,7 suggested speci¯c texture zeros in the mass matrices to explain both masses and mixing angles. Weinberg8 highlighted early on the challenge of maintaining mass relations under renormalization group evolution, an issue that our present work addresses through approximate one-loop RG invariance. Two-loop terms shift Rq by a few 103 over 12 decades, safely below current lattice sensitivity. Because Koide's formula for charged leptons1 stands out for its remarkable precision and simplicity, various attempts have been made to extend Koide-like relations to quarks,9,10 but these either required additional parameters or achieved less precision than the original lepton relation. Sumino11 proposed a family gauge symmetry to explain the stability of the Koide relation against quantum corrections, while Rodejohann and Zhang12 extended Koide-type relations to the neutrino sector. Recent years have seen increased interest in modular symmetries as the origin of °avor structure. Feruglio4 and others have developed frameworks where modular forms constrain Yukawa couplings, while Kobayashi et al.13 connected these ideas to intersecting D-brane models. The precision frontier has advanced with improved lattice calculations14 and precise determinations of running quark masses.15 These developments provide crucial empirical inputs against which theoretical proposals can be tested. The use of orbifold CFT in understanding particle physics has roots in early string theory work by Dixon et al.16,17 These constructions naturally encode discrete symmetries and their breaking patterns. Vafa18 established connections between 2650073-2 Uni¯ed fermion mass ratios conformal theories and the emergence of realistic particle spectra, while Hamidi and Vafa19 studied interactions on orbifolds that could give rise to observed Yukawa structures. Our work builds on these foundations by explicitly connecting the mathematical structure of orbifold CFTs to observed mass relations, suggesting a deeper geometric origin of °avor in the Standard Model. Unlike previous quark mass ratio models,9,10 which often require additional parameters or lack RG-invariance, our orbifold CFT Mod. Phys. Lett. A Downloaded from www.worldscientific.com approach identi¯es a parameter-free quark mass-ratio functional whose numerical value is consistent with 1=2, stable across energy scales. Speci¯cally, in this paper, we show that orbifold conformal ¯eld theory16 provides a natural framework that uni¯es both sectors. Orbifold CFT describes a quantum ¯eld theory with discrete symmetries, similar to compacti¯ed dimensions in string theory, where twisted sectors lead to fractional scaling of particle properties. We construct a master formula relating particle masses to their charges, demonstrate its approximate invariance under renormalization group °ow, and present a precise, parameter-free result for quark mass ratios that can be tested with upcoming lattice calculations.20,21 This approach reveals a deeper geometric origin of °avor in the Standard Model, linking low-energy mass patterns to high-energy conformal ¯eld theory structures, with potential implications for °avor mixing, and beyond-Standard-Model physics. Aim of this work. The orbifold/simple-current constraints are used to motivate and ¯x a minimal, parameter-free functional Rðk; qÞ of charged-fermion masses for by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. each charge family. The paper then tests whether the resulting specializations R‘ and Rq take simple rational values in Nature and studies their RG stability. We do not derive the numerical values 2=3 or 1=2 from ¯rst-principles CFT operator algebra; rather, we show that the CFT-motivated invariant, once ¯xed, evaluates to those values (within uncertainties) on current data and remains stable under running. Scope and limitations. We do not yet derive the exponents of the power sums on mass eigenvalues from ¯rst principles. Rather, the twist structure motivates two independent homogeneous \degrees" (chiral k, left{right 2k), and we choose the unique permutation-symmetric, degree-zero monomial consistent with charge selection. We test this ansatz against data and renormalization-group running. A full derivation would require proving that Yukawa matrices in the relevant orbifold class admit a simple-current grading whose spectral invariants collapse to these two degrees; we leave this operator-level derivation to future work. 1. Master Formula and Parameter Fixing 1.1. Assumptions We work within the following explicit framework: A0 (Diagonal modular invariant). We consider diagonal modular invariants,  and  ¼ 2h. This covers the standard orbifold vacua used in SM-like model so h ¼ h building. For non-diagonal invariants the second exponent depends on h þ h  in the relevant blocks; we defer that analysis. 2650073-3 S. Varma A1 (Simple-current Z k orbifold). For each charged-fermion sector considered, we assume the dominant Yukawa couplings are governed by a Zk simplecurrent selection rule, generated by J with J k ¼ 1.a A2 (U(1) current algebra factor). A U(1) current algebra factor at level L, with hQ ¼ Q2 =ð2LÞ. The minimal charge unit is q > 0; SM charges are Q ¼ nq. In what follows, q denotes the minimal charge quantum in the chosen U(1) normalization: q ¼ 1 for charged leptons (in units of jej), q ¼ 13 for quarks (in units of the Mod. Phys. Lett. A Downloaded from www.worldscientific.com down-quark charge magnitude). In Appendix B we work with integer-normalized hypercharge Q ~ ¼ 6Y at level L ¼ 2 with q ¼ 1 for clarity; this is equivalent to the main-text choice q ¼ 13 in electric-charge units. A3 (Permutation symmetry and homogeneity). The observable is symmetric under family permutations and homogeneous of degree zero under mi ! mi . Because the ratio is de¯ned purely on mass eigenvalues (not on gauge quantum numbers), it treats the masses as an unordered multiset. For the combined quark ratio Rq , the symmetry group is S6 : invariance under any permutation of the six quark mass eigenvalues fu; c; t; d; s; bg. A4 (Leading-twist dominance, UV-motivated). In weakly-coupled orbifolds, excited/higher twist blocks are suppressed by ðE=Ms Þ2n (oscillator insertions) and/or by eSws (worldsheet instantons). For typical compacti¯cations with Ms  1016 {1018 GeV and electroweak energies E  v  102 GeV, oscillator suppression gives ðE=Ms Þ2  1028 per insertion; worldsheet-instanton factors with by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. generic wrapped areas yield eSws . e10  104:5 for generic compacti¯cation choices. We therefore retain only the minimal twist orders in our parameter-free invariant. Justi¯cation. In simple-current Zk orbifolds, ¯elds in a twisted sector with twist P fraction r=k have minimal conformal weight hmin ðrÞ ¼ i 12 rki ð1  rki Þ. Consequently, contributions from subleading twists (r > 1) and oscillator excitations carry larger h and enter 4D couplings with suppressions / ðE=Ms Þh as well as moduli{overlap factors, making them negligible for E  Ms . Mixing between distinct twists arises through higher-dimensional operators and organizes into the same symmetric{homogeneous structure, so to ¯rst nontrivial order such e®ects correct the numerator and denominator of Rðk; qÞ coherently and are parametrically subleading. This motivates retaining only minimal twists in our parameter-free invariant, while treating subleading-twist e®ects as controlled corrections (see also Appendix A). A5 (Single Higgs VEV for the sector used in R). We restrict to the SM (one Higgs doublet) so that all charged-fermion masses are proportional to the same v. In multi-Higgs models, sectoral ratios (e.g. Ru for up-type, Rd for down-type) where the VEV cancels within each sector are the correct objects. a For classic constructions and ¯xed-point resolution in simple-current modular invariants see Refs. 22 and 23; for orbifold CFT foundations see Refs. 16 and 24. 2650073-4 Uni¯ed fermion mass ratios 1.2. Construction of the minimal invariant Mod. Phys. Lett. A Downloaded from www.worldscientific.com Two minimal homogeneous degrees from twist structure. In a Zk orbifold with diagonal invariant (A0), twisted sectors introduce two family-symmetric P 1=k building blocks labeled by the minimal twist orders: S1 ¼ i m i (chiral degree) P 1=ð2kÞ and S2 ¼ i m i (left{right combined degree, where the factor of 2 arises from  ¼ 2h). We do not claim these are raw CFT observables; rather, they are the building blocks of the unique invariant under the stated constraints (A0{A3). Construction. Any symmetric homogeneous observable constructed from S1 and S2 is of the form S 1 S 2 . Homogeneity (A3) requires =k þ =ð2kÞ ¼ 0, hence by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles.  ¼ 2. Minimal charge selection ¯xes the smallest rational  ¼ kq consistent with Eq. (5). Thus the unique parameter-free choice is P  1=k kq i mi Rðk; qÞ ¼ P  : ð2Þ 1=ð2kÞ 2kq i m i Substituting ðk; qÞ ¼ ð1; 1Þ for leptons and ð2; 13 Þ for quarks ¯xes the lepton and quark specializations of the invariant, R‘  Rð1; 1Þ and Rq  Rð2; 13 Þ. Their numerical values are then obtained by evaluating the resulting expressions on the measured mass eigenvalues at a speci¯ed renormalization scale. With current inputs, one ¯nds R‘ ’ 2=3 (Koide) and Rq ’ 1=2 within uncertainties. In this sense the orbifold argument ¯xes the form of the relation, while the appearance of a neat rational value is an empirical outcome of the parameter-free specialization. This is the CFT-motivated ansatz we test in what follows. 1.3. Monodromy constraint and parameter ¯xing Monodromy{twist identi¯cation. For the simple current J that shifts charge by q, the monodromy charge of a ¯eld of charge Q is Qq QJ ðQÞ ¼  ðmod 1Þ: ð3Þ L (The sign is conventional; modulo 1, Qq=L are equivalent. Derivation in Appendix A.) 2650073-5 S. Varma In a Zk orbifold one has kQJ ðQÞ 2 Z. Writing the twisted-sector label as s, nonvanishing correlators enforce Qq X s  k ðmod kÞ; si  0 ðmod kÞ: ð4Þ L i In particular, the global power that counts minimal charge insertions is ¯xed by the rational kq, which appears in Eq. (2). Mod. Phys. Lett. A Downloaded from www.worldscientific.com Modular-invariance constraint. Simple-current modular invariance is equivalent to monodromy-charge integrality for all ¯elds in the spectrum. For the U(1) level-L factor and the simple current that shifts charge by q, this gives the familyuniform rational kq2 Nrat ¼2 Q; ð5Þ L held ¯xed within a given charge family. For SM charges, the minimal ðk; qÞ are: leptons ð1; 1Þ, quarks ð2; 13 Þ.b Remark (universality). The construction only requires (A0{A5): a Zk simplecurrent action with diagonal invariant and a U(1)L factor. This covers the standard geometric Zk orbifolds and many intersecting-brane models at the level of their worldsheet selection rules. Other constructions may also realize the same constraints; our claim is universality within this class, not uniqueness of the UV completion. by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. 2. Renormalization-Group Stability This section demonstrates that the mass ratio Rq in Eq. (2) is approximately invariant under Standard-Model (SM) renormalization-group (RG) evolution at one loop and quanti¯es the size of two-loop corrections. Road map. (i) Show that the one-loop RGEs for fermion masses are °avor diagonal in each SM sector. (ii) Derive the logarithmic derivative  d d log R and identify coe±cients i . (iii) Show that one exact identity and one numerically small hypercharge combination make the one-loop R small. (iv) Evaluate the two-loop anomalous dimensions to bound the residual running. 2.1. One-loop RG equations for fermion masses At one loop the running of a diagonal fermion mass mf obeys25,26 ! d 1 X 2 X3 2  log mf ¼ c 0y 0  dfa g a ; ð6Þ d 162 f 0 ff f a¼1 b Historically, R‘ ¼ 2=3 is quoted using charged-lepton pole masses. Using MS masses at  ¼ MZ shifts R‘ by Oð103 Þ; unless stated otherwise our statements about R‘ refer to the pole-mass de¯nition. 2650073-6 Uni¯ed fermion mass ratios where yf 0 are Yukawa couplings and ga the gauge couplings. For each SM family the matrices are °avor diagonal, cff 0 ¼ cf ff 0 and dfa ¼ da ff 0 : 8 <3 ðquarksÞ; cf ¼ 3 dfa ¼ ðd1 Q 2f ; d2 ; d3 Þ ð7Þ : ðcharged leptonsÞ; 2 with explicit coe±cients in Ref. 27. Mod. Phys. Lett. A Downloaded from www.worldscientific.com 2.2. Beta function of the mass ratio For the universal ratio, P  1=k kq i mi Rðk; qÞ ¼ P  ; 1=ð2kÞ 2kq i mi its logarithmic derivative is d X d  log R ¼ i  log mi ; ð8Þ d i d 1=k 1=ð2kÞ m m i ¼ q P i 1=k  q P i 1=ð2kÞ : ð9Þ j mj j mj by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. 2.3. Cancellation mechanism The identity X f ¼ 0 ð10Þ f P follows directly from Eq. (9). The U(1)Y combination f f Q 2f is numerically small for the exponents 12 ; 14 , leaving a suppressed one-loop remainder. Numerically at MZ , P 2 3 f f Q f  Oð10 Þ, rendering the hypercharge term subleading. Substituting Eq. (6) gives 2 d 1 6 X  log R ¼ 4 cf f y 2f d 16 2 |{z} f flavorindep: 3 X 2 X 2 X 7  da ga f  d 1 g 1 f Q 2f 5: ð11Þ a |{z} f f flavorindep: P Gauge terms proportional to SU(2)SU(3) cancel exactly because f f ¼ 0. The remaining U(1)Y term is numerically small{both because g1  g3 and because P the up/down sector contributions to f f Q 2f nearly cancel for the exponents 12 and 2650073-7 S. Varma Mod. Phys. Lett. A Downloaded from www.worldscientific.com Fig. 1. The ratio Rq across energy scales from MZ to 1012 GeV with 1 uncertainty band, demonstrating stability around Rq ¼ 1=2. Despite individual quark masses spanning ¯ve orders of magnitude and running signi¯cantly, this dimensionless combination remains stable within a few 103 . Data from Ref. 15, evolved using two-loop SM RGEs. Table 1. Values of Rquark at di®erent energy scales using running quark masses from Ref. 15. by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. Scale (GeV) Theory Rquark 91:1876 (MZ ) Full SM 0:49954ð17Þ 125:1 (Mh ) Full SM 0:49970ð17Þ 173:1 (Mt ) Full SM 0:49981ð17Þ 105 Full SM 0:50234ð20Þ 108 Full SM 0:50398ð22Þ 1012 Full SM 0:50533ð24Þ 1 4 . The net one-loop  R is consistent with zero within our uncertainties; the full twoloop running shifts Rq by a few 103 up to 1014 GeV (Fig. 1; Table 1). Two-loop stability. The two-loop anomalous-dimension matrices of Ref. 27, whose °avor structure is again diagonal, introduce gauge-Yukawa cross terms proportional to the same charge combinations that vanish at one loop. Evaluated with these matrices, the residual contribution shifts Rq;‘ by a few 103 between MZ and 1014 GeV   safely below current lattice uncertainties. The one-loop invariance therefore su±ces for present phenomenology. 3. Numerical Test and High-Scale Prediction Unless stated otherwise, all quark masses are MS parameters quoted at  ¼ MZ ; central values and uncertainties are taken from Ref. 15 (see Table 2). We evolve with 2650073-8 Uni¯ed fermion mass ratios two-loop SM RGEs including thresholds at mc ; mb ; mt and propagate uncertainties using the Monte-Carlo procedure described in Sec. 4. To quantitatively test our ansatz, we examine the quark mass ratio Rq across energy scales. Using the six MS quark masses from Ref. 15, we compute Rq at the electroweak scale (MZ  91:2 GeV). The values mu ¼ 0:00123 GeV, md ¼ 0:00267 GeV, ms ¼ 0:05316 GeV, mc ¼ 0:620 GeV, mb ¼ 2:839 GeV, and mt ¼ 168:26 GeV yield: Mod. Phys. Lett. A Downloaded from www.worldscientific.com P 1=2 2=3 m Rq ðMZ Þ ¼ P ¼ 0:4995  0:0017; ð12Þ ð m1=4 Þ4=3 where the uncertainty is propagated via Monte Carlo (Sec. 4). The stability of Rq as shown in Fig. 1 is remarkable given that the individual quark masses span ¯ve orders of magnitude and run signi¯cantly across 10 orders of magnitude of energy scale. This stability provides strong evidence for the fundamental nature of our relation. To produce a de¯nitive high-scale prediction, we evolved the six quark masses to  ¼ 1014 GeV using RunDec v3.128 with the four-loop QCD -function and mass ð5Þ anomalous dimension in a six-°avor scheme. The calculation uses  s ðMZ Þ ¼ 0:1179  0:0009 and incorporates two-loop electroweak corrections. At this energy scale, we ¯nd (in GeV): mu ¼ 0:00049ð8Þ, md ¼ 0:00107ð8Þ, ms ¼ 0:0213ð17Þ, mc ¼ 0:248ð7Þ, mb ¼ 1:127ð11Þ, and mt ¼ 67:6ð4Þ, yielding: by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. Rq ð Þ ¼ 0:500  0:002: ð13Þ This result is statistically indistinguishable from 1=2 (deviation < 0:1 ).c The high-energy prediction in Eq. (13) can be tested with upcoming lattice determinations. Building on recent nonperturbative step-scaling work by the ALPHA Collaboration,20 ongoing Nf ¼ 2 þ 1 þ 1 programs aim to push quark-mass determinations to   1014 GeV with  0:2% precision, su±cient to con¯rm or falsify our prediction. 4. Monte Carlo Error Propagation and Raw Data We performed a comprehensive statistical analysis to assess the stability of our quark mass ratio Rq under uncertainties in the input quark masses. This section presents the methodology and detailed results. 4.1. Quark mass inputs and uncertainties Table 2 shows the running quark masses at di®erent energy scales used in our analysis, with their associated uncertainties. These values are based on Ref. 15. c Caveat. The evolution above assumes the Standard Model particle content between MZ and H . Additional thresholds or BSM interactions would perturb the near-invariance of Rq ; in such scenarios Eq. (8) implies controlled, model-dependent deviations from 12 . 2650073-9 S. Varma Table 2. Running quark masses (in GeV) at di®erent energy scales with uncertainties from Ref. 15. Scale (GeV) mu md ms mc mb mt 91:1876 0.00123(21) 0.00267(19) 0.053(5) 0.620(17) 2.839(26) 168.3(8) 125:1 0.00120(21) 0.00261(18) 0.052(5) 0.607(17) 2.768(26) 165.1(8) 173:1 0.00118(20) 0.00256(18) 0.051(4) 0.594(17) 2.702(25) 162.0(8) 105 0.00086(15) 0.00188(13) 0.0375(33) 0.435(13) 1.908(21) 123.8(9) 108 0.00069(12) 0.00152(11) 0.0303(27) 0.350(11) 1.502(18) 102.5(9) Mod. Phys. Lett. A Downloaded from www.worldscientific.com 1012 0.00056(10) 0.00124(9) 0.0248(22) 0.283(9) 1.194(15) 85.1(9) 4.2. Monte Carlo methodology We performed a 10,000-trial Monte Carlo simulation to propagate the uncertainties in quark masses to the ratio Rq . For each trial: (1) We randomly sampled each quark mass from a Gaussian distribution centered at its central value with P standard deviation equal to its uncertainty. ð m1=2 Þ2=3 (2) We calculated Rq ¼ ðP m1=4 Þ4=3 using these sampled masses. (3) We repeated this process for each energy scale. 4.3. Statistical results Table 3 summarizes the statistical results from our Monte Carlo analysis with data by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. sourced directly from Ref. 15. For the high-scale prediction at 1014 GeV, we evolved these masses using RunDec v3.1, as described in Sec. 3, yielding Rq ¼ 0:4998 0:0018. Figure 2 shows the distribution of Rq values at di®erent energy scales. 4.4. Analysis of stability Our Monte Carlo analysis reveals several important features: (1) Despite the relatively large uncertainties in individual quark masses (particularly for light quarks), the ratio Rq remains remarkably stable, with 95% of values falling within the range 0.497–0.508 across all energy scales. Table 3. Monte Carlo analysis results for Rq at di®erent energy scales. Scale (GeV) Mean Rq 91:1876 (MZ ) 0:4996ð13Þ 125:1 (Mh ) 0:4998ð13Þ 173:1 (Mt ) 0:4999ð13Þ 105 0:5024ð13Þ 108 0:5041ð14Þ 1012 0:5054ð14Þ 2650073-10 Uni¯ed fermion mass ratios Mod. Phys. Lett. A Downloaded from www.worldscientific.com Fig. 2. Histograms of Rq values from 10,000 Monte Carlo trials at di®erent energy scales. Despite individual quark masses spanning ¯ve orders of magnitude, the ratio Rq remains remarkably stable around 1=2 across all energy scales. (2) This stability is particularly impressive given that the quark masses themselves span ¯ve orders of magnitude and run signi¯cantly with energy. (3) The ratio Rq shifts by only about 1% when moving from the electroweak scale to 1014 GeV, consistent with the expected higher-order corrections beyond by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. one-loop. The robustness of Rq  1=2 under these large uncertainties strongly suggests that this is not a numerical coincidence but re°ects a deeper underlying principle consistent with the orbifold-CFT-motivated invariant. 5. Extensions 5.1. Neutrino case: Dirac versus Majorana Our invariant targets charged fermions. For neutrinos we make no nontrivial prediction; if one formally sets q ¼ 0 in Eq. (2), the expression reduces to R ¼ 1, but the construction is not intended for Majorana masses, which arise from a di®erent operator class with di®erent selection rules. 6. Phenomenological Outlook Our orbifold CFT approach to fermion mass relations has several signi¯cant implications for particle physics. The constructed relation provides a parameter-free theoretical framework that uni¯es lepton and quark mass patterns, suggesting a deeper geometric origin of °avor structure. 2650073-11 S. Varma The most immediate application is our prediction for the quark mass ratio at high energy scales. With Rq ð Þ ¼ 0:500  0:002 at  ¼ 1014 GeV, we provide a precise target for upcoming lattice QCD calculations. Recent advances in nonperturbative step-scaling techniques21 are bringing the required precision within reach. A veri¯cation of this prediction would provide strong evidence for the fundamental role of conformal structures in determining fermion mass patterns. The modular invariance constraint that ¯xes our parameters also implies correMod. Phys. Lett. A Downloaded from www.worldscientific.com lations in the °avor mixing sector. Because the o®-diagonal Yukawa elements transform in the same orbifold representations, the same orbifold reasoning that yields mass relations should generate constraints on CKM matrix elements when applied to o®-diagonal Yukawa couplings. A systematic study of these correlations is in progress and could potentially explain the observed hierarchical structure of °avor mixing. A preliminary exploration of the possible CKM-matrix correlations is given in Appendix C. Beyond the Standard Model, any exotic fermion with a di®erent charge, or any extended °avor symmetry that alters the modular-invariance condition, would modify our predicted ratios. This provides a clean discriminator between models with vector-like quarks or horizontal symmetries. The appearance of orbifold modular constraints in low-energy °avor observables hints at a geometric origin of fermion masses, reminiscent of how string compacti¯cations tie geometry to particle spectra. by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. 7. Conclusion We have presented a uni¯ed theoretical framework that bridges lepton and quark mass relations through a single, parameter-free formula motivated by orbifold conformal ¯eld theory, o®ering a novel perspective on °avor structure in the Standard Model. Our approach reveals that the well-known Koide relation for leptons and a newly identi¯ed quark mass relation emerge from the same underlying principle: modular invariance of twisted sectors in orbifold CFTs. The key ¯ndings are: . A single modular-covariant formula X kq . X 2kq R¼ m1=k m1=ð2kÞ applies to leptons and quarks once ðk; qÞ are ¯xed by their SM charge assignments. . For charged leptons the formula reproduces the known Koide result R‘ ¼ 2=3; for quarks it predicts Rq ¼ 0:500  0:002 at  ¼ 1014 GeV — consistent with 12 within uncertainties — which was not ¯t to data. . Flavor-independent one-loop coe±cients make the ratio approximately one-loop RG invariant; two-loop terms shift Rq by a few 103 over 12 decades. . Any exotic charge assignment that alters ðk; qÞ would shift R, providing a sharp discriminator for models with vector-like fermions or horizontal symmetries. 2650073-12 Uni¯ed fermion mass ratios Most importantly, this relation is falsi¯able. With current lattice methods trending toward sub-percent precision, a 0.5% discrepancy in Rq ð1014 Þ GeV would disqualify the minimal orbifold-CFT scenario presented here, assuming no °avornon-universal thresholds between the electroweak and test scales. Limitations. The map \twist degree ! mass-eigenvalue exponent" is not yet proven from ¯rst principles; higher-twist contamination is expected to be UVsuppressed (with ðE=Ms Þ2  1028 or eSws . 104:5 for generic compacti¯cation Mod. Phys. Lett. A Downloaded from www.worldscientific.com choices, A4) but could shift R at the percent level if special moduli enhance nexttwist couplings; multi-Higgs models require sectoral versions. The connection between orbifold modular invariance and fermion mass ratios presented here suggests a deeper geometric origin of °avor structure in the Standard Model. This framework opens new avenues for understanding °avor mixing patterns and provides a distinctive signature for models with exotic fermions or extended °avor symmetries. Appendix A. Monodromy-Charge Lemma For a U(1) factor at level L, hQ ¼ Q2 =ð2LÞ. For the simple current J of charge q, hJ ¼ q2 =ð2LÞ and hJQ ¼ hQþq ¼ hQ þ hJ þ QqL . Thus Qq QJ ðQÞ ¼ hJ þ hQ  hJQ ¼  ðmod 1Þ: ðA:1Þ L by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. The simple-current modular invariant exists i® kQJ ðQÞ 2 Z for all ¯elds kept in the spectrum, giving kQq=L 2 Z. For a ¯xed-charge family Q ¼ nq, the integer 2 2 equals n kqL , so kqL is a family-uniform rational that (once minimized) pins ðk; q; LÞ as used in Eq. (2). This is exactly the constraint quoted in Eq. (5). Appendix B. Worked Model: Z 2 Orbifold of U(1)L Setup. Work with integer-normalized charges ~ Q ¼ 6Y (so that all SM hypercharges become integers). Choose U(1) level L ¼ 2 and simple current J that shifts charge by q ¼ 1. Then h ~ ¼ Q~ 2 =ð2LÞ ¼ Q ~ 2 =4 and h ~ ¼ h ~ þ hq þ Qq=L ~ ~ ¼ h ~ þ h1 þ Q=2. Q Qþq Q Q Monodromy charge and twist label. The monodromy charge is ~ ~ ~ ¼  Qq ¼  Q ðmod 1Þ: QJ ðQÞ ðB:1Þ L 2 Since J 2 ¼ 1, modular invariance requires 2QJ ð ~ QÞ 2 Z, which holds for all integer ~ Q. The twisted-sector label is ~ Qq s  2  Q ~  ~Qðmod 2Þ: ðB:2Þ L Hence: . ~ ¼ þ4; dR : Q Right-handed singlets (uR : Q ~ ¼ 2; eR : Q ~ ¼ 6) have s ¼ 0 (even Q). ~ . EW doublets (QL : Q~ ¼ þ1; LL : Q~ ¼ 3; H: Q ~ ¼ þ3) have s ¼ 1 (odd ~ Q). 2650073-13 S. Varma Yukawa selection. The up-type Yukawa involves Q  L, H ~ , uR . With integernormalized charges ~Q ¼ 6Y , one has ~ QðQ L Þ ¼ 1, ~ QðH ~ Þ ¼ 3, ~ QðuR Þ ¼ þ4, so  1  3 þ 4 ¼ 0. Twist labels add sðQ L Þ þ sðH Þ þ sðuR Þ  1 þ 1 þ 0  0 ðmod 2Þ. ~ ~ Q The down-type case with H is analogous: Qð  L Þ þ QðHÞ ~ ~ R Þ ¼ 1 þ 3  2 ¼ 0 þ Qðd and 1 þ 1 þ 0  0 ðmod 2Þ. Thus both Yukawas satisfy U(1) and twist selection rules. Two minimal twist orders. For a Z2 orbifold, the holomorphic minimal twist Mod. Phys. Lett. A Downloaded from www.worldscientific.com gives a homogeneous degree corresponding to exponent 12 in symmetric sums; the left{right combined block gives exponent 14 . These de¯ne the two family-symmetric building sums across the quark sector: X 1=2 6 X6 1=4 S1 ¼ m i ; S2 ¼ mi : ðB:3Þ i¼1 i¼1 Charge-selection power and ratio. Global U(1) selection ¯xes the minimal power 2q ¼ 2 ð 13 Þ ¼ 23 in SM hypercharge units; the unique symmetric homogeneous monomial is 2=3 S1 Rq ¼ 4=3 : ðB:4Þ S2 This is precisely the k ¼ 2 specialization of the master expression Eq. (2), and is what we test in Sec. 3. by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. Worldsheet{to{4D bridge. In orbifold compacti¯cations, 4D Yukawas come from sphere three-point functions of vertex operators. The s-label split (doublet versus singlet) organizes which three-point couplings are allowed pre-EWSB; after EWSB and diagonalization, the physical mass eigenvalues lose memory of the pre-EWSB slabels. The ¯nal ratio acts on eigenvalues post-EWSB and does not track s labels. The two building sums S1 ; S2 arise because the twist structure organizes allowed Yukawa blocks, and these blocks produce the two homogeneous degrees we use. Relation to the lepton-sector specialization. This appendix is intended as an explicit illustration of how a k ¼ 2 (Z2 ) simple-current selection rule can underwrite the quark-sector building blocks used in Rq . The charged-lepton Koide specialization corresponds to the k ¼ 1 case and need not arise from the same single Z2 factor in this toy construction. In realistic internal CFTs (often products with multiple simple currents and multiple twisted sectors), di®erent Yukawa couplings can be controlled by di®erent e®ective simple-current gradings; in that sense the lepton and quark specializations can coexist within one overall background, even though this appendix isolates only the quark-sector Z2 ingredient. Appendix C. CKM Correlations Outline The orbifold CFT approach that predicts mass relations can be extended to the °avor mixing sector. Here, we outline how similar modular invariance constraints may generate relations among CKM matrix elements. 2650073-14 Uni¯ed fermion mass ratios The Yukawa couplings in the Standard Model can be written as  iH LYukawa ¼ y uij Q  i HdRj þ h:c: ~ uRj þ y dij Q ðC:1Þ When electroweak symmetry is broken, these couplings give rise to mass matrices: v v M iju ¼ y uij pffiffi ; M ijd ¼ y dij pffiffi : ðC:2Þ 2 2 Mod. Phys. Lett. A Downloaded from www.worldscientific.com The CKM matrix arises from diagonalizing these mass matrices: VCKM ¼ ðU Lu Þ† U Ld ; ðC:3Þ where U Lu and U Ld are the unitary matrices that diagonalize M u and M d . Extending our orbifold reasoning to o®-diagonal elements will require imposing the charge-uniformity condition on mixed Yukawa indices. Preliminary algebra shows that this produces constraints linking certain mass ratios to mixing angles. A detailed analysis is underway; initial results suggest the framework may shed light on: (1) The hierarchical structure of the CKM matrix (jVus j jVcb j jVub j). (2) The correlation between quark mass hierarchies and mixing angles. (3) The near-unity value of jVtb j: by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. This extension would provide a uni¯ed geometric framework for understanding both mass ratios and mixing patterns in the Standard Model. Acknowledgments The author received no speci¯c funding for this work. The author gratefully acknowledges the use of computational resources provided by VS Asset Management, LLC. The author has no competing interests to declare. All numerical data and code used in this work are available from the corresponding author upon reasonable request. ORCID Samir Varma https://orcid.org/0000-0001-6437-5395 References 1. Y. Koide, Phys. Lett. B 120, 161 (1983). 2. H. Georgi and C. Jarlskog, Phys. Lett. B 86, 297 (1979). 3. H. Fritzsch, Phys. Lett. B 73, 317 (1977). 4. F. Feruglio, Frascati Phys. Ser. 64, 174 (2017). 5. Z.-Z. Xing, Phys. Rep. 854, 1 (2020). 6. M. Gell-Mann, The eightfold way: A theory of strong interaction symmetry, California Institute of Technology Synchrotron Laboratory Report No. CTSL-20 (1961). 2650073-15 S. Varma 7. H. Fritzsch, Nucl. Phys. B 155, 189 (1979). 8. S. Weinberg, Trans. N. Y. Acad. Sci. 38, 185 (1977). 9. A. Kartavtsev, arXiv:1111.0480 [hep-ph]. 10. J.-M. Gerard, F. Go±net and M. Herquet, Phys. Lett. B 631, 125 (2005). 11. Y. Sumino, J. High Energy Phys. 2009, 075 (2009). 12. W. Rodejohann and H. Zhang, Phys. Lett. B 698, 152 (2011). 13. T. Kobayashi, S. Nagamoto and S. Uemura, Prog. Theor. Exp. Phys. 2017, 023B02 (2017). Mod. Phys. Lett. A Downloaded from www.worldscientific.com 14. Y. Aoki et al., Eur. Phys. J. C 82, 869 (2022). 15. G.-y. Huang and S. Zhou, Phys. Rev. D 103, 016010 (2021). 16. L. Dixon, J. A. Harvey, C. Vafa and E. Witten, Nucl. Phys. B 261, 651 (1985). 17. L. Dixon, J. A. Harvey, C. Vafa and E. Witten, Nucl. Phys. B 274, 285 (1986). 18. C. Vafa, Phys. Lett. B 206, 421 (1988). 19. S. J. Hamidi and C. Vafa, Nucl. Phys. B 279, 465 (1987). 20. ALPHA Collab. (M. Bruno et al.), Phys. Rev. Lett. 119, 102001 (2017). 21. I. Campos, P. Fritzsch, C. Pena, D. Preti, A. Ramos and A. Vladikas, Eur. Phys. J. C 78, 387 (2018). 22. A. N. Schellekens and S. Yankielowicz, Nucl. Phys. B 327, 673 (1989). 23. J. Fuchs, A. N. Schellekens and C. Schweigert, Commun. Math. Phys. 180, 39 (1996), Revised version (v3) February 1998. 24. L. Dixon and J. A. Harvey, Nucl. Phys. B 274, 93 (1986). 25. M. E. Machacek and M. T. Vaughn, Nucl. Phys. B 222, 83 (1983). 26. R. Alonso, E. E. Jenkins, A. V. Manohar and M. Trott, JHEP 04, 159 (2014). 27. T. van Ritbergen, J. Vermaseren and S. Larin, Phys. Lett. B 400, 379 (1997). 28. F. Herren and M. Steinhauser, Comput. Phys. Commun. 224, 333 (2018). by 74.108.22.145 on 03/13/26. Re-use and distribution is strictly not permitted, except for Open Access articles. 2650073-16