OCTOptimizeHarmonicSpectrum

OCTOptimizeHarmonicSpectrum

Section Calculation Modes::Optimal Control
Type block
Default no

(Experimental) If OCTTargetOperator = oct_tg_hhg, the target is the harmonic emission spectrum. In that case, you must supply an OCTOptimizeHarmonicSpectrum block in the inp file. The target is given, in general, by:

$J_1 = \int_0^\infty d\omega \alpha(\omega) H(\omega)$,

where $H(\omega)$ is the harmonic spectrum generated by the system, and $\alpha(\omega)$ is some function that determines what exactly we want to optimize. The role of the OCTOptimizeHarmonicSpectrum block is to determine this $\alpha(\omega)$ function. Currently, this function is defined as:

$\alpha(\omega) = \sum_{L=1}^{M} \frac{\alpha_L}{a_L} \sqcap( (\omega - L\omega_0)/a_L )$,

where $\omega_0$ is the carrier frequency. $M$ is the number of columns in the OCTOptimizeHarmonicSpectrum block. The values of L will be listed in the first row of this block; $\alpha_L$ in the second row, and $a_L$ in the third.

Example:

%OCTOptimizeHarmonicSpectrum
   7 | 9 | 11
   -1 | 1 | -1
   0.01 | 0.01 | 0.01
%

 if (parse_block(namespace, 'OCTOptimizeHarmonicSpectrum', blk) == 0) then