Skip to content

Commit e433f43

Browse files
cajoharecajohare
committed
updates for XENON, Fermi-Supernovae and solar axion basin
1 parent 6946bf7 commit e433f43

16 files changed

+423
-3465
lines changed

‎AxionElectron.ipynb‎

Lines changed: 1 addition & 9 deletions
Large diffs are not rendered by default.

‎AxionPhoton.ipynb‎

Lines changed: 3 additions & 10 deletions
Large diffs are not rendered by default.

‎PlotFuncs.py‎

Lines changed: 41 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -630,19 +630,19 @@ def AstroBounds(ax,projection=True,fs=15):
630630
y2 = ax.get_ylim()[1]
631631
### Astrophysical constraints
632632

633-
# SN-gamma rays arXiv:[1410.3747]
633+
# Fermi extragalactic SN gamma rays arXiv:[2006.06722]
634634
SNgamma_col = [0.05, 0.5, 0.06]
635635
SNgamma = loadtxt("limit_data/AxionPhoton/SN-gamma.txt")
636636
plt.plot(SNgamma[:,0],SNgamma[:,1],'k-',alpha=0.6,zorder=0.21,lw=2)
637637
plt.fill_between(SNgamma[:,0],SNgamma[:,1],y2=y2,edgecolor=None,facecolor=SNgamma_col,zorder=0.21)
638638
plt.text(3e-11,2e-11,r'{\bf SN}-$\gamma$',fontsize=fs,color='w',ha='left',va='top')
639639

640640
# M87 Limits from arXiv:[1703.07354]
641-
M87_col = [0.0, 0.66, 0.42]
641+
M87_col = 'seagreen'
642642
M87 = loadtxt("limit_data/AxionPhoton/M87.txt")
643-
plt.plot(M87[:,0],M87[:,1],'k-',lw=2,alpha=0.8,zorder=0.2)
643+
plt.plot(M87[:,0],M87[:,1],'k-',lw=2,alpha=1,zorder=0.2)
644644
plt.fill_between(M87[:,0],M87[:,1],y2=y2,edgecolor=None,facecolor=M87_col,zorder=0.2)
645-
plt.text(1.4e-12,4e-12,r'\quad {\bf M87}',fontsize=fs,color='w',ha='left',va='top')
645+
plt.text(1.4e-12,5e-12,r'\quad {\bf M87}',fontsize=fs,color='w',ha='left',va='top')
646646

647647
# HYDRA-A arXiv:[1304.0989]
648648
HYDRA_col = [0.24, 0.71, 0.54]
@@ -855,25 +855,43 @@ def g_x(C_ae,m_a):
855855

856856
return
857857

858-
def XENON1T(ax,col='m',fs=20):
859-
# XENON1T LDM Searches arXiv:[1907.11485]
858+
def XENON1T(ax,col='crimson',fs=20):
859+
# XENON1T S2 analysis arXiv:[1907.11485]
860860
y2 = ax.get_ylim()[1]
861-
dat = loadtxt("limit_data/AxionElectron/XENON1T.txt")
861+
dat = loadtxt("limit_data/AxionElectron/XENON1T_DM_S2.txt")
862862
plt.plot(dat[:,0],dat[:,1],'k-',alpha=0.6,zorder=0.51,lw=2)
863863
plt.fill_between(dat[:,0],dat[:,1],y2=y2,edgecolor=None,facecolor=col,zorder=0.51)
864-
plt.text(0.25e2,4e-14,r'{\bf XENON1T}',fontsize=fs,color=col,ha='left',va='top')
864+
865+
# XENON1T S1+S2 analysis arXiv:[2006.09721]
866+
dat = loadtxt("limit_data/AxionElectron/XENON1T_DM_S1S2.txt")
867+
plt.plot(dat[:,0],dat[:,1],'k-',alpha=0.6,zorder=0.51,lw=2)
868+
plt.fill_between(dat[:,0],dat[:,1],y2=y2,edgecolor=None,facecolor=col,zorder=0.51)
869+
870+
plt.text(1.2e2,4e-14,r'{\bf XENON1T}',fontsize=fs,color=col,ha='center',va='top')
871+
plt.text(1.2e2,2.5e-14,r'(DM)',fontsize=fs,color=col,ha='center',va='top')
872+
873+
# Solar axion basin arXiv:[2006.12431]
874+
col = 'royalblue'
875+
y2 = ax.get_ylim()[1]
876+
dat = loadtxt("limit_data/AxionElectron/XENON1T_S2_SolarAxionBasin.txt")
877+
plt.plot(dat[:,0],dat[:,1],'k-',alpha=0.6,zorder=0.6,lw=2)
878+
plt.fill_between(dat[:,0],dat[:,1],y2=y2,edgecolor=None,facecolor=col,zorder=0.6)
879+
plt.text(1.5e3,2e-11,r'{\bf XENON1T}',fontsize=fs,color='k',ha='center',va='top')
880+
plt.text(1.5e3,1.3e-11,r'(Solar axion',fontsize=fs,color='k',ha='center',va='top')
881+
plt.text(1.5e3,0.8e-11,r' basin)',fontsize=fs,color='k',ha='center',va='top')
882+
865883
return
866884

867-
def LUX(ax,col='darkorchid',fs=20):
885+
def LUX(ax,col='darkred',fs=20):
868886
# LUX arXiv:[1704.02297]
869887
y2 = ax.get_ylim()[1]
870888
dat = loadtxt("limit_data/AxionElectron/LUX.txt")
871889
plt.plot(dat[:,0],dat[:,1],'k-',alpha=0.6,zorder=0.52,lw=2)
872890
plt.fill_between(dat[:,0],dat[:,1],y2=y2,edgecolor=None,facecolor=col,zorder=0.52)
873-
plt.text(0.5e-8,6e-12,r'{\bf LUX}',fontsize=30,color='w',alpha=0.8,ha='left',va='top')
891+
plt.text(0.5e-8,6e-12,r'{\bf LUX} (Solar axions)',fontsize=30,color='w',alpha=0.8,ha='left',va='top')
874892
return
875893

876-
def PandaX(ax,col='mediumvioletred',fs=20):
894+
def PandaX(ax,col='firebrick',fs=20):
877895
# PandaX arXiv:[1707.07921]
878896
y2 = ax.get_ylim()[1]
879897
# Currently not using Solar pandaX limit
@@ -886,7 +904,7 @@ def PandaX(ax,col='mediumvioletred',fs=20):
886904
plt.text(1.2e3,4.5e-13,r'{\bf PandaX}',fontsize=fs-2,color='w',ha='left',va='top',rotation=20)
887905
return
888906

889-
def EDELWEISS(ax,col='crimson',projection=True,fs=20):
907+
def EDELWEISS(ax,col='darkred',projection=True,fs=20):
890908
# EDELWEISS arXiv:[1808.02340]
891909
y2 = ax.get_ylim()[1]
892910
dat = loadtxt("limit_data/AxionElectron/EDELWEISS.txt")
@@ -898,17 +916,17 @@ def EDELWEISS(ax,col='crimson',projection=True,fs=20):
898916
plt.text(9e0,7e-13,r'{\bf EDELWEISS}',fontsize=fs,color=col,ha='left',va='top')
899917
return
900918

901-
def SuperCDMS(ax,col='orchid',fs=20):
919+
def SuperCDMS(ax,col='maroon',fs=20):
902920
# SuperCDMS arXiv:[1911.11905]
903921
y2 = ax.get_ylim()[1]
904922
dat = loadtxt("limit_data/AxionElectron/SuperCDMS.txt")
905923
plt.fill_between(dat[:,0],dat[:,1],y2=y2,edgecolor=None,facecolor=col,zorder=0.58)
906924
plt.plot(dat[:,0],dat[:,1],'-',color='k',alpha=0.5,zorder=0.58,lw=3)
907-
plt.text(1e2,9e-12,r'{\bf SuperCDMS}',fontsize=fs,color='k',ha='left',va='top',alpha=0.8)
925+
plt.text(2e1,2.7e-11,r'{\bf SuperCDMS}',fontsize=fs-1,color='w',ha='left',va='top',alpha=0.8,rotation=-82)
908926
return
909927

910-
def DARWIN(ax,col='blueviolet',fs=20):
911-
# PandaX arXiv:[1606.07001]
928+
def DARWIN(ax,col='brown',fs=20):
929+
# DARWIN arXiv:[1606.07001]
912930
y2 = ax.get_ylim()[1]
913931
dat = loadtxt("limit_data/AxionElectron/Projections/DARWIN.txt")
914932
plt.plot(dat[:,0],dat[:,1],'--',color=col,alpha=1.0,zorder=0.1,lw=3)
@@ -927,17 +945,19 @@ def Magnon(ax,col='rebeccapurple',fs=20):
927945
# Axion-magnon conversion arXiv:[2005.10256]
928946
y2 = ax.get_ylim()[1]
929947
dat = loadtxt("limit_data/AxionElectron/Projections/Magnon.txt")
930-
plt.plot(dat[:,0],dat[:,1],'-',color=col,alpha=1.0,zorder=0.5,lw=3)
948+
plt.fill_between(dat[:,0],dat[:,1],y2=1,color=col,alpha=0.4,zorder=0.51)
949+
plt.plot(dat[:,0],dat[:,1],'--',color=col,alpha=1.0,zorder=0.5,lw=3)
931950
plt.text(2e-6,1e-14,r'{\bf Magnons \newline (YIT, NiSP$_3$)}',fontsize=fs,color=col,ha='left',va='top',rotation=0)
932951
return
933952

934953
def MagnonScan(ax,col='mediumvioletred',fs=20):
935954
# Axion-magnon conversion arXiv:[2005.10256 and 2001.10666]
936955
y2 = ax.get_ylim()[1]
937956
dat = loadtxt("limit_data/AxionElectron/Projections/MagnonScan.txt")
938-
plt.plot(dat[:,0],dat[:,1],'-',color=col,alpha=1.0,zorder=0.5,lw=3)
939-
plt.text(1.8e-5,1e-13,r'{\bf Magnons}',fontsize=fs-1,color=col,ha='center',va='top',rotation=0)
940-
plt.text(1.8e-5,0.7*1e-13,r'{\bf (Scanning)}',fontsize=fs-1,color=col,ha='center',va='top',rotation=0)
957+
plt.fill_between(dat[:,0],dat[:,1],y2=1,color=col,alpha=0.4,zorder=0.51)
958+
plt.plot(dat[:,0],dat[:,1],'--',color=col,alpha=1.0,zorder=0.5,lw=3)
959+
plt.text(1.8e-5,0.5e-13,r'{\bf Magnons}',fontsize=fs-1,color=col,ha='center',va='top',rotation=0)
960+
plt.text(1.8e-5,0.7*0.5e-13,r'{\bf (Scanning)}',fontsize=fs-1,color=col,ha='center',va='top',rotation=0)
941961

942962
return
943963

@@ -980,7 +1000,7 @@ def StellarBounds(ax,fs=30,Hint=True):
9801000
plt.text(0.5e-8,3.5e-13,r'{\bf Red giants}',fontsize=fs,color='w')
9811001

9821002
# Solar neutrinos arXiv:[0807.2926]
983-
SolarNu_col = [0.01, 0.75, 0.24]
1003+
SolarNu_col = 'seagreen'
9841004
SolarNu = loadtxt("limit_data/AxionElectron/SolarNu.txt")
9851005
plt.plot(SolarNu[:,0],SolarNu[:,1],color='k',lw=2,alpha=0.7,zorder=1)
9861006
plt.fill_between(SolarNu[:,0],SolarNu[:,1],y2=y2,edgecolor=None,facecolor=SolarNu_col,zorder=0.7)

‎README.md‎

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -66,7 +66,7 @@ Download: [png](https://github.com/cajohare/AxionLimits/raw/master/plots/plots_p
6666
* HESS (PKS 2155-304): [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionPhoton/HESS.txt), [reference](https://arxiv.org/abs/1311.3148)
6767
* Horizontal branch stars: [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionPhoton/HorizontalBranch.txt), [reference](https://arxiv.org/abs/1406.6053)
6868
* Solar neutrinos: [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionPhoton/SolarNu.txt), [reference](https://arxiv.org/abs/1501.01639)
69-
* SN1987A gamma rays: [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionPhoton/SN-gamma.txt), [reference](https://arxiv.org/abs/1410.3747)
69+
* Extragalacic supernovae gamma rays (Fermi): [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionPhoton/SN-gamma.txt), [reference](https://arxiv.org/abs/2006.06722)
7070
* SN1987A ALP decay: [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionPhoton/SN1986A_decay.txt), [reference](https://arxiv.org/abs/1702.02964)
7171
* SN1987A neutrinos: [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionPhoton/SN1986A_neutrinos.txt), [reference](https://arxiv.org/abs/1808.10136)
7272
* VLT (Abell 2667 and 2390): [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionPhoton/Telescopes.txt), [reference](https://arxiv.org/abs/astro-ph/0611502)
@@ -88,7 +88,9 @@ Download: [png](https://github.com/cajohare/AxionLimits/raw/master/plots/plots_p
8888
* PandaX-II: [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionElectron/PandaX.txt), [reference](https://arxiv.org/abs/1707.07921)
8989
* Semiconductors (absorption): [projection](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionElectron/Projections/SemiconductorAbsorption.txt), [reference](https://arxiv.org/abs/1608.02123)
9090
* SuperCDMS: [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionElectron/SuperCDMS.txt), [reference](https://arxiv.org/abs/1911.11905)
91-
* XENON1T: [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionElectron/XENON1T.txt), [reference](https://arxiv.org/abs/1907.11485)
91+
* XENON1T (ALP DM search, S2 only): [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionElectron/XENON1T_DM_S2.txt), [reference](https://arxiv.org/abs/1907.11485)
92+
* XENON1T (ALP DM search S1+S2): [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionElectron/XENON1T_DM_S1S2.txt), [reference](https://arxiv.org/abs/2006.09721)
93+
* XENON1T (Solar axion basin): [limit](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionElectron/XENON1T_DM_SolarAxionBasin.txt), [reference](https://arxiv.org/abs/2006.12431)
9294

9395
### Haloscopes:
9496
<!---* Electron spin magnetometers: [projection](https://github.com/cajohare/AxionLimits/raw/master/limit_data/AxionElectron/Projections/ElectronSpinMagnetometers.txt), (in prep.)-->

‎limit_data/AxionElectron/XENON1T_DM_S1S2.txt‎

Lines changed: 151 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,151 @@
1+
# XENON1T ALP DM search (S1+S2 analysis)
2+
# m_a [eV] g_ae [dimensionless]
3+
1009.3086652232158 2.112486108036773e-14
4+
1060.530742821423 2.253884958697497e-14
5+
1114.3523237468867 2.4047483141860204e-14
6+
1170.905331925094 2.5657096793095987e-14
7+
1222.738640075412 2.741996401113324e-14
8+
1284.7922176005948 2.935267814229211e-14
9+
1349.9949934560434 3.1369464159662915e-14
10+
1418.506788405017 3.3469173199710197e-14
11+
1481.3008482423143 3.576879847393283e-14
12+
1556.4763714586693 3.8162976108752464e-14
13+
1635.4670273656984 4.085291848732439e-14
14+
1718.4664326730035 4.365987171914255e-14
15+
1794.539021738994 4.658223682137684e-14
16+
1885.6112776223779 4.994852584609645e-14
17+
1981.3054200689458 5.302688384129695e-14
18+
2081.8560082778317 5.5275567641807465e-14
19+
2201.0877480823506 5.714298132269343e-14
20+
2327.148109905053 5.88775338252279e-14
21+
2460.428181544738 6.05640406060622e-14
22+
2601.3414491209737 6.261011851756162e-14
23+
2750.325079867309 6.451062379825979e-14
24+
2907.841278392418 6.63584865975844e-14
25+
3074.3787206168354 6.814597661439105e-14
26+
3250.4540698338874 6.90576973893902e-14
27+
3436.6135795982404 7.021452066893802e-14
28+
3633.434788414836 7.186670180778643e-14
29+
3841.5283114857893 7.246614107491433e-14
30+
4061.5397350739486 7.009797161451962e-14
31+
4241.3348345690765 6.47253205294023e-14
32+
4401.766540869343 5.956621435290134e-14
33+
4596.622704982557 5.44552621896985e-14
34+
4785.276171462834 4.9976192980798346e-14
35+
5012.594393102949 4.6892811903603445e-14
36+
5266.981884084281 4.4538979874928406e-14
37+
5568.631667287945 4.279811535614e-14
38+
5887.557490870973 4.167550757006132e-14
39+
6224.748785583221 4.112529481324146e-14
40+
6581.251648701564 4.167550757006132e-14
41+
6958.172089426554 4.139948713825137e-14
42+
7356.679460149023 3.978133381598392e-14
43+
7730.028485235956 3.7596626477147126e-14
44+
8122.324848627962 3.5121093208120324e-14
45+
8429.557945145049 3.2001100275380147e-14
46+
8802.715247130589 2.954837957377501e-14
47+
9192.391371684447 2.906155358624124e-14
48+
9658.902163058034 3.1006783076686725e-14
49+
10086.479729374543 3.3636394434328224e-14
50+
10307.31495104639 3.7471916815756524e-14
51+
10598.365206228636 4.1883729918059596e-14
52+
10863.968832018641 4.673726638683464e-14
53+
11205.35316197475 5.163597229819835e-14
54+
11701.388586480725 5.620151892393207e-14
55+
12295.230148368979 5.976445567743231e-14
56+
12919.208971148875 6.250619200719417e-14
57+
13491.112882935613 6.803286113541094e-14
58+
13915.051372001344 7.491431608274467e-14
59+
14396.78621617917 8.262897903194691e-14
60+
14987.655600110444 8.963654576969398e-14
61+
15651.12492593042 8.415282752124532e-14
62+
16093.069692847517 7.562288289281618e-14
63+
16293.474572855348 6.739541086175394e-14
64+
16496.375068472757 6.006305540695147e-14
65+
16857.54970449753 5.338042694985124e-14
66+
17495.19930506458 5.112383846704688e-14
67+
18383.074825485233 5.400480875069855e-14
68+
19256.338505129646 5.810021467329499e-14
69+
20233.591251640977 6.168096620365035e-14
70+
21326.320750033643 6.408336444917617e-14
71+
22547.71853957586 6.66900319270737e-14
72+
23256.24761441874 7.282820193331213e-14
73+
23618.81811602753 8.153808605913877e-14
74+
24061.371830002692 9.189827833158399e-14
75+
24361.004928376573 1.0357483194030279e-13
76+
24664.36931831096 1.1596185730415935e-13
77+
25282.478412618882 1.3040674227102533e-13
78+
26483.489143344246 1.3669168529416178e-13
79+
27570.418037022857 1.295430473237702e-13
80+
28348.930893500725 1.1615466238119864e-13
81+
29239.75458120567 1.032312696734139e-13
82+
30065.40491523978 1.1443089350999293e-13
83+
30439.805447855528 1.295430473237702e-13
84+
30723.662358594844 1.4535735552438523e-13
85+
31202.651654602767 1.6418967920260103e-13
86+
31787.30620030351 1.8464210281813085e-13
87+
32583.92186853278 2.049008439267139e-13
88+
33921.22323228361 2.2166340442860606e-13
89+
34343.63976320848 1.9151526238900976e-13
90+
34343.63976320848 1.682393657943967e-13
91+
34450.063272449326 1.4960383941380625e-13
92+
34556.81656512015 1.3244449090281905e-13
93+
34771.31659162848 1.175133037675352e-13
94+
35422.83840875893 1.0707397390143385e-13
95+
36990.926663837075 1.1596185730415935e-13
96+
37800.81273993419 1.2983030816965963e-13
97+
38271.541288677865 1.4632648830384718e-13
98+
38628.430609079274 1.6382639526066958e-13
99+
38868.20383512307 1.8525660712102295e-13
100+
39230.657145986755 2.0949011560307784e-13
101+
39596.49040214231 2.3689362132558843e-13
102+
40213.8092243823 2.678817931967696e-13
103+
41272.647197787475 2.939993108568314e-13
104+
42385.58461688223 2.7785348590754925e-13
105+
42913.406726730136 2.4982785616689494e-13
106+
43313.582315355896 2.2092813717993816e-13
107+
43852.960658474374 1.951552372842763e-13
108+
44674.64828538584 1.7392311553684984e-13
109+
45511.73214885879 1.5500096507264727e-13
110+
46364.500733339744 1.387507872100456e-13
111+
47379.61354493188 1.2351837289316054e-13
112+
48566.9850782951 1.1087537184043357e-13
113+
49630.3194827121 9.908653889733159e-14
114+
50248.35932287419 8.772145321270649e-14
115+
50874.09553186896 7.783213713191108e-14
116+
52149.041617337905 6.96722777966668e-14
117+
54289.32984766643 7.150942098722502e-14
118+
56692.59464332828 7.770294359901126e-14
119+
58293.43775487094 8.523149213093274e-14
120+
59019.35826733234 9.606096690549871e-14
121+
59939.48427761564 1.0850650992656199e-13
122+
60311.53883137397 1.2229329791709483e-13
123+
60498.431241567094 1.387507872100456e-13
124+
61251.810247874186 1.5638039936952845e-13
125+
62786.8302789329 1.7742512375005979e-13
126+
62942.432348253635 2.002346192442627e-13
127+
64360.31915664624 1.8220440148370026e-13
128+
64559.757861097394 1.616635092436319e-13
129+
64960.49123728906 1.4343830340072084e-13
130+
65363.71203357204 1.272677240444581e-13
131+
65973.24094150437 1.1292014196659139e-13
132+
67626.58385413409 1.010816841197818e-13
133+
70839.09536370885 9.563635069187179e-14
134+
72614.38053365429 1.0426539049097599e-13
135+
73518.63786276784 1.175133037675352e-13
136+
73974.98052791749 1.3244449090281905e-13
137+
74664.81113411573 1.4828418127337042e-13
138+
76063.83073400985 1.6712510216327152e-13
139+
79677.14250830266 1.7950043533171298e-13
140+
83720.730874267 1.893010402264567e-13
141+
87426.85665886912 2.0365651373393526e-13
142+
90733.84380733279 2.2314128816186317e-13
143+
94165.92025235867 2.357159194791977e-13
144+
98030.6549543512 1.97328632114584e-13
145+
98395.2987698412 2.137650911110506e-13
146+
99867.4882040226 1.829323131593249e-13
147+
100176.95594323197 2.574112815165452e-13
148+
100176.95594323197 2.376820230357642e-13
149+
100487.38265601033 2.1946491133389848e-13
150+
100487.38265601033 2.0264404809255733e-13
151+
100487.38265601033 1e0

limit_data/AxionElectron/XENON1T.txt renamed to limit_data/AxionElectron/XENON1T_DM_S2.txt

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -111,4 +111,3 @@
111111
5185.731498551633 8.052973387032824e-13
112112
5306.322345981944 9.372876582529338e-13
113113
5457.927107530393 1.0874803051256767e-12
114-
5506.584544044661 1.4069541080088371e-12

0 commit comments

Comments
 (0)