We briefly review the various proposed scenarios that may lead to nonthermal radio emissions from exoplanetary systems (planetary magnetospheres, magnetosphere-ionosphere and magnetosphere-satellite coupling, and star-planet interactions), and the physical information that can be drawn from their detection. The latter scenario is especially favorable to the production of radio emission above 70\,MHz. We summarize the results of past and recent radio searches, and then discuss FAST characteristics and observation strategy, including synergies. We emphasize the importance of polarization measurements and a high duty-cycle for the very weak targets that radio-exoplanets prove to be.

All magnetized planets are known to produce intense non thermal radio emissions through a mechanism known as Cyclotron Maser Instability (CMI), requiring the presence of accelerated electrons generally arising from magnetospheric current systems. In return, radio emissions are a good probe of these current systems and acceleration processes. The CMI generates highly anisotropic emissions and leads to important visibility effects, which have to be taken into account when interpreting the data. Several studies showed that modeling the radio source anisotropic beaming pattern can reveal a wealth of physical information about the planetary or exoplanetary magnetospheres that produce these emissions. We present a numerical tool, called ExPRES (Exoplanetary and Planetary Radio Emission Simulator), which is able to reproduce the occurrence in time-frequency plane of CMI-generated radio emissions from planetary magnetospheres, exoplanets or star-planet interacting systems. Special attention is given to the computation of the radio emission beaming at and near its source. We explain what physical information about the system can be drawn from such radio observations, and how it is obtained. These information may include the location and dynamics of the radio sources, the type of current system leading to electron acceleration and their energy and, for exoplanetary systems, the magnetic field strength, the orbital period of the emitting body and the rotation period, tilt and offset of the planetary magnetic field. Most of these parameters can be remotely measured only via radio observations. The ExPRES code provides the proper framework of analysis and interpretation for past (Cassini, Voyager, Galileo), current (Juno, groundbased radiotelescopes) and future (BepiColombo, Juice) observations of planetary radio emissions, as well as for future detection of radio emissions from exoplanetary systems.

The high latitude radio emissions produced by the Cyclotron Maser Instability (CMI) in Jupiter's magnetosphere extend from a few kHz to 40 MHz. Part of the decametric emissions is of auroral origin, and part is driven by the moons Io, Europa and Ganymede. After summarizing the method used to identify Jupiter-satellite radio emissions, which consists in comparing space- and ground-based radio observations to ExPRES simulations of CMI-driven emissions in the time-frequency plane, we present a parametric study of the free parameters required by the ExPRES code (electron distribution function and resonant energy, magnetic field model, lead angle, and altitude of the ionospheric cut-off) in order to assess the accuracy of our simulations in the Io-Jupiter case. We find that Io-DAM arcs are fairly modeled by loss-cone driven CMI with electrons of 1-10 keV energy, using the ISaAC, VIPAL or VIP4 magnetic field model and a simple sinusoidal lead angle model. The altitude of the ionospheric cut-off has a marginal impact on the simulations. We discuss the impact of our results on the identification of Europa-DAM and Ganymede-DAM emissions.

Voyager/PRA (Planetary Radio Astronomy) data from digitized tapes archived at CNES have been reprocessed and recalibrated. The data cover the Jupiter and Saturn flybys of both Voyager probes. We have also reconstructed goniopolarimetric datasets (flux and polarization) at full resolution. These datasets are currently not available to the scientific community, but they are of primary interest for the analysis of the Cassini data at Saturn, and the Juno data at Jupiter, as well as for the preparation of the JUICE mission. We present the first results derived from the re-analysis of this dataset.

Accurately predicting structural properties of dark matter halos is one of the fundamental goals of modern cosmology. We use the new suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. The MultiDark simulations cover a large range of masses 1e10-1e15Msun and volumes upto 50Gpc**3. The total number of dark matter halos in all the simulations exceeds 60 billion. We find that in order to understand the structure of dark matter halos and to make ~1% accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the traditional ratio of the virial radius to the core radius in the NFW profile. For massive halos the averge density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter alpha in the Einasto approximation. Combining results from different redshifts, masses and cosmologies, we show that halos progress through three stages of evolution. (1) They start as rare density peaks that experience very fast and nearly radial infall. This radial infall brings mass closer to the center producing a high concentrated halo. Here, the halo concentration increases with the increasing halo mass and the concentration is defined by the alpha parameter with nearly constant core radius. Later halos slide into (2) the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. Now the concentration does not depend on halo mass. (3) Once the rate of accretion slows down, halos move into the domain of declining concentration-mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. We provide accurate analytical fits to the numerical results for halo density profiles and concentrations.

This work investigates the impact of cosmic flows and density perturbations on Hubble constant $H_0$ measurements using nonlinear phase-space reconstructions of the Local Universe (LU). In particular, we rely on 25 constrained N-body simulations of the LU using the 2MRS galaxy sample within distances of 90 Mpc/h. These have been randomly extended up to volumes enclosing 360 Mpc/h with augmented Lagrangian perturbation theory (750 simulations in total), accounting in this way for gravitational mode coupling from larger scales, correcting for periodic boundary effects, and estimating systematics of missing attractors ($\sigma_{\rm large}=134$ km/s). We report on Local Group speed reconstructions, which are compatible with CMB-dipole measurements: $|v_{\rm LG}|=685\pm137$ km/s. The direction $(l,b)=(260.5\pm 13.3,39.1\pm 10.4)^\circ$ is compatible with observations after considering the variance of large scales. Accounting for large scales, our local bulk flow estimations assuming a $\Lambda$CDM model are compatible with estimates based on velocity data derived from the Tully-Fisher relation. We focus on low redshift supernova measurements ($0.01<z<0.025$), which have been found to disagree with probes at larger distances. Our analysis indicates that there are two effects contributing to this tension. First, the anisotropic distribution of supernovae aligns with the velocity dipole and induces a systematic boost in $H_0$. Second, a divergent region surrounding the Virgo Supercluster is responsible for an additional positive bias in $H_0$. Taking these effects into account yields a correction of $\Delta H_0=-1.76 \pm 0.21 $ km/s/Mpc, thereby reducing the tension between local probes and more distant probes. Effectively $H_0$ is lower by about $2\%$.

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

We present, for the first time, a Local Universe (LU) characterization using high precision constrained $N$-body simulations based on self-consistent phase-space reconstructions of the large-scale structure in the Two-Micron All-Sky Galaxy Redshift Survey. We analyse whether we live in a special cosmic web environment by estimating cosmic variance from a set of unconstrained $\Lambda$CDM simulations as a function of distance to random observers. By computing volume and mass filling fractions for voids, sheets, filaments and knots, we find that the LU displays a typical scatter of about $1\sigma$ at scales $r>15\,h^{-1}\,{\rm Mpc}$, in agreement with $\Lambda$CDM, converging to a fair unbiased sample when considering spheres of about $60\,h^{-1}\,{\rm Mpc}$ radius. Additionally, we compute the matter density profile of the LU and found a reasonable agreement with the estimates of Karachentsev (2012) only when considering the contribution of dark haloes. This indicates that observational estimates may be biased towards low density values. As a first application of our reconstruction, we investigate the likelihood of different galaxy morphological types to inhabit certain cosmic web environments. In particular, we find that, irrespective of the method used to define the web, either based on the density or the peculiar velocity field, elliptical galaxies show a clear tendency to preferentially reside in clusters as opposed to voids (up to a level of $5.3\sigma$ and $9.8\sigma$ respectively) and conversely for spiral galaxies (up to a level of $5.6\sigma$ and $5.4\sigma$ respectively). These findings are compatible with previous works, albeit at higher confidence levels.

The increasing precision in the determination of the Hubble parameter has reached a per cent level at which large-scale cosmic flows induced by inhomogeneities of the matter distribution become non-negligible. Here we use large-scale cosmological N-body simulations to study statistical properties of the local Hubble parameter as measured by local observers. We show that the distribution of the local Hubble parameter depends not only on the scale of inhomogeneities, but also on how one defines the positions of observers in the cosmic web and what reference frame is used. Observers located in random dark matter haloes measure on average lower expansion rates than those at random positions in space or in the centres of cosmic voids, and this effect is stronger from the halo rest frames compared to the CMB rest frame. We compare the predictions for the local Hubble parameter with observational constraints based on type Ia supernovae (SNIa) and CMB observations. Due to cosmic variance, for observers located in random haloes we show that the Hubble constant determined from nearby SNIa may differ from that measured from the CMB by 0.8 per cent at 1sigma statistical significance. This scatter is too small to significantly alleviate a recently claimed discrepancy between current measurements assuming a flat LCDM model. However, for observers located in the centres of the largest voids permitted by the standard LCDM model, we find that Hubble constant measurements from SNIa would be biased high by 5 per cent, rendering this tension inexistent in this extreme case.

Accurate predictions of the abundance and clustering of dark matter haloes play a key role in testing the standard cosmological model. Here, we investigate the accuracy of one of the leading methods of connecting the simulated dark matter haloes with observed galaxies -- the Halo Abundance Matching (HAM) technique. We show how to choose the optimal values of the mass and force resolution in large-volume N-body simulations so that they provide accurate estimates for correlation functions and circular velocities for haloes and their subhaloes -- crucial ingredients of the HAM method. At the 10% accuracy, results converge for 50 particles for haloes and 150 particles for progenitors subhaloes. In order to achieve this level of accuracy a number of conditions should be satisfied. The force resolution for the smallest resolved (sub)haloes should be in the range (0.1-0.3)rs, where rs is the scale radius of (sub)haloes. The number of particles for progenitors of subhaloes should be 150. We also demonstrate that the two-body scattering plays a minor role for the accuracy of N-body simulations thanks to the relatively small number of crossing-times of dark matter in haloes, and the limited force resolution of cosmological simulations.

In the standard model of cosmology, structure emerges out of non-rotational flow and the angular momentum of collapsing halos is induced by tidal torques. The growth of halo angular momentum in the linear and quasi-linear phases is associated with a shear, curl-free, flow and it is well described within the linear framework of tidal torque theory (TTT). However, TTT is rendered irrelevant as haloes approach turn around and virialization. At that stage the flow field around halos has non-zero vorticity. Using a cosmological simulation, we have examined the importance of the curl of the velocity field (vorticity) in determining halo spin, finding a strong alignment between the two. We have also examined the alignment of vorticity with the principle axes of the shear tensor, finding that it tends to be perpendicular to the axis along which material is collapsing fastest (e1). This behavior is independent of halo masses and cosmic web environment. Our results agree with previous findings on the tendency of halo spin to be perpendicular to e1, and of the spin of (simulated) halos and (observed) galaxies to be aligned with the large-scale structure. Our results imply that angular momentum growth proceeds in two distinct phases. In the first phase angular momentum emerges out of a shear, curl-free, potential flow, as described by TTT. In the second phase, in which haloes approach virialization, the angular momentum emerges out of a vortical flow and halo spin becomes strongly aligned with the vorticity of the ambient flow field.

Lagrangian smoothed particle hydrodynamics (SPH) is a well-established approach to model fluids in astrophysical problems, thanks to its geometric flexibility and ability to automatically adjust the spatial resolution to the clumping of matter. However, a number of recent studies have emphasized inaccuracies of SPH in the treatment of fluid instabilities. The origin of these numerical problems can be traced back to spurious surface effects across contact discontinuities, and to SPH's inherent prevention of mixing at the particle level. We here investigate a new fluid particle model where the density estimate is carried out with the help of an auxiliary mesh constructed as the Voronoi tessellation of the simulation particles instead of an adaptive smoothing kernel. This Voronoi-based approach improves the ability of the scheme to represent sharp contact discontinuities. We show that this eliminates spurious surface tension effects present in SPH and that play a role in suppressing certain fluid instabilities. We find that the new `Voronoi Particle Hydrodynamics' described here produces comparable results than SPH in shocks, and better ones in turbulent regimes of pure hydrodynamical simulations. We also discuss formulations of the artificial viscosity needed in this scheme and how judiciously chosen correction forces can be derived in order to maintain a high degree of particle order and hence a regular Voronoi mesh. This is especially helpful in simulating self-gravitating fluids with existing gravity solvers used for N-body simulations.

We present optical and X-ray data for a sample of serendipitous XMM-Newton sources that are selected to have 0.5-2 keV vs 2-4.5 keV X-ray hardness ratios which are harder than the X-ray background. The sources have 2-4.5 keV X-ray flux >= 10^-14 cgs, and in this paper we examine a subsample of 42 optically bright (r < 21) sources; this subsample is 100 per cent spectroscopically identified. All but one of the optical counterparts are extragalactic, and we argue that the single exception, a Galactic M star, is probably a coincidental association. The X-ray spectra are consistent with heavily absorbed power laws (21.8 < log NH < 23.4), and all of them appear to be absorbed AGN. The majority of the sources show only narrow emission lines in their optical spectra, implying that they are type-2 AGN. Only a small fraction of the sources (7/42) show broad optical emission lines, and all of these have NH < 10^23 cm^-2. This implies that ratios of X-ray absorption to optical/UV extinction equivalent to > 100 times the Galactic gas-to-dust ratio are rare in AGN absorbers (at most a few percent of the population), and may be restricted to broad absorption-line QSOs. Seven objects appear to have an additional soft X-ray component in addition to the heavily absorbed power law. We consider the implications of our results in the light of the AGN unified scheme. We find that the soft components in narrow-line objects are consistent with the unified scheme provided that > 4 per cent of broad-line AGN have ionised absorbers that attenuate their soft X-ray flux by >50 per cent. In at least one of the X-ray absorbed, broad-line AGN in our sample the X-ray spectrum requires an ionised absorber, consistent with this picture.

We present a 23 ks XMM-Newton observation of the ULIRG Mrk 273. The hard X-ray spectrum can be modeled by a highly absorbed power law plus an Fe Kalpha emission line. The iron line is broad, suggesting possible superposition of a neutral iron line at 6.4 keV, and a blend of ionized iron lines from Fe XXV and Fe XXVI. Given the relatively short exposure, the three line components can not be singularly resolved with high statistical significance: the neutral component is detected at ~2.5sigma and the Fe XXV line at ~2sigma c.l., while for the Fe XXVI line we can only estimate an upper limit. The broad band spectrum requires, in addition to a highly absorbed power law, at least three collisionally ionized plasma components, which may be associated with star-forming regions. The temperatures of the three plasmas are about 0.3, 0.8 and 6 keV, where the highest of the three is sufficient to produce ionized iron emission lines. An alternative interpretation for the origin of the soft emission might also be given in terms of reflection off some photoionized gas, as has been observed in a number of nearby Compton-thick Seyfert 2 galaxies. A hot gas, photoionized by the primary, continuum can also produce ionized iron lines. Unfortunately, given the limited statistics and the lack of high resolution spectroscopy, it is not possible to distinguish between the two models investigated. We further compare the XMM-Newton findings with the Chandra data obtaining consistent spectral results. The absorption corrected hard X-ray (2-10 keV) luminosity of Mrk 273 corresponds to ~0.2% of the far-IR luminosity, similar to typical values found in pure starbursts. The thermal contribution to the soft X-ray luminosity is ~0.2-0.7x10^42 erg s^-1, comparable to those found in NGC 6240 and other starburst dominated ULIRGs.