{"id":"24100d16-9d7e-477e-bc10-1336d4930785","arxiv_id":"2607.16180","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A dynamical phase-field model with an electronic-polarization field predicts spatially resolved refractive index and electro-optic response in BaTiO3 microstructures, including >4000 pm/V near domain walls.","lead":"A new computer model simulates how light bends and changes speed inside ferroelectric crystals by connecting the crystal's electrical domains to its optical response. The model predicts that electro-optic effects are dramatically stronger at domain walls, which could guide the design of faster photonic devices.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Domain-wall electro-optic enhancement rests on an untested local electronic-response assumption: electronic gradient terms are neglected and bulk-fitted couplings are applied at walls.","rationale":"The reader's weakest-assumption analysis identifies the same central vulnerability: the local, instantaneously equilibrated electronic polarization with neglected electronic gradient terms. I agree that this is the least secure load-bearing assumption, because it directly underlies the paper's most striking new predictions—nanoscale domain-wall electro-optic enhancements and transient switching response. The concern is not that the model is internally inconsistent; within its stated assumptions the derivations are coherent. Rather, the transferability of bulk-derived coupling constants and the neglect of electronic nonlocality at walls are unsupported by any independent spatially resolved calculation or measurement. The quantitative agreement with film-averaged experiments is suggestive but not decisive: the PsiQuantum point is a single unscaled comparison, and the Eltes temperature-dependence comparison requires a free scaling factor. A DFT calculation of the wall-region electronic dielectric response would directly test the local-closure assumption and would either support or undermine the central claim. This does not change the reader's CONDITIONAL verdict: the paper is a credible methodological advance, but the strongest quantitative and microscopic claims should be regarded as conditional on this test and on making the scaling argument explicit.","tokens_in":15440,"tokens_out":8546,"duration_ms":87299,"concrete_test":"Perform a first-principles (DFT or constrained-polarization) calculation of the electronic dielectric stiffness B_e for a BaTiO3 supercell containing a 90° domain wall with the same P_L(x) profile as in Fig. 4, comparing B_e at the wall center and 2 nm away with Eq. (19) evaluated from the local P_L and stress. If the DFT wall-region B_e differs from the local-model value by more than ~20%, recompute the wall r in Fig. 4d-f; if the >4000 pm/V peak does not survive, the neglect of electronic gradient terms and the transfer of bulk-fitted parameters to the wall are load-bearing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing part of the central claim is the spatially resolved electro-optic response, specifically the predicted >4000 pm/V enhancements near domain walls and the 6000 pm/V transient during switching. These values are computed from Eq. (15) using B_e(x) = ε0 ∂²f/∂P_e², where the electronic free-energy density is strictly local: it contains a quadratic term in P_e plus coupling to the local lattice polarization and stress (Eqs. 21, 22, 24), and the electronic gradient energy is explicitly neglected (text after Eq. 1 and Methods Eq. 16). At a domain wall, P_L varies by roughly 0.5 C/m² over a few nanometres; there is no evidence that the bulk-fitted g_ijkl^ee and B_e0 coefficients remain valid in that region, nor that nonlocal electronic response does not smear or suppress the wall peak. Because r is obtained from a difference of inverse B_e tensors, even a moderate error in the wall-region B_e can translate into a large error in the claimed enhancement. The paper provides no spatially resolved atomistic or experimental check of the wall coefficients; the only external validation is film-averaged, and the Eltes comparison in Fig. 6d uses an arbitrary constant scaling factor. Thus the microscopically novel prediction is the least secure part of the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a dynamical phase-field model for ferroelectrics that adds an electronic polarization field to the usual lattice polarization order parameter. The electronic polarization is assumed to equilibrate instantaneously to the local lattice polarization and stress, and its inverse dielectric stiffness B_e(x) is computed from a local free-energy curvature (Eqs. 12, 19, 21-22). This yields spatially resolved refractive-index and electro-optic maps that follow the evolving ferroelectric domain structure. Applying the model to BaTiO3 thin films, the authors report local electro-optic coefficients exceeding 4000 pm/V near domain walls, transient film-averaged responses up to ~6000 pm/V during switching, and temperature-dependent effective coefficients that are compared with experiments on BaTiO3-on-Si films (Eltes et al., PsiQuantum). The central claim is that the ferroelectric domain structure strongly modifies the local electro-optic response and that the simulations quantitatively reproduce film measurements.","tokens_in":15812,"tokens_out":5571,"duration_ms":54421,"significance":"If the predictions are reliable, this would be a useful mesoscale tool: it connects ferroelectric microstructure to optical properties, which existing phase-field models do not do, and it provides a mechanism for domain-wall and phase-boundary contributions to the electro-optic response. The time-scale separation between electronic and lattice polarization is physically sensible, the perturbation solution in Eq. (13)-(14) is straightforward, and the simulated microstructures in Figs. 3-5 are realistic. The paper also makes data available. However, the headline quantitative claims rest on parameters imported from bulk BaTiO3 and on a strictly local electronic-response assumption at domain walls, so the degree of independent prediction is smaller than the abstract suggests.","major_comments":[{"comment":"The central claim of local electro-optic enhancement >4000 pm/V (and ~6000 pm/V during switching) is computed from B_e(x)=ε0 ∂²f/∂P_e², with the electronic free energy treated as strictly local and electronic gradient terms explicitly neglected. The coefficients B_e0 and g_ijkl^ee are bulk-fitted values. At a domain wall P_L changes by roughly 0.5 C/m² over a few nanometres, and there is no evidence that the bulk coefficients transfer to that region or that nonlocal electronic response does not smear or suppress the wall peak. Since r is obtained from a difference of inverse B_e tensors (Eq. 15), even a moderate error in the wall-region B_e can dominate the claimed enhancement. Please provide a concrete test—e.g., a first-principles calculation of the electronic dielectric response at a 90° or 180° wall—or explicitly reframe the >4000 pm/V result as an untested model prediction.","section":"Methods, Eqs. (16), (21)-(22); Results, Fig. 4"},{"comment":"The model imports B_e0 and g_ijkl^ee from bulk BaTiO3, so the simulated film-averaged electro-optic coefficients inherit the bulk baseline they are compared with. The agreement with the PsiQuantum measurement (1080 vs 988 pm/V at 295 K) is therefore not an independent prediction of the intrinsic coefficient. The genuinely new output is the microstructural/domain-wall contribution, but the paper does not separate it from the bulk-imposed baseline. Please decompose the simulated r into bulk-intrinsic and microstructure-induced parts, or state clearly that the baseline is reproduced by construction and only the microstructural modulation is predicted.","section":"Methods, Table S1; Results, Fig. 6c"},{"comment":"The claim of quantitative agreement with Eltes et al. is obtained after multiplying the simulated data by a constant factor chosen to match the 0° measurement at 300 K. This scaling tests only the temperature dependence and relative orientation ratios, not the absolute magnitude. Please report the unscaled comparison and a quantitative metric of agreement. If interfacial dead layers or incomplete poling are invoked, include a parameterized physical model rather than an arbitrary rescaling factor.","section":"Results, Fig. 6d; Discussion"},{"comment":"The large transient values during switching (up to 6000 pm/V) are slopes of the refractive-index hysteresis loop at the coercive field and include domain nucleation, growth, and wall motion. These are not conventional linear electro-optic coefficients and depend on field history, sweep rate, and dynamic domain behavior. Comparing them with bulk r_51=1300 pm/V is misleading unless they are explicitly labeled as effective, history-dependent responses. Please report the small-field linear-regime values and separate intrinsic and extrinsic contributions.","section":"Results, Figs. 4c and 5c; Eq. (15)"}],"minor_comments":[{"comment":"'quanitative agreement' should be 'quantitative agreement'.","section":"Introduction"},{"comment":"The symbol B_ij is used for the inverse optical dielectric tensor, while B_e earlier denotes the electronic dielectric stiffness. This dual use of B is confusing and should be disambiguated.","section":"Eq. (15)"},{"comment":"The eigenstrain expression includes a contribution from electronic polarization via π_ijkl, but the text immediately after says the electronic contribution is ignored. Please reconcile these statements.","section":"Methods, Eq. (24)"},{"comment":"The abstract describes the model as predicting 'wavelength-dependent' optical properties, but simulations are shown only at 1550 nm. Either include a wavelength sweep or temper the claim.","section":"Abstract / Results"}],"recommendation":"major_revision","confidential_remarks":"The model is plausible and the microstructural framework is a genuine contribution, but the headline quantitative claims currently outrun the evidence: the absolute agreement is partly inherited from fitted bulk parameters, the Eltes comparison uses an arbitrary scale factor, and the domain-wall enhancement rests on an unvalidated locality assumption. These are fixable by recalibrating the claims and adding targeted validation, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news here is the added electronic-polarization field in a phase-field framework. That lets the model compute local refractive index and electro-optic response directly from the evolving ferroelectric microstructure, which neither the lattice-only phase-field work (ref. 27) nor the authors' earlier static thermodynamic theory did. The separation of timescales is physically sensible, the perturbation derivation is clean, and the temperature-dependent phase-competition picture is a nice payoff. I believe the central methodological claim: this is a new simulation route to local optical properties at device-relevant length scales.\n\nWhat the paper does well: it states the model clearly, gives the free-energy structure, and shows how the inverse electronic stiffness becomes the optical dielectric tensor. The BaTiO3 case study is appropriate, and the comparison with two experimental datasets gives the temperature dependence a real anchor. The DOI for figure data is a plus.\n\nWhere it gets soft, in proportion: the validation is partly circular. B_e0 and the g_ijkl^ee couplings are imported or fitted to reproduce bulk BaTiO3 n and r, so the film-averaged electro-optic coefficients inherit the bulk value they are compared with. The Eltes comparison in Fig. 6d requires an arbitrary constant scaling factor; calling that 'quantitative agreement' overstates what the data support. More importantly, the abstract's headline numbers—over 4000 pm/V near domain walls and a 6000 pm/V transient—rest on a strictly local electronic response with no electronic gradient term and no atomistic check of the bulk-fitted couplings inside a domain wall. That is a real weak spot, not a manufactured one. If the local-response assumption breaks down at the wall, those enhancements could be artifacts. The paper offers no spatially resolved experimental or first-principles check of that point.\n\nMinor: the code is proprietary (MuPRO), which limits reproducibility, though the data DOI helps. The paper also says 'quantitative' in a few places where the evidence is scaled or parameter-inherited.\n\nBottom line: the methodological advance is defensible and worth publishing, but the strongest claims need tempering and the wall-region electronic-response assumption needs either a justification or a caveat. I would send this to peer review rather than desk-reject, and I would ask the authors to address the validation circularity and the wall-region assumption explicitly. For a reader in ferroelectric photonics or phase-field modeling, this is worth a look; I'd probably cite it for the electronic-polarization formulation, but not for the 4000 pm/V number until that region is better supported.","headline":"A genuinely useful step: phase-field with an explicit electronic polarization gives spatially resolved refractive-index and electro-optic maps from evolving domain structure; the dramatic domain-wall numbers are the least secure part and the film-averaged validation is partly circular.","tokens_in":16303,"tokens_out":1481,"would_cite":true,"duration_ms":17234,"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":"This paper claims that a phase-field model with an added electronic polarization field can predict spatially resolved refractive-index and electro-optic maps in ferroelectric microstructures, reproducing measured BaTiO3 thin-film electro-op","keywords":["ferroelectrics","electro-optic coefficients","phase-field simulation","electronic polarization","refractive index","BaTiO3","domain walls","photonics"],"falsifier":"Measure the local electro-optic coefficient at a single a/c domain wall in a BaTiO3 thin film with near-field optical microscopy; the model predicts a wall-localized r coefficient above 4000 pm/V decaying over nanometers, so an absent or much weaker wall enhancement would falsify the central claim.","tokens_in":15316,"feed_emoji":"⚡","tokens_out":6845,"duration_ms":56600,"temperature":0.7,"pith_summary":"The paper establishes a computational method for predicting how the optical properties of a ferroelectric material change as its domain structure evolves. It does this by introducing an electronic polarization field—the part of the polarization that can respond at optical frequencies—into a phase-field simulation, so the local refractive index and electro-optic coefficients are computed directly from the evolving lattice polarization and stress. The payoff is that mesoscale phenomena that were previously inaccessible—domain-wall-localized enhancement, transient spikes during polarization switching, and phase-coexistence effects—become predictable. Applied to BaTiO3 thin films, the model reproduces experimentally measured electro-optic coefficients on silicon and explains temperature- and orientation-dependent trends that bulk single-crystal values cannot. If correct, it gives device designers a way to simulate photonic response alongside domain engineering, rather than splicing bulk optical tensors onto assumed domain patterns.","feed_headline":"Ferroelectric domain walls push electro-optic response past 4000 pm/V","feed_subtitle":"New phase-field model maps local refractive index and electro-optic response in BaTiO3 films, reproducing measured temperature trends.","key_machinery":"The engine of the model is the electronic polarization field P_e(x,t), the part of the polarization carried by distortion of electron orbitals that can respond at optical frequencies. The optical susceptibility at each point is set by an electronic dielectric stiffness tensor B_e(x) = ε0 ∂²f_e/∂P_e ∂P_e, which is the curvature of the free-energy landscape with respect to P_e. Because B_e is evaluated at the local lattice polarization and stress, the refractive-index map inherits the domain microstructure directly; the electro-optic coefficient is then obtained numerically as the finite difference of the inverse optical dielectric tensor between two applied fields. The separation of timescale","core_discovery":"The central claim is that the optical response of a ferroelectric is not set by the bulk electro-optic tensor assigned to artificial domains, but emerges from the coupled evolution of lattice polarization, stress, and an electronic polarization field. In this formulation, the local refractive index is determined by the inverse electronic dielectric stiffness, so every domain wall, phase boundary, and monoclinic bridging phase leaves a mark on the optical properties. The paper demonstrates this for BaTiO3 thin films: local electro-optic coefficients exceed 4000 pm/V near domain walls—about three times the bulk single-crystal r_51 = 1300 pm/V—and during polarization switching the film-average","pith_inferences":["A testable design corollary the paper leaves implicit: deliberately stabilizing low-symmetry monoclinic phases or phase coexistence—rather than maximizing the stability of a single phase—could be the most effective route to large electro-optic response in thin films.","The predicted wall-localized r > 4000 pm/V suggests that films engineered with dense, stable 90° domain-wall arrays could show bulk-averaged electro-optic coefficients far above single-crystal values without relying on transient switching states.","Because the model assumes instantaneous local electronic equilibrium, an extension that includes electronic gradient energy or nonlocal response at walls could either soften or sharpen the predicted wall enhancement; this remains an open question the author's approach does not settle.","The same framework could be used to design quasi-phase-matched nonlinear devices by simulating the domain pattern directly and computing the resulting nonlinear coefficient map, connecting microstructure simulation to device layout without manual domain assignment."],"forward_implications":["The effective electro-optic coefficient of a multidomain ferroelectric film can exceed the bulk single-crystal value, so device design cannot rely on volume-averaged bulk tensors.","Dense or movable domain-wall configurations are predicted to create local electro-optic hotspots; engineering wall populations is a lever for enhancing the response.","Near ferroelectric phase boundaries—where phases coexist and low-symmetry monoclinic phases mediate transitions—the average electro-optic response peaks at roughly 2.5 times the mid-tetragonal value.","The model quantitatively matches measured temperature and field-angle dependence of BaTiO3/Si films, with the lower experimental magnitudes attributable to interfacial dead layers or incomplete poling.","The same electronic-polarization-field machinery extends beyond the linear electro-optic effect to nonlinear susceptibilities, thermo-optic and piezo-optic coefficients, and frequency-dependent optical response."],"fun_headline_variants":["Domain walls triple electro-optic response in BaTiO3 films","New phase-field model maps local optical response in ferroelectrics","Coupled model predicts temperature-dependent optics in ferroelectric domains","Electro-optic effect near domain walls hits 4000 pm/V","Phase-field approach captures domain-wall-enhanced electro-optics"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that at optical frequencies the electronic polarization responds instantly and locally to the lattice polarization and stress, with the electronic gradient energy neglected; if bulk-derived coupling constants fail or nonlocal electronic effects matter at domain walls, the predicted >4000 pm/V enhancements and temperature dependence would not survive.","fun_headline_variants_meta":{"raw":{"variants":["Domain walls triple electro-optic response in BaTiO3 films","New phase-field model maps local optical response in ferroelectrics","Coupled model predicts temperature-dependent optics in ferroelectric domains","Electro-optic effect near domain walls hits 4000 pm/V","Phase-field approach captures domain-wall-enhanced electro-optics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000421,"raw_usage":{"total_tokens":2016,"prompt_tokens":774,"completion_tokens":1242,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":518,"completion_tokens_details":{"reasoning_tokens":1155}},"tokens_in":518,"tokens_out":1242,"duration_ms":9115,"temperature":1.0,"reasoning_tokens":1155,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T21:06:12.135714+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the local electro-optic coefficient at a single a/c domain wall in a BaTiO3 thin film with near-field optical microscopy; the model predicts a wall-localized r coefficient above 4000 pm/V decaying over nanometers, so an absent or much weaker wall enhancement would falsify the central claim.","supporting_citations":[],"review_version":1}