{"id":"1400dda3-d375-44d7-9904-736da5bb4368","arxiv_id":"2506.06177","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Mean force emission theory with the attractive Kelbg potential reproduces classical molecular dynamics bremsstrahlung spectra in electron-ion plasmas and links emission to absorption and conductivity.","lead":"This paper extends a mean force emission theory to electron-ion plasmas, where the attractive electron-ion interaction and short-range quantum effects are modeled with a smoothed interaction potential. It shows the theory reproduces computer-simulated bremsstrahlung spectra and connects emission to absorption and electrical conductivity, with implications for warm dense matter.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Validation is partially circular: the potential of mean force is extracted from the same MD trajectories that serve as the reference spectrum, so the reported agreement does not yet establish independent predictive power for electron-ion plasmas.","rationale":"The reader identified the same weakest assumption: the high-frequency spectrum is computed from a potential of mean force derived from the same MD trajectories used for validation. This is indeed the most load-bearing concern because the central claim is an empirical validation of mean force emission theory for electron-ion plasmas. The paper's comparisons are internally consistent and physically instructive, and the framework is a reasonable extension of prior work, but the circularity means the agreement does not yet demonstrate predictive power for conditions or methods not already contained in the simulation. The reader's conditional verdict correctly reflects this: the claim is plausible and well supported as a consistency check, but an independent test of the potential of mean force (e.g., from HNC closure or from a separate simulation with different particle number) would be needed to remove the circularity. I found no additional concern that would flip the verdict to accept or reject: the bound-state peaks are explicitly acknowledged as outside the two-body model, and the plasma-frequency peak at strong coupling is likewise transparently discussed as a limitation. The lack of error bars and absence of public data are reproducibility issues, but they are secondary to the input-dependence of the validation. Therefore the conditional verdict stands without change.","tokens_in":23113,"tokens_out":3834,"duration_ms":36669,"concrete_test":"Recompute the high-frequency emission spectrum j_h(omega) (Sec. IIIA) using a potential of mean force obtained independently of the reference MD run. For a first check, use the weakly coupled analytic form w_ei(r) = U_K(r) exp(-r/lambda_D) at Gamma=0.01, Theta=468, and compare the resulting j_h(omega) with the MD spectrum over omega/omega_pe > 1. If the spectrum changes by more than the visual scatter of the MD data, the agreement in Fig. 8 is partially attributable to using MD-derived input. A stronger test: compute w_ei from an Ornstein-Zernike/HNC closure with the same Kelbg potential and repeat the Gamma=1, Theta=0.468 case. If the high-frequency or low-frequency predictions shift substantially, an independent-source demonstration is required before the framework can be claimed validated for real plasmas.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that mean force emission theory, using the potential of mean force w_ei(r) = -k_B T ln g_ei(r) (Eq. 22), reproduces the MD bremsstrahlung spectrum. But in every comparison shown (Figs. 8-10), w_ei is obtained from the same MD simulation that provides the reference spectrum (Sec. IIA, Eq. 9), with a smoothing spline applied to the noisy g_ei(r). This makes the high-frequency model (Sec. IIIA, Eq. 24) a consistency test rather than an independent prediction: the theory is given the exact equilibrium pair structure of the very system it is then asked to describe. The agreement could partly reflect information already encoded in g_ei(r) — particularly the short-range force plateau that controls the high-frequency decay — rather than validating the two-body-in-mean-force reduction. The paper acknowledges this by listing Ornstein-Zernike theory as a possible alternative source of w_ei, but it never performs that independent test. The low-frequency Drude model (Sec. IIIB) uses a collision frequency derived from the same MD-based PMF, so the same caveat applies there. Because the conclusion 'the same framework can be applied to describe electron-ion plasmas' rests on this agreement, the lack of an independent PMF is the most load-bearing weakness. A secondary symptom is the sensitivity of the high-frequency spectrum to the smoothing of g_ei at small r, where the MD data are explicitly noisy; this could inject uncontrolled artifacts into the force that drives emission.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends mean force emission theory to attractively interacting electron-ion plasmas. Classical molecular dynamics simulations use a Kelbg pseudopotential for electron-ion interactions to avoid Coulomb collapse, and the bremsstrahlung emission coefficient is computed from the force autocorrelation function. The theory is built from two limits: high-frequency emission from isolated two-body collisions in the potential of mean force w_ei = -k_B T ln g_ei (Sec. IIIA, Eq. (24)) and low-frequency Drude-like emission with a mean-force kinetic-theory collision frequency (Sec. IIIB, Eq. (26)), combined into a comprehensive model by Eq. (27). The results are compared with MD spectra across Gamma = 0.01-1 and Theta = 0.468-468, including a discussion of bound-state peaks, classical and Born Gaunt factors, the plasma-frequency peak at strong coupling, and the relation between emission, absorption, and the real part of the dynamic conductivity via Eq. (14).","tokens_in":23344,"tokens_out":5609,"duration_ms":58734,"significance":"If validated, the framework provides a practical classical route to bremsstrahlung emission, absorption, and dynamic conductivity in strongly coupled plasmas, with the linear-response bridge in Eq. (14) making the conclusions transferable to transport coefficients. The paper's strengths include the explicit check of the autocorrelation formalism against the direct dipole formula (Fig. 7), the analytic screening correction derived in Appendix B, and a transparent qualitative account of Kelbg-scale and bound-state effects. The quantitative support is strongest for weakly coupled cases and for the high- and low-frequency limits. The current validation, however, is partly circular because the potential of mean force is taken from the same MD trajectories that provide the reference spectra; an independent source for w_ei would substantially strengthen the central claim.","major_comments":[{"comment":"The main validation is partially circular. In every comparison shown in Figs. 8-10, the potential of mean force used in Eqs. (22)-(24) and in the collision frequency of Sec. IIIB is extracted from the same MD runs that generate the reference spectrum, via Eq. (9) followed by a smoothing spline. Agreement between theory and MD is therefore a consistency test of the mean-force reduction rather than an independent prediction for electron-ion plasmas. The manuscript notes at Eq. (22) that g_ei can be obtained from Ornstein-Zernike theory or hypernetted-chain closure, but it does not perform that independent computation. Please either recompute w_ei from an OZ/HNC closure (or another independent route) at the simulated conditions and repeat the spectral comparison, or provide a sensitivity study showing that the spectra are insensitive to the spline smoothing and to the small-r noise in g_ei. Without one of these, the conclusion that 'the same framework can be applied to describe electron-ion plasmas' is not fully supported.","section":"Sec. III, Eq. (22)"},{"comment":"At strong coupling the comprehensive model does not reproduce the plasma-frequency peak, and the paper's argument that this peak is due to free-electron collective motion relies mainly on the comparison with the repulsive-Kelbg spectrum. The attractive MD case at Gamma = 1, Theta = 0.468 also contains bound-state contributions (Sec. IIB3), and no decomposition of the MD spectrum into free and bound parts is provided. The claim that the Drude correction 'cannot capture' this peak would be strengthened by isolating the free-electron contribution, for example by recomputing the spectrum from unbound trajectories only, or by a more quantitative comparison with a model that includes the oscillatory force autocorrelation at intermediate timescales.","section":"Sec. IIIC, Eq. (27), Fig. 10"},{"comment":"The piecewise model in Eq. (43) is used to compute the frequency-averaged Gaunt factor and the screening-induced reduction of total bremsstrahlung power, but the connection frequency x is never defined or determined. The numerical value of the screening correction in Fig. 15 depends on x, so the quantitative claim that screening 'slightly reduces the total bremsstrahlung power' is not fully specified. Please state how x is chosen (for example, by a continuity condition, a matching frequency, or a fit) and show the sensitivity of G-bar to x over a reasonable range.","section":"Sec. IV, Eq. (43), Fig. 15"}],"minor_comments":[{"comment":"The thermostat name should be 'Nosé-Hoover' rather than 'Nose-Hoover', and the text 'data.39.' appears to contain a typographical artifact from reference formatting.","section":"Sec. IIA"},{"comment":"The caption's scaling instructions ('multiplied by 10, 0.02, 10^-9') are easy to misread because the Gamma, Theta pairs are not aligned with the listed factors; a table or explicit per-curve labels would improve clarity.","section":"Fig. 2 caption"},{"comment":"The orbit-frequency argument in Fig. 5 is qualitative: the shaded 'important frequency' regions are based on a decrease in force and a decrease in bound-electron probability, but no occupancy weighting is given. Please state the weighting used to define the shaded regions or label them explicitly as heuristic ranges.","section":"Sec. IIB3, Eqs. (11)-(12)"},{"comment":"The reported high-frequency decay power law of approximately (hbar omega/k_B T)^-4.1 is given without an error estimate or a stated fitting range; please specify how this exponent was extracted from the MD data.","section":"Sec. IVD, Fig. 13"},{"comment":"The Barry approximation introduces h, h_infty, q, and G in an order that is slightly confusing; defining all symbols before displaying the approximation would improve readability.","section":"Appendix B, Eq. (B7)"}],"recommendation":"major_revision","confidential_remarks":"The central issue is the same-MD circularity for the potential of mean force. I do not recommend rejection because the paper is transparent about the limitation and the theoretical framework is valuable, but the requested independent PMF test is necessary to move the validation from consistency check to prediction. The undefined connection frequency x in Eq. (43) should also be fixed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, honest extension of the group's mean force emission theory to attractive electron-ion interactions, and it deserves a serious referee. The novel pieces are real: replacing Coulomb collapse with the Kelbg potential, identifying the bound-state peaks and their orbit-frequency origin, and showing the high-frequency decay and the strong-coupling plasma-frequency peak. The relation between emission, absorption, and real sigma via linear response is clean and useful. The authors also acknowledge the Kelbg potential's limits at very high frequency and give a plausible screening correction in the appendix. Credit where due: the MD comparisons span a useful range of Gamma and Theta, and the two-frequency-limit construction is physically transparent.\n\nThe soft spot the stress test flags is real and load-bearing. In every comparison, w_ei = -kT ln g_ei is extracted from the same MD trajectories that produce the reference spectrum. So the high-frequency model is fed the exact equilibrium pair structure of the system it is asked to reproduce. That makes the agreement a consistency check, not an independent validation of the mean-force reduction. The authors list Ornstein-Zernike and HNC as alternative sources for w_ei but never run that test. I would want that before believing the framework has predictive power for electron-ion plasmas. The low-frequency Drude model inherits the same issue because nu_ei is computed from the same potential.\n\nSecondary issues in proportion: the comprehensive model misses the bound-state peaks and the strong-coupling plasma-frequency peak, which the authors state; the smoothing spline on noisy g_ei at small r could inject artifacts; and there are no error bars or public code and data. None of these are fatal, but they all push in the same direction: the paper is a convincing demonstration of internal consistency and a promising step, rather than a fully independent validation.\n\nWho gains: people working on radiation transport in warm dense matter and on dynamic conductivity in strongly coupled plasmas. A serious referee should engage. The circularity can be addressed with an independent PMF test or by reframing the results as consistency checks. For the record, I think the authors are not hiding anything; they openly state where w_ei comes from. The limitation is in the design, not in their presentation. I would accept this for review.","headline":"Solid, honest extension of mean force emission theory to electron-ion plasmas, but the validation is partly circular because the potential of mean force comes from the same MD simulations used as the reference.","tokens_in":23983,"tokens_out":1748,"would_cite":true,"duration_ms":19002,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.25.Os","52.65.Yy"],"model":"deepseek-v4-flash","headline":"Classical bremsstrahlung in electron-ion plasmas, from weak to strong coupling, is computed from the potential of mean force plus a Drude collision-frequency model, and the same formula yields absorption and dynamic conductivity.","keywords":["mean force emission theory","bremsstrahlung","electron-ion plasma","strongly coupled plasma","Kelbg potential","potential of mean force","dynamic conductivity","Gaunt factor"],"falsifier":"A decisive check is to feed the theory a potential of mean force obtained independently of the spectrum it is meant to predict — for example from an Ornstein-Zernike solution with a hypernetted-chain closure, or from a differently seeded simulation — and see whether the high- and low-frequency emission predictions still match the MD spectra. A second, frequency-specific check: compute the spectrum for $\\hbar\\omega/k_B T$ well above the Kelbg decay frequency $\\omega_{\\max,K}$ with an exact quantum-mechanical Gaunt factor, since the paper reports a $\\propto(\\hbar\\omega/k_B T)^{-4.1}$ decay there that disagrees with the Born approximation's $\\propto(\\hbar\\omega/k_B T)^{-0.5}$.","tokens_in":22822,"feed_emoji":"⚡","tokens_out":19455,"duration_ms":158566,"temperature":0.7,"pith_summary":"Mean force emission theory, previously built for repulsively interacting positron-ion plasmas, is here extended to attractive electron-ion plasmas and claimed to reproduce the classical bremsstrahlung spectrum across coupling strengths. The radiation is computed from binary electron-ion collisions moving in the equilibrium potential of mean force $w_{ei}(r) = -k_B T \\ln g_{ei}(r)$ at high frequency, and from a Drude-like collision-frequency model at low frequency. Molecular dynamics simulations that avoid Coulomb collapse through the Kelbg potential confirm the predictions, except for peaks from classically bound electron-ion states and a strong-coupling peak near the plasma frequency that the standard Drude correction misses. Because the emission spectrum is the Fourier transform of a force autocorrelation function, the paper ties it to bremsstrahlung absorption and the real part of the dynamic conductivity through one linear-response relation. If the claim holds, the same framework gives emission, absorption, and conductivity in warm dense matter from a static pair distribution alone.","feed_headline":"Plasma bremsstrahlung reduced to one equilibrium force","feed_subtitle":"The same averaged electron-ion force reproduces simulated spectra and ties emission, absorption, and conductivity.","key_machinery":"The load-bearing object is the potential of mean force $w_{ei}(r) = -k_B T \\ln g_{ei}(r)$, the canonical average of the electron-ion force at fixed separation, obtained from the electron-ion radial distribution function $g_{ei}(r)$; the paper computes $g_{ei}(r)$ from the MD trajectories, though the definition itself admits experimental or integral-equation sources. At high frequency the machinery is a sum over isolated binary collisions: Newton's equation is integrated in this potential, and the single-collision radiation spectrum $W(\\omega, v, b)$ is averaged over impact parameters and a Maxwellian speed distribution. At low frequency it is a Drude-like spectrum built from an exponentially decaying velocity autocorrelation $Z(t) = (k_B T/m_e)\\exp(-\\nu_{ei} t)$, with the electron-ion collision frequency $\\nu_{ei}$ taken from mean force kinetic theory; the combined solution multiplies the high-frequency spectrum onto the low-frequency Drude form. The Kelbg potential $U_K(r)$, which is linear at short range and saturates the force inside the thermal de Broglie wavelength, is what lets the attractive MD simulations run without Coulomb collapse and sets the high-frequency cutoff. The identity that carries the argument to other observables is $j(\\omega) = \\omega^2 k_B T/(2\\pi^2 c^3 \\epsilon_0)\\,\\mathrm{Re}\\,\\sigma(\\omega)$, obtained by writing emission as the Fourier transform of the force autocorrelation function.","core_discovery":"The paper's central claim is that the same mean force emission theory validated on repulsive systems describes electron-ion plasmas: at frequencies above the plasma frequency the emission coefficient is obtained by averaging the single-collision radiation spectrum over impact parameters and speeds in the potential of mean force, and at low frequencies it is reproduced by a Drude form whose collision frequency comes from mean force kinetic theory. Compared with molecular dynamics, the combined model matches the simulated spectra across the coupling strengths $\\Gamma = 0.01$ to $1$ and the degeneracies studied, with exceptions the paper identifies: classically bound electron-ion states produce peaks that the two-body theory cannot capture (their locations are predicted by a circular-orbit model), and a peak near the plasma frequency at strong coupling is not produced by the standard Drude low-frequency form, though it also appears in repulsive, bound-state-free simulations and is attributed to correlated free-electron motion. The paper further argues that the emission coefficient, governed by the force autocorrelation function $\\langle \\dot{\\mathbf{J}}(t)\\cdot\\dot{\\mathbf{J}}(0)\\rangle$, is related to absorption and to $\\mathrm{Re}\\,\\sigma(\\omega)$ by $j(\\omega) = \\omega^2 k_B T/(2\\pi^2 c^3 \\epsilon_0)\\,\\mathrm{Re}\\,\\sigma(\\omega)$, so the framework extends to those transport coefficients. The Kelbg potential is used deliberately as a numerical proxy for quantum mechanics; the paper states it should not be read as an accurate description of physical warm dense matter, particularly at very high frequencies where the predicted decay disagrees with the Born approximation.","pith_inferences":["The framework suggests a practical route to opacities and electrical conductivities of warm dense matter from static structural data alone, if $g_{ei}(r)$ can be supplied by experiment or by an integral-equation theory; the paper itself only tests the case where the pair distribution comes from the same MD run as the spectrum.","A cleaner test of the two-body assumption would be to remove classically bound orbits from the MD spectra and check whether the residual free-electron spectrum matches mean force emission theory exactly over the full frequency range, not just where the bound-state peaks are absent.","The link between the emission spectrum and $\\mathrm{Re}\\,\\sigma(\\omega)$ implies an experimentally accessible signature: optical conductivity measurements on a strongly coupled plasma should show the same near-plasma-frequency peak structure as the radiation spectrum, since the paper argues both trace to the same correlated electron motion.","Replacing the Kelbg potential with a different short-range pseudopotential should shift the high-frequency cutoff but leave the low-frequency and plasma-frequency features unchanged; identifying which spectral features are insensitive to that choice would mark which predictions are stable signatures of the mean force concept."],"forward_implications":["The same pair distribution function that defines the potential of mean force supplies both the high-frequency binary-collision spectrum and the low-frequency electron-ion collision frequency, so no radiation-specific parameter is needed.","Screening enters through the potential of mean force and produces a plateau near the plasma frequency that lowers the frequency-averaged Gaunt factor; screening therefore slightly reduces total bremsstrahlung power even in weakly coupled, non-degenerate plasmas.","Because emission, absorption, and the real part of the dynamic conductivity share one linear-response formula, the conclusions of the study transfer to both transport coefficients, not just radiation.","The Kelbg potential reproduces the quantum cutoffs inserted by hand into classical Gaunt factors at low frequency, but deviates from the Born approximation at very high frequencies, so the quantum decay of the spectrum is only qualitatively captured.","The strong-coupling peak near the plasma frequency also appears in repulsive, bound-state-free simulations, indicating it comes from correlated free-electron motion; the standard Drude correction cannot produce it."],"supporting_citations":[{"why":"The prior mean force emission theory for repulsive positron-ion plasmas that this paper extends to attractive electron-ion plasmas.","marker":"[12]"},{"why":"Improved Kelbg potential for correlated Coulomb systems; with [28], the basis of the pseudopotential used in the simulations.","marker":"[14]"},{"why":"Supplies the specific temperature-dependent Kelbg potential form and fit parameter used as Eq. (1) and Eqs. (2)-(3).","marker":"[28]"},{"why":"Standard reference for the emission-coefficient formula (Eq. 5) and for the classical Gaunt-factor baselines the paper compares against.","marker":"[18]"},{"why":"Dawson-Oberman linear-dielectric calculation giving the screened low-frequency Gaunt-factor limit that mean force emission theory reproduces and generalizes.","marker":"[19]"},{"why":"Oster's radio-frequency emission, absorption, and conductivity model; the classical Gaunt-factor baseline that mean force emission theory is compared against and generalizes.","marker":"[20]"},{"why":"Earlier molecular dynamics study of optical conductivity with Kelbg potentials; the paper confirms its conclusion that very high frequency emission is not accurately captured.","marker":"[31]"},{"why":"Source of the potential-of-mean-force definition $w_{ei} = -k_B T \\ln g_{ei}$ (Eq. 22) and of the linear response background for Eq. (14).","marker":"[44]"},{"why":"Kubo's linear response relation for conductivity, the quantum formula whose classical limit ties emission to $\\mathrm{Re}\\,\\sigma(\\omega)$.","marker":"[45]"},{"why":"Mean force kinetic theory used to compute the electron-ion collision frequency entering the low-frequency Drude model.","marker":"[49]"}],"fun_headline_variants":["One force governs plasma bremsstrahlung, absorption, conductivity","Kelbg potential adds quantum effects to classical force model","Strong coupling peak eludes standard Drude correction","Bound states add peaks to plasma emission spectrum","Mean force model matches molecular dynamics with caveats"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that high-frequency emission can be treated as independent two-body electron-ion collisions moving in the equilibrium potential of mean force, with the surrounding plasma entering only through that equilibrium average; the paper tests this using a potential of mean force extracted from the very same molecular dynamics trajectories that produce the spectrum, not from an independent source, so collective dynamics beyond the average are not separately tested.","fun_headline_variants_meta":{"raw":{"variants":["One force governs plasma bremsstrahlung, absorption, conductivity","Kelbg potential adds quantum effects to classical force model","Strong coupling peak eludes standard Drude correction","Bound states add peaks to plasma emission spectrum","Mean force model matches molecular dynamics with caveats"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000551,"raw_usage":{"total_tokens":2720,"prompt_tokens":1125,"completion_tokens":1595,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":741,"completion_tokens_details":{"reasoning_tokens":1520}},"tokens_in":741,"tokens_out":1595,"duration_ms":12539,"temperature":1.0,"reasoning_tokens":1520,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:59:22.583332+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check is to feed the theory a potential of mean force obtained independently of the spectrum it is meant to predict — for example from an Ornstein-Zernike solution with a hypernetted-chain closure, or from a differently seeded simulation — and see whether the high- and low-frequency emission predictions still match the MD spectra. A second, frequency-specific check: compute the spectrum for $\\hbar\\omega/k_B T$ well above the Kelbg decay frequency $\\omega_{\\max,K}$ with an exact quantum-mechanical Gaunt factor, since the paper reports a $\\propto(\\hbar\\omega/k_B T)^{-4.1}$ decay there that disagrees with the Born approximation's $\\propto(\\hbar\\omega/k_B T)^{-0.5}$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The prior mean force emission theory for repulsive positron-ion plasmas that this paper extends to attractive electron-ion plasmas."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Improved Kelbg potential for correlated Coulomb systems; with [28], the basis of the pseudopotential used in the simulations."},{"cited_title":"Bekefi ,\\ title title Radiation processes in plasmas , \\ @noop journal journal Wiley series in plasma physics \\ ( year 1966 ) NoStop","cited_arxiv_id":null,"evidence_quote":"Standard reference for the emission-coefficient formula (Eq. 5) and for the classical Gaunt-factor baselines the paper compares against."},{"cited_title":"Dawson \\ and\\ author C","cited_arxiv_id":null,"evidence_quote":"Dawson-Oberman linear-dielectric calculation giving the screened low-frequency Gaunt-factor limit that mean force emission theory reproduces and generalizes."},{"cited_title":"Morozov , author H","cited_arxiv_id":null,"evidence_quote":"Earlier molecular dynamics study of optical conductivity with Kelbg potentials; the paper confirms its conclusion that very high frequency emission is not accurately captured."},{"cited_title":"\\ Hansen \\ and\\ author I","cited_arxiv_id":null,"evidence_quote":"Source of the potential-of-mean-force definition $w_{ei} = -k_B T \\ln g_{ei}$ (Eq. 22) and of the linear response background for Eq. (14)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Mean force kinetic theory used to compute the electron-ion collision frequency entering the low-frequency Drude model."}],"review_version":1}