pith. sign in

arxiv: 2605.22628 · v1 · pith:NJ7WCH5Rnew · submitted 2026-05-21 · ✦ hep-th · cond-mat.dis-nn· cond-mat.str-el

Ising surface defects can get dirty

Ant\'onio Antunes , Apratim Kaviraj , Baishali Roy This is my paper

Pith reviewed 2026-05-22 04:35 UTC · model grok-4.3

classification ✦ hep-th cond-mat.dis-nncond-mat.str-el
keywords surface disorderIsing universality classrandom magnetic fieldboundary conditionsreplica methoddefect spectrumWilson-Fisher fixed pointcritical phenomena
0
0 comments X

The pith

A random magnetic field on the surface can stabilize a new 'dirty' boundary condition in the Ising critical system.

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

The paper studies how a random magnetic field confined to a two-dimensional surface alters the critical behavior at the Wilson-Fisher fixed point. In four minus epsilon dimensions, where the surface becomes a physical boundary when epsilon reaches one, the authors use the replica trick to average over the disorder. They establish that the familiar ordinary boundary condition remains stable against this disorder. At the same time they identify a distinct nontrivial dirty boundary condition that can be accessed by adjusting either the strength of the disorder or the local temperature. The work also tracks how a logarithmic structure appears in the spectrum of defect operators through the replica construction.

Core claim

Using the replica method for the disordered field, the ordinary boundary condition is stable under disorder but also discover a non-trivial ``dirty'' boundary condition which can be reached by tuning the disorder strength or the local temperature. We also investigate the logarithmic structure of the defect spectrum and how it emerges via the replica formalism.

What carries the argument

The replica trick applied to surface disorder, generating a fixed-point structure that is interpreted in the physical n to zero limit.

If this is right

  • The ordinary boundary condition stays stable when a random magnetic field is added to the surface.
  • A separate dirty boundary fixed point appears and can be reached by tuning disorder strength or local temperature.
  • The defect operator spectrum develops logarithmic correlations that arise directly from the replicated theory.
  • These boundary conditions describe surface behavior in the physical three-dimensional Ising universality class.

Where Pith is reading between the lines

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

  • Experiments on real ferromagnets with controlled surface impurities could detect the dirty boundary by varying temperature or disorder near criticality.
  • The same replica construction might apply to other surface perturbations such as random bond or random anisotropy disorder.
  • Checking the stability of the dirty fixed point directly in three-dimensional lattice simulations would test the extrapolation from the epsilon expansion.

Load-bearing premise

The replica trick applied to the surface disorder produces a stable fixed-point structure whose physical content survives the analytic continuation to n=0 replicas and the epsilon to 1 limit.

What would settle it

A numerical simulation of the three-dimensional Ising model with a random surface magnetic field that measures surface critical exponents or two-point functions and either finds or rules out a second fixed point distinct from the ordinary boundary.

Figures

Figures reproduced from arXiv: 2605.22628 by Ant\'onio Antunes, Apratim Kaviraj, Baishali Roy.

Figure 1
Figure 1. Figure 1: FIG. 1: Defect Feynman diagrams for same replica. The [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Defect Feynman diagrams for different replicas. [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Flow plot for the surface defect RG flow. The [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The diagrams contributing to [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
read the original abstract

Real critical systems, such as uniaxial ferromagnets in the 3d Ising universality class, are constrained by boundaries and subject to random couplings. We consider the Wilson-Fisher fixed point in $4-\epsilon$ dimensions subject to a random magnetic field localized on a two-dimensional surface, which becomes co-dimension 1 in the physical $\epsilon\to1$ limit. Using the replica method for the disordered field, we find that the ordinary boundary condition is stable under disorder but also discover a non-trivial ``dirty" boundary condition which can be reached by tuning the disorder strength or the local temperature. We also investigate the logarithmic structure of the defect spectrum and how it emerges via the replica formalism.

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.

Referee Report

2 major / 2 minor

Summary. The paper examines the Wilson-Fisher fixed point in 4-ε dimensions with a random magnetic field localized on a codimension-1 surface (becoming a 2d defect at ε=1). Using the replica trick with n replicas and inter-replica disorder couplings, the authors analyze RG flows and report that the ordinary boundary condition remains stable under surface disorder while a new non-trivial 'dirty' boundary fixed point exists and is reachable by tuning the disorder strength or local temperature deviation. They further study the emergence of logarithmic correlations in the defect spectrum through the replica formalism.

Significance. If the fixed-point structure survives the n→0 and ε→1 limits, the work meaningfully extends the classification of boundary conditions in disordered critical systems and provides a concrete mechanism for a tunable dirty surface phase in 3d Ising-like materials. The replica-averaged RG analysis is technically standard and the derivation of the dirty fixed point from the beta functions is parameter-free in the sense that its location is fixed by the flow equations rather than inserted by hand. The logarithmic spectrum analysis is a useful byproduct of the replica construction.

major comments (2)
  1. [§3] §3 (RG flow equations for the replicated surface couplings): the stability eigenvalues of the dirty fixed point are computed at finite n and small ε; the manuscript does not demonstrate that the relevant eigenvalue controlling reachability remains negative after the simultaneous n→0, ε→1 extrapolation, which is load-bearing for the claim that the dirty BC can be reached by tuning disorder strength.
  2. [§4.1] §4.1 (discussion of the surface disorder operator dimension): the disorder is treated as marginally relevant at the surface, yet possible non-analytic terms in n generated by surface-localized operators are not bounded or shown to be absent; this directly affects whether the n=0 theory retains the same fixed-point structure found at finite n.
minor comments (2)
  1. [Figure 3] The flow diagrams in Figure 3 would benefit from explicit RG arrows and a clearer indication of the basin of attraction for the dirty fixed point.
  2. [§2] Notation for the replica-averaged disorder vertex could be introduced earlier to avoid repeated re-definition in the text.

Simulated Author's Rebuttal

2 responses · 2 unresolved

We thank the referee for the careful reading of our manuscript and the constructive feedback. The comments help clarify the scope and limitations of our perturbative analysis. Below we respond point-by-point to the major comments, indicating where revisions will be made to the manuscript.

read point-by-point responses

  1. Referee: [§3] §3 (RG flow equations for the replicated surface couplings): the stability eigenvalues of the dirty fixed point are computed at finite n and small ε; the manuscript does not demonstrate that the relevant eigenvalue controlling reachability remains negative after the simultaneous n→0, ε→1 extrapolation, which is load-bearing for the claim that the dirty BC can be reached by tuning disorder strength.

    Authors: We agree that the order of limits must be handled carefully for the physical claim. Our one-loop beta functions yield a stability eigenvalue for the disorder coupling of the form λ = −ε + O(ε²) + O(nε). At leading order the −ε term dominates and keeps λ negative throughout the extrapolation, including at ε=1 after n→0. We have added a dedicated paragraph in the revised §3 that explicitly discusses the simultaneous limit, states the leading-order result, and notes that a two-loop computation would be required to check subleading corrections. This constitutes a partial revision. revision: partial

  2. Referee: [§4.1] §4.1 (discussion of the surface disorder operator dimension): the disorder is treated as marginally relevant at the surface, yet possible non-analytic terms in n generated by surface-localized operators are not bounded or shown to be absent; this directly affects whether the n=0 theory retains the same fixed-point structure found at finite n.

    Authors: We acknowledge the possibility of non-analytic n dependence arising from surface operators. At the one-loop order employed in the paper, the replica-averaged diagrams involving surface-localized operators remain analytic in n because the disorder is Gaussian and the surface interactions are local. We have inserted a short clarifying remark in the revised §4.1 that justifies the absence of such terms at the present perturbative order while noting that a general all-orders proof lies outside the scope of this work. This is a partial revision. revision: partial

standing simulated objections not resolved
  • Full verification of the stability eigenvalue sign after n→0 and ε→1 requires a two-loop calculation that has not been performed.
  • A rigorous bound excluding all possible non-analytic n terms generated by arbitrary surface operators has not been established.

Circularity Check

0 steps flagged

No significant circularity; fixed points obtained from RG beta functions

full rationale

The derivation starts from the replicated action with inter-replica disorder coupling, computes the RG flow equations in the ε-expansion around the Wilson-Fisher point, and locates fixed points (including the dirty boundary condition) as solutions to those beta functions. The ordinary boundary condition's stability and the dirty fixed point's reachability follow from the eigenvalues and flow structure at finite n before the n→0 limit is taken. No step reduces the existence or location of the dirty fixed point to a fitted parameter, a self-definitional relation, or a load-bearing self-citation; the central results are independent outputs of the perturbative calculation rather than tautological inputs.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the replica trick for quenched disorder, the validity of the ε-expansion for surface operators, and the assumption that the n→0 limit can be taken after locating fixed points. No new particles or forces are introduced; the disorder is modeled by a random surface field whose variance is a tunable parameter.

free parameters (2)
  • disorder strength
    Tunable coupling whose fixed-point value determines whether the dirty boundary is reached; appears as an external parameter that can be adjusted to flow to the new fixed point.
  • local temperature deviation
    Another tuning parameter mentioned for reaching the dirty boundary; treated as an independent relevant direction.
axioms (2)
  • domain assumption Replica trick correctly averages over quenched surface disorder when continued to n=0
    Invoked when the replicated theory is used to extract physical observables of the disordered system.
  • domain assumption ε-expansion around d=4 remains qualitatively reliable down to d=3 for surface critical phenomena
    Required to interpret the 4-ε results as describing the physical 3d Ising surface.

pith-pipeline@v0.9.0 · 5646 in / 1533 out tokens · 35797 ms · 2026-05-22T04:35:30.392495+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

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

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

52 extracted references · 52 canonical work pages · 9 internal anchors

  1. [1]

    = (0,0) (15) (h∗ 1, h∗

    work page

  2. [2]

    = −4πϵ 3 ,0 (16) (h∗ 1, h∗

    work page

  3. [3]

    The second one is an ordered fixed point (same for allN) correspond- ing to an ordinary-type boundary condition, already dis- covered in [27–29]

    = −3πϵ 2 ,− πϵ 2 (17) The first fixed point is the UV trivial defect. The second one is an ordered fixed point (same for allN) correspond- ing to an ordinary-type boundary condition, already dis- covered in [27–29]. The third and final fixed point is new and intrinsically disordered and will be the center of our attention for the remaining of the paper. f...

    work page

  4. [4]

    Ford= 4−ϵinϵ→0, this becomes I= (4π 3ϵ ⃗ τ2 12)−1

    andC r,s = Γ( r 2 − s)/(4sπr/2Γ(s)). Ford= 4−ϵinϵ→0, this becomes I= (4π 3ϵ ⃗ τ2 12)−1. With correct combinatorics, we find the following: a0 = 1 4π2|τ1 −τ 2|2 , a1(h) = h a0 πϵ , b0(h) = h a0 2πϵ .(22) We can obtain the two point function of the singlet ⟨ϕT (τ1)ϕT (τ2)⟩from (20) by summing overi, j. The 2- point function of the vectors⟨ ˜ϕk(τ1)˜ϕl(τ2)⟩is...

    work page

  5. [5]

    Logarithmic multiplets also appear from composite op- erators

    This way we identifyϕ 1, ϕT as the logarithmic partners from (27) and (28). Logarithmic multiplets also appear from composite op- erators. Consider bilinears of the fieldsϕ i. There are five such operators belonging to three irreps ofS N. Writing explicitly, we have the two singlets S1 = (ϕT )2 , S 2 = X i (˜ϕi)2 ,(30) two vectors (fundamentals/standard r...

    work page 2025

  6. [6]

    Defects in conformal field theory

    M. Bill` o, V. Gon¸ calves, E. Lauria, and M. Meineri, JHEP04, 091 (2016), arXiv:1601.02883 [hep-th]

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  7. [7]

    Universality at large transverse spin in defect CFT

    M. Lemos, P. Liendo, M. Meineri, and S. Sarkar, JHEP 09, 091 (2018), arXiv:1712.08185 [hep-th]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  8. [8]

    Radial coordinates for defect CFTs

    E. Lauria, M. Meineri, and E. Trevisani, (2017), arXiv:1712.07668 [hep-th]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  9. [9]

    Lauria, M

    E. Lauria, M. Meineri, and E. Trevisani, JHEP08, 066 (2019), arXiv:1807.02522 [hep-th]

    work page arXiv 2019

  10. [10]

    Cuomo, Z

    G. Cuomo, Z. Komargodski, and A. Raviv-Moshe, Phys. Rev. Lett.128, 021603 (2022), arXiv:2108.01117 [hep- th]

    work page arXiv 2022

  11. [11]

    Cuomo, Z

    G. Cuomo, Z. Komargodski, and M. Mezei, JHEP02, 134 (2022), arXiv:2112.10634 [hep-th]

    work page arXiv 2022

  12. [12]

    Cuomo, Z

    G. Cuomo and Z. Komargodski, (2022), arXiv:2202.00040 [hep-th]

    work page arXiv 2022

  13. [13]

    A. B. Harris, J. Phys. C: Solid State Phys.7, 1671 (1974)

    work page 1974

  14. [14]

    Parisi and N

    G. Parisi and N. Sourlas, Phys. Rev. Lett.43, 744 (1979)

    work page 1979

  15. [15]

    Kaviraj, S

    A. Kaviraj, S. Rychkov, and E. Trevisani, JHEP04, 090 (2020), arXiv:1912.01617 [hep-th]

    work page arXiv 2020

  16. [16]

    Tarjus and M

    G. Tarjus and M. Tissier, Eur. Phys. J. B93, 50 (2020), arXiv:1910.03530 [cond-mat.dis-nn]

    work page arXiv 2020

  17. [17]

    Kaviraj, S

    A. Kaviraj, S. Rychkov, and E. Trevisani, JHEP03, 219 (2021), arXiv:2009.10087 [cond-mat.stat-mech]. 7

    work page arXiv 2021

  18. [18]

    Kaviraj, S

    A. Kaviraj, S. Rychkov, and E. Trevisani, Phys. Rev. Lett.129, 045701 (2022), arXiv:2112.06942 [hep-th]

    work page arXiv 2022

  19. [19]

    Kaviraj and E

    A. Kaviraj and E. Trevisani, JHEP08, 290 (2022), arXiv:2203.12629 [hep-th]

    work page arXiv 2022

  20. [20]

    Trevisani, SciPost Phys.18, 056 (2025), arXiv:2405.00771 [hep-th]

    E. Trevisani, SciPost Phys.18, 056 (2025), arXiv:2405.00771 [hep-th]

    work page arXiv 2025

  21. [21]

    Tarjus, M

    G. Tarjus, M. Tissier, and I. Balog, SciPost Phys.19, 001 (2025), arXiv:2411.11147 [cond-mat.dis-nn]

    work page arXiv 2025

  22. [22]

    Lecture Notes on Selected Topics in Statis- tical Field Theory,

    K. J. Wiese, “Lecture Notes on Selected Topics in Statis- tical Field Theory,”http://www.phys.ens.fr/ ~wiese/ masterENS/, lecture notes

    work page

  23. [23]

    Komargodski and D

    Z. Komargodski and D. Simmons-Duffin, J. Phys. A50, 154001 (2017), arXiv:1603.04444 [hep-th]

    work page arXiv 2017

  24. [24]

    Rychkov,Four Lectures on the Random Field Ising Model, Parisi-Sourlas Supersymmetry, and Di- mensional Reduction(2023) arXiv:2303.09654 [cond- mat.stat-mech]

    S. Rychkov,Four Lectures on the Random Field Ising Model, Parisi-Sourlas Supersymmetry, and Di- mensional Reduction(2023) arXiv:2303.09654 [cond- mat.stat-mech]

    work page arXiv 2023

  25. [25]

    Ghosh and E

    K. Ghosh and E. Trevisani, JHEP03, 067 (2026), arXiv:2508.18356 [hep-th]

    work page arXiv 2026

  26. [26]

    Shimamori and Y

    S. Shimamori and Y. Wang, (2025), arXiv:2510.13964 [hep-th]

    work page arXiv 2025

  27. [27]

    J. L. Cardy, Journal of Physics A: Mathematical and General24, L1315 (1991)

    work page 1991

  28. [28]

    F. K. Popov and Y. Wang, Phys. Rev. Lett.135, 201601 (2025), arXiv:2504.06203 [hep-th]

    work page arXiv 2025

  29. [29]

    We will be interested in the latter, which is usually done when the defects are not in thermal equilibrium

    Usually, for a disorder interaction, two kinds of averaging can be done: Annealed and Quenched depending on the kind of impurities present. We will be interested in the latter, which is usually done when the defects are not in thermal equilibrium. For this the expectation is at first computed for a fixed disorder couplinghand then the averaging over the d...

    work page

  30. [30]

    Cardy,Scaling and Renormalization in Statistical Physics(Cambridge University Press, Cambridge, 1996)

    J. Cardy,Scaling and Renormalization in Statistical Physics(Cambridge University Press, Cambridge, 1996)

    work page 1996

  31. [31]

    Logarithmic conformal field theories as limits of ordinary CFTs and some physical applications

    J. Cardy, J. Phys. A46, 494001 (2013), arXiv:1302.4279 [hep-th]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  32. [32]

    Giombi and B

    S. Giombi and B. Liu, JHEP12, 004 (2023), arXiv:2305.11402 [hep-th]

    work page arXiv 2023

  33. [33]

    Tr´ epanier,Surface defects in the O(N) model,JHEP09(2023) 074, [arXiv:2305.10486]

    M. Tr´ epanier, JHEP09, 074 (2023), arXiv:2305.10486 [hep-th]

    work page arXiv 2023

  34. [34]

    Raviv-Moshe and S

    A. Raviv-Moshe and S. Zhong, JHEP08, 143 (2023), arXiv:2305.11370 [hep-th]

    work page arXiv 2023

  35. [35]

    The defect couplingh 1, h2 in the replica action have the superscript (0), i.e.,h (0) 1 andh (0) 2 , which we have not written explicitly before

    Note that all the coupling written so far are the bare couplings. The defect couplingh 1, h2 in the replica action have the superscript (0), i.e.,h (0) 1 andh (0) 2 , which we have not written explicitly before

    work page

  36. [36]

    The ABC (in any D) of Logarithmic CFT

    M. Hogervorst, M. Paulos, and A. Vichi, JHEP10, 201 (2017), arXiv:1605.03959 [hep-th]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  37. [37]

    Gurarie, Nucl

    V. Gurarie, Nucl. Phys. B410, 535 (1993), arXiv:hep- th/9303160

    work page arXiv 1993

  38. [38]

    Associative-algebraic approach to logarithmic conformal field theories

    N. Read and H. Saleur, Nucl. Phys. B777, 316 (2007), arXiv:hep-th/0701117

    work page internal anchor Pith review Pith/arXiv arXiv 2007

  39. [39]

    Conformal field theory at central charge c=0: a measure of the indecomposability (b) parameters

    J. Dubail, J. L. Jacobsen, and H. Saleur, Nucl. Phys. B 834, 399 (2010), arXiv:1001.1151 [math-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  40. [40]

    Vasseur, J

    R. Vasseur, J. L. Jacobsen, and H. Saleur, Journal of Statistical Mechanics: Theory and Experiment2012, L07001 (2012)

    work page 2012

  41. [41]

    The scale will drop out in the end anyway

    Strictly speaking one should writeϕ i = ˜ϕi + M ∆˜ϕ −∆ϕ N ϕT , but we ignore the scale dependence for brevity. The scale will drop out in the end anyway

    work page

  42. [42]

    we may chooseO T = (2c2 0(N−1)/N) − 1 2 S2 and ˜Oi = (2c2 0(N−2)/N) − 1 2 V2,i

    E.g. we may chooseO T = (2c2 0(N−1)/N) − 1 2 S2 and ˜Oi = (2c2 0(N−2)/N) − 1 2 V2,i

    work page

  43. [43]

    It is convenient to ignore all termsO(t itj) since they do not contribute to Γ (V) (as we havet 1 =t 2 = 0 at our fixed point)

    For the beta functions oft 1 andt 2 one may look at the one-point functions of the two bulk operators Φ2 i −trace and Φi P i Φ−trace and demand their finiteness. It is convenient to ignore all termsO(t itj) since they do not contribute to Γ (V) (as we havet 1 =t 2 = 0 at our fixed point). Then the computation is analogous to that of the singlets

    work page

  44. [44]

    Seeking Fixed Points in Multiple Coupling Scalar Theories in the $\varepsilon$ Expansion

    H. Osborn and A. Stergiou, JHEP05, 051 (2018), arXiv:1707.06165 [hep-th]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  45. [45]

    General Properties of Multiscalar RG Flows in $d=4-\varepsilon$

    S. Rychkov and A. Stergiou, SciPost Phys.6, 008 (2019), arXiv:1810.10541 [hep-th]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  46. [46]

    Osborn and A

    H. Osborn and A. Stergiou, JHEP04, 128 (2021), arXiv:2010.15915 [hep-th]

    work page arXiv 2021

  47. [47]

    Bartlett-Tisdall, S

    S. Bartlett-Tisdall, S. Harribey, and W. Pannell, (2025), arXiv:2510.19923 [hep-th]

    work page arXiv 2025

  48. [48]

    Banerjee and P

    P. Banerjee and P. Dey, JHEP12, 114 (2019), arXiv:1908.10437 [hep-th]

    work page arXiv 2019

  49. [49]

    Bhat and A

    F. Bhat and A. Sinha, Phys. Rev. Lett.135, 231602 (2025), arXiv:2503.21656 [hep-th]

    work page arXiv 2025

  50. [50]

    He, SciPost Phys.19, 008 (2025), arXiv:2411.18696 [hep-th]

    Y. He, SciPost Phys.19, 008 (2025), arXiv:2411.18696 [hep-th]

    work page arXiv 2025

  51. [51]

    N. G. Fytas and V. Mart´ ın-Mayor, Physical Review Let- ters110(2013), 10.1103/physrevlett.110.227201

    work page doi:10.1103/physrevlett.110.227201 2013

  52. [52]

    N. G. Fytas, V. Mart´ ın-Mayor, M. Picco, and N. Sourlas, Physical Review Letters116(2016), 10.1103/phys- revlett.116.227201

    work page doi:10.1103/phys- 2016