Source author record

Constantinos Kalapotharakos

Constantinos Kalapotharakos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

The Fundamental Plane Relation for Gamma-Ray Pulsars Implied by 4FGL

We explore the validity of the recently reported fundamental plane (FP) relation of gamma-ray pulsars using 190 pulsars included in the latest 4FGL-DR3 catalog. This sample number is more than twice as large as that of the original study. The FP relation incorporates 4 parameters, i.e., the spin-down power, $\dot{\mathcal{E}}$, the surface magnetic field, $B_{\star}$, the total gamma-ray luminosity, $L_γ$, and a spectral cutoff energy, $ε_{\rm cut}$. The derivation of $ε_{\rm cut}$ is the most intriguing one because $ε_{\rm cut}$ depends on the proper interpretation of the available phase-averaged spectra. We construct synthetic phase-averaged spectra, guided by the few existing phase-resolved ones, to find that the best fit cutoff energy, $ε_{\rm c1}$, corresponding to a purely exponential cutoff (plus a power law) spectral form, is the parameter that optimally probes the maximum cutoff energy of the emission that originates from the core of the dissipative region, i.e., the equatorial current sheet. Computing this parameter for the 190 4FGL pulsars, we find that the resulting FP relation, i.e. the gamma-ray luminosity in terms of the other observables, reads $L_γ=10^{14.3\pm 1.3}(ε_{\rm c1}/{\rm MeV})^{1.39\pm0.17}(B_{\star}/{\rm G})^{0.12\pm 0.03}(\dot{\mathcal{E}}/{\rm erg\;s^{-1}})^{0.39\pm 0.05}{\rm ~erg\;s^{-1}}$; this is in good agreement with both the empirical relation reported by Kalapotharakos et al. (2019) and the theoretically predicted relation for curvature radiation. Finally, we revisit the radiation reaction limited condition, to find it is a sufficient but not necessary condition for the theoretical derivation of the FP relation. However, the assumption of the radiation reaction limited acceleration reveals the underlying accelerating electric field component and its scaling with $\dot{\mathcal{E}}$.

preprint2021arXiv

Probing the High-energy $γ$-ray Emission Mechanism in the Vela Pulsar via Phase-resolved Spectral and Energy-dependent Light Curve Modeling

Recent kinetic simulations sparked a debate regarding the emission mechanism responsible for pulsed GeV $γ$-ray emission from pulsars. Some models invoke curvature radiation, while other models assume synchrotron radiation in the current-sheet. We interpret the curved spectrum of the Vela pulsar as seen by H.E.S.S. II (up to $\sim$100 GeV) and the $Fermi$ Large Area Telescope (LAT) to be the result of curvature radiation due to primary particles in the pulsar magnetosphere and current sheet. We present phase-resolved spectra and energy-dependent light curves using an extended slot gap and current sheet model, invoking a step function for the accelerating electric field as motivated by kinetic simulations. We include a refined calculation of the curvature radius of particle trajectories in the lab frame, impacting the particle transport, predicted light curves, and spectra. Our model reproduces the decrease of the flux of the first peak relative to the second one, evolution of the bridge emission, near-constant phase positions of peaks, and narrowing of pulses with increasing energy. We can explain the first of these trends because we find that the curvature radii of the particle trajectories in regions where the second $γ$-ray light curve peak originates are systematically larger than those associated with the first peak, implying that the spectral cutoff of the second peak is correspondingly larger. However, an unknown azimuthal dependence of the $E$-field as well as uncertainty in the precise spatial origin of the GeV emission, precludes a simplistic discrimination of emission mechanisms.

preprint2021arXiv

Very-High-Energy Emission From Pulsars

Air-Cherenkov telescopes have detected pulsations at energies above 50 GeV from a growing number of Fermi pulsars. These include the Crab, Vela, PSR B1706-44 and Geminga, with the first two having pulsed detections above 1 TeV. In some cases, there appears to be very-high-energy (VHE) emission that is an extension of the Fermi spectra to high energies, while in other cases, additional higher-energy spectral components that require a separate emission mechanism may be present. We present results of broad-band spectral modeling using global magnetosphere fields and multiple emission mechanisms that include synchro-curvature (SC) and inverse Compton scattered (ICS) radiation from accelerated particles (primaries) and synchrotron-self Compton (SSC) emission from lower-energy pairs. Our models predict three distinct VHE components: SC from primaries whose high-energy tail can extend to 100 GeV, SSC from pairs that can extend to several TeV and ICS from primary particles accelerated in the current sheet, scattering pair synchrotron radiation, that appears beyond 10 TeV. Our models suggest that H.E.S.S.-II and MAGIC have detected the high-energy tail of the primary SC component that produces the Fermi spectrum in Vela, Geminga and PSR B1706-44. We argue that the ICS component peaking above 10 TeV from Vela has been seen by H.E.S.S. Detection of this emission component from the Crab and other pulsars is possible with HAWC and CTA, and directly measures the maximum particle energy in pulsars.

preprint2020arXiv

The Multipolar Magnetic Field of Millisecond Pulsar PSR J0030+0451

Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed to constrain the neutron star mass and radius, has inferred surface hot-spots (the magnetic polar caps) that imply significantly non-dipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole plus quadrupole components using static vacuum field and force-free global magnetosphere models. Taking into account the compactness and observer angle values provided by Miller et al. (2019) and Riley et al. (2019), we compute geodesics from the observer plane to the polar caps to compute the resulting X-ray light curve. We explore, through Markov chain Monte Carlo techniques, the detailed magnetic field configurations that can reproduce the observed X-ray light curve and have discovered degeneracies, i.e., diverse field configurations, which can provide sufficient descriptions to the NICER X-ray waveforms. Having obtained the force-free field structures, we then compute the corresponding synchronous gamma-ray light curves following Kalapotharakos et al. (2014) these we compare to those obtained by Fermi-LAT, to provide models consistent with both the X-ray and the gamma-ray data, thereby restricting further the multipole field parameters. An essential aspect of this approach is the proper computation of the relative phase between the synchronous X- and gamma-ray light curves. We conclude with a discussion of the broader implications of our study.

preprint2015arXiv

Synchrotron Self-Compton Emission from the Crab and Other Pulsars

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.

preprint2015arXiv

Testing dissipative magnetosphere model light curves and spectra with FERMI pulsars

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 (σ). 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 γ-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 σ 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.

preprint2014arXiv

Gamma-Ray Emission in Dissipative Pulsar Magnetospheres: From Theory to Fermi Observations

We compute the patterns of $γ$-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 $γ$-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 $γ$-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 $γ$-ray emission on the $γ$-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 $γ$-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.

preprint2012arXiv

Gamma-Ray Light Curves from Pulsar Magnetospheres with Finite Conductivity

We investigate the shapes of γ-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 σ. (b) We define realistic trajectories of radiating particles in magnetospheres of different σ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 σ, especially in the second approach. While still not self-consistent, these results are a step forward in understanding the physics of pulsar γ-radiation.

preprint2012arXiv

Toward a Realistic Pulsar Magnetosphere

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=ρ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.

preprint2011arXiv

Gamma-Ray Pulsar Light Curves in Vacuum and Force-Free Geometry

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.

preprint2011arXiv

Pulsar Emission Geometry and Accelerating Field Strength

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.

preprint2011arXiv

The Extended Pulsar Magnetosphere

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.

preprint2010arXiv

Asymptotic Orbits in Barred Spiral Galaxies

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.

preprint2010arXiv

Nonlinear force-free reconstruction of the global solar magnetic field: methodology

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.

preprint2010arXiv

Orbital structure in barred spiral galaxies

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.

preprint2009arXiv

The pulsar synchrotron in 3D: curvature radiation

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.

preprint2006arXiv

Special Features of Galactic Dynamics

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.