Fast radio bursts (FRBs) are short (millisecond) radio pulses originating from enigmatic sources at extragalactic distances so far lacking a detection in other energy bands. Magnetized neutron stars (magnetars) have been considered as the sources powering the FRBs, but the connection is controversial because of differing energetics and the lack of radio and X-ray detections with similar characteristics in the two classes. We report here the detection by the AGILE satellite on April 28, 2020 of an X-ray burst in coincidence with the very bright radio burst from the Galactic magnetar SGR 1935+2154. The burst detected by AGILE in the hard X-ray band (18-60 keV) lasts about 0.5 seconds, it is spectrally cutoff above 80 keV, and implies an isotropically emitted energy ~ $10^{40}$ erg. This event is remarkable in many ways: it shows for the first time that a magnetar can produce X-ray bursts in coincidence with FRB-like radio bursts; it also suggests that FRBs associated with magnetars may emit X-ray bursts of both magnetospheric and radio-pulse types that may be discovered in nearby sources. Guided by this detection, we discuss SGR 1935+2154 in the context of FRBs, and especially focus on the class of repeating-FRBs. Based on energetics, magnetars with fields B ~ $10^{15}$ G may power the majority of repeating-FRBs. Nearby repeating-FRBs offer a unique occasion to consolidate the FRB-magnetar connection, and we present new data on the X-ray monitoring of nearby FRBs. Our detection enlightens and constrains the physical process leading to FRBs: contrary to previous expectations, high-brightness temperature radio emission coexists with spectrally-cutoff X-ray radiation.

Measuring the cosmic ray flux over timescales comparable to the age of the solar system, $\sim 4.5\,$Gyr, could provide a new window on the history of the Earth, the solar system, and even our galaxy. We present a technique to indirectly measure the rate of cosmic rays as a function of time using the imprints of atmospheric neutrinos in paleo-detectors, natural minerals which record damage tracks from nuclear recoils. Minerals commonly found on Earth are $\lesssim 1\,$Gyr old, providing the ability to look back across cosmic ray history on timescales of the same order as the age of the solar system. Given a collection of differently aged samples dated with reasonable accuracy, this technique is particularly well-suited to measuring historical changes in the cosmic ray flux at Earth and is broadly applicable in astrophysics and geophysics.

FRB 180916 is a most intriguing source at 150 Mpc distance capable of producing repeating fast radio bursts with a periodic 16.35 day temporal pattern. We report on the X-ray and $\gamma$-ray observations of FRB 180916 obtained by AGILE and Swift. We focused on the recurrent 5-day time intervals of active radio bursting and present results obtained on Feb. 3 - 8; Feb. 25; Mar. 5 - 10; Mar. 22 - 28, 2020 during a multiwavelength campaign involving high-energy and radio observations. We also searched for temporal coincidences at millisecond timescales between all known radio bursts of FRB 180916 and X-ray and MeV events detectable by AGILE. We do not detect any simultaneous event or any extended X-ray and $\gamma$-ray emission on timescales of hours/days/weeks. Our cumulative X-ray (0.3-10 keV) flux upper limit of $5 \times\,10^{-14} \rm \, erg \, cm^{-2} s^{-1}$ (obtained during 5-day active intervals) translates into an isotropic luminosity upper limit of $L_{X,UL} \sim 1.5 \times\, 10^{41} \rm erg \, s^{-1}$. Observations above 100 MeV over a many-year timescale provide an average luminosity upper limit one order of magnitude larger. These results provide the so-far most stringent limits on high-energy emission from FRB 180916 and constrain the dissipation of magnetic energy from a magnetar-like source of radius $R_m$, internal magnetic field $B_m$ and dissipation timescale $\tau_d$ to satisfy the relation $R_{m,6}^3 B_{m,16}^2 \tau_{d,8}^{-1} \lesssim 1$, where $R_{m,6}$ is $R_m$ in units of $10^6$ cm, $B_{m,16}$ is $B_m$ in units of $10^{16}$ G, and $\tau_{d,8}$ in units of $10^8$ s.

The interactions of cosmic rays with the solar atmosphere produce secondary particle which can reach the Earth. In this work we present a comprehensive calculation of the yields of secondary particles as gamma-rays, electrons, positrons, neutrons and neutrinos performed with the FLUKA code. We also estimate the intensity at the Sun and the fluxes at the Earth of these secondary particles by folding their yields with the intensities of cosmic rays impinging on the solar surface. The results are sensitive on the assumptions on the magnetic field nearby the Sun and to the cosmic-ray transport in the magnetic field in the inner solar system.

We focus on two repeating fast radio bursts (FRBs) recently detected by the CHIME/FRB experiment in 2018--2019 (Source 1: 180916.J0158+65, and Source 2: 181030.J1054+73). These sources have low excess dispersion measures (DMs) ($ < 100 \rm \, pc \, cm^{-3}$ and $ < 20 \rm \, pc \, cm^{-3}$, respectively), implying relatively small maximal distances. They were repeatedly observed by AGILE in the MeV--GeV energy range. We do not detect prompt emission simultaneously with these repeating events. This search is particularly significant for the submillisecond and millisecond integrations obtainable by AGILE. The sources are constrained to emit a MeV-fluence in the millisecond range below $F'_{MeV} = 10^{-8} \, \rm erg \, cm^{-2}$ corresponding to an isotropic energy near $E_{MeV,UL} \simeq 2 \times 10^{46}\,$erg for a distance of 150 Mpc (applicable to Source 1). We also searched for $\gamma$-ray emission for time intervals up to 100 days, obtaining 3$\,\sigma$ upper limits (ULs) for the average isotropic luminosity above 50 MeV, $L_{\gamma,UL} \simeq \,$(5-10)$\,\times 10^{43} \rm \, erg \, s^{-1}$. For a source distance near 100 kpc (possibly applicable to Source 2), our ULs imply $E_{MeV,UL}\simeq10^{40} \rm erg$, and $L_{\gamma,UL} \simeq \,$2$\,\times 10^{37} \rm \, erg \, s^{-1}$. Our results are significant in constraining the high-energy emission of underlying sources such as magnetars, or other phenomena related to extragalactic compact objects, and show the prompt emission to be lower than the peak of the 2004 magnetar outburst of SGR 1806-20 for source distances less than about 100 Mpc.

In the era of big data, radio astronomical image reconstruction algorithms are challenged to estimate clean images given limited computing resources and time. This article is driven by the need for large scale image reconstruction for the future Square Kilometre Array (SKA), which will become in the next decades the largest low and intermediate frequency radio telescope in the world. This work proposes a scalable wideband deconvolution algorithm called MUFFIN, which stands for "MUlti Frequency image reconstruction For radio INterferometry". MUFFIN estimates the sky images in various frequency bands given the corresponding dirty images and point spread functions. The reconstruction is achieved by minimizing a data fidelity term and joint spatial and spectral sparse analysis regularization terms. It is consequently non-parametric w.r.t. the spectral behaviour of radio sources. MUFFIN algorithm is endowed with a parallel implementation and an automatic tuning of the regularization parameters, making it scalable and well suited for big data applications such as SKA. Comparisons between MUFFIN and the state-of-the-art wideband reconstruction algorithm are provided.

The origin of cosmic neutrinos is still largely unknown. Using data obtained by the gamma-ray imager on board of the AGILE satellite, we systematically searched for transient gamma-ray sources above 100 MeV that are temporally and spatially coincident with ten recent high-energy neutrino IceCube events. We find three AGILE candidate sources that can be considered possible counterparts to neutrino events. Detecting 3 gamma-ray/neutrino associations out of 10 IceCube events is shown to be unlikely due to a chance coincidence. One of the sources is related to the BL Lac source TXS 0506+056. For the other two AGILE gamma-ray sources there are no obvious known counterparts, and both Galactic and extragalactic origin should be considered.

In radio astronomy, accurate calibration is of crucial importance for the new generation of radio interferometers. More specifically, because of the potential presence of outliers which affect the measured data, robustness needs to be ensured. On the other hand, calibration is improved by taking advantage of these new instruments and exploiting the known structure of parameters of interest across frequency. Therefore, we propose in this paper an iterative robust multi-frequency calibration algorithm based on a distributed and consensus optimization scheme which aims to estimate the complex gains of the receivers and the directional perturbations caused by the ionosphere. Numerical simulations reveal that the proposed distributed calibration technique outperforms the conventional non-robust algorithm and per-channel calibration.

This paper investigates calibration of sensor arrays in the radio astronomy context. Current and future radio telescopes require computationally efficient algorithms to overcome the new technical challenges as large collecting area, wide field of view and huge data volume. Specifically, we study the calibration of radio interferometry stations with significant direction dependent distortions. We propose an iterative robust calibration algorithm based on a relaxed maximum likelihood estimator for a specific context: i) observations are affected by the presence of outliers and ii) parameters of interest have a specific structure depending on frequency. Variation of parameters across frequency is addressed through a distributed procedure, which is consistent with the new radio synthesis arrays where the full observing bandwidth is divided into multiple frequency channels. Numerical simulations reveal that the proposed robust distributed calibration estimator outperforms the conventional non-robust algorithm and/or the mono-frequency case.

In this paper, we aim to design robust estimation techniques based on the compound-Gaussian (CG) process and adapted for calibration of radio interferometers. The motivation beyond this is due to the presence of outliers leading to an unrealistic traditional Gaussian noise assumption. Consequently, to achieve robustness, we adopt a maximum a posteriori (MAP) approach which exploits Bayesian statistics and follows a sequential updating procedure here. The proposed algorithm is applied in a multi-frequency scenario in order to enhance the estimation and correction of perturbation effects. Numerical simulations assess the performance of the proposed algorithm for different noise models, Student's t, K, Laplace, Cauchy and inverse-Gaussian compound-Gaussian distributions w.r.t. the classical non-robust Gaussian noise assumption.

AGILE is a mission of the Italian Space Agency (ASI) Scientific Program dedicated to gamma-ray astrophysics, operating in a low Earth orbit since April 23, 2007. It is designed to be a very light and compact instrument, capable of simultaneously detecting and imaging photons in the 18 keV to 60 keV X-ray energy band and in the 30 MeV{50 GeV gamma-ray energy with a good angular resolution (< 1 deg at 1 GeV). The core of the instrument is the Silicon Tracker complemented with a CsI calorimeter and a AntiCoincidence system forming the Gamma Ray Imaging Detector (GRID). Before launch, the GRID needed on-ground calibration with a tagged gamma-ray beam to estimate its performance and validate the Monte Carlo simulation. The GRID was calibrated using a tagged gamma-ray beam with energy up to 500 MeV at the Beam Test Facilities at the INFN Laboratori Nazionali di Frascati. These data are used to validate a GEANT3 based simulation by comparing the data and the Monte Carlo simulation by measuring the angular and energy resolutions. The GRID angular and energy resolutions obtained using the beam agree well with the Monte Carlo simulation. Therefore the simulation can be used to simulate the same performance on-light with high reliability.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

We report on the measurement of the $^{7}$Be($n, p$)$^{7}$Li cross section from thermal to approximately 325 keV neutron energy, performed in the high-flux experimental area (EAR2) of the n_TOF facility at CERN. This reaction plays a key role in the lithium yield of the Big Bang Nucleosynthesis (BBN) for standard cosmology. The only two previous time-of-flight measurements performed on this reaction did not cover the energy window of interest for BBN, and showed a large discrepancy between each other. The measurement was performed with a Si-telescope, and a high-purity sample produced by implantation of a $^{7}$Be ion beam at the ISOLDE facility at CERN. While a significantly higher cross section is found at low-energy, relative to current evaluations, in the region of BBN interest the present results are consistent with the values inferred from the time-reversal $^{7}$Li($p, n$)$^{7}$Be reaction, thus yielding only a relatively minor improvement on the so-called Cosmological Lithium Problem (CLiP). The relevance of these results on the near-threshold neutron production in the p+$^{7}$Li reaction is also discussed.

The LIGO-Virgo Collaboration (LVC) detected, on 2017 August 17, an exceptional gravitational-wave (GW) event temporally consistent within $\sim\,1.7 \, \rm s$ with the GRB 1708117A observed by Fermi-GBM and INTEGRAL. The event turns out to be compatible with a neutron star-neutron star (NS-NS) coalescence that subsequently produced a radio/optical/X-ray transient detected at later times. We report the main results of the observations by the AGILE satellite of the GW170817 localization region (LR) and its electromagnetic (e.m.) counterpart. At the LVC detection time $T_0$, the GW170817 LR was occulted by the Earth. The AGILE instrument collected useful data before and after the GW-GRB event because in its spinning observation mode it can scan a given source many times per hour. The earliest exposure of the GW170817 LR by the gamma-ray imaging detector (GRID) started about 935 s after $T_0$. No significant X-ray or gamma-ray emission was detected from the LR that was repeatedly exposed over timescales of minutes, hours, and days before and after GW170817, also considering Mini-calorimeter and Super-AGILE data. Our measurements are among the earliest ones obtained by space satellites on GW170817 and provide useful constraints on the precursor and delayed emission properties of the NS-NS coalescence event. We can exclude with high confidence the existence of an X-ray/gamma-ray emitting magnetar-like object with a large magnetic field of $10^{15} \, \rm G$. Our data are particularly significant during the early stage of evolution of the e.m. remnant.

We present a new formulation of a family of proximity operators that generalize the projector step for phase retrieval. These proximity operators for noisy intensity measurements can replace the classical "noise free" projection in any projection-based algorithm. They are derived from a maximum likelihood formulation and admit closed form solutions for both the Gaussian and the Poisson cases. In addition, we extend these proximity operators to undersampled intensity measurements. To assess their performance, these operators are exploited in a classical Gerchberg Saxton algorithm. We present numerical experiments showing that the reconstructed complex amplitudes with these proximity operators perform always better than using the classical intensity projector while their computational overhead is moderate.

The LIGO/Virgo Collaboration (LVC) detected on 2017 January 4, a significant gravitational-wave (GW) event (now named GW170104). We report in this Letter the main results obtained from the analysis of hard X-ray and gamma-ray data of the AGILE mission that repeatedly observed the GW170104 localization region (LR). At the LVC detection time $T_0$ AGILE observed about 36% of the LR. The gamma-ray imaging detector did not reveal any significant emission in the energy range 50 MeV--30 GeV. Furthermore, no significant gamma-ray transients were detected in the LR that was repeatedly exposed over timescales of minutes, hours, and days. We also searched for transient emission using data near $T_0$ of the omnidirectional detector MCAL operating in the energy band 0.4--100 MeV. A refined analysis of MCAL data shows the existence of a weak event (that we call "E2") with a signal-to-noise ratio of $4.4\,\sigma$ lasting about 32 ms and occurring $0.46\,\pm\,0.05 \,\rm s$ before $T_0$. A study of the MCAL background and of the false-alarm rate of E2 leads to the determinination of a post-trial significance between $2.4\,\sigma$ and $2.7\,\sigma$ for a temporal coincidence with GW170104. We note that E2 has characteristics similar to those detected from the weak precursor of the short GRB 090510. The candidate event E2 is worth consideration for simultaneous detection by other satellites. If associated with GW170104, it shows emission in the MeV band of a short burst preceding the final coalescence by 0.46 sec and involving $\sim 10^{-7}$ of the total rest mass energy of the system.

The measured fluxes of secondary particles produced by the interactions of Cosmic Rays (CRs) with the astronomical environment play a crucial role in understanding the physics of CR transport. In this work we present a comprehensive calculation of the secondary hadron, lepton, gamma-ray and neutrino yields produced by the inelastic interactions between several species of stable or long-lived cosmic rays projectiles (p, D, T, 3He, 4He, 6Li, 7Li, 9Be, 10Be, 10B, 11B, 12C, 13C, 14C, 14N, 15N, 16O, 17O, 18O, 20Ne, 24Mg and 28Si) and different target gas nuclei (p, 4He, 12C, 14N, 16O, 20Ne, 24Mg, 28Si and 40Ar). The yields are calculated using FLUKA, a simulation package designed to compute the energy distributions of secondary products with large accuracy in a wide energy range. The present results provide, for the first time, a complete and self-consistent set of all the relevant inclusive cross sections regarding the whole spectrum of secondary products in nuclear collisions. We cover, for the projectiles, a kinetic energy range extending from $0.1 GeV/n$ up to $100 TeV/n$ in the lab frame. In order to show the importance of our results for multi-messenger studies about the physics of CR propagation, we evaluate the propagated spectra of Galactic secondary nuclei, leptons, and gamma rays produced by the interactions of CRs with the insterstellar gas, exploiting the numerical codes DRAGON and GammaSky. We show that, adopting our cross section database, we are able to provide a good fit of a complete sample of CR observables, including: leptonic and hadronic spectra measured at Earth, the local interstellar spectra measured by Voyager, and the gamma-ray emissivities from Fermi-LAT collaboration. We also show a set of gamma-ray and neutrino full-sky maps and spectra.

The measured fluxes of secondary particles produced by the interactions of cosmic rays with the astronomical environment represent a powerful tool to infer some properties of primary cosmic rays. In this work we investigate the production of secondary particles in inelastic hadronic interactions between several cosmic rays species of projectiles and different target nuclei of the interstellar medium. The yields of secondary particles have been calculated with the FLUKA simulation package, that provides with very good accuracy the energy distributions of secondary products in a large energy range. An application to the propagation and production of secondaries in the Galaxy is presented.

Magnetic reconnection is a plasma phenomenon where a topological rearrangement of magnetic field lines with opposite polarity results in dissipation of magnetic energy into heat, kinetic energy and particle acceleration. Such a phenomenon is considered as an efficient mechanism for energy release in laboratory and astrophysical plasmas. An important question is how to make the process fast enough to account for observed explosive energy releases. The classical model for steady state magnetic reconnection predicts reconnection times scaling as $S^{1/2}$ (where $S$ is the Lundquist number) and yields times scales several order of magnitude larger than the observed ones. Earlier two-dimensional MHD simulations showed that for large Lundquist number the reconnection time becomes independent of $S$ ("fast reconnection" regime) due to the presence of the secondary tearing instability that takes place for $S \gtrsim 1 \times 10^4$. We report on our 3D MHD simulations of magnetic reconnection in a magnetically confined cylindrical plasma column under either a pressure balanced or a force-free equilibrium and compare the results with 2D simulations of a circular current sheet. We find that the 3D instabilities acting on these configurations result in a fragmentation of the initial current sheet in small filaments, leading to enhanced dissipation rate that becomes independent of the Lundquist number already at $S \simeq 1\times 10^3$.

We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data.

We report the results of an extensive search in the AGILE data for a gamma-ray counterpart of the LIGO gravitational wave event GW150914. Currently in spinning mode, AGILE has the potential of covering with its gamma-ray instrument 80 % of the sky more than 100 times a day. It turns out that AGILE came within a minute from the event time of observing the accessible GW150914 localization region. Interestingly, the gamma-ray detector exposed about 65 % of this region during the 100 s time intervals centered at -100 s and +300 s from the event time. We determine a 2-sigma flux upper limit in the band 50 MeV - 10 GeV, $UL = 1.9 \times 10^{-8} \rm \, erg \, cm^{-2} \, s^{-1}$ obtained about 300 s after the event. The timing of this measurement is the fastest ever obtained for GW150914, and significantly constrains the electromagnetic emission of a possible high-energy counterpart. We also carried out a search for a gamma-ray precursor and delayed emission over timescales ranging from minutes to days: in particular, we obtained an optimal exposure during the interval -150 / -30 s. In all these observations, we do not detect a significant signal associated with GW150914. We do not reveal the weak transient source reported by Fermi-GBM 0.4 s after the event time. However, even though a gamma-ray counterpart of the GW150914 event was not detected, the prospects for future AGILE observations of gravitational wave sources are decidedly promising.

The advent of enhanced technologies in radio interferometry and the perspective of the SKA telescope bring new challenges in image reconstruction. One of these challenges is the spatio-spectral reconstruction of large (Terabytes) data cubes with high fidelity. This contribution proposes an alternative implementation of one such 3D prototype algorithm, MUFFIN (MUlti-Frequency image reconstruction For radio INterferometry), which combines spatial and spectral analysis priors. Using a recently proposed primal dual algorithm, this new version of MUFFIN allows a parallel implementation where computationally intensive steps are split by spectral channels. This parallelization allows to implement computationally demanding translation invariant wavelet transforms (IUWT), as opposed to the union of bases used previously. This alternative implementation is important as it opens the possibility of comparing these efficient dictionaries, and others, in spatio-spectral reconstruction. Numerical results show that the IUWT-based version can be successfully implemented at large scale with performances comparable to union of bases.

We describe a "spatio-spectral" deconvolution algorithm for wide-band imaging in radio interferometry. In contrast with the existing multi-frequency reconstruction algorithms, the proposed method does not rely on a model of the sky-brightness spectral distribution. This non-parametric approach can be of particular interest for the new generation of low frequency radiotelescopes. The proposed solution formalizes the reconstruction problem as a convex optimization problem with spatial and spectral regularizations. The efficiency of this approach has been already proven for narrow-band image reconstruction and the present contribution can be considered as its extension to the multi-frequency case. Because the number of frequency bands multiplies the size of the inverse problem, particular attention is devoted to the derivation of an iterative large scale optimization algorithm. It is shown that the main computational bottleneck of the approach, which lies in the resolution of a linear system, can be efficiently overcome by a fully parallel implementation w.r.t. the frequencies, where each processor reconstructs a narrow-band image. All the other optimization steps are extremely fast. A parallel implementation of the algorithm in Julia is publicly available at https://github.com/andferrari. Preliminary simulations illustrate the performances of the method and its ability to reconstruct complex spatio-spectral structures.

Current and future radio interferometric arrays such as LOFAR and SKA are characterized by a paradox. Their large number of receptors (up to millions) allow theoretically unprecedented high imaging resolution. In the same time, the ultra massive amounts of samples makes the data transfer and computational loads (correlation and calibration) order of magnitudes too high to allow any currently existing image reconstruction algorithm to achieve, or even approach, the theoretical resolution. We investigate here decentralized and distributed image reconstruction strategies which select, transfer and process only a fraction of the total data. The loss in MSE incurred by the proposed approach is evaluated theoretically and numerically on simple test cases.

Astronomical optical interferometers (OI) sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid atmospheric perturbations, the phases of the complex Fourier samples (visibilities) cannot be directly exploited , and instead linear relationships between the phases are used (phase closures and differential phases). Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic OI instruments are now paving the way to multiwavelength imaging. This paper presents the derivation of a spatio-spectral ("3D") image reconstruction algorithm called PAINTER (Polychromatic opticAl INTErferometric Reconstruction software). The algorithm is able to solve large scale problems. It relies on an iterative process, which alternates estimation of polychromatic images and of complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also from differential phases, which help to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm.

Astronomical optical interferometers sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid perturbations caused by atmospheric turbulence, the phases of the complex Fourier samples (visibilities) cannot be directly exploited. Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic optical interferometric instruments are now paving the way to multiwavelength imaging. This paper is devoted to the derivation of a spatio-spectral (3D) image reconstruction algorithm, coined PAINTER (Polychromatic opticAl INTErferometric Reconstruction software). The algorithm relies on an iterative process, which alternates estimation of polychromatic images and of complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also differential phases, which helps to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm and in particular the relevance of injecting a differential phases model in the reconstruction.

High-energy phenomena in the cosmos, and in particular processes leading to the emission of gamma- rays in the energy range 10 MeV - 100 GeV, play a very special role in the understanding of our Universe. This energy range is indeed associated with non-thermal phenomena and challenging particle acceleration processes. The technology involved in detecting gamma-rays is challenging and drives our ability to develop improved instruments for a large variety of applications. GAMMA-LIGHT is a Small Mission which aims at an unprecedented advance of our knowledge in many sectors of astrophysical and Earth studies research. The Mission will open a new observational window in the low-energy gamma-ray range 10-50 MeV, and is configured to make substantial advances compared with the previous and current gamma-ray experiments (AGILE and Fermi). The improvement is based on an exquisite angular resolution achieved by GAMMA-LIGHT using state-of-the-art Silicon technology with innovative data acquisition. GAMMA-LIGHT will address all astrophysics issues left open by the current generation of instruments. In particular, the breakthrough angular resolution in the energy range 100 MeV - 1 GeV is crucial to resolve patchy and complex features of diffuse sources in the Galaxy as well as increasing the point source sensitivity. This proposal addresses scientific topics of great interest to the community, with particular emphasis on multifrequency correlation studies involving radio, optical, IR, X-ray, soft gamma-ray and TeV emission. At the end of this decade several new observatories will be operational including LOFAR, SKA, ALMA, HAWK, CTA. GAMMA-LIGHT will "fill the vacuum" in the 10 MeV-10 GeV band, and will provide invaluable data for the understanding of cosmic and terrestrial high-energy sources.

In optical interferometry, the visibility squared modulus are generally assumed to follow a Gaussian distribution and to be independent of each other. A quantitative analysis of the relevance of such assumptions is important to help improving the exploitation of existing and upcoming multi-wavelength interferometric instruments. Analyze the statistical behaviour of both the absolute and the colour-differential squared visibilities: distribution laws, correlations and cross-correlations between different baselines. We use observations of stellar calibrators obtained with AMBER instrument on VLTI in different instrumental and observing configurations, from which we extract the frame-by-frame transfer function. Statistical hypotheses tests and diagnostics are then systematically applied. For both absolute and differential squared visibilities and under all instrumental and observing conditions, we find a better fit for the Student distribution than for the Gaussian, log-normal and Cauchy distributions. We find and analyze clear correlation effects caused by atmospheric perturbations. The differential squared visibilities allow to keep a larger fraction of data with respect to selected absolute squared visibilities and thus benefit from reduced temporal dispersion, while their distribution is more clearly characterized. The frame selection based on the criterion of a fixed SNR value might result in either a biased sample of frames or in a too severe selection.

We present a variability study of a sample of bright gamma-ray (30 MeV -- 50 GeV) sources. This sample is an extension of the first AGILE catalogue of gamma-ray sources (1AGL), obtained using the complete set of AGILE observations in pointing mode performed during a 2.3 year period from July 9, 2007 until October 30, 2009. The dataset of AGILE pointed observations covers a long time interval and its gamma-ray data archive is useful for monitoring studies of medium-to-high brightness gamma-ray sources. In the analysis reported here, we used data obtained with an improved event filter that covers a wider field of view, on a much larger (about 27.5 months) dataset, integrating data on observation block time scales, which mostly range from a few days to thirty days. The data processing resulted in a better characterized source list than 1AGL was, and includes 54 sources, 7 of which are new high galactic latitude (|BII| >= 5) sources, 8 are new sources on the galactic plane, and 20 sources from the previous catalogue with revised positions. Eight 1AGL sources (2 high-latitude and 6 on the galactic plane) were not detected in the final processing either because of low OB exposure and/or due to their position in complex galactic regions. We report the results in a catalogue of all the detections obtained in each single OB, including the variability results for each of these sources. In particular, we found that 12 sources out of 42 or 11 out of 53 are variable, depending on the variability index used, where 42 and 53 are the number of sources for which these indices could be calculated. Seven of the 11 variable sources are blazars, the others are Crab pulsar+nebula, LS I +61{\deg}303, Cyg X-3, and 1AGLR J2021+4030.

Context: AGILE is a gamma-ray astrophysics mission which has been in orbit since 23 April 2007 and continues to operate reliably. The gamma-ray detector, AGILE-GRID, has observed Galactic and extragalactic sources, many of which were collected in the first AGILE Catalog. Aims: We present the calibration of the AGILE-GRID using in-flight data and Monte Carlo simulations, producing Instrument Response Functions (IRFs) for the effective area A_eff), Energy Dispersion Probability (EDP), and Point Spread Function (PSF), each as a function of incident direction in instrument coordinates and energy. Methods: We performed Monte Carlo simulations at different gamma-ray energies and incident angles, including background rejection filters and Kalman filter-based gamma-ray reconstruction. Long integrations of in-flight observations of the Vela, Crab and Geminga sources in broad and narrow energy bands were used to validate and improve the accuracy of the instrument response functions. Results: The weighted average PSFs as a function of spectra correspond well to the data for all sources and energy bands. Conclusions: Changes in the interpolation of the PSF from Monte Carlo data and in the procedure for construction of the energy-weighted effective areas have improved the correspondence between predicted and observed fluxes and spectra of celestial calibration sources, reducing false positives and obviating the need for post-hoc energy-dependent scaling factors. The new IRFs have been publicly available from the Agile Science Data Centre since November 25, 2011, while the changes in the analysis software will be distributed in an upcoming release.

We investigate the dynamical propagation of the South-East jet from the Crab pulsar interacting with supernova ejecta by means of three-dimensional relativistic MHD numerical simulations with the PLUTO code. The initial jet structure is set up from the inner regions of the Crab Nebula. We study the evolution of hot, relativistic hollow outflows initially carrying a purely azimuthal magnetic field. Our jet models are characterized by different choices of the outflow magnetization ($\sigma$ parameter) and the bulk Lorentz factor ($\gamma_{j}$). We show that the jet is heavily affected by the growth of current-driven kink instabilities causing considerable deflection throughout its propagation length. This behavior is partially stabilized by the combined action of larger flow velocities and/or reduced magnetic field strengths. We find that our best jet models are characterized by relatively large values of $\sigma$ ($\gtrsim 1$) and small values of $\gamma_{j}\simeq 2$. Our results are in good agreement with the recent X-ray (\textit{Chandra}) data of the Crab Nebula South-East jet indicating that the jet changes direction of propagation on a time scale of the order of few years. The 3D models presented here may have important implications in the investigation of particle acceleration in relativistic outflows.

The $^{63}$Ni($n, \gamma$) cross section has been measured for the first time at the neutron time-of-flight facility n\_TOF at CERN from thermal neutron energies up to 200 keV. In total, capture kernels of 12 (new) resonances were determined. Maxwellian Averaged Cross Sections were calculated for thermal energies from kT = 5 keV to 100 keV with uncertainties around 20%. Stellar model calculations for a 25 M$_\odot$ star show that the new data have a significant effect on the $s$-process production of $^{63}$Cu, $^{64}$Ni, and $^{64}$Zn in massive stars, allowing stronger constraints on the Cu yields from explosive nucleosynthesis in the subsequent supernova.

Gamma-ray emission from the Crab Nebula has been recently shown to be unsteady. In this paper, we study the flux and spectral variability of the Crab above 100 MeV on different timescales ranging from days to weeks. In addition to the four main intense and day-long flares detected by AGILE and Fermi-LAT between Sept. 2007 and Sept. 2012, we find evidence for week-long and less intense episodes of enhanced gamma-ray emission that we call "waves". Statistically significant "waves" show timescales of 1-2 weeks, and can occur by themselves or in association with shorter flares. We present a refined flux and spectral analysis of the Sept. - Oct. 2007 gamma-ray enhancement episode detected by AGILE that shows both "wave" and flaring behavior. We extend our analysis to the publicly available Fermi-LAT dataset and show that several additional "wave" episodes can be identified. We discuss the spectral properties of the September 2007 "wave"/flare event and show that the physical properties of the "waves" are intermediate between steady and flaring states. Plasma instabilities inducing "waves" appear to involve spatial distances l \sim 10^{16} cm and enhanced magnetic fields B \sim (0.5 - 1) mG. Day-long flares are characterized by smaller distances and larger local magnetic fields. Typically, the deduced total energy associated with the "wave" phenomenon (E_w \sim 10^{42} erg, where E_w is the kinetic energy of the emitting particles) is comparable with that associated to the flares, and can reach a few percent of the total available pulsar spindown energy. Most likely, flares and waves are the product of the same class of plasma instabilities that we show acting on different timescales and radiation intensities.

In 2008 AGILE and Fermi detected gamma-ray flaring activity from the unidentified EGRET source 3EG J1236+0457, recently associated with a flat spectrum radio quasar GB6 J1239+0443 at z=1.762. The optical counterpart of the gamma-ray source underwent a flux enhancement of a factor 15-30 in 6 years, and of ~10 in six months. We interpret this flare-up in terms of a transition from an accretion-disk dominated emission to a synchrotron-jet dominated one. We analysed a Sloan Digital Sky Survey (SDSS) archival optical spectrum taken during a period of low radio and optical activity of the source. We estimated the mass of the central black hole using the width of the CIV emission line. In our work, we have also investigated SDSS archival optical photometric data and UV GALEX observations to estimate the thermal-disk emission contribution of GB6 J1239+0443. Our analysis of the gamma-ray data taken during the flaring episodes indicates a flat gamma-ray spectrum, with an extension of up to 15 GeV, with no statistically-relevant sign of absorption from the broad line region, suggesting that the blazar-zone is located beyond the broad line region. This result is confirmed by the modeling of the broad-band spectral energy distribution (well constrained by the available multiwavelength data) of the flaring activity periods and by the accretion disk luminosity and black hole mass estimated by us using archival data.

During its first 2 years of operation, the gamma-ray AGILE satellite accumulated an extensive dataset for the Galactic plane. The data have been monitored for transient sources and several gamma-ray sources were detected. Their variability and possible association were studied. In this talk we will focus on the results of extensive observations of the Carina Region during the time period 2007 July - 2009 January, for a total livetime of ~130 days. The region is extremely complex, hosting massive star formation, with the remarkable colliding wind binary Eta Carinae, massive star clusters and HII regions (e.g. NGC 3324, RCW49, Westerlund II) and a giant molecular cloud extending over 150 pc (between l=284.7 and l=289). The Carina Nebula itself is the largest and IR highest surface brightness nebula of the Southern emisphere. We monitored several gamma ray sources in the Carina Region. In particular we detect a gamma ray source (1AGL J1043-5931) consistent with the position of Eta Carinae and report a remarkable 2-days gamma-ray flaring episode from this source on 2008 Oct 11-13. If 1AGL J1043-5931 is associated with the Eta Car system, our data provides the long sought first detection above 100 MeV of a colliding wind binary.

At the core of the AGILE scientific instrument, designed to operate on a satellite, there is the Gamma Ray Imaging Detector (GRID) consisting of a Silicon Tracker (ST), a Cesium Iodide Mini-Calorimeter and an Anti-Coincidence system of plastic scintillator bars. The ST needs an on-ground calibration with a $\gamma$-ray beam to validate the simulation used to calculate the energy response function and the effective area versus the energy and the direction of the $\gamma$ rays. A tagged $\gamma$-ray beam line was designed at the Beam Test Facility (BTF) of the INFN Laboratori Nazionali of Frascati (LNF), based on an electron beam generating $\gamma$ rays through bremsstrahlung in a position-sensitive target. The $\gamma$-ray energy is deduced by difference with the post-bremsstrahlung electron energy \cite{prest}-\cite{hasan}. The electron energy is measured by a spectrometer consisting of a dipole magnet and an array of position sensitive silicon strip detectors, the Photon Tagging System (PTS). The use of the combined BTF-PTS system as tagged photon beam requires understanding the efficiency of $\gamma$-ray tagging, the probability of fake tagging, the energy resolution and the relation of the PTS hit position versus the $\gamma$-ray energy. This paper describes this study comparing data taken during the AGILE calibration occurred in 2005 with simulation.

We report the AGILE gamma-ray observations and the results of the multiwavelength campaigns on seven flaring blazars detected by the mission: During two multiwavelength campaigns, we observed gamma-ray activity from two Flat Spectrum Radio Quasars of the Virgo region, e.g. 3C 279 and 3C 273 (the latter being the first extragalactic source simultaneously observed with the gamma-ray telescope and the hard X ray imager of the mission). Due to the large FOV of the AGILE/GRID instrument, we achieved an almost continuous coverage of the FSRQ 3C 454.3. The source showed flux above 10E-6 photons/cm2/s (E > 100 MeV) and showed day by day variability during all the AGILE observing periods. In the EGRET era, the source was found in high gamma-ray activity only once. An other blazar, PKS 1510-089 was frequently found in high gamma-ray activity. S5 0716+71, an intermediate BL Lac object, exhibited a very high gamma-ray activity and fast gamma-ray variability during a period of intense optical activity. We observed high gamma-ray activity from W Comae, a BL Lac object, and Mrk 421, an high energy peaked BL Lac object. For this source, a multiwavelength campaign from optical to TeV has been performed.

The AGILE scientific instrument has been calibrated with a tagged $\gamma$-ray beam at the Beam Test Facility (BTF) of the INFN Laboratori Nazionali di Frascati (LNF). The goal of the calibration was the measure of the Point Spread Function (PSF) as a function of the photon energy and incident angle and the validation of the Monte Carlo (MC) simulation of the silicon tracker operation. The calibration setup is described and some preliminary results are presented.

We report on the extremely intense and fast gamma-ray are above 100 MeV detected by AGILE from the Crab Nebula in mid-April 2011. This event is the fourth of a sequence of reported major gamma-ray flares produced by the Crab Nebula in the period 2007/mid-2011. These events are attributed to strong radiative and plasma instabilities in the inner Crab Nebula, and their properties are crucial for theoretical studies of fast and efficient particle acceleration up to 10^15 eV. Here we study the very rapid flux and spectral evolution of the event that reached on April 16, 2011 the record-high peak flux of F = (26 +/- 5) x 10^-6 ph cm^-2 s^-1 with a risetime timescale that we determine to be in the range 6-10 hrs. The peak flaring gamma-ray spectrum reaches a distinct maximum near 500 MeV with no substantial emission above 1 GeV. The very rapid risetime and overall evolution of the Crab Nebula are strongly constrain the acceleration mechanisms and challenge MHD models. We briefly discuss the theoretical implications of our observations.

We consider the problem of convergence in homogeneous shearing box simula- tions of magneto-rotationally driven turbulence. When there is no mean magnetic flux, if the equations are non dimensionalized with respect to the diffusive scale, the only free parameter in the problem is the size of the computational domain. The problem of convergence then relates to the asymptotic form of the solutions as the computational box size becomes large. By using a numerical code with a high order of accuracy we show that the solutions become asymptotically inde- pendent of domain size. We also show that cases with weak magnetic flux join smoothly to the zero flux cases as the flux vanishes. These results are consistent with the operation of a subcritical small-scale dynamo driving the turbulence. We conclude that for this type of turbulence the angular momentum transport is a proportional to the diffusive flux and therefore has limited relevance in as- trophysical situations.

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all frequencies, and during the last four years it has exhibited more than one gamma-ray flare per year, becoming the most active gamma-ray blazar in the sky. We present for the first time the multi-wavelength AGILE, SWIFT, INTEGRAL, and GASP-WEBT data collected in order to explain the extraordinary gamma-ray flare of 3C 454.3 which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 reached a peak flux (E>100 MeV) of F_gamma(p) = (6.8+-1.0)E-5 ph/cm2/s on a time scale of about 12 hours, more than a factor of 6 higher than the flux of the brightest steady gamma-ray source, the Vela pulsar, and more than a factor of 3 brighter than its previous super-flare on 2009 December 2-3. The multi-wavelength data make a thorough study of the present event possible: the comparison with the previous outbursts indicates a close similarity to the one that occurred in 2009. By comparing the broadband emission before, during, and after the gamma-ray flare, we find that the radio, optical and X-ray emission varies within a factor 2-3, whereas the gamma-ray flux by a factor of 10. This remarkable behavior is modeled by an external Compton component driven by a substantial local enhancement of soft seed photons.

We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak.

Stars form in regions of the galaxy that are denser and cooler than the mean interstellar medium. These regions are called Giant Molecular Clouds. At the beginning of their life, up to $10^5-10^6$ years, stars accrete matter from their rich surrounding environment and are origin of a peculiar phenomenon that is the jet emission. Jets from Young Stellar Objects (YSOs) are intensively studied by the astrophysical community by observations at different wavelengths, analytical and numerical modeling and laboratory experiments. Indications about the jet propagation and its resulting morphologies are here obtained by means of a combined study of hypersonic jets carried out both in the laboratory and by numerical simulations.

The well known Crab Nebula is at the center of the SN1054 supernova remnant. It consists of a rotationally-powered pulsar interacting with a surrounding nebula through a relativistic particle wind. The emissions originating from the pulsar and nebula have been considered to be essentially stable. Here we report the detection of strong gamma-ray (100 MeV-10 GeV) flares observed by the AGILE satellite in September, 2010 and October, 2007. In both cases, the unpulsed flux increased by a factor of 3 compared to the non-flaring flux. The flare luminosity and short timescale favor an origin near the pulsar, and we discuss Chandra Observatory X-ray and HST optical follow-up observations of the nebula. Our observations challenge standard models of nebular emission and require power-law acceleration by shock-driven plasma wave turbulence within a ~1-day timescale.

Identification of gamma-ray-emitting Galactic sources is a long-standing problem in astrophysics. One such source, 1AGL J2022+4032, coincident with the interior of the radio shell of the supernova remnant Gamma Cygni (SNR G78.2+2.1) in the Cygnus Region, has recently been identified by Fermi as a gamma-ray pulsar, LAT PSR J2021+4026. We present long-term observations of 1AGL J2022+4032 with the AGILE gamma-ray telescope, measuring its flux and light curve. We compare the light curve of 1AGL J2022+4032 with that of 1AGL J2021+3652 (PSR J2021+3651), showing that the flux variability of 1AGL J2022+4032 appears to be greater than the level predicted from statistical and systematic effects and producing detailed simulations to estimate the probability of the apparent observed variability. We evaluate the possibility that the gamma-ray emission may be due to the superposition of two or more point sources, some of which may be variable, considering a number of possible counterparts. We consider the possibility of a nearby X-ray quiet microquasar contributing to the flux of 1AGL J2022+4032 to be more likely than the hypotheses of a background blazar or intrinsic gamma-ray variabilty of LAT PSR J2021+4026.

Terrestrial Gamma-Ray Flashes (TGFs) are very short bursts of high energy photons and electrons originating in Earth's atmosphere. We present here a localization study of TGFs carried out at gamma-ray energies above 20 MeV based on an innovative event selection method. We use the AGILE satellite Silicon Tracker data that for the first time have been correlated with TGFs detected by the AGILE Mini-Calorimeter. We detect 8 TGFs with gamma-ray photons of energies above 20 MeV localized by the AGILE gamma-ray imager with an accuracy of 5-10 degrees at 50 MeV. Remarkably, all TGF-associated gamma rays are compatible with a terrestrial production site closer to the sub-satellite point than 400 km. Considering that our gamma rays reach the AGILE satellite at 540 km altitude with limited scattering or attenuation, our measurements provide the first precise direct localization of TGFs from space.

We present the gamma-ray data of the extraordinary flaring activity above 100 MeV from the flat spectrum radio quasar 3C 454.3 detected by AGILE during the month of December 2009. 3C 454.3, that has been among the most active blazars of the FSRQ type since 2007, was detected in the gamma-ray range with a progressively rising flux since November 10, 2009. The gamma-ray flux reached a value comparable with that of the Vela pulsar on December 2, 2009. Remarkably, between December 2 and 3, 2009 the source more than doubled its gamma-ray emission and became the brightest gamma-ray source in the sky with a peak flux of F_{\gamma,p} = (2000 \pm 400) x 10^-8 ph cm^-2 s^-1 for a 1-day integration above 100 MeV. The gamma-ray intensity decreased in the following days with the source flux remaining at large values near F \simeq (1000 \pm 200) x 10^-8 ph cm^-2 s^-1 for more than a week. This exceptional gamma-ray flare dissipated among the largest ever detected intrinsic radiated power in gamma-rays above 100 MeV (L_{\gamma, source, peak} \simeq 3 x 10^46 erg s^-1, for a relativistic Doppler factor of {\delta} \simeq 30). The total isotropic irradiated energy of the month-long episode in the range 100 MeV - 3 GeV is E_{\gamma,iso} \simeq 10^56 erg. We report the intensity and spectral evolution of the gamma-ray emission across the flaring episode. We briefly discuss the important theoretical implications of our detection.

During the month of December, 2009 the blazar 3C 454.3 became the brightest gamma-ray source in the sky, reaching a peak flux F ~2000E-8 ph/cm2/s for E > 100 MeV. Starting in November, 2009 intensive multifrequency campaigns monitored the 3C 454 gamma-ray outburst. Here we report the results of a 2-month campaign involving AGILE, INTEGRAL, Swift/XRT, Swift/BAT, RossiXTE for the high-energy observations, and Swift/UVOT, KANATA, GRT, REM for the near-IR/optical/UV data. The GASP/WEBT provided radio and additional optical data. We detected a long-term active emission phase lasting ~1 month at all wavelengths: in the gamma-ray band, peak emission was reached on December 2-3, 2009. Remarkably, this gamma-ray super-flare was not accompanied by correspondingly intense emission in the optical/UV band that reached a level substantially lower than the previous observations in 2007-2008. The lack of strong simultaneous optical brightening during the super-flare and the determination of the broad-band spectral evolution severely constrain the theoretical modelling. We find that the pre- and post-flare broad-band behavior can be explained by a one-zone model involving SSC plus external Compton emission from an accretion disk and a broad-line region. However, the spectra of the Dec. 2-3, 2009 super-flare and of the secondary peak emission on Dec. 9, 2009 cannot be satisfactorily modelled by a simple one-zone model. An additional particle component is most likely active during these states.

Cygnus X-1 is the archetypal black hole (BH) binary system in our Galaxy. We report the main results of an extensive search for transient gamma-ray emission from Cygnus X-1 carried out in the energy range 100 MeV - 3 GeV by the AGILE satellite, during the period 2007 July - 2009 October. The total exposure time is about 300 days, during which the source was in the "hard" X-ray spectral state. We divided the observing intervals in 2 or 4 week periods, and searched for transient and persistent emission. We report an episode of significant transient gamma-ray emission detected on 2009, October 16 in a position compatible with Cygnus X-1 optical position. This episode, occurred during a hard spectral state of Cygnus X-1, shows that a 1-2 day time variable emission above 100 MeV can be produced during hard spectral states, having important theoretical implications for current Comptonization models for Cygnus X-1 and other microquasars. Except for this one short timescale episode, no significant gamma-ray emission was detected by AGILE. By integrating all available data we obtain a 2$\sigma$ upper limit for the total integrated flux of $F_{\gamma,U.L.} = 3 \times 10^{-8} \rm ph cm^{-2} s^{-1}$ in the energy range 100 MeV - 3 GeV. We then clearly establish the existence of a spectral cutoff in the energy range 1-100 MeV that applies to the typical hard state outside the flaring period and that confirms the historically known spectral cutoff above 1 MeV.

We report on 18 months of multiwavelength observations of the blazar 3C 454.3 (Crazy Diamond) carried out in July 2007-January 2009. We show the results of the AGILE campaigns which took place on May-June 2008, July-August 2008, and October 2008-January 2009. During the May 2008-January 2009 period, the source average flux was highly variable, from an average gamma-ray flux F(E>100MeV) > 200E-8 ph/cm2/s in May-June 2008, to F(E>100MeV)~80E-8 ph/cm2/s in October 2008-January 2009. The average gamma-ray spectrum between 100 MeV and 1 GeV can be fit by a simple power law (Gamma_GRID ~ 2.0 to 2.2). Only 3-sigma upper limits can be derived in the 20-60 keV energy band with Super-AGILE. During July-August 2007 and May-June 2008, RXTE measured a flux of F(3-20 keV)= 8.4E-11 erg/cm2/s, and F(3-20 keV)=4.5E-11 erg/cm2/s, respectively and a constant photon index Gamma_PCA=1.65. Swift/XRT observations were carried out during all AGILE campaigns, obtaining a F(2-10 keV)=(0.9-7.5)E-11 erg/cm2/s and a photon index Gamma_XRT=1.33-2.04. BAT measured an average flux of ~5 mCrab. GASP-WEBT monitored 3C 454.3 during the whole 2007-2008 period from the radio to the optical. A correlation analysis between the optical and the gamma-ray fluxes shows a time lag of tau=-0.4 days. An analysis of 15 GHz and 43 GHz VLBI core radio flux observations shows an increasing trend of the core radio flux, anti- correlated with the higher frequency data. The modeling SEDs, and the behavior of the long-term light curves in different energy bands, allow us to compare the jet properties during different emission states, and to study the geometrical properties of the jet on a time-span longer than one year.

The Supernova Remnant (SNR) IC 443 is an intermediate-age remnant well known for its radio, optical, X-ray and gamma-ray energy emissions. In this Letter we study the gamma-ray emission above 100 MeV from IC 443 as obtained by the AGILE satellite. A distinct pattern of diffuse emission in the energy range 100 MeV-3 GeV is detected across the SNR with its prominent maximum (source "A") localized in the Northeastern shell with a flux F = (47 \pm 10) 10^{-8} photons cm^{-2} s^{-1} above 100 MeV. This location is the site of the strongest shock interaction between the SNR blast wave and the dense circumstellar medium. Source "A" is not coincident with the TeV source located 0.4 degree away and associated with a dense molecular cloud complex in the SNR central region. From our observations, and from the lack of detectable diffuse TeV emission from its Northeastern rim, we demonstrate that electrons cannot be the main emitters of gamma-rays in the range 0.1-10 GeV at the site of the strongest SNR shock. The intensity, spectral characteristics, and location of the most prominent gamma-ray emission together with the absence of co-spatial detectable TeV emission are consistent only with a hadronic model of cosmic-ray acceleration in the SNR. A high-density molecular cloud (cloud "E") provides a remarkable "target" for nucleonic interactions of accelerated hadrons: our results show enhanced gamma-ray production near the molecular cloud/shocked shell interaction site. IC 443 provides the first unambiguous evidence of cosmic-ray acceleration by SNRs.

In the simple vs composite hypothesis test with a proper prior, the Bayes Factor (BF) is shown to be the posterior mean of the Likelihood Ratio (LR). Therefore, the posterior standard deviation of the LR or rather its posterior cumulative density function can be used to indicate the significativity of a detection by the BF and this detection procedure can be computed from a single Markov Chain. It is applied and compared for exoplanet detection. The previous statistics can be expressed from the Fractional BF (FBF) \cite{ohagan95} and the Probability distribution of the LR (PLR) \cite{aitkin97}. Two properties of the PLR related to the GLRT are noted and a procedure to optimize the PLR and the FBF two-parameters detectors according to their ROC curves is proposed. The performances of all tests are compared.

Pulsars are known to power winds of relativistic particles that can produce bright nebulae by interacting with the surrounding medium. These pulsar wind nebulae (PWNe) are observed in the radio, optical, x-rays and, in some cases, also at TeV energies, but the lack of information in the gamma-ray band prevents from drawing a comprehensive multiwavelength picture of their phenomenology and emission mechanisms. Using data from the AGILE satellite, we detected the Vela pulsar wind nebula in the energy range from 100 MeV to 3 GeV. This result constrains the particle population responsible for the GeV emission, probing multivavelength PWN models, and establishes a class of gamma-ray emitters that could account for a fraction of the unidentified Galactic gamma-ray sources.

BL Lac objects of the intermediate subclass (IBLs) are known to emit a substantial fraction of their power in the energy range 0.1--10 GeV. Detecting gamma-ray emission from such sources provides therefore a direct probe of the emission mechanisms and of the underlying powerhouse. The AGILE gamma-ray satellite detected the remarkable IBL S5 0716+714 (z \simeq 0.3) during a high state in the period from 2007 September - October, marked by two very intense flares reaching peak fluxes of 200\times10^{-8} ph / cm^2 s above 100 MeV, with simultaneous optical and X-ray observations. We present here a theoretical model for the two major flares and discuss the overall energetics of the source. We conclude that 0716+714 is among the brightest BL Lac's ever detected at gamma-ray energies. Because of its high power and lack of signs for ongoing accretion or surrounding gas, the source is an ideal candidate to test the maximal power extractable from a rotating supermassive black hole via the pure Blandford-Znajek (BZ) mechanism. We find that during the 2007 gamma-ray flares our source approached or just exceeded the upper limit set by BZ for a black hole of mass 10^9 M_sun

The study of relativistic particle acceleration is a major topic of high-energy astrophysics. It is well known that massive black holes in active galaxies can release a substantial fraction of their accretion power into energetic particles, producing gamma-rays and relativistic jets. Galactic microquasars (hosting a compact star of 1-10 solar masses which accretes matter from a binary companion) also produce relativistic jets. However, no direct evidence of particle acceleration above GeV energies has ever been obtained in microquasar ejections, leaving open the issue of the occurrence and timing of extreme matter energization during jet formation. Here we report the detection of transient gamma-ray emission above 100 MeV from the microquasar Cygnus X-3, an exceptional X-ray binary which sporadically produces powerful radio jets. Four gamma-ray flares (each lasting 1-2 days) were detected by the AGILE satellite simultaneously with special spectral states of Cygnus X-3 during the period mid-2007/mid-2009. Our observations show that very efficient particle acceleration and gamma-ray propagation out of the inner disk of a microquasar usually occur a few days before major relativistic jet ejections. Flaring particle energies can be thousands of times larger than previously detected maximum values (with Lorentz factors of 105 and 102 for electrons and protons, respectively). We show that the transitional nature of gamma-ray flares and particle acceleration above GeV energies in Cygnus X-3 is clearly linked to special radio/X-ray states preceding strong radio flares. Thus gamma-rays provide unique insight into the nature of physical processes in microquasars.

We report on the second AGILE multiwavelength campaign of the blazar 3C 454.3 during the first half of December 2007. This campaign involved AGILE, Spitzer, Swift,Suzaku,the WEBT consortium,the REM and MITSuME telescopes,offering a broad band coverage that allowed for a simultaneous sampling of the synchrotron and inverse Compton (IC) emissions.The 2-week AGILE monitoring was accompanied by radio to optical monitoring by WEBT and REM and by sparse observations in mid-Infrared and soft/hard X-ray energy bands performed by means of Target of Opportunity observations by Spitzer, Swift and Suzaku, respectively.The source was detected with an average flux of~250x10^{-8}ph cm^-2s^-1 above 100 MeV,typical of its flaring states.The simultaneous optical and gamma-ray monitoring allowed us to study the time-lag associated with the variability in the two energy bands, resulting in a possible ~1-day delay of the gamma-ray emission with respect to the optical one. From the simultaneous optical and gamma-ray fast flare detected on December 12, we can constrain the delay between the gamma-ray and optical emissions within 12 hours. Moreover, we obtain three Spectral Energy Distributions (SEDs) with simultaneous data for 2007 December 5, 13, 15, characterized by the widest multifrequency coverage. We found that a model with an external Compton on seed photons by a standard disk and reprocessed by the Broad Line Regions does not describe in a satisfactory way the SEDs of 2007 December 5, 13 and 15. An additional contribution, possibly from the hot corona with T=10^6 K surrounding the jet, is required to account simultaneously for the softness of the synchrotron and the hardness of the inverse Compton emissions during those epochs.

SuperAGILE is the hard X-ray monitor of the AGILE gamma ray mission, in orbit since 23$^{rd}$ April 2007. It is an imaging experiment based on a set of four independent silicon strip detectors, equipped with one-dimensional coded masks, operating in the nominal energy range 18-60 keV. The main goal of SuperAGILE is the observation of cosmic sources simultaneously with the main gamma-ray AGILE experiment, the Gamma Ray Imaging Detector (GRID). Given its $\sim$steradian-wide field of view and its $\sim$15 mCrab day-sensitivity, SuperAGILE is also well suited for the long-term monitoring of Galactic compact objects and the detection of bright transients. The SuperAGILE detector properties and design allow for a 6 arcmin angular resolution in each of the two independent orthogonal projections of the celestial coordinates. Photon by photon data are continuously available by the experiment telemetry, and are used to derive images and fluxes of individual sources, with integration times depending on the source intensity and position in the field of view. In this paper we report on the main scientific results achieved by SuperAGILE over its first two years in orbit, until April 2009.

Short gamma-ray bursts (GRBs), typically lasting less than 2 s, are a special class of GRBs of great interest. We report the detection by the AGILE satellite of the short GRB 090510 which shows two clearly distinct emission phases: a prompt phase lasting ~ 200 msec and a second phase lasting tens of seconds. The prompt phase is relatively intense in the 0.3-10 MeV range with a spectrum characterized by a large peak/cutoff energy near 3 MeV, in this phase, no significant high-energy gamma-ray emission is detected. At the end of the prompt phase, intense gamma-ray emission above 30 MeV is detected showing a power-law time decay of the flux of the type t^-1.3 and a broad-band spectrum remarkably different from that of the prompt phase. It extends from sub-MeV to hundreds of MeV energies with a photon index alpha ~ 1.5. GRB 090510 provides the first case of a short GRB with delayed gamma-ray emission. We present the timing and spectral data of GRB 090510 and briefly discuss its remarkable properties within the current models of gamma-ray emission of short GRBs.

We report the detection by the AGILE satellite of a rapid gamma-ray flare from the powerful gamma-ray quasar PKS 1510-089, during a pointing centered on the Galactic Center region from 1 March to 30 March 2008. This source has been continuosly monitored in the radio-to-optical bands by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). Moreover, the gamma-ray flaring episode triggered three ToO observations by the Swift satellite in three consecutive days, starting from 20 March 2008. In the period 1-16 March 2008, AGILE detected gamma-ray emission from PKS 1510-089 at a significance level of 6.2-sigma with an average flux over the entire period of (84 +/- 17) x 10^{-8} photons cm^{-2} s^{-1} for photon energies above 100 MeV. After a predefined satellite re-pointing, between 17 and 21 March 2008, AGILE detected the source at a significance level of 7.3-sigma, with an average flux (E > 100 MeV) of (134 +/- 29) x 10^{-8} photons cm^{-2} s^{-1} and a peak level of (281 +/- 68) x 10^{-8} photons cm^{-2} s^{-1} with daily integration. During the observing period January-April 2008, the source also showed an intense and variable optical activity, with several flaring episodes and a significant increase of the flux was observed at millimetric frequencies. Moreover, in the X-ray band the Swift/XRT observations seem to show an harder-when-brighter behaviour of the source spectrum. The spectral energy distribution of mid-March 2008 is modelled with a homogeneous one-zone synchrotron self Compton emission plus contributions from inverse Compton scattering of external photons from both the accretion disc and the broad line region. Indeed, some features in the optical-UV spectrum seem to indicate the presence of Seyfert-like components, such as the little blue bump and the big blue bump.

Relativistic magnetized jets are key elements in Active Galactic Nuclei and in other astrophysical environments. Their structure and evolution involves a complex nonlinear physics that can be approached by numerical studies only. Still, owing to a number of challenging computational aspects, only a few numerical investigations have been undertaken so far. In this paper, we present high-resolution three dimensional numerical simulations of relativistic magnetized jets carrying an initially toroidal magnetic field. The presence of a substantial toroidal component of the field is nowadays commonly invoked and held responsible for the process of jet acceleration and collimation. We find that the typical nose cone structures, commonly observed in axisymmetric two-dimensional simulations, are not produced in the 3D case. Rather, the toroidal field gives rise to strong current driven kink instabilities leading to jet wiggling. However, it appears to be able to maintain an highly relativistic spine along its full length. By comparing low and high resolution simulations, we emphasize the impact of resolution on the jet dynamical properties.

We present the first catalog of high-confidence gamma-ray sources detected by the AGILE satellite during observations performed from July 9, 2007 to June 30, 2008. Catalogued sources are detected by merging all the available data over the entire time period. AGILE, launched in April 2007, is an ASI mission devoted to gamma-ray observations in the 30 MeV - 50 GeV energy range, with simultaneous X-ray imaging capability in the 18-60 keV band. This catalog is based on Gamma-Ray Imaging Detector (GRID) data for energies greater than 100 MeV. For the first AGILE catalog we adopted a conservative analysis, with a high-quality event filter optimized to select gamma-ray events within the central zone of the instrument Field of View (radius of 40 degrees). This is a significance-limited (4 sigma) catalog, and it is not a complete flux-limited sample due to the non-uniform first year AGILE sky coverage. The catalog includes 47 sources, 21 of which are associated with confirmed or candidate pulsars, 13 with Blazars (7 FSRQ, 4 BL Lacs, 2 unknown type), 2 with HMXRBs, 2 with SNRs, 1 with a colliding-wind binary system, 8 with unidentified sources.

We present the results of extensive observations by the gamma-ray AGILE satellite of the Galactic region hosting the Carina nebula and the remarkable colliding wind binary Eta Carinae (Eta Car) during the period 2007 July to 2009 January. We detect a gamma-ray source (1AGL J1043-5931) consistent with the position of Eta Car. If 1AGL J1043-5931 is associated with the Eta Car system our data provide the long sought first detection above 100 MeV of a colliding wind binary. The average gamma-ray flux above 100 MeV and integrated over the pre-periastron period 2007 July to 2008 October is F = (37 +/- 5) x 10-8 ph cm-2 s-1 corresponding to an average gamma-ray luminosity of L = 3.4 x 10^34 erg s-1 for a distance of 2.3 kpc. We also report a 2-day gamma-ray flaring episode of 1AGL J1043-5931 on 2008 Oct. 11-13 possibly related to a transient acceleration and radiation episode of the strongly variable shock in the system.

This paper gives a simple and original presentation of various coronagraphs inherited from the Lyot coronagraph. We first present the Lyot and Roddier phase mask coronagraphs and study their properties as a function of the focal mask size. We show that the Roddier phase mask can be used to produce an apodization for the star. Optimal coronagraphy can be obtained from two main approaches, using prolate spheroidal pupil apodization and a finite-size focal mask, or using a clear aperture and an infinite mask of variable transmission.

The Mini-Calorimeter (MCAL) instrument on-board the AGILE satellite is a non-imaging gamma-ray scintillation detector sensitive in the 300keV-100MeV energy range with a total on-axis geometrical area of 1400cm^2. Gamma-Ray Bursts (GRBs) are one of the main scientific targets of the AGILE mission and the MCAL design as an independent self-triggering detector makes it a valuable all-sky monitor for GRBs. Furthermore MCAL is one of the very few operative instruments with microsecond timing capabilities in the MeV range. In this paper the results of GRB detections with MCAL after one year of operation in space are presented and discussed. A flexible trigger logic implemented in the AGILE payload data-handling unit allows the on-board detection of GRBs. For triggered events, energy and timing information are sent to telemetry on a photon-by-photon basis, so that energy and time binning are limited by counting statistics only. When the trigger logic is not active, GRBs can be detected offline in ratemeter data, although with worse energy and time resolution. Between the end of June 2007 and June 2008 MCAL detected 51 GRBs, with a detection rate of about 1 GRB/week, plus several other events at a few milliseconds timescales. Since February 2008 the on-board trigger logic has been fully active. Comparison of MCAL detected events and data provided by other space instruments confirms the sensitivity and effective area estimations. MCAL also joined the 3rd Inter-Planetary Network, to contribute to GRB localization by means of triangulation.

[Abridged] We report on a multiwavelength observation of the blazar 3C 454.3 (which we dubbed "crazy diamond") carried out on November 2007 by means of the astrophysical satellites AGILE, INTEGRAL, Swift, the WEBT Consortium, and the optical-NIR telescope REM. 3C 454.3 is detected at a $\sim 19-\sigma$ level during the 3-week observing period, with an average flux above 100 MeV of $F_{\rm E>100MeV} = (170 \pm 13) \times 10^{-8}$ \phcmsec. The gamma-ray spectrum can be fit with a single power-law with photon index $\Gamma_{\rm GRID} = 1.73 \pm 0.16$ between 100 MeV and 1 GeV. We detect significant day-by-day variability of the gamma-ray emission during our observations, and we can exclude that the fluxes are constant at the 99.6% ($\sim 2.9 \sigma$) level. The source was detected typically around 40 degrees off-axis, and it was substantially off--axis in the field of view of the AGILE hard X-ray imager. However, a 5-day long ToO observation by INTEGRAL detected 3C 454.3 at an average flux of about $F_{\rm 20-200 keV} = 1.49 \times 10^{-3}$ \phcmsec with an average photon index of $\Gamma_{\rm IBIS} = 1.75 \pm 0.24$ between 20--200 keV. Swift also detected 3C 454.3 with a flux in the 0.3--10 keV energy band in the range $(1.23-1.40) \times 10^{-2}$ \phcmsec{} and a photon index in the range $\Gamma_{\rm XRT} = 1.56-1.73$. In the optical band, both WEBT and REM show an extremely variable behavior in the $R$ band. A correlation analysis based on the entire data set is consistent with no time-lags between the gamma-ray and the optical flux variations. Our simultaneous multifrequency observations strongly indicate that the dominant emission mechanism between 30 MeV and 30 GeV is dominated by inverse Compton scattering of relativistic electrons in the jet on the external photons from the broad line region.

We report the gamma-ray activity from the intermediate BL Lac S5 0716+714 during observations acquired by the AGILE satellite in September and October 2007. These detections of activity were contemporaneous with a period of intense optical activity, which was monitored by GASP-WEBT. This simultaneous optical and gamma-ray coverage allows us to study in detail the light curves, time lags, gamma-ray photon spectrum, and Spectral Energy Distributions (SEDs) during different states of activity. AGILE observed the source with its two co-aligned imagers, the Gamma-Ray Imaging Detector (GRID) and the hard X-ray imager (Super-AGILE), which are sensitive to the 30 MeV-50 GeV and 18-60 keV energy ranges, respectively. Observations were completed in two different periods, the first between 2007 September 4-23, and the second between 2007 October 24-November 1. Over the period 2007 September 7-12, AGILE detected gamma-ray emission from the source at a significance level of 9.6-sigma with an average flux (E>100 MeV) of (97 +/- 15) x 10^{-8} photons cm^{-2} s^{-1}, which increased by a factor of at least four within three days. No emission was detected by Super-AGILE for the energy range 18-60 keV to a 3-sigma upper limit of 10 mCrab in 335 ksec. In October 2007, AGILE repointed toward S5 0716+714 following an intense optical flare, measuring an average flux of (47 +/- 11) x 10^{-8} photons cm^{-2} s^{-1} at a significance level of 6.0-sigma. The gamma-ray flux of S5 0716+714 detected by AGILE is the highest ever detected for this blazar and one of the most intense gamma-ray fluxes detected from a BL Lac object. The SED of mid-September appears to be consistent with the synchrotron self-Compton (SSC) emission model, but only by including two SSC components of different variabilities.

We investigate the linear stability properties of the plane interface separating two relativistic magnetized flows in relative motion. The two flows are governed by the (special) relativistic equations for a magnetized perfect gas in the infinite conductivity approximation. By adopting the vortex-sheet approximation, the relativistic magnetohydrodynamics equations are linearized around the equilibrium state and the corresponding dispersion relation is derived and discussed. The behavior of the configuration and the regimes of instability are investigated following the effects of four physical parameters, namely: the flow velocity, the relativistic and Alfv\'enic Mach numbers and the inclination of the wave vector on the plane of the interface. From the numerical solution of the dispersion relation, we find in general two separate regions of instability, associated respectively with the slow and fast magnetosonic modes. Modes parallel to the flow velocity are destabilized only for sufficiently low magnetization. For the latter case, stabilization is attained, additionally, at sufficiently large relativistic velocities between the two flows in relative motion. The relevance of these results to the study of the stability of astrophysical jets is briefly commented.

Strong observational evidence indicates that all extragalactic jets associated with AGNs move at relativistic speed up to 100 pc - 1 kpc scales from the nucleus. At larger distances, reflecting the Fanaroff-Riley radio source classification, we observe an abrupt deceleration in FR-I jets while relativistic motions persist up to Mpc scale in FR-II. Moreover, VLBI observations of some object like B2 1144+35, Mrk501 and M87 show limb brightening of the jet radio emission at the parsec scale. This effect is interpreted kinematically as due to the presence of a deboosted central spine at high Lorentz factor and of a weakly relativistic external layer. In this paper we investigate whether these effects can be interpreted by a breaking of the collimated flow by external medium entrainment favored by shear instabilities, namely Kelvin-Helmholtz instabilities. We examine in details the physical conditions under which significant deceleration of a relativistic flow is produced. We investigate the phenomenon by means of high-resolution three-dimensional relativistic hydrodynamic simulations using the PLUTO code for computational astrophysics. We find that the parameter of utmost importance in determining the instability evolution and the entrainment properties is the ambient/jet density contrast. We show that lighter jets suffer stronger slowing down in the external layer than in the central part and conserve a central spine at high Lorentz factor. Our model is verified by constructing synthetic emission maps from the numerical simulations that compare reasonably well with VLBI observations of the inner part of FR-I sources.

Aims: We study the changes in the properties of turbulence driven by the magnetorotational instability in a shearing box, as the computational domain size in the radial direction is varied relative to the height Methods: We perform 3D simulations in the shearing box approximation, with a net magnetic flux, and we consider computational domains with different aspect ratios Results: We find that in boxes of aspect ratio unity the transport of angular momentum is strongly intermittent and dominated by channel solutions in agreement with previous work. In contrast, in boxes with larger aspect ratio, the channel solutions and the associated intermittent behavior disappear. Conclusions: There is strong evidence that, as the aspect ratio becomes larger, the characteristics of the solution become aspect ratio independent. We conclude that shearing box calculations with aspect ratio unity or near unity may introduce spurious effects.

Primary GCR interact with the Earth's atmosphere originating atmospheric showers, thus giving rise to fluxes of secondary particles in the atmosphere. Electromagnetic and hadronic interactions interplay in the production of these particles, whose detection is performed by means of complementary techniques in different energy ranges and at different depths in the atmosphere, down to the Earth's surface. Monte Carlo codes are essential calculation tools which can describe the complexity of the physics of these phenomena, thus allowing the analysis of experimental data. However, these codes are affected by important uncertainties, concerning, in particular, hadronic physics at high energy. In this paper we shall report some results concerning inclusive particle fluxes and atmospheric shower properties as obtained using the FLUKA transport and interaction code. Some emphasis will also be given to the validation of the physics models of FLUKA involved in these calculations.

We investigate the stability, nonlinear development and equilibrium structure of vortices in a background shearing Keplerian flow. We make use of high-resolution global two-dimensional compressible hydrodynamic simulations. We introduce the concept of nonlinear adjustment to describe the transition of unbalanced vortical fields to a long-lived configuration. We discuss the conditions under which vortical perturbations evolve into long-lived persistent structures and we describe the properties of these equilibrium vortices. The properties of equilibrium vortices appear to be independent from the initial conditions and depend only on the local disk parameters. In particular we find that the ratio of the vortex size to the local disk scale height increases with the decrease of the sound speed, reaching values well above the unity. The process of spiral density wave generation by the vortex, discussed in our previous work, appear to maintain its efficiency also at nonlinear amplitudes and we observe the formation of spiral shocks attached to the vortex. The shocks may have important consequences on the long term vortex evolution and possibly on the global disk dynamics. Our study strengthens the arguments in favor of anticyclonic vortices as the candidates for the promotion of planetary formation. Hydrodynamic shocks that are an intrinsic property of persistent vortices in compressible Keplerian flows are an important contributor to the overall balance. These shocks support vortices against viscous dissipation by generating local potential vorticity and should be responsible for the eventual fate of the persistent anticyclonic vortices. Numerical codes have be able to resolve shock waves to describe the vortex dynamics correctly.

This paper is concerned with the theoretical properties of high contrast coronagraphic images in the context of exoplanet searches. We derive and analyze the statistical properties of the residual starlight in coronagraphic images, and describe the effect of a coronagraph on the speckle and photon noise. Current observations with coronagraphic instruments have shown that the main limitations to high contrast imaging are due to residual quasi-static speckles. We tackle this problem in this paper, and propose a generalization of our statistical model to include the description of static, quasi-static and fast residual atmospheric speckles. The results provide insight into the effects on the dynamic range of wavefront control, coronagraphy, active speckle reduction, and differential speckle calibration. The study is focused on ground-based imaging with extreme adaptive optics, but the approach is general enough to be applicable to space, with different parameters.

Statistical properties of the intensity in adaptive optics images are usually modeled with a Rician distribution. We study the central point of the image, where this model is inappropriate for high to very high correction levels. The central point is an important problem because it gives the Strehl ratio distribution. We show that the central point distribution can be modeled using a non-central Gamma distribution.

We present a new numerical code, PLUTO, for the solution of hypersonic flows in 1, 2 and 3 spatial dimensions and different systems of coordinates. The code provides a multi-physics, multi-algorithm modular environment particularly oriented towards the treatment of astrophysical flows in presence of discontinuities. Different hydrodynamic modules and algorithms may be independently selected to properly describe Newtonian, relativistic, MHD or relativistic MHD fluids. The modular structure exploits a general framework for integrating a system of conservation laws, built on modern Godunov-type shock-capturing schemes. Although a plethora of numerical methods has been successfully developed over the past two decades, the vast majority shares a common discretization recipe, involving three general steps: a piecewise polynomial reconstruction followed by the solution of Riemann problems at zone interfaces and a final evolution stage. We have checked and validated the code against several benchmarks available in literature. Test problems in 1, 2 and 3 dimensions are discussed.

Accretion disks and astrophysical jets are used to model many active astrophysical objects, viz., young stars, relativistic stars, and active galactic nuclei. In this paper we present self-consistent time-dependent simulations of supersonic jets launched from magnetized accretion disks, using high resolution numerical techniques. In particular we study the effects of the disk magnetic resistivity, parametrized through an alpha-prescription, in determining the properties of the inflow-outflow system. Moreover we analyze under which conditions steady state solutions of the type proposed in the self similar models of Blandford and Payne can be reached and maintained in a self consistent nonlinear stage. We use the resistive MHD FLASH code with adaptive mesh refinement, allowing us to follow the evolution of the structure for a time scale long enough to reach steady state. A detailed analysis of the initial configuration state is given. We obtain the expected solutions in the axisymmetric (2.5D) limit. Assuming a magnetic field around equipartition with the thermal pressure of the disk, we show how the characteristics of the disk jet system, as the ejection efficiency and the energetics, are affected by the anomalous resistivity acting inside the disk.

We derive an analytical solution to the computation of the output of a Lyot coronagraph for a given complex amplitude on the pupil plane. This solution, which does not require any simplifying assumption, relies on an expansion of the entrance complex amplitude on a Zernike base. According to this framework, the main contribution of the paper is the expression of the response of the coronagraph to a single base function. This result is illustrated by a computer simulation which describes the classical effect of propagation of a tip-tilt error in a coronagraph.

We perform a detailed analytical and numerical study of the dynamics of perturbations (vortex/aperiodic mode, Rossby and spiral-density waves) in 2D compressible disks with a Keplerian law of rotation. We draw attention to the process of spiral-density wave generation from vortices, discussing, in particular, the initial, most peculiar stages of wave emission. We show that the linear phenomenon of wave generation by vortices in smooth (without inflection points) shear flows found by using the so-called non-modal approach, is directly applicable to the present case. After an analytical non-modal description of the physics and characteristics of the spiral-density wave generation/propagation in the local shearing-sheet model, we follow the process of wave generation by small amplitude coherent circular vortex structures, by direct global numerical simulation, describing the main features of the generated waves.

X-Ray observations of groups and clusters of galaxies show that the Intra-Cluster Medium (ICM) in their cores is hotter than expected from cosmological numerical simulations of cluster formation which include star formation, radiative cooling and SN feedback. We investigate the effect of the injection of supersonic AGN jets into the ICM using axisymmetric hydrodynamical numerical simulations. A simple model for the ICM, describing the radial properties of gas and the gravitational potential in cosmological N-Body+SPH simulations of one cluster and three groups of galaxies at redshift z=0, is obtained and used to set the environment in which the jets are injected. We varied the kinetic power of the jet and the emission-weighted X-Ray temperature of the ICM. The jets transfer their energy to the ICM mainly by the effects of their terminal shocks. A high fraction of the injected energy can be deposited through irreversible processes in the cluster gas, up to 75% in our simulations. We show how one single, powerful jet can reconcile the predicted X-Ray properties of small groups, e.g. the Lx-Tx relation, with observations. We argue that the interaction between AGN jets and galaxy groups and cluster atmospheres is a viable feedback mechanism.

X-ray observations of extragalactic radiosources show strong evidences of interaction between the radio emitting plasma and the X-ray emitting ambient gas. In this paper we perform a detailed study of this interaction by numerical simulations. We study the propagation of an axisymmetric supersonic jet in an isothermal King atmosphere and we analyze the evolution of the resulting X-ray properties and their dependence on the jet physical parameters. We show the existence of two distinct and observationally different regimes of interaction, with strong and weak shocks. In the first case shells of enhanced X-ray emission are to be expected, while in the second case we expect deficit of X-ray emission coincident with the cocoon. By a comparison between analytical models and the results of our numerical simulations, we discuss the dependence of the transition between these two regimes on the jet parameters and we find that the mean controlling quantity results to be the jet kinetic power. We then discuss how the observed jets can be used to constrain the jet properties.

The 3-dimensional (3-D) calculation of the atmospheric neutrino flux by means of the FLUKA Monte Carlo model is here described in all details, starting from the latest data on primary cosmic ray spectra. The importance of a 3-D calculation and of its consequences have been already debated in a previous paper. Here instead the focus is on the absolute flux. We stress the relevant aspects of the hadronic interaction model of FLUKA in the atmospheric neutrino flux calculation. This model is constructed and maintained so to provide a high degree of accuracy in the description of particle production. The accuracy achieved in the comparison with data from accelerators and cross checked with data on particle production in atmosphere certifies the reliability of shower calculation in atmosphere. The results presented here can be already used for analysis by current experiments on atmospheric neutrinos. However they represent an intermediate step towards a final release, since this calculation does not yet include the bending of charged particles in atmosphere. On the other hand this last aspect, while requiring a considerable effort in a fully 3-D description of the Earth, if a high level of accuracy has to be maintained, does not affect in a significant way the analysis of atmospheric neutrino events.

Substantial fluxes of protons and leptons with energies below the geomagnetic cutoff have been measured by the AMS experiment at altitudes of 370-390 Km, in the latitude interval +/- 51.7 degrees. The production mechanisms of the observed trapped fluxes are investigated in detail by means of the FLUKA Monte Carlo simulation code. All known processes involved in the interaction of the cosmic rays with the atmosphere (detailed descriptions of the magnetic field and atmospheric density, as well as the electromagnetic and nuclear interaction processes) are included in the simulation. The results are presented and compared with the experimental data, indicating good agreement with the observed fluxes. The impact of secondary proton flux on particle production in atmosphere is briefly discussed.

We study, by numerical simulations, the entrainment process in a supersonic, radiative jet flow, during the evolution of Kelvin-Helmholtz instabilities, in the context of the the acceleration problem of molecular bipolar outflows, observed in Giant Molecular Clouds. Our results show that a large fraction of the initial jet momentum can be transferred to the ambient medium by this process. We therefore analyze in detail the instability evolution and compare some of the main observational properties of molecular outflows with those of the entrained material that we observe in our simulations. In particular, we find a good agreement for the mass vs. velocity distribution and for the outflow collimation structure, especially when a light jet is moving into a denser ambient medium. This is probably the case for (obscured) optical jets driving powerful molecular outflows in the denser environment of the inner regions of molecular clouds.

The analysis of the stability properties of astrophysical jets against Kelvin-Helmholtz (or shear-layer) instabilities plays a basic role in the understanding the origin and physical characteristics of these objects. Numerical simulations by Bodo et al. (1998) have shown that the three-dimensional non-linear evolution of KH instabilities in supersonic jets is substantially faster than in the two-dimensional case, leading to a cascade of modes towards smaller scales and a very effective mixing and momentum transfer to the ambient medium. On the other hand, Rossi et al. (1997) and Micono et al. (1998) found, in two dimensions, that radiative losses tend to reduce and delay mixing effects and momentum transfer to the ambient medium. In this paper, as a logical next step, we investigate the effects of radiative losses on the stability of 3D supersonic jets, assuming that the internal jet density is initially lower, equal and higher than the ambient medium, respectively. We find that light and equal-density radiative jets evolve in a qualitatively similar fashion with respect to the corresponding adiabatic ones. Conversely, we note substantial differences in the evolution of heavy jets: they remain more collimated and do not spread out, while the momentum gained by the ambient medium stays within ~ 5 jet radii.

Recent HST observations suggest that the NLR in Seyfert Galaxies can be the result of interaction between jet and external inhomogeneous medium; following this suggestion we perform numerical simulations considering the impact of a radiative jet on a dense cloud. We approach the problem adopting a hydrodynamical code, that consents us to study in detail the jet hydrodynamics, while we choose a more simplified treatment of radiative processes, in order to give a qualitatively good interpretation of the emission processes. Our three main purposes are: i) to reproduce in our simulations the physical conditions observed in the NLR of Seyfert Galaxies, ii) to obtain physical constraints of the jet parameters and iii) to study the jet capability to photoionize the surrounding medium.

We calculate the temporal evolution of distributions of relativistic electrons subject to synchrotron and adiabatic processes and Fermi-like acceleration in shocks. The shocks result from Kelvin-Helmholtz instabilities in the jet. Shock formation and particle acceleration are treated in a self-consistent way by means of a numerical hydrocode. We show that in our model the number of relativistic particles is conserved during the evolution, with no need of further injections of supra-thermal particles after the initial one. From our calculations, we derive predictions for values and trends of quantities like the spectral index and the cutoff frequency that can be compared with observations.

A detailed simulation of vertical showers in atmosphere produced by primary gammas and protons, in the energy range 1-100 TeV, has been performed by means of the FLUKA Monte Carlo code, with the aim of studying the time structure of the shower front at different detector heights. It turns out that the time delay distribution can be fitted using few parameters coincident with the distribution central moments. Such parameters exhibit a smooth behaviour as a function of energy. These results can be used both for detector design and for the interpretation of the existing measurements. Differences in the time structure between gamma and proton induced showers are found and explained in terms of the non-relativistic component of extensive air showers.

We present the results of a numerical analysis of the propagation and interaction of a supersonic jet with the external medium. We discuss the motion of the head of the jet into the ambient in different physical conditions, carrying out calculations with different Mach numbers and density ratios of the jet to the exteriors. Performing the calculation in a reference frame in motion with the jet head, we can follow in detail its long term dynamics. This numerical scheme allows us also to study the morphology of the cocoon for different physical parameters. We find that the propagation velocity of the jet head into the ambient medium strongly influences the morphology of the cocoon, and this result can be relevant in connection to the origin and structure of lobes in extragalactic radiosources.