Source author record

A. Costa

A. Costa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2025arXiv

How do recollimation-induced instabilities shape the propagation of hydrodynamic relativistic jets?

Recollimation is a phenomenon of particular importance in the dynamic evolution of jets and in the emission of high-energy radiation. Additionally, the full comprehension of this phenomenon provides insights into fundamental properties of jets in the vicinity of the Active Galactic Nucleus (AGN). Three-dimensional (magneto-)hydrodynamic simulations revealed that the jet conditions at recollimation favor the growth of strong instabilities, challenging the traditional view-supported from two-dimensional simulations-of confined jets undergoing a series of recollimation and reflection shocks. To investigate the stability of relativistic jets in AGNs at recollimation sites, we perform a set of long duration three-dimensional relativistic hydrodynamic simulations with the state-of-the-art PLUTO code, to focus on the development of hydrodynamical instabilities. We explore the non-linear growth of the instabilities and their effects on the physical jet properties as a function of the initial jet parameters: jet Lorentz factor, temperature, opening angle and jet-environment density-contrast. The parameter space is designed to describe low-power, weakly magnetized jets at small distances from the core (around the parsec scale). All collimating jets we simulated develop instabilities. Recollimation instabilities decelerate the jet, heat it, entrain external material, and move the recollimation point to shorter distances from the core. This is true for both conical and cylindrical jets. The instabilities, that are first triggered by the centrifugal instability, appear to be less disruptive in the case of narrower, denser, more relativistic, and warmer jets. These results provide valuable insights into the complex processes governing AGN jets and could be used to model the properties of low-power, weakly magnetized jetted AGNs.

preprint2022arXiv

A configurable computer simulation model for reducing patient waiting time in oncology departments

Nowadays, the increase in patient demand and the decline in resources are lengthening patient waiting times in many chemotherapy oncology departments. Therefore, enhancing healthcare services is necessary to reduce patient complaints. Reducing the patient waiting times in the oncology departments represents one of the main goals of healthcare manager. Simulation models are considered an effective tool for identifying potential ways to improve patient flow in oncology departments. This paper presents a new agent-based simulation model designed to be configurable and adaptable to the needs of oncology departments which have to interact with an external pharmacy. When external pharmacies are utilized, a courier service is needed to deliver the individual therapies from the pharmacy to the oncology department. An oncology department located in southern Italy was studied through the simulation model and different scenarios were compared with the aim of selecting the department configuration capable of reducing the patient waiting times.

preprint2022arXiv

ASTRI Mini-Array Core Science at the Observatorio del Teide

The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Project led by the Italian National Institute for Astrophysics (INAF) is developing and will deploy at the Observatorio del Teide a mini-array (ASTRI Mini-Array) composed of nine telescopes similar to the small-size dual-mirror Schwarzschild-Couder telescope (ASTRI-Horn) currently operating on the slopes of Mt. Etna in Sicily. The ASTRI Mini-Array will surpass the current Cherenkov telescope array differential sensitivity above a few tera-electronvolt (TeV), extending the energy band well above hundreds of TeV. This will allow us to explore a new window of the electromagnetic spectrum, by convolving the sensitivity performance with excellent angular and energy resolution figures. In this paper we describe the Core Science that we will address during the first four years of operation, providing examples of the breakthrough results that we will obtain when dealing with current open questions, such as the acceleration of cosmic rays, cosmology and fundamental physics and the new window, for the TeV energy band, of the time-domain astrophysics.

preprint2022arXiv

Extragalactic Observatory Science with the ASTRI Mini-Array at the Observatorio del Teide

The ASTRI Mini-Array is a next-generation system of nine imaging atmospheric Cherenkov telescopes that is going to be built at the Observatorio del Teide site. After a first phase, in which the instrument will be operated as an experiment prioritizing a schedule of primary science cases, an observatory phase is foreseen in which other significant targets will be pointed. We focus on the observational feasibility of extragalactic sources and on astrophysical processes that best complement and expand the ASTRI Mini-Array core science, presenting the most relevant examples that are at reach of detection over long-term time scales and whose observation can provide breakthrough achievements in the very-high energy extragalactic science. Such examples cover a wide range of $γ$-ray emitters, including the study of AGN low states in the multi-TeV energy range, the possible detection of Seyfert galaxies with long exposures and the searches of dark matter lines above 10 TeV. Simulations of the presented objects show that the instrument performance will be competitive at multi-TeV energies with respect to current arrays of Cherenkov telescopes.

preprint2022arXiv

Extreme blazars: the result of unstable recollimated jets?

Extreme BL Lacs (EHBL) form a subclass of blazars which challenge standard emission scenarios. In a recent study it has been argued that their peculiar properties can be explained if emitting electrons are accelerated in a series of oblique shocks induced by the recollimation of the relativistic jet. However, new 3D simulations of recollimated, weakly magnetized jets reveal that, in correspondence with the first recollimation shock, the flow develops a rapidly growing instability, becomes highly turbulent and decelerates, effectively hampering the formation of the multiple shock structure routinely observed in 2D simulations. Building on these new findings, we propose here a revised scenario for EHBL, in which the emission is produced by electrons accelerated at the recollimation shock and subsequently further energized through stochastic acceleration in the turbulent downstream flow. We apply a simple version of this scenario to the prototypical EHBL 1ES 0229-200, showing that the SED can be satisfactorily reproduced with standard values of the main physical parameters.

preprint2022arXiv

Galactic Observatory Science with the ASTRI Mini-Array at the Observatorio del Teide

The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Mini-Array will be composed of nine imaging atmospheric Cherenkov telescopes at the Observatorio del Teide site. The array will be best suited for astrophysical observations in the 0.3-200 TeV range with an angular resolution of few arc-minutes and an energy resolution of 10-15\%. A core-science programme in the first four years will be devoted to a limited number of key targets, addressing the most important open scientific questions in the very-high energy domain. At the same time, thanks to a wide field of view of about 10 degrees, ASTRI Mini-Array will observe many additional field sources, which will constitute the basis for the long-term observatory programme that will eventually cover all the accessible sky. In this paper, we review different astrophysical Galactic environments, e.g. pulsar wind nebulae, supernova remnants, and gamma-ray binaries, and show the results from a set of ASTRI Mini-Array simulations of some of these field sources made to highlight the expected performance of the array (even at large offset angles) and the important additional observatory science that will complement the core-science program.

preprint2022arXiv

Pseudostreamer influence on flux rope evolution

A critical aspect of solar activity is the coupling between eruptions and the surrounding coronal magnetic field, which determines the trajectory and morphology of the eruptive event. Pseudostreamers (PSs) are coronal magnetic structures formed by arcs of twin loops capped by magnetic field lines from coronal holes of the same polarity that meet at a central spine. They contain a single magnetic null point in the spine, potentially influencing the evolution of nearby flux ropes (FRs). To understand the net effect of the PS on FR eruptions is first necessary to study diverse and isolated FR-PS scenarios, which are not influenced by other magnetic structures. We performed numerical simulations in which a FR structure is in the vicinity of a PS magnetic configuration. The combined magnetic field of the PS and the FR results in the formation of two magnetic null points. We evolve this scenario by numerically solving the magnetohydrodynamic equations in 2.5D. The simulations consider a fully ionised compressible ideal plasma in the presence of a gravitational field and a stratified atmosphere. We find that the dynamic behaviour of the FR can be categorised into three different classes based on the FR trajectories and whether it is eruptive or confined. Our analysis indicates that the magnetic null points are decisive in the direction and intensity of the FR deflection and their hierarchy depends on the topological arrangement of the scenario. Moreover, the PS lobe acts as a magnetic cage enclosing the FR. We report that the total unsigned magnetic flux of the cage is a key parameter defining whether the FR is ejected or not.

preprint2021arXiv

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology.

preprint2021arXiv

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]"

preprint2020arXiv

First detection of the Crab Nebula at TeV energies with a Cherenkov telescope in a dual-mirror Schwarzschild-Couder configuration: the ASTRI-Horn telescope

We report on the first detection of very high-energy (VHE) gamma-ray emission from the Crab Nebula by a Cherenkov telescope in dual-mirror Schwarzschild-Couder (SC) configuration. The result has been achieved by means of the 4 m size ASTRI-Horn telescope, operated on Mt. Etna (Italy) and developed in the context of the Cherenkov Telescope Array Observatory preparatory phase. The dual-mirror SC design is aplanatic and characterized by a small plate scale, allowing us to implement large field of view cameras with small-size pixel sensors and a high compactness. The curved focal plane of the ASTRI camera is covered by silicon photo-multipliers (SiPMs), managed by an unconventional front-end electronics based on a customized peak-sensing detector mode. The system includes internal and external calibration systems, hardware and software for control and acquisition, and the complete data archiving and processing chain. The observations of the Crab Nebula were carried out in December 2018, during the telescope verification phase, for a total observation time (after data selection) of 24.4 h, equally divided into on- and off-axis source exposure. The camera system was still under commissioning and its functionality was not yet completely exploited. Furthermore, due to recent eruptions of the Etna Volcano, the mirror reflection efficiency was reduced. Nevertheless, the observations led to the detection of the source with a statistical significance of 5.4 sigma above an energy threshold of ~3 TeV. This result provides an important step towards the use of dual-mirror systems in Cherenkov gamma-ray astronomy. A pathfinder mini-array based on nine large field-of-view ASTRI-like telescopes is under implementation.

preprint2020arXiv

Scientific prospects for a mini-array of ASTRI telescopes: a gamma-ray TeV data challenge

ASTRI is a project aiming at the realization of a gamma-ray imaging Cherenkov telescope that observes the sky in the TeV band. Recently, the development of a mini-array (MA) of ASTRI telescopes has been funded by the Istituto Nazionale di Astrofisica. The ASTRI Comprehensive Data Challenge (ACDC) project aims at optimizing the scientific exploitation and analysis techniques of the ASTRI MA, by performing a complete end-to-end simulation of a tentative scientific program, from the generation of suitable instrument response functions to the proposal, selection, analysis, and interpretation of the simulated data. We assumed that the MA will comprise nine ASTRI telescopes arranged in a (almost) square geometry (mean distance between telescopes of ~250m). We simulated three years of observations, adopting a realistic pointing plan that takes into account, for each field, visibility constraints for an assumed site in Paranal (Chile) and observational time slots in dark sky conditions. We simulated the observations of nineteen Galactic and extragalactic fields selected for their scientific interest, including several classes of objects (such as pulsar wind nebulae, supernova remnants, gamma-ray binaries etc), for a total of 81 point-like and extended sources. Here we present an overview of the ACDC project, providing details on the different software packages needed to carry out the simulated three-years operation of the ASTRI MA. We discuss the results of a systematic analysis applied on the whole simulated data, by making use of prototype science tools widely adopted by the TeV astronomical community. Furthermore, particular emphasis is also given to some targets used as benchmarks.

preprint2020arXiv

Slow and sausage loop mode excitation due to local and global spontaneous perturbations

We analyse the capability of different type of perturbations -associated with usual environment energy fluctuations of the solar corona- to excite slow and sausage modes in solar flaring loops. We perform numerical simulations of the MHD ideal equations considering straight plasma magnetic tubes subject to local and global energy depositions. We find that local loop energy depositions of typical microflares [$\sim$($10^{27}$-$10^{30}$) erg] are prone to drive slow shock waves that induce slow mode patterns. The slow mode features are obtained for every tested local energy deposition inside the loop. Meanwhile, to obtain an observable sausage mode pattern a global perturbation, capable to modify instantaneously the internal loop temperature, is required, i.e. the characteristic conductive heating time must be much smaller than the radiative cooling one. Experiments carried out by varying parameter $β$ show us that the excitation of sausage modes does not depend significantly on the value of this parameter but on the global or local character of the energy source.

preprint2020arXiv

Solar slow magneto-acoustic-gravity waves: an erratum correction and a revisited scenario

Slow waves are commonly observed on the entire solar atmosphere. Assuming a thin flux tube approximation, the cut-off periods of slow-mode magneto-acoustic-gravity waves that travel from the photosphere to the corona were obtained in Costa et al. (2018). In that paper, however, a typo in the specific heat coefficient at constant pressure $c_{\mathrm{p}}$ value led to an inconsistency in the cut-off calculation, which is only significant at the transition region. Due to the abrupt temperature change in the region, a change of the mean atomic weight (by a factor of approximately two) also occurs, but is often overlooked in analytical models for simplicity purposes. In this paper, we revisit the calculation of the cut-off periods of magneto-acoustic-gravity waves in Costa et al. (2018) by considering an atmosphere in hydrostatic equilibrium with a temperature profile, with the inclusion of the variation of the mean atomic weight and the correction of the inconsistency aforementioned. In addition, we show that the cut-off periods obtained analytically are consistent with the corresponding periods measured in observations of a particular active region.

preprint2018arXiv

Slow magnetoacoustic gravity waves in an equilibrium stratified solar atmosphere: cut-off periods through the transition region

Assuming the thin flux tube approximation, we introduce an analytical model that contemplates the presence of: a non-isothermal temperature; a varying magnetic field and a non-uniform stratified medium in hydrostatic equilibrium due to a constant gravity acceleration. This allows the study of slow magnetoacoustic cut-off periods across the solar transition region, from the base of the solar chromosphere to the lower corona. The used temperature profile approaches the VAC solar atmospheric model. The periods obtained are consistent with observations. Similar to the acoustic cut-off periods, the resulting magnetoacoustic gravity ones follow the sharp temperature profile, but shifted towards larger heights; in other words, at a given height the magnetoacoustic cut-off period is significantly lower than the corresponding acoustic one. Along a given longitude of an inclined thin magnetic tube, the greater its inclination the softener the temperature gradient it crosses. Changes in the magnetic field intensity do not significantly modify the periods at the coronal level but modulate the values below the transition region within periods between $\sim [2\,- 6]\,$min. Within the limitations of our model, we show that monochromatic oscillations of the solar atmosphere are the atmospheric response at its natural frequency to random or impulsive perturbations, and not a consequence of the forcing from the photosphere.

preprint2017arXiv

MHD simulations of coronal supra-arcade downflows including anisotropic thermal conduction

Coronal supra-arcade downflows (SADs) are observed as dark trails descending towards hot turbulent fan shaped regions. Due to the large temperature values, and gradients in these fan regions the thermal conduction should be very efficient. While several models have been proposed to explain the triggering and the evolution of SADs, none of these scenarios address a systematic consideration of thermal conduction. Thus, we accomplish this task numerically simulating the evolution of SADs within this framework. That is, SADs are conceived as voided (subdense) cavities formed by non-linear waves triggered by downflowing bursty localized reconnection events in a perturbed hot fan. We generate a properly turbulent fan, obtained by a stirring force that permits control of the energy and vorticity input in the medium where SADs develop. We include anisotropic thermal conduction and consider plasma properties consistent with observations. Our aim is to study if it is possible to prevent SADs to vanish by thermal diffusion. We find that this will be the case, depending on the turbulence parameters. In particular, if the magnetic field lines are able to envelope the voided cavities, thermally isolating them from the hot environment. Velocity shear perturbations that are able to generate instabilities of the Kelvin-Helmholtz type help to produce magnetic islands, extending the life-time of SADs.

preprint2017arXiv

Prospects for CTA observations of the young SNR RX J1713.7-3946

We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX~J1713.7$-$3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very-high-energy (VHE) gamma rays. Special attention is paid to explore possible spatial (anti-)correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H{\sc i} emission. We present a series of simulated images of RX J1713.7$-$3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the non-thermal X-ray emission observed by {\it XMM-Newton}, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H{\sc i} observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic vs leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator.

preprint2016arXiv

Photo-ionization of planetary winds: case study HD209458b

Close-in hot Jupiters are exposed to a tremendous photon flux that ionizes the neutral escaping material from the planet leaving an observable imprint that makes them an interesting laboratory for testing theoretical models. In this work we present 3D hydrodynamic simulations with radiation transfer calculations of a close-in exoplanet in a blow-off state. We calculate the Ly-$α$ absorption and compare it with observations of HD 209458b an previous simplified model results.Our results show that the hydrodynamic interaction together with a proper calculation of the photoionization proccess are able to reproduce the main features of the observed Ly-$α$ absorption, in particular at the blue-shifted wings of the line. We found that the ionizing stellar flux produce an almost linear effect on the amount of absorption in the wake. Varying the planetary mass loss rate and the radiation flux, we were able to reproduce the $10\%$ absorption observed at $-100~\mathrm{km~s^{-1}}$.

preprint2015arXiv

3D MHD simulation of flare supra-arcade downflows in a turbulent current sheet medium

Supra-arcade downflows (SADs) are sunward, generally dark, plasma density depletions originated above posteruption flare arcades. In this paper using 3D MHD simulations we investigate if the SAD cavities can be produced by a direct combination of the tearing mode and Kelvin-Helmholtz instabilities leading to a turbulent current sheet (CS) medium or if the current sheet is merely the background where SADs are produced triggered by an impulsive deposition of energy. We find that to give account of the observational dark lane structures an addition of local energy, provided by a reconnection event, is required. We suggest that there may be a closed relation between characteristic SAD sizes and CS widths that must be satisfied to obtain an observable SAD.

preprint2015arXiv

A geometric interpretation of the Schützenberger group of a minimal subshift

The first author has associated in a natural way a profinite group to each irreducible subshift. The group in question was initially obtained as a maximal subgroup of a free profinite semigroup. In the case of minimal subshifts, the same group is shown in the present paper to also arise from geometric considerations involving the Rauzy graphs of the subshift. Indeed, the group is shown to be isomorphic to the inverse limit of the profinite completions of the fundamental groups of the Rauzy graphs of the subshift. A further result involving geometric arguments on Rauzy graphs is a criterion for freeness of the profinite group of a minimal subshift based on the Return Theorem of Berthé et. al.

preprint2015arXiv

Homological Domination in Large Random Simplicial Complexes

In this paper we state the homological domination principle for random multi-parameter simplicial complexes, claiming that the Betti number in one specific dimension (which is explicitly determined by the probability multi-parameter) significantly dominates the Betti numbers in all other dimensions. We also state and discuss evidence for two interesting conjectures which would imply a stronger version of the homological domination principle, namely that generically homology of a random simplicial complex coincides with that of a wedges of k-dimensional spheres. These two conjectures imply that under an additional assumption (specified in the paper) a random simplicial complex collapses to a k-dimensional complex homotopy equivalent to a wedge of spheres of dimension k.

preprint2015arXiv

Large random simplicial complexes, I

In this paper we develop further the multi-parameter model of random simplicial complexes. Firstly, we give an intrinsic characterisation of the multi-parameter probability measure. Secondly, we show that in multi-parameter random simplicial complexes the links of simplexes and their intersections are also multi-parameter random simplicial complexes. Thirdly, we find conditions under which a multi-parameter random simplicial complex is connected and simply connected.

preprint2015arXiv

Large Random Simplicial Complexes, II; the fundamental groups

In our recent work we described conditions under which a multi-parameter random simplicial complex is connected and simply connected. We showed that the Betti numbers of multi-parameter random simplicial complexes in one specific dimension dominate significantly the Betti numbers in all other dimensions. In this paper we focus mainly on the properties of fundamental groups of multi-parameter random simplicial complexes, which can be viewed as a new class of random groups. We describe thresholds for nontrivially and hyperbolicity (in the sense of Gromov) for these groups. Besides, we find domains in the multi-parameter space where these groups have 2-torsion. We also prove that these groups have never odd-prime torsion and their geometric and cohomological dimensions are either 0,1, 2 or infinity. Another result presented in this paper states that aspherical 2-dimensional subcomplexes of random complexes satisfy the Whitehead Conjecture, i.e. all their subcomplexes are also aspherical (with probability tending to one).

preprint2015arXiv

Large Random Simplicial Complexes, III; The Critical Dimension

In this paper we study the notion of critical dimension of random simplicial complexes in the general multi-parameter model described in our previous papers of this series. This model includes as special cases the Linial-Meshulam-Wallach model as well as the clique complexes of random graphs. We characterise the concept of critical dimension in terms of various geometric and topological properties of random simplicial complexes such as their Betti numbers, the fundamental group, the size of minimal cycles and the degrees of simplexes. We mention in the text a few interesting open questions.

preprint2015arXiv

Random Simplicial Complexes

In this paper we introduce a new model of random simplicial complexes depending on multiple probability parameters. This model includes the well-known Linial - Meshulam random simplicial complexes and random clique complexes as special cases. Topological and geometric properties of a multi-parameter random simplicial complex depend on the whole combination of the probability parameters and the thresholds for topological properties are convex sets rather than numbers (as in all previously known models). We discuss the containment properties, density domains and dimension of the random simplicial complexes.

preprint2015arXiv

Two step chromospheric Moreton wave excitation in a blast-wave scenario. A case study: the December 06, 2006 event

We examine the capability of a blast-wave scenario -associated to a coronal flare or to the expansion of CME flanks- to reproduce a chromospheric Moreton phenomenon. We also simulate the Moreton event of December 06, 2006 considering both the corona and the chromosphere. To obtain a sufficiently strong coronal shock -able to generate a detectable chromospheric Moreton wave- a relatively low magnetic field intensity is required, in comparison with the active region values. Employing reasonable coronal constraints, we show that a flare ignited blast-wave or the expansion of the CME flanks emulated as an instantaneous or a temporal piston model, respectively, are capable to reproduce the observations.

preprint2013arXiv

Muon tomography imaging algorithms for nuclear threat detection inside large volume containers with the Muon Portal detector

Muon tomographic visualization techniques try to reconstruct a 3D image as close as possible to the real localization of the objects being probed. Statistical algorithms under test for the reconstruction of muon tomographic images in the Muon Portal Project are here discussed. Autocorrelation analysis and clustering algorithms have been employed within the context of methods based on the Point Of Closest Approach (POCA) reconstruction tool. An iterative method based on the log-likelihood approach was also implemented. Relative merits of all such methods are discussed, with reference to full Geant4 simulations of different scenarios, incorporating medium and high-Z objects inside a container.

preprint2013arXiv

On the sensitivity of extrasolar mass-loss rate ranges: HD 209458b a case study

We present a 3D hydrodynamic study of the effects that different stellar wind conditions and planetary wind structures have on the calculated Ly-$α$ absorptions produced during the transit of HD 209458b. Considering a range of stellar wind speeds $\sim$[350-800] km s$^{-1}$, coronal temperature $\sim$[3-7] $\times10^{6}$ K and two values of the polytropic index $Γ$ $\sim$[1.01-1.13], while keeping fixed the stellar mass loss rate, we found a that a $\dot M_p$ range between $\sim$[3-5] $\times 10^{10}$g s$^{-1}$ give account for the observational absorption in Ly-$α$ measured for the planetary system. Also, several models with anisotropic evaporation profiles for the planetary escaping atmosphere were carried out, showing that both, the escape through polar regions and through the night side yields larger absorptions than an isotropic planetary wind.

preprint2012arXiv

Simulation of descending multiple supra--arcade reconnection outflows in solar flares

After recent AIA observations by Savage, McKenzie and Reeves we revisit the scenario proposed by us in previous papers. We have shown that sunward, generally dark plasma features originated above posteruption flare arcades are consistent with a scenario where plasma voids (which we identify as supra--arcade reconnection outflows, SAROs) generate the bouncing and interfering of shocks and expansion waves upstream of an initial localized deposition of energy which is collimated in the magnetic field direction. In this paper we analyze the multiple production and interaction of SAROs and their individual structure that make them relatively stable features while moving. We compare our results with observations and with the scenarios proposed by other authors.

preprint2011arXiv

Numerical simulation of fundamental trapped sausage modes

Context: We integrate the 2D MHD ideal equations of a straight slab to simulate observational results associated with fundamental sausage trapped modes. Aims: Starting from a non-equilibrium state with a dense chromospheric layer, we analyse the evolution of the internal plasma dynamics of magnetic loops, subject to line-tying boundary conditions, and with the coronal parameters described in Asai et al. (2001) and Melnikov et al. (2002) to investigate the onset and damping of sausage modes. Methods: To integrate the equations we used a high resolution shock-capturing (HRSC) method specially designed to deal appropriately with flow discontinuities. Results: Due to non-linearities and inhomogeneities, pure modes are difficult to sustain and always occur coupled among them so as to satisfy, e.g., the line-tying constraint. We found that, in one case, the resonant coupling of the sausage fundamental mode with a slow one results in a non-dissipative damping of the former. Conclusions: In scenarios of thick and dense loops, where the analytical theory predicts the existence of fundamental trapped sausage modes, the coupling of fast and slow quasi-periodic modes -with a node at the center of the longitudinal speed- occur contributing to the damping of the fast mode. If a discontinuity in the total pressure between the loop and the corona is assumed, a fundamental fast sausage transitory leaky regime is spontaneously produced and an external compressional Alfvén wave takes away the magnetic energy.

preprint2011arXiv

Simulation of dark lanes in post--flare supra--arcades

Using two simulations of 1.5D, for the first time, in Costa et al. (2009) and Shulz et al. (2010) we numerically reproduce the observational dark inflows described in Verwichte et al. (2005). We show that the dark tracks can be explained as hot plasma vacuums generated upstream of a slow magnetoacoustic shock wave produced by a localized deposition of energy. In this work, we show that the overall 2D results are in agreement with the observational behaviour. However they show a slight shift in the characteristic parameter with respect to those found previously. Also, we confirm qualitatively the behaviour found in the previous papers, i.e. for a given numerical domain the period of the kink--like structure is a function of the magnetic field intensity: larger periods are associated with lower magnetic field intensities. Contrary to the 1D result -where the sunward dynamic is independent of the magnetic field intensity due to its exclusively waveguide role- in the 2D simulation the sunward speed is larger for larger values of the magnetic field. This can be interpreted as the capability of the low coronal plasma to collimate the deposition of energy into the magnetic field direction. The moving features consistent of low--density and high--temperature plasma cavities have larger inside values of the structuring parameter beta than the neighboring media. Thus, the voids seem to be the emergence structures of a whole nonlinear interacting plasma context of shocks and waves more than voided plasma loops magnetically structured.

preprint2011arXiv

VisIVOWeb: A WWW Environment for Large-Scale Astrophysical Visualization

This article presents a newly developed Web portal called VisIVOWeb that aims to provide the astrophysical community with powerful visualization tools for large-scale data sets in the context of Web 2.0. VisIVOWeb can effectively handle modern numerical simulations and real-world observations. Our open-source software is based on established visualization toolkits offering high-quality rendering algorithms. The underlying data management is discussed with the supported visualization interfaces and movie-making functionality. We introduce VisIVOWeb Network, a robust network of customized Web portals for visual discovery, and VisIVOWeb Connect, a lightweight and efficient solution for seamlessly connecting to existing astrophysical archives. A significant effort has been devoted for ensuring interoperability with existing tools by adhering to IVOA standards. We conclude with a summary of our work and a discussion on future developments.

preprint2010arXiv

Discontinuous Transition in a Boundary Driven Contact Process

The contact process is a stochastic process which exhibits a continuous, absorbing-state phase transition in the Directed Percolation (DP) universality class. In this work, we consider a contact process with a bias in conjunction with an active wall. This model exhibits waves of activity emanating from the active wall and, when the system is supercritical, propagating indefinitely as travelling (Fisher) waves. In the subcritical phase the activity is localised near the wall. We study the phase transition numerically and show that certain properties of the system, notably the wave velocity, are discontinuous across the transition. Using a modified Fisher equation to model the system we elucidate the mechanism by which the the discontinuity arises. Furthermore we establish relations between properties of the travelling wave and DP critical exponents.

preprint2010arXiv

VisIVO - Integrated Tools and Services for Large-Scale Astrophysical Visualization

VisIVO is an integrated suite of tools and services specifically designed for the Virtual Observatory. This suite constitutes a software framework for effective visual discovery in currently available (and next-generation) very large-scale astrophysical datasets. VisIVO consists of VisiVO Desktop - a stand alone application for interactive visualization on standard PCs, VisIVO Server - a grid-enabled platform for high performance visualization and VisIVO Web - a custom designed web portal supporting services based on the VisIVO Server functionality. The main characteristic of VisIVO is support for high-performance, multidimensional visualization of very large-scale astrophysical datasets. Users can obtain meaningful visualizations rapidly while preserving full and intuitive control of the relevant visualization parameters. This paper focuses on newly developed integrated tools in VisIVO Server allowing intuitive visual discovery with 3D views being created from data tables. VisIVO Server can be installed easily on any web server with a database repository. We discuss briefly aspects of our implementation of VisiVO Server on a computational grid and also outline the functionality of the services offered by VisIVO Web. Finally we conclude with a summary of our work and pointers to future developments.

preprint2009arXiv

Numerical simulation of dark lanes in post-flare supra-arcade

We integrate the MHD ideal equations to simulate dark void sunwardly moving structures in post--flare supra--arcades. We study the onset and evolution of the internal plasma instability to compare with observations and to gain insight into physical processes and characteristic parameters of these phenomena. The numerical approach uses a finite-volume Harten-Yee TVD scheme to integrate the 1D1/2 MHD equations specially designed to capture supersonic flow discontinuities. The integration is performed in both directions, the sunward radial one and the transverse to the magnetic field. For the first time, we numerically reproduce observational dark voids described in Verwichte et al. (2005). We show that the dark tracks are plasma vacuums generated by the bouncing and interfering of shocks and expansion waves, upstream an initial slow magnetoacoustic shock produced by a localized deposition of energy modeled with a pressure perturbation. The same pressure perturbation produces a transverse to the field or perpendicular magnetic shock giving rise to nonlinear waves that compose the kink--like plasma void structures, with the same functional sunward decreasing phase speed and constancy with height of the period, as those determined by the observations.

preprint2009arXiv

Numerical simulation of the internal plasma dynamics of post-flare loops

We integrate the MHD ideal equations of a slender flux tube to simulate the internal plasma dynamics of coronal post-flare loops. We study the onset and evolution of the internal plasma instability to compare with observations and to gain insight into physical processes and characteristic parameters associated with flaring events. The numerical approach uses a finite-volume Harten-Yee TVD scheme to integrate the 1D1/2 MHD equations specially designed to capture supersonic flow discontinuities. We could reproduce the observational sliding down and upwardly propagating of brightening features along magnetic threads of an event occurred on October 1st, 2001. We show that high--speed downflow perturbations, usually interpreted as slow magnetoacoustic waves, could be better interpreted as slow magnetoacoustic shock waves. This result was obtained considering adiabaticity in the energy balance equation. However, a time--dependent forcing from the basis is needed to reproduce the reiteration of the event which resembles observational patterns -commonly known as quasi--periodic pulsations (QPPs)- which are related with large scale characteristic longitudes of coherence. This result reinforces the interpretation that the QPPs are a response to the pulsational flaring activity.

preprint2005arXiv

Viscosity of high crystal content melts: dependence on solid fraction

The rheological properties of suspensions containing high solid fractions are investigated. Attention is focused on viscosity of silicate and magmatic melt systems. A general empirical equation which describes the relative viscosity of suspensions as a function of suspended solid fraction is proposed. In the limit of very dilute solid concentrations it reduces to the Einstein equation. The proposed relationship is satisfactorily applied to reproduce available experimental data relative to silicate melts. Moreover, the extrapolation of the model to very high concentrations is compared with experimental observations on partially-melted granite.

preprint2003arXiv

Viscous heating in fluids with temperature-dependent viscosity: implications for magma flows

Viscous heating plays an important role in the dynamics of fluids with strongly temperature-dependent viscosity because of the coupling between the energy and momentum equations. The heat generated by viscous friction produces a local temperature increase near the tube walls with a consequent decrease of the viscosity which may dramatically change the temperature and velocity profiles. These processes are mainly controlled by the Peclét number, the Nahme number, the flow rate and the thermal boundary conditions. The problem of viscous heating in fluids was investigated in the past for its practical interest in the polymer industry, and was invoked to explain some rheological behaviours of silicate melts, but was not completely applied to study magma flows. In this paper we focus on the thermal and mechanical effects caused by viscous heating in tubes of finite lengths. We find that in magma flows at high Nahme number and typical flow rates, viscous heating is responsible for the evolution from Poiseuille flow, with a uniform temperature distribution at the inlet, to a plug flow with a hotter layer near the walls. When the temperature gradients induced by viscous heating are very pronounced, local instabilities may occur and the triggering of secondary flows is possible. For completeness, this paper also describes magma flow in infinitely long tubes both at steady state and in transient phase.