New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Bessel beams #209
Comments
/cc @stefaniagl |
Proof-of-principle is available at https://github.com/stefaniagl/adda |
6 Bessel beam types have been implemented - LE - linear electric field, LM - linear magnetic field, CS - circularly symmetric energy density, TEC, TMC - polarizations generating TE and TM Bessel beams, and M - generalized type. |
After the pull request (#304) will be merged we still need the following to finalize this issue (to be submitted as the new pull request):
|
- adds vector Bessel beams of LE, LM, CS, CS', TEM, TEL, and generalized types as incident field (#209). They rely on new Fortran file `bessel.f90`. - new command line option `-beam_center <x> <y> <z>`. Corresponding arguments to some older `-beam ...` are now obsolete (#245). - beam center in laboratory reference frame is now specified always (even when at origin) - calculation of Cext through F(0) is now used only for default position of the plane wave - adds beam vorticity as a general parameter. It affects the relation between two incident polarizations (rotation + phase change) and adds phase factor to computation of the amplitude matrix. The latter is ignored in some places, when only the Mueller matrix is computed at the end. - removes a lot of redundant code in calculation of plane incident wave in surface mode (probably was used long ago for alternative calculation of Cext) - adds auxiliary function ScanDouble3Error() in param.c - updates instructions for adding new beams Co-authored-by: stefanyagl <42984899+stefanyagl@users.noreply.github.com> Co-authored-by: Maxim Yurkin <yurkin@gmail.com>
Non-trivial application of Bessel beams for imitating circular polarization of a plane wave are discussed at https://groups.google.com/g/adda-discuss/c/k-ClVwvwQQk . In particular, interesting combinations are used in tests/equiv/ext_CD.py. One of them highlights the lack of robustness for M matrices corresponding to TE and TM beams for very small half-cone angles. |
Original issue reported on code.google.com by
yurkin
on 2 Jul 2015 at 8:59The text was updated successfully, but these errors were encountered: