19.2 old set

19.2.1 P=(0,0,0)

19.2.1.1 ../../momentum/out/G2t_T128_L24_msq0-4.900000_msq1-4.650000_l02.500000_l12.500000_mu5.000000_g0.000000_rep0_output

GEVP_phi0_phi03_phi1_meffl0(L24T128) = 0.129046(40) χ2/dof= 0.026755

GEVP_phi0_phi03_phi1_meffl1(L24T128) = 0.40693(27) χ2/dof= 0.58826

GEVP_phi0_phi03_phi1_meffl2(L24T128) = 0.446792(68) χ2/dof= 0.95693

E1_0(L24T128) = 0.129057(40) χ2/dof= 0.047287

E1_1(L24T128) = 0.446791(68) χ2/dof= 0.96327

E3_0_vev(L24T128) = 0.40752(82) 0.00048373(72) 0.26(27) 0.0005905(66) χ2/dof= 0.0026487

0102030405060−0.500.511.522.533.54
mylabelGEVP_phi0_phi03_phi1_meffl0(L24T128)GEVP_phi0_phi03_phi1_meffl1(L24T128)GEVP_phi0_phi03_phi1_meffl2(L24T128)E1_0(L24T128)E1_1(L24T128)E3_0_vev(L24T128)fitty

19.2.2 P=(1,0,0)

19.2.2.1 ../../momentum/out/G2t_T128_L24_msq0-4.900000_msq1-4.650000_l02.500000_l12.500000_mu5.000000_g0.000000_rep0_output

GEVP_phi0_phi03_phi1_p1_meffl0(L24T128) = 0.290220(31) χ2/dof= 1.6069

GEVP_phi0_phi03_phi1_p1_meffl1(L24T128) = 0.515128(60) χ2/dof= 0.61451

GEVP_phi0_phi03_phi1_p1_meffl2(L24T128) = 0.56689(19) χ2/dof= 1.1024

E1_0_p1(L24T128) = 0.290237(13) χ2/dof= 4.9672

E1_1_p1(L24T128) = 0.515041(44) χ2/dof= 0.8219

E3_0_p1_vev(L24T128) = 0.5713(13) 0.00018486(53) 0.381(73) 50(22) 0.9(0.2)e-4 χ2/dof= 0.01676

E3_0_p1_meff(L24T128) = 0.440918(53) χ2/dof= 37945

0102030405060−4−20246
mylabelGEVP_phi0_phi03_phi1_p1_meffl0(L24T128)GEVP_phi0_phi03_phi1_p1_meffl1(L24T128)GEVP_phi0_phi03_phi1_p1_meffl2(L24T128)E1_0_p1(L24T128)E1_1_p1(L24T128)E3_0_p1_meff(L24T128)fitty

19.2.2.2 ../../momentum/out/G2t_T256_L24_msq0-4.900000_msq1-4.650000_l02.500000_l12.500000_mu5.000000_g0.000000_rep0_output

GEVP_phi0_phi03_phi1_p1_meffl0(L24T256) = 0.29002(33) χ2/dof= 0.29721

GEVP_phi0_phi03_phi1_p1_meffl1(L24T256) = 0.4724(11) χ2/dof= 184.03

GEVP_phi0_phi03_phi1_p1_meffl2(L24T256) = 0.51572(51) χ2/dof= 1.3818

E1_0_p1(L24T256) = 0.29023(19) χ2/dof= 0.57503

E1_1_px(L24T256) = 0.51572(51) χ2/dof= 1.3828

E3_0_p1_vev(L24T256) = 0.5711(67) 0.0001847(13) 1588(149) -3.1(1.5)e+9 8.8(1.0)e-5 χ2/dof= 0.0028509

E2_0(L24T256) = 0.26551(85) χ2/dof= 0.91054

E2_0_p1(L24T256) = 0.4264(16) 0.0018022(22) -879(971) χ2/dof= 0.015153

E1_0(L24T256) = 0.12967(35) χ2/dof= 0.20384

E3_0_p1_meff(L24T256) = 0.43477(93) χ2/dof= 290.71

020406080100120−6−4−20246
mylabelGEVP_phi0_phi03_phi1_p1_meffl0(L24T256)GEVP_phi0_phi03_phi1_p1_meffl1(L24T256)GEVP_phi0_phi03_phi1_p1_meffl2(L24T256)E1_0_p1(L24T256)E1_1_px(L24T256)E3_0_p1_vev(L24T256)E3_0_p1_meff(L24T256)fitty
#> <div id="htmlwidget-b2551f2f9d7bcd218b79" style="width:900px;height:480px;" class="plotly html-widget"></div>
#> <script type="application/json" data-for="htmlwidget-b2551f2f9d7bcd218b79">{"x":{"data":[{"visible":false,"showlegend":false,"xaxis":null,"yaxis":null,"hoverinfo":"text","frame":null}],"layout":{"margin":{"t":29.6986301369863,"r":7.30593607305936,"b":17.351598173516,"l":10.958904109589},"plot_bgcolor":"rgba(235,235,235,1)","paper_bgcolor":"rgba(255,255,255,1)","font":{"color":"rgba(0,0,0,1)","family":"","size":14.6118721461187},"xaxis":{"domain":[0,1],"automargin":true,"type":"linear","autorange":true,"range":[],"tickmode":"auto","ticktext":[],"tickvals":[],"categoryorder":"array","categoryarray":[],"nticks":null,"ticks":"outside","tickcolor":"rgba(51,51,51,1)","ticklen":3.65296803652968,"tickwidth":0.66417600664176,"showticklabels":true,"tickfont":{"color":"rgba(77,77,77,1)","family":"","size":11.689497716895},"tickangle":-0,"showline":false,"linecolor":null,"linewidth":0,"showgrid":true,"gridcolor":"rgba(255,255,255,1)","gridwidth":0.66417600664176,"zeroline":false,"anchor":"y","title":"t","hoverformat":".2f","showexponent":"all","exponentformat":"e"},"yaxis":{"domain":[0,1],"automargin":true,"type":"linear","autorange":true,"range":[],"tickmode":"auto","ticktext":[],"tickvals":[],"categoryorder":"array","categoryarray":[],"nticks":null,"ticks":"outside","tickcolor":"rgba(51,51,51,1)","ticklen":3.65296803652968,"tickwidth":0.66417600664176,"showticklabels":true,"tickfont":{"color":"rgba(77,77,77,1)","family":"","size":11.689497716895},"tickangle":-0,"showline":false,"linecolor":null,"linewidth":0,"showgrid":true,"gridcolor":"rgba(255,255,255,1)","gridwidth":0.66417600664176,"zeroline":false,"anchor":"x","title":"y","hoverformat":".2f","showexponent":"all","exponentformat":"e"},"shapes":[{"type":"rect","fillcolor":null,"line":{"color":null,"width":0,"linetype":[]},"yref":"paper","xref":"paper","x0":0,"x1":1,"y0":0,"y1":1}],"showlegend":false,"legend":{"bgcolor":"rgba(255,255,255,1)","bordercolor":"transparent","borderwidth":1.88976377952756,"font":{"color":"rgba(0,0,0,1)","family":"","size":11.689497716895}},"hovermode":"closest","width":900,"height":480,"barmode":"relative","title":"fit"},"config":{"doubleClick":"reset","modeBarButtonsToAdd":["hoverclosest","hovercompare"],"showSendToCloud":false},"source":"A","attrs":{"4bc734b64e3":{"type":"scatter"}},"cur_data":"4bc734b64e3","visdat":{"4bc734b64e3":["function (y) ","x"]},"highlight":{"on":"plotly_click","persistent":false,"dynamic":false,"selectize":false,"opacityDim":0.2,"selected":{"opacity":1},"debounce":0},"shinyEvents":["plotly_hover","plotly_click","plotly_selected","plotly_relayout","plotly_brushed","plotly_brushing","plotly_clickannotation","plotly_doubleclick","plotly_deselect","plotly_afterplot","plotly_sunburstclick"],"base_url":"https://plot.ly"},"evals":[],"jsHooks":[]}</script>

19.2.3 plotting the exponential

|A0ϕ1|2(eEϕ1t+eEϕ1(Tt))+|A0ϕ0|2(eEϕ0t+eEϕ0(Tt))+|A03ϕ0|2(eE3ϕ0t+eE3ϕ0(Tt))+|Aϕ02ϕ0(p)|2eEϕ0T(e(E2ϕ0(p)Eϕ0)t+e(E2ϕ0(p)Eϕ0)(Tt))+|Aϕ0(p)2ϕ0|2eEϕ0(p)T(e(E2ϕ0Eϕ0(p))t+e(E2ϕ0Eϕ0(p))(Tt)).

library(ggplot2)
library(Rose)
T<-256
t<-c(1:100)*48/100
E1_0= 0.12967
E1_0_p1= 0.29023
E1_1_p1=0.51572
E2_0=0.26551
E2_0_p1=0.4264

e1<-log10( exp(-E1_1_p1*t)+exp(-E1_1_p1*(T-t)))
e2<-log10( exp(-E1_0_p1*t)+exp(-E1_0_p1*(T-t)))
e3<-log10( (exp(-(E2_0_p1-E1_0)*t)+exp(-(E2_0_p1-E1_0)*(T-t)) ) *exp(-E1_0*T))
e4<-log10( (exp(-(E2_0-E1_0_p1)*t)+exp(-(E2_0-E1_0_p1)*(T-t)) ) *exp(-E1_0_p1*T))

gg<-ggplot()
gg<- gg+ geom_line(aes(x=t, y=e1, color="e1"))
gg<- gg+ geom_line(aes(x=t, y=e2, color="e2"))
gg<- gg+ geom_line(aes(x=t, y=e3, color="e3"))
gg<- gg+ geom_line(aes(x=t, y=e4, color="e4"))

gg<- gg+ geom_line(aes(x=(corr_E3_0_p1[,1]),
                       y=log10(corr_E3_0_p1[,2]), color="corr_E3"))

myplotly(gg, to_print = FALSE, xrange=c(0,48))
010203040−30−25−20−15−10−50
coloure1e2e3e4corr_E3