{"id":"caf5f7d2-bf92-43a4-87a0-57f05992f69d","arxiv_id":"2411.14346","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Row-wise standardized daily load profiles lie on a mathematical sphere, PCA projects them onto an arc, and a principal curve plus a von Mises-Fisher distribution provides ordering and generative sampling.","lead":"Daily electrical load profiles from medium voltage transformers, when normalized and compressed with PCA, arrange themselves along a curved arc on a three-dimensional sphere. The authors use this geometric picture to detect faulty meters, order profiles by similarity, and generate synthetic load profiles.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's explicit simplification to a spherical (rather than ellipsoidal) projection in Section II-A is load-bearing and untested; because the PCA components are highly unequal in variance, the fitted sphere, marginal outlier thresholds, and single-kappa VMF generative model are likely…","rationale":"The reader's weakest_assumption already identifies the spherical-versus-ellipsoid simplification; I agree. The additional value here is pinning down why it is load-bearing: the anisotropy is not a minor perturbation but is visible in Table III (phi std about 50 degrees, theta std about 9 degrees) and follows directly from the unequal eigenvalues in Table II. Every contribution, outlier detection, ordering, and generative VMF modeling, depends on the geometry being spherical and isotropic. The ordering via principal curve is more robust because it uses the arc itself, but the outlier and generative claims require the isotropy that is assumed, not shown. Since the existing verdict is CONDITIONAL with moderate confidence, this concern strengthens the conditions but does not overturn the verdict; the paper can address it by fitting anisotropic directional distributions and reporting whether conclusions change.","tokens_in":16104,"tokens_out":5426,"duration_ms":58852,"concrete_test":"Fit by maximum likelihood both a 3D isotropic von Mises-Fisher distribution and a 5-parameter Kent (anisotropic Fisher-Bingham) distribution to the centered 3D PCA coordinates of the clean Municipality 1 data. If the Kent model wins by an information criterion such as BIC with a difference larger than 10, or a likelihood-ratio test rejects the isotropic restriction, then the spherical assumption underlying Eqs. (15), (16)-(18), and (23) is not supported; the paper would need to replace the isotropic sphere with an anisotropic directional model for outlier thresholds and generative sampling.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Row-wise standardization (Eqs. 2-8) indeed places every profile on a unit hypersphere, and the first three PCs retain roughly 90% variance, so the projected points lie close to the unit ball's surface. But the claim that the projected distribution is 'spherical-shaped' rather than ellipsoidal is assumed 'for simplicity' immediately after Eq. (14), and the paper never tests it. Table II shows PC1 explains about 69% of variance versus about 20% for PC2 and about 4% for PC3; after centering, the 3D coordinates have strong anisotropy. Table III quantifies this: the azimuth phi has std about 50 degrees, while the polar angle theta has std about 9 degrees; the data form an arc, i.e., a one-dimensional curve on the sphere, not a rotationally symmetric cloud. This matters for every downstream step. The sphere fit (Eq. 15) minimizes radial residuals and can absorb the anisotropy into a mean radius, but the residual angular structure remains. Outlier detection then fits independent von Mises marginals for phi and theta and a skew-normal for r and flags points outside 95% CIs; this ignores the strong phi-theta dependence along the arc, so it can miss off-arc points and flag legitimate arc endpoints. The generative model (Eq. 23) uses a single concentration kappa, which is rotationally symmetric around mu=f(s); it cannot represent the anisotropic local spread of a real arc without different concentrations for tangential versus normal directions. The paper takes kappa=7.1 from the theta marginal, which does not represent the much wider phi spread. The empirical validation in Section V-D shows only marginal distributions for four municipalities; Municipality 4 has only 81.8% CEV for three PCs and 16.9% outliers, so the spherical-shell approximation is weakest exactly where the paper claims generalization. Thus the central geometric assumption is plausible but unverified, and the claimed unification depends on it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that daily medium-voltage load profiles, after row-wise z-score standardization and PCA projection into three dimensions, lie on a sphere with an arc-shaped latent distribution. On this basis it proposes a unified pipeline: fitting a sphere and spherical-coordinate marginals for outlier detection, fitting a principal curve through the arc to order profiles and define clusters, and using a von Mises-Fisher distribution centered on the principal curve for generative modeling. The claims are illustrated on load-profile data from four Dutch municipalities.","tokens_in":1960,"tokens_out":2034,"duration_ms":74841,"significance":"If the spherical and arc-shaped structure were established, the proposed framework would be a genuinely useful unification of visualization, outlier detection, ordering, and generation for load-profile data, with a small parameter count and stable PCA projections. The paper is also honest about several simplifications and provides a reproducible principal-curve implementation. However, the main empirical claim is weakened in two ways. First, the fact that row-wise standardized profiles lie on a hypersphere is a mathematical consequence of the normalization, not an empirical discovery; the empirical content lies entirely in the shape of the three-dimensional PCA projection. Second, the paper explicitly assumes, rather than tests, that this projection is spherical rather than ellipsoidal, even though its own Tables II and III show strong anisotropy. The proposed outlier and generative models inherit this untested geometric assumption. The work is therefore a promising demonstration, but its central latent-structure claim needs substantially more quantitative support.","major_comments":[{"comment":"Equation (6) as printed is incorrect: the numerator should be the sum of squared deviations, not the sum of deviations, and the denominator should be that same sum divided by D. As written, the numerator is zero by the definition of the mean. The conclusion ||p_i|| = sqrt(D) is correct after this correction, but the proof needs to be fixed. In addition, the statement that standardization 'unveils' a latent spherical distribution overstates what the algebra shows: any row-wise z-score normalized matrix lies on a hypersphere by construction. The empirical claim is specifically that the three-dimensional PCA projection has spherical and arc-shaped structure, and that claim needs to be separated from the algebraic identity.","section":"Section II-A, Eqs. (4)-(7)"},{"comment":"Equation (15) is not a well-posed optimization problem as written. The objective is linear in rho with a negative coefficient, so for any center c one can make the objective arbitrarily small by taking rho arbitrarily large; no finite minimizer exists. Presumably the intended objective is a least-squares sphere fit such as sum_i ( ||z_i - c||^2 - rho )^2, where rho = r^2. Because the fitted center c and radius rho define the spherical-coordinate transformation used for all downstream outlier detection and generative modeling, this error is load-bearing and must be corrected.","section":"Section III, Eq. (15)"},{"comment":"The assumption that the three-dimensional PCA projection is 'spherical-shaped instead of an ellipsoid' is never tested, and the paper's own summary statistics indicate strong anisotropy. Table II shows PC1 explaining roughly 69% of variance versus about 20% and 4% for PC2 and PC3, and Table III reports a standard deviation of about 50 degrees for the azimuthal angle versus about 9 degrees for the polar angle. These numbers are more consistent with an arc-like, effectively one-dimensional distribution than with a rotationally symmetric spherical cloud. The paper should provide a direct test of sphericity, for example by comparing the sphere-fit residuals of Eq. (15) with an ellipsoid fit, or by evaluating the eigenvalue spectrum of the local covariance structure. Until such a test is provided, the use of isotropic directional-statistics tools is not justified.","section":"Section II-A, after Eq. (14); Tables II and III"},{"comment":"The outlier detection rule fits independent von Mises marginals for phi and theta and a skew-normal marginal for r, then flags points outside the 95% confidence interval of each marginal. This ignores the strong dependence between phi and theta along the arc: points at the arc endpoints can be flagged as outliers merely because one marginal is extreme, while points that deviate from the arc in the joint (phi, theta) plane can pass all three marginal tests. A joint directional outlier criterion, for example a density threshold under a fitted VMF or Kent distribution, would be more appropriate. In addition, the 95% confidence level is chosen without justification, and no quantitative evaluation of outlier-detection performance against labeled anomalies is provided.","section":"Section III, Fig. 3"},{"comment":"The generative model uses a single concentration parameter kappa = 7.1, taken from the theta marginal, for a von Mises-Fisher distribution that is rotationally symmetric around its mean direction f(s). Along a curved, anisotropic arc the local spread is generally different in the tangential direction and the two normal directions, so a single kappa cannot reproduce the observed distribution. The paper should either fit an anisotropic directional distribution (such as a Kent distribution) or demonstrate empirically that the generated points reproduce the joint distribution of the projected data, rather than only the marginal summaries in Table I.","section":"Section V-C, Eqs. (23)-(24)"},{"comment":"The cross-municipality validation is descriptive rather than quantitative. The paper shows that the first three eigenprofiles look similar across municipalities and that the spherical-coordinate marginals have comparable moments, but it does not test whether the three-dimensional projections in the other municipalities are spherical, nor whether the arc ordering is present. The case of Municipality 4 is particularly concerning: its three-component cumulative explained variance is only 81.8% and its outlier fraction is 16.9%, in contrast to the other municipalities. The claim that the spherical structure with latent ordering is a general property of MV load profiles requires a quantitative test, for example a goodness-of-fit measure for the spherical model or a comparison of principal-curve reconstruction error across all four datasets.","section":"Section V-D, Tables II and III, Fig. 8"}],"minor_comments":[{"comment":"The title contains a stray space in 'Voltage' ('V oltage'), and the Conclusion refers to the 'von Fisher-Misses distribution'; this should be 'von Mises-Fisher distribution'.","section":"Title and Conclusion"},{"comment":"There is a typo 'hypershpere' that should read 'hypersphere'.","section":"Section II-A, paragraph after Eq. (14)"},{"comment":"Equation (20a) appears to omit the square in the Gaussian exponent: it should be exp(-1/2 (t - mu)^2), not exp(-1/2 (t - mu)).","section":"Section IV-A, Eq. (20a)"},{"comment":"The phrase 'it is possible can identify outlier profiles' is ungrammatical; it should be 'it is possible to identify'.","section":"Section I, paragraph after 'It is critical to notice'"},{"comment":"The first two paragraphs of Section II-A repeat nearly the same statement about standardization of daily MV load profiles; one of the two paragraphs should be removed or merged.","section":"Section II-A, first two paragraphs"}],"recommendation":"major_revision","confidential_remarks":"The paper's central geometric claim is a mixture of a trivial algebraic fact (row-wise standardization places points on a hypersphere) and an interesting but untested empirical observation (the 3D PCA projection is spherical and arc-shaped). The current manuscript does not provide the quantitative evidence needed to separate these two components. The most serious technical issue is the ill-posed sphere-fitting objective in Eq. (15), which must be corrected before the results can be reproduced. If the authors add a proper sphericity test, replace the marginal outlier rule with a joint directional method, fit an anisotropic generative model, and provide quantitative cross-dataset validation, the paper could become a solid application-oriented contribution. In its present form, the claims outrun the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is a decent applied contribution that overclaims its theoretical novelty. The observation that row-wise z-scored profiles lie on a hypersphere is a mathematical identity—Eq. (6) is just the definition of variance. What's actually new is the empirical finding that the PCA projection of these standardized profiles forms a thin, arc-shaped band on the sphere, and that a principal curve along that arc gives a meaningful ordering of load profiles, from commercial to residential with PV. The four-municipality validation is a real plus; the eigenprofiles and spherical coordinate distributions are genuinely similar across datasets.\n\nThe soft spots are real but not fatal. The paper explicitly assumes the 3D projection is spherical rather than ellipsoidal, and never tests the anisotropy. Table III shows the azimuth has about 50 degrees std while the polar angle has about 9 degrees; that is not a rotationally symmetric cloud. This matters for the VMF generative model, which uses a single concentration parameter and on which the paper pins continuous cluster mixing. It also matters for outlier detection, where independent marginals on phi and theta will miss off-arc points or flag legitimate arc endpoints. Still, the method demonstrably catches a known faulty-meter pattern (absolute-value recording) that would be hard to catch otherwise, so it is not without value.\n\nThe bigger issue is the gap between the 'unification' claim and the actual workflow. The outlier threshold is hand-chosen 95% confidence intervals, the cluster bins are set post hoc to match AggHC, and the principal curve fitting has its own tuning. That is normal for applied work, but it undercuts the promise of a single, parameter-light representation.\n\nI would send this to review with a request for major revision: test the sphericity assumption explicitly (e.g., compare against an ellipsoidal or kappa-anisotropic model), report outlier detection precision/recall on labeled faults, and release code and data so the ordering and generative claims are reproducible. The authors are honest about the 'for simplicity' assumption, which makes me think they can address these.\n\nThis paper is for DNO data analysts and load profiling researchers. It will not reshape the field, but it is a useful, citable example of using directional statistics in power systems. A serious referee should engage with it.\n\nRecommendation: accept for peer review, with major revisions.","headline":"A useful practical workflow for MV load profile analysis, but the headline 'spherical structure' is a mathematical artifact of z-scoring; the real value is the arc-shaped ordering and the generative model, both of which need stronger validation.","tokens_in":17097,"tokens_out":2721,"would_cite":false,"duration_ms":26299,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Standardized daily load profiles sit on a hidden sphere, and a principal curve through them reveals a continuous order, unifying outlier detection, clustering, and generative modeling.","keywords":["load profiles","spherical representation","principal component analysis","outlier detection","generative modelling","principal curves","von Mises-Fisher distribution","medium voltage grids"],"falsifier":"For any of the four municipalities, fit a general ellipsoid to the first three PCA coordinates and compare residuals with the best-fit sphere from the paper's optimization; if the radius $r$ after centering shows a systematic trend with polar or azimuth angle, or its distribution is not tightly concentrated near 1, the spherical assumption and the von Mises-Fisher generative model lose support.","tokens_in":15940,"feed_emoji":"⚡","tokens_out":5614,"duration_ms":50825,"temperature":0.7,"pith_summary":"This paper argues that the shapes of daily medium-voltage load profiles, once each profile is standardized by its own mean and standard deviation, live on a hypersphere. Projecting those standardized profiles with PCA into three dimensions gives a spherical cloud with an arc-shaped latent structure, and the paper shows the same pattern appearing in data from four Dutch municipalities. The authors use this single representation to replace four separate toolchains: cluster visualization, outlier detection, profile ordering, and generative modeling all happen in the same spherical coordinate system. If the claim holds, load-profile analysis becomes a problem of directional statistics rather than a patchwork of per-task methods, and continuous transitions between consumption types become directly observable and sampleable.","feed_headline":"Daily power load profiles trace an arc on a hidden sphere","feed_subtitle":"A single spherical model unifies outlier detection, ordering, and generative modeling for medium-voltage grids.","key_machinery":"The central object is the row-wise standardized load-profile matrix $X = \\frac{1}{\\sqrt{D}}\\hat{P}$, whose rows have unit norm; dot products between rows therefore equal the cosines of the angles between profiles, giving the Gram matrix $S = XX^{\\top}$. The argument runs through three linked pieces: PCA projects the hypersphere data into $\\mathbb{R}^3$, producing an arc-shaped spherical cloud; a principal curve $f(s)$ parameterized by $s \\in [0,1]$ is fitted through the cloud to recover the latent ordering; and a von Mises-Fisher distribution $f_{\\mathrm{vmf}}(z; \\mu_{\\mathrm{vmf}}, \\kappa)$ with $\\mu_{\\mathrm{vmf}} = f(s)$ generates new profiles that can be selected to lie between clusters. The paper's key simplifying assumption is that the three-dimensional projection is spherical rather than ellipsoidal.","core_discovery":"The load-bearing discovery is mathematical: after row-wise standardization, each profile vector has squared norm exactly $D$, so the profiles lie on a hypersphere of radius $\\sqrt{D}$; after normalization and PCA, the first three principal components yield a three-dimensional latent distribution that is approximately spherical and arc-shaped. The paper shows this structure is not an artifact of one dataset: the same eigenprofiles, radius concentration near 1, and single-mode angle distributions appear across four municipalities. On this basis it constructs a von Mises-Fisher generative model centered on a principal curve fitted to the arc, so sampling along the curve produces synthetic profiles that interpolate continuously between clusters. It also uses the fitted probability models in spherical coordinates to flag outlier profiles, including meter failures that per-meter time-series methods miss.","pith_inferences":["The hypersphere property follows only from row-wise standardization, so the same latent-sphere claim should hold for any shape-focused profile dataset, including low-voltage or hourly data; checking the radius distribution on such data would test this extension directly.","The arc ordering is reminiscent of seriation and pseudotime problems in other fields, so the principal-curve parameter $s$ could plausibly be correlated with spatial or demographic variables, an analysis the paper does not perform.","The three-component reconstruction deliberately discards high-frequency structure; a practical extension would model the residual as independent noise per time step while keeping the spherical generative model for shape, improving 15-minute accuracy.","The 95% confidence-interval outlier rule depends on the fitted von Mises and skew-normal distributions; a nonparametric density on the sphere would provide a robustness check for the outlier labels."],"forward_implications":["One spherical representation replaces separate pipelines: the same three coordinates show cluster structure, expose outliers, and serve as the sampling space for synthetic profiles.","Outlier detection in radius and angle coordinates catches systematic meter faults, such as meters storing absolute power values, that smoothing or per-meter historical methods cannot.","Each profile receives a scalar $s$ from the principal curve; clusters become intervals on $s$, and the curve gives a ranked, continuous ordering of areas by consumption similarity.","Sampling a von Mises-Fisher distribution centered at $f(s)$ yields profiles with continuous mixtures between clusters, such as a specified share of residential and commercial consumption, something discrete mixture models cannot do.","Because the spherical structure and eigenprofiles repeat across four municipalities, the directional-statistics toolkit transfers to other medium-voltage datasets."],"supporting_citations":[{"why":"Supplies the principal-curve method used to fit the one-dimensional curve through the spherical cloud and define the latent ordering.","marker":"[31]"},{"why":"Provides the robust spline principal-curve algorithm for data concentrated on a spherical shell, used to parameterize the arc.","marker":"[32]"},{"why":"Serves as the multivariate Gaussian baseline generative model that the von Mises-Fisher model is compared against.","marker":"[33]"},{"why":"Justifies using the dot-product similarity matrix and principal coordinate analysis, since it yields the same lower-dimensional projection as PCA on standardized data.","marker":"[30]"},{"why":"Documents the horseshoe or arc effect in principal coordinates analysis, which the paper cites to explain the arc-shaped latent structure.","marker":"[26]"},{"why":"Supplies the seriation background connecting arc-shaped projections to an underlying latent ordering of samples.","marker":"[28]"}],"fun_headline_variants":["Load profiles found on a hidden sphere, shaped as an arc","Spherical arc structure unifies load profile analytics","Hidden sphere and arc simplify grid load modeling","One sphere model covers load clustering, outliers, generation","Load profiles trace an arc across a latent sphere"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole construction assumes the three-dimensional PCA projection is spherical rather than ellipsoidal; if the true shape is an ellipsoid, the fitted von Mises-Fisher distributions and radius-based outlier thresholds are misspecified.","fun_headline_variants_meta":{"raw":{"variants":["Load profiles found on a hidden sphere, shaped as an arc","Spherical arc structure unifies load profile analytics","Hidden sphere and arc simplify grid load modeling","One sphere model covers load clustering, outliers, generation","Load profiles trace an arc across a latent sphere"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000667,"raw_usage":{"total_tokens":3022,"prompt_tokens":906,"completion_tokens":2116,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":522,"completion_tokens_details":{"reasoning_tokens":2042}},"tokens_in":522,"tokens_out":2116,"duration_ms":14290,"temperature":1.0,"reasoning_tokens":2042,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:16:54.584649+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For any of the four municipalities, fit a general ellipsoid to the first three PCA coordinates and compare residuals with the best-fit sphere from the paper's optimization; if the radius $r$ after centering shows a systematic trend with polar or azimuth angle, or its distribution is not tightly concentrated near 1, the spherical assumption and the von Mises-Fisher generative model lose support.","supporting_citations":[{"cited_title":"Principal curves,","cited_arxiv_id":null,"evidence_quote":"Supplies the principal-curve method used to fit the one-dimensional curve through the spherical cloud and define the latent ordering."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the robust spline principal-curve algorithm for data concentrated on a spherical shell, used to parameterize the arc."},{"cited_title":"Conditional Multivariate Elliptical Copulas to Model Residential Load Profiles From Smart Meter Data,","cited_arxiv_id":null,"evidence_quote":"Serves as the multivariate Gaussian baseline generative model that the von Mises-Fisher model is compared against."},{"cited_title":"Multidimensional Scaling,","cited_arxiv_id":null,"evidence_quote":"Justifies using the dot-product similarity matrix and principal coordinate analysis, since it yields the same lower-dimensional projection as PCA on standardized data."},{"cited_title":"Uncovering the Horseshoe Effect in Microbial Analyses,","cited_arxiv_id":null,"evidence_quote":"Documents the horseshoe or arc effect in principal coordinates analysis, which the paper cites to explain the arc-shaped latent structure."},{"cited_title":"A mathematical approach to seriation,","cited_arxiv_id":null,"evidence_quote":"Supplies the seriation background connecting arc-shaped projections to an underlying latent ordering of samples."}],"review_version":1}