{"id":"b0f31924-1f9a-4aa9-b7ca-74089d6ca120","arxiv_id":"2507.16239","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"SDFT with a Kirkwood-factor renormalization quantitatively reproduces longitudinal and transverse polarization dynamics in the Stockmayer fluid, as verified against Brownian dynamics simulations.","lead":"This paper develops a stochastic density functional theory (SDFT) for the dynamics of polarized molecules in a dipolar fluid, and compares its predictions against Brownian dynamics simulations of the Stockmayer fluid. The theory captures longitudinal polarization relaxation well, and can be patched with the Kirkwood factor to also match the transverse fluctuations, pointing the way to cheap coarse-grained models for electrolyte and solvent dynamics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The transverse normalization is internally inconsistent (Eqs. 13, 35, 46 vs. Eq. 9 with nu_T=2), leaving the claimed quantitative match of S_T(q,omega) dependent on an unresolved factor-of-2 convention.","rationale":"The reader identified the Kivelson-Madden relation as the weakest assumption, but the paper actually provides a direct test of the g_K rescaling: the effective-SDFT relaxation times, constructed with tau_rs_tilde=tau_rs g_K, reproduce the simulated longitudinal and transverse relaxation times in Fig. 1(c)-(d), and the frequency-dependent spectra in Figs. 4-5 then validate the Lorentzian form. Thus that concern is substantially mitigated by the paper's own comparisons. The unresolved factor-of-2 normalization, by contrast, is a purely internal inconsistency in the published equations: Eq. (13), Eq. (35), and Eq. (46) cannot all be correct with the definitions in Eqs. (3) and (9). Because the transverse PSD amplitude is the central observable, a reader cannot reproduce Figs. 4-5 unambiguously from the text. This does not invalidate the physical idea of Kirkwood-factor renormalization, but it means the quantitative claim is conditional on a convention that must be stated and checked. The recommended verdict therefore remains CONDITIONAL, unchanged from the reader, because the work is otherwise solid and the issue is likely correctable; the test above would settle whether the plotted agreement is real or an artifact of normalization.","tokens_in":19193,"tokens_out":13813,"duration_ms":147258,"concrete_test":"Recompute S_T from the stored BD trajectories for the largest-dipole system (p*=6.15) using both conventions: the full transverse vector P_T=P-(P.qhat)qhat (two components) and a single Cartesian transverse component. Then evaluate Eq. (37) with the effective parameters from Eqs. (47)-(51) and the discrete-sampling correction Eq. (A1) at qmin. If the full-vector convention is used, Eq. (37) should match the simulated PSD; if only a single component is used, Eq. (37) overestimates by a factor of 2 and the plotted agreement must be renormalized. This single check decides whether the quantitative comparison in Fig. 4 is internally consistent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in Section V.E rests on the absolute amplitude of the transverse dynamic structure factor, but the manuscript uses two incompatible transverse normalizations. Equation (3) defines F_T with the full transverse vector, and Eq. (35) then gives F_T(0)=2p^2/3 for independent dipoles; Eq. (9) consistently has nu_T=2. Together with chi_L=1-1/epsilon_r and chi_T=epsilon_r-1, these relations imply S_L/S_T=1/(2 epsilon_r) and S_T(q->0)=2 p^2 g_K/3. In contrast, Section V.A and Eq. (46) state S_T(q->0)=p^2/3, and Eq. (13) gives S_L/S_T=1/epsilon_r, missing the transverse degeneracy factor of 2. The effective-dipole rescaling Eq. (47), p_tilde=p sqrt(g_K), is calibrated against the latter normalization; if the full-vector convention is used, the required rescaling for the transverse amplitude would instead be p_tilde=p sqrt(g_K/2). Since Eq. (37) contains an explicit factor 2 in the static amplitude, it is unclear whether the S_T(q,omega) curves in Figs. 4 and 5 were produced with the full-vector or per-component convention. This is an internal inconsistency, not a matter of external consensus, and it directly determines the headline quantity. The Kivelson-Madden assumption is a weaker concern because the relaxation-time comparisons in Fig. 1(d) and the spectra in Fig. 4 effectively test the g_K rescaling; the factor-of-2 ambiguity is not resolved anywhere in the text.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a linearized stochastic density functional theory (SDFT) for the polarization field of a Stockmayer fluid and obtains explicit expressions for the longitudinal and transverse intermediate scattering functions and dynamic structure factors. The predictions are compared with Brownian dynamics simulations for a range of dipole moments and for variations of the translational and rotational diffusion coefficients. The authors find that plain SDFT describes longitudinal polarization fluctuations well but underestimates transverse fluctuations, and they propose an 'effective SDFT' in which the dipole moment is rescaled by the square root of the Kirkwood factor and the diffusion coefficients are divided by that factor. They report that this effective theory reproduces the simulated dynamic structure factors, including their dependence on the diffusion coefficients.","tokens_in":19486,"tokens_out":13602,"duration_ms":147867,"significance":"If established, the result would provide a simple analytic coarse-grained description of polarization dynamics in dipolar fluids, with potential applications to electrolytes and confined polar solvents. The paper has clear strengths: the SDFT equations are derived from a microscopic Langevin description, the simulation campaign covers a meaningful parameter range, and the dynamic predictions of the effective theory contain no additional free parameters beyond the statically measured Kirkwood factor. The dynamic lineshape comparisons in Figs. 4 and 5 are a nontrivial test of the theory. However, as written, the manuscript contains an internal normalization inconsistency in the transverse structure factor and an unexplained amplitude mismatch at finite wavevector; these issues directly affect the headline quantitative claims.","major_comments":[{"comment":"The transverse polarization structure factor is defined with two different normalizations. In Eq. (3), P_T is a two-component vector and the dot product is used, which gives the full transverse variance; this is consistent with Eq. (35), which yields F_T(0)=2p^2/3, and with the prefactor nu_T=2 in Eq. (9). Section V.A and Eq. (46), however, state S_T(q->0)=p^2/3 and S_L/S_T=1/epsilon_r, which are the per-component normalizations. The two conventions differ by a factor of 2. This matters because Eq. (47), tilde p^2=p^2 g_K, is calibrated against Eq. (46), i.e. against S_T(q->0)=p^2 g_K/3, whereas Eq. (37) contains an explicit factor 2 in the transverse dynamic structure factor amplitude. If the full-vector convention is used, tilde p=p sqrt(g_K) correctly gives S_T=2p^2 g_K/3; if the per-component convention is used, Eq. (37) must lose its factor 2, or the rescaling must instead be tilde p=p sqrt(g_K/2). The plotted values in Fig. 1(b) appear to follow the per-component convention, but the text says the curves are computed from Eqs. (36)-(37). The authors must choose one convention and apply it consistently to the equations, the simulations, and the figures; otherwise the claimed quantitative match of S_T(q,omega) in Figs. 4 and 5 is not defined unambiguously.","section":"Eqs. (13), (35), (37), (46), (47); Section V.A"},{"comment":"Effective SDFT as used in Fig. 5 is stated to be Eqs. (36)-(37) with effective parameters, whose static amplitudes are q-independent. Yet Section V.D reports that for the largest dipole, p*=6.15, the transverse static structure factor follows S_T(q)=p^2 g_K/[3(1+q^2 xi_T^2)] with xi_T~1.28 sigma. At q sigma=0.5 this gives S_T(q) about 29% smaller than S_T(q->0). Since the low-frequency plateau of S_T(q,omega) is S_T(q), the effective SDFT lines plotted with the q->0 amplitude cannot quantitatively agree with the BD data unless an additional q-dependent prefactor is introduced into Eq. (37). The text acknowledges that S_T(q=0.5 sigma^{-1}) differs significantly from S_T(q->0) but does not explain how effective SDFT accounts for this difference. The authors should either include the q-dependent static amplitude in the effective SDFT prediction, or restrict the claim of quantitative agreement to the q->0 limit and present the finite-q comparison as showing only the lineshape/frequency dependence.","section":"Section V.D, Section V.E, Fig. 5"},{"comment":"The dynamic rescaling tilde D_r^s = D_r^s/g_K is justified by the Kivelson-Madden relation tau_D = tau_1 g_K, imported from Refs. 53, 72, 73. This relation is not directly tested for the present Stockmayer fluid, and it is a central ingredient of the effective theory: without it, the frequency scale of the transverse spectra would be wrong. The agreement in Fig. 4 provides indirect evidence, but a more direct validation would strengthen the paper: the authors can compute the single-dipole reorientation time tau_1 from the same BD trajectories and check the relation tau_T(q->0) = tau_1 g_K. If that check is not feasible, the text should state explicitly that Eq. (50) is assumed rather than measured, so that the reader can weigh the evidence accordingly.","section":"Section V.C, Eq. (50)"}],"minor_comments":[{"comment":"In the first sentence of Section V.A, 'the longitudinal and static structure factors' should read 'the longitudinal and transverse structure factors'.","section":"Section V.A"},{"comment":"The word 'meau sed' in the caption is a typo and should be 'measured'.","section":"Fig. 1 caption"},{"comment":"The line tilde C_s = C_s = C_s g_K^0 is notationally awkward; since the exponent 0 is trivial, the statement that the concentration is unchanged should be made directly.","section":"Eqs. (47)-(49)"},{"comment":"Reference 50 is cited as a 2024 preprint without a journal reference; if it has been published or accepted, the full citation should be provided.","section":"Reference [50]"},{"comment":"The data availability statement is an incomplete placeholder; it should be completed before publication if the journal requires it.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":"The factor-of-2 normalization issue is genuine and currently affects the interpretation of the central plots. The finite-q amplitude concern is also specific and testable. Both are fixable with careful revision, and the underlying approach remains promising, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely new: it takes the polarization SDFT from Ref. 50, derives the longitudinal and transverse intermediate scattering functions and dynamic structure factors for the Stockmayer fluid, and confronts them with Brownian Dynamics simulations. The failure of bare SDFT on transverse fluctuations and the fix via a Kirkwood-factor rescaling are clearly identified. The effective SDFT matches the simulated dynamic structure factors over a wide frequency range, across dipole strengths, and for varied D_s and D_r. That is a practical and useful result for anyone wanting a tractable coarse-grained description of polar solvent dynamics.\n\nThe main soft spot is a real internal inconsistency in the transverse normalization. Eq. (3) defines the transverse ISF with the full two-component transverse vector, so F_T(0)=2p^2/3 in Eq. (35), and Eq. (9) consistently uses nu_T=2. But Section V.A and Eq. (46) quote S_T(q->0)=p^2 g_K/3, and Eq. (13) says S_L/S_T=1/eps_r, which only holds for a single transverse component. The effective dipole rescaling p_tilde=p sqrt(g_K) is calibrated against the per-component version, yet Eq. (37) contains an explicit factor 2 in the static amplitude. As written, the theory curve for S_T is twice as large as the BD data in Fig. 1(b) if that figure is per-component, which it appears to be. The authors do not state which convention they used for the comparison in Figs. 4 and 5. This is not a matter of taste; it affects the headline quantitative claim. It is fixable by carefully defining S_T as either the full transverse vector or a single component, then adjusting Eqs. (9), (13), (35), (37), and (46) consistently.\n\nTwo other points, both secondary. First, the rescaling of D_r through the Kivelson-Madden relation tau_D = tau_1 g_K is imported rather than independently verified, though the dynamic spectra comparisons in Figs. 4 and 5 do exercise that assumption and it survives. Second, the static agreement is partly circular because g_K is measured from the same simulations; the dynamic predictions are not fitted, so the central claim does not collapse. The data availability section is still a placeholder, and for a comparison paper that matters.\n\nSerious referee time: yes. The science is sound enough that a referee's work will improve the paper rather than bury it. Send it out, ask specifically about the factor of 2, and require the authors to state the normalization explicitly. After that is resolved, this is a solid contribution worth publishing.","headline":"Useful SDFT-to-simulation confrontation with a workable Kirkwood rescaling, but the transverse structure-factor normalization is internally inconsistent and needs cleaning before the quantitative claim holds.","tokens_in":20081,"tokens_out":4901,"would_cite":true,"duration_ms":48589,"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":"Effective stochastic density functional theory, with its four parameters rescaled by the Kirkwood factor, quantitatively reproduces the dynamic polarization structure factors of the Stockmayer fluid across dipole strengths and frequencies.","keywords":["Stockmayer fluid","stochastic density functional theory","polarization dynamics","intermediate scattering function","dynamic structure factor","Kirkwood factor","dielectric relaxation","Brownian dynamics simulation"],"falsifier":"In the same Brownian dynamics simulations, compute the single-dipole reorientation time $\\tau_1=\\int_0^\\infty \\langle\\hat{\\mathbf{u}}_i(t)\\cdot\\hat{\\mathbf{u}}_i(0)\\rangle\\,dt$ at each dipole strength and test whether $\\tau_T(q_{\\min})$ equals $g_K\\tau_1$; a departure of $\\tau_T/(g_K\\tau_1)$ from unity beyond statistical error across the $p^*$ range would falsify the diffusion-coefficient rescaling. A second check is to measure $S_T(q,\\omega)$ at wavevectors with $q\\xi_T>1$ for the strongest dipoles, where the $q$-independent static structure factor underlying the effective SDFT should break down.","tokens_in":18933,"feed_emoji":"⚛","tokens_out":13178,"duration_ms":126937,"temperature":0.7,"pith_summary":"This paper asks whether stochastic density functional theory (SDFT), a fluctuating mean-field description built from overdamped Langevin dynamics, can describe the time-dependent polarization fluctuations of a dipolar liquid. Using the Stockmayer fluid—Lennard-Jones particles carrying point dipoles—as a microscopic benchmark, the authors derive closed-form expressions for the intermediate scattering functions and dynamic structure factors of the longitudinal and transverse polarization components. They find that the bare theory captures the longitudinal channel but underestimates the transverse channel, and that both defects trace to neglected short-range orientational correlations. Rescaling the dipole moment and the two diffusion coefficients by the Kirkwood factor, a static measure of orientational order, brings the theory into quantitative agreement with Brownian dynamics simulations for both channels over the studied frequency range. If the result holds more generally, it gives a cheap, analytical route to dielectric relaxation spectra of polar fluids and their role in electrolytes.","feed_headline":"One rescaling fixes a mean-field theory of dipolar fluid dynamics","feed_subtitle":"Folding the Kirkwood factor into mean-field theory reproduces full polarization spectra.","key_machinery":"The central object is the Fourier-decomposed polarization density $\\mathbf{P}(\\mathbf{r},t)=\\sum_i p\\,\\hat{\\mathbf{u}}_i(t)\\,\\delta(\\mathbf{r}-\\mathbf{r}_i(t))$, split into longitudinal and transverse components with respect to the wavevector $\\mathbf{q}$. Linearized SDFT gives each component an Ornstein-Uhlenbeck relaxation with a single time scale, so the theory's output is two Lorentzian dynamic structure factors $S_{L,T}(q,\\omega)$. The load-bearing correction is the Kirkwood factor $g_K=\\langle M^2\\rangle/(N p^2)$, a static measure of local orientational correlations: it rescales the effective dipole as $p\\sqrt{g_K}$ and the rotational and translational diffusion coefficients by $1/g_K$, injecting the missing dipole correlations without changing the analytical form of the theory.","core_discovery":"Starting from the overdamped Langevin equations for positions and orientations, linearized SDFT yields a polarization field whose longitudinal and transverse Fourier components each relax as a single exponential, with rates $1/\\tau_L(q)=(1+3y)(1+q^2 a^2)/\\tau_s^r$ and $1/\\tau_T(q)=(1+q^2 a^2)/\\tau_s^r$; the intermediate scattering functions and dynamic structure factors are therefore exponential and Lorentzian. Against Brownian dynamics simulations of the Stockmayer fluid, the bare theory matches the longitudinal static structure factor and relaxation time but underestimates the transverse ones once dipole correlations develop. Replacing the bare parameters by effective ones, $\\tilde{p}=p\\,g_K^{1/2}$, $\\tilde{y}=y\\,g_K$, $\\tilde{D}_r=D_r/g_K$ and $\\tilde{D}_s=D_s/g_K$, with the Kirkwood factor $g_K=\\langle M^2\\rangle/(N p^2)$ measured from the same simulations, reproduces the static structure factors, relaxation times, and full frequency-dependent dynamic structure factors of both components for all dipole strengths studied, including at finite wavevector $q\\sigma=0.5$.","pith_inferences":["A direct measurement of $\\tau_1$ in the same simulations would convert the imported Debye-time relation into a testable input; if the relation is only approximate, the scheme could keep its structure while replacing $g_K$ with the simulated ratio $\\tau_T/\\tau_1$.","Because the effective parameters are fixed by a static average, one could compute $g_K$ from a short equilibrium run or from a static theory and then predict the full polarization spectrum, making SDFT a one-input dynamical theory.","At the strongest dipoles the transverse structure factor acquires a wavevector dependence with a correlation length of about one particle diameter, so effective SDFT should fail once $q\\xi_T$ approaches unity; locating this crossover in simulations would delimit the theory's range.","The same Kirkwood-factor rescaling could be imported into the SDFT description of ions in a polar solvent, replacing bare solvent parameters before computing electrolyte response, which the paper identifies as a future direction but does not carry out."],"forward_implications":["Effective SDFT predicts that polarization intermediate scattering functions remain single exponentials, so the dynamic structure factors keep their Lorentzian shape over the whole studied frequency window.","Agreement in the longitudinal channel is not a stringent test by itself: errors cancel between the numerator and denominator of $S_L(q\\to 0)$, while the transverse channel exposes the missing correlations.","The four-parameter rescaling ties the full frequency-dependent spectrum to one static equilibrium number, $g_K$, measured from the same simulations, so no dynamical fitting is required.","In the hydrodynamic limit the rescaled theory gives a Debye permittivity with relaxation time $\\tau_D=g_K\\tau_s^r$, connecting the microscopic rescaling to dielectric relaxation measurements.","The same effective parameters also describe the dynamics at finite wavevector $q\\sigma=0.5$ even when $g_K\\approx 12.6$, indicating the mapping extends beyond the $q\\to 0$ limit."],"supporting_citations":[{"why":"Foundational Langevin equation for the density of interacting Brownian particles, from which the SDFT evolution equations derive.","marker":"[31]"},{"why":"Earlier SDFT treatment of a polar solvent that supplies the linearized polarization evolution equation and the analytical framework extended to the Stockmayer fluid.","marker":"[50]"},{"why":"Source of the Kirkwood g-factor and the approximate relation between the Debye relaxation time and single-dipole reorientation time used to rescale the rotational diffusion coefficient.","marker":"[53,72,73]"},{"why":"Boundary-condition-dependent relations between permittivity and dipole fluctuations, used to connect the Kirkwood factor to the static permittivity via tin-foil boundary conditions.","marker":"[59,60]"},{"why":"Molecular density functional theory study of the Stockmayer fluid whose parameters set the reference simulation system and its expected permittivity.","marker":"[24]"},{"why":"Provides the Lennard-Jones interaction parameters used in the Brownian dynamics simulations.","marker":"[63]"},{"why":"P3M algorithm used to compute long-range dipolar interactions in the simulations that produce the reference data.","marker":"[64]"}],"fun_headline_variants":["Kirkwood rescaling fixes transverse polar dynamics in SDFT","One factor extends mean-field dipole theory to simulations","Rescaling with Kirkwood factor matches full dipole spectra","SDFT plus Kirkwood factor nails polar fluid relaxation","Simple rescaling makes mean-field theory exact for dipoles"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The dynamical rescaling rests on an imported approximate relation, $\\tau_D=\\tau_1 g_K$, stating that the collective Debye relaxation time equals the single-dipole reorientation time multiplied by the Kirkwood factor; the paper does not independently measure $\\tau_1$ in the Stockmayer fluid it simulates, so if the relation is inaccurate for this system the renormalized dynamics would fail even though the static agreement would survive.","fun_headline_variants_meta":{"raw":{"variants":["Kirkwood rescaling fixes transverse polar dynamics in SDFT","One factor extends mean-field dipole theory to simulations","Rescaling with Kirkwood factor matches full dipole spectra","SDFT plus Kirkwood factor nails polar fluid relaxation","Simple rescaling makes mean-field theory exact for dipoles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000148,"raw_usage":{"total_tokens":1212,"prompt_tokens":987,"completion_tokens":225,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":603,"completion_tokens_details":{"reasoning_tokens":145}},"tokens_in":603,"tokens_out":225,"duration_ms":2759,"temperature":1.0,"reasoning_tokens":145,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:16:22.010443+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In the same Brownian dynamics simulations, compute the single-dipole reorientation time $\\tau_1=\\int_0^\\infty \\langle\\hat{\\mathbf{u}}_i(t)\\cdot\\hat{\\mathbf{u}}_i(0)\\rangle\\,dt$ at each dipole strength and test whether $\\tau_T(q_{\\min})$ equals $g_K\\tau_1$; a departure of $\\tau_T/(g_K\\tau_1)$ from unity beyond statistical error across the $p^*$ range would falsify the diffusion-coefficient rescaling. A second check is to measure $S_T(q,\\omega)$ at wavevectors with $q\\xi_T>1$ for the strongest dipoles, where the $q$-independent static structure factor underlying the effective SDFT should break down.","supporting_citations":[{"cited_title":"Stochastic density functional theory for ions in a polar solvent","cited_arxiv_id":"2407.17232","evidence_quote":"Earlier SDFT treatment of a polar solvent that supplies the linearized polarization evolution equation and the analytical framework extended to the Stockmayer fluid."},{"cited_title":"Ramirez , author R","cited_arxiv_id":null,"evidence_quote":"Provides the Lennard-Jones interaction parameters used in the Brownian dynamics simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"P3M algorithm used to compute long-range dipolar interactions in the simulations that produce the reference data."}],"review_version":1}