HOMERUN: Multi-Domain Models in Science & Sports

• HOMERUN in astrophysics models emission-line spectra with a multi-cloud photoionization framework using non-negative least squares for accurate, stratified ISM analyses.
• In particle physics, the HOMER method infers effective Lund string fragmentation functions from hadronic observables, achieving reconstruction errors as low as 1–5% under varying complexities.
• Baseball analytics leverages HOMERUN to normalize home-run statistics via era-adjusted and hierarchical models, redefining player rankings through refined probabilistic assessments.

Searching arXiv for papers on “HOMERUN” to ground the article in published work. HOMERUN appears in current technical literature in several distinct senses. In astrophysics it denotes “Highly Optimized Multi-cloud Emission-line Ratios Using photo-ionizatioN,” a multi-cloud photoionisation framework that models an observed spectrum as a non-negative combination of CLOUDY single-cloud calculations and is used to infer metallicities, abundance ratios, density structure, ionisation parameters, attenuation, and line-to-mass conversion factors (Marconi et al., 2024). A related but separate usage appears in particle physics through the HOMER or H method, which addresses a restricted inverse problem of hadronization by extracting an effective Lund string fragmentation function from hadron-level observables (Assi et al., 7 Mar 2025). In baseball analytics and sports physics, home runs are the direct object of statistical and dynamical analysis, including era detrending, hierarchical prediction, run-to-win valuation, and ballistic reconstruction (Petersen et al., 2010).

1. Terminological scope and research domains

The same string, “HOMERUN,” therefore spans at least three research regimes: photoionisation modelling, hadronization inference, and quantitative baseball analysis. The usages are not historically continuous and do not share a common formalism.

Usage	Field	Core object
HOMERUN	Astrophysical photoionisation modelling	Multi-cloud fitting of emission-line spectra
HOMER or H method	Hadronization / event generators	Extraction of $f(z)$ from observable hadronic data
Home run analytics	Baseball statistics and physics	Era-adjusted counts, component rates, run value, trajectory dynamics

The astrophysical usage is the most explicit acronymic one in the supplied literature, and it is the only case where HOMERUN is formally expanded as “Highly Optimized Multi-cloud Emission-line Ratios Using photo-ionizatioN” (Marconi et al., 2024). The hadronization paper explicitly states that “HOMERUN” does not appear in that manuscript; the operative label there is the HOMER or H method, with “HOMERUN” only a conceptual extension in the accompanying explanation (Assi et al., 7 Mar 2025). In baseball work, by contrast, the term is literal rather than acronymic.

2. HOMERUN as a multi-cloud photoionisation framework

In its astrophysical sense, HOMERUN replaces the single-zone approximation with a non-negative linear combination of constant-density single-cloud models. In the notation used across the method papers and applications, an observed line flux or luminosity is approximated as

$$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$

or equivalently

$$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$

The weights are solved by non-negative least squares, and the fit quality is evaluated with a $\chi^2$-like loss function built from observed line fluxes and adopted uncertainties (Marconi et al., 2024).

The original methodological paper constructs large CLOUDY grids spanning ionisation parameter $\log U$, density $\log N_H$, stellar ionising continua from BPASS, dust/no-dust variants, and gas-phase metallicity expressed through $12+\log(\mathrm{O/H})$. A central refinement is that N/O and S/O are allowed to vary through scaling parameters applied to nitrogen and sulphur lines. This flexibility is not ancillary: when nitrogen and sulphur are not rescaled, the minimum $\chi^2$ can worsen by factors up to $100$, and inferred metallicities can shift by up to $\pm 0.4$ dex (Marconi et al., 2024).

The framework is explicitly designed to fit broad ionization ranges simultaneously. In nearby H II regions, all lines are reproduced with an accuracy better than $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$0, including the difficult combination of [O I] $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$1, [O II] $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$2, [O III] $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$3, [S II] $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$4, and [S III] $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$5 (Marconi et al., 2024). This is the point at which HOMERUN most clearly departs from constant-pressure or single-cloud models: the improvement is not only in goodness of fit, but in the ability to reproduce low-, intermediate-, and high-ionisation zones within one self-consistent spectral model.

The same architecture also provides derived physical quantities beyond line fitting. The weighted cloud distribution yields effective densities and ionisation parameters, while the best-fitting grid determines $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$6, abundance ratios such as N/O and S/O, and, in later AGN applications, line-luminosity-to-mass conversion factors for ionized outflows (Marconcini et al., 11 Dec 2025). The acceptable-model envelope is defined through

$$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$7

which is used repeatedly to quote uncertainties on fitted quantities (Marconi et al., 2024).

3. Astrophysical applications: stratified ISM, high-redshift galaxies, and AGN outflows

HOMERUN’s later applications are dominated by systems in which single-zone modelling fails. In the $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$8 galaxy GHZ2, multi-zone photoionisation modelling shows that the spectrum cannot be reproduced by single-density spectro-photometric models. The fitted solutions require a strongly stratified ISM in which low-/intermediate-density gas and high-density regions with $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$9 coexist. The line spectrum is consistent either with a composite star-formation plus AGN scenario or with star formation in a combination of radiation-bounded and matter-bounded regions. Purely radiation-bounded stellar models fail to reproduce the observed He II emission, making an additional hard ionising component unavoidable. The same fit yields $$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$0, corresponding to an N/O ratio about $$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$1 times solar (Castellano et al., 9 Dec 2025).

In local AGN, HOMERUN is increasingly coupled to the 3D kinematic code MOKA$$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$2. In NGC 1068, the integrated optical-plus-mid-IR line set is fitted by two AGN-ionized components: a dust-rich phase with $$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$3, $$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$4, and $$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$5, and a dust-poor phase with $$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$6, $$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$7, and $$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$8. The best-fit metallicity is $$L_{l,\mathrm{mod}} = \sum_{j=1}^{m} w_j\,L_{l,\mathrm{mod},j}.$$9. When combined with MOKA$\chi^2$0, the analysis finds that [O IV] traces an outflow $\chi^2$1 faster than [O III], and that the mid-IR-revealed dusty component carries a significantly larger ionized-gas mass than optical lines alone would imply (Marconcini et al., 11 Dec 2025).

In NGC 1365, HOMERUN fits $\chi^2$2 optical-to-mid-IR emission lines in a $\chi^2$3 aperture using a dust-free AGN component and a dusty star-formation component. The best metallicity is $\chi^2$4, with AGN and star-forming gas characterized by $\chi^2$5, $\chi^2$6 and $\chi^2$7, $\chi^2$8, respectively. Coupled to MOKA$\chi^2$9, the modelling yields consistent outflow masses from [Ne V] and [O III], whereas classical single-zone methods differ by orders of magnitude between those tracers (Ceci et al., 10 Jul 2025).

The broader diagnostic lesson is stated most sharply in “One cloud is not enough.” Across three case studies at $\log U$0, even a small fraction of unresolved high-density clumps can contribute more than half of the observed flux of auroral lines while contributing negligibly to standard optical density tracers. Under those conditions, $\log U$1-method metallicities can be underestimated by $\log U$2 dex. The same study shows that disagreements between UV- and optical-based N/O estimates do not necessarily imply chemical inhomogeneities: in RXCJ2248-ID they can arise from ionisation and density structure alone, whereas in the Sunburst Arc the data favour genuine chemical stratification, with an N-enriched component coexisting with a chemically normal one (Moreschini et al., 13 Jan 2026).

4. The HOMER method in hadronization and the restricted inverse problem

A separate usage, closely adjacent in name but different in domain, appears in particle-physics work on hadronization. The updated HOMER method addresses a restricted inverse problem: assuming the Lund string model is correct, infer an effective fragmentation function from hadron-level observables by reweighting simulated string breaks (Assi et al., 7 Mar 2025).

The underlying fragmentation function is the symmetric Lund form

$\log U$3

with $\log U$4 the light-cone momentum fraction and $\log U$5 the transverse mass. The method proceeds in three stages. First, a classifier distinguishes synthetic “data” from “simulation” using event-level observables, with event weights

$\log U$6

Second, per-string-break weights are inferred so that their product over a chain reproduces the event weights in expectation. Third, new simulations are reweighted to reconstruct an effective $\log U$7 (Assi et al., 7 Mar 2025).

The technical innovation required for gluon-rich strings is an observable-space smearing over neighbouring events,

$\log U$8

which approximates the average over invisible fragmentation histories. The classifier uses $\log U$9 high-level observables, including event-shape variables $\log N_H$0, $\log N_H$1, $\log N_H$2, $\log N_H$3, and $\log N_H$4, multiplicities $\log N_H$5 and $\log N_H$6, and moments of $\log N_H$7 with $\log N_H$8 (Assi et al., 7 Mar 2025).

The paper studies four increasingly complex $\log N_H$9 scenarios at $12+\log(\mathrm{O/H})$0 GeV: a fixed $12+\log(\mathrm{O/H})$1 string with no gluons; a fixed $12+\log(\mathrm{O/H})$2 string; a variable one-gluon string; and a full parton-shower string with an unrestricted number of gluons. Across these cases, the extracted fragmentation function remains accurate, with degradation described as relatively modest: the reconstruction error rises from the $12+\log(\mathrm{O/H})$3–$12+\log(\mathrm{O/H})$4 level in simpler cases to about $12+\log(\mathrm{O/H})$5 in the full $12+\log(\mathrm{O/H})$6 case (Assi et al., 7 Mar 2025). This suggests that global event shapes and multiplicities retain substantial information about longitudinal string fragmentation even after shower complexity is introduced.

5. Home runs in sabermetric modelling and statistical inference

In baseball research, home runs are analysed both as outcomes to be normalized across eras and as component rates embedded in larger batting and team-level models. One line of work detrends seasonal and career totals by comparing a player’s home-run rate to the league-wide weighted seasonal baseline. For player $12+\log(\mathrm{O/H})$7 in season $12+\log(\mathrm{O/H})$8,

$12+\log(\mathrm{O/H})$9

with $\chi^2$0 the number of home runs and $\chi^2$1 the number of at-bats. Detrended seasonal home runs are then

$\chi^2$2

and detrended career totals are sums over seasons. This procedure removes league-wide shifts due to factors such as performance-enhancing drugs, expansion, equipment, and rule changes without assigning causal responsibility to any single factor (Petersen et al., 2010).

Applied to career distributions from 1920–2009, raw and detrended home-run totals have essentially the same right-skewed functional form, well approximated by a Gamma distribution or truncated power law. For raw career HR, the maximum-likelihood parameters are $\chi^2$3 and $\chi^2$4; for detrended HR they are $\chi^2$5 and $\chi^2$6. The average ratio $\chi^2$7 is $\chi^2$8 and bimodal, reflecting inflation of early-era totals and deflation of late-era totals. The paper also quantifies the “steroids era” shift: average league home-run prowess over 1994–2009 is about $\chi^2$9, versus $100$0 over 1978–1993, a statistically significant $100$1 increase (Petersen et al., 2010).

This era adjustment reorders historical rankings. Babe Ruth’s raw total of $100$2 HR ranks third, behind Barry Bonds ($100$3) and Hank Aaron ($100$4), but his detrended total of $100$5 ranks first by a wide margin; Barry Bonds falls to eighth with $100$6, and Hank Aaron to fifth with $100$7. The same framework defines objective benchmarks for extraordinary careers through extreme-value thresholds of the fitted Gamma law, giving about $100$8 detrended HR as a top-tail benchmark at the $100$9 level (Petersen et al., 2010).

A second line of work embeds home runs in hierarchical component models of batting. Batting average is decomposed as

$\pm 0.4$0

where

$\pm 0.4$1

Because the multinomial likelihood factorizes into three independent binomial likelihoods, strikeout probability, conditional home-run probability, and hit-in-play probability can be estimated separately with exchangeable beta random-effects models (Albert, 2015).

For the 2011 season, using players with at least $\pm 0.4$2 AB, the fitted home-run hyperparameters are $\pm 0.4$3 and $\pm 0.4$4, implying $\pm 0.4$5 across players. The posterior-mean estimator

$\pm 0.4$6

implements partial pooling, and expected future home-run totals follow

$\pm 0.4$7

In a $\pm 0.4$8-year prediction contest for batting average, the component model is generally superior to a single beta-binomial model, especially in 1963–1980 and 1995–2012 (Albert, 2015).

A third line of work translates runs into wins. The generalized Pythagorean framework models runs scored and allowed as independent three-parameter Weibull variables with possibly different shape parameters, estimates their parameters by Method of Moments, and computes winning percentage as $\pm 0.4$9 through a numerical integral when $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$00. Over the last $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$01 MLB seasons, this Method-of-Moments model narrowly outperforms Pythag$$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$02 and clearly outperforms least-squares Weibull fits in mean squared error of predicted wins (Almeida et al., 2023). A plausible implication is that the value of a home-run hitter is not only a shift in mean runs scored but also, potentially, a change in scoring variance, which the unequal-shape Weibull model can accommodate (Almeida et al., 2023).

6. Home-run trajectory physics and the Mantle reconstruction problem

Home runs are also a computational-physics problem. The reconstruction of Mickey Mantle’s 22 May 1963 Yankee Stadium shot treats the trajectory as a two-dimensional boundary-value problem with gravity, quadratic drag, and Magnus lift. The equations of motion use a state vector $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$03, with drag coefficient

$$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$04

and Magnus coefficient

$$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$05

The ball is required to start at $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$06, $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$07 and reach $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$08, $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$09, corresponding to the Yankee Stadium facade (Warren, 2024).

Using historical weather, the authors adopt $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$10, $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$11, $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$12, and $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$13. With no wind and a “still rising” constraint at impact, the minimum exit speed consistent with the geometry occurs for $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$14 and $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$15, yielding

$$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$16

This is above any Statcast-recorded home run. When the “still rising” constraint is relaxed so the ball may be slightly falling at impact, and a tailwind is included, the minimum becomes

$$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$17

with an unobstructed range of about $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$18 or $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$19 (Warren, 2024).

The comparison set is explicit. Since 2015, Statcast’s hardest-hit home run is Giancarlo Stanton’s $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$20 ($$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$21), and the longest is Nomar Mazara’s $$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$22 ($$F_i^{\rm obs} \approx \sum_j w_j F_{i,j}^{\rm model}, \qquad w_j \ge 0,$$23). Mantle’s reconstructed shot is therefore either beyond modern records under strict eyewitness constraints or, under more conservative assumptions, within the uppermost modern exit-velocity range while still exceeding current distance records (Warren, 2024). The paper’s final conclusion is deliberately qualified: there is no single definitive answer, but Mantle’s 1963 home run remains physically compatible with the upper edge of human capability and compares favorably with the most powerful home runs of the Statcast era (Warren, 2024).