{"id":"7f234aec-b04f-48d3-9aa8-ce5c80b63205","arxiv_id":"2504.21246","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"Atomistic simulations and Kramers theory connect activation entropy to a predicted increase of critical resolved shear stress with temperature in copper, but the paper's own fitted equations contradict the claimed increase for edge dislocations.","lead":"Simulations of copper dislocations show that the rate at which dislocation pairs bypass each other changes with thermostat friction, and that an entropy term that grows with temperature can make flow stress rise instead of fall. The paper offers a candidate explanation for a long-standing anomaly in copper and other metals, but the supporting equations inside the paper do not consistently reproduce the claimed trend.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (15) contradicts the central claim: edge CRSS falls from 11.97 MPa at 0 K to 4.04 MPa at 300 K, and Table 6 disagrees with Table 5 by roughly 9 eV at 100 MPa.","rationale":"The reader's verdict of REJECT is well supported, but the most load-bearing weakness is not the unvalidated Kramers mass in Eqs. (10)-(11); it is the internal contradiction between the fitted CRSS polynomials (Table 6 and Eqs. 15-16) and the directly computed activation free energies (Table 5). The central claim requires the CRSS to increase with temperature, but the edge expression that the paper itself presents decreases monotonically over the anomalous regime and becomes negative near 387 K. The fits in Table 6 also disagree with Table 5 by an order of magnitude in the energy scale, which means the zeroing procedure that defines CRSS is operating on polynomials that do not represent the computed barriers. A corrected fit through the Table 5 edge data yields a much higher 0 K CRSS and a different temperature trend, so the entropy-based explanation of the flow-stress anomaly is not established. The MD-versus-Kramers friction-rate comparison may still be a useful contribution, but it does not rescue the macroscopic CRSS conclusion. I therefore agree with the rejection while identifying a different, more direct internal inconsistency as the decisive concern.","tokens_in":16616,"tokens_out":6251,"duration_ms":63377,"concrete_test":"Recompute CRSS(T) directly from Table 5: fit Delta G(T, tau) to the reported values (screw: 0.718 eV at 30 MPa and 0.331 eV at 40 MPa; edge: 1.468 eV at 100 MPa and 0.483 eV at 120 MPa) with a linear stress term, extrapolate each fit to Delta G = 0, and compare the resulting tau_CRSS(T) with Fig. 9. Also evaluate Eq. (15) numerically at 0, 150, and 300 K; if the values are 11.97, about 9.3, and 4.04 MPa as written, the edge CRSS decreases with temperature and the claimed entropic increase is not present.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline claim that activation entropy raises the CRSS with temperature is derived by setting the fitted activation free energies in Table 6 to zero, yielding Eqs. (15) and (16). The edge formula Eq. (15) is monotonically decreasing over the experimentally anomalous range: 11.97 MPa at 0 K, about 9.3 MPa at 150 K, 4.04 MPa at 300 K, and it crosses zero near 387 K. The screw formula decreases slightly from 0 K to about 300 K before rising, and the average therefore does not reproduce the claimed entropic increase at low-to-intermediate temperatures. More seriously, Table 6 is internally inconsistent with Table 5, which reports the directly computed activation free energies. For the edge dipole at 100 MPa and 0 K, Table 5 gives about 1.468 eV, whereas the Table 6 expression gives 1.005 - 8.394 = -7.389 eV; at 120 MPa the discrepancy is similar. A two-point linear fit through the edge entries of Table 5 gives dG/dtau about -0.049 eV/MPa and a 0 K CRSS near 130 MPa, not the 11.97 MPa from Eq. (15). Thus the polynomial used to produce the central CRSS curve contradicts the paper's own computed barriers. The Kramers mass assumption flagged by the reader is also unvalidated, but the CRSS conclusion is already unsupported by the paper's own data.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports atomistic molecular dynamics (MD) and nudged-elastic-band (NEB) calculations of thermally activated edge and screw dislocation dipole bypass in face-centered cubic copper, using the Mishin potential. Rates obtained from 465 MD simulations are compared with Kramers rate theory at several Langevin friction coefficients, temperatures, and shear stresses. The authors also compute activation entropies along the NEB paths using Schoeck's formalism and construct temperature- and stress-dependent activation free energies. The paper's headline claim is that activation entropy produces a critical resolved shear stress (CRSS) that increases with temperature, thereby explaining the low-to-intermediate temperature flow-stress anomaly in copper and other FCC metals.","tokens_in":16978,"tokens_out":8317,"duration_ms":82004,"significance":"If correct, the observation that Langevin friction changes only the Kramers prefactor and not the activation barrier would be a useful contribution, and the entropy-based explanation of the flow-stress anomaly would be significant because it would provide a mechanism distinct from impurities, cross-slip, and dynamic strain aging. The direct MD/Kramers comparison across friction coefficients (Figs. 7 and 8, Table 7) is a genuine strength, and the use of Schoeck's entropy with temperature-dependent elastic constants is methodologically interesting. However, the central CRSS claim is not supported by the manuscript's own equations: the fitted activation free energies in Table 6 are inconsistent with the directly computed barriers in Table 5, and Eq. (15) for the edge dislocation gives a CRSS that decreases with temperature over the anomalous range. The rate-versus-friction portion of the paper is separable and more robust, but the advertised explanation of the anomaly fails in its present form.","major_comments":[{"comment":"Equation (15) contradicts the stated conclusion of a slight increase of CRSS with temperature. Evaluating the edge expression gives 11.97 MPa at 0 K, 10.54 MPa at 100 K, 7.83 MPa at 200 K, and 4.04 MPa at 300 K, and the expression crosses zero near 387 K. The screw expression, Eq. (16), decreases slightly from 39.71 MPa at 0 K to about 38.6 MPa at 300 K before rising slowly. The average of the edge and screw CRSS is therefore monotonically decreasing over the experimentally anomalous range below 300 K, in direct contradiction to the text and to Fig. 9, which claim that the theoretical CRSS follows the same increasing trend as the experimental flow stress.","section":"§4.1, Eq. (15) and Fig. 9"},{"comment":"Table 6 is internally inconsistent with the directly computed activation free energies in Table 5. For the edge dipole at τ = 100 MPa and T = 0 K, Table 5 gives ΔG_eff ≈ 1.468 eV, whereas the Table 6 expression gives 1.005 − 8.394×10⁻² × 100 = −7.389 eV. At 120 MPa the discrepancy is similar. Since the anharmonic activation free energy ΔG_nh at T = 0 reduces to the activation enthalpy, a negative value is unphysical. The polynomial used to derive Eqs. (15) and (16) therefore does not represent the paper's own NEB/Schoeck data, and the CRSS curves built on it are not a valid prediction.","section":"§4.1, Table 5 versus Table 6"},{"comment":"The temperature derivatives implied by Table 6 do not support the proposed entropy mechanism. For the edge dipole at 100 MPa, ∂ΔG/∂T at T = 0 is approximately −9.5×10⁻⁵ eV/K, and for the screw dipole at 30 MPa it is approximately −6.6×10⁻⁴ eV/K. Thus the fitted activation free energy decreases with temperature at these conditions, while the text states that negative activation entropy increases the activation free energy and leads to an increasing τCRSS. The fitted polynomials and the qualitative mechanism asserted in Section 4.1 are mutually inconsistent.","section":"§4.1, Table 6 temperature derivatives"},{"comment":"The CRSS is defined by setting the fitted activation free energy polynomial to zero and extrapolating far outside the simulated stress window (30–120 MPa) and to temperatures below those used in the MD and NEB calculations. No uncertainty estimate, validation against direct barrier calculations, or sensitivity analysis is provided for this extrapolation. This is particularly problematic because the polynomial already fails to reproduce the barriers in Table 5 at the stress values on which it is ostensibly based.","section":"§4.1, CRSS extrapolation"}],"minor_comments":[{"comment":"The caption reads 'Activation free energy per as a function of temperature' and is missing a word; it should say what the free energy is normalized by.","section":"Table 5 caption"},{"comment":"The labels '1st jump' and '2nd jump' for the edge dislocation are not defined in the text; please explain what these two barriers correspond to physically.","section":"Tables 3 and 4"},{"comment":"The figure mixes computed CRSS values with experimental flow-stress data from different mechanisms and alloys; the caption should state the plotted quantity, the stress units, and whether the comparison is quantitative or only qualitative.","section":"Fig. 9"},{"comment":"The Kramers prefactor uses the total mass of all atoms in the simulation cell as the effective mass along the reaction coordinate; this assumption is not tested or justified, and the text should at least discuss its sensitivity, even though it does not affect the CRSS determination.","section":"§2.4, Eqs. (10) and (11)"}],"recommendation":"reject","confidential_remarks":"The rate-versus-friction comparison and the entropy calculations are potentially useful, but the paper's advertised central claim is contradicted by its own equations and tables. I would be willing to reconsider a substantially revised manuscript in which the CRSS analysis is corrected or the anomaly claim is removed; in its current form the internal inconsistency between Tables 5 and 6 and Eq. (15) makes the main conclusion unsupported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's headline claim—that activation entropy raises the critical resolved shear stress with temperature in copper—does not survive contact with its own equations. Eq. (15), which the authors use to define τ_CRSS for edge dislocations, gives 11.97 MPa at 0 K, 10.54 MPa at 100 K, 7.83 MPa at 200 K, and 4.04 MPa at 300 K. That is a monotonic decrease, the opposite of the abstract's promise. Figure 9 shows an increasing trend, but it appears to be drawn from the other root of the fit or from an averaging artifact, not from Eq. (15). The stress-test note also catches a serious internal inconsistency: plugging T=0, τ=100 MPa into the Table 6 expression for the edge gives −7.39 eV, while Table 5 reports 1.468 eV for the same conditions. A ~9 eV discrepancy means the polynomial used to derive the CRSS curve is not consistent with the paper's own computed free energies. The argument that activation entropy explains the anomalous flow stress is therefore not supported by the data as presented.\n\nThat said, the paper has real value in its secondary contributions. The MD simulations of edge and screw dipole bypass under different Langevin friction coefficients, the comparison with Kramers rate theory, and the demonstration that friction shifts the prefactor without changing the barrier are solid and reasonably convincing. The Kramers rates match the MD trends to within an order of magnitude for most cases, and the non-Arrhenius behavior at high temperature near τ_CRSS is an interesting observation, even if the temperature-dependent entropy improvement is mixed (the MSLE at 40 MPa is 11–27, which is poor). The Schoeck entropy application to dipole interactions is a legitimate extension of prior work, and the paper is honest about the limitations of a single-dipole configuration.\n\nThe soft spots beyond the CRSS issue are the unvalidated choice of total cell mass in the Kramers prefactor, the narrow friction range, and the fact that the CRSS result is an extrapolation of the same fitted polynomial that is internally inconsistent. These are fixable in principle, but as it stands the central claim is wrong on its own terms.\n\nWho should read this? Anyone working on atomistic rate calculations for dislocation-obstacle interactions will find the friction-dependence study useful, and the non-Arrhenius behavior is worth knowing about. But the paper should not be published with the current conclusions. It deserves a serious referee because the underlying data and methods are mostly sound, but the framing needs major revision and the CRSS analysis must be corrected. My recommendation: reject in current form, but encourage a resubmission that fixes the fitting inconsistency and either retracts the CRSS claim or presents it accurately. This is a case where the paper is not a lost cause, but the headline is wrong.","headline":"The main CRSS claim is contradicted by the paper's own Eq. (15) and the Table 5/6 inconsistency; the friction-rate study is salvageable, but the central explanation of the flow stress anomaly fails.","tokens_in":858,"tokens_out":850,"would_cite":false,"duration_ms":38162,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Activation entropy alone could explain why copper's flow stress rises with temperature.","keywords":["activation entropy","flow stress anomaly","critical resolved shear stress","dislocation dipole","Kramers rate theory","Schoeck entropy formalism","Langevin friction","copper"],"falsifier":"Perform the same dipole-bypass rate calculations with the effective mass extracted from the actual normal mode at the initial and saddle points instead of the total cell mass; if the resulting rates and the temperature dependence of the critical resolved shear stress differ substantially from the paper's predictions, the central claim is falsified.","tokens_in":16384,"feed_emoji":"⚛️","tokens_out":7039,"duration_ms":65874,"temperature":0.7,"pith_summary":"This paper claims that activation entropy, not impurities or cross-slip, can explain the long-standing anomaly in which copper's flow stress increases with temperature at low to intermediate temperatures. Using atomistic simulations of edge and screw dislocation dipoles, the authors compute the rate at which dislocations overcome their mutual interaction barrier and compare it with Kramers rate theory. They find that the temperature dependence of the elastic constants makes the activation entropy strongly negative and growing with temperature, so the activation free energy barrier rises with temperature when the enthalpy barrier is small. This entropic effect raises the critical resolved shear stress with temperature, matching the experimentally observed trend. The paper further shows that Langevin friction modifies only the attempt frequency, not the barrier, a result confirmed by molecular dynamics.","feed_headline":"Entropy may explain why copper's flow stress rises with temperature","feed_subtitle":"Atomistic simulations trace the anomaly to temperature-dependent activation entropy in pure Cu.","key_machinery":"The load-bearing machinery is Schoeck's entropy formalism, which expresses the change in entropy induced by internal strains as $\\Delta S = \\alpha_V K \\int_\\Omega V_{ii}\\,d\\Omega - \\frac{1}{2}\\int_\\Omega \\frac{\\partial C_{iklm}}{\\partial T} V_{ik} V_{lm}\\,d\\Omega$, in combination with nudged elastic band (NEB) enthalpy profiles along the minimum energy path. This yields the activation free energy $\\Delta G = \\Delta H - T\\Delta S$ at each stress and temperature. The rate is then evaluated with Kramers' expression, whose attempt frequency depends on the curvature of the enthalpy at the initial and saddle points and on the Langevin friction coefficient. The temperature dependence of the elastic constants (and the thermal expansion term, found negligible here) makes the entropy temperature dependent, and that is the mechanism driving the anomalous increase of $\\tau_{\\mathrm{CRSS}}$ with temperature.","core_discovery":"The central discovery is that for both edge and screw dislocation dipoles in face-centered cubic copper, the activation entropy computed with Schoeck's elasticity formalism depends on temperature and stress, is negative, and grows in magnitude with temperature. When this temperature-dependent entropy is inserted into the activation free energy $\\Delta G = \\Delta H - T\\Delta S$, the barrier that a dislocation must overcome increases with temperature, especially at stresses approaching the critical resolved shear stress at 0 K. Consequently, the critical resolved shear stress $\\tau_{\\mathrm{CRSS}}$ increases with temperature, reproducing the anomalous flow-stress behavior observed experimentally. Kramers' transition state theory, with harmonic prefactors derived from the enthalpy curvature at the minimum and saddle points, reproduces the rates measured by molecular dynamics and correctly captures the effect of Langevin friction, which acts only on the prefactor. The non-Arrhenius rate behavior at high temperature is captured only when the entropy is treated as temperature dependent, not when it is held constant.","pith_inferences":["If activation entropy is the driver, then alloying or pressure that changes $\\partial C_{iklm}/\\partial T$ should systematically shift the magnitude and temperature range of the flow-stress anomaly, a prediction that goes beyond the paper's pure-copper dipoles.","The same Schoeck-plus-Kramers machinery could be applied to other face-centered cubic metals (Al, Ni, Ag) to test whether the anomaly's magnitude scales with the temperature sensitivity of their elastic constants.","A testable extension is to compute the rates with a different thermostat (e.g., Nose-Hoover) or with explicit heat baths of varying coupling; if the friction dependence vanishes or changes sign, the Kramers prefactor identification would need revision."],"forward_implications":["Kramers rate theory reproduces the MD-computed rates and their dependence on Langevin friction, so friction enters only through the attempt frequency and leaves the energy barrier unchanged.","At stresses near the zero-temperature critical resolved shear stress, the bypass rate can become non-Arrhenius: increasing temperature can reduce the rate, because the entropic barrier grows faster than $k_B T$.","The critical resolved shear stress for dislocation dipole bypass in pure copper increases with temperature, matching the direction of the anomaly seen experimentally in Cu and Cu-based alloys.","No impurities, cross-slip, or dynamic strain aging are needed to produce the anomalous increase; activation entropy alone is sufficient in the simulated configurations."],"supporting_citations":[{"why":"Supplies the Kramers rate expression used to compute bypass rates and their dependence on Langevin friction.","marker":"[1]"},{"why":"Provides the previous methodology for computing activation entropy along the minimum energy path that this work builds on.","marker":"[22]"},{"why":"Gives the elasticity-based entropy formula (Schoeck's formalism) central to computing the activation entropy.","marker":"[39]"},{"why":"Provides the nudged elastic band method used to determine the minimum energy path and enthalpy barriers.","marker":"[40]"},{"why":"Supplies the interatomic potential for copper used in the MD simulations and for the temperature-dependent elastic constants.","marker":"[35]"},{"why":"Provides experimental data on the increase of critical resolved shear stress with temperature below 200 K in Cu-Co alloys that the results are compared against.","marker":"[11]"},{"why":"Supplies the comprehensive Kramers turnover theory and the high/low-friction limits used in the rate expressions.","marker":"[32]"}],"fun_headline_variants":["Copper's odd strength rise linked to activation entropy","Temperature-dependent entropy flips copper's hardening trend","Why copper gets stronger when heated: entropy's role","Entropy explains copper's anomalous thermal hardening","Activation entropy: culprit behind copper's stress-temperature puzzle"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The Kramers rate calculation treats the dislocation bypass as a one-dimensional crossing whose attempt frequencies use the total mass of all atoms in the simulation cell, so if the true effective mass along the reaction path differs, every predicted rate and the extracted critical stress change.","fun_headline_variants_meta":{"raw":{"variants":["Copper's odd strength rise linked to activation entropy","Temperature-dependent entropy flips copper's hardening trend","Why copper gets stronger when heated: entropy's role","Entropy explains copper's anomalous thermal hardening","Activation entropy: culprit behind copper's stress-temperature puzzle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000157,"raw_usage":{"total_tokens":1197,"prompt_tokens":897,"completion_tokens":300,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":513,"completion_tokens_details":{"reasoning_tokens":226}},"tokens_in":513,"tokens_out":300,"duration_ms":3554,"temperature":1.0,"reasoning_tokens":226,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:10:08.448640+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform the same dipole-bypass rate calculations with the effective mass extracted from the actual normal mode at the initial and saddle points instead of the total cell mass; if the resulting rates and the temperature dependence of the critical resolved shear stress differ substantially from the paper's predictions, the central claim is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Kramers rate expression used to compute bypass rates and their dependence on Langevin friction."},{"cited_title":"Schoeck, The entropy of internal stress fields, physica status solidi (b) 97 (1) (1980) 345–354","cited_arxiv_id":null,"evidence_quote":"Gives the elasticity-based entropy formula (Schoeck's formalism) central to computing the activation entropy."},{"cited_title":"Fusenig, E","cited_arxiv_id":null,"evidence_quote":"Provides experimental data on the increase of critical resolved shear stress with temperature below 200 K in Cu-Co alloys that the results are compared against."}],"review_version":1}