{"id":"28dd544a-b84c-4170-a3ab-2a3f005b9252","arxiv_id":"2506.05054","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"At high laser fluence, (Ba,Sr)TiO3 films expand less than heat diffusion predicts, a deviation the authors attribute to charge injection across the Schottky barrier at the SrRuO3 interface.","lead":"This paper compares the pyroelectric response of a barium strontium titanate thin film measured by standard hysteresis loops with picosecond structural deformations seen in time-resolved X-ray diffraction after a laser pulse. It finds that at higher laser powers the film expands less than a heat diffusion model predicts, suggesting a transient charge effect at the metal/ferroelectric interface that could hide part of the pyroelectric response.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The high-fluence 'frustrated expansion' claim depends on a constant-material-parameter heat diffusion prediction, but Section 5 concedes temperature-dependent parameters cannot be neglected; a temperature-dependent thermal expansion coefficient could explain the apparent frustration.","rationale":"The reader's weakest assumption correctly identifies that the heat diffusion simulation, after fluence scaling, must correctly predict the BST temperature and strain for the central claim to hold. My stress-test focuses on a more specific and manuscript-supported failure mode: the conversion of simulated temperature to lattice shift via Eq. (3) uses a single thermal expansion coefficient, while Section 5 explicitly states that temperature-dependent material parameters cannot be neglected for the ~100 K temperature rise. The low-fluence agreement does not constrain this nonlinearity, and the static p(T) data in Section 2 show strong temperature dependence near the diffuse transition. The paper's own limitation statement thus undercuts the premise that the high-fluence BST strain prediction is reliable enough to attribute the deficit to charge injection. This is not a disagreement with consensus; it is an internal consistency concern between the simulation's constant-parameter treatment and the paper's stated limitation. The proposed concrete test would settle the question by replacing the constant alpha with the measured temperature-dependent alpha and checking whether the apparent frustration disappears. Because the reader already rendered a CONDITIONAL verdict based largely on this class of concern, my analysis does not change the verdict; it sharpens the condition that should be met before the charge mechanism is accepted. I agree with the reader's weakest assumption and find no need to adjust the verdict.","tokens_in":12715,"tokens_out":5469,"duration_ms":74121,"concrete_test":"Re-run the published udkm1dsim heat diffusion model with temperature-dependent alpha_BST(T) measured on the same film by static XRD over 20-150 degrees C, and with literature c_p(T) if available, while keeping the fluence scaling fixed by the SrRuO3 sensor fit. Compare predicted and measured BST (002) shifts at 10 and 15 mJ/cm2. If the revised model matches the measured data, the frustration is explained by thermal nonlinearity and the charge-injection mechanism is unnecessary; if the discrepancy persists, the charge-injection interpretation survives this test.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that at 10 and 15 mJ/cm2 the (Ba,Sr)TiO3 film expands less than the heat diffusion prediction, with the deficit assigned to charge injection across the SrRuO3/(Ba,Sr)TiO3 Schottky barrier. The load-bearing premise is that the simulation's BST strain prediction is thermally accurate after scaling the fluence to match the SrRuO3 sensor. Section 5 explicitly concedes this premise is unsafe: for the ~100 K temperature rise, 'the temperature dependence cannot be neglected' for the material parameters, and no quantitative assessment is made. Eq. (3) converts simulated temperature to lattice shift with a single thermal expansion coefficient alpha; the heat diffusion input likewise uses literature or single-temperature values for c_p and k. The low-fluence (6 mJ/cm2) agreement calibrates one operating point and does not constrain behavior over a ~100 K rise in a clamped film with a diffuse transition. If the true alpha at high temperature is smaller than the value used, the constant-parameter model will overpredict the BST expansion, producing exactly the reported frustration without any charge mechanism. The Schottky interpretation is also unsupported by measured barrier height or injected charge density, as the authors acknowledge by proposing DLTS; the 'frustration' evidence is therefore not yet distinguished from a thermal nonlinearity artifact.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports static pyroelectric measurements of epitaxial (Ba,Sr)TiO3 films on NdScO3 with a SrRuO3 bottom electrode, deriving the pyroelectric coefficient from temperature-dependent polarization hysteresis loops and an indirect electrocaloric temperature change. The authors then perform time-resolved x-ray diffraction after femtosecond laser excitation of the SrRuO3 layer and compare the transient lattice expansion of BST and SrRuO3 to one-dimensional heat diffusion simulations. At low fluence (6 mJ/cm2) the simulation matches the BST expansion after scaling the effective fluence to 4.5 mJ/cm2, whereas at higher fluences (10 and 15 mJ/cm2, scaled to 65% and 45% respectively) the BST expands less than the constant-parameter heat diffusion model predicts. This 'frustrated thermal expansion' is tentatively attributed to charge injection across the SrRuO3/BST Schottky barrier, which would counteract thermal expansion via the inverse piezoelectric effect. The authors explicitly note that quantitative assessment of the secondary pyroelectric effect is currently impossible due to unknown temperature dependence of material parameters.","tokens_in":12990,"tokens_out":3114,"duration_ms":39368,"significance":"If the central claim were fully supported, the paper would be significant because it proposes a non-contact, time-resolved route to disentangle secondary pyroelectric contributions in clamped ferroelectric films, and it identifies a possible transient charge-screening mechanism at the electrode/ferroelectric interface. The static pyroelectric and electrocaloric characterization of Ba0.7Sr0.3TiO3 is useful and follows standard procedures. The time-resolved methodology, including the use of SrRuO3 as an embedded temperature sensor, is also of interest. However, the load-bearing 'frustration' claim rests on a constant-parameter heat diffusion model whose quantitative accuracy over the ~100 K temperature rise is explicitly conceded in Section 5 to be inadequate, and the Schottky mechanism is proposed without measured barrier heights or injected charge densities. The paper is honest about these limitations, but as it stands the main new finding is not yet established beyond a plausible alternative (temperature-dependent thermal expansion coefficients).","major_comments":[{"comment":"The 'frustrated thermal expansion' claim is defined as a deviation of the measured BST strain from a heat diffusion simulation that uses constant (or literature single-temperature) values of the thermal expansion coefficient, specific heat, and thermal conductivity. Section 5 explicitly states that for the ~100 K temperature rise in the experiment 'the temperature dependence cannot be neglected' for the relevant material parameters. If the true thermal expansion coefficient of BST decreases with temperature over this range (as is common near a diffuse transition), the constant-parameter model will overpredict the expansion and produce exactly the reported frustration without any charge mechanism. To make the central claim load-bearing, the authors must either implement temperature-dependent parameters in the simulation, or quantitatively bound the effect of their temperature dependence and show that the observed deviation exceeds this bound.","section":"Section 3, Fig. 3 and Section 5"},{"comment":"The effective fluence is separately tuned to 65% (for 10 mJ/cm2) and 45% (for 15 mJ/cm2) so that the simulated SrRuO3 sensor response matches the measurement. The same calibrated model is then used to predict the BST expansion. This introduces a two-parameter fit that can absorb not only the stated reduction in energy transfer but also any systematic error in the thermal model (e.g., incorrect absorption, interface resistance, or parameter temperature dependence). No physical mechanism is given for why the scaling factor changes with fluence, and no independent constraint is provided for the BST temperature prediction. Please provide an uncertainty analysis that propagates the fluence scaling uncertainty into the predicted BST strain, or an alternative validation (e.g., a second sensor layer or independent temperature probe).","section":"Section 3, fluence scaling"},{"comment":"The measured transient qz shifts are shown without error bars, and the simulated curves are presented as deterministic. The magnitude of the claimed BST deviation at 10 and 15 mJ/cm2 is not quantified relative to experimental precision (e.g., peak position uncertainty from the reciprocal space map fits, laser fluence stability, x-ray time jitter). Please provide error bars or confidence intervals on at least the final plateau values of the BST shift at each fluence, and state whether the deviation is statistically significant.","section":"Section 3, Figs. 3a-c"},{"comment":"The Schottky-barrier charge injection mechanism is presented as an explanation for the frustrated expansion, but no quantitative estimate is given of the space charge density, barrier height, or injected charge required to produce the observed strain compensation. The manuscript even notes that no quantitative values for the relevant (Ba,Sr)TiO3/SrRuO3 barrier parameters are reported in the literature. Without an order-of-magnitude demonstration that plausible charge injection can generate the required counteracting strain (via inverse piezoelectric coupling), the mechanism remains only a speculation. A quantitative estimate, even with wide uncertainty ranges, is needed to support the interpretation over the thermal nonlinearity alternative.","section":"Section 4, Eq. (4) and proposed mechanism"}],"minor_comments":[{"comment":"The abstract and conclusion use 'frustration of the thermal expansion' whereas the body text in Section 3 and 6 sometimes says 'reduced thermal expansion.' Please adopt a single term for consistency, preferably one that does not presume a non-thermal origin.","section":"Abstract and Conclusion"},{"comment":"Equation (3) reads '∆qz = qz,0 (α∆T + 1)', which is dimensionally inconsistent because ∆qz and qz,0 have the same units but the right-hand side is qz,0 times a dimensionless factor plus qz,0. This should presumably be qz = qz,0(1 + α∆T) or ∆qz/qz,0 = α∆T. Please correct.","section":"Eq. (3)"},{"comment":"The text for the 6 mJ/cm2 measurement states the simulation reproduces the data with an excitation fluence of 4.5 mJ/cm2, while the Figure 3d caption says the simulation was performed for 6.5 mJ/cm2. These numbers are inconsistent; please clarify which fluence values correspond to which curves and why they differ from the nominal experimental fluences.","section":"Section 3, Fig. 3 caption"},{"comment":"The caption defines the substrate as NdGaO3 (NGO) while the text consistently uses NdScO3 (NSO). Please correct the figure caption or the text so the substrate material is reported consistently.","section":"Figure 1 caption"},{"comment":"There is a typo in the text: 'SrRu03' should be 'SrRuO3' in the sentence on work functions. Also, in the same section the phrase 'permeability in (Ba,Sr)TiO3' likely means 'permittivity' or 'dielectric properties'; please revise the wording.","section":"Section 4"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a promising experimental approach but the central claim is not yet supported because the heat diffusion model's constant-parameter approximation is acknowledged in Section 5 to be inadequate over the relevant temperature range. The authors will need to either supply temperature-dependent simulations or a quantitative bounding argument. The ad hoc fluence scaling factors also weaken the model-based interpretation. If the authors can address the thermal nonlinearity alternative with a quantitative sensitivity analysis, the manuscript could become suitable for publication. I also note that the static pyroelectric and electrocaloric measurements appear solid but are not themselves the main novelty; the paper's fit to Journal of Alloys and Compounds may hinge on the strength of the time-resolved result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe short version: this paper combines routine pyroelectric hysteresis measurements with TR-XRD on the same BST/SRO/NSO stack and reports a fluence-dependent under-expansion of the BST film relative to a heat diffusion model. If the observation holds, it is a nice new transient interface effect; the Schottky screening explanation, however, is hand-waving, and the paper admits as much.\n\nWhat is new: the TR-XRD study on the same material used for static pyroelectric characterization, and the specific deviation at 10 and 15 mJ/cm². The low-fluence data match the simulation, which gives some confidence in the model mechanics. The static pyroelectric coefficient of about -170 µC/m²K and the derived electrocaloric ΔT of 0.37 K are plausible and consistent with literature. The paper is transparent about its limits: Section 5 explicitly says no quantitative assessment of the secondary pyroelectric effect is possible because temperature-dependent piezoelectric and elastic stiffnesses are unknown, and DLTS is suggested as a follow-up.\n\nThe soft spot is the load-bearing assumption that the heat diffusion model predicts the BST strain accurately at high fluence. The model uses constant α, c_p, and k. The fluence is scaled to 65% and 45% of the nominal values to match the SRO sensor. That the deviation in BST appears despite this scaling is suggestive, but the model's prediction for BST depends on extrapolating material parameters over a ~100 K temperature rise. The stress-test concern is on point: if the true α_BST drops with temperature, the constant-α model would overpredict expansion and produce exactly this 'frustration' without any charge effect. The paper's own Section 5 concedes the temperature dependence cannot be neglected. So the central observation is real, but the mechanism is not distinguished from a thermal nonlinearity artifact. There are also no error bars on the TR-XRD shifts or the static p values, which makes the strength of the deviation hard to judge.\n\nWhere does that leave us? The paper is worth a serious referee. The measurements are careful, the writing is honest, and the observation, even if it turns out to be thermal, is a useful data point for people modeling ultrafast heating in ferroelectric stacks. The referee should push for error bars, a temperature-dependent α test, and ideally a direct probe of charge dynamics. I would send it to review rather than desk reject, but the authors should be asked to tone down the Schottky interpretation or support it with more evidence.\n\nI'd probably bring it to a reading group if the topic is ferroelectrics; for a general audience it is a bit narrow. I'd cite it if I were working on transient lattice response in oxide heterostructures.\n\nBest,\n[You]","headline":"Solid, honest TR-XRD study of BST/SRO with a fluence-dependent expansion anomaly, but the Schottky explanation is under-supported and the heat diffusion calibration needs scrutiny.","tokens_in":13512,"tokens_out":3522,"would_cite":true,"duration_ms":40451,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["77.70.-a"],"model":"deepseek-v4-flash","headline":"Picosecond x-ray diffraction shows a (Ba,Sr)TiO3 film expanding less than heat diffusion predicts, a deficit the authors attribute to charge injection across the SrRuO3/(Ba,Sr)TiO3 Schottky barrier.","keywords":["ferroelectric thin film","time-resolved x-ray diffraction","pyroelectric effect","ultrafast lattice dynamics","Schottky barrier","heat diffusion simulation","barium strontium titanate","pump-probe"],"falsifier":"Repeat the high-fluence pump-probe measurement on an identical (Ba,Sr)TiO3 stack with a non-Schottky bottom electrode (or with a DC bias applied to suppress charge injection) and check whether the expansion deficit disappears; if it does not disappear, the charge-screening mechanism is not the cause.","tokens_in":1808,"feed_emoji":"⚛️","tokens_out":3077,"duration_ms":154204,"temperature":0.7,"pith_summary":"This paper is trying to establish that picosecond x-ray diffraction can separate the thermal from the non-thermal parts of the pyroelectric response in a ferroelectric thin film. The authors heat a (Ba,Sr)TiO3 film through an underlying SrRuO3 electrode with femtosecond laser pulses and watch the lattice expand on a picosecond timescale. At low pump fluence the expansion matches a heat diffusion simulation; at higher fluence the film expands less than the model predicts, and the paper attributes that 'frustrated thermal expansion' to electrons injected across the Schottky barrier at the SrRuO3/(Ba,Sr)TiO3 interface, whose screening field opposes the expansion through the inverse piezoelectric effect. If correct, this gives a contactless way to access pyroelectric contributions that ordinary hysteresis measurements average together.","feed_headline":"X-rays show ferroelectric film expanding less than heat predicts","feed_subtitle":"Strain snapshots hint hot electrons cross the electrode barrier, revealing a hidden pyroelectric term.","key_machinery":"The central object is the Schottky barrier at the SrRuO3/(Ba,Sr)TiO3 interface, a contact barrier formed because the metal and ferroelectric work functions differ, with a space-charge region of width $d_n = \\sqrt{2\\epsilon_r\\epsilon_0(\\phi_M-\\phi_S)/(e n_D)}$ extending into the ferroelectric film. The argument is carried by comparing measured transient Bragg-peak shifts to a one-dimensional heat diffusion simulation of the heterostructure. At low fluence the simulation matches both layers and serves as a calibration; at high fluence the SrRuO3 sensor is matched only after rescaling the input fluence, while the (Ba,Sr)TiO3 film still falls short of the predicted thermal expansion. The shortfall is then interpreted as an inverse piezoelectric strain generated by the transient electric field of an electron-screened, extended space-charge region.","core_discovery":"The paper's central claim is that the picosecond deformation of (Ba,Sr)TiO3 after femtosecond heating of the SrRuO3 electrode is not purely thermal: at excitation fluences of 10 and 15 mJ/cm2, the ferroelectric film expands less than a heat diffusion simulation predicts, even though the same simulation reproduces the SrRuO3 temperature sensor after scaling the fluence to 65% and 45%. The authors propose that this 'frustrated thermal expansion' arises because hot electrons from the laser-heated SrRuO3 cross the Schottky barrier, screen the positive space charge, and thereby change the electric field in the space-charge region. The field change produces an inverse piezoelectric strain that counteracts thermal expansion, and the injected carriers also alter the permittivity and polarization. The mechanism is explicitly tentative: the paper states that a quantitative separation of pyroelectric contributions requires temperature-dependent piezoelectric and elastic constants and better characterization of the interface, for example by deep-level transient spectroscopy.","pith_inferences":["If the charge-injection picture is right, the strain deficit is a direct time-resolved measure of hot-electron transfer across an oxide interface, so the same experiment could probe Schottky-barrier dynamics in other metal/ferroelectric stacks without electrical contacts.","The empirical fluence rescaling (65% and 45%) could also be explained by reduced optical absorption at high intensity; measuring the pump reflectivity and absorbed energy at each fluence would discriminate between a real charge effect and an absorption artifact.","Because the inverse piezoelectric strain should be proportional to the space-charge field, applying an external bias across the stack in a future experiment would modulate the magnitude or sign of the frustrated expansion if the mechanism is correct."],"forward_implications":["At low fluence the heat diffusion model reproduces both layers, so the experiment supplies a calibrated thermal baseline; any strain shortfall at higher fluence is a direct signature of non-thermal response.","Static hysteresis measurements give only the sum of all pyroelectric contributions, whereas the transient strain data isolate the strain-related secondary channel, even though a quantitative value is not yet reached.","The frustration appears only above a threshold fluence and is absent in static measurements, implying the extra dissipation channel is a transient non-equilibrium process.","A quantitative separation would require temperature-dependent piezoelectric and elastic stiffness constants plus a characterization of the Schottky barrier, which the paper identifies as the necessary next step."],"supporting_citations":[{"why":"Supplies the finite-element heat diffusion simulation used to compute transient layer temperatures.","marker":"[26]"},{"why":"Models the ultrafast laser heating regime and justifies the delta-like thermal excitation.","marker":"[8]"},{"why":"Supports neglecting in-plane heat diffusion in the first 10 ns after excitation.","marker":"[9]"},{"why":"Provides the relation between Bragg-peak shift and temperature used to read the SrRuO3 thermometer.","marker":"[10]"},{"why":"Supplies the growth and structural characterization of the (Ba,Sr)TiO3 films.","marker":"[12]"},{"why":"Provides the decomposition of the pyroelectric coefficient into primary, secondary, tertiary, and quaternary parts and the secondary contribution formula.","marker":"[5]"},{"why":"Gives the sub-500 fs electron-phonon coupling time in SrRuO3 that supports impulsive heating.","marker":"[30]"},{"why":"Supplies the space-charge region width formula used to describe the Schottky barrier.","marker":"[32]"},{"why":"Provides the work-function values for SrRuO3 needed to estimate the barrier height.","marker":"[34]"},{"why":"Reports ballistic energy and charge transport across a metal/absorber interface, supporting the proposed hot-electron injection.","marker":"[40]"}],"fun_headline_variants":["Picosecond strain reveals frustrated thermal expansion in ferroelectric","Hot electrons from electrode curb ferroelectric expansion","Laser-heated film expands less than heat model predicts","Time-resolved X-rays spot hidden pyroelectric contribution","X-ray snapshots show ferroelectric expanding less than heat predicts"],"cache_read_input_tokens":15616,"weakest_assumption_plain":"The interpretation stands or falls on the assumption that the heat diffusion simulation, after scaling the laser fluence to match the SrRuO3 temperature transient, correctly predicts the temperature inside the (Ba,Sr)TiO3 layer, so the expansion shortfall must come from a non-thermal process rather than from an error in the model or in the temperature-dependent material parameters.","fun_headline_variants_meta":{"raw":{"variants":["Picosecond strain reveals frustrated thermal expansion in ferroelectric","Hot electrons from electrode curb ferroelectric expansion","Laser-heated film expands less than heat model predicts","Time-resolved X-rays spot hidden pyroelectric contribution","X-ray snapshots show ferroelectric expanding less than heat predicts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000312,"raw_usage":{"total_tokens":1783,"prompt_tokens":960,"completion_tokens":823,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":742}},"tokens_in":576,"tokens_out":823,"duration_ms":9872,"temperature":1.0,"reasoning_tokens":742,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:26:26.131685+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the high-fluence pump-probe measurement on an identical (Ba,Sr)TiO3 stack with a non-Schottky bottom electrode (or with a DC bias applied to suppress charge injection) and check whether the expansion deficit disappears; if it does not disappear, the charge-screening mechanism is not the cause.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports neglecting in-plane heat diffusion in the first 10 ns after excitation."},{"cited_title":"Shayduk, H","cited_arxiv_id":null,"evidence_quote":"Provides the relation between Bragg-peak shift and temperature used to read the SrRuO3 thermometer."},{"cited_title":"Wawra, K","cited_arxiv_id":null,"evidence_quote":"Supplies the growth and structural characterization of the (Ba,Sr)TiO3 films."},{"cited_title":"Mattern, J.-E","cited_arxiv_id":null,"evidence_quote":"Gives the sub-500 fs electron-phonon coupling time in SrRuO3 that supports impulsive heating."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the space-charge region width formula used to describe the Schottky barrier."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the work-function values for SrRuO3 needed to estimate the barrier height."},{"cited_title":"Plech, P","cited_arxiv_id":null,"evidence_quote":"Reports ballistic energy and charge transport across a metal/absorber interface, supporting the proposed hot-electron injection."}],"review_version":1}