Group Abstract Group Abstract

Message Boards Message Boards

0
|
341 Views
|
4 Replies
|
2 Total Likes
View groups...
Share
Share this post:

Quantum optics. How to analytically solve the equations below?

Recently, I have read a paper about Dicke Model, named" Dissipative Phase Transition in the Two-Photon Dicke Model".

enter image description here

"These equations can be solved analytically to determine the steady ..."but i cannot use "Solve" function to solve these analytically. By the way, the left of all equation is zero because of steady states. My Notebook is below, i want to know how to get the analytical solve.

POSTED BY: Tree Clean
4 Replies

Like previous comments said that system of equations cannot solved analytically, that is why they say they do numerical evolution. So the approach to be preferred should be using NSolve or NDSolve.

Also relevant can be exploring the second quantization functions of the Wolfram QuantumFramework https://wolfr.am/1tscRdMNP. Depending on the specific problem you want to simulate some accurate numerical simulations can be done in a truncated Fock space.

POSTED BY: Bruno Tenorio

There are 7 equations in 10 variables. I have assigned a numeric value to three of them, and then found a numerical solution with respect to the others using FindRoot. Unfortunately, k2 has come out negative:

eqa1 = -(\[Kappa]1/2 + I Subscript[\[Omega], c]) a - 
    2 I \[Lambda] jx a\[Conjugate] -
    \[Kappa]2 (2 n a + A a\[Conjugate] - 2 a^2 a\[Conjugate]) == 0;
eqa2 = -(\[Kappa]1/2 - I Subscript[\[Omega], c]) a\[Conjugate] + 
    2 I \[Lambda] jx a -
    \[Kappa]2 (2 n a\[Conjugate] + A\[Conjugate] a - 
       2 a\[Conjugate]^2 a) == 0;
eqa3 = -(\[Kappa]1 + 2 I Subscript[\[Omega], c]) A - 
    4 I \[Lambda] jx n -
    2 \[Kappa]2 (3 n A - 2 a\[Conjugate] a^3) == 0;
eqa4 = -(\[Kappa]1 - 2 I Subscript[\[Omega], c]) A\[Conjugate] + 
    4 I \[Lambda] jx n -
    2 \[Kappa]2 (3 n A\[Conjugate] - 2 a a\[Conjugate]^3) == 0;
eqa5 = 2 I \[Lambda] jx (A - A\[Conjugate]) - \[Kappa]1 n -
    \[Kappa]2 (4 n^2 - 4 a\[Conjugate]^2 a^2 + 2 A\[Conjugate] A) == 0;
eqa6 = 2 Subscript[\[Omega], a] jx - 
    2 \[Lambda] jz (A + A\[Conjugate]) == 0;
eqa7 = jx^2 + jz^2 - 1 == 0;
eqs = {eqa1, eqa2, eqa3, eqa4, eqa5, eqa6, eqa7};
assigned = {Subscript[\[Omega], a] -> 1, 
   Subscript[\[Omega], c] -> 1, \[Lambda] -> 1};
startingValues = {{a, 2}, {A, 2}, {jx, 2}, {jz, 2},
   {n, 2}, {\[Kappa]1, 2}, {\[Kappa]2, 2}};
sol = FindRoot[eqs /. assigned, startingValues]
Chop[sol, 10^-7]
eqs[[All, 1]] /. sol // Chop
POSTED BY: Gianluca Gorni

I cannot solve the problem, but one issue is the use of Subscript[] in algebraic equations. Consider your $\omega_a$. The actual code representing what you see is the following:

Subscript[\[Omega], a]

Note the a. That is the same variable a that appears multiple times in each equation. Subscript[\[Omega], a] is treated algebraically as a function Subscript[] with two arguments, \[Omega] and a. Since Subscript[] has no definition, it is impossible to an equation containing Subscript[\[Omega], a] for the variable a.

Note 1: Fixing that problem did not make the system solvable in the time I was willing to wait.

Note 2: I am unfamiliar with the notation and subject matter of the paper. Therefore, I cannot comment on the correctness of your system beyond the remarks about syntax above.

POSTED BY: Michael Rogers

The article also says that a closed-form exact solution was not possible.

POSTED BY: Gianluca Gorni
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard