{"id":"5bfa2923-1dd9-40a9-bc3e-4548d5291284","arxiv_id":"2501.12070","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A Hamiltonian and open-quantum-systems framework is constructed for dissipative polarizable media, unifying classical polarizable models with quantum dynamics.","lead":"This paper builds a quantum description of polarizable media by turning the classical equations of a damped harmonic oscillator into a non-Hermitian Hamiltonian with auxiliary variables. It then derives a self-consistent semiclassical electric field and a master equation for a quantum system coupled to such a medium, with applications to plasmonic nanostructures.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The master equation for the QPM is derived by grafting the Breuer-Petruccione construction onto a non-Hermitian bath, but the rate matrix gamma(omega) is never shown to be positive semidefinite; without that, the dissipator in Eq.","rationale":"The reader's weakest_assumption identifies exactly the most load-bearing concern: the master-equation derivation assumes the standard Born-Markov-Lindblad machinery applies to a non-Hermitian environment Hamiltonian H0, but the required positivity of the rate matrix gamma(omega) is not established. This is load-bearing because the master equation (Eq. 95-96) is one of the three principal results advertised in the abstract, and if gamma is not positive semidefinite, the dissipator is not completely positive and the reduced dynamics is not a physical quantum channel. The paper itself concedes the lack of Hermitian symmetries in the correlation functions, so the gap is not an artifact of the referee; it is a genuine missing proof. A concrete one-dimensional test can settle the issue numerically. I also considered the initial-condition subtlety in Eqs. (7b)-(7c), where the stated dot x(0) omits a possible -f(0) term, making the first-order system exactly equivalent to Eq. (1) only when f(0)=0; however, for typical physical driving (fields switched on after t=0) this is a secondary concern, whereas the master equation positivity gap directly undermines the paper's central open-quantum-system claim. The reader's conditional verdict is therefore appropriate: the framework is promising, but the master equation, and hence the 'open quantum dynamics in QPM' headline, is not yet established. No change to the verdict is needed; the conditionality already flags this issue.","tokens_in":36155,"tokens_out":10682,"duration_ms":111926,"concrete_test":"Take the simplest QPM: n = 1, K = omega_0^2, Gamma = gamma, with A = 1 (as in the symmetric, undamped-limit case). Compute the 2x2 environment correlation matrix from Eq. (109) (or Eq. (106)) for a thermal Gaussian state at inverse temperature beta. Numerically evaluate gamma_{alpha beta}(omega) via Eq. (98) for a range of omega, e.g., omega in [-5 omega_0, 5 omega_0]. Check whether the 2x2 matrix gamma(omega) is Hermitian and positive semidefinite at every omega. Then propagate a two-level system with the dissipator in Eq. (96) and verify whether the reduced density matrix remains positive semidefinite. If any eigenvalue of gamma is negative, or if rho_S(t) develops a negative eigenvalue, the master equation is not completely positive and the central open-system claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section V.B, the paper writes the dissipator in Lindblad form, Eq. (96), with rates gamma_{alpha beta}(omega) = Xi_{alpha beta}(omega) + Xi*_{beta alpha}(omega), Eq. (98), where Xi is a half-Fourier transform of <x_alpha(t) x_beta(0)>. For a Hermitian bath, x†(t) = x(t) and stationarity make gamma the full Fourier transform of a positive-definite kernel, so gamma(omega) is positive semidefinite and the GKSL theorem applies. Here H0 is non-Hermitian, so x(t) is not x†(t); the paper explicitly states 'the correlation functions will lack the symmetries present in the Hermitian case' but does not prove that gamma(omega) is Hermitian positive semidefinite. If gamma has a negative eigenvalue, the dissipator is not completely positive and can generate unphysical (negative-probability) reduced states. The QPM also has a finite number of modes, so the correlation functions are finite sums of exponentials and do not decay; the Markovian limit invoked by using the Breuer-Petruccione result is not justified for a finite environment. Since the master equation is one of the three headline results and the only concrete application to open quantum systems, this gap is load-bearing: the claimed 'master equation describing the open dynamics' may not describe any physical open quantum system. The paper's own note in Section V.B acknowledges the asymmetry but stops short of establishing the necessary positivity condition.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a Hamiltonian and quantum-mechanical formulation of dissipative classical polarizable media governed by Eq. (1). It introduces auxiliary variables v=u˙, a matrix √K whose nonzero spectrum matches the poles of (ω^2+2iωΓ−K)^{-1}, and a quadratic Hamiltonian (Eq. 11) that reproduces Eq. (1) as Ehrenfest mean dynamics after canonical quantization. The resulting Quantum Polarizable Medium (QPM) is analyzed via pseudo-boson operators and bi-coherent states (Sec. II.C), applied to reproduce the polarizability spectra of Ag147 and graphene-disk nanostructures using eigenpair filtering (Sec. III.B), and used to derive a self-consistent semiclassical field equation (Eq. 84) and a Born-Markov master equation for a quantum system embedded in the QPM (Eqs. 95-98), with correlation functions for Gaussian states (Eqs. 104-110). The paper also derives algebraic decompositions of the polarizability in terms of the eigenvalues and eigenvectors of √K.","tokens_in":36549,"tokens_out":6324,"duration_ms":64985,"significance":"The paper is ambitious and addresses a genuine need: putting classical polarizable embedding models on a quantum footing for open-system simulations. The spectral decomposition of the polarizability in terms of √K is original and potentially useful, and the phase-space treatment of quadratic non-Hermitian Hamiltonians is carefully developed in the appendices. The authors are explicit about several limitations, including the unresolved gauge dependence of the interaction term and the loss of Hermitian correlation-function symmetries. However, the central open-system result, the master equation, rests on an unproven extension of the standard GKSL derivation to a non-Hermitian finite bath, and the thermal-state covariance formula is used without proof of validity for the indefinite B that appears here. These gaps are load-bearing rather than presentational. If the positivity of the rates and the validity of the thermal state can be established, the framework would be a substantial contribution to the theory of dissipative quantum polarizable media.","major_comments":[{"comment":"The derivation of the master equation assumes that the standard Born-Markov-Lindblad construction (Ref. [114]) extends to a non-Hermitian environment Hamiltonian H0. The paper defines the rates γαβ(ω)=Ξαβ(ω)+Ξ*βα(ω) in Eq. (98) but does not prove that the matrix γ(ω) is positive semidefinite. For a Hermitian bath, this follows from stationarity and x†(t)=x(t), which makes γ the full Fourier transform of a positive-definite kernel. The paper explicitly states before Sec. V.A that \"the correlation functions will lack the symmetries present in the Hermitian case,\" but it does not address the positivity consequence. Without positivity, the dissipator (96) is not in Lindblad form and can generate non-positive reduced states. Additionally, the QPM has a finite number of modes, so the correlation functions in Eq. (104) are finite sums of exponentials and do not decay; the Markovian limit invoked by the Breuer-Petruccione construction is not justified for such a finite environment, and the \"limit η→0+\" in Eq. (106) is a distributional identity rather than a physical Markov approximation. This is load-bearing because the master equation is one of the three headline results and the only concrete open-system application.","section":"Section V.B, Eqs. (95)-(99)"},{"comment":"The thermal Gaussian state covariance M0 = -ℏ/2 cot(ℏβ J B /2) J^T is used as an input for the correlation functions in Eqs. (108)-(110). In the standard phase-space derivation, this formula requires the quadratic Hamiltonian to be Hermitian and the matrix JB to have appropriate spectral properties so that M0 is a legitimate covariance matrix (positive definite). Here B is only symmetric and JB is generally non-Hermitian, and the paper provides no proof that the coth formula yields a valid Gaussian state with a positive Wigner function (Eq. 61). If M0 fails to be positive definite, the \"arbitrary Gaussian states\" claim and the correlation functions derived from it are not physically meaningful. A proof of validity under explicit conditions on K and Γ, or a restriction to parameter regimes where the formula is known to hold, is needed.","section":"Section IV.A, Eq. (107)"},{"comment":"The Hamiltonian quantization is defined only up to the similarity matrix A, and the paper states that \"the situation is, however, less clear for the interaction term H1(t), an issue that will remain open in this work.\" This unresolved gauge ambiguity propagates into the later results that use H1(t): the self-consistent electric field equation (84) and the open-system interaction Hamiltonian in Sec. V.A (Eq. 88) both depend on the interaction term. Furthermore, the claim that H0 is identically zero on-shell relies on the complex relation π(t)=iA^{-1}√K x(t); for a non-Hermitian H0, the physical meaning of \"zero on-shell\" as an energy is not established. The authors acknowledge the H1 issue, but it should either be resolved or explicitly stated as a condition restricting the validity of the subsequent derivations.","section":"Section II.B, Eq. (11)"}],"minor_comments":[{"comment":"The text describes a \"spectral equivalence\" between √K and Eq. (1), but Eq. (10) states an inclusion Sp{−√K} ⊆ {ω : det(ω^2+2iωΓ−K)=0} ∪ {0}. Since the applications rely on the eigenvalues of √K for filtering, please clarify whether equality is intended and how zero eigenvalues are handled.","section":"Section II.A, Eq. (10)"},{"comment":"The red vertical lines in Fig. 2 are described as the real part of the eigenvalues of √K, but the caption does not state the normalization or the binning/density of the eigenvalue distribution; please clarify.","section":"Section III.B, Fig. 2 caption"},{"comment":"The dissipator is written with a summation over α,β=1,...,N, but the operators O_α(ω) are defined for α up to 2n=N in Eq. (97); please ensure the index ranges are consistent.","section":"Section V.B, Eq. (96)"},{"comment":"Many key derivations are relegated to the Supplemental Material (e.g., Refs. [83], [98], [105], [111]-[113]); as a standalone journal article, at least the main steps of these derivations should be summarized in the main text or the supplement made openly available.","section":"General"},{"comment":"The notation for the tensor product in Eq. (36)-(37) alternates between f(ω)⊗R_s and f(ω)R_s^T; please unify the notation and explicitly define the dimensions of the matrices involved.","section":"Section III.A, Eqs. (36)-(37)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious theoretical contribution with a genuine need and a well-developed phase-space formalism. The main obstacle to acceptance is the master-equation derivation in Sec. V.B: the positivity of the rate matrix γ(ω) is not established, and the Markovian approximation for a finite non-Hermitian bath is not justified. I recommend requesting a rigorous positivity proof or an explicit restriction to parameter regimes where the Lindblad form holds, and similarly a justification for the thermal covariance formula in Eq. (107). If these gaps can be closed, the paper would be suitable for publication in a journal that welcomes non-Hermitian and open-quantum-system methods."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The main thing to know: this paper is a real attempt to give a finite-dimensional quantum Hamiltonian formulation for dissipative polarizable media, and it does a lot of work well. The lift to sqrt(K), the pseudo-boson diagonalization, and the phase-space derivation of a self-consistent emitted field are competently done. The spectral filtering analysis on realistic Ag147 and graphene disk models is a useful practical tool. I also credit the authors for flagging some open issues themselves, notably the gauge dependence of the interaction term H1(t).\n\nThe soft spots are real, and one is load-bearing. The master equation in Section V.B is derived by grafting the Breuer-Petruccione construction onto a non-Hermitian bath. The rates gamma_{alpha beta}(omega) = Xi_{alpha beta}(omega) + Xi*_{beta alpha}(omega) are never shown to be Hermitian positive semidefinite. Because x(t) is not x^dagger(t) under non-Hermitian H0, the standard symmetries fail, and the paper explicitly says the correlation functions will lack the Hermitian symmetries, but it does not provide the necessary positivity check. Without that, the dissipator may not be completely positive, and the claimed master equation may not describe any physical open quantum system. That is not a footnote; it is one of the three headline results. The Markovian limit is also assumed for a finite set of modes whose correlation functions are finite sums of exponentials and do not decay, so the time-scale separation is not justified. The H1 gauge issue is secondary by comparison.\n\nThe applications are honest but somewhat circular: the QPM is constructed so its mean dynamics coincide with the classical equation of motion, and the spectra plotted re-express the same classical response functions. That is a consistency check, not new predictive content. No code or data are shipped, so independent reproduction of the figures is not possible.\n\nWho gets value? Theorists working on non-Hermitian quantum mechanics and open quantum systems, and quantum-chemistry researchers who want a formal bridge between classical polarizable models and genuine quantum environments. It deserves a serious referee: the program is worth engaging, but the master-equation section needs substantial revision before the claims stand. I would send it to peer review with a clear request that the authors either prove complete positivity (or state conditions under which it holds) or present the master equation as a formal, unproven construction.","headline":"A serious non-Hermitian quantization of classical polarizable media, but the derived master equation lacks a complete-positivity proof and the applications re-derive known spectra.","tokens_in":37065,"tokens_out":2129,"would_cite":false,"duration_ms":24885,"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":"The paper shows that the classical equation of motion for dissipative polarizable media can be quantized exactly as a finite quadratic Hamiltonian, reproducing the classical dynamics as Ehrenfest means and yielding a master equation for…","keywords":["quantum polarizable medium","damped harmonic oscillator","pseudo-bosons","open quantum systems","master equation","plasmonics","polarizability","Gaussian states"],"falsifier":"Compute the rate matrix $\\gamma_{\\alpha\\beta}(\\omega)$ from Eq. (98) for the simplest QPM, a single damped harmonic oscillator with $\\Gamma>0$ in a thermal state, and check whether $\\gamma(\\omega)$ is positive semidefinite for all $\\omega$; a negative eigenvalue anywhere would show that the derived dissipator is not of the standard completely positive form, so the master equation could not describe a physical open system.","tokens_in":35919,"feed_emoji":"⚛️","tokens_out":12540,"duration_ms":123149,"temperature":0.7,"pith_summary":"The paper sets out to show that the dissipative polarizable media used throughout computational chemistry and plasmonics, normally treated as classical equations of motion with friction, are the mean-field limit of a genuine quantum system with the same finite number of degrees of freedom. It constructs a quadratic Hamiltonian in an extended phase space, quantizes it under the condition that expectation values follow the classical damped equation, and calls the result the Quantum Polarizable Medium (QPM). From there it derives a spectral decomposition of the medium's polarizability, a self-consistent semiclassical equation for the re-emitted electric field, and a master equation for a quantum emitter embedded in the medium, with correlation functions valid for arbitrary Gaussian states. If these derivations are right, standard polarizable models inherit a fully quantum description of decoherence and energy exchange without invoking an infinite bath or external noise.","feed_headline":"Damped polarizable media quantized without an infinite bath","feed_subtitle":"Classical damping becomes Ehrenfest means of a finite quantum system, yielding field formulas and a master equation.","key_machinery":"The load-bearing object is the matrix square root of the extended dynamical matrix, $\\sqrt{\\mathcal{K}} = i\\begin{pmatrix}0&-I\\\\K&2\\Gamma\\end{pmatrix}$, whose eigenvalues coincide, up to a possible null mode, with the poles of the polarizability kernel $(\\omega^2 + 2i\\omega\\Gamma - K)^{-1}$. This matrix converts the second-order damped equation into a first-order system and supplies the spectrum used to decompose absorption and dispersion. On top of it the paper builds a quadratic Hamiltonian $H = \\frac{1}{2}\\pi^T A\\pi + \\frac{1}{2}x^T A^{-1}\\mathcal{K}x + x^T A^{-1}F(t)$, with $A$ the symmetric matrix realizing $\\mathcal{K} = A\\mathcal{K}^T A^{-1}$; quantization of $x,\\pi$ under the canonical commutation relation makes the classical equation the Ehrenfest limit. Pseudo-boson operators $b,\\tilde b$ diagonalize $H_0$ and generate bi-coherent states, while the symplectic evolution $\\Lambda_t = e^{JBt}$ in phase space carries the Gaussian-state calculations that lead to the self-consistent field equation and the correlation functions.","core_discovery":"The central claim is that the damped oscillator equation $\\ddot u + 2\\Gamma\\dot u + Ku + f=0$, which underlies fluctuating-charge, fluctuating-dipole, discrete-interaction, and continuum polarizable models, can be embedded in an extended phase space whose quadratic Hamiltonian $H = \\frac{1}{2}\\pi^T A\\pi + \\frac{1}{2}x^T A^{-1}\\mathcal{K}x + x^T A^{-1}F(t)$ is finite and time-independent when no external field acts. Quantizing this Hamiltonian under the canonical commutation relations defines the Quantum Polarizable Medium; the classical damped equation is recovered as the Ehrenfest mean-value dynamics, and the time-independent part $H_0$ is a constant of motion that is identically zero on-shell, so the isolated medium's energy comes only from the interaction term. The spectrum of the extended matrix $\\sqrt{\\mathcal{K}} = i\\begin{pmatrix}0&-I\\\\K&2\\Gamma\\end{pmatrix}$ is shown to coincide with the poles of the polarizability kernel $(\\omega^2 + 2i\\omega\\Gamma - K)^{-1}$, and $H_0$ is diagonalized by pseudo-boson operators whose bi-coherent states evolve by a simple phase-space rule. On this basis the paper obtains a self-consistent semiclassical equation for the emitted electric field and a Markovian master equation for an embedded quantum system, with correlation functions computed over arbitrary Gaussian states.","pith_inferences":["The paper does not check whether the rate matrix $\\gamma_{\\alpha\\beta}(\\omega)$ is positive semidefinite; a single-oscillator calculation would settle whether the derived master equation is physically admissible, a question the authors leave open.","A direct extension would use the same finite-dimensional non-Hermitian Hamiltonian as a reservoir in numerically exact methods such as quantum trajectories or hierarchical equations of motion, rather than stopping at the Markovian master equation; this would test the Markovian approximation when the environment itself has non-Hermitian spectral features.","The intercept filter suggests a practical model-order reduction: keeping only the eigenpairs with the largest intercept coefficients should reproduce time-domain polarization dynamics with far fewer variables than the full $\\sqrt{\\mathcal{K}}$ matrix, a numerical strategy the paper illustrates but does not propose as a general algorithm."],"forward_implications":["Any classical polarizable medium of the form $\\ddot u + 2\\Gamma\\dot u + Ku + f=0$ acquires a finite-dimensional quantum counterpart, so dissipation no longer requires an infinite bath or an added noise term.","The polarizability of a plasmonic structure decomposes into absorptive and dispersive eigenmode contributions of $\\sqrt{\\mathcal{K}}$; the intercept filter reconstructs the main spectral features with very few modes, up to 4 for the regular graphene disk.","The first-order solution of the self-consistent field equation expresses electric-field enhancement as a sum of Gaussian modes in $k$-space weighted by the resonances of $(\\omega E + iJB)^{-1}$, connecting the medium's quantum spectrum to near-field enhancement.","A quantum emitter coupled to the QPM obeys the derived master equation, with rates determined by the $K$ and $\\Gamma$ matrices and correlation functions for thermal or arbitrary Gaussian states.","Because $H_0$ vanishes on-shell, the average energy of the isolated medium is exactly $E(t)=\\mathrm{tr}\\{H_1(t)\\rho\\}$; only the interaction with external fields contributes to the energetics."],"supporting_citations":[{"why":"Establishes the open-quantum-system treatment of polarizable continuum models and supplies the $A_\\alpha u_\\alpha$ interaction form that Section V generalizes to auxiliary variables.","marker":"[32]"},{"why":"The reference for the standard Markovian master-equation derivation that Section V.B follows.","marker":"[114]"},{"why":"Provides the pseudo-boson operator framework used to diagonalize $H_0$ and construct bi-coherent states.","marker":"[87]"},{"why":"Supplies the pseudo-boson coherent-state and number-state formalism, including the non-orthogonality relations used in the dynamics.","marker":"[88]"},{"why":"Establishes the similarity transformation $\\mathcal{K} = A\\mathcal{K}^T A^{-1}$ whose symmetric matrix $A$ enters the Hamiltonian.","marker":"[54]"},{"why":"Earlier recognition of the square-root matrix $\\sqrt{\\mathcal{K}}$, the object at the center of the paper's spectral analysis.","marker":"[75]"},{"why":"Recent Lagrangian derivation of the same square-root matrix, which the paper extends to a Hamiltonian and quantum formulation.","marker":"[76]"},{"why":"Provides the identity $[q,H] = -i\\hbar J\\nabla H$ used to derive the Ehrenfest dynamics from the quantized Hamiltonian.","marker":"[84]"},{"why":"Gives the Wigner-function evolution for quadratic systems, the phase-space method used to compute Gaussian-state expectation values.","marker":"[100]"}],"fun_headline_variants":["No bath needed for damped quantum media","Damped media get finite quantum description","Quantizing dissipative media without infinite baths","Exact quantum fields from damped media","Finite quantum phase space for dissipative media"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument rests on the assumption that the standard open-quantum-system master-equation derivation remains valid when the environment Hamiltonian is non-Hermitian; if the correlation rates $\\gamma_{\\alpha\\beta}(\\omega)$ from Eq. (98) can be negative, the master equation is not completely positive and would not describe a physical system.","fun_headline_variants_meta":{"raw":{"variants":["No bath needed for damped quantum media","Damped media get finite quantum description","Quantizing dissipative media without infinite baths","Exact quantum fields from damped media","Finite quantum phase space for dissipative media"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000556,"raw_usage":{"total_tokens":2673,"prompt_tokens":1000,"completion_tokens":1673,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":616,"completion_tokens_details":{"reasoning_tokens":1607}},"tokens_in":616,"tokens_out":1673,"duration_ms":12693,"temperature":1.0,"reasoning_tokens":1607,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:32:57.896965+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the rate matrix $\\gamma_{\\alpha\\beta}(\\omega)$ from Eq. (98) for the simplest QPM, a single damped harmonic oscillator with $\\Gamma>0$ in a thermal state, and check whether $\\gamma(\\omega)$ is positive semidefinite for all $\\omega$; a negative eigenvalue anywhere would show that the derived dissipator is not of the standard completely positive form, so the master equation could not describe a physical open system.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The reference for the standard Markovian master-equation derivation that Section V.B follows."},{"cited_title":"Quantum dynamics in the self-consistent quadratic approximation","cited_arxiv_id":"2403.11327","evidence_quote":"Provides the pseudo-boson operator framework used to diagonalize $H_0$ and construct bi-coherent states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier recognition of the square-root matrix $\\sqrt{\\mathcal{K}}$, the object at the center of the paper's spectral analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Recent Lagrangian derivation of the same square-root matrix, which the paper extends to a Hamiltonian and quantum formulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the identity $[q,H] = -i\\hbar J\\nabla H$ used to derive the Ehrenfest dynamics from the quantized Hamiltonian."},{"cited_title":"[12], Eqs","cited_arxiv_id":null,"evidence_quote":"Gives the Wigner-function evolution for quadratic systems, the phase-space method used to compute Gaussian-state expectation values."}],"review_version":1}