We investigate the three-dimensional motion of a test particle in the gravitational field generated by a non-spherical compact object endowed with a mass quadrupole moment, described by the Erez-Rosen metric, and a radiation field, including the general relativistic Poynting-Robertson effect, coming from a rigidly rotating spherical emitting source located outside of the compact object. We derive the equations of motion for test particles influenced by such radiation field, recovering the two-dimensional description, and the weak-field approximation. This dynamical system admits the existence of a critical hypersurface, region where gravitational and radiation forces balance. Selected test particle orbits for different set of input parameters are displayed. The possible configurations on the critical hypersurfaces can be either latitudinal drift towards the equatorial ring or suspended orbits. We discuss about the existence of multiple hypersurface solutions through a simple method to perform the calculations. We graphically prove also that the critical hypersurfaces are stable configurations within the Lyapunov theory.

Objectives: A systematic study on the general relativistic Poynting-Robertson effect has been developed so far by introducing different complementary approaches, which can be mainly divided in two kinds: (1) improving the theoretical assessments and model in its simple aspects, and (2) extracting mathematical and physical information from such system with the aim to extend methods or results to other similar physical systems of analogue structure. Methods/Analysis: We use these theoretical approaches: relativity of observer splitting formalism; Lagrangian formalism and Rayleigh potential with a new integration method; Lyapunov theory os stability. Findings: We determined the three-dimensional formulation of the general relativistic Poynting-Robertson effect model. We determine the analytical form of the Rayleigh potential and discuss its implications. We prove that the critical hypersurfaces (regions where there is a balance between gravitational and radiation forces) are stable configurations. Novelty /Improvement: Our new contributions are: to have introduced the three-dimensional description; to have determined the general relativistic Rayleigh potential for the first time in the General Relativity literature; to have provided an alternative, general and more elegant proof of the stability of the critical hypersurfaces.

We derive the equations of motion of a test particle in the equatorial plane around a static and spherically symmetric wormhole influenced by a radiation field including the general relativistic Poynting-Robertson effect. From the analysis of this dynamical system, we develop a diagnostic to distinguish a black hole from a wormhole, which can be timely supported by several and different observational data. This procedure is based on the possibility of having some wormhole metrics, which smoothly connect to the Schwarzschild metric in a small transition surface layer very close to the black hole event horizon. To detect such a metric-change, we analyse the emission proprieties from the critical hypersurface (stable region where radiation and gravitational fields balance) together with those from an accretion disk in the Schwarzschild spacetime toward a distant observer. Indeed, if the observational data are well fitted within such model, it immediately implies the existence of a black hole; while in case of strong departures from such description it means that a wormhole could be present. Finally, we discuss our results and draw the conclusions.

We present a new method to measure the radius-to-mass ratio (R/M) of weakly magnetic, disc-accreting neutron stars by exploiting the occultation of parts of the inner disc by the star itself. This occultation imprints characteristic features on the X-ray line profile that are unique and are expected to be present in low mass X-ray binary systems seen under inclinations higher than ~65 degrees. We analyse a NuSTAR observation of a good candidate system, 4U 1636-53, and find that X-ray spectra from current instrumentation are unlikely to single out the occultation features owing to insufficient signal-to-noise. Based on an extensive set of simulations we show that large-area X-ray detectors of the future generation could measure R/M to ~2{\div}3% precision over a range of inclinations. Such is the precision in radius determination required to derive tight constraints on the equation of state of ultradense matter and it represents the goal that other methods too aim to achieve in the future.

We determine for the first time in the literature the analytic form of the Rayleigh potential of the general relativistic Poynting-Robertson effect. The employed procedure is based on the use of an integrating factor and a new integration strategy where the test particle's dissipated energy represents the fundamental variable. The obtained results and their implications are discussed. Finally, concluding remarks and future projects are drawn.

In this paper we introduce a new procedure, termed by us \emph{energy formalism}, to deal with dissipative systems in metric theories of gravity. This approach aims at determining the analytic expression of Rayleigh dissipation function in the context of the inverse problem in the calculus of variations. We describe our method in detail, presenting a simple example. After, we consider as first extensive application the general relativistic Poynting-Robertson effect. The obtained results and future implications are discussed.

We prove the stability of the critical hypersurfaces associated with the three-dimensional general relativistic Poynting-Robertson effect. The equatorial ring configures to be as a stable attractor and the whole critical hypersurface as a basin of attraction for this dynamical system. We introduce a new, simpler (in terms of calculations), and more physical approach within the Lyapunov theory. We propose three different Lyapunov functions, each one carrying important information and very useful for understanding such phenomenon under different aspects.

In my doctoral thesis, I have focussed my attention on radiation processes in high-energy astrophysics connected with the accretion flow physics around compact objects. Generally, a radiation field beside to exert an outward radiation pressure, there is also the presence of a radiation drag force, which both can drastically change or even halt the motion of the surrounding matter. The radiation drag force, known as Poynting-Robertson effect, acts as a dissipative force against the matter's orbital motion, removing very efficiently angular momentum and energy from it. The thesis is organised in three parts: (1) for ray-tracing purposes, I have developed a mathematical method for deriving a set of high-accurate approximate polynomial formulae to easily integrate photon geodesics in a Schwarzschild spacetime; (2) I gave two fundamental contributions in the field of the general relativistic treatment of the Poynting-Robertson effect (Lagrangian formulations and extension of the model in three dimensions); (3) I reduced the data of three accreting millisecond X-ray pulsars: IGR J00291+5934, IGR J18245-2452, and SAX J1748.9-2021. This thesis offers innovative ideas in the field of radiation processes involving the Poynting-Robertson effect in high-energy astrophysics, opening thus up future interesting perspectives both in theoretical and observational physics.

SAX J1748.9-2021 is a transiently accreting X-ray millisecond pulsar. It is also known as an X-ray burster source discovered by Beppo-SAX. We analysed the persistent emission and type-I X-ray burst properties during its 2015 outburst. The source varied from hard to soft state within half day. We modeled the broad-band spectra of the persistent emission in the 1 - 250 keV energy band for both spectral states using the quasi-simultaneous INTEGRAL and Swift data. The broad-band spectra are well fitted by an absorbed thermal Componization model, compps, in a slab geometry. The best-fits for the two states indicate significantly different plasma temperature of 18 and 5 keV and the Thomson optical depth of 3 and 4, respectively. In total, 56 type-I X-ray bursts were observed during the 2015 outburst, of which 26 detected by INTEGRAL in the hard state, 25 by XMM-Newton in the soft state, and 5 by Swift in both states. As the object transited from the hard to the soft state, the recurrence time for X-ray bursts decreased from $\approx 2$ to $\approx1$ hr. The relation between the recurrence time, $ \Delta t_{\rm rec}$, and the local mass accretion rate per unit area onto the compact object, $\dot m$, is fitted by a power-law model, and yielded as best fit at $\Delta t_{\rm rec} \sim \langle \dot{m} \rangle^{-1.02\pm0.03}$ using all X-ray bursts. In both cases, the observed recurrence times are consistent with the mixed hydrogen/helium bursts. We also discuss the effects of type-I X-ray bursts prior to the hard to soft transition.

We propose a polynomial approximation of the global Lense-Thirring rigid precession frequency to study low frequency quasi-periodic oscillations around spinning black holes. This high-performing approximation allows to determine the expected frequencies of a precessing thick accretion disc with fixed inner radius and variable outer radius around a black hole with given mass and spin. We discuss the accuracy and the applicability regions of our polynomial approximation, showing that the computational times are reduced by a factor of $\approx70$ in the range of minutes.

IGR~J18245--2452/PSR J1824--2452I is one of the rare transitional accreting millisecond X-ray pulsars, showing direct evidence of switches between states of rotation powered radio pulsations and accretion powered X-ray pulsations, dubbed transitional pulsars. IGR~J18245--2452 is the only transitional pulsar so far to have shown a full accretion episode, reaching an X-ray luminosity of $\sim10^{37}$~erg~s$^{-1}$ permitting its discovery with INTEGRAL in 2013. In this paper, we report on a detailed analysis of the data collected with the IBIS/ISGRI and the two JEM-X monitors on-board INTEGRAL at the time of the 2013 outburst. We make use of some complementary data obtained with the instruments on-board XMM-Newton and Swift in order to perform the averaged broad-band spectral analysis of the source in the energy range 0.4 -- 250~keV. We have found that this spectrum is the hardest among the accreting millisecond X-ray pulsars. We improved the ephemeris, now valid across its full outburst, and report the detection of pulsed emission up to $\sim60$ keV in both the ISGRI ($10.9 \sigma$) and Fermi/GBM ($5.9 \sigma$) bandpass. The alignment of the ISGRI and Fermi GBM 20 -- 60 keV pulse profiles are consistent at a $\sim25\ \mu$s level. We compared the pulse profiles obtained at soft X-rays with \xmm\ with the soft \gr-ray ones, and derived the pulsed fractions of the fundamental and first harmonic, as well as the time lag of the fundamental harmonic, up to $150\ \mu$s, as a function of energy. We report on a thermonuclear X-ray burst detected with \Integ, and using the properties of the previously type-I X-ray burst, we show that all these events are powered primarily by helium ignited at a depth of $y_{\rm ign} \approx 2.7\times10^8$ g cm${}^{-2}$. For such a helium burst the estimated recurrence time of $\Delta t_{\rm rec}\approx5.6$ d is in agreement with the observations.

IGR J00291+5934 is the fastest-known accretion-powered X-ray pulsar, discovered during a transient outburst in 2004. In this paper, we report on Integral and Swift observations during the 2015 outburst, which lasts for $\sim25$ d. The source has not been observed in outburst since 2008, suggesting that the long-term accretion rate has decreased by a factor of two since discovery. The averaged broad-band (0.1 - 250 keV) persistent spectrum in 2015 is well described by a thermal Comptonization model with a column density of $N_{\rm H} \approx4\times10^{21}$ cm$^{-2}$, a plasma temperature of $kT_{\rm e} \approx50$ keV, and a Thomson optical depth of $\tau_{\rm T}\approx1$. Pulsations at the known spin period of the source are detected in the Integral data up to the $\sim150$ keV energy band. We also report on the discovery of the first thermonuclear burst observed from IGR J00291+5934, which lasts around 7 min and occurs at a persistent emission level corresponding to roughly $1.6\%$ of the Eddington accretion rate. The properties of the burst suggest it is powered primarily by helium ignited at a depth of $y_{\rm ign}\approx1.5\times10^9$ g cm$^{-2}$ following the exhaustion by steady burning of the accreted hydrogen. The Swift/BAT data from the first $\sim20$ s of the burst provide indications of a photospheric radius expansion phase. Assuming this is the case, we infer a source distance of $d = 4.2 \pm 0.3$ kpc.

We develop a method for deriving approximate analytical formulae to integrate photon geodesics in a Schwarzschild spacetime. Based on this, we derive the approximate equations for light bending and propagation delay that have been introduced empirically. We then derive for the first time an approximate analytical equation for the solid angle. We discuss the accuracy and range of applicability of the new equations and present a few simple applications of them to known astrophysical problems.