Meep Tutorial/Multilevel-atomic susceptibility

From AbInitio

(Difference between revisions)
Jump to: navigation, search
Revision as of 03:12, 17 December 2016 (edit)
Ardavan (Talk | contribs)

← Previous diff
Revision as of 19:39, 17 December 2016 (edit)
Ardavan (Talk | contribs)

Next diff →
Line 15: Line 15:
(set! pml-layers (list (make pml (thickness dpml) (side High)))) (set! pml-layers (list (make pml (thickness dpml) (side High))))
(define-param freq-21 (/ 40 (* 2 pi))) ; emission frequency (units of 2\pia/c) (define-param freq-21 (/ 40 (* 2 pi))) ; emission frequency (units of 2\pia/c)
- (define-param gamma-21 (/ 4 (* 2 pi))) ; emission gamma (units of 2\pia/c)+ (define-param gamma-21 (/ 4 (* 2 pi))) ; emission linewidth (units of 2\pia/c)
- (define-param sigma-21 8e-23) ; emission sigma+ (define-param sigma-21 8e-23) ; dipole coupling strength
(set! sigma-21 (/ sigma-21 (sqr freq-21))) (set! sigma-21 (/ sigma-21 (sqr freq-21)))
(define-param rate-21 0.005) ; non-radiative rate (units of c/a) (define-param rate-21 0.005) ; non-radiative rate (units of c/a)
- (define-param N0 1e23) ; initial population of ground state+ (define-param N0 1e23) ; initial population density of ground state
- (define-param Rp 0.015) ; pumping rate of ground to highest state+ (define-param Rp 0.015) ; pumping rate of ground to excited state
(define two-level (make medium (index ncav) (define two-level (make medium (index ncav)
(E-polarizations (make multilevel-atom (sigma 1) (E-polarizations (make multilevel-atom (sigma 1)
Line 34: Line 34:
(define print-field (lambda () (print "field:, " (meep-time) ", " (define print-field (lambda () (print "field:, " (meep-time) ", "
(real-part (get-field-point Ex (vector3 0 0 (+ (* -0.5 sz) (* 0.5 Lcav))))) "\n"))) (real-part (get-field-point Ex (vector3 0 0 (+ (* -0.5 sz) (* 0.5 Lcav))))) "\n")))
- (define-param endt 50000)+ (define-param endt 20000)
(run-until endt (after-time (- endt 250) print-field)) (run-until endt (after-time (- endt 250) print-field))
-The choice of the atomic media parameters requires some care. For example, choosing a pumping rate <code>Rp</code> that lies far beyond threshold will cause large inversion which is not physical and thus produces meaningless results. For these types of problems involving steady-state phenomena, it is important to check that the fields used in the analysis do not contain relaxation oscillations. The choice of the simulation time is also important when operating near the lasing threshold of a particular mode as there the modal intensity is typically weak. For this particular cavity, we can compute the lasing thresholds analytically using the steady-state ab-initio laser theory (SALT) developed by Prof. A. Douglas Stone and his group at Yale. Based on the default parameters of the simulation script above, two modes, <math>m</math>=18 and 19, should begin to lase within the pumping-rate interval of [0.0050,0.0096]. For comparison, note that this is just above the relaxation rate of <code>rate-21</code>=0.005. We also need to choose a large resolution to ensure numerical stability of the polarization fields which is also a feature of linear absorbers. A plot of the variation of the modal intensity with pumping rate is shown below.+The choice of the atomic media parameters requires some care. For example, choosing a pumping rate <code>Rp</code> that lies far beyond threshold will cause large inversion which is not physical and thus produces meaningless results. The choice of the simulation time is also important when operating near the lasing threshold of a particular mode as the fields contain relaxation oscillations and require some time to reach steady state. For this particular two-level cavity, we can compute the lasing thresholds analytically using the steady-state ab-initio laser theory (SALT) developed by Prof. A. Douglas Stone and his group at Yale. Based on the default parameters of the simulation script above, two modes, <math>m</math>=18 and 19, should begin to lase very close to the relaxation rate. We also need to choose a small timestep to ensure that the data is smooth and continuous which therefore requires a large resolution. A plot of the variation of the modal intensity with pumping rate is shown below.
Note that the slopes of each curve is decreasing with increasing pumping rate. This occurs because the onset of lasing from additional modes means there is less gain available to the other modes. The modal intensities reach an asymptote in the limit of large pumping rates. We can convert Meep's dimensionless parameters into real units by specifying the units of the cavity length <math>L_{cav}</math> and then multiplying the rate terms by <math>L_{cav}/c</math>. Note that the slopes of each curve is decreasing with increasing pumping rate. This occurs because the onset of lasing from additional modes means there is less gain available to the other modes. The modal intensities reach an asymptote in the limit of large pumping rates. We can convert Meep's dimensionless parameters into real units by specifying the units of the cavity length <math>L_{cav}</math> and then multiplying the rate terms by <math>L_{cav}/c</math>.

Revision as of 19:39, 17 December 2016

Meep 1.4 introduced a feature to model saturable absorption/gain via multilevel-atomic susceptibility. This is based on a generalization of the Maxwell-Bloch equations which involves the interaction of a quantized system having an arbitrary number of levels with the electromagnetic fields. We will demonstrate this feature by computing the lasing thresholds of a two-level, multimode cavity in 1d. This is similar to the example used in A. Cerjan et al (Fig. 2).

The cavity consists of a high-index medium with a perfect-metallic mirror on one end and an abrupt termination in air on the other. We will specify an initial population density for the ground state of the two-level system. The field within the cavity is initialized to an arbitrary profile and a fictitious source is used to pump the cavity at a fixed rate. The fields are time stepped until reaching steady state. Near the end of the time stepping, we will output the electric field in the middle of the cavity and then, in post processing, compute its Fourier transform to obtain the spectra.

The angular frequency of the passive cavity modes can be computed analytically using the equation ωcav = (m + 0.5)π / (ncavLcav) where ncav and Lcav are the cavity index and length, and m is an integer. Similar to the paper referenced above, we focus on the angular-frequency range of [36,44] which contains four modes corresponding to m=17-20. In the continuum limit, these modes would appear as Dirac delta functions in the cavity spectra. The discretized model, however, produces peaks with finite width. We therefore need to integrate a fixed number of points around each peak to smoothen out the mode intensity.

(set-param! resolution 500)
(define-param ncav 1.5)                          ; cavity refractive index
(define-param Lcav 1)                            ; cavity length
(define-param dpad 1)                            ; padding thickness
(define-param dpml 1)                            ; PML thickness
(define-param sz (+ Lcav dpad dpml))
(set! geometry-lattice (make lattice (size no-size no-size sz)))
(set! dimensions 1)
(set! pml-layers (list (make pml (thickness dpml) (side High))))
(define-param freq-21 (/ 40 (* 2 pi)))           ; emission frequency  (units of 2\pia/c)
(define-param gamma-21 (/ 4 (* 2 pi)))           ; emission linewidth  (units of 2\pia/c)
(define-param sigma-21 8e-23)                    ; dipole coupling strength
(set! sigma-21 (/ sigma-21 (sqr freq-21)))
(define-param rate-21 0.005)                     ; non-radiative rate  (units of c/a)
(define-param N0 1e23)                           ; initial population density of ground state
(define-param Rp 0.015)                          ; pumping rate of ground to excited state
(define two-level (make medium (index ncav)
 (E-polarizations (make multilevel-atom (sigma 1)
  (transitions (make transition (from-level 1) (to-level 2) (pumping-rate Rp)
                                (frequency freq-21) (gamma gamma-21) (sigma sigma-21))
               (make transition (from-level 2) (to-level 1) (transition-rate rate-21)))
  (initial-populations N0)))))
(set! geometry (list (make block (center 0 0 (+ (* -0.5 sz) (* 0.5 Lcav)))
                          (size infinity infinity Lcav) (material two-level))))
(init-fields)
(meep-fields-initialize-field fields Ex 
         (lambda (p) (if (= (vector3-z p) (+ (* -0.5 sz) (* 0.5 Lcav))) 1 0)))
(define print-field (lambda () (print "field:, " (meep-time) ", "
      (real-part (get-field-point Ex (vector3 0 0 (+ (* -0.5 sz) (* 0.5 Lcav))))) "\n")))
(define-param endt 20000)
(run-until endt (after-time (- endt 250) print-field))

The choice of the atomic media parameters requires some care. For example, choosing a pumping rate Rp that lies far beyond threshold will cause large inversion which is not physical and thus produces meaningless results. The choice of the simulation time is also important when operating near the lasing threshold of a particular mode as the fields contain relaxation oscillations and require some time to reach steady state. For this particular two-level cavity, we can compute the lasing thresholds analytically using the steady-state ab-initio laser theory (SALT) developed by Prof. A. Douglas Stone and his group at Yale. Based on the default parameters of the simulation script above, two modes, m=18 and 19, should begin to lase very close to the relaxation rate. We also need to choose a small timestep to ensure that the data is smooth and continuous which therefore requires a large resolution. A plot of the variation of the modal intensity with pumping rate is shown below.

Note that the slopes of each curve is decreasing with increasing pumping rate. This occurs because the onset of lasing from additional modes means there is less gain available to the other modes. The modal intensities reach an asymptote in the limit of large pumping rates. We can convert Meep's dimensionless parameters into real units by specifying the units of the cavity length Lcav and then multiplying the rate terms by Lcav / c.

Personal tools