The vast majority of pulsars detected by the Fermi Large Area Telescope (LAT) display exponentially cutoff spectra with cutoffs falling in a narrow band around a few GeV. Early spectral modelling predicted spectral cutoffs at energies of up to 100 GeV, assuming curvature radiation. It was therefore not expected that pulsars would be visible in the very-high energy (VHE) regime (>100 GeV). The VERITAS announcement of the detection of pulsed emission from the Crab pulsar at energies up to 400 GeV (and now up to 1.5 TeV as detected by MAGIC) therefore raised important questions about our understanding of the electrodynamics and local environment of pulsars. H.E.S.S. has now detected pulsed emission from the Vela pulsar down to tens of GeV, making this the second pulsar detected by a ground-based Cherenkov telescope. Deep upper limits have also been obtained by VERITAS and MAGIC for the Geminga pulsar. We will review the latest developments in VHE pulsar science, including an overview of the latest observations, refinements, and extensions to radiation models and magnetic field structures, and the implementation of new radiation mechanisms. This will assist us in understanding the VHE emission detected from the Crab pulsar, and predicting the level of VHE emission expected from other pulsars, which is very important for the upcoming CTA.

The vast majority of the pulsars detected by the Fermi Large Area Telescope (LAT) display spectra with exponential cutoffs falling in a narrow range around a few GeV. Early spectral modelling predicted spectral cutoff energies of up to 100 GeV. More modern studies estimated spectral cutoff energies in the 1-20 GeV range. It was therefore not expected that pulsars would be visible in the very-high-energy (VHE; >100 GeV) regime. The VERITAS detection (confirmed by MAGIC) of pulsed emission from the Crab pulsar up to 400 GeV (and now possibly up to 1 TeV) therefore raised important questions about our understanding of the electrodynamics and local environment of pulsars. H.E.S.S. has now detected pulsed emission from the Vela pulsar in the 20-120 GeV range, making this the second pulsar detected by a ground-based Cherenkov telescope. We will review the latest developments in VHE pulsar science, including an overview of recent observations and refinements to radiation models and magnetic field structures. This will assist us in interpreting the VHE emission detected from the Crab and Vela pulsars, and predicting the level of VHE emission expected from other pulsars, which will be very important for the upcoming CTA.

Over the past six years, the Fermi Large Area Telescope has detected more than 150 gamma-ray pulsars, discovering a variety of light curve trends and classes. Such diversity hints at distinct underlying magnetospheric and/or emission geometries. We implemented an offset-dipole magnetic field, with an offset characterised by parameters epsilon and magnetic azimuthal angle, in an existing geometric pulsar modelling code which already includes static and retarded vacuum dipole fields. We use these different magnetic field solutions in conjunction with standard emission geometries, namely the two-pole caustic and outer gap models (the latter only for non-offset dipoles), and construct intensity maps and light curves for several pulsar parameters. We compare our model light curves to the Vela data from the second pulsar catalogue of Fermi. We use a refined chi-square grid search method for finding best-fit light curves for each of the different models. Our best fit is for the retarded vacuum dipole field and the outer gap model.

The second pulsar catalogue of the Fermi Large Area Telescope (LAT) will contain in excess of 100 gamma-ray pulsars. The light curves (LCs) of these pulsars exhibit a variety of shapes, and also different relative phase lags with respect to their radio pulses, hinting at distinct underlying emission properties (e.g., inclination and observer angles) for the individual pulsars. Detailed geometric modelling of the radio and gamma-ray LCs may provide constraints on the B-field structure and emission geometry. We used different B-field solutions, including the static vacuum dipole and the retarded vacuum dipole, in conjunction with an existing geometric modelling code, and constructed radiation sky maps and LCs for several different pulsar parameters. Standard emission geometries were assumed, namely the two-pole caustic (TPC) and outer gap (OG) models. The sky maps and LCs of the various B-field and radiation model combinations were compared to study their effect on the resulting LCs. As an application, we compared our model LCs with Fermi LAT data for the Vela pulsar, and inferred the most probable configuration in this case, thereby constraining Vela's high-altitude magnetic structure and system geometry.

The light curves of gamma-ray pulsars detected by the Fermi Large Area Telescope show great variety in profile shape and position relative to their radio profiles. Such diversity hints at distinct underlying magnetospheric and/or emission geometries for the individual pulsars. We implemented an offset-dipole magnetic field in an existing geometric pulsar modelling code which already includes static and retarded vacuum dipole fields. In our model, this offset is characterised by a parameter epsilon (with epsilon = 0 corresponding to the static dipole case). We constructed sky maps and light curves for several pulsar parameters and magnetic fields, studying the effect of an offset dipole on the resulting light curves. A standard two-pole caustic emission geometry was used. As an application, we compared our model light curves with Fermi data for the bright Vela pulsar.