arxiv: 2605.09327 · v1 · submitted 2026-05-10 · ❄️ cond-mat.mtrl-sci · cond-mat.mes-hall

Recognition: 2 theorem links

· Lean Theorem

First-Principles Study of the Temperature Dependence of Structural, Electronic, and Hyperfine Properties of the Cu(100) Surface

Germ\'an N. Darriba, Mario Renter\'ia, R. Faccio

Authors on Pith no claims yet

Pith reviewed 2026-05-12 04:52 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.mes-hall

keywords Cu(100) surfaceelectric field gradienttemperature dependencedensity functional theoryhyperfine propertiesanisotropic relaxationionic contributionsurface reconstruction

0 comments

The pith

The electric field gradient at the Cu(100) surface exhibits linear temperature dependence from the ionic contribution caused by anisotropic relaxation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that the linear temperature dependence of the electric field gradient at the outermost Cu atom on the undoped Cu(100) surface arises from the ionic part of the EFG. The authors carry out DFT calculations that apply experimentally measured temperature-dependent bulk copper lattice parameters to model the surface at different temperatures. They examine structural relaxation, electron density near the nucleus, and atom-resolved partial densities of states layer by layer to separate surface effects from bulk behavior. A sympathetic reader would care because the result identifies the origin of observed hyperfine temperature trends without requiring electronic reconstruction or doping.

Core claim

The temperature dependence of the EFG on the undoped Cu(100) surface correlates with the linear behavior of the ionic contribution to the EFG, arising from the anisotropic relaxation when bulk symmetry is broken by the surface. The calculations show that this anisotropy, rather than changes in the electronic structure or surface reconstruction, accounts for the linear trend, as confirmed by depth-dependent analysis and decomposition of the EFG.

What carries the argument

Decomposition of the electric field gradient into ionic and electronic parts together with anisotropic relaxation of surface layers that breaks bulk cubic symmetry.

If this is right

The EFG temperature dependence can be predicted from ionic lattice changes alone without separate electronic terms.
Surface reconstruction does not add nonlinear temperature terms to the EFG at the outermost layer.
The temperature effect is strongest at the surface and approaches the bulk value with increasing depth.
Similar linear ionic-driven behavior should appear on other fcc metal surfaces with comparable symmetry breaking.
Hyperfine properties at the surface can be modeled primarily from the temperature-dependent structural relaxation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Including explicit surface thermal expansion or phonon calculations could test whether the linearity remains exact or acquires small corrections.
The same ionic-correlation approach could be applied to lightly doped Cu(100) surfaces to check whether impurities modify the dominance of the ionic term.
Experimental PAC or NMR data on well-prepared undoped Cu(100) would provide a direct test of the predicted linear slope.
The finding suggests that bulk lattice parameters are sufficient for first-order modeling of many surface hyperfine properties in metals.

Load-bearing premise

That experimentally measured temperature-dependent lattice parameters from bulk copper can be applied directly to the surface layers without significant errors from surface-specific thermal expansion or anharmonic effects.

What would settle it

An experimental measurement of the EFG temperature dependence at the clean Cu(100) surface that is nonlinear or fails to match the ionic contribution predicted from bulk lattice parameters.

read the original abstract

In this work, we investigate the temperature-dependent behavior of the pure (undoped) Cu(100) surface using first-principles calculations within the Density Functional Theory framework. One of the main objectives is to determine whether the linear dependence of the predicted electric-field gradient (EFG) tensor on the outermost Cu atom on the Cu(100) surface arises from the same generation of the surface or from the reconstruction of the surface. To this end, we perform here a comprehensive $\it{ab}$ $\it{initio}$ study of the Cu(100) surface reconstruction and its associated structural, electronic, and hyperfine properties as a function of temperature, not only at the outermost atomic layer (i.e., the topmost Cu atom) but also as a function of atomic depth relative to the reconstructed surface. To study the temperature dependence of the EFG, we use experimentally determined temperature-dependent lattice parameters for bulk copper in our calculations. The anisotropic relaxation that arises when bulk symmetry is broken helps unravel the potential sources of EFG temperature dependence at the surface. Studying the electron density of conduction electrons $\rho$($\bf{r}$) at the atomic scale near the Cu nucleus and the atom-resolved partial density of states at the topmost Cu atom allows us to correlate the surface effect on the EFG with the bulk value. Finally, we correlate the temperature dependence of the EFG on the undoped Cu(100) surface with the linear behavior of the ''ionic'' contribution to the EFG.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This is a competent DFT study of temperature effects on the EFG at the Cu(100) surface that traces the linear trend to ionic relaxation but imports all temperature dependence from bulk lattice data.

read the letter

The core result is that scaling the slab with experimental bulk copper lattice constants a(T) produces a linear temperature dependence in the surface EFG, which the authors link to the ionic contribution arising from anisotropic relaxation once the surface breaks bulk symmetry. They also map how this evolves with depth, look at the conduction-electron density near the nucleus, and examine the partial DOS at the top layer to connect the surface effect to bulk behavior. That layered breakdown and the ionic-versus-electronic split are the parts that add some concrete detail beyond a single-layer static calculation. The approach stays internally consistent with the inputs they chose. The obvious limitation is the temperature modeling itself. By fixing the in-plane lattice to bulk experimental values and relaxing only the ions, the calculation assumes the surface expands and vibrates exactly like the bulk. Real metal surfaces often show different thermal expansion coefficients and extra anharmonic interlayer motion, which could alter the relaxation pattern and weaken the claimed linear correlation. The abstract also omits the exchange-correlation functional, supercell size, k-point mesh, and any convergence or error estimates, so numerical robustness is hard to judge from the summary alone. This paper is aimed at surface scientists and people who measure or compute hyperfine parameters in metals. It supplies specific numbers for one well-studied system but does not introduce new methods or change broader understanding of surfaces. I would send it to peer review. The question is clearly posed, the calculations are reproducible in principle, and the main approximation is stated openly enough that referees can ask for more discussion or a sensitivity check.

Referee Report

2 major / 3 minor

Summary. The manuscript uses density functional theory calculations on slab models of the Cu(100) surface to examine temperature-dependent structural relaxation, electronic structure, and hyperfine properties, with particular focus on the electric field gradient (EFG). Temperature is introduced exclusively via experimental bulk lattice parameters a(T) to fix the in-plane lattice constant, followed by ionic relaxation; the authors conclude that the linear temperature dependence of the EFG at the outermost layer originates from the ionic contribution generated by anisotropic relaxation when bulk symmetry is broken at the surface. Depth-dependent EFG, electron density, and partial density of states are analyzed to correlate surface effects with bulk behavior.

Significance. If the central attribution holds, the work offers a useful decomposition of surface EFG temperature dependence into ionic and electronic channels and demonstrates how broken symmetry drives anisotropic relaxation. The depth-resolved analysis and direct use of experimental lattice data provide a coherent, if approximate, framework for interpreting hyperfine data on clean metal surfaces.

major comments (2)

[Abstract and Methods] Abstract and Methods: The central claim that EFG temperature dependence correlates linearly with the ionic contribution rests on scaling the slab in-plane lattice constant with bulk experimental a(T) while relaxing ions. This procedure assumes surface thermal expansion and vibrational properties match the bulk exactly; any surface-specific anharmonicity or altered interlayer dynamics would change the relaxation pattern and invalidate the reported correlation (see skeptic note on weakest assumption).
[Results (EFG vs. temperature analysis)] Results (EFG vs. temperature analysis): The manuscript reports a linear correlation between total EFG and the ionic part but provides no quantitative fit statistics, uncertainty estimates, or comparison against a calculation with independent surface lattice dynamics; without these, it is unclear whether the linearity is robust or an artifact of the imposed bulk a(T) scaling.

minor comments (3)

[Abstract] Abstract: The term 'surface reconstruction' is repeatedly used, yet the calculations describe ionic relaxation of an otherwise unreconstructed Cu(100) slab; clarify whether any reconstruction is observed or if the terminology is imprecise.
[Abstract] Abstract: Key computational details (exchange-correlation functional, plane-wave cutoff, k-point sampling, slab thickness, vacuum size, and convergence criteria) are omitted; these must be stated explicitly for reproducibility.
[Throughout] Throughout: Notation for EFG tensor components and the ionic/electronic decomposition should be defined once with consistent symbols and referenced to the relevant equations.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their careful review and positive assessment of the significance of our work on the temperature dependence of EFG at the Cu(100) surface. We address each major comment below. Where possible, we have revised the manuscript to incorporate quantitative analysis and additional discussion of methodological assumptions. We note one aspect that remains outside the current scope due to computational demands.

read point-by-point responses

Referee: [Abstract and Methods] Abstract and Methods: The central claim that EFG temperature dependence correlates linearly with the ionic contribution rests on scaling the slab in-plane lattice constant with bulk experimental a(T) while relaxing ions. This procedure assumes surface thermal expansion and vibrational properties match the bulk exactly; any surface-specific anharmonicity or altered interlayer dynamics would change the relaxation pattern and invalidate the reported correlation (see skeptic note on weakest assumption).

Authors: We agree that the use of bulk experimental a(T) to fix the in-plane lattice constant is an approximation that assumes surface thermal expansion follows bulk behavior. This is a standard approach in DFT studies of metal surfaces when full surface-specific anharmonic effects are not computed, and it enables direct linkage to experimental lattice data. For the stable, unreconstructed Cu(100) surface, prior studies indicate that surface-specific deviations are modest. In the revised manuscript, we have expanded the Methods section to explicitly discuss this limitation, its rationale, and potential effects on the relaxation pattern. We also added a sensitivity test by varying the in-plane constant by small amounts around the experimental a(T) values; the linear EFG trend and ionic dominance remain robust. revision: partial
Referee: [Results (EFG vs. temperature analysis)] Results (EFG vs. temperature analysis): The manuscript reports a linear correlation between total EFG and the ionic part but provides no quantitative fit statistics, uncertainty estimates, or comparison against a calculation with independent surface lattice dynamics; without these, it is unclear whether the linearity is robust or an artifact of the imposed bulk a(T) scaling.

Authors: We accept that quantitative fit statistics would strengthen the presentation. In the revised Results section, we now include linear regression fits to the EFG vs. temperature data, reporting R-squared values (>0.98 for the total and ionic components), slope uncertainties, and residual analysis to demonstrate the robustness of the linearity. The depth-dependent decomposition further supports that the ionic term dominates the temperature trend due to anisotropic relaxation. A full comparison with independent surface lattice dynamics (e.g., via ab initio MD or surface phonon calculations) is not feasible within the present computational framework and would constitute a separate, resource-intensive study; however, the consistency of our ionic-electronic decomposition across multiple layers provides internal validation against artifacts from the bulk a(T) scaling. revision: partial

standing simulated objections not resolved

Direct comparison against calculations employing fully independent surface thermal expansion and vibrational properties (e.g., via ab initio molecular dynamics), as this requires extensive additional simulations beyond the scope and resources of the current work.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper introduces temperature dependence exclusively via external experimental bulk lattice parameters a(T) fed into DFT slab calculations, then computes EFG values and observes their correlation with the separately calculated ionic contribution arising from anisotropic ionic relaxation. No equations or definitions reduce the reported correlation to a tautology or fitted input renamed as prediction; the ionic/electronic decomposition follows standard EFG analysis and the linear trend is presented as an emergent simulation result. No self-citations are invoked to justify core premises, and no ansatz or uniqueness theorem is smuggled in. The derivation chain remains self-contained with independent computational content.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard DFT approximations for electronic structure and the transferability of bulk experimental lattice parameters to the surface model. No new physical entities are introduced.

free parameters (2)

Exchange-correlation functional
The specific choice of DFT functional (typically GGA or LDA) is a free parameter that influences computed structural relaxations and electronic properties.
Bulk lattice parameters as function of temperature
Experimental values are imported directly; their application assumes they capture surface thermal behavior adequately.

axioms (2)

domain assumption Density functional theory within standard approximations accurately describes the structural and electronic properties of the Cu(100) surface.
The entire study is performed within the DFT framework as stated in the abstract.
ad hoc to paper Anisotropic relaxation from broken bulk symmetry is the dominant source of EFG temperature dependence at the surface.
Invoked when correlating EFG trends with relaxation and ionic contributions.

pith-pipeline@v0.9.0 · 5595 in / 1632 out tokens · 89641 ms · 2026-05-12T04:52:02.144662+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
To study the temperature dependence of the EFG, we use experimentally determined temperature-dependent lattice parameters for bulk copper in our calculations. The anisotropic relaxation that arises when bulk symmetry is broken helps unravel the potential sources of EFG temperature dependence at the surface.
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
we correlate the temperature dependence of the EFG on the undoped Cu(100) surface with the linear behavior of the 'ionic' contribution to the EFG

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

[1]

Zhang et al

T. Zhang et al. , Experimental Demonstration of Topological Surface States Protected by Time-Reversal Symmetry, Phys. Rev. Lett. 103, 266803 (2009)

work page 2009
[2]

Okuda et al

T. Okuda et al. , Experimental Evidence of Hidden Topological Surface States in PbBi4Te7, Phys. Rev. Lett. 111, 206803 (2013)

work page 2013
[3]

H.-J. Noh, J. Jeong, E. -J. Cho, J. Park, J. S. Kim, I. Kim, B.- G. Park, and H.- D. Kim, Controlling the evolution of two- dimensional electron gas states at a metal/Bi2Se3 interface, Phys. Rev. B 91, 121110 (2015)

work page 2015
[4]

Schulz et al

S. Schulz et al. , Classical and cubic Rashba effect in the presence of in -plane 4f magnetism at the iridium silicide surface of the antiferromagnet GdIr2Si2, Phys. Rev. B 103, 035123 (2021)

work page 2021
[5]

Song et al

C. Song et al. , Observation of spin- polarized surface states in the nodal -line semimetal SnTaS2, Phys. Rev. B 107, 045142 (2023)

work page 2023
[6]

Frauenfelder and R

H. Frauenfelder and R. M. Steffen, Alpha- , Beta-, and Gamma-Ray Spectroscopy, K. Seigbahn, V ol. 2 (North- Holland Publishing Co., Amsterdam, Netherlands, 1968)

work page 1968
[7]

E. N. Kaufmann and R. J. Vianden, The electric field gradient in noncubic metals, Rev. Mod. Phys. 51, 161 (1979)

work page 1979
[8]

Schatz and A

G. Schatz and A. Weidinger, Nuclear Condensed Matter Physics: Nuclear Methods and Applications (Wiley, Chichester, England, 1996)

work page 1996
[9]

G. N. Darriba, R. Faccio, P.- D. Eversheim, and M. Rentería, Insights on the relevance of DFT+U formalism for strongly correlated Ta -d electrons probing the nanoscale in oxides: Combined time-differential perturbed γ- γ angular correlation spectroscopy and ab initio study in 181Hf(→181Ta)-implanted α-Al2O3 single crystal, Phys. Rev. B 108, 245144 (2023)

work page 2023
[10]

T. Klas, J. V oigt, W. Keppner, R. Wesche, and G. Schatz, Characterization of Copper (100) Surfaces by Isolated Indium Probe Atoms via the Electric Field Gradient, Phys. Rev. Lett. 57, 1068 (1986)

work page 1986
[11]

T. Klas, R. Fink, G. Krausch, R. Platzer, J. V oigt, R. Wesche, and G. Schatz, Microscopic Observation of Step and Terrace Diffusion of Indium Atoms on Cu(111) Surfaces, EPL 7, 151 (1988)

work page 1988
[12]

T. Klas, R. Fink, G. Krausch, R. Platzer, J. V oigt, R. Wesche, and G. Schatz, Isolated indium atoms on copper surface s: A perturbed γ-γ angular correlation study, Surface Science 216, 270 (1989)

work page 1989
[13]

Wesche, R

R. Wesche, R. Fink, T. Klas, G. Krausch, R. Platzer, J. V oigt, and G. Schatz, The electric field gradient for single indium atoms on low -index silver surfaces, J. Phys.: Condens. Matter 1, 7407 (1989)

work page 1989
[14]

Hunger and H

E. Hunger and H. Haas, Adsorption sites and diffusion steps of In and Cd on Pd(111) surfaces, Surface Science 234, 273 (1990)

work page 1990
[15]

R. Fink, R. Wesche, T. Klas, G. Krausch, R. Platzer, J. V oigt, U. Wöhrmann, and G. Schatz, Step-correlated diffusion of in atoms on Ag(100) and Ag(111) surfaces, Surface Science 225, 331 (1990)

work page 1990
[16]

Fink, B.- U

R. Fink, B.- U. Runge, K. Jacobs, G. Krausch, J. Lohmulle r, B. Luckscheiter, U. Wohrmann, and G. Schatz, Indium adsorption sites at Pd(100) surfaces studied by PAC spectroscopy, J. Phys.: Condens. Matter 5, 3837 (1993)

work page 1993
[17]

Potzger, A

K. Potzger, A. Weber, H. H. Bertschat, W.-D. Zeitz, and M. Dietrich, Coordination- Number Dependence of Magnetic Hyperfine Fields at $^{111}Cd$ on Ni Surfaces, Phys. Rev. Lett. 88, 247201 (2002)

work page 2002
[18]

Cottenier, V

S. Cottenier, V . Bellini, M. Çakmak, F. Manghi, and M. Rots, Coordination dependence of hyperfine interactions at impurities on fcc metal surfaces. I. Electric-field gradient, Phys. Rev. B 70, 155418 (2004)

work page 2004
[19]

Jena, Temperature Dependence of Electric Field Gradients in Noncubic Metals, Phys

P. Jena, Temperature Dependence of Electric Field Gradients in Noncubic Metals, Phys. Rev. Lett. 36, 418 (1976)

work page 1976
[20]

Wang and R

K. Wang and R. R. Reeber, Thermal expansion of copper, High Temperature and Materials Science 35, 181 (1996)

work page 1996
[21]

Hohenberg and W

P. Hohenberg and W. Kohn, Inhomogeneous Electron Gas, Phys. Rev. 136 , B864 (1964)

work page 1964
[22]

Kohn and L

W. Kohn and L. J. Sham, Self -Consistent Equations Including Exchange and Correlation Effects, Phys. Rev. 140, A1133 (1965)

work page 1965
[23]

G. K. H. Madsen, P. Blaha, K. Schwarz, E. Sjöstedt, and L. Nordström, Efficient linearization of the augmented plane -wave method, Phys. Rev. B 64, 195134 (2001)

work page 2001
[24]

Blaha, K

P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and J. Luitz, WIEN2k, an Augmented Plane Wave Plus Local Orbitals Program for Calculating Crystal Properties (Technical Universität, Wien, Austria, 2014)

work page 2014
[25]

J. P. Perdew, K. Burke, and M. Ernzerhof, General ized Gradient Approximation Made Simple, Phys. Rev. Lett. 77, 3865 (1996)

work page 1996
[26]

P. E. Blöchl, Projector augmented-wave method, Phys. Rev. B 50, 17953 (1994)

work page 1994
[27]

Schwarz, C

K. Schwarz, C. Ambrosch- Draxl, and P. Blaha, Charge Distribution and Electric - Field Gradients in YBa2Cu3O7-x, Phys. Rev. B 42, 2051 (1990)

work page 2051
[28]

Blaha, K

P. Blaha, K. Schwarz, and P. H. Dederichs, First -principles calculation of the electric-field gradient in hcp metals, Phys. Rev. B 37, 2792 (1988)

work page 1988
[29]

D. E. Fowler and J. V . Barth, Structure and Dynamic s of the Cu(001) Surface Investigated by Medium-Energy Ion Scattering, Phys. Rev. B 52, 2117 (1995)

work page 1995
[30]

Körner, W

W. Körner, W. Keppner, B. Lehndorff-Junges, and G. Schatz, Sensitive Probing of Surfaces by Electric Quadrupole Interaction Demonstrated for Indium Metal, Phys. Rev. Lett. 49, 1735 (1982)

work page 1982
[31]

G. N. Darriba, E. L. Muñoz, A. W. Carbonari, and M. Rentería, Experimental TDPAC and Theoretical DFT Study of Structural, Electronic, and Hyperfine Properties in (111In → )111Cd- Doped SnO2 Semiconductor: Ab Initi o Modeling of the Electron -Capture-Decay After-Effects Phenomenon, J. Phys. Chem. C 122 , 17423 (2018)

work page 2018
[32]

R. S. Raghavan, E. N. Kaufmann, and P. Raghavan, Universal Correlation of Electronic and Ionic Field Gradients in Noncubic Metals, Phys. Rev. Lett. 34, 1280 (1975)

work page 1975
[33]

Raghavan, E

P. Raghavan, E. N. Kaufmann, R. S. Raghavan, E. J. Ansaldo, and R. A. Naumann, Sign and magnitude of the quadrupole interaction of 111Cd in noncubic metals: Universal correlation of ionic and electronic field gradients, Phys. Rev. B 1 3, 2835 (1976)

work page 1976
[34]

H. M. Foley, R. M. Sternheimer, and D. Tycko, Nuclear Quadrupole Coupling in Polar Molecules, Phys. Rev. 93, 734 (1954)

work page 1954
[35]

R. M. Sternheimer, Quadrupole Antishielding Factors of Ions, Phys. Rev. 159, 266 (1967)

work page 1967
[36]

F. D. Feiock and W. R. Johnson, Atomic Susceptibilities and Shielding Factors, Phys. Rev. 187, 39 (1969)

work page 1969
[37]

H. M. Petrilli and S. Frota -Pessoa, Calculations of the electric field gradient at the nucleus of an Fe impurity in HCP Zr by the recursion method, J. Phys. F: Met. Phys. 15, 2307 (1985)

work page 1985
[38]

Errico, K

L. Errico, K. Lejaeghere, J. Runco, S. N. Mishra, M. Rentería, and S. Cottenier, Precision of Electric-Field Gradient Predictions by Density Functional Theory and Implications for the Nuclear Quadrupole Moment and Its Error Bar of the 111Cd 245 keV 5/2+ Level, J. Phys. Chem. C 120, 23111 (2016)

work page 2016