4
$\begingroup$

I somehow know that when simplifying the expression $$ \frac{2}{3} \left(2 \sqrt{3 B^2+4 J^2} \cos \left(\frac{1}{3} \cos ^{-1}\left(-\frac{J \left(9 B^2+8 J^2\right)}{\sqrt{\left(3 B^2+4 J^2\right)^3}}\right)\right)+J\right) $$ it is equal $2 \sqrt{B^2+J^2}$, as they are identical

J = -1;
Plot[{2 Sqrt[B^2 + J^2], 
  2/3 (J + 2 Sqrt[3 B^2 + 4 J^2]
       Cos[1/3 ArcCos[-((J (9 B^2 + 8 J^2))/
          Sqrt[(3 B^2 + 4 J^2)^3])]])}, {B, 0, 5}, 
 PlotStyle -> {{Red, Dashed}, {Yellow, Thin}}]

enter image description here

How can I perform this simplifying using Mathematica? where $B$ and $J$ are reals.

$\endgroup$
1
  • 1
    $\begingroup$ It seems Mathematica has difficulty simplifying because ArcCos has multiple inverses. You need the Cos[...] part to equal (-2 J + 3 Sqrt[B^2 + J^2])/(4 Sqrt[3 B^2 + 4 J^2]) however, I haven't got a result out of Reduce[Cos[ArcCos[-((J (9 B^2 + 8 J^2))/Sqrt[(3 B^2 + 4 J^2)^3])]/3] == (-2 J + 3 Sqrt[B^2 + J^2])/(4 Sqrt[3 B^2 + 4 J^2]), {B, J}, Reals] $\endgroup$
    – flinty
    Jul 9, 2020 at 13:10

2 Answers 2

4
$\begingroup$
expr = 2/3 (J+2 Sqrt[3 B^2+4 J^2] Cos[1/3 ArcCos[-((J (9 B^2+8 J^2))/Sqrt[(3 B^2+4 J^2)^3])]]);
expr2 = Simplify[expr /. {B->r Cos[t]/√3, J->r Sin[t]/2 },{r>0, 0<t<2Pi}]
expr3 = expr2 /. e:Cos[ArcCos[_]/3 ] :> First[y/.Solve[Reduce[{y==e, 0<t<2Pi},t,Reals]/.c:_Cos:>TrigExpand[c],y,Reals]//Simplify//Normal//Union]//Simplify

Simplify[TrigExpand //@ expr3 /. {r->Sqrt[3 B^2 + 4 J^2], t->ArcTan[√3B, 2J]}, _∈ Reals]

Output

$\frac{1}{3} r \left(4 \cos \left(\frac{1}{3} \cos ^{-1}\left(-\frac{1}{4} (5+\cos (2 t)) \sin (t)\right)\right)+\sin (t)\right)$

$\frac{r \sqrt{\cos (2 t)+7}}{\sqrt{6}}$

$2 \sqrt{B^2+J^2}$

$\endgroup$
3
$\begingroup$

Another way to reduce expression.

expr = 2/3 (J + 
  2 Sqrt[3 B^2 + 4 J^2] Cos[
    1/3 ArcCos[-((J (9 B^2 + 8 J^2))/Sqrt[(3 B^2 + 4 J^2)^3])]]) //
Expand;

rule = Cos[
1/3 ArcCos[-((J (9 B^2 + 8 J^2))/Sqrt[(3 B^2 + 4 J^2)^3])]
] -> cc;

Solve for unknown simplified expression x

sol = Solve[x == expr /. rule, cc, Reals][[1, 1]] // Normal

(*   cc -> (-2 J + 3 x)/(4 Sqrt[3 B^2 + 4 J^2])   *)

eq1 = Equal @@ sol /. Reverse[rule]

(*   Cos[1/3 ArcCos[-((J (9 B^2 + 8 J^2))/
 Sqrt[(3 B^2 + 4 J^2)^3])]] == (-2 J + 3 x)/(
 4 Sqrt[3 B^2 + 4 J^2])   *)

eq2 = ArcCos[#] & /@ eq1 // 
   Simplify[#, 
0 < ArcCos[-((J (9 B^2 + 8 J^2))/Sqrt[(3 B^2 + 4 J^2)^3])] < Pi] &

(*    ArcCos[-((J (9 B^2 + 8 J^2))/Sqrt[(3 B^2 + 4 J^2)^3])] == 
      3 ArcCos[(-2 J + 3 x)/(4 Sqrt[3 B^2 + 4 J^2])]   *)

eq3 = TrigExpand[Cos /@ eq2]

(*   -((9 B^2 J Sqrt[(3 B^2 + 4 J^2)^3])/(3 B^2 + 4 J^2)^3) - (
8 J^3 Sqrt[(3 B^2 + 4 J^2)^3])/(3 B^2 + 4 J^2)^3 == -(J^3/(
2 (3 B^2 + 4 J^2)^(3/2))) + (3 J)/(2 Sqrt[3 B^2 + 4 J^2]) + (
9 J^2 x)/(4 (3 B^2 + 4 J^2)^(3/2)) - (9 x)/(
4 Sqrt[3 B^2 + 4 J^2]) - (27 J x^2)/(8 (3 B^2 + 4 J^2)^(3/2))  + (
27 x^3)/(16 (3 B^2 + 4 J^2)^(3/2))   *)

Solve[eq3, x, Reals] // Simplify

(*   {{x -> ConditionalExpression[2 J, B != 0]}, 
      {x -> ConditionalExpression[-2 Sqrt[B^2 + J^2], B != 0]}, 
      {x -> ConditionalExpression[2 Sqrt[B^2 + J^2], B != 0]}}   *)

I didn't examine further, why 2 J is given as solution for B!=0. It is the solution for B==0.

expr /. B -> 0 // Simplify[#, J \[Element] Reals] &
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.