We present an inverse mapping approach to determining the emission height of the optical photons from pulsars, which is directly constrained by empirical data. The model discussed is for the case of the Crab pulsar. Our method, using the optical Stokes parameters, determines the most likely geometry for emission including magnetic field inclination angle ($\alpha$), observers line of sight angle ($\chi$) and emission height. We discuss the computational implementation of the approach, along with any physical assumptions made. We find that the most likely emission altitude is at 20% of the light cylinder radius above the stellar surface, in the open field region. We also present a general treatment of the expected polarisation from synchrotron source with a truncated power law spectrum of particles.

This work presents a highly optimized computational framework for the Discrete Dipole Approximation, a numerical method for calculating the optical properties associated with a target of arbitrary geometry that is widely used in atmospheric, astrophysical and industrial simulations. Core optimizations include the bit-fielding of integer data and iterative methods that complement a new Discrete Fourier Transform (DFT) kernel, which efficiently calculates the matrix vector products required by these iterative solution schemes. The new kernel performs the requisite 3-D DFTs as ensembles of 1-D transforms, and by doing so, is able to reduce the number of constituent 1-D transforms by 60% and the memory by over 80%. The optimizations also facilitate the use of parallel techniques to further enhance the performance. Complete OpenMP-based shared-memory and MPI-based distributed-memory implementations have been created to take full advantage of the various architectures. Several benchmarks of the new framework indicate extremely favorable performance and scalability. OpenDDA is available following the usual open source regulations from http://www.opendda.org

The magnetospheric plasma distribution in the vicinity of a pulsar at various inclination angles is investigated using a new relativistic, parallel 3D Particle-In-Cell (PIC) code DYMPHNA3D. DYMPHNA3D uses a superposition of the electromagnetic fields associated with a rotating magnetised conducting sphere in a vacuum (the pulsar fields) and the electromagnetic fields due to the presence of the magnetospheric plasma surrounding the pulsar (the plasma fields), as the total fields. The plasma is moved self-consistently through the magnetosphere using the PIC methodology. Our initial simulation results are presented here. These show similar solutions to those obtained from previous numerical simulations, which show the fundamental plasma distribution in the vicinity of an aligned rotating neutron star to consist of two polar domes and an equatorial torus of trapped non-neutral plasma of opposite sign. The aligned case being the case in which the inclination angle between the magnetic dipole moment and the rotation axis of the star is zero. Furthermore, our code allows for off-axis simulations and we have found that this plasma distribution collapses into a Quad-Lobe charge-separated non-neutral magnetospheric plasma distribution in the case of an orthogonal rotator, i.e., the case in which the magnetic dipole moment is oriented at right angles to the rotation axis of the neutron star, with the plasma remaining trapped close to the stellar surface by the force-free surfaces. We find that if initialised with a Goldreich-Julian type distribution, the system is seen to collapse rapidly into these stable Dome-Torus structures.