{"id":"da7f10f9-cd40-41f4-93f6-1bd01debd1c2","arxiv_id":"2608.05968","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For two tunnel-coupled Luttinger liquids, entanglement dynamics is solved analytically in the Gaussian regime: early-time growth exponent is set by the protocol's first non-vanishing derivative, and long-time averages are exact at finite temperature.","lead":"This paper derives exact formulas for how quantum entanglement grows between two coupled one-dimensional quantum fluids when a tunneling interaction is switched on over time. It predicts that the early growth rate is set by the first non-zero derivative of the switching protocol, and it gives exact long-time averages at finite temperature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Gaussian expansion of the tunneling term is least controlled exactly where the early-time power law is derived (J(t) near zero), so the headline scaling may not describe the physical system.","rationale":"The reader's weakest assumption was the harmonic approximation of the tunneling term, and that is exactly the load-bearing concern here. Within the Gaussian model, the derivations appear correct: the Ermakov factorization, the symplectic eigenvalue formulas, the early-time expansion, and the long-time averages are mutually consistent and numerically cross-checked. The paper is honest about the approximation's limits in Sec. VII, but it does not establish validity in the small-J regime that controls the early-time power law. Because the central new claim is a universal exponent derived from a Taylor expansion around J=0, the physical relevance of the result is not yet settled. The reader's CONDITIONAL verdict is appropriate: the paper should be accepted only if the Gaussian predictions are benchmarked against the full sine-Gordon dynamics or the claim is explicitly restricted to the Gaussian model as a mathematical statement. No internal inconsistency was found, and the zero-mode issue is less severe than the harmonic approximation. Therefore the stress-test pass does not change the reader's verdict.","tokens_in":55100,"tokens_out":20741,"duration_ms":226262,"concrete_test":"Simulate the full time-dependent sine-Gordon Hamiltonian with a truncated-Wigner or MPS method for the linear-ramp protocol of Fig. 3 using the same experimental parameters, and extract the early-time growth of mutual information in the window t/t*_1 << 1. Compare the fitted exponent and prefactor with Eq. (43). If the full model does not reproduce I(A:B) ~ t^2 with the predicted amplitude, or if non-Gaussian correlations make the leading growth differ, the early-time scaling law is not a physical consequence of the tunnel-coupled Luttinger liquid but only of the Gaussian approximation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The analytical machinery inside the Gaussian model is internally consistent: the Ermakov solution, the symplectic-eigenvalue formulas of Sec. III, and the long-time averages of Sec. V all follow from the stated harmonic Hamiltonian. The load-bearing risk is the mapping from the physical tunnel-coupled Luttinger liquids to that Hamiltonian. Eq. (6) expands cos(phi_A - phi_B) to quadratic order around a phase minimum, which is justified only when the tunneling-induced mass is large enough to suppress relative-phase fluctuations. The new early-time result, I(A:B) ~ t^{n+1} in Eq. (43), is derived by Taylor-expanding J(t) around t=0, i.e. in the window where J(t) is arbitrarily small and the relative phase is essentially unpinned. This is the regime where the harmonic approximation is least controlled, not merely the adiabatic slow-ramp regime acknowledged in Sec. VII. Consequently, the universal exponent fixed by the first non-null derivative may be an artifact of the quadratic expansion rather than a property of the sine-Gordon dynamics. The zero-mode omission is secondary: it can shift the result by O(1) or O(log L), but it does not threaten the extensive thermodynamic-limit scaling. The central claim therefore needs an external validity check against the full time-dependent sine-Gordon model before it can be presented as a physical prediction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a Gaussian-level theory of entanglement dynamics between two identical Luttinger liquids coupled by a time-dependent tunneling amplitude J(t). After a harmonic expansion of the cosine tunneling term, the model decouples into a gapless symmetric sector and a time-dependent gapped antisymmetric sector, with all mode dynamics encoded in solutions gamma_q(t) of an Ermakov equation. The authors derive analytical expressions for mutual information, logarithmic negativity, and Rényi entropies in terms of gamma_q(t), and from these obtain three main results: an early-time power law I(A:B) ~ t^{n+1} whose exponent is fixed by the first nonvanishing derivative of J(t) at t=0 (Eq. 43); exact long-time averages after saturation controlled by J_f and {gamma_q(t_f), dot_gamma_q(t_f)} (Tables II and III); and a characterization of the adiabatic limit. The analytical results are benchmarked against numerical solutions of the Ermakov equation for several protocols and temperatures.","tokens_in":55327,"tokens_out":14062,"duration_ms":151902,"significance":"Within the quadratic model, the paper is thorough and internally consistent: the appendices contain complete derivations, the symplectic-eigenvalue formulas are explicit, and Figs. 3-6 show good agreement between the approximate closed forms and numerical integration of the Ermakov equations. The reduction of all information-theoretic quantities to the single function F[gamma_q(t)] is elegant, and the early-time exponent and the L sqrt(4 omega J_f/v) scaling of the sudden-quench averages are concrete, falsifiable predictions of the Gaussian model. The central caveat, acknowledged in Sec. VII but not resolved there, is that the mapping from the physical sine-Gordon tunneling term to the quadratic Hamiltonian of Eq. (6) is not controlled in some regimes where predictions are made; if that mapping can be justified or tested, the framework would be a valuable benchmark for field-space entanglement in 1D quantum simulators.","major_comments":[{"comment":"The early-time power law is the headline physical prediction, but it is derived exactly in the regime where the harmonic approximation of the tunneling cosine is least controlled. Equation (6) is obtained by expanding cos(phi_A - phi_B) to quadratic order around a phase minimum; this is standard when the tunneling-induced mass is large relative to fluctuations. However, the scaling result (43) follows from Taylor-expanding J(t) around t=0 (Eq. (36)), i.e. for J(t) arbitrarily small, where the relative phase is essentially unpinned and the omitted non-Gaussian terms of the cosine are not small. Section VII concedes a possible breakdown in the adiabatic slow-ramp regime, but not in the early-time window. Since Eq. (43) is the central claim, the manuscript should either (i) provide a quantitative condition under which the quadratic expansion remains valid at early times, (ii) test the prediction against the full time-dependent sine-Gordon dynamics using, for example, the truncated-Wigner or self-consistent harmonic methods cited in Sec. VII, or (iii) explicitly limit the claim to the Gaussian model and avoid presenting the exponent as a property of tunnel-coupled Luttinger liquids.","section":"Sec. II A and Sec. IV, Eqs. (6) and (43)"}],"minor_comments":[{"comment":"The definition of M_q in Table III appears garbled and dimensionally inconsistent; it should read M_q = 1/2 [gamma_q^2(t_f) + (dot_gamma_q(t_f)/Omega_qf)^2 + (zeta_q/gamma_q(t_f))^2], consistent with Eq. (51), but the printed expression contains an extra factor of gamma_q^2(t_f).","section":"Table III"},{"comment":"The crossover time t_cut^(0) is written in a way that looks dimensionful; the statement should be made as the dimensionless ratio t_cut^(0)/t_0^* = v pi/(sqrt(omega J_f) L), as in the caption of Fig. 3.","section":"Sec. IV (after Eq. (43))"},{"comment":"The removal of the zero mode q=0 is justified only heuristically by saying it is one out of a continuum of modes; for the finite-size expressions such as Eq. (42), an estimate of the q=0 contribution would be needed before the finite-L quantitative claims can be taken as physical.","section":"Sec. II B"},{"comment":"The convergence statements 'tf ~ 10 L/v for MI' and 'tf ~ L/v for LN' are based on visual overlap with the adiabatic solution; a quantitative error measure or a bound on the deviation would make these claims more robust.","section":"Sec. VI"}],"recommendation":"major_revision","confidential_remarks":"The main risk is external validity: the Gaussian expansion of the cosine is uncontrolled precisely in the early-time window that produces the headline scaling law. I would encourage the editor to request either a comparison with sine-Gordon numerics or an explicit and prominent restriction of the claims to the quadratic model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely new: it gives a complete analytic solution for mutual information, logarithmic negativity, and Rényi entropies in a Gaussian model of two tunnel-coupled Luttinger liquids, and it organizes the whole dynamics around the Ermakov factor. The early-time scaling law — I grows like t^{n+1} with n fixed by the first non-null derivative of J(t) — is new, and so are the exact long-time averages in terms of J_f and the Ermakov factors at saturation. The derivations are careful and internally consistent, and the numerical checks against exact Ermakov solutions in Figs. 3–6 are the right kind of evidence. The threshold temperature for entanglement is a nice byproduct. The special-function results in Tables II and III are substantial work, and the appendices look complete.\n\nThe main soft spot is exactly what the stress-test note flags: the harmonic approximation of the cosine is load-bearing, and the early-time result is derived for J(t) near zero, where the mass is small and the relative phase is effectively unpinned. That is where the quadratic expansion is least controlled. The paper acknowledges the limitation for slow adiabatic ramps in Sec. VII, but not for the early-time regime that produces its headline exponent. So the universal power law may be a property of the Gaussian model rather than the physical sine-Gordon dynamics. I don't think this kills the paper: the abstract and Sec. II are explicit that the results are within a Gaussian approximation, so the internal math holds up. But the physical claim needs a caveat, and ideally a sanity check against the full time-dependent sine-Gordon model (truncated Wigner or self-consistent harmonic methods) for at least one smooth ramp.\n\nLesser points: the zero-mode omission in Sec. II B is heuristic but probably minor, shifting results by O(1) or O(log L). The abstract also slightly overstates analytic completeness — early-time LN is left open, and finite-ramp LN is an integral. Those are framing issues, not mathematical errors.\n\nBottom line: this is a serious piece of theory, clearly thought through, and the central Gaussian machinery holds together. It deserves a full referee pass, with the referee asked to push on the early-time validity condition before publication.","headline":"A technically careful Gaussian solution with new early-time and long-time results, but the headline scaling law sits exactly where the harmonic approximation is least controlled; worth a serious referee who pushes on the regime of validity.","tokens_in":55890,"tokens_out":2128,"would_cite":true,"duration_ms":24562,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Within a Gaussian approximation, field-space entanglement between two tunnel-coupled Luttinger liquids is fully encoded in the Ermakov factors $\\gamma_q(t)$, so mutual information, logarithmic negativity, and Rényi entropies are known…","keywords":["field-space entanglement","Luttinger liquid","Ermakov equation","logarithmic negativity","mutual information","Rényi entropies","time-dependent tunneling","Gaussian approximation"],"falsifier":"Solve the full time-dependent sine-Gordon model for a sudden quench at zero temperature and extract the early-time mutual information in the thermodynamic limit; the Gaussian theory predicts $I(A{:}B)\\propto t$, so a robustly different power law would falsify the claim that the Ermakov factors control the entanglement dynamics.","tokens_in":54897,"feed_emoji":"⚛️","tokens_out":14994,"duration_ms":135599,"temperature":0.7,"pith_summary":"Two gapless one-dimensional quantum fluids, initially independent and at zero or finite temperature, are brought into contact by a tunneling term acting uniformly along their whole length; the paper asks how entanglement between the fields themselves grows in time. The answer it establishes is that, within a Gaussian approximation to the tunneling interaction, every information-theoretic measure of field-space entanglement — mutual information, logarithmic negativity, and Rényi entropies — is a known function of the mode-dependent Ermakov factors $\\gamma_q(t)$, one nonlinear ordinary differential equation per momentum mode. From that single ingredient the paper derives a zero-temperature early-time power law $I(A{:}B)\\propto t^{n+1}$ in the thermodynamic limit, with $n$ the order of the first non-null derivative of the tunneling protocol at $t=0$, and exact long-time averages after the coupling saturates that depend only on the final coupling and the terminal values of $\\gamma_q$ and its derivative. The reason to care is that the construction converts a nonequilibrium many-body question into one forced harmonic oscillator per mode, and all predictions are expressed in quantities that 1D Bose-gas simulators can in principle measure.","feed_headline":"Coupled 1D fluids: entanglement grows by a t^(n+1) law","feed_subtitle":"One nonlinear equation drives mutual information, negativity, and Rényi entropies for any tunneling ramp.","key_machinery":"The central object is the set of Ermakov factors $\\gamma_q(t)$: positive mode-dependent scale factors solving $\\ddot{\\gamma}_q+\\Omega_q(t)^2\\gamma_q=\\Omega_{q0}^2/\\gamma_q^3$, with $\\gamma_q(0)=1$ and $\\dot{\\gamma}_q(0)=0$, where $\\Omega_q(t)^2=(vq)^2+4\\omega J(t)$ is the instantaneous frequency of the antisymmetric sector. The two coupled fluids are first exchanged for symmetric and antisymmetric phase combinations; the symmetric sector stays a gapless Luttinger liquid, while the antisymmetric sector becomes a time-dependent gapped oscillator at each momentum. The mechanism that generates field-space entanglement is the difference between those two sectors, and the Ermakov factor exactly encodes it: it factorizes the antisymmetric symplectic propagator into a squeezing-and-shearing part built from $(\\gamma_q,\\dot\\gamma_q)$ and a rotation part, and every entropic quantity depends only on $F[\\gamma_q]$. Solving one ordinary differential equation per mode therefore replaces the full field-theoretic dynamics.","core_discovery":"The central claim is that field-space entanglement dynamics between two identical tunnel-coupled Luttinger liquids is solvable at the Gaussian level: for each positive momentum mode $q$, everything follows from the Ermakov equation $\\ddot{\\gamma}_q(t)+\\Omega_q(t)^2\\gamma_q(t)=\\Omega_{q0}^2/\\gamma_q^3(t)$ with $\\gamma_q(0)=1$, $\\dot{\\gamma}_q(0)=0$, and $\\Omega_q(t)^2=(vq)^2+4\\omega J(t)$. Once the Ermakov factors are known, the symplectic eigenvalues that feed all correlation measures are fixed by the composite function $F[\\gamma_q(t)]=\\frac{1}{4}[(\\dot{\\gamma}_q/\\Omega_{q0})^2+(\\gamma_q-1/\\gamma_q)^2]$, so mutual information, logarithmic negativity, and Rényi entropies are known at every time and temperature for any tunneling profile. The paper reads three generic behaviours off these expressions: zero-temperature early-time MI grows as $t^{n+1}$ in the thermodynamic limit (after a finite-size transient $t^{2n+2}$ with logarithmic corrections), long-time averages after the ramp saturates are completely fixed by $J_f$ and $\\{\\gamma_q(t_f),\\dot\\gamma_q(t_f)\\}$, and finite-temperature logarithmic negativity is monotonically suppressed and vanishes above a protocol-dependent threshold. It also gives the adiabatic limit $\\gamma_q^{\\rm ad}(t)=\\sqrt{\\Omega_{q0}/\\Omega_q(t)}$ for slow ramps. These statements are made within the harmonic expansion of the tunneling cosine, which the paper treats as the defining regime of its analysis.","pith_inferences":["Not claimed by the paper: the power-law exponent $n+1$ offers a direct experimental fingerprint, since protocols whose first non-null derivative has order 0, 1, 2, ... should produce measured early-time MI exponents 1, 2, 3, ... in a 1D Bose-gas simulator.","An editorial stress test of the long-time reduction would compare two ramps with identical $J_f$, $\\gamma_q(t_f)$, and $\\dot\\gamma_q(t_f)$ but different interiors; equal time-averaged correlations would confirm that the terminal Ermakov data are the only memory of the protocol.","The paper drops the zero mode $q=0$; at short system sizes or under near-uniform initial conditions this single mode could contribute visibly, so a calculation retaining the zero mode would bound the size of that correction.","The slow-ramp regime is where the harmonic approximation is least guaranteed; analyzing the same field-space entanglement measures in the full sine-Gordon dynamics would reveal whether the adiabatic predictions survive beyond the Gaussian regime."],"forward_implications":["At zero temperature and in the thermodynamic limit, early-time mutual information is a pure power law $I(A{:}B)\\propto t^{n+1}$: linear for a sudden quench, quadratic for a linear ramp, cubic for a quadratic ramp, and so on.","After the coupling has saturated, the long-time averages of mutual information, logarithmic negativity, and Rényi entropies carry no memory of the ramp interior: only $J_f$ and the terminal pair $\\{\\gamma_q(t_f),\\dot\\gamma_q(t_f)\\}$ enter.","For a zero-temperature sudden quench, the thermodynamic-limit averaged mutual information, logarithmic negativity, and Rényi-2 entropy all scale as $L\\sqrt{4\\omega J_f}/v$.","At finite temperature the averaged logarithmic negativity decreases monotonically and drops to zero at a threshold $T^*$; the averaged mutual information also decreases with temperature but remains positive.","For ramp durations much longer than the slowest-mode period $L/v$, mutual information and logarithmic negativity converge to the adiabatic evolution of the same Gaussian equations."],"supporting_citations":[{"why":"Motivates and grounds the field-space partition between two tunnel-coupled Tomonaga-Luttinger liquids and its experimental detection.","marker":"[44]"},{"why":"Supplies the nonlinear differential-equation solution (Ermakov/Pinney) on which the mode-by-mode scale factors rest.","marker":"[53]"},{"why":"Provides the Gaussian-state framework: covariance matrices, symplectic eigenvalues, entropies, and logarithmic negativity.","marker":"[54]"},{"why":"Gives the exact sudden-quench Ermakov solution used as the basis of the early-time and long-time quench results.","marker":"[66]"},{"why":"Gives the analytic linear-ramp Ermakov solution used to verify the finite-ramp long-time averages.","marker":"[67]"},{"why":"Defines the Luttinger-liquid low-energy description that underlies the phase and density field decomposition.","marker":"[46]"},{"why":"Justifies the sine-Gordon form of the tunneling coupling between the two one-dimensional condensates.","marker":"[41]"}],"fun_headline_variants":["Field-space entanglement: Luttinger liquids obey t^(n+1) growth","Ermakov equation controls entanglement between tunnel-coupled 1D fluids","Tunneling ramp sets entanglement scaling: t^(n+1) for Luttinger liquids","Exact entanglement dynamics of two coupled Luttinger liquids from one ODE","Luttinger liquids: field-space entanglement follows a t^(n+1) law exactly"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the cosine of the relative phase in the tunneling coupling may be replaced by its quadratic expansion around a phase minimum; if this harmonic/Gaussian approximation fails during the ramp, the Ermakov-factor equations and all derived power laws and long-time averages no longer describe the physical model.","fun_headline_variants_meta":{"raw":{"variants":["Field-space entanglement: Luttinger liquids obey t^(n+1) growth","Ermakov equation controls entanglement between tunnel-coupled 1D fluids","Tunneling ramp sets entanglement scaling: t^(n+1) for Luttinger liquids","Exact entanglement dynamics of two coupled Luttinger liquids from one ODE","Luttinger liquids: field-space entanglement follows a t^(n+1) law exactly"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000639,"raw_usage":{"total_tokens":3045,"prompt_tokens":1151,"completion_tokens":1894,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":767,"completion_tokens_details":{"reasoning_tokens":1787}},"tokens_in":767,"tokens_out":1894,"duration_ms":14744,"temperature":1.0,"reasoning_tokens":1787,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T20:10:25.432324+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the full time-dependent sine-Gordon model for a sudden quench at zero temperature and extract the early-time mutual information in the thermodynamic limit; the Gaussian theory predicts $I(A{:}B)\\propto t$, so a robustly different power law would falsify the claim that the Ermakov factors control the entanglement dynamics.","supporting_citations":[{"cited_title":"Murtadho, M","cited_arxiv_id":null,"evidence_quote":"Motivates and grounds the field-space partition between two tunnel-coupled Tomonaga-Luttinger liquids and its experimental detection."},{"cited_title":"Pinney, The nonlinear differential equation y′′ + p(x) y + c y−3 = 0, Proc","cited_arxiv_id":null,"evidence_quote":"Supplies the nonlinear differential-equation solution (Ermakov/Pinney) on which the mode-by-mode scale factors rest."},{"cited_title":"Serafini, Quantum continuous variables: A primer of theoretical meth ods (1st ed.), CRC Press (2017)","cited_arxiv_id":null,"evidence_quote":"Provides the Gaussian-state framework: covariance matrices, symplectic eigenvalues, entropies, and logarithmic negativity."},{"cited_title":"Dupays, D","cited_arxiv_id":null,"evidence_quote":"Gives the exact sudden-quench Ermakov solution used as the basis of the early-time and long-time quench results."},{"cited_title":"Dupays, B","cited_arxiv_id":null,"evidence_quote":"Gives the analytic linear-ramp Ermakov solution used to verify the finite-ramp long-time averages."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Luttinger-liquid low-energy description that underlies the phase and density field decomposition."},{"cited_title":"Gritsev, A","cited_arxiv_id":null,"evidence_quote":"Justifies the sine-Gordon form of the tunneling coupling between the two one-dimensional condensates."}],"review_version":1}