We present a full-fledged analysis of Brans-Dicke cosmology with a cosmological constant and cold dark matter (BD-$\Lambda$CDM for short). We extend the scenarios where the current cosmological value of the BD-field is restricted by the local astrophysical domain to scenarios where that value is fixed only by the cosmological observations, which should be more natural in view of the possible existence of local screening mechanims. Our analysis includes both the background and perturbations equations in different gauges. We find that the BD-$\Lambda$CDM is favored by the overall cosmological data as compared to the concordance GR-$\Lambda$CDM model, namely data on distant supernovae, cosmic chronometers, local measurements of the Hubble parameter, baryonic acoustic oscillations, Large-Scale Structure formation and the cosmic microwave background under full Planck 2018 CMB likelihood. We also test the impact of Strong and Weak-Lensing data on our results, which can be significant. We find that the BD-$\Lambda$CDM can mimic effective quintessence with a significance of about $3-3.5\sigma$ c.l. (depending on the lensing datasets). The fact that the BD-$\Lambda$CDM behaves effectively as a Running Vacuum Model (RVM) when viewed from the GR perspective helps to alleviate some of the existing tensions with the data, such as the $\sigma_8$ excess predicted by GR-$\Lambda$CDM. On the other hand, the BD-$\Lambda$CDM model has a crucial bearing on the acute $H_0$-tension with the local measurements, which is rendered virtually harmless owing to the small increase of the effective value of the gravitational constant with the expansion. The simultaneous alleviation of the two tensions is a most remarkable feature of BD-gravity with a cosmological constant in the light of the current observations, and hence goes in support of BD-$\Lambda$CDM against GR-$\Lambda$CDM

In this work a weak-turbulence closure is used to determine the structure of the two-time power spectrum of weak magnetohydrodynamic (MHD) turbulence from the nonlinear equations describing the dynamics. The two-time energy spectrum is a fundamental quantity in turbulence theory from which most statistical properties of a homogeneous turbulent system can be derived. A closely related quantity, obtained via a spatial Fourier transform, is the two-point two-time correlation function describing the space-time correlations arising from the underlying dynamics of the turbulent fluctuations. Both quantities are central in fundamental turbulence theories as well as in the analysis of turbulence experiments and simulations. However, a first-principles derivation of these quantities has remained elusive due to the statistical closure problem, in which dynamical equations for correlations at order $n$ depend on correlations of order $n+1$. The recent launch of the Parker Solar Probe (PSP), which will explore the near-Sun region where the solar wind is born, has renewed the interest in the scientific community to understand the structure, and possible universal properties of space-time correlations. The weak MHD turbulence regime that we consider in this work allows for a natural asymptotic closure of the two-time spectrum, which may be applicable to other weak turbulence regimes found in fluids and plasmas. An integro-differential equation for the scale-dependent temporal correlation function is derived for incompressible Alfv\'enic fluctuations whose nonlinear dynamics is described by the reduced MHD equations.

The Eulerian space-time correlation of strong Magnetohydrodynamic (MHD) turbulence in strongly magnetized plasmas is investigated by means of direct numerical simulations of Reduced MHD turbulence and phenomenological modeling. Two new important results follow from the simulations: 1) counter-propagating Alfv\'enic fluctuations at a each scale decorrelate in time at the same rate in both balanced and imbalanced turbulence; and 2) the scaling with wavenumber of the decorrelation rate is consistent with pure hydrodynamic sweeping of small-scale structures by the fluctuating velocity of the energy-containing scales. An explanation of the simulation results is proposed in the context of a recent phenomenological MHD model introduced by Bourouaine and Perez 2019 (BP19) when restricted to the strong turbulence regime. The model predicts that the two-time power spectrum exhibits an universal, self-similar behavior that is solely determined by the probability distribution function of random velocities in the energy-containing range. Understanding the scale-dependent temporal evolution of the space-time turbulence correlation as well as its associated universal properties is essential in the analysis and interpretation of spacecraft observations, such as the recently launched Parker Solar Probe (PSP).

In classical mechanics, the Kepler potential and the Harmonic potential share the following remarkable property: in either of these potentials, a bound test particle orbits with a radial period that is independent of its angular momentum. For this reason, the Kepler and Harmonic potentials are called \it{isochrone}. In this paper, we solve the following general problem: are there any other isochrone potentials, and if so, what kind of orbits do they contain? To answer these questions, we adopt a geometrical point of view initiated by Michel H\'enon in 1959, in order to explore and classify exhaustively the set of isochrone potentials and isochrone orbits. In particular, we provide a geometric generalization of Kepler's third law, and give a similar law for the apsidal angle, for any isochrone orbit. We also relate the set of isochrone orbits to the set of parabolae in the plane under linear transformations, and use this to derive an analytical parameterization of any isochrone orbit. Along the way we compare our results to known ones, pinpoint some interesting details of this mathematical physics problem, and argue that our geometrical methods can be exported to more generic orbits in potential theory.

In this paper we apply a methodology, recently proposed by Bourouaine and Perez (2019) (BP19), to interpret solar-wind turbulent power spectra beyond Taylor approximation (TA). The turbulent power spectra were measured using \emph{Helios} spacecraft data near 0.6 au. We use the model proposed in BP19 to reproduce the field-perpendicular power spectrum $E(k_\perp)$ of anti-sunward Alfv\'enic fluctuations in the plasma frame (where $k_\perp$ is the field-perpendicular wavenumber) from the corresponding measured frequency power spectrum $P_{\rm sc}(\omega,\theta_b)$ along the sampling angle $\theta_b$, which is the angle between the local magnetic field and the sampling direction. Here $\omega=2\pi f$ and $f$ is the frequency of the time signal. Interestingly enough, we found that for all corresponding measured frequency power spectrum $P_{\rm sc}(\omega,\theta_b)$ the reproduced field-perpendicular power spectrum $E(k_\perp)$ is the same and independent of the considered sampling angle $\theta_b$. This finding is consistent with the fact that the analyzed turbulence is strong and highly anisotropic with $k_\| \ll k_\perp$ (where $k_\|$ is the field-parallel wavenumber). Furthermore, for this specific time signal we found that the commonly used TA is still approximately valid with the important difference that a broadening in $k_\perp$ for each angular frequency $\omega$ is present. This broadening can be described in the context of the methodology proposed in BP19.

Comet Interceptor is an ESA F-class mission expected to launch in 2028 on the same launcher as ESA's ARIEL mission. Comet Interceptor's science payload consists of three spacecraft, a primary spacecraft that will carry two smaller probes to be released at the target. The three spacecraft will fly-by the target along different chords, providing multiple simultaneous perspectives of the comet nucleus and its environment. Each of the spacecraft will be equipped with different but complementary instrument suites designed to study the far and near coma environment and surface of a comet or interstellar object (ISO). The primary spacecraft will perform a fly-by at ~1000 km from the target. The two smaller probes will travel deeper into the coma, closer to the nucleus. The mission is being designed and launched without a specific comet designated as its main target. Comet Interceptor will travel to the Sun-Earth L2 Lagrangian point with ARIEL and wait in hibernation until a suitable long-period comet (LPC) is found that will come close enough to the Sun for the spacecraft to maneuver to an encounter trajectory. To prepare for all eventualities, the science team has assembled a preliminary set of backup targets from the known Jupiter family comets, where a suitable fly-by trajectory can be achieved during the nominal mission timeline (including the possibility of some launch delay). To better prioritize this list, we are releasing our potential backup targets in order to solicit the planetary community's help with observations of these objects over future apparitions and to encourage publication of archival data on these objects.

We analyze Brans-Dicke gravity with a cosmological constant, $\Lambda$, and cold dark matter (BD-$\Lambda$CDM for short) in the light of the latest cosmological observations on distant supernovae, Hubble rate measurements at different redshifts, baryonic acoustic oscillations, large scale structure formation data, gravitational weak-lensing and the cosmic microwave background under full Planck 2015 CMB likelihood. Our analysis includes both the background and perturbations equations. We find that BD-$\Lambda$CDM is observationally favored as compared to the concordance $\Lambda$CDM model, which is traditionally defined within General Relativity (GR). In particular, some well-known persisting tensions of the $\Lambda$CDM with the data, such as the excess in the mass fluctuation amplitude $\sigma_8$ and specially the acute $H_0$-tension with the local measurements, essentially disappear in this context. Furthermore, viewed from the GR standpoint, BD-$\Lambda$CDM cosmology mimics quintessence at $\gtrsim3\sigma$ c.l. near our time.

We present 3D numerical simulations and an analytic model of reflection-driven MHD turbulence in the solar wind. Our simulations describe transverse, non-compressive MHD fluctuations within a narrow magnetic flux tube that extends from the photosphere out to a heliocentric distance $r$ of 21 solar radii $(R_s)$. We launch outward-propagating "$z^+$ fluctuations" into the simulation domain by imposing a randomly evolving photospheric velocity field. As these fluctuations propagate away from the Sun, they undergo partial reflection, producing inward-propagating "$z^-$ fluctuations." Counter-propagating fluctuations subsequently interact, causing fluctuation energy to cascade to small scales and dissipate. Our analytic model incorporates alignment, allows for strongly or weakly turbulent nonlinear interactions, and divides the $z^+$ fluctuations into two populations with different characteristic radial correlation lengths. The inertial-range power spectra in our simulations evolve toward a $k_\perp^{-3/2}$ scaling at $r>10 R_s$, where $k_\perp$ is the wave-vector component perpendicular to the background magnetic field. In two of our simulations, the $z^+$ power spectra are much flatter between the coronal base and $r \simeq 4 R_s$. We argue that these spectral scalings are caused by: (1) high-pass filtering in the upper chromosphere; (2) the anomalous coherence of inertial-range $z^-$ fluctuations in a reference frame propagating outwards with the $z^+$ fluctuations; and (3) the change in the sign of the radial derivative of the Alfv\'en speed at $r=r_m \simeq 1.7 R_s$, which disrupts this anomalous coherence between $r=r_m$ and $r\simeq 2r_m$. At $r>1.3 R_s$, the turbulent heating rate in our simulations is comparable to the heating rate in a previously developed solar-wind model that agreed with a number of observational constraints.

This paper presents a survey of microwave front-end receivers installed at radio telescopes throughout the World. This unprecedented analysis was conducted as part of a review of front-end developments for Italian radio telescopes, initiated by the Italian National Institute for Astrophysics in 2016. Fifteen international radio telescopes have been selected to be representative of the instrumentation used for radio astronomical observations in the frequency domain from 300 MHz to 116 GHz. A comprehensive description of the existing receivers is presented and their characteristics are compared and discussed. The observing performances of the complete receiving chains are also presented. An overview of on-going developments illustrates and anticipates future trends in front-end projects to meet the most ambitious scientific research goals.

In this letter we propose a practical methodology to interpret future Parker Solar Probe (PSP) turbulent time signals even when Taylor's hypothesis is not valid. By extending Kraichnan's sweeping model used in hydrodynamics we derive the Eulerian spacetime correlation function in magnetohydrodynamics (MHD) turbulence. It is shown that in MHD, the temporal decorrelation of small-scale fluctuations arises from a combination of hydrodynamic sweeping induced by large-scale fluid velocity $\delta u_0$ and by the Alfv\'enic propagation along the local magnetic field. The resulting temporal part of the space-time correlation function is used to determine the wavenumber range $\Delta k_\perp=[k_{\rm min},k_{\rm max}]$ of the turbulent fluctuations that contribute to the power of a given frequency $\omega$ of the time signal measured in the spacecraft frame. Our analysis also shows that the shape of frequency power spectrum $P_{sc}(\omega)$ of the time signal will follow the same power-law of the reduced power spectrum $E(k_\perp)\sim k^{-\alpha}_\perp$ in the plasma frame, where $\alpha$ is the spectral index. The proposed framework for the analysis of PSP time signals entirely relies on two simple dimensionless parameters that can be empirically obtained from PSP measurements, namely, $\epsilon=\delta u_0/\sqrt 2 V_\perp$ (where $V_\perp$ is the perpendicular velocity of PSP seen in the plasma frame) and the spectral index $\alpha$.

Investigations on dark energy (DE) are currently inconclusive about its time evolution. Hints of this possibility do however glow now and then in the horizon. Herein we assess the current status of dynamical dark energy (DDE) in the light of a large body of updated $SNIa+H(z)+BAO+LSS+CMB$ observations, using the full Planck 2015 CMB likelihood. The performance of the $\Lambda$CDM model (with equation of state $w=-1$ for $\Lambda$) is confronted with that of the general XCDM and CPL parametrizations, as well as with the traditional $\Phi$CDM model based on the scalar field potential $V\sim \Phi^{-\alpha}$. In particular, we gauge the impact of the bispectrum in the LSS and BAO parts, and show that the subset of $CMB+BAO+LSS$ observations may contain the bulk of the DDE signal. The departure from $w=-1$ is significant: roughly $2.6\sigma$ for XCDM and $2.9\sigma$ for $\Phi$CDM. In both cases the full Bayesian evidence is found to be positive even for a prior range of the DDE parameters extending over several standard deviations from the mean when the bispectrum is taken into account. Positive signs follow as well from the preliminary results of Planck 2018 data using the compressed CMB likelihood. As a bonus we also find that the $\sigma_8$-tension becomes reduced with DDE.

We consider the impact of dynamical dark energy (DDE) in the possible solution of the existing tensions in the $\Lambda$CDM. We test both interacting and non-interacting DE models with dark matter (DM). Among the former, the running vacuum model (RVM) interacting with DM appears as a favored option. The non-interacting scalar field model based on the potential $V\sim \phi^{-\alpha}$, and the generic XCDM parametrization, also provide consistent signs of DDE. The important novelty of our analysis with respect to the existing ones in the literature is that we use the matter bispectrum, together with the power spectrum. Using a complete and updated set of cosmological observations on $SNIa+BAO+H(z)+LSS+CMB$, we find that the crucial triad $BAO+LSS+CMB$ (i.e. baryonic acoustic oscillations, large scale structure formation data and the cosmic microwave background) provide the bulk of the signal. The bispectrum data is instrumental to get hold of the DDE signal, as our analysis shows. If the bispectrum is not included, the DDE signal could not be currently perceived at a significant confidence level.

This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.

In the potential theory, isochrony was introduced by Michel H\'enon in 1959 to characterize astrophysical observations of some globular clusters. Today, Michel Henon's isochrone potential is mainly used for his integrable property in numerical simulations, but is generally not really known. In a recent paper [29], we have presented new fundamental and theoretical results about isochrony that have particular importance in self-gravitating dynamics and which are detailed in this paper. In particular, new characterization of the isochrone state has been proposed which are investigated in order to analyze the product of the fast relaxation of a self-gravitating system. The general paradigm consists in considering that this product is a lowered isothermal sphere (King Model). By a detailed numerical study we show that this paradigm fails when the isochrone model succeeds in reproducing the quasi-equilibrium state obtained just after fast relaxation.

Brans-Dicke (BD) cosmology is reconsidered from an approach in which the model can be formulated in $\Lambda$CDM form, but at the expense of replacing the rigid cosmological constant $\Lambda$ with a dynamical quasivacuum component which depends on a small parameter $\epsilon$. The product of $\epsilon$ times the BD-parameter $\omega_{\rm BD}$ becomes exactly determined in terms of the ordinary cosmological parameters. The GR limit is recovered for $\omega_{\rm BD}\to\infty$ when $\epsilon\to 0$. We solve the background cosmology of the model as well as the perturbations equations. When fitted to the cosmological data, we find that the BD-cosmology, transcribed in such an effective GR form, appears to be competitive with the $\Lambda$CDM and emulates the running vacuum model.

We study a generic class of time-evolving vacuum models which can provide a better phenomenological account of the overall cosmological observations as compared to the $\Lambda$CDM. Among these models, the running vacuum model (RVM) appears to be the most motivated and favored one, at a confidence level of $\sim 3\sigma$. We further support these results by computing the Akaike and Bayesian information criteria. Our analysis also shows that we can extract fair signals of dynamical dark energy (DDE) by confronting the same set of data to the generic XCDM and CPL parametrizations. In all cases we confirm that the combined triad of modern observations on Baryonic Acoustic Oscillations, Large Scale Structure formation, and the Cosmic Microwave Background, provide the bulk of the signal sustaining a possible vacuum dynamics. In the absence of any of these three crucial data sources, the DDE signal can not be perceived at a significant confidence level. Its possible existence could be a cure for some of the tensions existing in the $\Lambda$CDM when confronted to observations.

Revisiting and extending an old idea of Michel H\'enon, we geometrically and algebraically characterize the whole set of isochrone potentials. Such potentials are fundamental in potential theory. They appear in spherically symmetrical systems formed by a large amount of charges (electrical or gravitational) of the same type considered in mean-field theory. Such potentials are defined by the fact that the radial period of a test charge in such potentials, provided that it exists, depends only on its energy and not on its angular momentum. Our characterization of the isochrone set is based on the action of a real affine subgroup on isochrone potentials related to parabolas in the $\mathbb{R}^2$ plane. Furthermore, any isochrone orbits are mapped onto associated Keplerian elliptic ones by a generalization of the Bohlin transformation. This mapping allows us to understand the isochrony property of a given potential as relative to the reference frame in which its parabola is represented. We detail this isochrone relativity in the special relativity formalism. We eventually exploit the completeness of our characterization and the relativity of isochrony to propose a deeper understanding of general symmetries such as Kepler's Third Law and Bertrand's theorem.

After about two decades of the first observational papers confirming the accelerated expansion of the universe, we are still facing the question whether the cause of it is a rigid cosmological constant $\Lambda$-term or a mildly evolving dynamical dark energy (DDE). While studies focusing mainly on CMB measurements do not perceive signs of physics beyond the $\Lambda$CDM, in this work we show that if we take a large string $SNIa+BAO+H(z)+LSS+CMB$ of modern cosmological observations, in which not only the CMB but also a rich sample of large scale structure formation data are included, one can extract $\sim 3.3\sigma$ signs of DDE using a simple XCDM parameterization. These signs can be enhanced up to near $3.8\sigma$ in the context of the running vacuum model (RVM), in which the vacuum energy density is in interaction with dark matter. Recently the RVM has been shown to provide an efficient and economical solution to the $\sigma_8$-tension, which is one of the intriguing phenomenological problems that has not been possible to solve within the $\Lambda$CDM so far. This fact contributes to strengthen the possibility that dynamical vacuum energy, or in general DDE, could be presently favored by the observations.

The Atacama Millimeter/submillimeter Array (ALMA) Phasing Project (APP) has developed and deployed the hardware and software necessary to coherently sum the signals of individual ALMA antennas and record the aggregate sum in Very Long Baseline Interferometry (VLBI) Data Exchange Format. These beamforming capabilities allow the ALMA array to collectively function as the equivalent of a single large aperture and participate in global VLBI arrays. The inclusion of phased ALMA in current VLBI networks operating at (sub)millimeter wavelengths provides an order of magnitude improvement in sensitivity, as well as enhancements in u-v coverage and north-south angular resolution. The availability of a phased ALMA enables a wide range of new ultra-high angular resolution science applications, including the resolution of supermassive black holes on event horizon scales and studies of the launch and collimation of astrophysical jets. It also provides a high-sensitivity aperture that may be used for investigations such as pulsar searches at high frequencies. This paper provides an overview of the ALMA Phasing System design, implementation, and performance characteristics.

In this year, in which we celebrate 100 years of the cosmological term, $\Lambda$, in Einstein's gravitational field equations, we are still facing the crucial question whether $\Lambda$ is truly a fundamental constant or a mildly evolving dynamical variable. After many theoretical attempts to understand the meaning of $\Lambda$, and in view of the enhanced accuracy of the cosmological observations, it seems now mandatory that this issue should be first settled empirically before further theoretical speculations on its ultimate nature. In this work, we summarize the situation of some of these studies. Devoted analyses made recently show that the $\Lambda=$const. hypothesis, despite being the simplest, may well not be the most favored one. The overall fit to the cosmological observables $SNIa+BAO+H(z)+LSS+CMB$ singles out the class RVM of the "running" vacuum models, in which $\Lambda=\Lambda(H)$ is an affine power-law function of the Hubble rate. It turns out that the performance of the RVM as compared to the "concordance" $\Lambda$CDM model (with $\Lambda=$const.) is much better. The evidence in support of the RVM may reach $\sim 4\sigma$ c.l., and is bolstered with Akaike and Bayesian criteria providing strong evidence in favor of the RVM option. We also address the implications of this framework on the tension between the CMB and local measurements of the current Hubble parameter.

Despite the outstanding achievements of modern cosmology, the classical dispute on the precise value of $H_0$, which is the first ever parameter of modern cosmology and one of the prime parameters in the field, still goes on and on after over half a century of measurements. Recently the dispute came to the spotlight with renewed strength owing to the significant tension (at $>3\sigma$ c.l.) between the latest Planck determination obtained from the CMB anisotropies and the local (distance ladder) measurement from the Hubble Space Telescope (HST), based on Cepheids. In this work, we investigate the impact of the running vacuum model (RVM) and related models on such a controversy. For the RVM, the vacuum energy density $\rho_{\Lambda}$ carries a mild dependence on the cosmic expansion rate, i.e. $\rho_{\Lambda}(H)$, which allows to ameliorate the fit quality to the overall $SNIa+BAO+H(z)+LSS+CMB$ cosmological observations as compared to the concordance $\Lambda$CDM model. By letting the RVM to deviate from the vacuum option, the equation of state $w=-1$ continues to be favored by the overall fit. Vacuum dynamics also predicts the following: i) the CMB range of values for $H_0$ is more favored than the local ones, and ii) smaller values for $\sigma_8(0)$. As a result, a better account for the LSS structure formation data is achieved as compared to the $\Lambda$CDM, which is based on a rigid (i.e. non-dynamical) $\Lambda$ term.

Recent analyses in the literature suggest that the concordance $\Lambda$CDM model with rigid cosmological term, $\Lambda=$const., may not be the best description of the cosmic acceleration. The class of "running vacuum models", in which $\Lambda=\Lambda(H)$ evolves with the Hubble rate, has been shown to fit the string of $SNIa+BAO+H(z)+LSS+CMB$ data significantly better than the $\Lambda$CDM. Here we provide further evidence on the time-evolving nature of the dark energy (DE) by fitting the same cosmological data in terms of scalar fields. As a representative model we use the original Peebles & Ratra potential, $V\propto\Phi^{-\alpha}$. We find clear signs of dynamical DE at $\sim 4\sigma$ c.l., thus reconfirming through a nontrivial scalar field approach the strong hints formerly found with other models and parametrizations.

The energy released in solar flares derives from a reconfiguration of magnetic fields to a lower energy state, and is manifested in several forms, including bulk kinetic energy of the coronal mass ejection, acceleration of electrons and ions, and enhanced thermal energy that is ultimately radiated away across the electromagnetic spectrum from optical to X-rays. Using an unprecedented set of coordinated observations, from a suite of instruments, we here report on a hitherto largely overlooked energy component -- the kinetic energy associated with small-scale turbulent mass motions. We show that the spatial location of, and timing of the peak in, turbulent kinetic energy together provide persuasive evidence that turbulent energy may play a key role in the transfer of energy in solar flares. Although the kinetic energy of turbulent motions accounts, at any given time, for only $\sim (0.5-1)$\% of the energy released, its relatively rapid ($\sim$$1-10$~s) energization and dissipation causes the associated throughput of energy (i.e., power) to rival that of major components of the released energy in solar flares, and thus presumably in other astrophysical acceleration sites.

Despite the fact that a rigid $\Lambda$-term is a fundamental building block of the concordance $\Lambda$CDM model, we show that a large class of cosmological scenarios with dynamical vacuum energy density $\rho_{\Lambda}$ and/or gravitational coupling $G$, together with a possible non-conservation of matter, are capable of seriously challenging the traditional phenomenological success of the $\Lambda$CDM. In this paper, we discuss these "running vacuum models" (RVM's), in which $\rho_{\Lambda}=\rho_{\Lambda}(H)$ consists of a nonvanishing constant term and a series of powers of the Hubble rate. Such generic structure is potentially linked to the quantum field theoretical description of the expanding Universe. By performing an overall fit to the cosmological observables $SNIa+BAO+H(z)+LSS+BBN+CMB$ (in which the WMAP9, Planck 2013 and Planck 2015 data are taken into account), we find that the class of RVM's appears significantly more favored than the $\Lambda$CDM, namely at an unprecedented level of $\gtrsim4.2\sigma$. Furthermore, the Akaike and Bayesian information criteria confirm that the dynamical RVM's are strongly preferred as compared to the conventional rigid $\Lambda$-picture of the cosmic evolution.

We report on electrostatic measurements made on board the European Space Agency mission LISA Pathfinder. Detailed measurements of the charge-induced electrostatic forces exerted on free-falling test masses (TMs) inside the capacitive gravitational reference sensor are the first made in a relevant environment for a space-based gravitational wave detector. Employing a combination of charge control and electric-field compensation, we show that the level of charge-induced acceleration noise on a single TM can be maintained at a level close to 1.0 fm/s^2/sqrt(Hz) across the 0.1-100 mHz frequency band that is crucial to an observatory such as LISA. Using dedicated measurements that detect these effects in the differential acceleration between the two test masses, we resolve the stochastic nature of the TM charge build up due to interplanetary cosmic rays and the TM charge-to-force coupling through stray electric fields in the sensor. All our measurements are in good agreement with predictions based on a relatively simple electrostatic model of the LISA Pathfinder instrument.

How can effective barotropic matter emerge from the interaction of cosmological fluids in an isotropic and homogeneous cosmological model ? The dynamics of homogeneous and isotropic Friedmann-Lema\^itre-Robertson-Walker universes is a natural special case of generalized Lotka-Volterra systems where each of the universe's fluid components can be seen as a competitive species in a predator-prey model. (Jungle universe : arXiv:1306.1037) In addition to numerical simulations illustrating this behaviour among the barotropic fluids filling the universe, we analytically pinpoint that effective time-dependent barotropic indices can arise from a physical coupling between those fluids whose dynamics could then look like that of another type of cosmic fluid, such as a cosmological constant. Since the nature of dark energy is still unknown, this dynamical approach could help understanding some of the properties of dark matter and dark energy at large cosmological scales.

OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.

In this work we study the process of energy dissipation triggered by a slow large scale motion of a magnetized conducting fluid. Our consideration is motivated by the problem of heating the solar corona, which is believed to be governed by fast reconnection events set off by the slow motion of magnetic field lines anchored in the photospheric plasma. To elucidate the physics governing the disruption of the imposed laminar motion and the energy transfer to small scales, we propose a simplified model where the large-scale motion of magnetic field lines is prescribed not at the footpoints but rather imposed volumetrically. As a result, the problem can be treated numerically with an efficient, highly-accurate spectral method, allowing us to use a resolution and statistical ensemble exceeding those of the previous work. We find that, even though the large-scale deformations are slow, they eventually lead to reconnection events that drive a turbulent state at smaller scales. The small-scale turbulence displays many of the universal features of field-guided MHD turbulence like a well developed inertial range spectrum. Based on these observations, we construct a phenomenological model that gives the scalings of the amplitude of the fluctuations and the energy dissipation rate as functions of the input parameters. We find a good agreement between the numerical results and the predictions of the model.

The interaction between Alfven-wave turbulence and the background solar wind affects the cross helicity in two ways. Non-WKB reflection converts outward-propagating Alfven waves into inward-propagating Alfven waves and vice versa, and the turbulence transfers momentum to the background flow. When both effects are accounted for, the total cross helicity is conserved. In the special case that the background density and flow speed are independent of time, the equations of cross-helicity conservation and total-energy conservation can be combined to recover a well-known equation derived by Heinemann and Olbert that has been interpreted as a non-WKB generalization of wave-action conservation. This latter equation (in contrast to cross-helicity and energy conservation) does not hold when the background varies in time.

Recently there have been claims on model-independent evidence of dynamical dark energy. Herein we consider a fairly general class of cosmological models with a time-evolving cosmological term of the form $\Lambda(H)=C_0+C_H H^2+C_{\dot{H}} \dot{H}$, where $H$ is the Hubble rate. These models are well motivated from the theoretical point of view since they can be related to the general form of the effective action of quantum field theory in curved spacetime. Consistency with matter conservation can be achieved by letting the Newtonian coupling $G$ change very slowly with the expansion. We solve these dynamical vacuum models and fit them to the wealth of expansion history and linear structure formation data. The results of our analysis show a significantly better agreement as compared to the concordance $\Lambda$CDM model, thus supporting the possibility of a dynamical cosmic vacuum.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

The Solar Probe Plus (SPP) spacecraft will explore the near-Sun environment, reaching heliocentric distances less than $10 R_{\odot}$. Near Earth, spacecraft measurements of fluctuating velocities and magnetic fields taken in the time domain are translated into information about the spatial structure of the solar wind via Taylor's "frozen turbulence" hypothesis. Near the perihelion of SPP, however, the solar-wind speed is comparable to the Alfv\'en speed, and Taylor's hypothesis in its usual form does not apply. In this paper, we show that, under certain assumptions, a modified version of Taylor's hypothesis can be recovered in the near-Sun region. We consider only the transverse, non-compressive component of the fluctuations at length scales exceeding the proton gyroradius, and we describe these fluctuations using an approximate theoretical framework developed by Heinemann and Olbert. We show that fluctuations propagating away from the Sun in the plasma frame obey a relation analogous to Taylor's hypothesis when $V_{\rm sc,\perp} \gg z^-$ and $z^+ \gg z^-$, where $V_{\rm sc,\perp}$ is the component of the spacecraft velocity perpendicular to the mean magnetic field and $\bm{z}^+$ ($\bm{z}^-$) is the Elsasser variable corresponding to transverse, non-compressive fluctuations propagating away from (towards) the Sun in the plasma frame. Observations and simulations suggest that, in the near-Sun solar wind, the above inequalities are satisfied and $\bm{z}^+$ fluctuations account for most of the fluctuation energy. The modified form of Taylor's hypothesis that we derive may thus make it possible to characterize the spatial structure of the energetically dominant component of the turbulence encountered by SPP.

The Large Magellanic Cloud, a satellite galaxy of the Milky Way, has been observed with the High Energy Stereoscopic System (H.E.S.S.) above an energy of 100 billion electron volts for a deep exposure of 210 hours. Three sources of different types were detected: the pulsar wind nebula of the most energetic pulsar known N 157B, the radio-loud supernova remnant N 132D and the largest non-thermal X-ray shell - the superbubble 30 Dor C. The unique object SN 1987A is, surprisingly, not detected, which constrains the theoretical framework of particle acceleration in very young supernova remnants. These detections reveal the most energetic tip of a gamma-ray source population in an external galaxy, and provide via 30 Dor C the unambiguous detection of gamma-ray emission from a superbubble.

Very high energy (VHE, $E>$100 GeV) $\gamma$-ray flaring activity of the high-frequency peaked BL Lac object \pg\ has been detected by the \hess\ telescopes. The flux of the source increased by a factor of 3 during the nights of 2012 April 26 and 27 with respect to the archival measurements with hint of intra-night variability. No counterpart of this event has been detected in the \fla\ data. This pattern is consistent with VHE $\gamma$ ray flaring being caused by the injection of ultrarelativistic particles, emitting $\gamma$ rays at the highest energies. The dataset offers a unique opportunity to constrain the redshift of this source at \bestz\ using a novel method based on Bayesian statistics. The indication of intra-night variability is used to introduce a novel method to probe for a possible Lorentz Invariance Violation (LIV), and to set limits on the energy scale at which Quantum Gravity (QG) effects causing LIV may arise. For the subluminal case, the derived limits are $\textrm{E}_{\rm QG,1}>4.10\times 10^{17}$ GeV and $\textrm{E}_{\rm QG,2}>2.10\times 10^{10}$ GeV for linear and quadratic LIV effects, respectively.

The LISA Pathfinder mission will demonstrate the technology of drag-free test masses for use as inertial references in future space-based gravitational wave detectors. To accomplish this, the Pathfinder spacecraft will perform drag-free flight about a test mass while measuring the acceleration of this primary test mass relative to a second reference test mass. Because the reference test mass is contained within the same spacecraft, it is necessary to apply forces on it to maintain its position and attitude relative to the spacecraft. These forces are a potential source of acceleration noise in the LISA Pathfinder system that are not present in the full LISA configuration. While LISA Pathfinder has been designed to meet it's primary mission requirements in the presence of this noise, recent estimates suggest that the on-orbit performance may be limited by this `suspension noise'. The drift-mode or free-flight experiments provide an opportunity to mitigate this noise source and further characterize the underlying disturbances that are of interest to the designers of LISA-like instruments. This article provides a high-level overview of these experiments and the methods under development to analyze the resulting data.

We propose an association between the phase-space mixing level of a self-gravitating system and the indistinguishability of its constituents (stars or dark matter particles). This represents a refinement in the study of systems exhibiting incomplete violent relaxation. Within a combinatorial analysis similar to that of Lynden-Bell, we make use of this association to obtain a distribution function that deviates from the Maxwell-Boltzmann distribution, increasing its slope for high energies. Considering the smallness of the occupation numbers for large distances from the center of the system, we apply a correction to Stirling's approximation which increases the distribution slope also for low energies. The distribution function thus obtained presents some resemblance to the "S" shape of distributions associated with cuspy density profiles (as compared to the distribution function obtained from the Einasto profile), although it is not quite able to produce sharp cusps. We also argue how the association between mixing level and indistinguishability can provide a physical meaning to the assumption of particle-permutation symmetry in the N-particle distribution function, when it is used to derive the one-particle Vlasov equation, which raises doubts about the validity of this equation during violent relaxation.

A survey of Michel Henon contributions to the study of globular cluster systems.

The region around the supernova remnant (SNR) W41 contains several TeV sources and has prompted the H.E.S.S. Collaboration to perform deep observations of this field of view. This resulted in the discovery of the new very high energy (VHE) source HESS J1832-093, at the position $\rm RA=18^h 32^m 50^s \pm 3^s_{stat} \pm 2^s_{syst}, \rm Dec=-9^\circ 22' 36'' \pm 32''_{stat} \pm 20''_{syst} (J2000)$, spatially coincident with a part of the radio shell of the neighboring remnant G22.7-0.2. The photon spectrum is well described by a power-law of index $\Gamma = 2.6 \pm 0.3_{\rm stat} \pm 0.1_{\rm syst}$ and a normalization at 1 TeV of $\Phi_0=(4.8 \pm 0.8_{\rm stat}\pm 1.0_{\rm syst})\,\times\,10^{-13}\,\rm{cm} ^{-2}\,s^{-1}\,TeV^{-1}$. The location of the gamma-ray emission on the edge of the SNR rim first suggested a signature of escaping cosmic-rays illuminating a nearby molecular cloud. Then a dedicated XMM-Newton observation led to the discovery of a new X-ray point source spatially coincident with the TeV excess. Two other scenarios were hence proposed to identify the nature of HESS J1832-093. Gamma-rays from inverse Compton radiation in the framework of a pulsar wind nebula scenario or the possibility of gamma-ray production within a binary system are therefore also considered. Deeper multi-wavelength observations will help to shed new light on this intriguing VHE source.

Dwarf spheroidal galaxies of the Local Group are close satellites of the Milky Way characterized by a large mass-to-light ratio and are not expected to be the site of non-thermal high-energy gamma-ray emission or intense star formation. Therefore they are amongst the most promising candidates for indirect dark matter searches. During the last years the High Energy Stereoscopic System (H.E.S.S.) of imaging atmospheric Cherenkov telescopes observed five of these dwarf galaxies for more than 140 hours in total, searching for TeV gamma-ray emission from annihilation of dark matter particles. The new results of the deep exposure of the Sagittarius dwarf spheroidal galaxy, the first observations of the Coma Berenices and Fornax dwarves and the re-analysis of two more dwarf spheroidal galaxies already published by the H.E.S.S. Collaboration, Carina and Sculptor, are presented. In the absence of a significant signal new constraints on the annihilation cross-section applicable to Weakly Interacting Massive Particles (WIMPs) are derived by combining the observations of the five dwarf galaxies. The combined exclusion limit depends on the WIMP mass and the best constraint is reached at 1-2 TeV masses with a cross-section upper bound of ~3.9x10-24 cm^3 s-1 at a 95% confidence level.

The inertial-interval energy spectrum of strong magnetohydrodynamic (MHD) turbulence with a uniform background magnetic field was observed numerically to be close to $k^{-3/2}$ by a number of independent groups. A dissenting opinion has been voiced by Beresnyak, A. 2011, PRL, 106, 075001-. 2012, MNRAS, 422, 3495-. 2014, ApJ, 784, L20 that the spectral scaling is close to $k^{-5/3}$. The conclusions of these papers are however incorrect as they are based on numerical simulations that are drastically unresolved, so that the discrete numerical scheme does not approximate the physical solution at the scales where the measurements are performed. These results have been rebutted in our more detailed papers Perez, J. C., Mason, J., Boldyrev, S., & Cattaneo, F. 2012, PRX, 2, 041005-. 2014, ApJL, 793, L13; here, by popular demand, we present a brief and simple explanation of our major criticism of Beresnyak's work.

A variety of studies of magnetised plasma turbulence invoke theories for the advection of a passive scalar by turbulent fluctuations. Examples include modelling the electron density fluctuations in the interstellar medium, understanding the chemical composition of galaxy clusters and the intergalactic medium, and testing the prevailing phenomenological theories of magnetohydrodynamic turbulence. While passive scalar turbulence has been extensively studied in the hydrodynamic case, its counterpart in MHD turbulence is significantly less well understood. Herein we conduct a series of high-resolution direct numerical simulations of incompressible, field-guided, MHD turbulence in order to establish the fundamental properties of passive scalar evolution. We study the scalar anisotropy, establish the scaling relation analogous to Yaglom's law, and measure the intermittency of the passive scalar statistics. We also assess to what extent the pseudo Alfven fluctuations in strong MHD turbulence can be modelled as a passive scalar. The results suggest that the dynamics of a passive scalar in MHD turbulence is considerably more complicated than in the hydrodynamic case.

We investigate the intermittency of energy dissipation in magnetohydrodynamic (MHD) turbulence by identifying dissipative structures and measuring their characteristic scales. We find that the probability distribution of energy dissipation rates exhibits a power law tail with index very close to the critical value of -2.0, which indicates that structures of all intensities contribute equally to energy dissipation. We find that energy dissipation is uniformly spread among coherent structures with lengths and widths in the inertial range. At the same time, these structures have thicknesses deep within the dissipative regime. As the Reynolds number is increased, structures become thinner and more numerous, while the energy dissipation continues to occur mainly in large-scale coherent structures. This implies that in the limit of high Reynolds number, energy dissipation occurs in thin, tightly packed current sheets which nevertheless span a continuum of scales up to the system size, exhibiting features of both coherent structures and nanoflares previously conjectured as a coronal heating mechanism.

Magnetohydrodynamic (MHD) turbulence in the majority of natural systems, including the interstellar medium, the solar corona, and the solar wind, has Reynolds numbers far exceeding the Reynolds numbers achievable in numerical experiments. Much attention is therefore drawn to the universal scaling properties of small-scale fluctuations, which can be reliably measured in the simulations and then extrapolated to astrophysical scales. However, in contrast with hydrodynamic turbulence, where the universal structure of the inertial and dissipation intervals is described by the Kolmogorov self-similarity, the scaling for MHD turbulence cannot be established based solely on dimensional arguments due to the presence of an intrinsic velocity scale -- the Alfven velocity. In this Letter, we demonstrate that the Kolmogorov first self-similarity hypothesis cannot be formulated for MHD turbulence in the same way it is formulated for the hydrodynamic case. Besides profound consequences for the analytical consideration, this also imposes stringent conditions on numerical studies of MHD turbulence. In contrast with the hydrodynamic case, the discretization scale in numerical simulations of MHD turbulence should decrease faster than the dissipation scale, in order for the simulations to remain resolved as the Reynolds number increases.

In this paper we report on the analysis of all the available optical and very high-energy $\gamma$-ray ($>$200 GeV) data for the BL Lac object PKS 2155$-$304, collected simultaneously with the ATOM and H.E.S.S. telescopes from 2007 until 2009. This study also includes X-ray (RXTE, Swift) and high-energy $\gamma$-ray (Fermi-LAT) data. During the period analysed, the source was transitioning from its flaring to quiescent optical states,and was characterized by only moderate flux changes at different wavelengths on the timescales of days and months. A flattening of the optical continuum with an increasing optical flux can be noted in the collected dataset, but only occasionally and only at higher flux levels. We did not find any universal relation between the very high-energy $\gamma$-ray and optical flux changes on the timescales from days and weeks up to several years. On the other hand, we noted that at higher flux levels the source can follow two distinct tracks in the optical flux-colour diagrams, which seem to be related to distinct $\gamma$-ray states of the blazar. The obtained results therefore indicate a complex scaling between the optical and $\gamma$-ray emission of PKS 2155$-$304, with different correlation patterns holding at different epochs, and a $\gamma$-ray flux depending on the combination of an optical flux and colour rather than a flux alone.

Previous observations with HESS have revealed the existence of an extended very-high-energy (VHE; E>100 GeV) gamma-ray source, HESS J1834-087, coincident with the SNR W41. The origin of the gamma-ray emission has been further investigated with HESS and the Fermi-LAT. The gamma-ray data provided by 61h (HESS) and 4 yrs (Fermi LAT) of observations cover over 5 decades in energy (1.8GeV - 30TeV). The morphology and spectrum of the TeV and GeV sources have been studied and multi-wavelength data have been used to investigate the origin of the observed emission. The TeV source can be modeled with a sum of two components: one point-like and one significantly extended (sig_TeV = 0.17{\deg}), both centered on SNR W41 and exhibiting spectra described by a power law of index 2.6. The GeV source detected with Fermi is extended (sig_GeV =0.15{\deg}) and morphologically matches the VHE emission. Its spectrum can be described by a power-law with index 2.15 and joins smoothly the one of the whole TeV source. A break appears in the spectra around 100 GeV. Two main scenarios are proposed to explain the emission: a pulsar wind nebula (PWN) or the interaction of SNR W41 with a molecular cloud. X-ray observations suggest the presence of a point-like source (pulsar candidate) near the center of the SNR and non-thermal X-ray diffuse emission which could arise from a potential PWN. The PWN scenario is supported by the match of of the TeV and GeV positions with the putative pulsar. However, the overall spectrum is reproduced by a 1-zone leptonic model only if an excess of low-energy electrons is injected by a high spin-down power pulsar. This low-energy component is not needed if the point-like TeV source is unrelated to the extended GeV and TeV sources. The interacting SNR scenario is supported by the spatial coincidence between the gamma-ray sources, the detection of OH maser lines and the hadronic modeling.

The long gamma-ray burst (GRB) 100621A, at the time the brightest X-ray transient ever detected by Swift-XRT in the $0.3\textrm{--}10$ keV range, has been observed with the H.E.S.S. imaging air Cherenkov telescope array, sensitive to gamma radiation in the very-high-energy (VHE, $>100$ GeV) regime. Due to its relatively small redshift of $z\sim0.5$, the favourable position in the southern sky and the relatively short follow-up time ($<700 \rm{s}$ after the satellite trigger) of the H.E.S.S. observations, this GRB could be within the sensitivity reach of the H.E.S.S. instrument. The analysis of the H.E.S.S. data shows no indication of emission and yields an integral flux upper limit above $\sim$380 GeV of $4.2\times10^{-12} \rm cm^{-2}s^{-1}$ (95 % confidence level), assuming a simple Band function extension model. A comparison to a spectral-temporal model, normalised to the prompt flux at sub-MeV energies, constraints the existence of a temporally extended and strong additional hard power law, as has been observed in the other bright X-ray GRB 130427A. A comparison between the H.E.S.S. upper limit and the contemporaneous energy output in X-rays constrains the ratio between the X-ray and VHE gamma-ray fluxes to be greater than 0.4. This value is an important quantity for modelling the afterglow and can constrain leptonic emission scenarios, where leptons are responsible for the X-ray emission and might produce VHE gamma rays.

The non-thermal nature of the X-ray emission from the shell-type supernova remnants (SNRs) G1.9+0.3 and G330.2+1.0 is an indication of intense particle acceleration in the shock fronts of both objects. This suggests that the SNRs are prime candidates for very-high-energy (VHE; E $>$ 0.1 TeV) {\gamma}-ray observations. G1.9+0.3, recently established as the youngest known SNR in the Galaxy, also offers a unique opportunity to study the earliest stages of SNR evolution in the VHE domain. The purpose of this work is to probe the level of VHE {\gamma}-ray emission from both SNRs and use this to constrain their physical properties. Observations were conducted with the H.E.S.S. (High Energy Stereoscopic System) Cherenkov telescope array over a more than six-year period spanning 2004-2010. The obtained data have effective livetimes of 67 h for G1.9+0.3 and 16 h for G330.2+1.0. The data are analyzed in the context of the multi-wavelength observations currently available and in the framework of both leptonic and hadronic particle acceleration scenarios. No significant {\gamma}-ray signal from G1.9+0.3 or G330.2+1.0 was detected. Upper limits (99% confidence level) to the TeV flux from G1.9+0.3 and G330.2+1.0 for the assumed spectral index {\Gamma} = 2.5 were set at 5.6 $\times$ 10$^{-13}$ cm$^{-2}$ s$^{-1}$ above 0.26 TeV and 3.2 $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ above 0.38 TeV, respectively. In a one-zone leptonic scenario, these upper limits imply lower limits on the interior magnetic field to B$_{\mathrm{G1.9}}$ $\gtrsim$ 11 {\mu}G for G1.9+0.3 and to B$_{\mathrm{G330}}$ $\gtrsim$ 8 {\mu}G for G330.2+1.0. In a hadronic scenario, the low ambient densities and the large distances to the SNRs result in very low predicted fluxes, for which the H.E.S.S. upper limits are not constraining.

Context. On March 4, 2013, the Fermi-LAT and AGILE reported a flare from the direction of the Crab Nebula in which the high-energy (HE; E > 100 MeV) flux was six times above its quiescent level. Simultaneous observations in other energy bands give us hints about the emission processes during the flare episode and the physics of pulsar wind nebulae in general. Aims. We search for variability of the emission of the Crab Nebula at very-high energies (VHE; E > 100 GeV), using contemporaneous data taken with the H.E.S.S. array of Cherenkov telescopes. Methods. Observational data taken with the H.E.S.S. instrument on five consecutive days during the flare were analysed concerning the flux and spectral shape of the emission from the Crab Nebula. Night-wise light curves are presented with energy thresholds of 1 TeV and 5 TeV. Results. The observations conducted with H.E.S.S. on 2013 March 6 to March 10 show no significant changes in the flux. They limit the variation on the integral flux above 1 TeV to less than 63% and the integral flux above 5 TeV to less than 78% at a 95% confidence level.

Composite supernova remnants (SNRs) constitute a small subclass of the remnants of massive stellar explosions where non-thermal radiation is observed from both the expanding shell-like shock front and from a pulsar wind nebula (PWN) located inside of the SNR. These systems represent a unique evolutionary phase of SNRs where observations in the radio, X-ray, and $\gamma$-ray regimes allow the study of the co-evolution of both these energetic phenomena. In this article, we report results from observations of the shell-type SNR G15.4+0.1 performed with the High Energy Stereoscopic System (H.E.S.S.) and XMM-Newton. A compact TeV $\gamma$-ray source, HESSJ1818-154, located in the center and contained within the shell of G15.4+0.1 is detected by H.E.S.S. and featurs a spectrum best represented by a power-law model with a spectral index of $-2.3 \pm 0.3_{stat} \pm 0.2_{sys}$ and an integral flux of F$(>$0.42$\,\mathrm{TeV}$)=($0.9 \pm 0.3_{\mathrm{stat}} \pm 0.2_{\mathrm{sys}}) \times 10^{-12}$cm$^{-2}$s$^{-1}$. Furthermore, a recent observation with XMM-Newton reveals extended X-ray emission strongly peaked in the center of G15.4+0.1. The X-ray source shows indications of an energy-dependent morphology featuring a compact core at energies above 4 keV and more extended emission that fills the entire region within the SNR at lower energies. Together, the X-ray and VHE $\gamma$-ray emission provide strong evidence of a PWN located inside the shell of G15.4+0.1, and this SNR can therefore be classified as a \emph{composite} based on these observations. The radio, X-ray, and $\gamma$-ray emission from the PWN is compatible with a one-zone leptonic model that requires a low average magnetic field inside the emission region. An unambiguous counterpart to the putative pulsar, which is thought to power the PWN, has been detected neither in radio nor in X-ray observations of G15.4+0.1.

The results of follow-up observations of the TeV gamma-ray source HESSJ 1640-465 from 2004 to 2011 with the High Energy Stereoscopic System (H.E.S.S.) are reported in this work. The spectrum is well described by an exponential cut-off power law with photon index Gamma=2.11 +/- 0.09_stat +/- 0.10_sys, and a cut-off energy of E_c = (6.0 +2.0 -1.2) TeV. The TeV emission is significantly extended and overlaps with the north-western part of the shell of the SNR G338.3-0.0. The new H.E.S.S. results, a re-analysis of archival XMM-Newton data, and multi-wavelength observations suggest that a significant part of the gamma-ray emission from HESS J1640-465 originates in the SNR shell. In a hadronic scenario, as suggested by the smooth connection of the GeV and TeV spectra, the product of total proton energy and mean target density could be as high as W_p n_H ~ 4 x 10^52 (d/10kpc)^2 erg cm^-3.

About 40% of the observation time of the High Energy Stereoscopic System (H.E.S.S.) is dedicated to studying active galactic nuclei (AGN), with the aim of increasing the sample of known extragalactic very-high-energy (VHE, E>100 GeV) sources and constraining the physical processes at play in potential emitters. H.E.S.S. observations of AGN, spanning a period from April 2004 to December 2011, are investigated to constrain their gamma-ray fluxes. Only the 47 sources without significant excess detected at the position of the targets are presented. Upper limits on VHE fluxes of the targets were computed and a search for variability was performed on the nightly time scale. For 41 objects, the flux upper limits we derived are the most constraining reported to date. These constraints at VHE are compared with the flux level expected from extrapolations of Fermi-LAT measurements in the two-year catalog of AGN. The H.E.S.S. upper limits are at least a factor of two lower than the extrapolated Fermi-LAT fluxes for 11 objects. Taking into account the attenuation by the extragalactic background light reduces the tension for all but two of them, suggesting intrinsic curvature in the high-energy spectra of these two AGN. Compilation efforts led by current VHE instruments are of critical importance for target-selection strategies before the advent of the Cherenkov Telescope Array, CTA.

Context: Very-high-energy (VHE; E>100 GeV) {\gamma}-ray emission from blazars inevitably gives rise to electron-positron pair production through the interaction of these {\gamma}-rays with the Extragalactic Background Light (EBL). Depending on the magnetic fields in the proximity of the source, the cascade initiated from pair production can result in either an isotropic halo around an initially beamed source or a magnetically broadened cascade flux. Aims: Both extended pair halo (PH) and magnetically broadened cascade (MBC) emission from regions surrounding the blazars 1ES 1101-232, 1ES 0229+200 and PKS 2155-304 were searched for, using VHE {\gamma}-ray data taken with the High Energy Stereoscopic System (H.E.S.S.), and high energy (HE; 100 MeV<E<100 GeV) {\gamma}-ray data with the Fermi Large Area Telescope (LAT). Methods: By comparing the angular distributions of the reconstructed gamma-ray events to the angular profiles calculated from detailed theoretical models, the presence of PH and MBC was investigated. Results: Upper limits on the extended emission around 1ES 1101-232, 1ES 0229+200 and PKS 2155-304 are found to be at a level of few percent of the Crab nebula flux above 1 TeV, depending on the assumed photon index of the cascade emission. Assuming strong Extra-Galactic Magnetic Field (EGMF) values, > 10$^{-12}$G, this limits the production of pair halos developing from electromagnetic cascades. For weaker magnetic fields, in which electromagnetic cascades would result in magnetically broadened cascades, EGMF strengths in the range (0.3 - 3)$\times 10^{-15}$G were excluded for PKS 2155-304 at the 99% confidence level, under the assumption of a 1 Mpc coherence length.

In this paper, we exploit the fact that the dynamics of homogeneous and isotropic Friedmann-Lemaitre universes is a special case of generalized Lotka-Volterra system where the competitive species are the barotropic fluids filling the Universe. Without coupling between those fluids, Lotka-Volterra formulation offers a pedagogical and simple way to interpret usual Friedmann-Lemaitre cosmological dynamics. A natural and physical coupling between cosmological fluids is proposed which preserve the structure of the dynamical equations. Using the standard tools of Lotka-Volterra dynamics, we obtain the general Lyapunov function of the system when one of the fluids is coupled to dark energy. This provides in a rigorous form a generic asymptotic behavior for cosmic expansion in presence of coupled species, beyond the standard de Sitter, Einstein-de Sitter and Milne cosmologies. Finally, we conjecture that chaos can appear for at least four interacting fluids.

HESS J0632+057 is the only gamma-ray binary known so far whose position in the sky allows observations with ground-based observatories both in the northern and southern hemispheres. Here we report on long-term observations of HESS J0632+057 conducted with the VERITAS and H.E.S.S. Cherenkov Telescopes and the X-ray Satellite Swift, spanning a time range from 2004 to 2012 and covering most of the system's orbit. The VHE emission is found to be variable, and is correlated with that at X-ray energies. An orbital period of $315 ^{+6}_{-4}$ days is derived from the X-ray data set, which is compatible with previous results, $P = (321 \pm 5$) days. The VHE light curve shows a distinct maximum at orbital phases close to 0.3, or about 100 days after periastron passage, which coincides with the periodic enhancement of the X-ray emission. Furthermore, the analysis of the TeV data shows for the first time a statistically significant ($> 6.5 \sigma$) detection at orbital phases 0.6--0.9. The obtained gamma-ray and X-ray light curves and the correlation of the source emission at these two energy bands are discussed in the context of the recent ephemeris obtained for the system. Our results are compared to those reported for other gamma-ray binaries.

Axionlike particles (ALPs) are hypothetical light (sub-eV) bosons predicted in some extensions of the Standard Model of particle physics. In astrophysical environments comprising high-energy gamma rays and turbulent magnetic fields, the existence of ALPs can modify the energy spectrum of the gamma rays for a sufficiently large coupling between ALPs and photons. This modification would take the form of an irregular behavior of the energy spectrum in a limited energy range. Data from the H.E.S.S. observations of the distant BL Lac object PKS 2155-304 (z = 0.116) are used to derive upper limits at the 95% C.L. on the strength of the ALP coupling to photons, $g_{\gamma a} < 2.1\times 10^{-11}$ GeV$^{-1}$ for an ALP mass between 15 neV and 60 neV. The results depend on assumptions on the magnetic field around the source, which are chosen conservatively. The derived constraints apply to both light pseudoscalar and scalar bosons that couple to the electromagnetic field.

The active galactic nucleus PKS 0301-243 (z=0.266) is a high-synchrotron-peaked BL Lac object that is detected at high energies (HE, 100 MeV < E < 100 GeV) by Fermi/LAT. This paper reports on the discovery of PKS 0301-243 at very high energies (E>100 GeV) by the High Energy Stereoscopic System (H.E.S.S.) from observations between September 2009 and December 2011 for a total live time of 34.9 hours. Gamma rays above 200 GeV are detected at a significance of 9.4{\sigma}. A hint of variability at the 2.5{\sigma} level is found. An integral flux I(E > 200 GeV) = (3.3 +/- 1.1_stat +/- 0.7_syst)e-12 ph cm^-2s^-1 and a photon index {\Gamma} = 4.6 +/- 0.7_stat +/- 0.2_syst are measured. Multi-wavelength light curves in HE, X-ray and optical bands show strong variability, and a minimal variability timescale of eight days is estimated from the optical light curve. A single-zone leptonic synchrotron self-Compton scenario satisfactorily reproduces the multi-wavelength data. In this model, the emitting region is out of equipartition and the jet is particle dominated. Because of its high redshift compared to other sources observed at TeV energies, the very high energy emission from PKS 0301-243 is attenuated by the extragalactic background light (EBL) and the measured spectrum is used to derive an upper limit on the opacity of the EBL.

We present direct numerical simulations of inhomogeneous reduced magnetohydrodynamic (RMHD) turbulence between the Sun and the Alfv\'en critical point. These are the first such simulations that take into account the solar-wind outflow velocity and the radial inhomogeneity of the background solar wind without approximating the nonlinear terms in the governing equations. RMHD turbulence is driven by outward-propagating Alfv\'en waves ($z^+$ fluctuations) launched from the Sun, which undergo partial non-WKB reflection to produce sunward-propagating Alfv\'en waves ($z^-$ fluctuations). We present ten simulations with different values of the correlation time $\tau_{\rm c\,\sun}^+$ and perpendicular correlation length $L_{\perp \sun}$ of outward-propagating Alfv\'en waves (AWs) at the coronal base. We find that between 15% and 33% of the $z^+$ energy launched into the corona dissipates between the coronal base and Alfv\'en critical point. Between 33% and 40% of this input energy goes into work on the solar-wind outflow, and between 22% and 36% escapes as $z^+$ fluctuations through the simulation boundary at $r=r_{\rm A}$. The $z^\pm$ power spectra scale like $k_\perp^{-\alpha^\pm}$, where $k_\perp$ is the wavenumber in the plane perpendicular to $\vec{B}_0$. In our simulation with the smallest value of $\tau_{\rm c\,\sun}^+$ ($\sim 2 \mbox{min}$) and largest value of $L_{\perp \sun}$ ($2\times 10^4 \mbox{km}$), we find that $\alpha^+$ decreases approximately linearly with increasing $\ln(r)$, reaching a value of 1.3 at $r=11.1 R_{\sun}$. Our simulations with larger values of $\tau_{\rm c\,\sun}^+$ exhibit alignment between the contours of constant $\phi^+$, $\phi^-$, $\Omega_0^+$, and $\Omega_0^-$, where $\phi^\pm$ are the Els\"asser potentials and $\Omega_0^\pm$ are the outer-scale parallel Els\"asser vorticities.

Recent theoretical studies argue that the rate of stochastic ion heating in low-frequency Alfv\'en-wave turbulence is given by $Q_\perp = c_1 [(\delta u)^3 /\rho] \exp(-c_2/\epsilon)$, where $\delta u$ is the rms turbulent velocity at the scale of the ion gyroradius $\rho$, $\epsilon = \delta u/v_{\perp \rm i}$, $v_{\perp \rm i}$ is the perpendicular ion thermal speed, and $c_1$ and $c_2$ are dimensionless constants. We test this theoretical result by numerically simulating test particles interacting with strong reduced magnetohydrodynamic (RMHD) turbulence. The heating rates in our simulations are well fit by this formula. The best-fit values of $c_1$ are $\sim 1$. The best-fit values of $c_2$ decrease (i.e., stochastic heating becomes more effective) as the grid size and Reynolds number of the RMHD simulations increase. As an example, in a $1024^2 \times 256$ RMHD simulation with a dissipation wavenumber of order the inverse ion gyroradius, we find $c_2 = 0.21$. We show that stochastic heating is significantly stronger in strong RMHD turbulence than in a field of randomly phased Alfv\'en waves with the same power spectrum, because coherent structures in strong RMHD turbulence increase orbit stochasticity in the regions where ions are heated most strongly. We find that $c_1$ increases by a factor of $\sim 3$ while $c_2$ changes very little as the ion thermal speed increases from values $\ll v_{\rm A}$ to values $\sim v_{\rm A}$, where $v_{\rm A}$ is the Alfv\'en speed. We discuss the importance of these results for perpendicular ion heating in the solar wind.

We present a set of hydrodynamical/Nbody controlled simulations of isolated gas rich galaxies that self-consistently include SN feedback and a detailed chemical evolution model, both tested in cosmological simulations. The initial conditions are motivated by the observed star forming galaxies at z ~ 2-3. We find that the presence of a multiphase interstellar media in our models promotes the growth of disc instability favouring the formation of clumps which in general, are not easily disrupted on timescales compared to the migration time. We show that stellar clumps migrate towards the central region and contribute to form a classical-like bulge with a Sersic index, n > 2. Our physically-motivated Supernova feedback has a mild influence on clump survival and evolution, partially limiting the mass growth of clumps as the energy released per Supernova event is increased, with the consequent flattening of the bulge profile. This regulation does not prevent the building of a classical-like bulge even for the most energetic feedback tested. Our Supernova feedback model is able to establish a self-regulated star formation, producing mass-loaded outflows and stellar age spreads comparable to observations. We find that the bulge formation by clumps may coexit with other channels of bulge assembly such as bar and mergers. Our results suggest that galactic bulges could be interpreted as composite systems with structural components and stellar populations storing archaeological information of the dynamical history of their galaxy.

LISA Pathfinder is a mission planned by the European Space Agency to test the key technologies that will allow the detection of gravitational waves in space. The instrument on-board, the LISA Technology package, will undergo an exhaustive campaign of calibrations and noise characterisation campaigns in order to fully describe the noise model. Data analysis plays an important role in the mission and for that reason the data analysis team has been developing a toolbox which contains all the functionalities required during operations. In this contribution we give an overview of recent activities, focusing on the improvements in the modelling of the instrument and in the data analysis campaigns performed both with real and simulated data.

A deep observation campaign carried out by the High Energy Stereoscopic System (H.E.S.S.) on Centaurus A enabled the discovery of gamma rays from the blazar 1ES 1312-423, two degrees away from the radio galaxy. With a differential flux at 1 TeV of (1.9 +/-0.6(stat) +/-0.4(sys)) x 10^{-13} /cm^2 /s /TeV corresponding to 0.5% of the Crab nebula differential flux and a spectral index of 2.9 +/- 0.5 (stat) +/- 0.2 (sys), 1ES 1312-423 is one of the faintest sources ever detected in the very high energy (E>100 GeV) extragalactic sky. A careful analysis using three and a half years of Fermi-LAT data allows the discovery at high energies (E>100 MeV) of a hard spectrum (index of 1.4 +/- 0.4 (stat) +/- 0.2 (sys)) source coincident with 1ES 1312-423. Radio, optical, UV and X-ray observations complete the spectral energy distribution of this blazar, now covering 16 decades in energy. The emission is successfully fitted with a synchrotron self Compton model for the non-thermal component, combined with a black-body spectrum for the optical emission from the host galaxy.

The detection of the high-frequency peaked BL Lac object (HBL) SHBL J001355.9-185406 ($z$=0.095) at high (HE; 100 MeV$<$E$<$300 GeV) and very high-energy (VHE; $E>100\,{\rm GeV}$) with the \fer\ Large Area Telescope (LAT) and the High Energy Stereoscopic System (H.E.S.S.) is reported. Dedicated observations have been performed with the H.E.S.S. telescopes, leading to a detection at the $5.5\,\sigma$ significance level. The measured flux above 310 GeV is $(8.3 \pm 1.7_{\rm{stat}}\pm 1.7_{\rm{sys}})\times 10^{-13}$ photons \cms\ (about 0.6% of that of the Crab Nebula), and the power law spectrum has a photon index of \indexHESS. Using 3.5 years of publicly available \fla\ data, a faint counterpart has been detected in the LAT data at the $5.5\,\sigma$ significance level, with an integrated flux above 300 MeV of $(9.3 \pm 3.4_{\rm stat} \pm 0.8_{\rm sys})\times 10^{-10}$ photons \cms\ and a photon index of $\Gamma = 1.96 \pm 0.20_{\rm stat} \pm 0.08_{\rm sys}$. X-ray observations with \textit{Swift}-XRT allow the synchrotron peak energy in $\nu F_\nu$ representation to be located at $\sim 1.0\,{\rm keV}$. The broadband spectral energy distribution is modelled with a one-zone synchrotron self-Compton (SSC) model and the optical data by a black-body emission describing the thermal emission of the host galaxy. The derived parameters are typical for HBLs detected at VHE, with a particle dominated jet.

The nature of subproton scale fluctuations in the solar wind is an open question, partly because two similar types of electromagnetic turbulence can occur: kinetic Alfven turbulence and whistler turbulence. These two possibilities, however, have one key qualitative difference: whistler turbulence, unlike kinetic Alfven turbulence, has negligible power in density fluctuations. In this Letter, we present new observational data, as well as analytical and numerical results, to investigate this difference. The results show, for the first time, that the fluctuations well below the proton scale are predominantly kinetic Alfven turbulence, and, if present at all, the whistler fluctuations make up only a small fraction of the total energy.

Observations reveal nearly power-law spectra of magnetic and density plasma fluctuations at subproton scales in the solar wind, which indicates the presence of a turbulent cascade. We discuss the three-field and two-field models for micro-scale plasma fluctuations, and then present the results of numerical simulations of a two-field model of kinetic-Alfven turbulence, which models plasma motion at sub-proton scales.

The quasar PKS 1510-089 (z=0.361) was observed with the H.E.S.S. array of imaging atmospheric Cherenkov telescopes during high states in the optical and GeV bands, to search for very high energy (VHE, defined as E >= 0.1 TeV) emission. VHE \gamma-rays were detected with a statistical significance of 9.2 standard deviations in 15.8 hours of H.E.S.S. data taken during March and April 2009. A VHE integral flux of I(0.15 TeV < E < 1.0 TeV) = (1.0 +- 0.2 (stat) +- 0.2 (sys) x 10^{-11} cm^{-2}s^{-1} is measured. The best-fit power law to the VHE data has a photon index of \Gamma=5.4 +- 0.7 (stat) +- 0.3 (sys). The GeV and optical light curves show pronounced variability during the period of H.E.S.S. observations. However, there is insufficient evidence to claim statistically significant variability in the VHE data. Because of its relatively high redshift, the VHE flux from PKS 1510-089 should suffer considerable attenuation in the intergalactic space due to the extragalactic background light (EBL). Hence, the measured \gamma-ray spectrum is used to derive upper limits on the opacity due to EBL, which are found to be comparable with the previously derived limits from relatively-nearby BL Lac objects. Unlike typical VHE-detected blazars where the broadband spectrum is dominated by non-thermal radiation at all wavelengths, the quasar PKS 1510-089 has a bright thermal component in the optical to UV frequency band. Among all VHE detected blazars, PKS 1510-089 has the most luminous broad line region (BLR). The detection of VHE emission from this quasar indicates a low level of \gamma-\gamma absorption on the internal optical to UV photon field.

Very high-energy gamma-ray emission from PKS 0447-439 was detected with the H.E.S.S. Cherenkov telescope array in December 2009. This blazar is one of the brightest extragalactic objects in the Fermi Bright Source List and has a hard spectrum in the MeV to GeV range. In the TeV range, a photon index of 3.89 +- 0.37 (stat) +- 0.22 (sys) and a flux normalisation at 1 TeV, Phi(1 TeV) = (3.5 +- 1.1 (stat) +- 0.9 (sys)) x 10^{-13} cm^{-2} s^{-1} TeV^{-1}, were found. The detection with H.E.S.S. triggered observations in the X-ray band with the Swift and RXTE telescopes. Simultaneous UV and optical data from Swift UVOT and data from the optical telescopes ATOM and ROTSE are also available. The spectrum and light curve measured with H.E.S.S. are presented and compared to the multi-wavelength data at lower energies. A rapid flare is seen in the Swift XRT and RXTE data, together with a flux variation in the UV band, at a time scale of the order of one day. A firm upper limit of z < 0.59 on the redshift of PKS 0447-439 is derived from the combined Fermi-LAT and H.E.S.S. data, given the assumptions that there is no upturn in the intrinsic spectrum above the Fermi-LAT energy range and that absorption on the Extragalactic Background Light (EBL) is not weaker than the lower limit provided by current models. The spectral energy distribution is well described by a simple one-zone Synchrotron Self-Compton (SSC) scenario, if the redshift of the source is less than z <~ 0.4.

We develop a framework for studying the statistical properties of current sheets in numerical simulations of 3D magnetohydrodynamic (MHD) turbulence. We describe an algorithm that identifies current sheets in a simulation snapshot and then determines their geometrical properties (including length, width, and thickness) and intensities (peak current density and total energy dissipation rate). We then apply this procedure to simulations of reduced MHD turbulence and perform a statistical analysis on the obtained population of current sheets. We evaluate the role of reconnection by separately studying the populations of current sheets which contain magnetic X-points and those which do not. We find that the statistical properties of the two populations are different in general. We compare the scaling of these properties to phenomenological predictions obtained for the inertial range of MHD turbulence. Finally, we test whether the reconnecting current sheets are consistent with the Sweet-Parker model.

Aim. In this paper we present very high energy (VHE; E>100 GeV) data from the \gamma-ray binary system PSR B1259-63/LS 2883 taken around its periastron passage (15th of December 2010) with the High Energy Stereoscopic System (H.E.S.S.) of Cherenkov Telescopes. We aim to search for a possible TeV counterpart of the GeV flare detected by the Fermi LAT. In addition, we aim to study the current periastron passage in the context of previous observations taken at similar orbital phases, testing the repetitive behavior of the source. Methods. Observations at VHE were conducted with H.E.S.S. from 9th to 16th of January 2011. The total dataset amounts to around 6 h of observing time. Results. The source is detected in the 2011 data at a significance level of 11.5\sigma\ revealing an averaged integral flux above 1 TeV of (1.01 \pm 0.18_{stat} \pm 0.20_{sys}) \times 10^{-12} cm^{-2}s^{-1}. The differential energy spectrum follows a power-law shape with a spectral index \Gamma = 2.92 \pm 0.30_{stat} \pm 0.20_{sys} and a flux normalisation at 1 TeV of N_{0} = 1.95 \pm 0.32_{stat} \pm 0.39_{sys}) \times 10^{-12} TeV^{-1} cm^{-2} s^{-1}. The measured lightcurve does not show any evidence for variability of the source on the daily scale. Conclusions. The measured integral flux and the spectral shape of the 2011 data are compatible with the results obtained around previous periastron passages. The absence of variability in the H.E.S.S. data indicates that the GeV flare observed by Fermi LAT in the time period covered also by H.E.S.S. observations originates in a different physical scenario than the TeV emission. Additionaly, new results compared to those obtained in the observations which were performed in 2004 at a similar orbital phase, further support the hypothesis of the repetitive behavior of the source.

Globular clusters (GCs) are established emitters of high-energy (HE, 100 MeV<E<100 GeV) \gamma-ray radiation which could originate from the cumulative emission of the numerous millisecond pulsars (msPSRs) in the clusters' cores or from inverse Compton (IC) scattering of relativistic leptons accelerated in the GC environment. These stellar clusters could also constitute a new class of sources in the very-high-energy (VHE, E>100 GeV) \gamma-ray regime, judging from the recent detection of a signal from the direction of Terzan 5 with the H.E.S.S. telescope array. We searched for point-like and extended VHE \gamma-ray emission from 15 GCs serendipitously covered by H.E.S.S observations and also performed a stacking analysis combining the data from all GCs to investigate the hypothesis of a population of faint emitters. Assuming IC emission as the origin of the VHE \gamma-ray signal from the direction of Terzan 5, we calculated the expected \gamma-ray flux from each of the 15 GCs, based on their number of millisecond pulsars, their optical brightness and the energy density of background photon fields. We did not detect significant VHE \gamma-ray emission from any of the 15 GCs in either of the two analyses. Given the uncertainties related to the parameter determinations, the obtained flux upper limits allow to rule out the simple IC/msPSR scaling model for NGC 6388 and NGC 7078. The upper limits derived from the stacking analyses are factors between 2 and 50 below the flux predicted by the simple leptonic scaling model, depending on the assumed source extent and the dominant target photon fields. Therefore, Terzan 5 still remains exceptional among all GCs, as the VHE \gamma-ray emission either arises from extra-ordinarily efficient leptonic processes, or from a recent catastrophic event, or is even unrelated to the GC itself.

Gamma-ray line signatures can be expected in the very-high-energy (VHE; E_\gamma > 100 GeV) domain due to self-annihilation or decay of dark matter (DM) particles in space. Such a signal would be readily distinguishable from astrophysical \gamma-ray sources that in most cases produce continuous spectra which span over several orders of magnitude in energy. Using data collected with the H.E.S.S. \gamma-ray instrument, upper limits on line-like emission are obtained in the energy range between ~500 GeV and ~25 TeV for the central part of the Milky Way halo and for extragalactic observations, complementing recent limits obtained with the Fermi-LAT instrument at lower energies. No statistically significant signal could be found. For monochromatic \gamma-ray line emission, flux limits of (2x10^-7 - 2x10^-5) m^-2 s^-1 sr^-1 and (1x10^-8 - 2x10^-6) m^-2 s^-1 sr^-1 are obtained for the central part of the Milky Way halo and extragalactic observations, respectively. For a DM particle mass of 1 TeV, limits on the velocity-averaged DM annihilation cross section < \sigma v >(\chi\chi -> \gamma\gamma) reach ~10^-27 cm^3 s^-1, based on the Einasto parametrization of the Galactic DM halo density profile.

The extragalactic background light (EBL) is the diffuse radiation with the second highest energy density in the Universe after the cosmic microwave background. The aim of this study is the measurement of the imprint of the EBL opacity to gamma-rays on the spectra of the brightest extragalactic sources detected with the High Energy Stereoscopic System (H.E.S.S.). The originality of the method lies in the joint fit of the EBL optical depth and of the intrinsic spectra of the sources, assuming intrinsic smoothness. Analysis of a total of ~10^5 gamma-ray events enables the detection of an EBL signature at the 8.8 std dev level and constitutes the first measurement of the EBL optical depth using very-high energy (E>100 GeV) gamma-rays. The EBL flux density is constrained over almost two decades of wavelengths (0.30-17 microns) and the peak value at 1.4 micron is derived as 15 +/- 2 (stat) +/- 3 (sys) nW / m^2 sr.

The previously unidentified very high-energy (VHE; E > 100 GeV) \gamma-ray source HESS J1303-631, discovered in 2004, is re-examined including new data from the H.E.S.S. Cherenkov telescope array. Archival data from the XMM-Newton X-ray satellite and from the PMN radio survey are also examined. Detailed morphological and spectral studies of VHE \gamma-ray emission as well as of the XMM-Newton X-ray data are performed. Significant energy-dependent morphology of the \gamma-ray source is detected with high-energy emission (E > 10 TeV) positionally coincident with the pulsar PSR J1301-6305 and lower energy emission (E <2 TeV) extending \sim 0.4^{\circ} to the South-East of the pulsar. The spectrum of the VHE source can be described with a power-law with an exponential cut-off N_{0} = (5.6 \pm 0.5) X 10^{-12} TeV^-1 cm^-2 s^-1, \Gamma = 1.5 \pm 0.2) and E_{\rm cut} = (7.7 \pm 2.2) TeV. The PWN is also detected in X-rays, extending \sim 2-3' from the pulsar position towards the center of the \gamma-ray emission region. The spectral energy distribution (SED) is well described by a one zone leptonic scenario which, with its associated caveats, predicts a very low average magnetic field for this source. Significant energy-dependent morphology of this source, as well as the identification of an associated X-ray PWN from XMM-Newton observations enable identification of the VHE source as an evolved PWN associated to the pulsar PSR J1303-6305. However, the large discrepancy in emission region sizes and the low level of synchrotron radiation suggest a multi-population leptonic nature. The low implied magnetic field suggests that the PWN has undergone significant expansion. This would explain the low level of synchrotron radiation and the difficulty in detecting counterparts at lower energies, the reason this source was originally classified as a "dark" VHE \gamma-ray source.

Vela X is a region of extended radio emission in the western part of the Vela constellation: one of the nearest pulsar wind nebulae (PWNe), and associated with the energetic Vela pulsar (PSR B0833-45). Extended very-high-energy (VHE) $\gamma$-ray emission (HESS $\mathrm{J0835\mhyphen 455}$) was discovered using the H.E.S.S. experiment in 2004. The VHE $\gamma$-ray emission was found to be coincident with a region of X-ray emission discovered with ${\it ROSAT}$ above 1.5 keV (the so-called \textit{Vela X cocoon}): a filamentary structure extending southwest from the pulsar to the centre of Vela X. A deeper observation of the entire Vela X nebula region, also including larger offsets from the cocoon, has been performed with H.E.S.S. This re-observation was carried out in order to probe the extent of the non-thermal emission from the Vela X region at TeV energies and to investigate its spectral properties. In order to increase the sensitivity to the faint $\gamma$-ray emission from the very extended Vela X region, a multivariate analysis method combining three complementary reconstruction techniques of Cherenkov-shower images is applied for the selection of $\gamma$-ray events. The analysis is performed with the On/Off background method, which estimates the background from separate observations pointing away from Vela X; towards regions free of $\gamma$-ray sources but with comparable observation conditions. The $\gamma$-ray surface brightness over the large Vela X region reveals that the detection of non-thermal VHE $\gamma$-ray emission from the PWN HESS $\mathrm{J0835\mhyphen 455}$ is statistically significant over a region of radius 1.2$^{\circ}$ around the position $\alpha$ = 08$^{\mathrm{h}}$ 35$^{\mathrm{m}}$ 00$^{\mathrm{s}}$, $\delta$ = -45$^{\circ}$ 36$^{\mathrm{\prime}}$ 00$^{\mathrm{\prime}\mathrm{\prime}}$ (J2000).

The energy spectrum of magnetohydrodynamic turbulence attracts interest due to its fundamental importance and its relevance for interpreting astrophysical data. Here we present measurements of the energy spectra from a series of high-resolution direct numerical simulations of MHD turbulence with a strong guide field and for increasing Reynolds number. The presented simulations, with numerical resolutions up to 2048^3 mesh points and statistics accumulated over 30 to 150 eddy turnover times, constitute, to the best of our knowledge, the largest statistical sample of steady state MHD turbulence to date. We study both the balanced case, where the energies associated with Alfv\'en modes propagating in opposite directions along the guide field, E^+ and $E^-, are equal, and the imbalanced case where the energies are different. In the balanced case, we find that the energy spectrum converges to a power law with exponent -3/2 as the Reynolds number is increased, consistent with phenomenological models that include scale-dependent dynamic alignment. For the imbalanced case, with E^+>E^-, the simulations show that E^- ~ k_{\perp}^{-3/2} for all Reynolds numbers considered, while E^+ has a slightly steeper spectrum at small Re. As the Reynolds number increases, E^+ flattens. Since both E^+ and E^- are pinned at the dissipation scale and anchored at the driving scales, we postulate that at sufficiently high Re the spectra will become parallel in the inertial range and scale as E^+ ~ E^- ~ k_{\perp}^{-3/2}. Questions regarding the universality of the spectrum and the value of the "Kolmogorov constant" are discussed.

We present the significant detection of the first extragalactic pulsar wind nebula (PWN) detected in gamma rays, N157B, located in the large Magellanic Cloud (LMC). Pulsars with high spin-down luminosity are found to power energised nebulae that emit gamma rays up to energies of several tens of TeV. N157B is associated with PSRJ0537-6910, which is the pulsar with the highest known spin-down luminosity. The High Energy Stereoscopic System telescope array observed this nebula on a yearly basis from 2004 to 2009 with a dead-time corrected exposure of 46 h. The gamma-ray spectrum between 600 GeV and 12 TeV is well-described by a pure power-law with a photon index of 2.8 \pm 0.2(stat) \pm 0.3(syst) and a normalisation at 1 TeV of (8.2 \pm 0.8(stat) \pm 2.5(syst)) \times 10^-13 cm^-2s^-1TeV^-1. A leptonic multi-wavelength model shows that an energy of about 4 \times 10^49erg is stored in electrons and positrons. The apparent efficiency, which is the ratio of the TeV gamma-ray luminosity to the pulsar's spindown luminosity, 0.08% \pm 0.01%, is comparable to those of PWNe found in the Milky Way. The detection of a PWN at such a large distance is possible due to the pulsar's favourable spin-down luminosity and a bright infrared photon-field serving as an inverse-Compton-scattering target for accelerated leptons. By applying a calorimetric technique to these observations, the pulsar's birth period is estimated to be shorter than 10 ms.

Very-high-energy (VHE; E >100 GeV) and high-energy (HE; 100 MeV < E < 100 GeV) data from \gamma-ray observations performed with the H.E.S.S. telescope array and the Fermi-LAT instrument, respectively, are analysed in order to investigate the non-thermal processes in the starburst galaxy NGC 253. The VHE \gamma-ray data can be described by a power law in energy with differential photon index \Gamma=2.14 \pm 0.18_stat \pm 0.30_sys and differential flux normalisation at 1 TeV of F_0 = (9.6 \pm 1.5_stat (+5.7,-2.9)_sys) x 10^{-14} TeV^{-1} cm^{-2} s^{-1}. A power-law fit to the differential HE \gamma-ray spectrum reveals a photon index of \Gamma=2.24 \pm 0.14_stat \pm 0.03_sys and an integral flux between 200 MeV and 200 GeV of F(0.2-200 GeV) = (4.9 \pm 1.0_stat \pm 0.3_sys) x 10^{-9} cm^{-2} s^{-1}. No evidence for a spectral break or turnover is found over the dynamic range of both the LAT instrument and the H.E.S.S. experiment: a combined fit of a power law to the HE and VHE \gamma-ray data results in a differential photon index \Gamma=2.34 \pm 0.03 with a p-value of 30%. The \gamma-ray observations indicate that at least about 20% of the energy of the cosmic rays capable of producing hadronic interactions is channeled into pion production. The smooth alignment between the spectra in the HE and VHE \gamma-ray domain suggests that the same transport processes dominate in the entire energy range. Advection is most likely responsible for charged particle removal from the starburst nucleus from GeV to multiple TeV energies. In a hadronic scenario for the \gamma-ray production, the single overall power-law spectrum observed would therefore correspond to the mean energy spectrum produced by the ensemble of cosmic-ray sources in the starburst region.