Gabriel Turinici

We will proceed in 4 steps:

Note: the input example provided below can be recovered by pressing the "Reload/Refresh" button of your browser. The evolution equation is:
i d/dt Psi(t) = A Psi(t) + e1(t)*(Bx Psi(t)) + e2(t)*(By Psi(t)) +e3(t)*(Bz Psi(t))
where i=sqrt(-1), Psi(t)=wavefunction, A=internal hamiltonian, Bx,By,Bz=dipole moment matrices, e1(t),e2(t),e3(t)=laser intensities.

1. Input the system dimension (click here for help)
2. Input the system internal hamiltonian as a matrix prefixed by "complex" and suffixed by its name (one word) (click here for help) complex 0 0 0 0 0 2 0 0 0 0 2 0 0 0 0 2 A
3. Input the number of dipole matrices to be used (click here for help)
4. Input the system dipole matrices one by one, prefixed by "complex" and suffixed by its name (one word) as in exemple below (click here for help) :complex -0.000 -0.892 -0.000 -0.000 -0.892 +0.000 +0.000 +0.000 -0.000 -0.000 +0.000 +0.000 -0.000 +0.000 +0.000 +0.000 Bx complex +0.000 +0.000 +0.282 +0.000 +0.000 -0.000 -0.000 +0.000 +0.282 -0.000 +0.000 -0.000 +0.000 -0.000 -0.000 -0.000 By complex -0.000 -0.000 -0.000 -0.892 -0.000 +0.000 +0.000 +0.000 -0.000 +0.000 +0.000 +0.000 -0.892 +0.000 -0.000 +0.000 Bz

Back to the introduction page.

 

---

System dimension: the input should be the number of eigenstates used to describe the system. It is equal to the dimension of the hamiltonian and dipole matrices.

Back to input

 

---

System internal hamiltonian: the input should be the hamiltonian matrix. This should be a skew-hermitian matrix. However, ussualy this is a real symmetric matrix multiplied by "i"(=sqrt(-1)). Then the matrix is prefixed by keyword "complex" and the real symmetric matrix is provided. Note that all (not only under or upper diagonal!!!) entries are required. Entries are followed by a name to give to this matrix, for instance "A". Format for entries and name is free.

Back to input

 

---

The number of dipole matrices: input here the number of dipol moment matrices. Usually this equals 1, but any general cases may be considered.

Back to input

 

---

Dipole matrices: input here all dipol moment matrices. Their number should equal the previous input. Convention is as for internal hamiltonian matrix. Names are at users's choice, for instance "Bx","By","Bz". Format is free.

Back to input

(c) Copyright Gabriel Turinici, 2001-2010. All rights reserved. Notice: This page is work in progress intended to be released with public academic licence. As a developpement version, its content may NOT be copied or distributed in any way without prior permission of the author and of ALL other copyright owners.