Quantiﬁcation of Sub-Solar Star Ages from the Symmetry of Conjugate Histograms of Spin Period and Angular Velocity

: Empirical laws proposed for the decline in star spin with time have heretofore been tested using ambiguous ﬁtting models. We develop an analytical inverse model that uses histogram data to unequivocally determine the physical law governing how dwarf star spin depends on time ( t ) and mass ( M ). We analyze shapes of paired histograms of axial rotation period ( Π ) and angular velocity ( ω = 2 π / Π ) to utilize the fact that a variable and its reciprocal are governed by the same physics. Copious data on open clusters are used to test the formula ∂ ω / ∂ t ∝ − ω n where n is unrestricted, and thus covers diverse possibilities. Histogram conjugates for each of 15 clusters with 120 to 812 measurements provide n = 1.13 ± 0.19. Results are independent of initial spin rate, bin size, cluster parameters, and star mass. Notably, 11 large clusters with mostly M-types yield ﬁts with n = 1.07 ± 0.12. Associations behave similarly. Only exponential decay ( n = 1) explains the similar shapes of the conjugate histograms for the spin period and angular velocity, despite the asymmetric (inverse) relationship of these variables. This rate law is consistent with viscous dissipation. Forward modeling conﬁrms that n is near unity and further shows that coeval formation of all stars in a cluster does not occur. We therefore explore a constant rate of star production, which is reasonable for tiny stars. Inverse models show that episodic production increases with mass, but is unimportant below ~0.55 M Sun . We infer star and cluster ages, and ﬁnd that star production becomes less regular with time, as interstellar gas and dust are progressively depleted. Our new analytical approach of extracting a physical law from conjugate histograms is general and widely applicable.


Introduction
Star ages are very important to astronomy as they are intimately linked to stellar processes and evolution (e.g., [1][2][3]). Of all the stars, only our Sun has a known age, which is based on consistency among isotopic data for diverse meteorite types [4]. Ages of other stars are inferred from various models of stellar evolution where cluster ages are sometimes incorporated, but these too are derived from models.
One ubiquitous and measurable indicator of a star's age is its axial spin rate, which declines with time: young stars have periods (Π) as short as a few hours [5,6], whereas our ancient, 4.54 Ga (4.54 × 10 9 year-old) Sun spins slowly with Π~24.5 d at its equator. Different rate laws for spin decay with time have been proposed [7,8], but have not been validated [9]. The currently popular approach of gyrochronology is detailed below along with its known problems [10,11].
Given the lack of isotopic data on star ages, a statistical test of spin decay models is appropriate. Robust tests are possible because the database on spin rates now includes 11 open clusters that each have period measurements for >243 stars, plus several clusters and associations with at least 120 measurements (Table 1). Statistics were once used to study stellar rotation [12], but this approach was abandoned, in part because the methods used to measure spin rates have changed. Prior to~1996, Doppler methods provided apparent Symmetry 2021, 13, 1519 2 of 31 tangential velocities (vsini) of rotating stars. Converting these data to rotation periods (Π) requires estimating both stellar radius (R) and orientation (i) of the spin axis, because Π = 2πR/v equator . Moreover, period is undefined (not measureable) at an inclination of 0 • .   [39,40]. Open squares and fine curve = Pleiades [41,42]. Grey triangles and grey curve = Hyades [43]. Filled diamonds and dashed = Praesepe [44,45]; (b) Density (red color and left axis) and blackbody temperature (blue color and right axis) of main sequence stars. Symbols are the values listed in the summary tables of Zombeck [38]. Lines are least square fits of each dataset assuming a power law dependence on mass: each line is labeled with its fits. Horizontal lines show ambient densities for select materials. Main sequence star classifications (labels M8 to O5) are based on color which encapsulates all properties, except spin, of main sequence stars.
Recent studies furthermore characterize the stars by a color index, with reference to various filters. Frequently, the B-V filter is used (Table 1). Color index is important because, except for spin period, it links to all major star properties: surface temperature, luminescence, radius, and mass (M), through simple proportionalities (e.g., [37,38]; Figure 1b). Thus, starspot studies provide both angular velocities (ω) and equatorial tangential velocity via: Given the interconnected physical properties, the interrelationships of star spin, age, and mass are of great interest.

Purpose
We use a new tactic to establish the rate law for spin-down, which is to examine conjugate pairs of histograms of the same data. Our inverse method (Section 3) is mathematically novel but straightforward, and makes very few assumptions. Examining both Πand ω-histograms is essential because histograms of angular velocity place emphasis on old stars whereas histograms of period emphasize young stars. We also develop a new method to analyze histograms that are incomplete or under-sampled, and utilize forward modeling to test assumptions.

Organization
Section 2 provides background information on analytical inverse models because this approach is uncommon. Section 2 also summarizes the existing database on star spin, and the currently popular methodologies. Section 3 presents our analytical mathematics. Section 3.1 constructs histograms for the endmember case of coeval production and of stars forming with random initial spin rates, to evaluate the feasibility of these popular scenarios. Section 3.2 develops our basic model. Sections 3.3 and 3.4 show that assuming steady-state production of stars with time has a minor effect in our approach. Section 3.5 develops a new method for removing possible effects of bin size in histogram analysis. Section 4 Symmetry 2021, 13, 1519 4 of 31 demonstrates that exponential decay best represents star spin data and that other proposed rates of decay are inconsistent with observed histogram shapes. Section 5 evaluates whether possible mechanisms are consistent with exponential spin loss. Section 6 calibrates our formula for the dependence of spin rate (or period) on stellar age and M/M Sun , using the age of the Sun and some limiting values.
Astronomical implications of exponential spin-down for main sequence dwarfs are covered in Section 7. This concluding section also addresses wide applicability of our analytical inverse methods to statistical data bases that can be represented as histograms.

Inverse vs. Forward Modeling
A mainstay of modern science is using forward models with many parameters to explain data. Fitting has utility in that measurements can be reproduced and trends extrapolated, but the drawback is that an incorrect physical model can appear to match data via trade-offs among unconstrained input parameters, particularly when multiplied [46]. In contrast, analytical inverse models provide unique solutions [47,48]. Many analytical inverse approaches are historic. Groetsch [47] provides diverse examples from many fields. The astronomical example most closely related to the goals and approach of the present paper is the extraction of the proper motions of stars from observed radial velocities [49]. Recently, an inverse approach was used to extract mass and density of spiral galaxies from measurements of their rotational velocities as a function of radius [50].

Useful Mathematical Attributes of Histograms
Histograms are commonly used to visualize data. A histogram is described by a density function that has a corresponding distribution function that embodies the underlying physics [51]. Histograms of a variable and of its reciprocal typically have different shapes because their density functions differ in mathematical form, yet the underlying physical law is the same. On this basis, we develop an inverse method to extract a physical law from histogram shapes. Our approach assumes that one independent variable largely controls the dependent variable embodied by any given histogram. Consequently, the reciprocated dependent variable and its histogram are also controlled by the same independent variable. Our method is general, but as is common with inverse modeling, the mathematical model is deduced through solving a specific problem [47].

The Database on Spin Periods of Dwarf Stars
Star periods from analyzing light curves (intensity vs. time) are independent of the axial orientation of the star. Data quality is covered in many publications on light curves of cluster stars. Issues such as biases, measures of completeness, and minimal contamination from stars outside the cluster are addressed in detail by Hartman et al. [27]. Briefer discussions exist in Irwin et al. [16,19,25], and in most studies since 2015.
Importantly, spin values do not depend strongly on the research group: c.f. Irwin et al. [25] to Jackson and Jeffries [52]. An example is the Hyades, where 90% of the stars examined in more than one study have periods that match to better than 10% [53]. For Praesepe, periods for half of the 214 stars measured multiple times agree within 2%, and

The Da
Star p axial orien of cluster s tion from discussion Impo al. [25] to examined Praesepe, ¾ of the re than 10%, ogy seem m measurem measured of the repeats agree within 5% [29]. Thus, the available data should be accurate to better than 10%, with some outliers. Differences among studies prior to 2018 due to methodology seem moot. Notably, more recent studies use a different approach. Statistically, period measurements are no different, as shown by Cui et al. [54] via comparison with 1107 stars measured by McQuillen et al. [55]: the majority of periods agree within 1%. However, mass and luminosity are determined in a different manner in the most recent studies, and so these were not used to ascertain the effect of mass, unless mass of the stars was included in the tabulated data. Typically, temperatures were listed.
A few star spot studies of clusters were not analyzed because they either cover the same cluster with a similar mass range (cf. [14] to [56]) and provide very similar histograms, Symmetry 2021, 13, 1519 5 of 31 or fewer than 100 stars were measured (e.g., [57,58]). Neither are histograms of aggregated data from Kepler missions presented because, unlike open clusters, star production began at assorted times, which renders extraction of a rate law uncertain.
Importantly, two systematic errors exist in star spot studies.
• An upper limit on Π (and a corresponding lower limit on ω) exists and varies among the studies because starspot lifetimes must exceed the observational duration (e.g., [59]). Hence, stars with periods >25 days are difficult to measure and therefore underrepresented in data from the Kepler campaigns. The reliability study of 38 G-types in M67 found that 75% of the measurements of slow rotations were reliable, but only for Π < 32 days [60]. • For all studies, a lower limit on Π (and an upper limit on ω) is affected by factors such as the time between formation and the time at which a detectable (luminous) star exists. Under-representation of slow-and fast-spinning stars is accounted for in our approach (Section 3).
In summary, Section 4 focuses on post-2008 datasets from the Kepler 2 mission on individual open clusters, which are purported to be more accurate and sample more stars. Most studies yield >243 period measurements on individual open clusters, and describe properties of the open cluster studied (Table 1); also see Finlay [61]. Star spot studies of various groupings of stars near the Sun, associations, field M stars, and ground-based measurements are also analyzed. The reference dataset of Mt. Wilson stars includes sufficient information to convert tangential velocity to periods [36].
Data sets were downloaded from published supplementary material associated with the cited papers. Very few datasets were tabulated in the papers.

Relationships of Mass and Color Needed to Compare Star Spin Studies
This study mostly considers individual clusters. Towards the end of the report, data from the clusters listed in Table 1 are merged to explore of the effect of mass on spin-down.
To merge data, we use M/M Sun rather than color-temperature because some studies used filters other than B-V. Uncertainties in mass were estimated as~8% by Hartman et al. [20]. Moreover, the precise manner in which M/M Sun is computed is not important to our study because mass is only used to separate the stars into broad categories (e.g., M-types vs. K-types).
Reported mass is used, if available. If not, we used results from Zombeck [38] to construct the following empirical formulae for average mass: where R-I is a color index differing from B-V. Using Equation (2) or (3) yielded the same histograms as using the authors' M/M Sun values. Tabulated data for the V-K S index [37] were fit via splining to best cover the wide type ranges of the Pleiades and NCG 2516.

Formulae Previously Applied to Period Data
The currently popular approach, known as gyrochronology, stems from Skumanich's [7] proposal that equatorial rotational velocity (v) degrades in proportion to t −1/2 , where t is star age. The basis was Ca II emission strengths decreasing as rotational velocities decrease [62], plus model ages for the Pleiades and Hyades star clusters. Stars more massive than the Sun were used in the 1972 proposal, which remains a concern because large stars have very short lifetimes [1] and their properties, including chemical composition, are unlike those of the Sun (Figure 1b) and other dwarf stars. Furthermore, angular momentum losses from our Sun measured from 1994-2019 are a factor of two lower than expected from gyrochronology models [10]. Clearly problems exist with t −1/2 decline, yet alternatives such as exponential decay with time, proposed by Terndrup et al. [8], are little explored. This section therefore focuses on recent modifications of Skumanich's 1972 idea. Figure 2a illustrates application of gyrochronology forward modeling to the wellstudied Pleiades. Table 2 lists model assumptions of gyrochronology. In summary, data on plots of period vs. B-V are roughly bracketed by two ad hoc curves referred to as the I-and C-trends. The I-trend is alleged to represent slowly rotating stars. The original formula [63]: with 171 measurements all have n ~ 1.5 (Figure 5k), which is further from unity than values found for younger clusters. Peaks in their distributions cause comparatively poor fits. NGC 681 lacks small M stars. Fitting only the M-types of the other two clusters provides better fits that yield n values closer to unity ( Figure 8; Table 3), as is also the case for the Pleiades (Figure 3). Furthermore, peaks in period increase in height with the modeled age of the clusters for the older clusters (Figure 5g-l), whereas young clusters (Figure 5a-f) lack peaks. Middle-aged Pleiades ( Figure 3) having a very weak peak is consistent with this finding. Existence of peaks in П, weak or strong, is linked to the presence of old, white dwarfs ( Table 3). The above findings are consistent with star production becoming less regular with time as gas is depleted and mass is stored in white dwarfs, which are not included in period measurements. Notably, peaks in the distribution of M-type stars are only associated with the oldest clusters. These findings show that strong departures from steadystate production affect our fits, but not weak departures. Importantly, n is roughly 1 for all clusters, so for sub-solar dwarfs, star mass does not affect the law for spin decay.
In summary, fits in Figures 3, 5, 6 and 7 rule out n < ½ and n > 3/2, and thus unequivocally rule out Skumanich's rule (n = 3) for MKG stars. Crucially, rate laws with such values for n require that the appearances of the П-and ω-histograms differ substantially (details in Sections 4.4-4.6), contrary to the copious observations. Exponential decay, proposed by Terndrup et al. [8], best represents the available spin data on dwarf stars in individual clusters. Moreover, production of small stars is continuous in the youngest clusters, but this condition becomes less appropriate for larger stars as the clusters evolve, whereas production of M-types departs from continuous behavior only for the oldest clusters.
[(B-V − 0. 5) with 171 measurements all have n ~ 1.5 (Figure 5k), which is further from unity than values found for younger clusters. Peaks in their distributions cause comparatively poor fits. NGC 681 lacks small M stars. Fitting only the M-types of the other two clusters provides better fits that yield n values closer to unity ( Figure 8; Table 3), as is also the case for the Pleiades (Figure 3). Furthermore, peaks in period increase in height with the modeled age of the clusters for the older clusters (Figure 5g-l), whereas young clusters (Figure 5a-f) lack peaks. Middle-aged Pleiades ( Figure 3) having a very weak peak is consistent with this finding. Existence of peaks in П, weak or strong, is linked to the presence of old, white dwarfs ( Table 3).
The above findings are consistent with star production becoming less regular with time as gas is depleted and mass is stored in white dwarfs, which are not included in period measurements. Notably, peaks in the distribution of M-type stars are only associated with the oldest clusters. These findings show that strong departures from steadystate production affect our fits, but not weak departures. Importantly, n is roughly 1 for all clusters, so for sub-solar dwarfs, star mass does not affect the law for spin decay.
In summary, fits in Figures 3, 5, 6 and 7 rule out n < ½ and n > 3/2, and thus unequivocally rule out Skumanich's rule (n = 3) for MKG stars. Crucially, rate laws with such values for n require that the appearances of the Пand ω-histograms differ substantially (details in Sections 4.4-4.6), contrary to the copious observations. Exponential decay, proposed by Terndrup et al. [8], best represents the available spin data on dwarf stars in individual clusters. Moreover, production of small stars is continuous in the youngest clusters, but this condition becomes less appropriate for larger stars as the clusters evolve, whereas production of M-types departs from continuous behavior only for the oldest clusters.
was subsequently modified to: [64]. Other empirical relationships with a few free parameters are similar [65]. The narrow range in exponents for time of tasets.
Star types of two clusters differ substantially from those of the other 13, but these two still follow exponential decay. Cepheus OB3b (Figure 5c) has some young stars possibly not on the main sequence, and considerable data on comparatively large stars (types F and G-type), see [73] (their Figure 13). Similarly, NGC 869 (Figure 5e) consists of fast (young) rotators and is dominated by type-G stars [15].
Two older clusters with much data (M37 and NGC 2516, Figure 5j,l) plus NGC 6811 with 171 measurements all have n ~ 1.5 (Figure 5k), which is further from unity than values found for younger clusters. Peaks in their distributions cause comparatively poor fits. NGC 681 lacks small M stars. Fitting only the M-types of the other two clusters provides better fits that yield n values closer to unity ( Figure 8; Table 3), as is also the case for the Pleiades ( Figure 3). Furthermore, peaks in period increase in height with the modeled age of the clusters for the older clusters ( Figure   5g-l), whereas young clusters (Figure 5a-f) lack peaks. Middle-aged Pleiades ( Figure 3) having a very weak peak is consistent with this finding. Existence of peaks in П, weak or strong, is linked to the presence of old, white dwarfs ( Table 3).
The above findings are consistent with star production becoming less regular with time as gas is depleted and mass is stored in white dwarfs, which are not included in period measurements. Notably, peaks in the distribution of M-type stars are only associ-ated with the oldest clusters. These findings show that strong departures from steady-state production affect our fits, but not weak departures. Importantly, n is roughly 1 for all clusters, so for sub-solar dwarfs, star mass does not affect the law for spin decay.
In summary, fits in Figures 3, 5, 6 and 7 rule out n < ½ and n > 3/2, and thus unequivocally rule out Skumanich's rule (n = 3) for MKG stars. Crucially, rate laws with such values for n require that the appearances of the П-and ω-histograms differ substantially (details in Sections 4.4-4.6), contrary to the copious observations. Exponential decay, pro-posed by Terndrup et al. [8], best represents the available spin data on dwarf stars in in-dividual clusters. Moreover, production of small stars is continuous in the youngest clus-ters, but this condition becomes less appropriate for larger stars as the clusters evolve, whereas production of M-types departs from continuous behavior only for the oldest clus-ters.
to 0.57 in current use (e.g., [66]) makes little difference ( Figure 2a). Importantly, the square root shape for the I-trend in Figure 2a refers to the ad hoc B-V proportionality of Equation (4), rather than to any dependence on time.  [64] are not supported by larger, subsequent data sets [20,21]. Upper x-axis shows star mass. This graph has no connection with star age; (b) Comparison of data on periods near 0.3 M Sun to the I-and C-trends, calculated using Equations (5) and (6) with model cluster ages. Various symbols show individual clusters, see Table 1 for details. Pleiades stars are black cross symbols; (c) As above, but for 1 solar mass. These plots provide no support for Equations (5) and (6).

Model Assumptions
Gyrochronology Coeval (catastrophic) formation of dwarf stars occurs in any open cluster. The initial spin rate of stars is random. A physical law for spin decay is independent of star mass but depends on spin values. Star age can be inferred from a rate law that describes only the slowest dwarf stars (the I-trend), whereas the rate law describing the fast stars (the C-trend) is not relevant.
A t measurements all have n ~ 1.5 (Figure 5k), which is further from unity than valfor younger clusters. Peaks in their distributions cause comparatively poor fits. lacks small M stars. Fitting only the M-types of the other two clusters provides that yield n values closer to unity ( Figure 8; Table 3), as is also the case for the ( Figure 3). Furthermore, peaks in period increase in height with the modeled age sters for the older clusters (Figure 5g-l), whereas young clusters (Figure 5a-f) s. Middle-aged Pleiades ( Figure 3) having a very weak peak is consistent with ng. Existence of peaks in П, weak or strong, is linked to the presence of old, white able 3). above findings are consistent with star production becoming less regular with as is depleted and mass is stored in white dwarfs, which are not included in easurements. Notably, peaks in the distribution of M-type stars are only associthe oldest clusters. These findings show that strong departures from steadyduction affect our fits, but not weak departures. Importantly, n is roughly 1 for rs, so for sub-solar dwarfs, star mass does not affect the law for spin decay. mmary, fits in Figures 3, 5, 6 and 7 rule out n < ½ and n > 3/2, and thus unequivle out Skumanich's rule (n = 3) for MKG stars. Crucially, rate laws with such r n require that the appearances of the П-and ω-histograms differ substantially Sections 4.4-4.6), contrary to the copious observations. Exponential decay, pro-Terndrup et al. [8], best represents the available spin data on dwarf stars in inclusters. Moreover, production of small stars is continuous in the youngest clusthis condition becomes less appropriate for larger stars as the clusters evolve, production of M-types departs from continuous behavior only for the oldest clus-dependence is assumed for the I-trend, which rests on model ages of clusters. The I-and C-trends can be discerned visually from period vs. a temperature-color index scatter plot. That most stars fall between these two trends and are not described by either rate law is not germane.

Analytical inverse model (this paper)
Steady-state production approximates formation of sub-solar mass stars. A physical law for spin decay exists, and may or may not vary with star mass.
The "C-trend" was proposed [64] to describe fast-spinning stars: This equation is also ad hoc, but is uncommonly used to ascertain age. Problems in applying the gyrochronology model are recognized [67]. Concerns are as follows: • Equations (4)-(6) are fit to only a selected subset of the data, so data on many stars are ignored. From Figure 2, the strong curvature of the I-trend towards low B-V (large mass) differs among the various studies. Thus, the I-trend depends on the particular stars sampled. The axes of Figure 2 do not depict time, and so Skumanich's proposal of t e "gap" of Figure 2, in contrast to previous work (Section 2.5). In addition, avfits for ω and П to ascertain n reduces the importance of details.
lts of Fitting Individual Clusters dividual clusters define fits with n near unity, with a range of 0.83 to 1.55. Fits usters yield n between 1.0 and 1.2. The different values of n in the various fits 3, 5 and 6 are in part connected with zero bin size being impossible to attain. e important in fits of the period, since any slow rotators provide small ω. Adapplying the inverse relationship of ω with П converts the already broad peaks rams to very broad and weak peaks in ω-histograms, or adds a peak close to hich region of the histogram also depends on the maximum period measured. easons, angular velocity histograms adhere to power laws more closely than П-. Clearly, investigating pairs of histograms offsets imperfections in the dapes of two clusters differ substantially from those of the other 13, but these two exponential decay. Cepheus OB3b (Figure 5c) has some young stars possibly main sequence, and considerable data on comparatively large stars (types F e), see [73] (their Figure 13). Similarly, NGC 869 ( Figure 5e) consists of fast tators and is dominated by type-G stars [15]. lder clusters with much data (M37 and NGC 2516, Figure 5j,l) plus NGC 6811 easurements all have n ~ 1.5 (Figure 5k), which is further from unity than valfor younger clusters. Peaks in their distributions cause comparatively poor fits. acks small M stars. Fitting only the M-types of the other two clusters provides hat yield n values closer to unity ( Figure 8; Table 3), as is also the case for the igure 3). Furthermore, peaks in period increase in height with the modeled age ters for the older clusters ( Figure  bove findings are consistent with star production becoming less regular with s is depleted and mass is stored in white dwarfs, which are not included in surements. Notably, peaks in the distribution of M-type stars are only associthe oldest clusters. These findings show that strong departures from steadyction affect our fits, but not weak departures. Importantly, n is roughly 1 for , so for sub-solar dwarfs, star mass does not affect the law for spin decay. mary, fits in Figures 3, 5, 6 and 7 rule out n < ½ and n > 3/2, and thus unequivout Skumanich's rule (n = 3) for MKG stars. Crucially, rate laws with such n require that the appearances of the Пand ω-histograms differ substantially Sections 4.4-4.6), contrary to the copious observations. Exponential decay, proerndrup et al. [8], best represents the available spin data on dwarf stars in inusters. Moreover, production of small stars is continuous in the youngest clusis condition becomes less appropriate for larger stars as the clusters evolve, oduction of M-types departs from continuous behavior only for the oldest clus-cannot be evaluated from such scatter plots.

•
The C-trend also depends on the particular characteristics of the stars which were sampled ( Figure 2a), but the histograms do not, if a sufficient number of stars were measured ( Figure 3). • A mechanism for catastrophic (coeval) production of all stars in a cluster has not been explained. Counterevidence exists (e.g., [68,69]).

•
Evidence also exists that young stars start with spin that is fast (e.g., [5,6]), rather than from random rates. • Magnetic braking models fail to reproduce the distributions of Figure 2a when gyrochronology ages are utilized (e.g., [70]), even though this rate law purportedly describes this loss mechanism.
Four additional problems exist, but have not been discussed in the literature: • The postulated trends for fast and slow rotators cross near solar mass ( Figure 2a) for which B-V = 0.61 [37], showing gyrochronology models unrealistically require G2 stars to occupy a very narrow range of periods, and provide impossible, imaginary periods for more massive stars. • Gyrochronology calculations for 0.3M Sun (Equations (5) and (6)) based on model cluster ages only describe data on NGC 2516 ( Figure 2b). No solar mass stars in any of the clusters of Table 1 agree with this empirical formulation (Figure 2c).

•
The I-trend reflects two factors other than aging of the star: namely, that Kepler missions have a maximum reliable period, and that M-type stars are volumetrically much more abundant than more massive stars ( Figure 1a) and so statistically are likely to reveal a broader range of periods plus this upper limit. • Equations (4)-(6) require peculiar dimensions for their numerical constants, which results from color-period diagrams not depicting elapsed time.  [20]; (b) Data obtained during a recent mission [21], where the points are plotted at the top center of each bin; (c) A study focused on small, fast-spinning stars [22]. Histogram bins are shown as a symbol at the top of the bin in its middle. Black curve = least squares fit with the fit as indicated and n from the model developed in Section 3. If stated, the de-populated first bin was omitted in fitting.
To address known problems with gyrochronology models, recent studies propose more complicated multi-parameter models and "trends" in addition to C-and I-proposals (e.g., [11]). The age for our Sun is utilized, whereas other ages are gleaned from models. However, adding more equations and more parameters does not remedy the above flaws in gyrochronology, exemplified in Figure 2. These plots show that the slowest rotators detected depend on the maximum period that can be measured in finite duration missions (~25 days), whereas the fastest rotators detected (generally 0.1 to 0.5 days) depend on resolving periodicity from noise in the signal. Thus, the C-and I-trends largely reflect limitations of the measurements, not behavior of stars.
The remainder of our paper focuses on extracting a mass and time dependence from the measurements using an inverse approach. Skumanich's simple proposal for the timedependence of spin decay, which differs from the problematic gyrochronology formulations (Figure 2), is evaluated below.

Quantitative, Mathematical Analysis of Star Spin Histograms
Loss of spin of rotating stars with time is described by the general formula: where the exponent n can be non-integer. The proportionality constant α has units that depend on n. For the case of exponential decay (n = 1), τ = 1/α is the mean lifetime. If n = 0, then spin loss is controlled by factors other than spin, whereas n = 3 corresponds to Skumanich's "rule". Solutions to Equation (7) are: Symmetry 2021, 13, 1519 As in all previous studies, we assume that a star of a given mass starts its ignited life with a certain initial spin rate (ω 0 ) corresponding to a certain period (Π 0 ) and then slows down with time. Below, we use available data to decipher the best value of n, for each of many clusters, that controls the rate at which spin is lost on average, and then evaluate possible effects of stellar mass on the rate law. Consistency of the results is used to establish the physical law governing spin-down.

Histogram Shapes Do Not Support Coeval Star Production with Random Initial Spins
This section analytically tests a key assumption underlying previous models for star spin-down ( Table 2), namely that of simultaneous star formation in any given cluster. To probe the initial conditions, both ω and ω 0 must be considered as variables. Equations (8) and (9) lead to: Let us define a parameter G in terms of the entire histogram with a large number N total of stars: i.e., G = N total × bin width/range of ω 0 . Assuming random ω 0 and defining N as the number of stars in any one bin gives a relationship for incremental changes: Nδω = Gδω 0 . Rearranging and changing the increment to a derivative and then inserting Equation (10) gives: For constant G, differentiating Equation (11) and using Equation (10) provides: The quantity in square brackets is zero for n = 1, positive for n > 1, and negative for n < 1.
The slope ∂N/∂ω is thus zero for the exponential form, whereas slopes are positive for n > 1 and n < 0, but negative for 0 < n < 1. Figure 4 shows numerically generated synthetic histograms that represent coeval star production. First, at t = 0, random initial angular velocities are assumed to exist from 0 up to some limit (ω 0 = 0 to 1 in arbitrary units). How this synthetic cluster would appear at two finite times is illustrated for two cases: n = 1 (exponential decay) and n = 3 (Skumanich's proposal). The trends in the numerically generated histograms are fully consistent with Equations (11) and (12), confirming that these relationships are correct.
The synthetic shapes computed at finite times by assuming coeval production of stars with random initial spins ( Figure 4) are unlike histograms of real data (Figures 3 and 5). The synthetic histogram for n = 3 has a positive slope, which is completely at odds with the strong negative slopes observed for all large clusters that we examined (Figures 3 and 5). Figure 5 represents datasets with more than 243 stars. Similarly, the synthetic histogram for n = 1 has a null slope at all times (Figure 4), so that clusters of all ages should have "boxcar" histograms if production were coeval, which is likewise unsupported. Clearly, all cluster stars are not produced at the same time; the differences between the data in Figure 3 (or Figure 5) and Figure 4 are so extreme that no statistical test is needed. Therefore, we begin by assuming steady-state production rate as an endmember case, because this is the simplest remaining possibility. Sections 3.2 and 3.3 discuss the effect of relaxing this assumption on the rate law determined from the data. Figure 4, coupled with (12), demonstrates theoretically that the general shape of the histogram does not change with time if, and only if, decay is exponential. That similar histogram shapes are found for clusters of various presumed ages (Figures 3 and 5) points to exponential decay controlling spin-down, in agreement with the general discussion of Apostol [51]. This finding is irrespective of the details of fitting (Sections 3.5 and 4) and also of assuming steady-state star-production, which permits an analytical approach (Section 3.3).

Analytical Approach to Deciphering the Rate Law for Constant Production
To derive the number of spins in a counting bin, we begin with the general rate law of Equation (7). Recasting the general solution for ω and for n = 1 gives: Any bin on the histogram is bounded by two frequencies, ω a and ω b , and described by: For several reasons, we assume a constant rate (γ) of production for the stars, i.e., a constant rate at which stars become ignited and their periods can be measured. (1) This assumption provides the philosophical alternative to episodic production, akin to the historical viewpoints of Uniformitarianism and Catastrophism; (2) This assumption allows us to provide an analytical formula and to develop a method to quantitatively analyze paired histogram data. (3) Because the consensus view of catastrophic (coeval) production was ruled out above, we posit steady-state production as a feasible alternative. (4) Most of the data is on young clusters (Table 1), for which steady-state production of tiny stars is reasonable. (5) This assumption can be relaxed and others can be considered via numerical analysis in future work, as the database on measured periods expands. Below, episodic production is considered qualitatively.    Figure 3 for symbols and Pleiades data. All but two clusters (parts f,l) have >243 data points. Clusters are ordered from (a-l) using previous models for cluster age, listed in Table 1, along with data sources of spin. Panels (a,d,i,l) show the histograms in gray.
For constant production, the number of stars included in a given bin is the product of the production rate × the time interval, and thus: The width of the bin is ∆ω = ω a − ω b . From Equations (14) and (15), combining the constants γ and α to provide a new constant (β = α/γ) provides: We now expand the series assuming that that the bin size is small: Dropping the subscript b provides the remarkably simple solution: Although the above derivation considered n = 1, the same result holds for the special case of n = 1. Specifically, for n = 1, Utilizing a Taylor series and following the steps described above provides Equation (18). For periods, the starting point is Π n−1 = Π 0 n−1 + (2π) n−1 (n−1)αt for n = 1. Note that the power is (n−1). Solving for t and following the procedure as above leads to: For the special case of n = 1, the procedure described for angular velocity, analogous to (19), is followed, which results in Equation (20). In summary, if the stars constituting the histogram are produced at a constant rate, then for any power in the generalized rate law (Equation (7)), the number of stars in each bin of width ∆ω for the ω histogram is given by (18) while the number of stars for each bin of width ∆Π for the Π-histograms is given by (20). For both formulae, the constant β = α/γ describes the competition of spin loss to star production. For n = 1 only, both N(ω) and N(Π) inversely depend on their variables, due to the special properties of exponential functions. Hence, similar shapes will exist in the Πand ω-histograms if and only if n = 1, consistent with the general analysis of Apostol [51]. Section 4 uses Equations (18) and (20) to directly determine n from data on well-studied clusters.
Note that for small ω and small Π, the first bins in the histograms are under-populated because the null values for spin cannot exist in real datasets. Additionally, Equation (18) represents the bulk of the data for which ω is large compared to the increment ∆ω. Therefore, in comparing real data to the ideal curves, we omit the first bin, unless the cluster has a large number of stars, in which case under-population of the first bin is unimportant. For almost all clusters, we compare the data above~2∆ω and~2∆Π to Equations (18) and (20). Section 4.4 presents an alternative, numerical analysis which avoids potential problems with the variables being smaller than the increments and with an under-populated first bin. Virtually identical results for n are obtained.

What If Stars Are Produced with Different Spin Values or Decay at Different Rates?
Importantly, Equations (18) and (20), which describe the number of stars having given spins as a function of spin period or angular velocity, do not depend on the initial period Π 0 or initial angular velocity ω 0 , even though these are prominent factors in the rate law solutions, Equations (8) and (9). Instead, the only requirement is that the stars are produced at steady state, with a distribution of initial spin values that likewise does not change with time. Therefore, for steady state production, the starting point (i.e., the specific initial spin for any given star) is immaterial to the histogram shape. Our analysis only presumes that once a star has ignited and is spinning, that it slows down according to a certain rate law. The possibility that this rate law may be mass-dependent is considered below.
The exponent n in the rate law (Equation (7)) is revealed by the shape of the histogram curves. It is only the histogram shape that is important to ascertaining n. Neither the interval, decay constant, production rate, nor initial spin rate are particularly influential.
These findings are important for several reasons. First, our analysis is independent of processes possible for very young stars, such as disk formation or star contraction. Likewise, very old stars, which can experience significant mass loss, also have negligible influence. Second, different clusters can be compared together, as long as these are governed by the same rate law. Third, this behavior permits comparison of stars with different starting points; the only requirement is that stars in any given histogram pair are produced at a steady state, and afterward follow the same rate law ( Table 2).

Effect of Episodic Star Production
If stars in a cluster were produced over a narrow time interval and also with a restricted range in period, then the histograms would display a strong peak at some finite value of Π (or ω = 2π/Π). The range of initial distributions would control the breadth and shape of such a peak. Such behavior, if sufficiently strong, could interfere with ascertaining the rate law. However, only for the two oldest clusters do strong peaks exist (Figure 5k,l). A moderate peak does exist for the old Hyades cluster, but few stars were measured ( Figure 6). In all other well-populated histograms (Figures 3 and 5a-j), peaks are weak, if they exist and are associated with the larger stars. This is clearly the case for M34, which also has few data points (Figure 6b). Nevertheless, most fits are compatible with n = 1. For the sub-solar stars emphasized here, the histograms are consistent with nearly constant production, except possibly for the most ancient of clusters. Details are below. Figure 6. Paired histograms of clusters with~125 members. Symbols as in Figure 5: (a) Hyades, for which data from [28] are highly accurate and cover a wide range of periods; (b) M34 data from [26], which sampled K-and G-types, with few periods exceeding 15 days.

Bin Size Has a Negligible Effect When Analyzing Paired Histograms
To test whether our choice of bin size affects our results, paired histograms for 4 large clusters with different model ages were fit to a wide range of bin sizes. Large datasets were considered because they permit the test to be conducted over a large range in bin size. The maximum number of bins allowed in our plotting program is 400, which permitted reaching very tiny bin sizes (e.g., 0.1 days). Regarding wide bins, we typically stopped at bin sizes of 5 days or 5 radians per day, which yielded only 10 wide bins (e.g.,~5 day widths). Figure 7 shows very consistent trends in the powers derived from fitting, as in Figures 3, 5 and 6. The average value of n from histogram pairs is always near unity.  Table 1 were explored as these allow examination of the largest range in bin size. Least squares fits are shown. All evaluations point to n~1: (a) Data from [31]; (b) Data from [21]; (c) Data from [19]; (d) Data from [29].
As bin size approaches zero, the value of n for the angular velocity histogram decreases whereas the value of n for the period histograms increases. For the power fits shown, n varies from~0 to~2. The extracted values of n occupy a narrower range for all bin sizes. When the bin size is small, but not so small as to distort the shapes, n for each histogram type is near unity. Even the old cluster Praesepe, which has a strong peak (Figure 5k), behaves in this manner. Hence, the average n-value for the two fits is used to represent n. The bin size chosen has little effect on the determination of n for the rate law, as demonstrated by the crossover in Figure 7 near 0.2 to 1 widths. Our histograms typically used widths of 0.5 (Figures 3, 5 and 6).

Analysis of Histogram Data on Dwarf Star Spin
The primary goal is to determine the integer that best describes the underlying physical law (7). The secondary goal is to ascertain whether the same spin-down law holds for different star types (i.e., different mass ranges). The focus is on open clusters, most of which are active sites of star formation, as indicated by presence of nebulosity, dust, and large, short-lived B-types.

Individual Clusters
More than 243 stars have been measured in each of 11 clusters (Table 1). The next largest samplings (NGC 2547 and NGC 6811) have 176 and 171 data points, respectively. Figure 5 shows power law fits for the bin count as a function of both ω and Π for 12 of these 13 clusters, which are ordered from youngest to oldest in terms of previously modeled cluster ages. Figure 3 shows similarly analyzed Pleiades data, which would place after M50 (Figure 5g). Figure 6 shows two clusters with 120-129 data points. The value of n for each cluster (Table 3) is ascertained by averaging the two fits to ω and Π (Equations (8) and (9)).

Fits Are Slightly Affected by Star Mass and Sample Size
The number of measurements weakly affects the fits, showing that statistical significance is obtained for associations with~200 stars. However, another factor exists: NCG 2547 and NGC 6811 similarly have~173 data points, but only the former is fit well: the latter, which has only three M-type stars, is fit poorly, due to the presence of a peak in the period distribution (Figure 5l). That this peak is connected with larger stars, rather than with the number of measurements, is supported by comparing M34 to the Hyades ( Figure 6), which have 120 and 129 data points, respectively: As above, M34 has a peak, no data on M-types, and does not conform well to a simple power law. Star mass is explored further below.
Notably, the Hyades measurements are highly accurate and yield tightly constrained fits. This finding is consistent with random errors being reduced either as accuracy improves or as sample size increases [71].

Bin Size and Depopulation of Bins near the Origin Have Little Effect
Behavior at short and long ω or Π is not important to the analysis. In particular, we found that omitting the first bin does not significantly change the fit (i.e., the variation in n is~0.1). Figure 7 supports this deduction, since these fits omitted no stars, and further shows that bin size has a negligible effect. Determination of n does not hinge on the values at high and low ω or Π, mainly because the result for n depends on the curvature of the histogram. That is, fits mostly depend on bin counts for intermediate values of the variables. Hence, our results depend on the behavior of the majority of the stars, which are those in the "gap" of Figure 2, in contrast to previous work (Section 2.5). In addition, averaging the fits for ω and Π to ascertain n reduces the importance of details.

Results of Fitting Individual Clusters
All individual clusters define fits with n near unity, with a range of 0.83 to 1.55. Fits for most clusters yield n between 1.0 and 1.2. The different values of n in the various fits of Figures 3, 5 and 6 are in part connected with zero bin size being impossible to attain. This is more important in fits of the period, since any slow rotators provide small ω. Additionally, applying the inverse relationship of ω with Π converts the already broad peaks in Π-histograms to very broad and weak peaks in ω-histograms, or adds a peak close to ω = 0, in which region of the histogram also depends on the maximum period measured. For these reasons, angular velocity histograms adhere to power laws more closely than Π-histograms. Clearly, investigating pairs of histograms offsets imperfections in the datasets.
Star types of two clusters differ substantially from those of the other 13, but these two still follow exponential decay. Cepheus OB3b (Figure 5c) has some young stars possibly not on the main sequence, and considerable data on comparatively large stars (types F and G-type), see [73] (their Figure 13). Similarly, NGC 869 ( Figure 5e) consists of fast (young) rotators and is dominated by type-G stars [15].
Two older clusters with much data (M37 and NGC 2516, Figure 5j,l) plus NGC 6811 with 171 measurements all have n~1.5 (Figure 5k), which is further from unity than values found for younger clusters. Peaks in their distributions cause comparatively poor fits. NGC 681 lacks small M stars. Fitting only the M-types of the other two clusters provides better fits that yield n values closer to unity ( Figure 9; Table 3), as is also the case for the Pleiades (Figure 3). Furthermore, peaks in period increase in height with the modeled age of the clusters for the older clusters (Figure 5g-l), whereas young clusters (Figure 5a-f) lack peaks. Middle-aged Pleiades ( Figure 3) having a very weak peak is consistent with this finding. Existence of peaks in Π, weak or strong, is linked to the presence of old, white dwarfs ( Table 3).
The above findings are consistent with star production becoming less regular with time as gas is depleted and mass is stored in white dwarfs, which are not included in period measurements. Notably, peaks in the distribution of M-type stars are only associated with the oldest clusters. These findings show that strong departures from steady-state production affect our fits, but not weak departures. Importantly, n is roughly 1 for all clusters, so for sub-solar dwarfs, star mass does not affect the law for spin decay.
In summary, fits in Figure 3, Figure 5, Figure 6, Figure 7 rule out n < Symmetry 2021, 13, x FOR PEER REVIEW shows that bin size has a negligible effect. Determination of n does not hinge on t at high and low ω or П, mainly because the result for n depends on the curvat histogram. That is, fits mostly depend on bin counts for intermediate values of bles. Hence, our results depend on the behavior of the majority of the stars, w those in the "gap" of Figure 2, in contrast to previous work (Section 2.5). In add eraging the fits for ω and П to ascertain n reduces the importance of details.

Results of Fitting Individual Clusters
All individual clusters define fits with n near unity, with a range of 0.83 to for most clusters yield n between 1.0 and 1.2. The different values of n in the va of Figures 3, 5 and 6 are in part connected with zero bin size being impossible This is more important in fits of the period, since any slow rotators provide sma ditionally, applying the inverse relationship of ω with П converts the already bro in П-histograms to very broad and weak peaks in ω-histograms, or adds a pea ω = 0, in which region of the histogram also depends on the maximum period m For these reasons, angular velocity histograms adhere to power laws more closel histograms. Clearly, investigating pairs of histograms offsets imperfections i tasets.
Star types of two clusters differ substantially from those of the other 13, but still follow exponential decay. Cepheus OB3b (Figure 5c) has some young stars not on the main sequence, and considerable data on comparatively large stars and G-type), see [73] (their Figure 13). Similarly, NGC 869 ( Figure 5e) consis (young) rotators and is dominated by type-G stars [15].
Two older clusters with much data (M37 and NGC 2516, Figure 5j,l) plus N with 171 measurements all have n ~ 1.5 (Figure 5k), which is further from unity ues found for younger clusters. Peaks in their distributions cause comparatively NGC 681 lacks small M stars. Fitting only the M-types of the other two clusters better fits that yield n values closer to unity (Figure 8; Table 3), as is also the ca Pleiades ( Figure 3). Furthermore, peaks in period increase in height with the mo of the clusters for the older clusters (Figure 5g Table 3).
The above findings are consistent with star production becoming less reg time as gas is depleted and mass is stored in white dwarfs, which are not in period measurements. Notably, peaks in the distribution of M-type stars are on ated with the oldest clusters. These findings show that strong departures from state production affect our fits, but not weak departures. Importantly, n is roug all clusters, so for sub-solar dwarfs, star mass does not affect the law for spin de In summary, fits in Figures 3, 5, 6 and 7 rule out n < ½ and n > 3/2, and thus ocally rule out Skumanich's rule (n = 3) for MKG stars. Crucially, rate laws w values for n require that the appearances of the П-and ω-histograms differ sub (details in Sections 4.4-4.6), contrary to the copious observations. Exponential de posed by Terndrup et al. [8], best represents the available spin data on dwarf s dividual clusters. Moreover, production of small stars is continuous in the youn ters, but this condition becomes less appropriate for larger stars as the cluster whereas production of M-types departs from continuous behavior only for the ol ters. and n > 3/2, and thus unequivocally rule out Skumanich's rule (n = 3) for MKG stars. Crucially, rate laws with such values for n require that the appearances of the Πand ω-histograms differ substantially (details in Section 3.2), contrary to the copious observations. Exponential decay, proposed by Terndrup et al. [8], best represents the available spin data on dwarf stars in individual clusters. Moreover, production of small stars is continuous in the youngest clusters, but this condition becomes less appropriate for larger stars as the clusters evolve, whereas production of M-types departs from continuous behavior only for the oldest clusters.

Associations and the Solar Neighborhood
Star associations are less well-bound gravitationally than clusters. The Scorpious-Centaurus OB association [31] includes many young, large B and O stars. Many pre-main sequence stars were sampled: some of which appear to be accreting. Because the sample size (162) is marginally large for a statistical analysis and the data are limited to periods <9 days, the fits are less robust than those discussed above and similar to open clusters with inadequate sampling. Nonetheless, n near 1.1 is inferred (Figure 8a). A larger dataset on similarly young stars (Figure 8b) provides n = 1, within uncertainty. A long period for one star in µ Tau [32] tightly constrains n. The contrasting case of the Solar Neighborhood, which was culled to fast rotators with Π < 4 days [33], provides a reasonable fit due to the large number of points (481). This culling removes the old, slow stars and apparently reduces the effect of episodic production.
Thus, two different associations and the Solar Neighborhood behave equivalently to the open clusters, and corroborate that spin decays exponentially with time ( Table 3). The period measurements in these recent studies are considered more accurate (e.g., as in the Hyades [28]), which improves reliability of statistical analysis.

Ground-Based Measurements and a Comparison to Kepler Missions
Ground-based datasets [35,36] are accurate, can include slow rotators, and also provide n near unity ( Figure 10). However, the stars sampled need not be related. Thus, observing M-types in the Southern Hemisphere (Figure 10a) provides a generalized sample of the Solar Neighborhood and reveals a peak in Π near 100 days, which Kepler missions cannot probe (Figure 8c). This long-period peak is consistent with these M-types being on the order of 4.5 Ga old, like the Sun. Interestingly, the distribution is bimodal, consisting of fast and slow rotating M-types, yet histogram analysis still indicates exponential decay. Ground-based measurements differ from the remotely probed clusters and associations, where the curvature is controlled by the intermediate values, since bins for slow and fast rotators are underpopulated.  Figure 5: (a) Pre-main sequence stars of the Scorpious-Centaurus OB association from [31]; (b) A wide range of masses in the µ Tau association, which includes a B2 star and pre-main sequence stars [32]. Arrow points to the slow rotator; (c) Fast rotators in the Solar Neighborhood [33]. Adding the Sun makes little difference to the fits. Figure 9. Fits for M-type stars. Symbols described in Figure 5: (a) M37 (data from [27]) has a peak at 12 days, but is fit reasonably well; (b) NGC 2516 (data from [24]) has more M-types, in number and proportion, but this distribution also have a peak; (c) Small field stars close to the Sun: data from [34]. Note that the angular momentum panel is on a log-log scale. These stars are not related, need not be produced continuously, and their distribution is poorly represented by a power law. A very large sampling of nearby M-type field stars [34] from Kepler missions (Figure 9c) shows a strong peak at 30 days. As above, this dataset yields 1< n < 1.5, due to the peak being present. Plus, longer periods exist, but are not observed because the durations of the missions setting a cap on maximum Π. In contrast, ground-based observations from the Northern hemisphere (Figure 10b) include more distant stars of various masses, but have only a few M-types. For larger stars, a constant production rate is less likely. The peak at 40 days influences the Π fit, particularly as long periods are not present in this 100 member sampling. However, on average n = 1.
The datasets in Figures 9c, 8c and 10 thus represent nearby stars in complementary manners. All support exponential spin decay holding for stars with mass near solar and below, with perceived variations of n from unity and accuracy of fitting resulting from non-continuous production and the limitations of the various datasets.

Histograms of Stars from Open Clusters with Similar Mass
Because mass appears to govern the initial spin rate [74,75], mass may affect the decay constant α. Hence, we next analyze histograms where stars with similar mass in the various clusters are grouped together. The aggregates for M and K use data available prior to 2019, which were processed similarly: we used 362 measurements of Irwin et al. [26] for NGC 2516; 108 stars from Delorme et al. [76] plus Douglas et al. [53] for Hyades; and 71 stars from Meibom et al. [77] for NGC 6811. For the G stars, all studies of Table 1 were used to provide a larger dataset, exclusive of duplicate measurements.

Aggregated Data on M-Types with a Numerical Test
The distribution of M-type stars (i.e., with 0.03 < M/M Sun < 0.48) from diverse clusters ( Figure 11) exhibits the same strongly skewed pattern as observed for most individual clusters (Figures 3, 5, 6 and 9). Integer values for n were assumed in least squares fitting. For type-M stars, n = 1 provides the best fit. Moreover, Figure 11 clearly shows that only the exponential rate law provides similar skewed shapes for both the ωand Π-histograms. Skumanich's rule is clearly refuted by its unrealized prediction of a positive slope in the Π-histogram.  (17) and (19), as labeled. Counts for bins near the origin not included in fitting are indicated. Exponential decay (n = 1) provides the best fit for both visualizations.
To avoid assuming that ∆ω is smaller than the smallest measured ω, we now evaluate Equations (18) and (20) numerically. Figure 12a presents the ω histogram as N vs. 1/ω to place emphasis on the well-populated bins. By plotting against 1/ω, the under-populated initial bins (open symbols, which have values near the bin size in Figure 11) can be excluded from the fit, and young stars are still emphasized. Similarly, plotting N vs. 1/Π emphasizes old stars and well-populated bins (Figure 12b). Figure 11 on type-M stars shows that the analytical relationship for periods of (20) is valid for all but the first two bins. Examining Equation (7) suggests that N should be initially a constant, which is consistent with the fits. Numerical analysis (Figure 12) corroborates that n = 1, exponential decay, holds over the lifetime of M-type stars.  (18) and (20). Dots = count in Figure 11 for M-types: (a) Analysis of ω-histograms using bin size of 0.5 radians/day. By plotting count against 1/ω, the under-populated initial bins (open symbols) can be excluded from the fit, and young stars are emphasized; (b) Analysis of Π-histograms using a bin size of 0.5 days. Bin #1 is off scale with N = 8 and 1/Π = 11 d −1 . This approach does not assume that ω (or Π) is large compared to bin size. The linear trends show that spin attenuation is exponential (n = 1).
The patterns shown in Figures 11 and 12 are not caused by consolidating cluster data for two reasons.

•
Recent studies utilizing the Kepler 2 mission [21,29] largely describe M-types, and their histograms are likewise compatible only with n near 1 (Figures 3a and 5k).

•
Periods of M stars in individual clusters behave similarly to the cluster as a whole (Figure 9a,b). Only M37 shows a peak at Π = 12 d, which may represent episodic production. M37 is considered to be old and thus is most likely to depart from continuous production, since gas and dust are consumed with time.
In short, exponential spin-down describes M-type stars for any and all clusters, as well as M-types in star associations and in the Solar Neighborhood.

Type-K Stars
The most commonly sampled stars in missions previous to Kepler 2 are type-K's. Histograms of K stars ( Figure 13) have a secondary peak at Π = 7 d. The corresponding broad peak expected at 0.9 radian d −1 in ω-histograms is not resolvable. Despite the moderate peak existing at Π = 7 d, the overall best fit for both histograms on K-type stars is exponential decay (n~1). Cases n < shows that bin size has a negligible effect. Determination of n does not hinge on the values at high and low ω or П, mainly because the result for n depends on the curvature of the histogram. That is, fits mostly depend on bin counts for intermediate values of the variables. Hence, our results depend on the behavior of the majority of the stars, which are those in the "gap" of Figure 2, in contrast to previous work (Section 2.5). In addition, averaging the fits for ω and П to ascertain n reduces the importance of details.

Results of Fitting Individual Clusters
All individual clusters define fits with n near unity, with a range of 0.83 to 1.55. Fits for most clusters yield n between 1.0 and 1.2. The different values of n in the various fits of Figures 3, 5 and 6 are in part connected with zero bin size being impossible to attain. This is more important in fits of the period, since any slow rotators provide small ω. Additionally, applying the inverse relationship of ω with П converts the already broad peaks in П-histograms to very broad and weak peaks in ω-histograms, or adds a peak close to ω = 0, in which region of the histogram also depends on the maximum period measured. For these reasons, angular velocity histograms adhere to power laws more closely than Пhistograms. Clearly, investigating pairs of histograms offsets imperfections in the datasets.
Star types of two clusters differ substantially from those of the other 13, but these two still follow exponential decay. Cepheus OB3b (Figure 5c) has some young stars possibly not on the main sequence, and considerable data on comparatively large stars (types F and G-type), see [73] (their Figure 13). Similarly, NGC 869 (Figure 5e) consists of fast (young) rotators and is dominated by type-G stars [15].
Two older clusters with much data (M37 and NGC 2516, Figure 5j,l) plus NGC 6811 with 171 measurements all have n ~ 1.5 (Figure 5k), which is further from unity than values found for younger clusters. Peaks in their distributions cause comparatively poor fits. NGC 681 lacks small M stars. Fitting only the M-types of the other two clusters provides better fits that yield n values closer to unity (Figure 8; Table 3), as is also the case for the Pleiades (Figure 3). Furthermore, peaks in period increase in height with the modeled age of the clusters for the older clusters (Figure 5g-l), whereas young clusters (Figure 5a-f) lack peaks. Middle-aged Pleiades ( Figure 3) having a very weak peak is consistent with this finding. Existence of peaks in П, weak or strong, is linked to the presence of old, white dwarfs ( Table 3).
The above findings are consistent with star production becoming less regular with time as gas is depleted and mass is stored in white dwarfs, which are not included in period measurements. Notably, peaks in the distribution of M-type stars are only associated with the oldest clusters. These findings show that strong departures from steadystate production affect our fits, but not weak departures. Importantly, n is roughly 1 for all clusters, so for sub-solar dwarfs, star mass does not affect the law for spin decay.
In summary, fits in Figures 3, 5, 6 and 7 rule out n < ½ and n > 3/2, and thus unequivocally rule out Skumanich's rule (n = 3) for MKG stars. Crucially, rate laws with such values for n require that the appearances of the П-and ω-histograms differ substantially (details in Sections 4.4-4.6), contrary to the copious observations. Exponential decay, proposed by Terndrup et al. [8], best represents the available spin data on dwarf stars in individual clusters. Moreover, production of small stars is continuous in the youngest clusters, but this condition becomes less appropriate for larger stars as the clusters evolve, whereas production of M-types departs from continuous behavior only for the oldest clusters. and n > 2 are ruled out. To further explore the connection of histogram shape to mass, we constructed a histogram of small K-type stars (0.48 to 0.65 M Sun ; not shown). The appearance of the histogram pair for small type-K stars is intermediate to those of the M stars ( Figure 11) and K stars (Figure 13), but with the weak and very broad peak at 7 days being split into two even weaker peaks at 5 and 10 days. For small K stars, decay is close to exponential, but as mass increases, a secondary peak develops which suggests a component of episodic production. Such episodic production causes departures of the measurements from n~1 behavior, since our mathematical analysis (Section 3) assumed constant production.
Because episodic production increases with mass, the largest K-types were not included in Figure 13. Instead, the largest K stars were grouped with type-G, discussed further below.  (18) and (20), as labeled. Counts for bins near the origin not included in fitting are indicated. No curves provide optimal fits due to episodic production, particularly for M37. Nonetheless, the two plots constrain~1/2 < n < 2.

Type-G Stars
The clusters of Table 1 have varying numbers of type-G stars. To examine masses near solar and a large number of stars, Figure 14 includes the largest K's and the smallest F's in the clusters. The range of masses above and below solar is balanced in this manner, although very few F0-types were sampled. Despite production of these G stars having some episodic character, exponential decay of spin is supported.  (18) and (20), as labeled. Counts for bins near the origin not included in fitting are indicated. Coefficients are shown for the exponential case (red curves and text. None of the fits are optimal due to episodic production. However, n ≥ 2 can be ruled out based on period data. For angular velocity, n from 1 to 2 provides the best fits. The best fits to a floating power are: n = 0.74 for Π and n = 1.37 for ω. Episodicity in production of large stars in the young cluster NGC 869 exists, but is slight. Unlike other young clusters of Table 1, solar mass stars constitute the largest fraction of measurements, and a statistically significant number of F-types were measured ( Figure 15). Exponential decay describes the subgroups, recognizing that fitting period data with a peak and few data points typically gives n slightly larger than unity, and <1 for ω-histograms (not shown).

Laws Governing Spin Decay for MKG, and Probably F, Stars
Both ω-histograms and Π-histograms for aggregated type-M and -K stars of all clusters in Table 1 combined are strongly skewed. This behavior also exists for type-M stars in individual clusters with many data points, except for old clusters (cf. Figures 3,5,7 and 9), which have a peak in period that indicates some episodic production. The companion peak in ω-histograms is not resolved due to depopulation of the initial bins and/or the inverse relationship. Existence of a peak in Π correlates with large stellar mass (Figure 14), suggesting that star production either becomes more episodic as mass increases or that production of the large stars slows with time. The former explanation is consistent with older cluster populations having a greater fraction of the larger type-G stars.
Small M stars do not show evidence of episodic production, except for the oldest cluster and for members of the Solar Neighborhood, which is circa 4.5 Ga, based on meteorite ages. For this reason, the best fits are found for clusters with many measurements of smaller stars (Figures 3 and 5; Table 3).
The better match existing for small stars is attributed to their continuous production, consistent with their formation requiring only tiny volumes of gas and dust. The match for individual clusters, plus aggregated data on type-M and small K (Figures 11 and 12) unequivocally rules out the exponents of n < 1/2 and n > 1.8 in Equations (18) and (20), and points to exponential decay of spin for the smallest dwarf stars. Even for older clusters or those with larger stars (Figure 5e,j,k), exponential decay is suggested by the data whereas n = 3 is ruled out. Skumanich's 1972 proposal of n = 3 is inconsistent with the ω-histograms and completely incompatible with the Π-histograms.

Fundamental Problems with the Magnetic Braking Hypothesis
Skumanich's rule has found support in models of magnetic braking associated with ionized stellar winds (e.g., [78,79]). Regardless of whether flux or field strengths are considered, the principle underlying magnetic braking is that a moving charge experiences a sideways force (F) in a magnetic field (B) according to Lorentz's force law: where v is the velocity vector and q is charge. For two independent reasons, the net force (or torque) should be negligible on average: • Our Sun's magnetic field reverses polarity every 11 years, which is accompanied by its sunspot cycle. Forces for the two polarities would act in opposite directions (21), thus providing alternate braking and acceleration. Hence, magnetic braking cancels magnetic acceleration over very short-term of 22 yr. This 22 yr cycle is unchanged over~400 years of observations [80,81].

•
The sideways force depends on charge q, not on the mass of the charged particles. The stellar wind should be composed of roughly equal numbers of negatively and positively charged particles. Otherwise, the star would develop a net charge with time which then would preferentially retain the oppositely charged particles and restore charge balance. No net torques would arise.
Magnetic braking is not a feasible cause of spin-down if charge is conserved or if stars have reversals. Both behaviors are likely.

Viscous Dissipation as the Mechanism of Spin-Down
Viscous stress in a Newtonian fluid creates a force, F = Aµ∂v x /∂y where A = area, µ = the coefficient of viscosity, and flow is along the x direction [82]. The equatorial region dominates, so y~R, v x~R ω, and A is some function of R 2 , where R = radius. Using torque = R × F = I∂ω/∂t, where I = JMR 2 and J is the reduced moment of inertia, and ∂ω/∂R = 0 for a solid body gives: which has the same form as Equation (7) for n = 1 and f = α. The function f combines the mass dependence of all variables. Shear existing in stars is evident in the long-known latitudinal dependence of the rotation rate of our Sun, and also in the super-rotation of its core [83]. Differential rotation is also known in other stars: notably, part of the K2 mission was intended to explore differential rotation (e.g., [84]). Thus, viscous dissipation is a straightforward and logical explanation for progressive loss of stellar spin.

New Isochrones from Histogram Analysis
Quantifying spin-down using (7) requires both an initial angular velocity (or initial period) and a time constant for at least two stellar masses. We set limits, as follows.
The maximum initial spin rate is set by of the Virial Theorem, which relates gravitational potential energy (U g ) to rotational kinetic energy (R.E.): Equation (23) restricts the spin at formation for homogeneous density as: shows that bin size has a negligible effect. Determination of n does not hinge on the values at high and low ω or П, mainly because the result for n depends on the curvature of the histogram. That is, fits mostly depend on bin counts for intermediate values of the variables. Hence, our results depend on the behavior of the majority of the stars, which are those in the "gap" of Figure 2, in contrast to previous work (Section 2.5). In addition, averaging the fits for ω and П to ascertain n reduces the importance of details.

Results of Fitting Individual Clusters
All individual clusters define fits with n near unity, with a range of 0.83 to 1.55. Fits for most clusters yield n between 1.0 and 1.2. The different values of n in the various fits of Figures 3, 5 and 6 are in part connected with zero bin size being impossible to attain. This is more important in fits of the period, since any slow rotators provide small ω. Additionally, applying the inverse relationship of ω with П converts the already broad peaks in П-histograms to very broad and weak peaks in ω-histograms, or adds a peak close to ω = 0, in which region of the histogram also depends on the maximum period measured. For these reasons, angular velocity histograms adhere to power laws more closely than Пhistograms. Clearly, investigating pairs of histograms offsets imperfections in the datasets.
Star types of two clusters differ substantially from those of the other 13, but these two still follow exponential decay. Cepheus OB3b (Figure 5c) has some young stars possibly not on the main sequence, and considerable data on comparatively large stars (types F and G-type), see [73] (their Figure 13). Similarly, NGC 869 ( Figure 5e) consists of fast (young) rotators and is dominated by type-G stars [15].
Two older clusters with much data (M37 and NGC 2516, Figure 5j,l) plus NGC 6811 with 171 measurements all have n ~ 1.5 (Figure 5k), which is further from unity than values found for younger clusters. Peaks in their distributions cause comparatively poor fits. NGC 681 lacks small M stars. Fitting only the M-types of the other two clusters provides better fits that yield n values closer to unity (Figure 8; Table 3), as is also the case for the Pleiades (Figure 3). Furthermore, peaks in period increase in height with the modeled age of the clusters for the older clusters (Figure 5g-l), whereas young clusters (Figure 5a-f) lack peaks. Middle-aged Pleiades ( Figure 3) having a very weak peak is consistent with this finding. Existence of peaks in П, weak or strong, is linked to the presence of old, white dwarfs ( Table 3).
The above findings are consistent with star production becoming less regular with time as gas is depleted and mass is stored in white dwarfs, which are not included in period measurements. Notably, peaks in the distribution of M-type stars are only associated with the oldest clusters. These findings show that strong departures from steadystate production affect our fits, but not weak departures. Importantly, n is roughly 1 for all clusters, so for sub-solar dwarfs, star mass does not affect the law for spin decay.
In summary, fits in Figures 3, 5, 6 and 7 rule out n < ½ and n > 3/2, and thus unequivocally rule out Skumanich's rule (n = 3) for MKG stars. Crucially, rate laws with such values for n require that the appearances of the П-and ω-histograms differ substantially (details in Sections 4.4-4.6), contrary to the copious observations. Exponential decay, proposed by Terndrup et al. [8], best represents the available spin data on dwarf stars in individual clusters. Moreover, production of small stars is continuous in the youngest clusters, but this condition becomes less appropriate for larger stars as the clusters evolve, whereas production of M-types departs from continuous behavior only for the oldest clusters.
, (24) where G N is Newton's gravitational constant [75]. Although internal densification has some effect (see [85] for examples including polytropes), constant density suffices for the present purposes. Equation (24) provides Π 0 = 0.095 d for M = M Sun which agrees with the fastest spinning type-G stars found in four different clusters ( Figure 16). The t = 0 line must include the origin, since for null mass, the period is zero initially and remains zero for all time.
Using the origin provides a linear equation for initial period of Π 0 = 0.095M/M Sun days ( Figure 16) which agrees with the fastest-spinning stars known [5,6]. The initial spin is also confirmed by Solar System trends [74,75]. Our t = 0 trend is further supported by the observation that many spin at slightly slower rates, denoted as the breakup velocity = (MG/R) Symmetry 2021, 13, x FOR PEER REVIEW shows that bin size has a negligible effect. Determination of n does n at high and low ω or П, mainly because the result for n depends o histogram. That is, fits mostly depend on bin counts for intermedia bles. Hence, our results depend on the behavior of the majority o those in the "gap" of Figure 2, in contrast to previous work (Sectio eraging the fits for ω and П to ascertain n reduces the importance o

Results of Fitting Individual Clusters
All individual clusters define fits with n near unity, with a ran for most clusters yield n between 1.0 and 1.2. The different values of Figures 3, 5 and 6 are in part connected with zero bin size being This is more important in fits of the period, since any slow rotators ditionally, applying the inverse relationship of ω with П converts th in П-histograms to very broad and weak peaks in ω-histograms, o ω = 0, in which region of the histogram also depends on the maxim For these reasons, angular velocity histograms adhere to power law histograms. Clearly, investigating pairs of histograms offsets imp tasets.
Star types of two clusters differ substantially from those of the o still follow exponential decay. Cepheus OB3b (Figure 5c) has some not on the main sequence, and considerable data on comparative and G-type), see [73] (their Figure 13). Similarly, NGC 869 (Figu (young) rotators and is dominated by type-G stars [15].
Two older clusters with much data (M37 and NGC 2516, Figu with 171 measurements all have n ~ 1.5 (Figure 5k), which is furthe ues found for younger clusters. Peaks in their distributions cause co NGC 681 lacks small M stars. Fitting only the M-types of the other better fits that yield n values closer to unity ( Figure 8; Table 3), as Pleiades ( Figure 3). Furthermore, peaks in period increase in height of the clusters for the older clusters (Figure 5g-l), whereas young lack peaks. Middle-aged Pleiades ( Figure 3) having a very weak p this finding. Existence of peaks in П, weak or strong, is linked to the dwarfs ( Table 3).
The above findings are consistent with star production becom time as gas is depleted and mass is stored in white dwarfs, whic period measurements. Notably, peaks in the distribution of M-type ated with the oldest clusters. These findings show that strong de state production affect our fits, but not weak departures. Importan all clusters, so for sub-solar dwarfs, star mass does not affect the law In summary, fits in Figures 3, 5, 6 and 7 rule out n < ½ and n > 3 ocally rule out Skumanich's rule (n = 3) for MKG stars. Crucially values for n require that the appearances of the П-and ω-histogram (details in Sections 4.4-4.6), contrary to the copious observations. Ex posed by Terndrup et al. [8], best represents the available spin data dividual clusters. Moreover, production of small stars is continuou ters, but this condition becomes less appropriate for larger stars a whereas production of M-types departs from continuous behavior o ters. , e.g., [86]. Black symbols are reference points: square in circle = the Sun and a "young" Herbig object [6]. Solid lines = proposed isochrones, ages as labeled, where the upper panels consider p = 1 in (24) while the lower use p = 0.57. Dotted line = initial spin from the Virial theorem (see text) which is confirmed by the Herbig object and the fastest known O-star, with M = 22M sun and a period of 0.88 days [5]. Heavy red arrow on the upper axis indicates a typical observational cutoff. Stars in any open cluster have many ages: (a) Open clusters considered to be old, plus 29 stars in Coma Berenices (red-brown dots) [72] and 30 slow rotators in M67 (light pink dots) [60]; (b) Clusters with moderate model ages. Most stars are~3 Ga old, which a substantial portion of young stars, and some evidence for stellar demise in NGC 2516, which lacks fast rotators; (c) Clusters with young ages, plus data on 20 fast rotators in the young cluster IC4665 (pale blue-gray triangles) [57]. The average age is near 4 Ga, and no evidence exists for stellar demise; (d) Comparison of one large dataset for each of old, moderate, and young model ages to field M-types (grey dots). Note that the x-axis is shifted from that in the other panels. The trend for field stars [34] indicates that periods of M-types in the missions are limited to 40-60 days. Large green dots = data on 3 Be stars in Pleiades, which are consistent with the young ages of smaller members.
The only information available to establish lifetime is the Sun's current age (4.54 Ga) and period (24.5 d at the equator), which provide τ = 0.82 × 10 9 yr. Consequently, we suggest the following dependence of period on age: Combining the fit for density in Figure 1b with (24) suggests that power p = 0.57 (Figure 16c,d). This value assumes uniform density inside a star. We also consider p = 1 (Figure 16a,b). When M is near solar, p is immaterial. Thus, both cases for Equation (25) represent stars similar to our type-G2 Sun. As a first approximation, we assume (25) holds for main sequence M to F stars ( Figure 16). Isochrones from (25) are thus defined by t = τ ln(Π/Π 0 ).

Implications of Our Isochrones on Ages of Small Stars
Comparing Equation (25) to period data ( Figure 16) shows that ages of main sequence dwarf stars in clusters mostly range from~2 to 4.5 Ga. This finding is influenced by Kepler mission measurements of G-types being impacted by an upper limit of 25 to 32 days for reliable periods [60]. Ground-based measurements reach much larger Π of 225 days [35] and indicate ages approaching 8 Ga for M-types if p = 1 or 7 Ga if p = 0.57, whereas the maximum ages for K and G stars observed from Mount Wilson (e.g., [36]) are less,~5.5 Ga for either choice of p, as masses are near solar. Nearby stars are thus only a little older than the Sun. This finding is consistent with progressively larger stars having shorter life spans, as is well-known. Thus, our isochrones are supported by ground-based data.
Data on old open clusters are also consistent with the isochrones (Figure 16a). First, young stars are present only in Praesepe, which has the most measurements. Moreover, nebular gas is also absent, evolved stars exist, and the number of stars is few (Table 3). Finally, trends of period with mass for the oldest stars in evolved clusters which also have many measurements (Hyades, Praesepe, and NGC 6811) are parallel to the trend for ground-based data on the oldest, slowest stars (Figure 16a).
These parallel trends for the slowest star at any given mass describe destabilization of main sequence stars at old age. Destabilization at a certain period is explained qualitatively by fluid mechanics. Patterns of fluid flow in large rotating bodies depend on the relative effectiveness of inertial and Coriolis forces [82] so decreasing spin diminishes the latter, possibly triggering events which move the star off of the main sequence. Density and temperature both monotonically changing with mass ( Figure 1b) are consistent with destabilization depending on mass. Clusters with moderate model ages, all of which have white dwarfs (Table 1), also possess many young stars and rather few old stars (Figure 16b). Except for NGC 2516, a pile-up near the Sun's position is not seen. This cluster also has fewer young stars than the others with similar model ages.
Clusters with young model ages and nebulae have many young stars, but also have ancient members of 4 Ga (Figure 16c). The trends of individual clusters are nearly parallel to the isochrons, particularly if p = 0.57, and so stellar demise is not evident in the period statistics. None of these clusters are described as having white dwarfs. Greater ages in Figure 16c than 16b is consistent with their population becoming older, but very few members have left the main sequence trend. Figure 16d compares well-measured clusters with different ages to nearby field Mtypes and to Equation (25) with p = 0.57. The maximum age of 5.8 Ga for field M's is greater than that of Praesepe members, with a parallel trend. Because the number of measurements in a restricted mass range is high (1570), McQuillan et al.'s [34] data reveal the maximum reliable periods for the Kepler missions. This maximum of~50 days for 0.4 M Sun is consistent with the reliability study of solar masses (<32 days [60]) and larger stars having shorter lifetimes. Ground-based measurements have large ages of 7 Ga for M-types (Figure 16a), corroborating that the maximum periods in the oldest clusters from Kepler mission datasets are underestimated. Maximum ages for the clusters in Figure 16b,c are not affected, because very few members reach the limits of measurable periods.
Because star ages are most accurate near the solar mass, Table 3 lists mean and maximum ages for the G-types. Because very fast rotators are under-represented in the measurements, Table 3 includes the minimum age for all stars in each group.

Implications of Our Isochrones on Ages and Evolution of Clusters
Cluster ages are estimated from trends in type-G ages (Table 3) as follows. The cluster must be at least as old as its oldest stars. Even the youngest, actively producing clusters have type-G stars with average ages near 3.3 Ga, as well as some old stars approaching 5 Ga. However, measurements cannot reliably resolve long periods and thus large ages. The ages of the oldest stars of any type in any cluster correlate roughly with age of its youngest star (Figure 17a). Because the upper age of stars involves the former demise of ancient population, the cluster age is set by adding 3 Ga to the age of the average G-type. The additive constant is based on the mean of 3.3 Ga for all clusters being close to half of the age of the oldest type-G star (5.3 Ga in the Mt Wilson group: Table 3).  Table 3, with the addition of the maximum for M-types using p = 1: (a) Dependence of various star ages on the age of the youngest star in each group or cluster. Convergence of the average (black dots) and oldest (open squares) type-G stars at~4.5 Ga is consistent with constant production until nebula gas runs out and G stars die. M stars are typically the oldest members of any cluster (diamonds). The green line is the average age of the G star in the cluster plus 3 Ga (see text); Agreement of this line with the oldest stars provides the best estimate for the cluster ages; (b) Dependence of star ages on the distance of the group from the Sun. Trends are controlled by NGC 869, whose database is limited to large M's, and is the most distant cluster; (c) Dependence of star ages on previous models for cluster age (Table 1). These trends rest solely on the 38 reliable measurements of M67 [60].
The resulting age for the cluster is slightly greater than the oldest stars, which are almost always M-types (Figure 17a). This offset is consistent with periods being limited by the duration of the Kepler campaigns.
The clusters have similar ages (Table 3) because all are located in essentially the same part of the immense Milky Way ( Table 1). The oldest clusters have little gas and few remaining stars, while the youngest clusters have nebulosity or molecular clouds and lack white dwarfs. These findings are supported by the weak correlations of star age with distance from the Sun (Figure 17b) and the scatter plot of star age with modeled cluster age (Figure 17c). Previous inferences of cluster age (Table 1) depend on modeled ages of their younger, large main sequence stars.

Conclusions
This paper provides an inverse method to extract the physical law governing the loss of star spin with time from histograms of spin period and its reciprocal, angular velocity. All available data on large open clusters that each have >120 measured light curves provide histogram pairs with similar shapes, indicating that exponential decay best describes spin-down of stars smaller than three Solar masses. Superimposed peaks in the ωand Π-histograms are attributed to episodic production. Skumanich's formula and the previous gyrochronology model are contra-indicated. Our deduction of exponential spin-down is independent of the initial spin for the stars and is robust against the measurements being underpopulated for very short and very long periods.
If coeval production of all stars in open clusters is assumed, our forward models yield histogram shapes that are unlike any observed, regardless of what rate law is considered. The replacement of a coeval theory for star production with one assuming steady state production of small stars is parallel with the catastrophism vs. uniformitarianism dilemma faced by pioneer geologists (ca. 1785 to 1830), with the latter explanation ultimately prevailing for numerous processes.
Importantly, period measurements are made of small stars on the main sequence, which describes these dwarfs over most of their lifetime, while the star is stable. The exponential rate law likely applies to all stable stars of any size. The associated time constant probably depends on mass, but ascertaining this is beyond the scope of the present report.
Diverse phenomena are controlled by exponential attenuation, such as absorption of light by matter, radioactive decay, heat flow (Newton's law of cooling), capacitor discharge, or vibrational damping. Based on exponential spin-down, we constructed isochrones using the known age of the Sun and its estimated initial spin ( Figure 16). Our isochrones are most accurate for stars with ages of several billion years and masses near or slightly below solar. Figure 16d suggests that type-B6 stars can attain much greater ages than previously considered. Additional statistical analysis of data on large stars would be useful, and is possible using Doppler measurements. Nonetheless, the trend suggests a resolution to the puzzle as to how a young, large B8-O0 star could orbit an old black hole [87]. Specifically, the orbiter could be far older than assumed.

Astrophysical Implications
The Sun spins slowly compared to those near-solar-mass stars in open clusters and nearby regions of the Milky Way (Table 1 lists the many studies). This fact, combined with histogram analyses of the present report, points to it not being middle-aged as commonly assumed (Figures 3-9, 14 and 16). The Sun numbers among the oldest type-G stars ( Figure 17). From the Mt. Wilson study, maximum ages of G-types increase from 5 to 5.6 Ga as mass decreases (Figure 16a). Our finding, combined with the 4.5 Ga time-span on Earth for humanoids to develop, offers a resolution to Fermi's paradox. Namely, the requirement for star longevity coupled with planet habitability significantly reduces the frequency of possible civilizations.

General Applicability
Histograms are commonly used to represent a key variable of a population. Our mathematical formulation can be applied to other phenomena. The variable does not need not be angular momentum, or even to have a meaningful reciprocal (e.g., period), but can be any physical attribute. The limitation is that the primary variable must be described by a power law similar to Equation (7), which covers many possibilities, since n can be non-integer. Note that the independent variable does not need to be time.