{"id":"9f3f4f21-e8a1-48ff-afd6-310137f89800","arxiv_id":"2608.12264","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In a closed memristor model with mobile charged vacancies, the stationary endpoint concentrations and limiting resistances are independent of the vacancy-vacancy electrostatic interaction strength.","lead":"This paper derives analytical stationary vacancy distributions for a one-dimensional memristor model that includes electrostatic repulsion between charged vacancies. It shows that the limiting on and off resistances do not depend on the strength of that interaction, at least in the leading linear order.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Printed quadrature Eq. (9) is inconsistent with Eq. (8) because it omits a square root; the exact-solution formula is wrong as written, even though the alpha-independence result for endpoint concentrations is sound.","rationale":"The reader's verdict was CONDITIONAL, citing the missing square root in Eq. (9), the omitted derivation of Eq. (12), and the unproved bound 0<υ<1, while regarding the alpha-independence claim as well supported. My stress-test agrees that Eq. (9) is wrong as printed: differentiating the printed form contradicts Eq. (8), and even the α=0 limit fails unless p=1. This is a genuine load-bearing issue for the paper's 'exact solution' contribution and for the quantitative curves in Figs. 1 and 2. However, it is a repairable typographical/algebraic error rather than a fatal flaw: the corrected quadrature with the square root is consistent. More importantly, the alpha-independence of c(0) and c(1) can be proven directly from Eq. (8) and its first integral without using Eq. (9) at all. Integrating Eq. (8) yields Λ_x(1)=Λ_x(0); the first integral yields Λ_x² - αV = const; hence V(Λ0)=V(Λ0+p). Combined with Λ(1)-Λ(0)=p, this gives Eq. (12) independently of α. Thus the central physical claim identified by the reader is secure. The ideal-contact boundary condition Φ(0)=Φ(1)=0 is an explicit modeling assumption; it limits the scope to ideal contacts but does not create an internal inconsistency. I therefore see no reason to move the verdict. The CONDITIONAL verdict stands, primarily because the printed exact solution needs correction and the derivation of Eq. (12) should be stated explicitly, but the load-bearing physical conclusion is not in jeopardy.","tokens_in":12320,"tokens_out":13712,"duration_ms":123362,"concrete_test":"Set α=0 in Eqs. (9) and (10) with υ=1. The printed Eq. (9) gives Λ(x)=Λ0+A x, so the boundary condition Φ(1)=Φ(0)=0 forces A=p, whereas Eq. (10) gives A=p²; the two agree only for p=1. Replacing the integrand by 1/√(A+αV) removes the contradiction and reproduces the Burgers profile. Separately, verify that V(Λ0)=V(Λ0+p) with V(t)=log(1+e^t)-r t and Λ(1)-Λ(0)=p yields Eq. (12) with no α dependence, confirming that the endpoint-independence claim is robust despite the quadrature typo.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative contribution is the exact stationary solution, embodied in the quadrature Eq. (9). As printed, Eq. (9) cannot be correct. Let V(t)=log(1+e^t)-r t. If x = ∫_{Λ0}^{Λ} dt/(A+αV(t)), then dx/dΛ = 1/(A+αV), so Λ_x = A+αV, and differentiating again gives Λ_xx = α(c-r)(A+αV). Equation (8) requires Λ_xx = (α/2)(c-r), which would force A+αV = 1/2 identically, impossible. The α=0 limit exposes the same error explicitly: Eq. (9) yields Λ = Λ0 + A x, and the boundary condition Φ(1)=Φ(0)=0 gives Λ(1)-Λ(0)=p, so A=p, while Eq. (10) with υ=1 gives A=p², a contradiction for p≠1. The correct quadrature is ∫ dΛ/√(A+αV) = x, consistent with Eq. (10) and with the first integral of Eq. (8). This typo undermines the claimed analytical expressions, the numerical determination of υ from substituting x=1 into Eq. (9), the profiles in Fig. 1, and the asymptotics in Appendix A. By contrast, the headline alpha-independence claim does not depend on the quadrature: integrating Eq. (8) gives Λ_x(1)-Λ_x(0)=0, and the first integral Λ_x² - αV = const gives V(Λ0)=V(Λ0+p); together with Λ(1)-Λ(0)=p this yields Eq. (12) for Λ0 with no α dependence. That argument is exact and internally sound. The ideal-contact boundary condition is a stated modeling assumption that limits scope, not an inconsistency.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a one-dimensional closed memristor with mobile charged vacancies, adding electrostatic interaction to the authors' earlier nonlinear model. It derives stationary-state equations for the vacancy concentration and the electrostatic potential, introduces the logit variable Λ, and reduces the problem to an autonomous second-order ODE. The authors claim an exact quadrature solution, give an explicit expression for the endpoint logit Λ0, and conclude that the endpoint concentrations c(0) and c(1), and hence the limiting resistances in the on and off states, are independent of the electrostatic interaction strength α in the linear approximation (3). Figures show vacancy profiles for α=0 and α=1000, and an appendix provides asymptotic and rational approximations for a derivative scaling factor υ.","tokens_in":12653,"tokens_out":7025,"duration_ms":57335,"significance":"If correct, the alpha-independence of the limiting resistances is a nontrivial exact result that is robust to the strength of the electrostatic interaction, and the analytical stationary profiles are a useful contribution to memristor modeling. A strength of the derivation is that the independence of the endpoint concentrations follows directly from the first integral of Eq. (8) and the boundary conditions, without relying on the printed quadrature. However, the printed central quadrature contains an error that affects the quantitative results, and some supporting derivations are omitted, so the manuscript requires revision before the claims can be fully accepted.","major_comments":[{"comment":"The printed quadrature is dimensionally inconsistent with Eq. (8) and is missing a square root in the denominator. Differentiating the printed Eq. (9) gives Λ_x = A + αV(Λ) with V(t) = log(1+e^t)−rt, and differentiating again yields Λ_xx = α(c−r)(A+αV(Λ)), which can agree with Eq. (8) only if A+αV(Λ) = 1/2 identically, which is impossible. The α=0 limit of Eq. (9) gives A = p from the boundary condition Λ(1)−Λ(0)=p, whereas Eq. (10) with υ=1 gives A = p², a direct contradiction. The correct first integral is Λ_x² − αV(Λ) = A, so the quadrature should be ∫_{Λ0}^{Λ} dt / sqrt(A + αV(t)) = x. This error undermines the numerical determination of υ from x=1, the profiles in Fig. 1, and the asymptotics in Appendix A, all of which must be recomputed with the corrected formula.","section":"Section III, Eq. (9)"},{"comment":"The central claim that Λ0, and hence c(0) and c(1), are independent of α is based entirely on Eq. (12), but the derivation of Eq. (12) from the boundary conditions and particle-number conservation is not shown in the manuscript. Since this is the load-bearing step for the main result, a short derivation should be provided so that the reader can verify the claimed independence.","section":"Section III, Eq. (12)"}],"minor_comments":[{"comment":"The parameter υ is defined analytically in Eq. (11) as the ratio of derivatives at c=r, but the caption of Fig. 1 defines it geometrically as υ = tanφ/tanφ_B without introducing this geometric definition in the text; these definitions should be reconciled.","section":"Section III, Eq. (11) and Fig. 1 caption"},{"comment":"The statement 'It can be shown that always 0 < υ < 1' is made without proof or reference; a proof or a citation should be provided.","section":"Section III, after Eq. (11)"},{"comment":"Equation (A5) contains a misplaced parenthesis in 'h(576 +ap(1 + (r−1)r)))' and the notation 'a/b=' is ambiguous; write the expressions for a and b separately.","section":"Appendix A, Eq. (A5)"},{"comment":"The full Russian translation of the manuscript is appended after the English references, which is unconventional for a journal submission and should be removed or moved to supplementary material.","section":"General"},{"comment":"Please specify that all logarithms are natural logarithms, and clarify the notation 'cx|c=r' and 'bx|b=r' in Eq. (11) as derivatives with respect to x evaluated where c=r and b=r, respectively.","section":"Throughout"},{"comment":"The sentence 'It seems (although the authors did not emphasize this) ...' refers to the authors in the third person; rephrase to 'we did not emphasize' or use another suitable first-person construction.","section":"Section III, paragraph on transition region"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is an incremental extension of the authors' prior memristor model, but the central alpha-independence result is nontrivial and potentially interesting. The Russian duplicate of the full text is unusual for an arXiv submission and should be removed or made supplementary. The main technical error in Eq. (9) is a fixable typo, but it affects several quantitative claims, so major revision is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the headline physical claim is right, and the centerpiece formula as printed is wrong. The stress-test note is correct about Eq. (9). With V(t)=log(1+e^t)-rt, the printed form x = integral of dt/(A+alpha V) gives Lambda_x = A+alpha V, so Lambda_xx = alpha(c-r)(A+alpha V), which cannot match Eq. (8)'s (alpha/2)(c-r) unless A+alpha V = 1/2 identically. The alpha=0 limit makes the same point: Eq. (9) would give Lambda = Lambda0 + A x, and the boundary conditions force A=p, while Eq. (10) at upsilon=1 gives A=p^2. The missing square root is not cosmetic; it undermines the profiles in Fig. 1, the numerical determination of upsilon, and the asymptotics in Appendix A as printed.\n\nWhat is actually new survives this. The model extends the authors' earlier Burgers-equation memristor by adding the Poisson coupling, and the stationary profiles plus the intermediate neutral region are new. More importantly, the endpoint-independence result is real. It does not use the quadrature at all: integrating Eq. (8) gives Lambda_x(1)=Lambda_x(0), the first integral gives Lambda_x^2 - alpha V = const, and with Lambda(1)-Lambda(0)=p this yields Lambda0 independent of alpha. That argument is exact and internally sound, and it is the paper's most valuable contribution. The paper is also honest about the ideal-contact boundary condition as a modeling assumption; the scope limitation is stated, not hidden.\n\nSoft spots, in proportion: Eq. (9) is a load-bearing typo that must be corrected, but it is repairable and the central claim does not depend on it. The bound 0 < upsilon < 1 is asserted rather than proved, and the derivation of Eq. (12) is omitted; both should be supplied. The \"exact solution\" still needs a numerical root-find for upsilon, so it is exact in the quadrature sense, not fully closed-form; that is a fair characterization but worth stating precisely. The flux model in Eq. (1) is carried over from prior work without microscopic derivation; I do not see that as a flaw here since this paper is about the stationary consequences of that model.\n\nWho is this for? People working on analytical memristor device theory, especially on vacancy drift-diffusion models. It will not reshape the field, but the endpoint-independence result is a clean, surprising simplification that deserves to be in the literature. I would send it to peer review with a directive to fix Eq. (9) and supply the missing derivations; after that, I would accept it. I would also cite it if I worked on these models.","headline":"The alpha-independence claim is correct and survives the typo in Eq. (9), but the printed quadrature is broken and must be fixed before this is citable.","tokens_in":13213,"tokens_out":2087,"would_cite":true,"duration_ms":19307,"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":"For a closed memristor with mobile charged vacancies, the stationary end concentrations and the limiting on/off resistances are independent of the strength of electrostatic interaction among the vacancies, even though the interior vacancy…","keywords":["memristor","charged vacancies","electrostatic interaction","stationary states","Poisson equation","Burgers equation","logit transform","resistance switching"],"falsifier":"Compute stationary solutions of Eqs. (4) with the potential boundary conditions changed to $\\varphi(0)=V_0$ and $\\varphi(1)=V_1$, or to a Robin condition representing contact capacitance, holding $p$, $r$, and $\\alpha$ fixed, and check whether $\\Lambda_0$ and the endpoint concentrations stay independent of $\\alpha$; any $\\alpha$ variation would falsify the claim for that boundary class. Experimentally, one could vary $\\alpha$ through temperature or through the dielectric constant of the host while measuring $R_{\\rm on}$ and $R_{\\rm off}$ in a closed memristor with high-quality contacts: the claim predicts no change in the limiting resistances.","tokens_in":12076,"feed_emoji":"⚡","tokens_out":7540,"duration_ms":66845,"temperature":0.7,"pith_summary":"This paper extends a nonlinear model of a closed memristor—one that exchanges only electrons, not atoms—to include electrostatic repulsion among mobile charged vacancies, and asks what that interaction does to the stationary vacancy profiles that set the device's on and off resistances. The authors obtain exact analytical expressions for these stationary states, showing that the interaction reshapes the interior profile and creates an intermediate electrically neutral region between the vacancy-rich and vacancy-poor zones. Despite this strong reshaping, the vacancy concentrations at the two ends, $c(0)$ and $c(1)$, come out independent of the interaction strength $\\alpha$, and so do the limiting on/off resistances computed in the linear approximation. The reason is that the logit boundary constant $\\Lambda_0$ is fixed by the potential boundary conditions and vacancy conservation alone, with no $\\alpha$ dependence. If this is right, the static resistance contrast of such a memristor is a boundary and contact property, while the Coulomb interaction only rearranges the internal charge distribution.","feed_headline":"Electrostatic repulsion leaves memristor on/off resistance unchanged","feed_subtitle":"Vacancy-vacancy interactions reshape the profile but not the limiting resistances in this exact stationary solution.","key_machinery":"The load-bearing object is the logit transform $\\Lambda=\\log(c/(1-c))$, which linearizes the vacancy flux equation, reduces the stationary problem to the autonomous ODE $\\Lambda_{xx} = \\frac{\\alpha}{2}\\big(e^\\Lambda/(1+e^\\Lambda)-r\\big)$, and yields the solution as a quadrature with two integration constants $\\Lambda_0$ and $\\upsilon$. The ideal-contact potential conditions translate into $\\Lambda(1)-\\Lambda(0)=p$, and the explicit formula for $\\Lambda_0$ shows why $\\alpha$ drops out of the endpoint concentrations. The derivative-scaling factor $\\upsilon$, defined as the slope of the profile at $c=r$ normalized to the $\\alpha=0$ Burgers solution, carries the electrostatic reshaping of the interior profile.","core_discovery":"The paper's central claim is that in the stationary state of a closed memristor with mobile charged vacancies, the endpoint concentrations $c(0)$ and $c(1)$ are independent of the electrostatic interaction strength $\\alpha$, and therefore the limiting on and off state resistances are also $\\alpha$-independent to leading linear order. The argument runs through the logit transformation $\\Lambda=\\log(c/(1-c))$, which turns the zero-flux condition into $\\Lambda = p x - \\alpha\\varphi/2 + \\Lambda_0$ and the Poisson equation into an autonomous second-order ODE solvable by quadrature. The ideal-contact conditions $\\varphi(0)=\\varphi(1)=0$ force $\\Lambda(1)-\\Lambda(0)=p$, and conservation of the total vacancy number then fixes $\\Lambda_0 = \\log\\big((e^{pr}-1)/(e^p-e^{pr})\\big)$, which contains no $\\alpha$. Consequently the end concentrations, and hence the resistance combination $\\sigma = (c(0)/C_{\\max})\\cos^2\\theta + (c(1)/C_{\\max})\\sin^2\\theta$, do not change when the Coulomb interaction is turned on, even though the full profile develops an intermediate neutral region and the interior slope at the filling point is reduced by a factor $\\upsilon$ that lies strictly between 0 and 1.","pith_inferences":["My inference: if the $\\alpha$-independence is understood as a boundary effect, then engineering the contacts—changing $\\theta$, $R_1$, or $R_2$—is the only static lever on the on/off resistance ratio, while tuning the dielectric environment should affect switching dynamics but not the limiting resistance values.","My expectation: repeating the stationary calculation with non-ideal contacts, for example finite surface capacitance or a Schottky barrier replacing $\\varphi(0)=\\varphi(1)=0$, will make $c(0)$ and $c(1)$ depend on $\\alpha$; the reported cancellation is specific to ideal contacts.","A testable extension: the same quadrature method could be adapted to slow time-dependent switching around the stationary state, connecting the $\\alpha$-independent endpoints to the switching kinetics studied in the earlier linear-resistance model."],"forward_implications":["The static on/off resistance contrast of such a closed memristor is fixed by the contact asymmetry and the total vacancy filling, not by the strength of vacancy-vacancy repulsion.","Strong Coulomb interaction can be diagnosed through the interior profile: a neutral plateau forms between the depleted and enriched zones, and its width grows with $\\alpha$.","The interior slope at the filling point is always reduced by electrostatic interaction ($0<\\upsilon<1$), and it approaches the Burgers value in the large-current limit $p\\to\\infty$.","The quadrature solution plus the rational approximation for $\\upsilon$ give analytic estimates for stationary profiles over the whole parameter range without full numerics."],"supporting_citations":[{"why":"Supplies the nonlinear vacancy-flux model (thermal drift with saturation) and the Burgers-equation stationary solution that this paper generalizes by adding electrostatics.","marker":"[18]"},{"why":"Provides the linear-order resistance formula, Eq. (3), through which the endpoint concentrations determine the on and off state resistances.","marker":"[19]"},{"why":"Reports the first TiO2 memristor realization and the widely used two-phase model that the nonlinear model is contrasted with.","marker":"[3]"}],"fun_headline_variants":["Memristor on/off resistance immune to vacancy repulsion","Coulomb repulsion fails to shift memristor limit resistances","Stationary memristor state defies electrostatic interaction","Vacancy interactions alter profile, not on/off resistance"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire $\\alpha$-independence result rests on ideal contacts that pin the electrostatic potential to zero at both ends, together with the global-neutrality background; if real contacts hold a finite surface charge, the cancellation no longer follows.","fun_headline_variants_meta":{"raw":{"variants":["Memristor on/off resistance immune to vacancy repulsion","Coulomb repulsion fails to shift memristor limit resistances","Stationary memristor state defies electrostatic interaction","Vacancy interactions alter profile, not on/off resistance"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000182,"raw_usage":{"total_tokens":1302,"prompt_tokens":926,"completion_tokens":376,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":542,"completion_tokens_details":{"reasoning_tokens":307}},"tokens_in":542,"tokens_out":376,"duration_ms":3203,"temperature":1.0,"reasoning_tokens":307,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:10:47.829842+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute stationary solutions of Eqs. (4) with the potential boundary conditions changed to $\\varphi(0)=V_0$ and $\\varphi(1)=V_1$, or to a Robin condition representing contact capacitance, holding $p$, $r$, and $\\alpha$ fixed, and check whether $\\Lambda_0$ and the endpoint concentrations stay independent of $\\alpha$; any $\\alpha$ variation would falsify the claim for that boundary class. Experimentally, one could vary $\\alpha$ through temperature or through the dielectric constant of the host while measuring $R_{\\rm on}$ and $R_{\\rm off}$ in a closed memristor with high-quality contacts: the claim predicts no change in the limiting resistances.","supporting_citations":[{"cited_title":"Boilo and K","cited_arxiv_id":null,"evidence_quote":"Provides the linear-order resistance formula, Eq. (3), through which the endpoint concentrations determine the on and off state resistances."}],"review_version":1}