Materials in Meep

From AbInitio

Revision as of 07:53, 22 October 2005; Stevenj (Talk | contribs)
(diff) ←Older revision | Current revision | Newer revision→ (diff)
Jump to: navigation, search
Meep
Download
Release notes
FAQ
Meep manual
Introduction
Installation
Tutorial
Reference
C++ Tutorial
C++ Reference
Acknowledgements
License and Copyright

The material structure in Maxwell's equations is determined by the dielectric function ε(x), but ε is not only a function of position — in general, it also depends on frequency (material dispersion) and on the electric field E itself (nonlinearity). Material dispersion, in turn, is generally associated with absorption loss in the material, or possibly gain. All of these effects can be simulated in Meep, with certain restrictions.

In this section, we describe the form of the equations and material properties that Meep can simulate. The actual interface with which you specify these properties is described in the Meep reference.

Contents

Material dispersion

Physically, material dispersion arises because the polarization of the material does not respond instantaneously to an applied field E, and this is essentially the way that it is implemented in FDTD. In particular, \mathbf{D} = \varepsilon\mathbf{E} is expanded to:

\mathbf{D} = \varepsilon_\infty \mathbf{E} + \mathbf{P}

where \varepsilon_\infty is the instantaneous dielectric function (the infinite-frequency response) and P is the polarization density in the material. P, in turn, has its own time-evolution equation, and the exact form of this equation determines the frequency-dependence ε(ω). In particular, Meep supports any material dispersion of the form of a sum of harmonic resonances:

\varepsilon(\omega,\mathbf{x}) = \varepsilon_\infty(\mathbf{x}) + \sum_n \frac{\sigma_n(\mathbf{x})}{\omega_n^2 - \omega^2 - i\omega\gamma_n} ,

where ωn and γn are user-specified constants and \sigma_n(\mathbf{x}) is a user-specified function of position. This corresponds to evolving P via the equations:

\mathbf{P} = \sum_n \mathbf{P_n}
\frac{d^2\mathbf{P}_n}{dt^2} + \gamma_n \frac{d\mathbf{P}_n}{dt} +  \omega_n^2 \mathbf{P}_n = \sigma_n(\mathbf{x}) \mathbf{E}

That is, we must store and evolve a set of auxiliary fields \mathbf{P}_n along with the electric field in order to keep track of the polarization P. Essentially any ε(ω) could be modeled by including enough of these polarization fields — Meep allows you to specify any number of these, limited only by computer memory and time (which must increase with the number of polarization terms you require).

Loss and gain

If γ above is nonzero, then the dielectric function ε(ω) becomes complex, where the imaginary part is associated with absorption loss in the material if it is positive, or gain if it is negative.

If you look at Maxwell's equations, then d\mathbf{P}/dt plays exactly the same role as a current \mathbf{J}. Just as \mathbf{J} \cdot \mathbf{E} is the rate of change of mechanical energy (the power expended by the electric field on moving the currents), therefore, the rate at which energy is lost to absorption is given by:

absorption rate \sim \frac{d\mathbf{P}}{dt} \cdot \mathbf{E}

Meep can keep track of this energy, which for gain gives the amount of energy expended in amplifying the field. Using this energy, Meep supports the idea of a saturable gain (e.g. a situation in which there is a depletable population inversion causing the gain). For more information, see saturable gain in Meep.

Nonlinearity

In general, ε can be changed anisotropically by the E field itself, with:

\Delta\varepsilon_{ij} = \sum_{k} \chi_{ijk}^{(1)} E_k + \sum_{k\ell} \chi_{ijk\ell}^{(2)} E_k E_\ell  + \cdots

where the ij is the index of the change in the 3×3 ε tensor and the χ terms are the nonlinear susceptibilities. The χ(1) sum is the Pockels effect and the χ(2) sum is the Kerr effect. (If the above expansion is frequency-independent, then the nonlinearity is instantaneous; more generally, Δε would depend on some average of the fields at previous times.)

Currently, Meep supports instantaneous, isotropic Kerr nonlinearities, corresponding to a frequency-independent \chi_{ijk\ell}^{(2)} = \chi^{(2)} \cdot \delta_{ij} \delta_{k\ell}. Thus,

\mathbf{D} = \left( \varepsilon_\infty(\mathbf{x}) + \chi^{(2)}(\mathbf{x}) \cdot |\mathbf{E}|^2 \right) \mathbf{E} + \mathbf{P}

Normally, for nonlinear systems you will want to use real fields E. (This is normally the default. However, Meep uses complex fields if you have Bloch-periodic boundary conditions with a non-zero Bloch wavevector k, or in cylindrical coordinates. In the C++ interface, real fields must be explicitly specified.) For complex fields in nonlinear systems, the physical interpretration of the above equations is more subtle because one cannot simply obtain the physical solution by taking the real part any more.

Personal tools