We show that the $\gamma$-ray pulsar observables, i.e., their total $\gamma$-ray luminosity, $L_{\gamma}$, spectral cut-off energy, $\epsilon_{\rm cut}$, stellar surface magnetic field, $B_{\star}$, and spin-down power $\dot{\mathcal{E}}$, obey a relation of the form $L_{\gamma}=f(\epsilon_{\rm cut},B_{\star},\dot{\mathcal{E}})$, which represents a 3D plane in their 4D log-space. Fitting the data of 88 pulsars of the second Fermi pulsar catalog, we show this relation to be $L_{\gamma}\propto \epsilon_{\rm cut}^{1.18\pm 0.24}B_{\star}^{0.17\pm 0.05}\dot{\mathcal{E}}^{0.41\pm 0.08}$, a pulsar fundamental plane (FP). We show that the observed FP is remarkably close to the theoretical relation $L_{\gamma}\propto \epsilon_{\rm cut}^{4/3}B_{\star}^{1/6}\dot{\mathcal{E}}^{5/12}$ obtained assuming that the pulsar $\gamma$-ray emission is due to curvature radiation by particles accelerated at the pulsar equatorial current sheet just outside the light cylinder. Interestingly, the FP seems incompatible with emission by synchrotron radiation. The corresponding scatter about the FP is $\sim 0.35$dex and can only partly be explained by the observational errors while the rest is probably due to the variation of the inclination and observer angles. We predict also that $\epsilon_{\rm cut}\propto \dot{\mathcal{E}}^{7/16}$ toward low $\dot{\mathcal{E}}$ for both young and millisecond pulsars implying that the observed death-line of $\gamma$-ray pulsars is due to $\epsilon_{\rm cut}$ dropping below the Fermi-band. Our results provide a comprehensive interpretation of the observations of $\gamma$-ray pulsars, setting requirement for successful theoretical modeling.

A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics.

Pulsars emit pulsed emission across the entire electromagnetic spectrum and their light curve phenomenology is strongly dependent on energy. This is also true for the gamma-ray waveband. Continued detections by Fermi Large Area Telescope in the GeV band and ground-based Cherenkov telescopes in the TeV band (e.g., Crab and Vela above 1 TeV) raise important questions about our understanding of the electrodynamics and local environment of pulsar magnetospheres. We model energy-dependent light curves (as a function of geometry, e.g., pulsar inclination and observer angle) in the curvature radiation domain using a full emission code. We will discuss our refined calculation of the curvature radius of the particle trajectory and the effect thereof on the expected light curve shapes, as well as the origin of the light curve peaks in the magnetosphere. Our modelling should aid in differentiating between different emission mechanisms, as well as constraining the emission geometry by comparing our predictions to multi-wavelength data.

Pulsed emission from the Vela pulsar at energies above 3 TeV has recently been detected by the H.E.S.S. II air-Cherenkov telescope. We present a model for the broad-band spectrum of Vela from infra-red (IR) to beyond 10 TeV. Recent simulations of the global pulsar magnetosphere have shown that most of the particle acceleration occurs in the equatorial current sheet outside the light cylinder and that the magnetic field structure is nearly force-free for younger pulsars. We adopt this picture to compute the radiation from both electron-positron pairs produced in polar cap cascades and from primary particles accelerated in the separatrix and current sheet. The synchrotron spectrum from pairs resonantly absorbing radio photons at relatively low altitude can account for the observed IR-optical emission. We set the parallel electric field in the current sheet to produce the Fermi GeV emission through curvature radiation, producing particles with energies of 30-60 TeV. These particles then produce Very-High-Energy emission up to around 30 TeV through inverse-Compton scattering of the IR-Optical emission. We present model spectra and light curves that can match the IR-Optical, GeV and make predictions for the multi-TeV emission.

We present 3D global kinetic pulsar magnetosphere models, where the charged particle trajectories and the corresponding electromagnetic fields are treated self-consistently. For our study, we have developed a cartesian 3D relativistic particle-in-cell code that incorporates the radiation reaction forces. We describe our code and discuss the related technical issues, treatments, and assumptions. Injecting particles up to large distances in the magnetosphere, we apply arbitrarily low to high particle injection rates and get an entire spectrum of solutions from close to the Vacuum-Retarded-Dipole to close to the Force-Free solution, respectively. For high particle injection rates (close to FF solutions) significant accelerating electric field components are confined only near the equatorial current sheet outside the light-cylinder. A judicious interpretation of our models allows the calculation of the particle emission and consequently the derivation of the corresponding realistic high-energy sky-maps and spectra. Using model parameters that cover the entire range of spin-down powers of Fermi young and millisecond pulsars, we compare the corresponding model $\gamma$-ray light-curves, cutoff energies, and total $\gamma$-ray luminosities with those observed by Fermi to discover a dependence of the particle injection-rate, $\mathcal{F}$, on the spin-down power, $\dot{\mathcal{E}}$, indicating an increase of $\mathcal{F}$ with $\dot{\mathcal{E}}$. Our models guided by Fermi observations provide field-structures and particle distributions that are not only consistent with each other but also able to reproduce a broad range of the observed $\gamma$-ray phenomenology of both young and millisecond pulsars.

We performed ab-initio Particle-In-Cell (PIC) simulations of a pulsar magnetosphere with electron-positron plasma produced only in the regions close to the neutron star surface. We study how the magnetosphere transitions from the vacuum to a nearly force-free configuration. We compare the resulting force-free like configuration with ones obtained in a PIC simulation where particles are injected everywhere as well as with macroscopic force-free simulations. We found that although both PIC solutions have similar structure of electromagnetic fields and current density distributions, they have different particle density distribution. In fact in the injection from the surface solution, electrons and positrons counterstream only along parts of the return current regions and most of the particles leave the magnetosphere without returning to the star. We also found that pair production in the outer magnetosphere is not critical for filling the whole magnetosphere with plasma. We study how the current density distribution supporting the global electromagnetic configuration is formed by analyzing particle trajectories. We found that electrons precipitate to the return current layer inside the light cylinder and positrons precipitate to the current sheet outside the light cylinder by crossing magnetic field lines contributing to the charge density distribution required by the global electrodynamics. Moreover, there is a population of electrons trapped in the region close to the Y-point. On the other hand the most energetic positrons are accelerated close to the Y-point. These processes can have observational signatures that, with further modeling efforts, would help to distinguish this particular magnetosphere configuration from others.

A sub-population of energetic rotation-powered pulsars show high fluxes of pulsed non-thermal hard X-ray emission. While this MeV pulsar population includes some radio-loud pulsars like the Crab, a significant number have no detected radio or GeV emission, a mystery since gamma- ray emission is a common characteristic of pulsars with high spin-down power. Their steeply rising hard X-ray spectral energy distributions (SEDs) suggest peaks at 0.1 - 1 MeV but they have not been detected above 200 keV. Several upcoming and planned telescopes may shed light on the MeV pulsars. The Neutron star Interior Composition ExploreR (NICER) will observe pulsars in the 0.2 - 12 keV band and may discover additional MeV pulsars. Planned telescopes, such as All-Sky Medium-Energy Gamma-Ray Observatory (AMEGO) and e-ASTROGAM, will detect emission above 0.2 MeV and polarization in the 0.2 - 10 MeV band. We present a model for the spectrum and polarization of MeV pulsars where the X-ray emission comes from electron- positron pairs radiating in the outer magnetosphere and current sheet. This model predicts that the peak of the SED increases with surface magnetic field strength if the pairs are produced in polar cap cascades. For small inclination angles, a range of viewing angles can miss both the radio pulse and the GeV pulse from particles accelerating near the current sheet. Characterizing the emission and geometry of MeV pulsars can thus provide clues to the source of pairs and acceleration in the magnetosphere.

In recent years, surprise discoveries of pulsed emission from the Crab and Vela pulsars above 100 GeV have drawn renewed attention to this largely unexplored region of the energy range. In this paper, we discuss example light curves due to curvature emission, with good resolution in the different energy bands. Continued light curve modelling may help to discriminate between different emission mechanisms, as well as constrain the location where emission is produced within the pulsar magnetosphere, including regions beyond the light cylinder.

Based on the Fermi observational data we reveal meaningful constraints for the dependence of the macroscopic conductivity $(\sigma)$ of dissipative pulsar magnetosphere models on the corresponding spin-down rate, $\dot{\mathcal{E}}$. Our models are refinements of the FIDO (Force-Free Inside, Dissipative Outside) models whose dissipative regions are restricted on the equatorial current-sheet outside the light-cylinder. Taking into account the observed cutoff-energies of all the Fermi-pulsars and assuming that a) the corresponding $\gamma-$ray pulsed emission is due to curvature radiation at the radiation-reaction-limit regime and b) this emission is produced at the equatorial current-sheet near the light-cylinder, we show that the \emph{Fermi}-data provide clear indications about the corresponding accelerating electric-field components. A direct comparison between the \emph{Fermi} cutoff-energies and the model ones reveals that $\sigma$ increases with $\dot{\mathcal{E}}$ for high $\dot{\mathcal{E}}$-values while it saturates for low ones. This comparison indicates also that the corresponding gap-width increases toward low $\dot{\mathcal{E}}$-values. Assuming the Goldreich-Julian flux for the emitting particles we calculate the total $\gamma-$ray luminosity $(L_{\gamma})$. A comparison between the dependence of the Fermi $L_{\gamma}$-values and the model ones on $\dot{\mathcal{E}}$ indicates an increase of the emitting particle multiplicity with $\dot{\mathcal{E}}$. Our modeling guided by the \emph{Fermi}-data alone, enhances our understanding of the physical mechanisms behind the high energy emission in pulsar magnetospheres.

Polarization measurements provide strong constraints on models for emission from rotation-powered pulsars. We present multiwavelength polarization predictions showing that measurements over a range of frequencies can be particularly important for constraining the emission location, radiation mechanisms and system geometry. The results assume a generic model for emission from the outer magnetosphere and current sheet in which optical to hard X-ray emission is produced by synchrotron radiation from electron-positron pairs and gamma-ray emission is produced by curvature radiation or synchrotron radiation from accelerating primary electrons. The magnetic field structure of a force-free magnetosphere is assumed and the phase-resolved and phase-averaged polarization is calculated in the frame of an inertial observer. We find that large position angle swings and deep depolarization dips occur during the light curve peaks in all energy bands. For synchrotron emission, the polarization characteristics are strongly dependent on photon emission radius with larger, nearly $180^\circ$, position angle swings for emission outside the light cylinder as the line-of-sight crosses the current sheet. The phase-averaged polarization degree for synchrotron radiation is less that 10% and around 20% for emission starting inside and outside the light cylinder respectively, while the polarization degree for curvature radiation is much larger, up to 40% - 60%. Observing a sharp increase in polarization degree and a change in position angle at the transition between X-ray and gamma-ray spectral components would indicate that curvature radiation is the gamma-ray emission mechanism.

Results of a simulation of synchrotron-self Compton (SSC) emission from a rotation-powered pulsar are presented. The radiating particles are assumed to be both accelerated primary electrons and a spectrum of electron-positron pairs produced in cascades near the polar cap. They follow trajectories in a slot gap using 3D force-free magnetic field geometry, gaining pitch angles through resonant cyclotron absorption of radio photons, radiating and scattering synchrotron emission at high altitudes out to and beyond the light cylinder. Full angular dependence of the synchrotron photon density is simulated in the scattering and all processes are treated in the inertial observer frame. Spectra for the Crab and Vela pulsars as well as two energetic millisecond pulsars, B1821-24 and B1937+21 are simulated using this model. The simulation of the Crab pulsar radiation can reproduce both the flux level and the shape of the observed optical to hard X-ray emission assuming a pair multiplicity of $M_+ = 3 \times 10^5$, as well as the very-high-energy emission above 50 GeV detected by MAGIC and VERITAS, with both the synchrotron and SSC components reflecting the shape of the pair spectrum. Simulations of Vela, B1821$-$24 and B1937+21, for $M_+$ up to $10^5$, do not produce pair SSC emission that is detectable by current telescopes, indicating that only Crab-like pulsars produce significant SSC components. The pair synchrotron emission matches the observed X-ray spectrum of the millisecond pulsars and the predicted peak of this emission at 1 - 10 MeV would be detectable with planned Compton telescopes.

We explore the emission properties of a dissipative pulsar magnetosphere model introduced by Kalapotharakos et al. (2014), comparing its high energy light curves and spectra, due to curvature radiation, with data collected by the Fermi LAT. The magnetosphere structure is assumed to be near the force-free solution. The accelerating electric field, inside the light-cylinder, is assumed to be negligible, while outside the light-cylinder it rescales with a finite conductivity ({\sigma}). In our approach we calculate the corresponding high energy emission by integrating the trajectories of test particles that originate from the stellar surface, taking into account both the accelerating electric field components and the radiation reaction forces. First we explore the parameter space assuming different value sets for the stellar magnetic field, stellar period, and conductivity. We show that the general properties of the model are in a good agreement with observed emission characteristics of young {\gamma}-ray pulsars, including features of the phase resolved spectra. Second we find model parameters that fit each pulsar belonging to a group of eight bright pulsars that have a published phase-resolved spectrum. The {\sigma} values that best describe each of the pulsars in this group show an increase with the spin-down rate $(\dot{E})$ and a decrease with the pulsar age, expected if pair cascades are providing the magnetospheric conductivity. Finally, we explore the limits of our analysis and suggest future directions for improving such models.

Using both numerical and analytical approaches, we demonstrate the existence of an effective power-law relation $L\propto m^p$ between the mean Lyapunov exponent $L$ of stellar orbits chaotically scattered by a supermassive black hole in the center of a galaxy and the mass parameter $m$, i.e. ratio of the mass of the black hole over the mass of the galaxy. The exponent $p$ is found numerically to obtain values in the range $p \approx 0.3$--$0.5$. We propose a theoretical interpretation of these exponents, based on estimates of local `stretching numbers', i.e. local Lyapunov exponents at successive transits of the orbits through the black hole's sphere of influence. We thus predict $p=2/3-q$ with $q\approx 0.1$--$0.2$. Our basic model refers to elliptical galaxy models with a central core. However, we find numerically that an effective power law scaling of $L$ with $m$ holds also in models with central cusp, beyond a mass scale up to which chaos is dominated by the influence of the cusp itself. We finally show numerically that an analogous law exists also in disc galaxies with rotating bars. In the latter case, chaotic scattering by the black hole affects mainly populations of thick tube-like orbits surrounding some low-order branches of the $x_1$ family of periodic orbits, as well as its bifurcations at low-order resonances, mainly the Inner Lindbland resonance and the 4/1 resonance. Implications of the correlations between $L$ and $m$ to determining the rate of secular evolution of galaxies are discussed.

We compute the patterns of $\gamma$-ray emission due to curvature radiation in dissipative pulsar magnetospheres. Our ultimate goal is to construct macrophysical models that are able to reproduce the observed $\gamma$-ray light-curve phenomenology recently published in the Second Fermi Pulsar Catalog. We apply specific forms of Ohm's law on the open field lines using a broad range for the macroscopic conductivity values that result in solutions ranging, from near-vacuum to near Force-Free. Using these solutions, we generate model $\gamma$-ray light curves by calculating realistic trajectories and Lorentz factors of radiating particles, under the influence of both the accelerating electric fields and curvature radiation-reaction. We further constrain our models using the observed dependence of the phase-lags between the radio and $\gamma$-ray emission on the $\gamma$-ray peak-separation. We perform a statistical comparison of our model radio-lag vs peak-separation diagram and the one obtained for the Fermi standard pulsars. We find that for models of uniform conductivity over the entire open magnetic field line region, agreement with observations favors higher values of this parameter. We find, however, significant improvement in fitting the data with models that employ a hybrid form of conductivity; specifically, infinite conductivity interior to the light-cylinder and high but finite conductivity on the outside. In these models the $\gamma$-ray emission is produced in regions near the equatorial current sheet but modulated by the local physical properties. These models have radio-lags near the observed values and statistically best reproduce the observed light-curve phenomenology. Additionally, these models produce GeV photon cut-off energies.

We investigate the shapes of \gamma-ray pulsar light curves using 3D pulsar magnetosphere models of finite conductivity. These models, covering the entire spectrum of solutions between vacuum and force-free magnetospheres, for the first time afford mapping the GeV emission of more realistic, dissipative pulsar magnetospheres. To this end we generate model light curves following two different approaches: (a) We employ the emission patterns of the slot and outer gap models in the field geometries of magnetospheres with different conductivity \sigma. (b) We define realistic trajectories of radiating particles in magnetospheres of different \sigma and compute their Lorentz factor under the influence of magnetospheric electric fields and curvature radiation-reaction; with these at hand we then calculate the emitted radiation intensity. The light curves resulting from these prescriptions are quite sensitive to the value of \sigma, especially in the second approach. While still not self-consistent, these results are a step forward in understanding the physics of pulsar \gamma-radiation.

We present the magnetic and electric field structures as well as the currents and charge densities of pulsar magnetospheres which do not obey the ideal condition, ${\bf E \cdot B =0}$. Since the acceleration of particles and the production of radiation requires the presence of an electric field component parallel to the magnetic field, ${\bf E}_\parallel$, the structure of non-Ideal pulsar magnetospheres is intimately related to the production of pulsar radiation. Therefore, knowledge of the structure of non-Ideal pulsar magnetospheres is important because their comparison (including models for the production of radiation) with observations will delineate the physics and the parameters underlying the pulsar radiation problem. We implement a variety of prescriptions that support nonzero values for ${\bf E}_\parallel$ and explore their effects on the structure of the resulting magnetospheres. We produce families of solutions that span the entire range between the vacuum and the (ideal) Force-Free Electrodynamic solutions. We also compute the amount of dissipation as a fraction of the Poynting flux for pulsars of different angles between the rotation and magnetic axes and conclude that this is at most 20-40% (depending on the non-ideal prescription) in the aligned rotator and 10% in the perpendicular one. We present also the limiting solutions with the property $J=\rho c$ and discuss their possible implication on the determination of the "on/off" states of the intermittent pulsars. Finally, we find that solutions with values of $J$ greater than those needed to null ${\bf E}_\parallel$ locally produce oscillations, potentially observable in the data.

Recent studies have shown that gamma-ray pulsar light curves are very sensitive to the geometry of the pulsar magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magnetospheres have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve modeling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profiles. We find that, compared to the profiles derived from symmetric caps, the flux in the pulse peaks, which are caustics formed along the trailing magnetic field lines, increases significantly relative to the off-peak emission, formed along leading field lines. The enhanced contrast produces improved slot gap model fits to Fermi pulsar light curves like Vela, with vacuum dipole fits being more favorable.

The high-quality Fermi LAT observations of gamma-ray pulsars have opened a new window to understanding the generation mechanisms of high-energy emission from these systems. The high statistics allow for careful modeling of the light curve features as well as for phase resolved spectral modeling. We modeled the LAT light curves of the Vela and CTA 1 pulsars with simulated high-energy light curves generated from geometrical representations of the outer gap and slot gap emission models, within the vacuum retarded dipole and force-free fields. A Markov Chain Monte Carlo maximum likelihood method was used to explore the phase space of the magnetic inclination angle, viewing angle, maximum emission radius, and gap width. We also used the measured spectral cutoff energies to estimate the accelerating parallel electric field dependence on radius, under the assumptions that the high-energy emission is dominated by curvature radiation and the geometry (radius of emission and minimum radius of curvature of the magnetic field lines) is determined by the best fitting light curves for each model. We find that light curves from the vacuum field more closely match the observed light curves and multiwavelength constraints, and that the calculated parallel electric field can place additional constraints on the emission geometry.

We present the structure of the 3D ideal MHD pulsar magnetosphere to a radius ten times that of the light cylinder, a distance about an order of magnitude larger than any previous such numerical treatment. Its overall structure exhibits a stable, smooth, well-defined undulating current sheet which approaches the kinematic split monopole solution of Bogovalov 1999 only after a careful introduction of diffusivity even in the highest resolution simulations. It also exhibits an intriguing spiral region at the crossing of two zero charge surfaces on the current sheet, which shows a destabilizing behavior more prominent in higher resolution simulations. We discuss the possibility that this region is physically (and not numerically) unstable. Finally, we present the spiral pulsar antenna radiation pattern.

We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday's equation, give rise to a respective normal field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modelled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition - the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.

We study the formation of the spiral structure of barred spiral galaxies, using an $N$-body model. The evolution of this $N$-body model in the adiabatic approximation maintains a strong spiral pattern for more than 10 bar rotations. We find that this longevity of the spiral arms is mainly due to the phenomenon of stickiness of chaotic orbits close to the unstable asymptotic manifolds originated from the main unstable periodic orbits, both inside and outside corotation. The stickiness along the manifolds corresponding to different energy levels supports parts of the spiral structure. The loci of the disc velocity minima (where the particles spend most of their time, in the configuration space) reveal the density maxima and therefore the main morphological structures of the system. We study the relation of these loci with those of the apocentres and pericentres at different energy levels. The diffusion of the sticky chaotic orbits outwards is slow and depends on the initial conditions and the corresponding Jacobi constant.

We study the stellar response in a spectrum of potentials describing the barred spiral galaxy NGC 1300. These potentials have been presented in a previous paper and correspond to three different assumptions as regards the geometry of the galaxy. For each potential we consider a wide range of $\Omega_p$ pattern speed values. Our goal is to discover the geometries and the $\Omega_p$ supporting specific morphological features of NGC 1300. For this purpose we use the method of response models. In order to compare the images of NGC 1300 with the density maps of our models, we define a new index which is a generalization of the Hausdorff distance. This index helps us to find out quantitatively which cases reproduce specific features of NGC 1300 in an objective way. Furthermore, we construct alternative models following a Schwarzschild type technique. By this method we vary the weights of the various energy levels, and thus the orbital contribution of each energy, in order to minimize the differences between the response density and that deduced from the surface density of the galaxy, under certain assumptions. We find that the models corresponding to $\Omega_p\approx16$\ksk and $\Omega_p\approx22$\ksk are able to reproduce efficiently certain morphological features of NGC 1300, with each one having its advantages and drawbacks.

In a series of papers we study the stellar dynamics of the grand design barred-spiral galaxy NGC~1300. In the first paper of this series we estimate the gravitational potential and we give it in a form suitable to be used in dynamical studies. The estimation is done directly from near-infrared observations. Since the 3D distribution of the luminous matter is unknown, we construct three different general models for the potential corresponding to three different assumptions for the geometry of the system, representing limiting cases. A pure 2D disc, a cylindrical geometry (thick disc) and a third case, where a spherical geometry is assumed to apply for the major part of the bar. For the potential of the disc component on the galactic plane a Fourier decomposition method is used, that allows us to express it as a sum of trigonometric terms. Both even and odd components are considered, so that the estimated potential accounts also for the observed asymmetries in the morphology. For the amplitudes of the trigonometric terms a smoothed cubic interpolation scheme is used. The total potential in each model may include two additional terms (Plummer spheres) representing a central mass concentration and a dark halo component, respectively. In all examined models, the relative force perturbation points to a strongly nonlinear gravitational field, which ranges from 0.45 to 0.8 of the axisymmetric background with the pure 2D being the most nonlinear one. We present the topological distributions of the stable and unstable Lagrangian points as a function of the pattern speed $(\Omega_p)$. In all three models there is a range of $\Omega_p$ values, where we find multiple stationary points whose stability affects the overall dynamics of the system.

We present the orbital analysis of four response models, that succeed in reproducing morphological features of NGC 1300. Two of them assume a planar (2D) geometry with $\Omega_p$=22 and 16 \ksk respectively. The two others assume a cylindrical (thick) disc and rotate with the same pattern speeds as the 2D models. These response models reproduce most successfully main morphological features of NGC 1300 among a large number of models, as became evident in a previous study. Our main result is the discovery of three new dynamical mechanisms that can support structures in a barred-spiral grand design system. These mechanisms are presented in characteristic cases, where these dynamical phenomena take place. They refer firstly to the support of a strong bar, of ansae type, almost solely by chaotic orbits, then to the support of spirals by chaotic orbits that for a certain number of pat tern revolutions follow an n:1 (n=7,8) morphology, and finally to the support of spiral arms by a combination of orbits trapped around L$_{4,5}$ and sticky chaotic orbits with the same Jacobi constant. We have encountered these dynamical phenomena in a large fraction of the cases we studied as we varied the parameters of our general models, without forcing in some way their appearance. This suggests that they could be responsible for the observed morphologies of many barred-spiral galaxies. Comparing our response models among themselves we find that the NGC 130 0 morphology is best described by a thick disc model for the bar region and a 2D disc model for the spirals, with both components rotating with the same pattern speed $\Omega_p$=16 \ksk !. In such a case, the whole structure is included inside the corotation of the system. The bar is supported mainly by regular orbits, while the spirals are supported by chaotic orbits.

We study the orbital structure in a series of self-consistent $N$-body configurations simulating rotating barred galaxies with spiral and ring structures. We perform frequency analysis in order to measure the angular and the radial frequencies of the orbits at two different time snapshots during the evolution of each $N$-body system. The analysis is done separately for the regular and the chaotic orbits. We thereby identify the various types of orbits, determine the shape and percentages of the orbits supporting the bar and the ring/spiral structures, and study how the latter quantities change during the secular evolution of each system. Although the frequency maps of the chaotic orbits are scattered, we can still identify concentrations around resonances. We give the distributions of frequencies of the most important populations of orbits. We explore the phase space structure of each system using projections of the 4D surfaces of section. These are obtained via the numerical integration of the orbits of test particles, but also of the real $N$-body particles. We thus identify which domains of the phase space are preferred and which are avoided by the real particles. The chaotic orbits are found to play a major role in supporting the shape of the outer envelope of the bar as well as the rings and the spiral arms formed outside corotation.

The pattern speeds of spiral galaxies are closely related to the flow of material in their disks. Flows that follow the `precessing ellipses' paradigm (see e.g., Kalnajs 1973) are likely associated with slowly rotating spirals, which have corotation beyond their end. Such a flow can be secured by material trapped around stable, elliptical, x_1 periodic orbits precessing as their Jacobi constant varies. Contrarily, if part of the spiral arms is located at a corotation region then the spiral structure has to `survive' in chaotic regions. Barred-spiral systems with a single pattern speed and a bar ending before, but close to, corotation are candidates for having spirals supported by stars in chaotic motion. In this work we review the flows we have found in response models for various types of spiral potentials and indicate the cases, where order or chaos shapes the observed morphologies.

We investigate the strong electric current sheet that develops at the tip of the pulsar closed line region through time dependent three-dimensional numerical simulations of a rotating magnetic dipole. We show that curvature radiation from relativistic electrons and positrons in the current sheet may naturally account for several features of the high-energy pulsar emission. We obtain light curves and polarization profiles for the complete range of magnetic field inclination angles and observer orientations, and compare our results to recent observations from the Fermi gamma-ray telescope.

The unstable invariant manifolds of the short-period family of periodic orbits around the unstable Lagrangian points $L_1$ and $L_2$ of a barred galaxy define loci in the configuration space which take the form of a trailing spiral pattern. In the present paper we investigate this association in the case of the self-consistent models of Kaufmann & Contopoulos (1996) which provide an approximation of real barred-spiral galaxies. We also examine the relation of `response' models of barred-spiral galaxies with the theory of the invariant manifolds. Our main results are the following: The invariant manifolds yield the correct form of the imposed spiral pattern provided that their calculation is done with the spiral potential term turned on. We provide a theoretical model explaining the form of the invariant manifolds that supports the spiral structure. The azimuthal displacement of the Lagrangian points with respect to the bar's major axis is a crucial parameter in this modeling. When this is taken into account, the manifolds necessarily develop in a spiral-like domain of the configuration space, delimited from below by the boundary of a banana-like non-permitted domain, and from above either by rotational KAM tori or by cantori forming a stickiness zone. We construct `spiral response' models on the basis of the theory of the invariant manifolds and examine the connection of the latter to the `response' models (Patsis 2006) used to fit real barred-spiral galaxies, explaining how are the manifolds related to a number of morphological features seen in such models.

We investigate the three-dimensional structure of the pulsar magnetosphere through time-dependent numerical simulations of a magnetic dipole that is set in rotation. We developed our own Eulerian finite difference time domain numerical solver of force-free electrodynamics and implemented the technique of non-reflecting and absorbing outer boundaries. This allows us to run our simulations for many stellar rotations, and thus claim with confidence that we have reached a steady state. A quasi-stationary corotating pattern is established, in agreement with previous numerical solutions. We discuss the prospects of our code for future high-resolution investigations of dissipation, particle acceleration, and temporal variability.

We study a series of $N-$body simulations representing elliptical galaxies with central masses. Starting from two different systems with smooth centres, which have initially a triaxial configuration and are in equilibrium, we insert to them central masses of various values. Immediately after such an insertion a system presents a high fraction of particles moving in chaotic orbits, a fact causing a secular evolution towards a new equilibrium state. The chaotic orbits responsible for the secular evolution are identified. Their typical Lypaunov exponents are found to scale with the central mass as a power law $L\propto m^s$ with $s$ close to 1/2. The requirements for an effective secular evolution within a Hubble time are examined. These requirements are quantified by introducing a quantity called \emph{effective chaotic momentum} $\mathscr{L}$. This quantity is found to correlate well with the rate of the systems' secular evolution. In particular, we find that when $\mathscr{L}$ falls below a threshold value (0.004 in our $N-$body units) a system does no longer exhibit significant secular evolution.

We address the question of an appropriate choice of basis functions for the self-consistent field (SCF) method of simulation of the N-body problem. Our criterion is based on a comparison of the orbits found in N-body realizations of analytical potential-density models of triaxial galaxies, in which the potential is fitted by the SCF method using a variety of basis sets, with those of the original models. Our tests refer to maximally triaxial Dehnen gamma-models for values of $\gamma$ in the range 0<=gamma<=1. When an N-body realization of a model is fitted by the SCF method, the choice of radial basis functions affects significantly the way the potential, forces, or derivatives of the forces are reproduced, especially in the central regions of the system. We find that this results in serious discrepancies in the relative amounts of chaotic versus regular orbits, or in the distributions of the Lyapunov characteristic exponents, as found by different basis sets. Numerical tests include the Clutton-Brock and the Hernquist-Ostriker (HO) basis sets, as well as a family of numerical basis sets which are `close' to the HO basis set. The family of numerical basis sets is parametrized in terms of a quantity $\epsilon$ which appears in the kernel functions of the Sturm-Liouville (SL) equation defining each basis set. The HO basis set is the $\epsilon=0$ member of the family. We demonstrate that grid solutions of the SL equation yielding numerical basis sets introduce large errors in the variational equations of motion. We propose a quantum-mechanical method of solution of the SL equation which overcomes these errors. We finally give criteria for a choice of optimal value of $\epsilon$ and calculate the latter as a function of the value of gamma.

The present lecture notes are an introduction to selected topics of {\it Galactic Dynamics}. The focus is on topics that we consider more relevant to the main theme of this workshop, {\it Celestial Mechanics}. This is not intended to be a review article. In fact, any of the topics below could be the subject of a separate review. Only the main ideas and notions are introduced, as well as some important currently open problems in each topic. Some relevant results from our own research are also presented. We discuss topics related mostly to the so-called {\it ellipsoidal components} of galaxies. These are \textbf{a)} the dark halos of both elliptical and disk galaxies, \textbf{b)} the luminous matter in elliptical galaxies, and \textbf{c)} the bulges of disk galaxies. We shall only occasionally refer to the dynamics of disks, bars or spiral structure.

Dissipationless N-body models of rotating galaxies, iso-energetic to a non-rotating model, are examined as regards the mass in regular and in chaotic motion. The values of their spin parameters $\lambda$ are near the value $\lambda=0.22$ of our Galaxy. We obtain the distinction between the sets of particles moving in regular and in chaotic orbits and we show that the spatial distribution of these two sets of particles is much different. The rotating models are characterized by larger fractions of mass in chaotic motion ($\thickapprox 65%$) compared with the fraction of mass in chaotic motion in the non-rotating iso-energetic model ($\thickapprox 32%$). Furthermore, the Lyapunov numbers of the chaotic orbits in the rotating models become by about one order of magnitude larger than in the non-rotating model. Chaotic orbits are concentrated preferably in values of the Jacobi integral around the value of the effective potential at the corotation radius. We find that density waves form a central rotating bar embedded in a thin and a thick disc with exponential surface density profile. A surprising new result is that long living spiral arms are exited on the disc, composed almost completely by chaotic orbits. The bar excites an $m=2$ mode of spiral waves on the surface density of the disc, emanating from the corotation radius. These spiral waves are deformed, fade, or disappear temporarily, but they grow again re-forming a well developed spiral pattern. Spiral arms are discernible up to 20 or 30 rotations of the bar (lasting for about a Hubble time).