Inkstone simulates the electromagnetic properties of 3D and 2D multi-layered structures with in-plane periodicity, such as gratings, photonic-crystal slabs, metasurfaces, vertical-cavity or photonic-crystal surface-emitting lasers (VCSEL, PCSEL), (patterned) solar cells, nano-antennas, and more.
Internally, Inkstone implements rigorous coupled-wave analysis (RCWA), a. k. a. Fourier Modal Method (FMM).
Inkstone can calculate:
- the reflection, transmission, and absorption of the structure
- the total and by-order power fluxes of the propagating and the evanescent waves in each layer
- electric and magnetic field amplitudes at any locations in the structure,
- band-structures based on the determinant of the scattering matrix of the structure.
Features of Inkstone:
- It supports efficient and flexible parameter-scanning. You can change part of your structure such as the shapes and sizes of some patterns, or some material parameters. Inkstone only recalculates the modified parts and produces the final results efficiently.
- It allows both tensorial permittivities and tensorial permeabilities, such as in anisotropic, magneto-optical, or gyromagnetic materials.
- It can calculate the determinant of the scattering matrix on the complex frequency plane.
- Pre-defined shapes of patterns can be used, including rectangular, parallelogram, disk, ellipse, 1D, and polygons. Closed-form Fourier transforms and corrections for Gibbs phenomena are implemented.
- It is fully 3D.
- It is written in pure python, with heavy-lifting done in numpy and scipy.
Quick Start
Installation:
$ pip install inkstone
Or,
$ git clone git://github.com/alexysong/inkstone
$ pip install .
Usage
The examples folder contains various self-explaining examples to get you started.
Dependencies
- python 3.6+
- numpy
- scipy
Units, conventions, and definitions
Unit system
We adopt a natural unit system, where vacuum permittivity, permeability, and light speed are $\varepsilon_0=\mu_0=c_0=1$.
Sign convention
Sign conventions in electromagnetic waves:
$$e^{i(kx-\omega t)}$$
where $k$ is the wavevector, $x$ is spatial location, $\omega$ is frequency, $t$ is time.
By this convention, a permittivity of $\varepsilon_r + i\varepsilon_i$ with $\varepsilon_i>0$ means material loss, and $\varepsilon_i<0$ means material gain.
Coordinates and incident angles
(Inkstone, Incident $\bm{k}$ on stacked periodic nano electromagnetic structures.)
Citing
If you find Inkstone useful for your research, we would apprecite you citing our paper. For your convenience, you can use the following BibTex entry:
@article{song2018broadband,
title={Broadband Control of Topological Nodes in Electromagnetic Fields},
author={Song, Alex Y and Catrysse, Peter B and Fan, Shanhui},
journal={Physical review letters},
volume={120},
number={19},
pages={193903},
year={2018},
publisher={American Physical Society}
}
40 Dec 22, 2022
112 Jan 2, 2023
![[ICCV'21] PlaneTR: Structure-Guided Transformers for 3D Plane Recovery](https://github.com/IceTTTb/PlaneTR3D/raw/master/misc/network.jpg)
73 Nov 30, 2022
258 Jan 2, 2023
353 Dec 27, 2022
63 Jan 6, 2023
77 Jan 5, 2023
6 Jul 7, 2022
38 Oct 11, 2022
86 Jan 04, 2023
0 Jan 02, 2022
568 Jan 04, 2023
0 Feb 06, 2022
1.1k Jan 02, 2023
2 May 14, 2022
120 Dec 12, 2022
0 Dec 12, 2021
62 Dec 21, 2022
1.8k Jan 08, 2023
32 Oct 24, 2022
23 Mar 09, 2022
10 Dec 20, 2022
110 Dec 06, 2022
0 Dec 18, 2021
7 Sep 23, 2022
49 Dec 10, 2022
4 Oct 20, 2022
328 Dec 29, 2022
26 Dec 07, 2022