Reference.Programs.qm_balance
qm_balance: Perform balancing transformation
If the dimensionality of the problems under investigation is growing (particularly so for density matrix description in LvNE simulations), efficient methods for dimension reduction become mandatory. Instead of empirical truncation of the state space, we implemented here the method of balanced truncation, established as a powerful technique for the control of linear systems. This procedure is based on states identified as being simultaneously easily controllable and observable and on restricting the dynamics to the subspace spanned by such states.
As we have already demonstrated for some test cases, the procedure can greatly reduce the number of equations to be solved without notably changing the overall response of the system to the external control fields, see [1].
Source code
The MATLAB function qm_balance.m can be found here
Input parameters
The following input parameters are required when calling the function qm_balance.m
| parameter | description |
|---|---|
| filename | Typically used for equation of motion: 'tdse' or 'lvne' but other choices are also possible |
File I/O
As input this function requires matrices A,B,N,C,D along with vectors xi, xe which are typically generated using our function qm_abncd. These data are read from input data file tdse.mat or lvne.mat for TDSE or LvNE simulations, respectively,.
As output this function produces transformed matrices A,B,N,C,D along with vectors xi, xe which are saved as output data file tdse_b.mat or lvne_b.mat for TDSE or LvNE simulations, respectively. These matrices/vectors can then be truncated by running our function qm_truncate ...
Algorithms
For more details on the balancing transformation, see this page of our WavePacket central wiki.
-
Two options for stabilizing A-matrix: EVS and SSU
-
Two options for solving generalized Lyapunov equations: by iterative method or by biconjugate gradient method
-
Two options for performing the balancing transform: SRBT and MRMR
Variables
Note that function qm_init should be run previously in order to initialize all variables in use. Of particular interest are the following structures:
References
- B. Schäfer-Bung, C. Hartmann, B. Schmidt, and Ch. Schütte: J. Chem. Phys. 135, 014112 (2011)
