{"id":"a840d318-4529-416c-9346-fa4e4ee398e5","arxiv_id":"2504.17512","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Sweep frequency response analysis, step excitation, and eigensystem realization give consistent dq admittance estimates for a simulated grid-forming inverter between 1 Hz and 100 Hz, with step-based methods being faster.","lead":"This paper compares three ways to measure the electrical 'admittance' of a grid-forming inverter from the outside, without access to its internal control code. Time-domain step tests are shown to match frequency-sweep results in the low-frequency range, which could let engineers identify stability models faster.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed cross-method consistency is partly enforced by construction: ERA's system order was iteratively chosen until its Bode plot matched SEM/SFRA, so the agreement in Fig. 7 is not an independent confirmation.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the ERA system order is selected to match the other frameworks, so the consistency between ERA and SEM/SFRA is partly enforced by construction. This stress-test pass agrees with that assessment and does not find a more serious flaw. The paper's qualitative conclusion, that step-based methods can provide low-frequency admittance estimates comparable to frequency scanning, is plausible and practically useful, but the evidence as presented is weakened by the circular order selection and by the shared step-response data between ERA and SEM. A secondary issue, the apparent dimensional inconsistency in Eq. (6), also needs attention but does not by itself overturn the central claim. The concrete test proposed here, using independent order-selection criteria or a known reference model, would settle whether the claimed consistency is genuine. Given the absence of code, data, or an external baseline, the conditional verdict is appropriate: the central idea merits publication only after the comparison methodology is made non-circular and the derivation is clarified. No ad hominem is intended; the concern is strictly about the argument's structure and evidence.","tokens_in":8503,"tokens_out":3082,"duration_ms":31441,"concrete_test":"Re-run the identification with pre-specified, independent order-selection rules: for ERA, choose the order from the singular-value spectrum of the Hankel matrix (e.g., retain singular values above a threshold or use AIC); for SEM and SFRA, use cross-validation or an information criterion on the respective time- and frequency-domain data, without reference to the other methods. Then overlay the three Bode plots over 1–100 Hz for all four admittance elements. If the curves still agree within, say, 1 dB in magnitude and 5 degrees in phase, the consistency claim is supported. A complementary test is to apply all three methods to a known linear time-invariant model (e.g., the small-signal linearization of the GFM) and compare each identified admittance against the true model; this removes the circularity entirely and also checks whether the d2c conversion in Eqs. (4)–(5) is accurate.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim, that all three identification frameworks produce consistent dq admittance estimates over 1–100 Hz, rests on the comparison shown in Fig. 7. In Section IV, the ERA model order is selected by an iterative procedure that explicitly uses the other two frameworks as the target: the text states that truncation is performed 'starting from the lowest order and increasing until the Bode plot of the admittance matrix approximates the results from the other two frameworks,' and the final sixth-order model 'best matches the other frameworks.' This makes the ERA result dependent on SEM and SFRA, so the observed overlap in Fig. 7 cannot be read as independent corroboration. Furthermore, ERA and SEM are both fed the same step-response data (Figs. 3 and 4), so their mutual agreement is expected and provides no additional evidence. The SEM and SFRA orders are chosen by fit-quality thresholds, not by matching each other, but the ERA order selection is the most direct circular step. A secondary concern is Eq. (6): if the i(s) entries are the identified transfer functions P(s) from Eq. (5), then multiplying by s/g is dimensionally incorrect (P(s) already has units of admittance); if i(s) denotes Laplace transforms of step responses, the notation conflicts with Eqs. (7)–(9), where i(s) are transfer functions. This needs clarification but is not the primary threat to the consistency claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares three black-box identification approaches for extracting the dq admittance of a grid-forming inverter in a simulated EMT testbed: SFRA (frequency-domain sinusoidal injection with tfest fitting), SEM (step excitation with tfest fitting), and ERA (step-response Hankel/SVD realization). The authors derive a discrete-to-continuous conversion for ERA, collect step responses from 1% voltage perturbations, fit transfer functions, and report Bode plots of the four admittance elements. The central claim is that all three frameworks give consistent admittance estimates in the 1-100 Hz range, with ERA and SEM providing faster identification than SFRA for low-frequency EMT stability studies.","tokens_in":8732,"tokens_out":4321,"duration_ms":38759,"significance":"If the consistency claim holds, the practical payoff is substantial: two step perturbations could replace dozens of frequency scans for low-frequency admittance characterization, which is valuable when inverter firmware is inaccessible. The paper is clearly organized and provides a useful side-by-side comparison of the three techniques on a realistic GFM testbed, with transparent fit metrics (NRMSE) for SEM and SFRA. However, the evidence for consistency is weakened by the ERA system order being selected so that its Bode plot matches the other frameworks, and by a dimensional inconsistency in Eq. (6). These issues are fixable within the manuscript's scope, so the result is not fatally compromised but the central claim needs stronger support.","major_comments":[{"comment":"The ERA model order is selected iteratively 'starting from the lowest order and increasing until the Bode plot of the admittance matrix approximates the results from the other two frameworks,' and the sixth-order model is chosen because it 'best matches the other frameworks.' This makes the ERA-SEM/SFRA agreement in Fig. 7 partly enforced by construction rather than independently discovered. Furthermore, ERA and SEM use the same step-response data, so their mutual overlap in Fig. 7 is expected and provides no additional evidence of consistency. Please choose the ERA order using an independent criterion, such as the singular-value decay of the Hankel matrix or validation on a held-out portion of the step response, and then re-evaluate the comparison.","section":"Section IV, ERA truncation"},{"comment":"Equation (6) contains a dimensional inconsistency. If the entries i(1)od(s), i(1)oq(s), i(2)od(s), i(2)oq(s) are the identified transfer functions P(s) from Eq. (5), they already have units of admittance, so the extra factor s/g makes the right-hand side dimensionally incorrect. If these entries instead denote Laplace transforms of step responses, then the notation conflicts with Eqs. (7)-(9), where i(s) are transfer functions. Please clarify the meaning of i(s) and correct Eq. (6) accordingly.","section":"Equation (6)"},{"comment":"The discrete-to-continuous conversion d2c(WERA(z)) = gP(s)/s is exact only under specific sampling and hold assumptions, but the paper does not state which d2c method is used or validate its accuracy over 1-100 Hz. Because the ERA admittance estimate depends directly on this conversion, please report the d2c method and validate it, for example by comparing the identified continuous-time model against an independent SFRA scan or by checking its prediction on a separate time-domain experiment.","section":"Equations (4)-(5)"}],"minor_comments":[{"comment":"The phrase 'Laplace transform as shown in (6)' should refer to Eq. (5), where P(s) is recovered from WERA(z); Eq. (6) is the admittance matrix construction, not a Laplace transform definition.","section":"Section III.A"},{"comment":"The notation Z^{-1}[.] applied to a continuous-time function is ambiguous; please define the sampling period and state explicitly that the equality holds at sampling instants under an assumed zero-order hold.","section":"Equations (1)-(3)"},{"comment":"Please add a shaded vertical band or markers indicating the 1-100 Hz range to which the abstract's consistency claim refers, so that readers can verify the claim directly from the figure.","section":"Figure 7"},{"comment":"The NRMSE fit percentages are reported without defining the normalization; please state the formula or cite the MATLAB documentation for tfest's fit metric.","section":"Figures 5 and 6"},{"comment":"The text notes that SFRA 'shows a different response' above 100 Hz; please explicitly reconcile this observation with the abstract's 1-100 Hz consistency claim in the conclusions.","section":"Section IV, SFRA discussion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the scope of eess.SY, but the central consistency claim is currently supported by a partly circular ERA order-selection procedure. I recommend major revision because the issue is fixable by re-running the ERA identification with an independent order-selection rule, correcting Eq. (6), and validating the d2c conversion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful thing here is the side-by-side: ERA, SEM, and SFRA applied to the same simulated grid-forming inverter, with a clear testbed and a careful re-derivation of ERA for step inputs. That exact comparison is not in the prior literature, and the practical message—two step injections can give you low-frequency admittance comparable to frequency scanning—is worth having on record.\n\nThe paper does some things well. The ERA derivation that converts a step response into an impulse-response realization via an extended discrete-time system is clearly laid out. The discussion of SNR and time/computation trade-offs is sensible. The simulation results in the 1–100 Hz band are plausible, and the SFRA vs SEM/ERA divergence above 100 Hz is exactly what you would expect given the excitation.\n\nThe main soft spot is the circularity in the ERA order selection. The text says the order was increased until the Bode plot matched the other two frameworks. So the agreement in Fig. 7 is partly enforced by construction, not independently discovered. Also, ERA and SEM share the same step-response data, so their overlap is expected. This does not collapse the paper, but it changes the claim from 'three independent methods agree' to 'two methods that use the same data and one of them is tuned to match can look alike.' A revision should either use an independent order-selection criterion (e.g., Hankel singular values, prediction error) or explicitly frame the result as a demonstration that such matching is possible.\n\nThe secondary issue is Eq. (6). If the i(s) entries are already the identified transfer functions from Eq. (5), the extra s/g factor is dimensionally wrong. If they are meant to be Laplace transforms of raw step responses, then the notation conflicts with Eqs. (7)–(9), where i(s) are transfer functions. This needs clarification, but it is minor and fixable.\n\nAlso minor: no code or data, but the parameters in Table I are detailed enough that the simulation is reproducible in principle. The citation pattern is fine.\n\nBottom line: this is a modest, practical contribution for engineers who need to choose an identification method for black-box admittance extraction. It is not a breakthrough, but it is honest work with a real flaw in the comparison methodology. It deserves a serious referee and a conditional acceptance after revision, not a desk reject.","headline":"A useful comparative benchmark of dq admittance identification methods whose central consistency claim is partly enforced by construction; deserves a revised round, not a desk reject.","tokens_in":9347,"tokens_out":1351,"would_cite":false,"duration_ms":13900,"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":"Three identification methods agree on grid-forming inverter admittance from 1 to 100 Hz.","keywords":["admittance identification","dq-frame admittance","grid-forming inverter","sweep frequency response analysis","step excitation method","eigensystem realization algorithm","impedance-based stability","small-signal stability"],"falsifier":"Run the three identification methods on a different grid-forming inverter (or a hardware testbed with known admittance), choosing ERA's order from a pre-specified criterion such as the drop-off of Hankel singular values rather than by matching SEM/SFRA; if ERA's Bode curves then depart from SFRA inside 1-100 Hz, the paper's consistency claim fails. A simpler check is to compare step-derived transfer functions against a dense sine sweep in a case where the true admittance is known from the controller parameters.","tokens_in":8256,"feed_emoji":"⚡","tokens_out":9096,"duration_ms":85103,"temperature":0.7,"pith_summary":"The paper sets out to show that three black-box identification techniques -- sweep frequency response analysis, step excitation, and eigensystem realization -- yield consistent estimates of the dq admittance of a grid-forming inverter over 1 Hz to 100 Hz. The stakes are practical: inverter firmware is proprietary, so stability studies often have to build admittance models from external measurements, and a frequency sweep demands many sinusoidal injections, while the two step-based methods need only two perturbations. The paper demonstrates the consistency on a simulated grid-forming inverter by comparing the Bode plots of the four admittance entries, finding that ERA and SEM overlap closely with each other and with SFRA in the low-frequency band. Above 100 Hz the curves diverge, which the paper attributes to the higher signal-to-noise ratio of frequency scanning at high frequencies.","feed_headline":"Two step tests match a full frequency sweep for inverter admittance","feed_subtitle":"Grid-forming inverter dq admittance matches across three identification methods from 1 to 100 Hz.","key_machinery":"The operational core is the dq-frame admittance matrix $Y(s)$, defined with a negative sign so that positive current flows from the inverter to the grid, and its entries are extracted from the transfer functions between dq voltage perturbations and dq current responses. For ERA, the key identity is that feeding a step response to the algorithm is equivalent to identifying a discrete-time extended system $W_{ERA}(z)=W_{int}(z)W_{DD}(z)$, where $W_{int}=z/(z-1)$ is a discrete integrator; the continuous plant is then recovered through $P(s)=d2c(W_{ERA}(z))s/g$. This step-to-impulse equivalence lets ERA work with a simple voltage step, from which it builds a Hankel matrix, applies singular-value decomposition, and truncates at a chosen system order (six in this study). SEM fits continuous-time transfer functions directly from the same step-response data with a low-order model (four poles), and SFRA fits a transfer function to the measured frequency-response points; all three feed their estimates into the same $2\\times2$ admittance matrix structure.","core_discovery":"The central claim is that the dq admittance matrix of a grid-forming inverter can be identified equivalently in the time domain and the frequency domain, with agreement across methods from 1 Hz to 100 Hz. For each element of the admittance matrix, the inverter is perturbed separately on the d-axis and q-axis voltages: ERA and SEM apply a 1% step change to $v_{gd}$ and then to $v_{gq}$, record the resulting dq currents, and convert those time traces into transfer functions; SFRA instead injects 0.1 V peak-to-peak sinusoids at 100 discrete frequencies, computes phasors with a Fourier transform, and fits a transfer function through the points. The four fitted entries -- $Y_{dd}$, $Y_{dq}$, $Y_{qd}$, and $Y_{qq}$ -- overlap in Bode magnitude and phase within the 1-100 Hz band. The paper treats the divergence above 100 Hz as a regime difference: time-domain step methods lose signal-to-noise ratio at high frequencies, while SFRA's discrete-point excitation maintains accuracy there.","pith_inferences":["An extension the paper does not make: if the 1-100 Hz agreement survives on other inverter topologies and on hardware-in-the-loop testbeds, step-based identification could become the default screening tool for EMT studies, with frequency sweeps reserved for verifying the high-frequency tail.","The ERA order (six) is selected iteratively until its Bode plot matches the other two methods; a stronger test of the paper's consistency claim would fix the order in advance by an independent rule, such as the decay of Hankel singular values, and then compare the resulting curves.","The discrete-to-continuous conversion used for ERA is itself a candidate source of the high-frequency mismatch; sampling faster or using a different conversion could reveal whether the divergence above 100 Hz is a physical property of step excitation or an artifact of the conversion."],"forward_implications":["A grid-forming inverter's dq admittance in the 1-100 Hz band can be obtained from two step perturbations instead of a 100-point frequency sweep.","Time-domain identification is therefore a faster route to admittance models for low-frequency stability studies, with lower computational demand than frequency scanning.","Above 100 Hz, frequency scanning remains the more trustworthy of the three approaches, because it excites each discrete frequency with a strong, narrow-band signal.","ERA and SEM produce overlapping results from the same step data, so the choice between them can rest on convenience and noise handling rather than on the physics captured."],"supporting_citations":[{"why":"Defines the dq-frame admittance matrix and the stacking procedure used to assemble the ERA admittance estimate from the four transfer-function measurements.","marker":"[6]"},{"why":"Supplies the eigensystem realization algorithm itself, including the Hankel-matrix and singular-value-decomposition procedure that ERA applies to the step-response data.","marker":"[12]"},{"why":"Provides the step-excitation method for extracting dq admittance, including the transfer-function fitting approach that SEM applies to the time-domain measurements.","marker":"[14]"}],"fun_headline_variants":["Three methods agree on inverter admittance from 1 to 100 Hz","Step and sweep tests yield same admittance for grid inverters","Time-domain steps match frequency sweep for dq admittance","Inverter admittance identical via three identification techniques","Methods agree: grid-former admittance from 1 to 100 Hz"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The agreement between ERA and the other two methods rests on the selected system order (six), which the paper chooses iteratively until the Bode plots match, so the consistency is partly built into the comparison rather than independently discovered.","fun_headline_variants_meta":{"raw":{"variants":["Three methods agree on inverter admittance from 1 to 100 Hz","Step and sweep tests yield same admittance for grid inverters","Time-domain steps match frequency sweep for dq admittance","Inverter admittance identical via three identification techniques","Methods agree: grid-former admittance from 1 to 100 Hz"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000187,"raw_usage":{"total_tokens":1312,"prompt_tokens":912,"completion_tokens":400,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":313}},"tokens_in":528,"tokens_out":400,"duration_ms":3922,"temperature":1.0,"reasoning_tokens":313,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:38:15.568917+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the three identification methods on a different grid-forming inverter (or a hardware testbed with known admittance), choosing ERA's order from a pre-specified criterion such as the drop-off of Hankel singular values rather than by matching SEM/SFRA; if ERA's Bode curves then depart from SFRA inside 1-100 Hz, the paper's consistency claim fails. A simpler check is to compare step-derived transfer functions against a dense sine sweep in a case where the true admittance is known from the controller parameters.","supporting_citations":[{"cited_title":"Time-domain measurement-based dq-frame ad- mittance model identification for inverter-based resources,","cited_arxiv_id":null,"evidence_quote":"Defines the dq-frame admittance matrix and the stacking procedure used to assemble the ERA admittance estimate from the four transfer-function measurements."},{"cited_title":"An eigensystem realization algorithm for modal parameter identification and model reduction,","cited_arxiv_id":null,"evidence_quote":"Supplies the eigensystem realization algorithm itself, including the Hankel-matrix and singular-value-decomposition procedure that ERA applies to the step-response data."},{"cited_title":"Stability analysis of real-world subsyn- chronous oscillations via black-box emt models,","cited_arxiv_id":null,"evidence_quote":"Provides the step-excitation method for extracting dq admittance, including the transfer-function fitting approach that SEM applies to the time-domain measurements."}],"review_version":1}