Main Page

Computational Nanophotonics Wiki

Welcome to the Computational Nanophotonics Wiki! This Wiki is only accessible by those involved with the Computational Nanophotonics project headed by Stephen Gray (Chemistry, ANL). A wiki is a website that can be edited and expanded by many people. On this site we can maintain information about research issues, software architecture, who is doing what, and relevant references.

The official project page: http://www.mcs.anl.gov/Nanophotonics/.

Presentations and Relevant Publications

Bibliography

Papers

• Stephen K. Gray and Teobald Kupka, "Propagation of light in metallic nanowire arrays: Finite-difference time-domain studies of silver cylinders", PHYSICAL REVIEW B 68, 045415 (2003).

• M. S. Min, Q. Y. Chen, and Y. Maday, "Spectral method for 2D photonic band structures", Proc. of SPIE, Photonic Crystal Materials and Device II, Vol.5360, pp. 44-51 (2004).

• T.-W. Lee and S. K. Gray, "Controlled spatiotemporal excitation of metal nanoparticles with picosecond optical pulses",PHYSICAL REVIEW B 71, 035423 (2005).

• T.-W. Lee and S. K. Gray, "Subwavelength light bending by metal slit structures", Opt. Express 13, 9652 (2005).

• M. S. Min, "Discontinuous Galerkin method based on quadrilateral mesh for Maxwell's equations", Proc. of ACES/IEEE on Wireless Communications and Applied Computational Electromagnetics, April (2005).

• M.S. Min, T.-W. Lee, P.F. Fischer and S.K. Gray, "Fourier Spectral Simulations and Gegenbauer Reconstructions for Electromagnetic Waves in the Presence of a Metal Nanoparticle", Journal of Computational Physics, 213 (2) 730-747 (2006).

• Z. Meglicki, S. K. Gray and B. Norris, "Multigrid FDTD with Chombo", Computer Physics Communications, 172(2), 109-120, 15 January 2007.

Books

• Allen Taflove and Susan C. Hagness, Computational Electrodynamics: The Finite Difference Time-Domain Method, second edition, Artech House, Inc., Boston, 2000. (ISBN 1-58053-076-1)

• Dennis M. Sullivan, Electromagnetic Simulation Using the FDTD Method, IEEE Press, New York, 2000. (ISBN 0-7803-4747-1)

• Guiseppe Pelosi, Roberto Coccioli and Stefano Selleri, Quick Finite Elements for Electromagnetic Waves, Artech House, Inc., Boston, 1998. (ISBN 0-89006-848-8)

Reports

• Dan Martin and Gustav Meglicki, "Chombo Class PiecewiseLinearFillPatchFace", a code weave, August 22, 2006
• Dan Martin and Gustav Meglicki, "Chombo Class LevelFluxRegisterEdge", a code weave, November 29, 2006
• Gustav Meglicki, "Forms", code weave, last version
  • "Injected Signal", Gnuplot of signal injected in the (0,0,1) direction. The signal is (1,-1,0) polarized.
  • "Forms input file" for the simulation.
  • "Isosurface animation", the above signal scattering on a glazed metal ball.
  • "Cross section animation", same simulation, (y,z) plane cross-section.
  • "Forward scattering, Poynting flux", same simulation, forward scattered Poynting flux through the (x,y) plane in the scattered field region--animation.
  • "Back-scattered Poynting flux", same simulation, back-scattered Poynting flux through the (x,y) plane in the scattered field region--animation.

Presentations

• M.S. Min, Q. Y. Chen, Y. Maday, "Spectral methods for 2D photonic band structures", Photonic Band gap Materials and Devices II (OE12) as a part of SPIE's International Symposium on Integrated Optoelectronic Devices, San Jose, CA, Jan. 26-29, 2004.

• M. S. Min, T. W. Lee, P. F. Fischer, S. K. Gray, "Spectral methods for electromagnetic waves in discontinuous media", The 6th International Conference on Spectral and High-Order Methods (ICOSAHOM), Recent Advances in Spectral Methods and High Order Finite Difference Schemes, Providence, RI, Jun. 21-25, 2004.

• M.S. Min, P. F. Fischer, S. K. Gray, "Spectral element discontinuous Galerkin simulations for photonic crystals", Photonic Band Gap Materials and Devices II (OE12) as a part of SPIE's International Symposium on Integrated Optoelectronic Devices, San Jose, CA, Jan. 22-26, 2006.

• M.S. Min, P. F. Fischer, S. K. Gray, "Spectral element discontinuous Galerkin simulations for metallic nanoparticles", PIERS 2006 in Cambridge - Progress in Electromagnetics Research Symposium, Cambridge, MA, Mar 26-29, 2006.

• M.S. Min, "Spectral element discontinuous Galerkin simulations for waveguiding structures", PIERS 2006 in Tokyo - Progress in Electromagnetics Research Symposium, Tokyo, Japan, Aug 2-5, 2006.

• M.S. Min, P. F. Fischer, S. K. Gray, "Spectral-element discontinuous Galerkin simulations with exponential time integrator for nanophotonics: NEKCEM software", International Conference On Spectral and High Order Methods (ICOSAHOM), Beijing, China, Jun. 18-22, 2007.

Programs and Simulations

3D FDTD - "FORMS"

The Forms code solves Maxwell equations in vacuum and in the Lorentz media in three dimensions using a combination of Finite Difference Time Domain (FDTD) and static multigrid methods. The program uses Fortran stubs for basic FDTD computations, Chombo C++ wrapper for general logic and multigrid operations, and Guile, which is a dialect of Scheme, for multigrid and media layout specifications, for media properties specifications, for injected signal specifications, for some output specifications, and for steering the iteration process.

The standard Chombo library is additionally enhanced by new templates and classes that implement multigrid logic for face and edge mounted data in parallel computing context. These borrow from some functions of the staggeredChombo package mercifully contributed to the project by Dan Martin of the Berkeley National Lab. The staggeredChombo classes implement divergence-free averaging from fine to coarse grids, curl-refluxing corrections on the coarse side of the fine/coarse grid boundaries, and piecewise linear interpolation filling of cells on the fine side of the fine/coarse grid boundaries for face-mounted data.

• "FORMS User Manual" Code weave with introduction, user and administrator guides, and discussion of numerics, PDF, 552 pages (the user guide portion is about 200 pages), about 3MB--June 21, 2010;
• "Ex_yz_movie.gif" Animation: harmonic wave scattering on a dielectric ball--generated numerically with Forms, E_x data for the (x = x_0, y, z) slice;
• "mie_yz_movie.gif" Animation: harmonic wave scattering on a dielectric ball--generated by the Mie analytical solution, E_x data for the (x = x_0, y, z) slice;
• "Ex_xz_movie.gif" Animation: harmonic wave scattering on a dielectric ball--generated numerically with Forms, E_x data for the (x, y = y_0, z) slice;
• "mie_xz_movie.gif" Animation: harmonic wave scattering on a dielectric ball--generated by the Mie analytical solution, E_x data for the (x, y = y_0, z) slice;
• "Input.scm" An input file used to produce the above Forms animation;
• "ball-test.scm" A Guile script used to produce the Mie animation;
• "mie-movie for n=1" Animation: harmonic wave scattering on a dielectric ball with n=1--generated by the Mie analytical solution, E_x data for the (x = x_0, y, z) slice. This animation tests the correctness of the Mie code itself: should be a plane harmonic wave throughout the whole total field region.
• "FDTD scattering on a silver ball at 0.89eV, zy slice" GIF Animation, 180^3 resolution, 6.5MB
• "FDTD scattering on a silver ball at 0.89eV, zy slice" GIF Animation, 256^3 resolution, 9.8MB
• "Drude-Input.scm" A Guile script used to produce the above on the IU's Quarry cluster (x86_64 GNU/Linux 8-way SMP nodes)
• "Parallel-Drude-Input.scm" A Guile script to produce the above on the IU's Quarry cluster (x86_64 GNU/Linux 8-way SMP nodes), running on 8 cores under MPI
• "Analytical scattering on a silver ball at 0.89eV, zy slice" GIF Animation, 180^3 resolution, 12.5MB
• "Mie-Drude-Input.scm" A Guile script used to produce the above under CYGWIN_NT-6.0-WOW64 on a 4-core i686 desktop
• "FDTD scattering on a silver ball at 0.89eV, zx slice" GIF Animation, 180^3 resolution, 6.0MB
• "FDTD scattering on a silver ball at 0.89eV, zx slice" GIF Animation, 256^3 resolution, 9.0MB
• "Analytical scattering on a silver ball at 0.89eV, zx slice" GIF Animation, 180^3 resolution, 11.3MB

AMR-Enhanced 2D FDTD - "Shapes"

• "Shapes" Documentation:
  • shapes(1) - UNIX style man page that briefly describes the program.
  • shapes(5) - UNIX style man page that discusses in gruelling detail what should go into the "Shapes" input file.
  • "Shapes-2.1" User Guide This document explains in great detail how "Shapes" works and how to use it. Numerous hands-on illustrations are provided.
  • An animated GIF example This is an animation discussed in the "Shapes" User Guide, section 5, page 73, "Working with Multigrid"
  • An animated GIF example Same example, but this time level 1 propagation is shown.
  • An animated GIF example Same example, level 2.
• Visualizing Shapes Datasets How to build and use visualization software (email norris@mcs.anl.gov to obtain a user login for accessing this and other Wiki pages).

NekCEM - Spectral element discontinuous Galerkin code

• NekCEM instructions:
  • Quick Start this gives short instructions how NekCEM works and how to use it.
• Spectral element discontinuous Galerkin (SEDG) simulations:
  • An animated GIF example nanocylinder: scattered field energy profiles.
  • An animated GIF example nanocylinder: scattered field amplitude for an electric field component with spectral element mesh.
  • An animated GIF example nanocylinder: scattered field amplitude and PML boundary for electric field component with contour lines.
• Spectral element meshes:
  • 3D mesh (top view) for funnel arrays of nanospheres.
  • 3D simulations: electric field enhancement around the surface of nanospheres (2d cross-section view).
  • 3D plasmonic crystal structure.
  • 3D isolated hole/disk structure.

Internal

• Intraproject Pages