{"id":"3a6b843a-ffb4-47b6-a953-3d321d023a21","arxiv_id":"2608.06759","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"DFPT-based DFT+U+V with a carefully chosen Hubbard manifold reproduces birefringence in BaTiO3, LiNbO3, KNbO3 and PbTiO3, and predicts strong anisotropy enhancement in Bi-Mn doped BaTiO3.","lead":"Using density-functional theory with extended Hubbard corrections, the authors show that intersite U+V terms restore the optical anisotropy of ferroelectric perovskite oxides that plain onsite DFT+U suppresses. The result matters because a cheap, non-empirical functional that gets birefringence right could be used to screen crystals for photonic and quantum-optics applications.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"PbTiO3 geometry error undermines benchmark: relaxed c/a=1.150 vs 1.063, internal displacements ~50% too large, and the claimed birefringence agreement depends on choosing Ref [43] over Ref [37].","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the DFPT set II PbTiO3 geometry is over-distorted, so the birefringence agreement may be partly accidental. My reading sharpens this in two ways. First, Table 2 makes the magnitude concrete: c/a = 1.150 versus 1.063 and oxygen displacements about 50% too large. Second, the manuscript itself documents conflicting experimental birefringence data for PbTiO3, with DFPT following the lower curve of Ref. [37] while the higher curve of Ref. [43] is invoked as consistent with the calculation. A method claim that depends on choosing one of two experimental reports, at a geometry that is far outside the experimental structure, cannot support 'systematically validated' predictive performance for this compound. The independent-particle dielectric approximation (Eq. 11) is also a limitation, but it is a shared approximation across all four compounds and the anisotropy comparison may partly cancel common errors; the PbTiO3 geometry issue is a direct, table-level inconsistency with the paper's own accuracy claim. The concern is addressable rather than fatal: a recomputation at the experimental geometry would show whether the electronic-structure method or the relaxed structure is responsible for the reported agreement. If the birefringence remains high at the experimental geometry, the central claim survives this test; if it drops, the benchmark needs revision. Because the reader already assigned CONDITIONAL, and this test is exactly the kind of condition that should be required, the verdict should remain conditional pending this check.","tokens_in":14773,"tokens_out":4510,"duration_ms":43430,"concrete_test":"Recompute the PbTiO3 birefringence with DFPT set II Hubbard parameters but at the experimental structure (a = 3.904 Å, c/a = 1.063, internal coordinates from Ref. [41]) instead of the relaxed structure listed in Table 2. If the computed Δn drops from Ref. [43]-like values toward the Ref. [37] Sellmeier fit, the apparent agreement is an artifact of the overestimated ferroelectric distortion. A complementary run at the relaxed c/a = 1.150 but with experimental internal displacements would isolate whether the overestimated c/a or the overestimated oxygen off-centering is responsible for the birefringence overshoot.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that DFPT set II relaxations are accurate enough for birefringence to be predictive. Table 2 violates this for PbTiO3: c/a = 1.150 versus the experimental 1.063 (~8% error), and the internal displacements are overestimated by roughly 50% (for example, ΔzO1 = −0.671 Å versus −0.456 Å, and ΔzO3 = −0.656 Å versus −0.456 Å). Birefringence in a uniaxial ferroelectric is a strong function of tetragonal distortion, so agreement between DFPT set II and the PbTiO3 data of Ref. [43] does not discriminate between an accurate electronic-structure method and an over-polarized structure whose exaggerated anisotropy is compensated by other errors. The paper itself acknowledges two conflicting experimental reports: the DFPT results closely follow the lower-birefringence Sellmeier fit of Ref. [37] below 500 nm, while Ref. [43] reports an approximately fourfold larger birefringence in that region. Selecting the dataset that matches a calculation performed at a substantially wrong geometry leaves the predictive claim unsupported for one of the four benchmark compounds. This is load-bearing because the paper's stated conclusion is that the DFPT set II functional 'reproduces the experimental crystal geometries and optical anisotropies with good overall accuracy' across all four materials; PbTiO3 is a direct counterexample for geometry, and the optical agreement is not independent of that geometry error. The KNbO3 exception is an additional caveat, but the PbTiO3 case is more fundamental because the metric being validated (optical anisotropy) is computed at a structure that is far from the experimental one.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript benchmarks DFT+U and DFT+U+V against the anisotropic optical response of BaTiO3, LiNbO3, KNbO3, and PbTiO3, with Hubbard parameters computed non-empirically using DFPT and, for comparison, ACBN0, for three orbital manifolds: {U_TM,U_O}, {U_TM,V_TM-O}, and {U_TM,U_O,V_TM-O}. It reports that the DFPT set-II functional (TM onsite plus TM-O intersite) gives the best refractive indices and birefringence, proposes two scalar phase-matching metrics, and applies the method to Zn- and Bi/Mn-substituted BaTiO3, predicting enhanced birefringence and dichroism from Mn-induced mid-gap states. The paper's central claim is that DFPT set II reproduces both the experimental crystal geometries and optical anisotropies of the four benchmark oxides at mean-field cost.","tokens_in":15096,"tokens_out":6748,"duration_ms":66585,"significance":"The work is significant in its methodological comparison and in providing concrete, falsifiable predictions, notably the blue-shifted Type-I SHG phase-matching frequency and enhanced visible-range dichroic ratio in Ba7BiTi7MnO24. Its main strength is that U and V are not fitted to optical data: the DFPT and ACBN0 parameter sets follow well-defined linear-response prescriptions, and the manifold comparison is systematic. The LiNbO3 birefringence agreement and the phase-matching landscape analysis are convincing. However, the benchmark evidence is not as uniform as the conclusion claims: the PbTiO3 geometry errors are large, the BaTiO3 internal displacements are also overestimated by roughly 50%, and the optical response omits excitonic and local-field effects. Because these issues bear directly on the central predictive claim, the paper needs major revision rather than acceptance.","major_comments":[{"comment":"Table 2 shows that DFPT set II relaxes PbTiO3 to c/a = 1.150 versus the experimental 1.063, an error of about 8%, and internal displacements are overestimated by roughly 44% (e.g., Delta z_O1 = -0.671 Å versus -0.456 Å, and Delta z_O3 = -0.656 Å versus -0.456 Å). Birefringence in a uniaxial ferroelectric is highly sensitive to the tetragonal distortion, so agreement with the birefringence data of Ref. [43] does not independently validate the electronic-structure treatment: it may reflect compensation between an over-polarized structure and other errors. The claim in the Assessment section that DFPT set II reproduces the experimental crystal geometries and optical anisotropies with good overall accuracy across the four benchmark oxides is directly contradicted by this row. Please add a sensitivity test, for example computing birefringence at the experimental geometry, or substantially weaken the geometry part of the claim for PbTiO3.","section":"Table 2, PbTiO3 row; Assessment of optical anisotropy calculations"},{"comment":"The geometry problem is not confined to PbTiO3. For BaTiO3, DFPT set II gives internal oxygen displacements Delta z_O1 = -0.145 Å and Delta z_O3 = -0.090 Å versus experimental values of -0.093 Å and -0.056 Å, i.e., roughly 50-60% too large, even though c/a happens to be close. This pattern suggests that the good BaTiO3 birefringence agreement may also involve compensation between an overestimated polar distortion and other errors. Please quantify the sensitivity of the computed birefringence to the ferroelectric distortion, for instance by repeating the optical calculation at the experimental lattice parameters and internal coordinates.","section":"Table 2, BaTiO3 row"},{"comment":"The dielectric tensor is evaluated in the independent-particle approximation with a damping parameter, and Eqs. (12)-(14) derive the refractive index and extinction coefficient from it; no excitonic or local-field corrections are included. The comparisons in Fig. 3 extend into spectral regions near and above the optical gap (for example PbTiO3 below 500 nm), where excitonic and local-field effects can change both the magnitude and the anisotropy of the response. To support the predictive claim, either estimate these corrections for at least one benchmark compound or explicitly restrict the claim to the transparency region where independent-particle response is expected to be reliable.","section":"Computational Methods, dielectric response"},{"comment":"The recommendation of set II {U_TM,V_TM-O} is made after comparing sets I-III against the same experimental birefringence curves shown in Fig. 3. The Hubbard U and V values themselves are non-empirical, but the choice of the orbital manifold is benchmark-informed. This is a form of model selection on the target data and should be stated as a limitation. The predictive claim would be strengthened by demonstrating transferability on a held-out compound or property, or by a physical criterion that selects set II without reference to the optical benchmarks.","section":"Influence of the selected projection manifold"},{"comment":"The paper acknowledges that KNbO3 is the only exception, with ACBN0 agreeing better with experiment because of compensating errors, yet the conclusion still states that DFPT set II reproduces the measured anisotropic optical properties across the benchmark compounds. A claim of consistency across four materials should be weakened to three of four, or the KNbO3 failure should be analyzed more quantitatively, for example by decomposing the birefringence error into contributions from the gap error and the distortion error.","section":"Assessment of optical anisotropy calculations; Conclusions"}],"minor_comments":[{"comment":"The code name is misspelled as 'Qunatum ESPRESSO'; it should read 'Quantum ESPRESSO'.","section":"Computational Methods"},{"comment":"The dielectric-response formula in the Computational Methods is not numbered, but Eqs. (12)-(14) refer to its real and imaginary parts and to the refractive index; please number it consistently.","section":"Computational Methods, dielectric response"},{"comment":"The phase-matching metrics in Eqs. (10) and (11) use a Gaussian broadening parameter sigma that is not defined in the table caption for Table 4; please state the value and the integration domain explicitly at the point where the metrics are introduced.","section":"Birefringent phase-matching landscapes"},{"comment":"The text discusses two conflicting experimental birefringence datasets for PbTiO3 (Ref. [37] and Ref. [43]) but does not clearly state which one is plotted in Fig. 3. Please specify the dataset in the figure caption and justify the selection.","section":"Assessment of optical anisotropy calculations; Figure 3"},{"comment":"The data availability statement only offers data 'on a reasonable request'; for a benchmarking study of this type, depositing structures, Hubbard parameters, and computed dielectric tensors in a public repository would significantly improve reproducibility.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the stress-test concern about PbTiO3 lands in full. The manuscript's central claim is overstated relative to the evidence: one of the four benchmark compounds has a large geometry error, a second has overestimated internal displacements, and the optical response is computed without excitonic or local-field corrections. I do not see grounds for rejection, because the methodology is well-defined, the Hubbard parameters are genuinely non-empirical, the LiNbO3 result is convincing, and the solid-solution predictions are falsifiable. A revision that weakens the four-material claim, adds a geometry-sensitivity analysis, and addresses the independent-particle approximation would make the paper publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Name],\n\nHere's my take on arXiv:2608.06759. The paper does something useful: it benchmarks a non-empirical DFT+U+V scheme against the optical anisotropy of four ferroelectric perovskites and shows that the choice of Hubbard manifold matters a lot. The DFPT-computed U and V are not fitted to optical data, and the comparison across three parameter sets is systematic. The two phase-matching metrics (fraction and renormalized Dirichlet energy) are a reasonable addition, and the Bi-Mn doped BaTiO3 prediction is interesting, even if speculative.\n\nThe soft spot is PbTiO3. The relaxed c/a is 1.150 against the experimental 1.063, and the oxygen off-centerings are roughly 50% too large. Birefringence in a uniaxial ferroelectric is extremely sensitive to that distortion, so the agreement with the Zekri et al. data (Ref. [43]) is not a clean test. The paper acknowledges that Ref. [37] gives a much lower birefringence, and then argues the DFPT results match the fourfold-larger values from Ref. [43]. That choice of experiment needs an independent argument, not just \"consistent with our predictions.\" Without it, the PbTiO3 benchmark carries little weight. KNbO3 is also an acknowledged exception, so two of the four compounds don't clearly support the \"good overall accuracy\" claim.\n\nThe computational setup is standard and likely reproducible, but the data are not deposited, which is a minor issue. The independent-particle dielectric function without excitonic or local-field corrections is a known limitation, and the paper doesn't overclaim about it.\n\nOverall, the paper is a solid methodological contribution, but the central conclusion is overstated. A serious referee should ask for a discussion of the PbTiO3 geometry error and a more careful justification of the experimental data selection, or a softened claim. Worth sending to peer review.\n\nBest,\n[Your name]","headline":"Useful benchmark of DFT+U+V for optical anisotropy, but the PbTiO3 geometry error and a cherry-picked experimental comparison undercut the central claim.","tokens_in":15673,"tokens_out":3703,"would_cite":false,"duration_ms":34900,"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":"A low-cost DFT fix reproduces four ferroelectrics' birefringence.","keywords":["DFT+U+V","optical anisotropy","birefringence","perovskite oxides","Hubbard parameters","density-functional perturbation theory","phase matching","ferroelectric"],"falsifier":"Fix the relaxed geometry to the experimental lattice parameters for PbTiO3 (c/a = 1.063 instead of the calculated 1.150) and recompute birefringence with the same DFT+U+V set-II parameters; if the large birefringence that matches Ref. [43] collapses, the central agreement is an artifact of an over-distorted geometry rather than of the functional.","tokens_in":14576,"feed_emoji":"🔬","tokens_out":4557,"duration_ms":41953,"temperature":0.7,"pith_summary":"The paper tries to establish that the extended Hubbard functional DFT+U+V, with Hubbard parameters computed by density-functional perturbation theory and applied only to transition-metal onsite plus transition-metal–oxygen intersite terms, predicts the anisotropic optical response of BaTiO3, LiNbO3, KNbO3, and PbTiO3 in agreement with measured refractive indices and birefringence. The result matters because it offers a computationally cheap route to predictive optical anisotropy and phase-matching calculations without hybrid functionals or many-body perturbation theory. The paper also shows that plain DFT overestimates refractive indices while DFT+U alone suppresses anisotropy, and that the choice of orbital manifold is decisive. As a demonstration, it predicts that Bi–Mn co-substituted BaTiO3 develops a fourfold larger birefringence and a large visible dichroic ratio, with implications for polarization-sensitive photodetectors and integrated photonics.","feed_headline":"A low-cost DFT fix reproduces four ferroelectrics' birefringence","feed_subtitle":"Adding transition-metal–oxygen Hubbard terms with DFPT parameters matches measured refractive indices and flags doped BaTiO3 for photonics.","key_machinery":"The load-bearing object is the extended Hubbard functional DFT+U+V, whose correction energy adds an onsite U term and an intersite V term between atom-centered orbital occupations. The parameters are obtained non-empirically as second derivatives of the total energy with respect to orbital occupancies through density-functional perturbation theory, which includes electronic screening. The decisive choice is the projection manifold: the paper finds that set II, with onsite Hubbard terms on transition metals plus intersite transition-metal–oxygen terms but no onsite oxygen U, reproduces the measured anisotropy. The optical response is then evaluated from the independent-particle momentum-space dielectric function, and phase-matching behavior is summarized by two new dimensionless metrics, the phase-matching fraction and the renormalized Dirichlet energy of the refractive-index-mismatch function.","core_discovery":"The central claim is that correcting self-interaction errors with DFT+U alone is insufficient for non-centrosymmetric perovskite oxides because it localizes electrons too strongly and suppresses optical anisotropy, whereas adding intersite Hubbard V terms between transition-metal and oxygen orbitals restores the covalent hybridization that controls birefringence. Using DFPT-derived Hubbard parameters within the manifold {UTM, V_TM-O}, the paper reproduces the measured refractive-index dispersion and birefringence of BaTiO3, LiNbO3, KNbO3, and PbTiO3, identifies LiNbO3 as having the most favorable SHG phase-matching landscape among these, and extends the same parameter-free framework to substituted BaTiO3 solid solutions. The paper further claims that ACBN0, an alternative parameterization, systematically overestimates ferroelectric distortions and birefringence for these oxides, and that including onsite Hubbard U on oxygen within DFT+U+V destroys the anisotropy because oxygen orbitals respond weakly to the applied perturbation, making their U parameters spuriously large.","pith_inferences":["The paper's PbTiO3 agreement leans on the birefringence data of Ref. [43] rather than the Sellmeier fit of Ref. [37]; if experimental consensus shifts toward the fit, the claimed agreement for PbTiO3 would need revisiting.","Because birefringence in these crystals is extremely sensitive to the ferroelectric distortion, the same functional should be tested on other phases and on compounds with smaller distortions to see whether the set-II manifold remains optimal.","The independent-particle dielectric function omits excitonic and local-field effects; comparing these DFT+U+V spectra against BSE-level spectra would show whether the predicted anisotropy survives those corrections."],"forward_implications":["If the central claim is right, DFT+U+V with DFPT parameters and the set-II manifold can be used to screen ferroelectric oxides for birefringence and phase matching at a fraction of the cost of hybrid or GW calculations.","Plain DFT overestimates refractive indices and DFT+U underestimates anisotropy in these materials, so uncorrected or U-only calculations should not be trusted for optical-device design in this class.","The phase-matching fraction and renormalized Dirichlet energy give dimensionless scalar descriptors for ranking materials by phase-matching bandwidth in screening workflows.","Bi–Mn co-substituted BaTiO3 is predicted to combine a fourfold birefringence increase with Type-I SHG phase matching over the 0.75–1.0 eV range, making it a concrete candidate for polarization-sensitive photodetection and integrated photonics."],"supporting_citations":[{"why":"Supplies the DFPT linear-response method used to compute non-empirical Hubbard U and V parameters.","marker":"[21]"},{"why":"Provides the hp code used to carry out the DFPT Hubbard-parameter calculations.","marker":"[25]"},{"why":"Gives the ACBN0 extended-Hubbard benchmark whose birefringence results are compared against DFPT.","marker":"[30]"},{"why":"Establishes the linear-response definition of the onsite Hubbard U that underlies the DFPT derivation.","marker":"[33]"},{"why":"Supports the attribution of oxygen overlocalization to weak oxygen linear response, explaining why U on oxygen suppresses anisotropy.","marker":"[23]"},{"why":"Provides experimental refractive-index and birefringence data for BaTiO3 used as a benchmark.","marker":"[34]"},{"why":"Provides experimental refractive-index data for LiNbO3 used as a benchmark.","marker":"[35]"},{"why":"Provides the PbTiO3 Sellmeier fit used for comparison of refractive-index dispersion.","marker":"[37]"},{"why":"Provides the alternative PbTiO3 birefringence measurements that the DFPT predictions follow.","marker":"[43]"}],"fun_headline_variants":["Intersite Hubbard V restores perovskite birefringence","DFT+U+V reproduces four oxides' optical anisotropy","Low-cost correction nails perovskite birefringence","Hubbard V term fixes DFT+U's anisotropy failure","Accurate perovskite optics from DFT+U+V, cheaply"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole predictive chain assumes the Hubbard-corrected relaxed ferroelectric geometry is close enough to experiment and that omitting excitonic and local-field corrections does not change the anisotropy; PbTiO3's overestimated c/a of 1.150 versus the experimental 1.063 shows how tight that assumption is.","fun_headline_variants_meta":{"raw":{"variants":["Intersite Hubbard V restores perovskite birefringence","DFT+U+V reproduces four oxides' optical anisotropy","Low-cost correction nails perovskite birefringence","Hubbard V term fixes DFT+U's anisotropy failure","Accurate perovskite optics from DFT+U+V, cheaply"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000336,"raw_usage":{"total_tokens":1898,"prompt_tokens":1022,"completion_tokens":876,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":802}},"tokens_in":638,"tokens_out":876,"duration_ms":8716,"temperature":1.0,"reasoning_tokens":802,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:12:55.596288+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fix the relaxed geometry to the experimental lattice parameters for PbTiO3 (c/a = 1.063 instead of the calculated 1.150) and recompute birefringence with the same DFT+U+V set-II parameters; if the large birefringence that matches Ref. [43] collapses, the central agreement is an artifact of an over-distorted geometry rather than of the functional.","supporting_citations":[],"review_version":1}