ABSTRACT (original article): The site-decorated Ising model is introduced to advance the understanding and experimental realization of the recently discovered one-dimensional finite-temperature ultranarrow phase crossover in an external magnetic field, while mitigating the geometric complexities of traditional bond-decorated models. Furthermore, although higher-dimensional Ising models in an external field remain unsolved, an exact solution for a novel spin-reversal transition -- driven by an exotic, hidden “half-ice, half-fire’’ state induced by site decoration -- is derived. This transition, triggered by a slight variation in temperature or magnetic field even in the weak-field limit, offers a promising route toward energy-efficient applications such as data storage and processing. The results establish site decoration as a compelling new avenue for materials and device design, particularly in systems such as mixed - compounds, optical lattices, and neural networks. CITATION (original article): Weiguo Yin (2025), Site-Decorated Model for Unconventional Frustrated Magnets: Ultranarrow Phase Crossover and Spin Reversal Transition, arXiv:2502.11270. https://doi.org/10.48550/arXiv.2502.11270
Related: for context see previous works
◼
v1 .0 Single - chain bond - decorated Ising model (M = 1)
Weiguo Yin and A . M . Tsvelik, Phase Switch Driven by the Hidden Half - Ice, Half - Fire State in a Ferrimagnet, Phys . Rev . Lett . 133, 266701 (2024) https://community.wolfram.com/groups/-/m/t/3530595
Weiguo Yin and A . M . Tsvelik, Phase Switch Driven by the Hidden Half - Ice, Half - Fire State in a Ferrimagnet, Phys . Rev . Lett . 133, 266701 (2024) https://community.wolfram.com/groups/-/m/t/3530595
◼
v1 .1 2025 - 01 - 22 Magnetocaloric effect
◼
v2 .0 2025 - 02 - 05 Site - decorated Ising model (M = 2) submitted to PRB .
◼
v2 .1 2025 - 06 - 22 The hybrid - exact - mean - field theory
The unconventional frustration and novel physics are clarified by exactly mapping the one-dimensional site-decorated Ising model in a magnetic field onto a zero-field bond-decorated J1–J2 Ising model with conventional geometrical frustration. Furthermore, although higher-dimensional Ising models in an external field remain unsolved exactly, an exact solution for a novel spin-reversal transition—driven by an exotic, hidden ‘half-ice, half-fire’ state induced by site decoration—is derived. This transition, triggered by a slight variation in temperature or magnetic field—without changing its direction—even in the weak-field limit, offers a promising route toward energy-efficient applications such as data storage and processing. The results establish site decoration as a compelling new avenue for materials and device design, particularly in systems such as mixed d-f compounds, optical lattices, and neural networks. In addition, the site-decorated model offers a rigorous test ground for AI in science, as the analytic derivation of the present results was not only validated but also improved by a general-purpose large language model, inspiring the use of AI as scientific discoverer.
Set up output directory, color, precision
Set up output directory, color, precision
I. Methods
I. Methods
Step 1 : Set up the 4 states of the two backbone a - spins on a rung for the ladder or on a bond for the single chain.
Don’ t change this part
Don’ t change this part
In[]:=
Clear[S,outer1,outer2,states,parent,nstatesp]S=1/2;outer1=1;outer2=2;(*states=Array[{#1-(S+1),#2-(S+1),#3-(S+1)}&,{sdim,sdim,sdim}];*)states[S_,M_]:=Module[{sdim,matDim,tmp},sdim=2*S+1;matDim=ConstantArray[sdim,M];tmp=Array[{##}&,matDim]-(sdim-S);tmp=Reverse[Flatten[tmp,M-1]];tmp=If[OddQ[2S],2*tmp,If[EvenQ[2S],tmp]]]parent=states[S,2];nstatesp=Length[parent];Print["parent states= ",parent//MatrixForm]Print[nstatesp," states in total"]
parent states=
1 | 1 |
1 | -1 |
-1 | 1 |
-1 | -1 |
4 states in total
Step 2 : Set up the states for the M decorating b - spins : M = 1 for bond decoration; M = 2 for site decoration, each attached to one end backbone spin of a bond
In[]:=
Clear[Sc,M,child,nstatesc]M=2;Sc=1/2;Print["M= ",M]child=states[Sc,M];nstatesc=Length[child];Print["decorated states= ",child//MatrixForm]Print[nstatesc," states in total"]
M= 2
decorated states=
1 | 1 |
1 | -1 |
-1 | 1 |
-1 | -1 |
4 states in total
Step 3 : Sum over the decorating b - spins
In[]:=
Clear[childAllConf,childOneConf,gMat,yMat,muMat,J1,J2,a,b,h,pp,mm,pm,mp,SiteBond]SiteBond=1/2;(*1forbonddecoration,1/2forsitedecoration*)J2=0;gMat=Table[If[m1≠m2,J2,0],{m1,1,M},{m2,1,M}];(*yMat=Table[ConstantArray[J1,M],{m,1,2}];*)yMat={{J1,0},{0,J1}};muMat=ConstantArray[b,M];Print["yMat= ",yMat//MatrixForm]Print["gMat= ",gMat//MatrixForm]Print["muMat= ",muMat//MatrixForm]childOneConf[k_,i_]:=Module[{sumM,sumMM,sumSite},(*k_=1...4:parent'sspin1&2configuration,i_=1...2^M:children'sspinconfiguration*)sumM=0;sumMM=0;Do[sumM=sumM+child[[i,m1]]*(yMat[[outer1,m1]]*parent[[k,outer1]]+yMat[[outer2,m1]]*parent[[k,outer2]]+hmuMat[[m1]]),{m1,1,M}];Do[sumMM=sumMM+gMat[[m1,m2]]*child[[i,m1]]*child[[i,m2]],{m1,1,M},{m2,1,M}];Exp[sumM+sumMM/2]]childAllConf[k_]:=Module[{sum,Mdim},If[M==0,Return[1]];Mdim=nstatesc;sum=0;Do[sum=sum+childOneConf[k,i],{i,Mdim}];sum]pp=childAllConf[1]^SiteBond;pm=childAllConf[2]^SiteBond;mp=childAllConf[3]^SiteBond;mm=childAllConf[4]^SiteBond;pp=PowerExpand[FullSimplify[pp]];pm=PowerExpand[FullSimplify[pm]];mp=PowerExpand[FullSimplify[mp]];mm=PowerExpand[FullSimplify[mm]];Print["The decoration's contribution= ",{{pp,pm},{mp,mm}}//MatrixForm]Print[">>>>>Checkpoint 1: pm==mp? ",pmmp]Print[">>>>>Checkpoint 2: pp*mm==pm*mp? ",Simplify[pp*mm]==Simplify[pm*mp]]Clear[Ym,Yp,aa,bb,cc,transferM]Ym=Exp[ah]pp-Exp[-ah]mm;Yp=Exp[ah]pp+Exp[-ah]mm;Print["The frustration function= ",Ym]aa=Exp[J+ah]pp;bb=Exp[J-ah]mm;cc=Exp[-J]pm;transferM={{aa,cc},{cc,bb}};Print["transfer Matrix= ",transferM//MatrixForm]
yMat=
J1 | 0 |
0 | J1 |
gMat=
0 | 0 |
0 | 0 |
muMat=
b |
b |
The decoration's contribution=
2Cosh[bh+J1] | 2Cosh[2bh]+2Cosh[2J1] |
2Cosh[2bh]+2Cosh[2J1] | 2Cosh[bh-J1] |
>>>>>Checkpoint 1: pm==mp? True
>>>>>Checkpoint 2: pp*mm==pm*mp? True
The frustration function= -2Cosh[bh-J1]+2Cosh[bh+J1]
-ah
ah
transfer Matrix=
2 ah+J | -J 2Cosh[2bh]+2Cosh[2J1] |
-J 2Cosh[2bh]+2Cosh[2J1] | 2 -ah+J |
Step 4 : Solve the transfer Matrix using Method I
In[]:=
Clear[λ,paraT,free,entropy,Cv,magnetization,Chi,Caa,Cbb,ThisTc,OP,dT,maa,mbb,C12,C13,t,J,J1,J2,h,a,b]paraT[t_]:={JJ/t,J1J1/t,J2J2/t,hh/t};λ=1/2Exp[J](Yp+Sqrt[Ym^2+4Exp[-4J]pm^2]);free[t_]:=-tLog[λ]/.paraT[t]Print["Free energy free[t]= ",free[t]]entropy[t_]:=-D[free[ttt],ttt]/.{tttt}entropyDiff[t_,hh_]:=(entropy[t]/.{h->hh})-(entropy[t]/.{h->0})Cv[t_]:=tD[entropy[ttt],ttt]/.{tttt}magnetization[t_]:=-D[free[t],h]Chi[t_]:=D[magnetization[t],h]Chih0[t_]:=Exp[2J/t]/2(Exp[2J1/t](a+b)^2+Exp[-2J1/t](a-b)^2)/t-Sinh[2J/t]b^2/(2Cosh[2J1/t]+2Exp[-2J/t])/tCaa[t_]:=-D[free[t],a]Cbb[t_]:=-D[free[t],b]maa[t_]:=-D[free[t],a]/hmbb[t_]:=-D[free[t],b]/hC12[t_]:=-D[free[t],J]C13[t_]:=-D[free[t],J1]/2ThisTc[t_]:=Ym/.paraT[t];Print["ThisTc[t]= ",ThisTc[t]]OP[t_]:=Ym/Sqrt[Ym^2+4Exp[-4J]pm^2]/.paraT[t];dT[t_]:=1/D[OP[ttt],ttt]/.{tttt}Print["order parameter OP[t]= ",OP[t]]
Step 5 : Solve the transfer Matrix using Method II - the effective interaction method
Test low - T heff (T = 0.001)
Step 6 : Use the two methods to calculate T0 and dT
Use two methods to calculate T0 and dT . Confirmed
Test GetTc[r, h], GetdT[J, r, h] for h = 0. Passed
A. More tests
A. More tests
Test two methods : OP vs OPeff
Comparison plot testing equivalence between the direct order parameter (OP) and effective field order parameter (OPeff) for backbone type-a spins at multiple magnetic field values, verifying both computational methods yield identical results.
Test two methods : maa vs maaeff
Comparison plot verifying that the normalized magnetization of type-a spins computed via both methods (maa vs maaeff) produces consistent results across different field strengths.
Test two methods : mbb vs mbbeff
Comparison testing the magnetization of decorated type-b spins calculated using direct summation (mbb) versus the effective field approach (mbbeff), confirming consistency between methods.
Test two methods : mbbeff vs mbbeff2
Validation comparing mbbeff with an alternative analytical expression (mbbeff2) derived from the effective field and order parameter, demonstrating mathematical equivalence for site-decorated models.
II. Results
II. Results
A. One dimension
A. One dimension
1. Ultranarrow phase crossover
1. Ultranarrow phase crossover
Fig. 2 (a) : heff
Fig. 2 (b) : Definition of dT. The slope is verified.
[Placeholder 7]
Fig. 3 (c): entropy for J = 0
Fig. 3 (f): entropy for J = 20
Fig. 3 (f): entropy for J = 20
2 The weak h limit
2 The weak h limit
Fig. 4 (a) Weak h limit: T0 as a function of μa/μb.
Fig. 3 inset
Fig. 4 (d) Weak h limit: magnetization
B. Higher dimensions
B. Higher dimensions
1. Exact results
1. Exact results
Functions for magnetization at h=0 and heff=0 in a square lattice
Defines functions for computing exact spontaneous magnetization at the effective field condition (hₑff = 0) in 2D square lattices, enabling analysis of spin-reversal transitions in higher-dimensional systems.
Plots magnetization of decorated type-b spins in 2D at the spin-reversal transition point, showing distinct behavior for different backbone coupling strengths relative to the critical temperature.
Fig. 5(c): magnetization
Computes and visualizes total system magnetization (combining both type-a and type-b contributions) at the spin-reversal transition, demonstrating the magnetocaloric effect potential for applications.
2. Hybrid-exact-mean-field theory
2. Hybrid-exact-mean-field theory
Implements a hybrid theoretical framework combining exact 1D transfer matrix solutions with mean-field treatment of inter-chain coupling for 2D systems, enabling prediction of phase behavior beyond purely exact or purely mean-field approaches.
Parameters
Theory
0. Exact zero-field magnetization in a 2D square lattice
1. Purely exact 1 D (transfer - matrix) solution
3. Hybrid exact/mean - field
CITE THIS NOTEBOOK
CITE THIS NOTEBOOK
Site-decorated model for unconventional frustrated magnets: ultranarrow phase crossover and spin reversal transition
by Weiguo Yin
Wolfram Community, STAFF PICKS, December 4, 2025
https://community.wolfram.com/groups/-/m/t/3586246
by Weiguo Yin
Wolfram Community, STAFF PICKS, December 4, 2025
https://community.wolfram.com/groups/-/m/t/3586246