I have the following code
Clear[Rg3, SetRg3, Ig3, SetIg3, κ1, Setκ1, Γ, SetΓ, κ2, Setκ2, g1, Setg1, g2, Setg2, r1]
Setκ1 = 1; SetΓ = 0.01; Setκ2 = 20;
r1 = RegionPlot[SetRg3 = 0; SetIg3 = 0;
NMG3 = {{-Γ/
2, -I*g1, -I*Rg3 + Ig3}, {-I*g1, -κ1/
2, -I*g2}, {-I*Rg3 - Ig3, - I*g2, -κ2/2}} /. {Rg3 ->
SetRg3, Ig3 -> SetIg3, κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2, g1 -> Setg1,
g2 -> Setg2};
EigensysNMG3 = Eigensystem[NMG3, Cubics -> True]; {Chop[Min[{((Abs[
Normalize[
EigensysNMG3[[2]][[1]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}]), ((Abs[
Normalize[
EigensysNMG3[[2]][[2]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}]), ((Abs[
Normalize[
EigensysNMG3[[2]][[3]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}])}] >=
RankedMax[{((Abs[
Normalize[
EigensysNMG3[[2]][[1]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}])[[
1]], ((Abs[
Normalize[
EigensysNMG3[[2]][[2]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}])[[
1]], ((Abs[
Normalize[
EigensysNMG3[[2]][[3]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}])[[1]]},
2]]}, {Setg1, 0.01, 10}, {Setg2, 0.01, 12}, PlotRange -> Full, PlotLegends -> Automatic, PlotPoints -> Automatic, PlotStyle -> Directive[Red, Opacity[0.35]], PlotRangePadding -> None, BoundaryStyle -> {Black, Thick}]
which generates the following region plot
My goal is to only plot the line bounding the red and white region. How should I go about going that? I tried using ContourPlot and threw in the following in addition to the code above:
ContourPlot[{Chop[Min[{((Abs[
Normalize[
EigensysNMG3[[2]][[1]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}]), ((Abs[
Normalize[
EigensysNMG3[[2]][[2]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}]), ((Abs[
Normalize[
EigensysNMG3[[2]][[3]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}])}] ==
RankedMax[{((Abs[
Normalize[
EigensysNMG3[[2]][[1]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}])[[1]], ((Abs[
Normalize[
EigensysNMG3[[2]][[2]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}])[[1]], ((Abs[
Normalize[
EigensysNMG3[[2]][[3]]] /. {κ1 ->
Setκ1, Γ ->
SetΓ, κ2 -> Setκ2,
g3 -> Setg3}]).Transpose[{{1, 0, 0}}])[[1]]},
2]]}, {Setg1, 0.01, 10}, {Setg2, 0.01, 12}, PlotPoints -> Automatic]
but I was returned with a white blank plot. The inequality for the RegionPlot was a >=
and I tried doing the ContourPlot for =
. Appreciate any constructive help that I can take. Thanks for reading.
PlotStyle -> None
inRegionPlot
? $\endgroup$