DOI: 10.5281/zenodo.10598256
title: Source Code and Simulation Results: Chiral and directional optical emission from a dipole source coupled to a helical plasmonic antenna

description: Summary

This publication supplements the article "Chiral and directional optical emission from a dipole source coupled to a helical plasmonic antenna" with tabulated data and Matlab code that allows the reproduction of the results. Within the article, the chiral behavior of single and double plasmonic nano antennas made from silver is numerically investigated with a focus on the coupling of a linear polarized dipole as an excitation source to the helix.

Simulation Setup - FEM Simulations

The script "run_wavlengthscan.m" allows to reproduce all simulations of the article. It can be chosen between the single and double helices, by specifying the keys parameter "keys.doppelhelix" where 0 gives a single and 1 a double helix. The number of turns can be specified by choosing "keys.case". The dipol is located within a 20nm thick hBN substrate layer, on glass (BK7). Results of the Purcell enhancement can be plotted using the scripts "display_results_single_helix.m" and "display_results_doublehelix.m" in the folder "results". The far-field plots can be reproduced using the scripts "display_farfiel_polarization_single_helix.m" and "display_farfiel_polarization_double_helix.m" of the folder "FunctionsAndScripts".

The template for the mesh is contained in the folder "generate_grid_file", where the parameters of the helix (for example: radius, tube radius, and pitch height) can be modified.

Within the folder "project3D" all required .jcm files are stored. Copy the "grid.jcm" file with the geometry of interest to this folder to perform simulations.

All required keys parameters for the JCM template files (.jcmt, jcmpt) are set within the functions "set_numerical_parameter.m" and "set_physical_parameters.m", contained in the folder "FunctionsAndScripts". Therein, the function "set_sources.m" specifies the parameters for the dipole excitation, such as the position, and the strength (equivalent to the polarization).

Material

The material data has been taken from the refractiveindex.info database. For silver the data is taken from tabulated data from Johnson and Christy [1] . The dispersion relation for hBN comes from [2] and tabulated data for glass (BK7) from [3]. The MATLAB script "material_properties_plot.m" plots the material fits above the wavelengths of interest. The required tabulated data is given in the folder "material_data".

With 'material_properties_plot.m' the fits to the material data can be reproduced and plotted.

Usage

The .zip folder Helix_FEM contains all data and scripts to reproduce the plots from the 3D FEM simulations.

Requirements



JCMsuite (at least 5.4.0)

MATLAB (tested with version R2023b)


To run the simulations with JCMsuite you must replace corresponding placeholders with a path to your installation of JCMsuite. Free trial licenses are available, please refer to the homepage of JCMwave.

References

[1] P. B. Johnson and R.-W. Christy, “Optical constants of the noble metals,” Phys. Rev. B 6, 4370 (1972).

[2] S.-Y. Lee, T.-Y. Jeong, S. Jung, and K.-J. Yee, “Refractive index dispersion of hexagonal boron nitride in the visible and near-infrared,” Phys. Status Solidi B 256,  1800417 (2019).[3] “SCHOTT Zemax catalogue 2017-01-20b,” (2017).
description: Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany ́s Excellence Strategy – The Berlin Mathematics Research Center MATH+ (EXC-2046/1, project ID: 390685689) and also by DFG under project ID HO5461/3-1 "chiralFEBID". This project (20FUN05 SEQUME) has received funding from the EMPIR programme co-financed by the Participating States and from the European Union’s Horizon 2020 research and innovation programme. This project has received funding from the German Federal Ministry of Education and Research (BMBF Forschungscampus MODAL, project number 05M20ZBM).