<<peeters` ;

((peeters`setGitDir[ "blogit" ] )^\[Placeholder])

ps4 p2 (d)

Module[ {omega, sp, sm, s},

omega[qa_, pm_, deltaOverK0_] =  Sqrt[1 + pm Sqrt[Cos[ qa /2]^2 + deltaOverK0^2 Sin[qa/2]^2] ]
;

sp = Plot[ Evaluate@Table[omega[qa, sign, dOk],{sign, {(*-1,*)1}},
{dOk, {0.9, 0.1, 0}}
],
 {qa, 0, 2 Pi} ,
Ticks->{{0,Pi,2 Pi},Automatic},
AxesLabel -> { q a, ω /Sqrt[2 (Subscript[k, 0 ])/m]},
PlotLegends->Placed[{0.9, 0.1, 0},{Right,Bottom}]
] ;

sm = Plot[ Evaluate@Table[omega[qa, sign, dOk],{sign, {-1}},
{dOk, {0.9, 0.1, 0}}
],
 {qa, 0, 2 Pi} ,
Ticks->{{0,Pi,2 Pi},Automatic},
AxesLabel -> { q a, ω /Sqrt[2 (Subscript[k, 0 ])/m]},
PlotLegends->Placed[{0.9, 0.1, 0},{Right,Bottom}]
] ;

s = Plot[ Evaluate@Table[omega[qa, sign, dOk],{sign, {-1, 1}},
{dOk, {0.9, 0.1, 0}}
],
 {qa, 0, 2 Pi} ,
Ticks->{{0,Pi,2 Pi},Automatic},
AxesLabel -> { q a, ω /("2 \!\(\*SubscriptBox[\(k\), \(\(0\)\(\\ \)\)]\)/m")^(1/2)}
] ;

(*{sp, sm,*)
peeters`lblPlot[s]
]

Module[ {omega, sp, sm, s, p},

omega[qa_, pm_, deltaOverK0_] =  Sqrt[1 + pm Sqrt[Cos[ qa /2]^2 + deltaOverK0^2 Sin[qa/2]^2] ]
;

p = Plot[ Evaluate@Table[omega[qa, sign, dOk],{sign, {(*-1*),1}},
{dOk, {0}}
],
 {qa, 0, 2 Pi} ,
Ticks->{{0,Pi,2 Pi},Automatic},
AxesLabel -> { q a, ω /("2 \!\(\*SubscriptBox[\(k\), \(\(0\)\(\\ \)\)]\)/m")^(1/2)},
(*PlotLegends->Placed[{0.9, 0.1, 0},{Right,Bottom}],*)
PlotRange -> Full
]  ;
{p, peeters`exportForLatex["qmSolidsPs4OpticalDeltaZeroFig1", p]}
]
Module[ {omega, sp, sm, s, p},

omega[qa_, pm_, deltaOverK0_] =  Sqrt[1 + pm Sqrt[Cos[ qa /2]^2 + deltaOverK0^2 Sin[qa/2]^2] ]
;

p =Plot[ Evaluate@Table[omega[qa, sign, dOk],{sign, {-1}},
{dOk, {0}}
],
 {qa, 0, 2 Pi} ,
Ticks->{{0,Pi,2 Pi},Automatic},
AxesLabel -> { q a, ω /("2 \!\(\*SubscriptBox[\(k\), \(\(0\)\(\\ \)\)]\)/m")^(1/2)},
(*PlotLegends->Placed[{0.9, 0.1, 0},{Right,Bottom}],*)
PlotRange -> Full
] 
;
{p, peeters`exportForLatex["qmSolidsPs4AccousticDeltaZeroFig1", p]}
]