Search

We perform a large-scale hydrodynamic simulation of a massive star cluster whose stellar population mimics that of the Cygnus OB2 association. The main-sequence stars are first simulated during 1.6 Myr, until a quasi-stationary state is reached. At this time the three Wolf-Rayet stars observed in Cygnus OB2 are added to the simulation, which continues to 2 Myr. Using a high-resolution grid in the centre of the domain, we can resolve the most massive stars individually, which allows us to probe the kinetic structures at small (parsec) scales. We find that, although the cluster excavates a spherical "superbubble" cavity, the stellar population is too loosely distributed to blow a large-scale cluster wind termination shock, and that collective effects from wind-wind interactions are much less efficient than usually assumed. This challenges our understanding of the ultra-high energy emission observed from the region., Comment: 15 pages, 12 figures

We consider scenarios for non-thermal particle acceleration and re-acceleration in the central cores of compact massive star clusters, aided by insights from high resolution hydrodynamic simulations. We show that i) particles are unlikely to interact with many shocks during their lifetimes in the core; ii) colliding flows do not produce hard spectra; iii) turbulent re-acceleration in the core is suppressed. Inefficient re-acceleration mechanisms are not expected to produce hard components nor to increase the maximum energy within the cores of massive star clusters. Models in which the observed ultra-high energy gamma rays originate in the core of massive stellar clusters are thus disfavoured., Comment: 16 pages, 14 figures, accepted in MNRAS

Despite of various similar features, Functional Data Analysis and High-Dimensional Data Analysis are two major fields in Statistics that grew up recently almost independently one from each other. The aim of this paper is to propose a survey on methodological advances for variable selection in functional regression, which is typically a question for which both functional and multivariate ideas are crossing. More than a simple survey, this paper aims to promote even more new links between both areas., Comment: 22 pages

A new sparse semiparametric model is proposed, which incorporates the influence of two functional random variables in a scalar response in a flexible and interpretable manner. One of the functional covariates is included through a single-index structure, while the other is included linearly through the high-dimensional vector formed by its discretised observations. For this model, two new algorithms are presented for selecting relevant variables in the linear part and estimating the model. Both procedures utilise the functional origin of linear covariates. Finite sample experiments demonstrated the scope of application of both algorithms: the first method is a fast algorithm that provides a solution (without loss in predictive ability) for the significant computational time required by standard variable selection methods for estimating this model, and the second algorithm completes the set of relevant linear covariates provided by the first, thus improving its predictive efficiency. Some asymptotic results theoretically support both procedures. A real data application demonstrated the applicability of the presented methodology from a predictive perspective in terms of the interpretability of outputs and low computational cost., Comment: 33 pages, 6 figures, 10 tables

A fast and flexible $k$NN procedure is developed for dealing with a semiparametric functional regression model involving both partial-linear and single-index components. Rates of uniform consistency are presented. Simulated experiments highlight the advantages of the $k$NN procedure. A real data analysis is also shown., Comment: 14 pages, 1 figure, 6 tables

This paper aims to front with dimensionality reduction in regression setting when the predictors are a mixture of functional variable and high-dimensional vector. A flexible model, combining both sparse linear ideas together with semiparametrics, is proposed. A wide scope of asymptotic results is provided: this covers as well rates of convergence of the estimators as asymptotic behaviour of the variable selection procedure. Practical issues are analysed through finite sample simulated experiments while an application to Tecator's data illustrates the usefulness of our methodology., Comment: 40 pages, 7 figures, 5 tables

This paper develops a new automatic and location-adaptive procedure for estimating regression in a Functional Single-Index Model (FSIM). This procedure is based on $k$-Nearest Neighbours ($k$NN) ideas. The asymptotic study includes results for automatically data-driven selected number of neighbours, making the procedure directly usable in practice. The local feature of the $k$NN approach insures higher predictive power compared with usual kernel estimates, as illustrated in some finite sample analysis. As by-product we state as preliminary tools some new uniform asymptotic results for kernel estimates in the FSIM model., Comment: 31 pages, 5 figures, 4 tables

Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black holes and neutron star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with such pair plasmas. Their behaviour is quite different from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components and their role in the dynamics of such compact objects is believed to be fundamental. So far, our experimental inability to produce large yields of positrons in quasi-neutral beams has restricted the understanding of electron-positron pair plasmas to simple numerical and analytical studies which are rather limited. We present first experimental results confirming the generation of high-density, quasi-neutral, relativistic electron-positron pair beams using the 440 GeV/c beam at CERN's Super Proton Synchrotron (SPS) accelerator. The produced pair beams have a volume that fills multiple Debye spheres and are thus able to sustain collective plasma oscillations. Our work opens up the possibility of directly probing the microphysics of pair plasmas beyond quasi-linear evolution into regimes that are challenging to simulate or measure via astronomical observations., Comment: 14 pages, 8 figures

DOLCE, the first top-level (foundational) ontology to be axiomatized, has remained stable for twenty years and today is broadly used in a variety of domains. DOLCE is inspired by cognitive and linguistic considerations and aims to model a commonsense view of reality, like the one human beings exploit in everyday life in areas as diverse as socio-technical systems, manufacturing, financial transactions and cultural heritage. DOLCE clearly lists the ontological choices it is based upon, relies on philosophical principles, is richly formalized, and is built according to well-established ontological methodologies, e.g. OntoClean. Because of these features, it has inspired most of the existing top-level ontologies and has been used to develop or improve standards and public domain resources (e.g. CIDOC CRM, DBpedia and WordNet). Being a foundational ontology, DOLCE is not directly concerned with domain knowledge. Its purpose is to provide the general categories and relations needed to give a coherent view of reality, to integrate domain knowledge, and to mediate across domains. In these 20 years DOLCE has shown that applied ontologies can be stable and that interoperability across reference and domain ontologies is a reality. This paper briefly introduces the ontology and shows how to use it on a few modeling cases., Comment: 25 pages, 7 figures

Context: Young massive star clusters (YMCs) have come increasingly into the focus of discussions on the origin of galactic cosmic rays (CRs). The proposition of CR acceleration inside superbubbles (SBs) blown by the strong winds of these clusters avoids issues faced by the standard paradigm of acceleration at supernova remnant shocks. Aims: We provide an interpretation of the latest TeV $\gamma$-ray observations of the region around the YMC Westerlund 1 taken with the High Energy Stereoscopic System (H.E.S.S.) in terms of diffusive shock acceleration at the cluster wind termination shock, taking into account the spectrum and morphology of the emission. As Westerlund 1 is a prototypical example of a YMC, such a study is relevant to the general question about the role of YMCs for the Galactic CR population. Methods: We generate model $\gamma$-ray spectra, characterise particle propagation inside the SB based on the advection, diffusion, and cooling timescales, and constrain key parameters of the system. We consider hadronic emission from proton-proton interaction and subsequent pion decay and leptonic emission from inverse Compton scattering on all relevant photon fields, including the CMB, diffuse and dust-scattered starlight, and the photon field of Westerlund 1 itself. The effect of the magnetic field on cooling and propagation is discussed. Klein-Nishina effects are found to be important in determining the spectral evolution of the electron population. Results: A leptonic origin of the bulk of the observed $\gamma$-rays is preferable. The model is energetically plausible, consistent with the presence of a strong shock, and allows for the observed energy-independent morphology. The hadronic model faces two main issues: confinement of particles to the emission region and an unrealistic energy requirement., Comment: Accepted for publication in Astronomy & Astrophyics

We demonstrate that supernova remnant (SNR) shocks embedded within massive star clusters can reproduce both the cosmic-ray proton and all-particle spectra measured in the vicinity of the Earth up to hundreds of peta-electronvolts (PeV). We model two classes of massive star clusters. The first population are "loose clusters" which do not power a collective wind termination shock. SNR shocks then expand in a low-density and weakly magnetised medium, and this population mainly contributes up to the "knee" of the CR spectrum around 1 PeV. The second population are young compact clusters, which are powerful and compact enough to sustain a collective wind outflow. SNR shocks then expand from the cluster into the strongly magnetised wind and accelerate nuclei up to ultra-high energies. This population, representing only about 15% of all Galactic massive star clusters, nevertheless dominates the spectrum between ~ 1 and 100 PeV. While these two components alone can reproduce the shape of the CR spectrum up to hundreds of PeV, adding a light sub-ankle extragalactic component motivated by composition and anisotropy measurements, allows to reproduce the spectrum up to the highest energies. Fitting parameters are systematically linked to physical variables whose values are in line with theoretical expectations., Comment: 12 pages, 7 figures, accepted in MNRAS

We critically assess limits on the maximum energy of protons accelerated within superbubbles around massive stellar clusters, considering a number of different scenarios. In particular, we derive under which circumstances acceleration of protons above peta-electronvolt (PeV) energies can be expected. While the external forward shock of the superbubble may account for acceleration of particles up to 100 TeV, internal primary shocks such as supernova remnants expanding in the low density medium or the collective wind termination shock which forms around a young compact cluster provide more favourable channels to accelerate protons up to 1 PeV, and possibly beyond. Under reasonable conditions, clustered supernovae launching powerful shocks into the magnetised wind of a young and compact massive star cluster are found to be the most promising systems to accelerate protons above 10 PeV. On the other hand, stochastic re-acceleration in the strongly turbulent plasma is found to be much less effective than claimed in previous works, with a maximum proton energy of at most a few hundred TeV., Comment: 10 pages, 5 figures

We compute the production of cosmic rays in the dynamical superbubble produced by a cluster of massive stars. Stellar winds, supernova remnants and turbulence are found to accelerate particles so efficiently that the nonlinear feedback of the particles must be taken into account in order to ensure that the energy balance is not violated. High energy particles do not scatter efficiently on the turbulence and escape quickly after each supernova explosion, which makes both their intensity inside the bubble and injection in the interstellar medium intermittent. On the other hand, the stochastic acceleration of low energy particles hardens the spectra at GeV energies. Because cosmic rays damp the turbulence cascade, this hardening is less pronounced when nonlinearities are taken into account. Nevertheless, spectra with hard components extending up to 1 to 10 GeV and normalised to an energy density of 1 to 100 eV cm$^{-3}$ are found to be typical signatures of cosmic rays produced in superbubbles. Efficient shock reacceleration within compact clusters is further shown to produce hard, slightly concave spectra, while the presence of a magnetised shell is shown to enhance the confinement of cosmic rays in the bubble and therefore the collective plasma effects acting on them. We eventually estimate the overall contribution of superbubbles to the galactic cosmic ray content and show typical gamma-ray spectra expected from hadronic interactions in superbuble shells. In both cases, a qualitative agreement with observations is obtained., Comment: 20 pages, 14 figures

When the pressure of particles accelerated at shock waves is no longer negligible compared to the kinetic pressure of the gas, the linear theory of diffusive shock acceleration breaks down. This is expected in particular when the shock sweeps up preexisting cosmic rays, or when multiple shocks reaccelerate successively the same particles. To describe these systems, one has to account for the nonlinear backreaction of the particles on the magnetohydrodynamic flow. Using an up-to-date semi-analytical model of particle reacceleration at nonlinear shocks, we show that the presence of prexisting energetic particles strongly affects the shock profile, in such a way that the reacceleration of non thermal particles or the acceleration of particles from the thermal bath becomes less efficient. We further describe the evolution of the distribution of particles after several shocks and study the properties of the asymptotic solution. We detail the case of identical shocks as well as more realistic scenarios, including the heating of the medium or superbubble environments. When the particles are efficiently confined in the acceleration region, it is generally found that the spectrum converges toward a concave solution after a few tens of shocks, with a spectral index around 3.5 at the highest energy. The postshock cosmic ray pressure reaches an asymptotic value of about 4-5% of the ram pressure of one shock. Most of the shock pressure is transferred to escaping particles., Comment: 9 pages, 8 figures

We model the diffusive shock acceleration of particles in a system of two colliding shock waves and present a method to solve the time-dependent problem analytically in the test-particle approximation and high energy limit. In particular, we show that in this limit the problem can be analysed with the help of a self-similar solution. While a number of recent works predict hard ($E^{-1}$) spectra for the accelerated particles in the stationary limit, or the appearance of spectral breaks, we found instead that the spectrum of accelerated particles in a time-dependent collision follows quite closely the canonical $E^{-2}$ prediction of diffusive shock acceleration at a single shock, except at the highest energy, where a hardening appears, originating a bumpy feature just before the exponential cutoff. We also investigated the effect of the reacceleration of pre-existing cosmic rays by a system of two shocks, and found that under certain conditions spectral features can appear in the cutoff region. Finally, the mathematical methods presented here are very general and could be easily applied to a variety of astrophysical situations, including for instance standing shocks in accretion flows, diverging shocks, backward collisions of a slow shock by a faster shock, and wind-wind or shock-wind collisions., Comment: 11 pages, 8 figures

The first-order electroweak phase transition in the early universe could occur in multiple steps leading to specific multi-peaked signatures in the primordial gravitational wave (GW) spectrum. We argue that these signatures are generic phenomena in multi-scalar extensions of the Standard Model. In a simple example of such an extension, we have studied the emergence of reoccurring and nested vacuum bubble configurations and their role in the formation of multiple peaks in the GW spectrum. The conditions for potential detectability of these features by the forthcoming generation of interferometers have been studied., Comment: 7 pages, 3 captioned figures, 1 table; numerical analysis extended, conclusions unchanged; contributed talk by R. Pasechnik at the EPS-HEP2019, 10-17 July, 2019, Ghent, Belgium

The rich vacuum structure of multi-Higgs extensions of the Standard Model (SM) may have interesting cosmological implications for the electroweak phase transition (EWPT). As an important example of such class of models, we consider a particularly simple low-energy SM-like limit of a recently proposed Grand-Unified Trinification model with the scalar sector composed of two Higgs doublets and a complex singlet and with a global $\mathrm{U}(1)$ family symmetry. The fermion sector of this model is extended with a family of vector-like quarks which enhances CP violation. With the current study, we aim at exploring the generic vacuum structure and uncovering the features of the EWPT in this model relevant for cosmology. We show the existence of different phase transition patterns providing strong departure from thermal equilibrium. Most of these observations are not specific to the considered model and may generically be expected in other multi-Higgs extensions of the SM., Comment: 48 pages, 13 figures

We study static distributions of ferrofluid submitted to non-uniform magnetic fields. We show how the normal-field instability is modified in the presence of a weak magnetic field gradient. Then we consider a ferrofluid droplet and show how the gradient affects its shape. A rich phase transitions phenomenology is found. We also investigate the creation of droplets by successive splits when a magnet is vertically approached from below and derive theoretical expressions which are solved numerically to obtain the number of droplets and their aspect ratio as function of the field configuration. A quantitative comparison is performed with previous experimental results, as well as with our own experiments, and yields good agreement with the theoretical modeling., Comment: 41 pages, 18 figures