-
Space Plasma Physics Science Opportunities for the Lunar Orbital Platform -Gateway
Authors:
Iannis Dandouras,
Matt G G T Taylor,
Johan de Keyser,
Yoshifumi Futaana,
Ruth A Bamford,
Graziella Branduardi-Raymont,
Jean-Yves Chaufray,
Dragos Constantinescu,
Elisabetta de Angelis,
Pierre Devoto,
Jonathan Eastwood,
Marius Echim,
Philippe Garnier,
Benjamin Grison,
David Hercik,
Helmut Lammer,
André Laurens,
François Leblanc,
Anna Milillo,
Rumi Nakamura,
Lubomír Přech,
Elias Roussos,
Štěpán Štverák,
Julien Forest,
Arnaud Trouche
, et al. (4 additional authors not shown)
Abstract:
The Lunar Orbital Platform-Gateway (LOP-Gateway, or simply Gateway) is a crewed platform that will be assembled and operated in the vicinity of the Moon by NASA and international partner organizations, including ESA, starting from the mid-2020s. It will offer new opportunities for fundamental and applied scientific research. The Moon is a unique location to study the deep space plasma environment.…
▽ More
The Lunar Orbital Platform-Gateway (LOP-Gateway, or simply Gateway) is a crewed platform that will be assembled and operated in the vicinity of the Moon by NASA and international partner organizations, including ESA, starting from the mid-2020s. It will offer new opportunities for fundamental and applied scientific research. The Moon is a unique location to study the deep space plasma environment. Moreover, the lunar surface and the surface-bounded exosphere are interacting with this environment, constituting a complex multi-scale interacting system. This paper examines the opportunities provided by externally mounted payloads on the Gateway in the field of space plasma physics, heliophysics and space weather, but also examines the impact of the space environment on an inhabited platform in the vicinity of the Moon. It then presents the conceptual design of a model payload, required to perform these space plasma measurements and observations. It results that the Gateway is very well-suited for space plasma physics research. It allows a series of scientific objectives with a multidisciplinary dimension to be addressed.
△ Less
Submitted 20 December, 2022;
originally announced January 2023.
-
Detecting exoplanets with FAST?
Authors:
Philippe Zarka,
Di Li,
Jean-Mathias Grießmeier,
Laurent Lamy,
Julien N. Girard,
Sébastien L. G. Hess,
T. Joseph W. Lazio,
Gregg Hallinan
Abstract:
We briefly review the various proposed scenarios that may lead to nonthermal radio emissions from exoplanetary systems (planetary magnetospheres, magnetosphere-ionosphere and magnetosphere-satellite coupling, and star-planet interactions), and the physical information that can be drawn from their detection. The latter scenario is especially favorable to the production of radio emission above 70\,M…
▽ More
We briefly review the various proposed scenarios that may lead to nonthermal radio emissions from exoplanetary systems (planetary magnetospheres, magnetosphere-ionosphere and magnetosphere-satellite coupling, and star-planet interactions), and the physical information that can be drawn from their detection. The latter scenario is especially favorable to the production of radio emission above 70\,MHz. We summarize the results of past and recent radio searches, and then discuss FAST characteristics and observation strategy, including synergies. We emphasize the importance of polarization measurements and a high duty-cycle for the very weak targets that radio-exoplanets prove to be.
△ Less
Submitted 2 April, 2019;
originally announced April 2019.
-
ExPRES: a Tool to Simulate Exoplanetary and Planetary Radio Emissions
Authors:
C. K. Louis,
S. L. G. Hess,
B. Cecconi,
P. Zarka,
L. Lamy,
S. Aicardi,
A. Loh
Abstract:
All magnetized planets are known to produce intense non thermal radio emissions through a mechanism known as Cyclotron Maser Instability (CMI), requiring the presence of accelerated electrons generally arising from magnetospheric current systems. In return, radio emissions are a good probe of these current systems and acceleration processes. The CMI generates highly anisotropic emissions and leads…
▽ More
All magnetized planets are known to produce intense non thermal radio emissions through a mechanism known as Cyclotron Maser Instability (CMI), requiring the presence of accelerated electrons generally arising from magnetospheric current systems. In return, radio emissions are a good probe of these current systems and acceleration processes. The CMI generates highly anisotropic emissions and leads to important visibility effects, which have to be taken into account when interpreting the data. Several studies showed that modeling the radio source anisotropic beaming pattern can reveal a wealth of physical information about the planetary or exoplanetary magnetospheres that produce these emissions. We present a numerical tool, called ExPRES (Exoplanetary and Planetary Radio Emission Simulator), which is able to reproduce the occurrence in time-frequency plane of CMI-generated radio emissions from planetary magnetospheres, exoplanets or star-planet interacting systems. Special attention is given to the computation of the radio emission beaming at and near its source. We explain what physical information about the system can be drawn from such radio observations, and how it is obtained. These information may include the location and dynamics of the radio sources, the type of current system leading to electron acceleration and their energy and, for exoplanetary systems, the magnetic field strength, the orbital period of the emitting body and the rotation period, tilt and offset of the planetary magnetic field. Most of these parameters can be remotely measured only via radio observations. The ExPRES code provides the proper framework of analysis and interpretation for past (Cassini, Voyager, Galileo), current (Juno, groundbased radiotelescopes) and future (BepiColombo, Juice) observations of planetary radio emissions, as well as for future detection of radio emissions from exoplanetary systems.
△ Less
Submitted 31 January, 2019;
originally announced January 2019.
-
Simulating Jupiter-satellite decametric emissions with ExPRES: a parametric study
Authors:
C. K. Louis,
L. Lamy,
P. Zarka,
B. Cecconi,
S. L. G. Hess,
X. Bonnin
Abstract:
The high latitude radio emissions produced by the Cyclotron Maser Instability (CMI) in Jupiter's magnetosphere extend from a few kHz to 40 MHz. Part of the decametric emissions is of auroral origin, and part is driven by the moons Io, Europa and Ganymede. After summarizing the method used to identify Jupiter-satellite radio emissions, which consists in comparing space- and ground-based radio obser…
▽ More
The high latitude radio emissions produced by the Cyclotron Maser Instability (CMI) in Jupiter's magnetosphere extend from a few kHz to 40 MHz. Part of the decametric emissions is of auroral origin, and part is driven by the moons Io, Europa and Ganymede. After summarizing the method used to identify Jupiter-satellite radio emissions, which consists in comparing space- and ground-based radio observations to ExPRES simulations of CMI-driven emissions in the time-frequency plane, we present a parametric study of the free parameters required by the ExPRES code (electron distribution function and resonant energy, magnetic field model, lead angle, and altitude of the ionospheric cut-off) in order to assess the accuracy of our simulations in the Io-Jupiter case. We find that Io-DAM arcs are fairly modeled by loss-cone driven CMI with electrons of 1-10 keV energy, using the ISaAC, VIPAL or VIP4 magnetic field model and a simple sinusoidal lead angle model. The altitude of the ionospheric cut-off has a marginal impact on the simulations. We discuss the impact of our results on the identification of Europa-DAM and Ganymede-DAM emissions.
△ Less
Submitted 27 April, 2018;
originally announced April 2018.
-
Refurbishing Voyager 1 & 2 Planetary Radio Astronomy (PRA) Data
Authors:
Baptiste Cecconi,
Anaïs Pruvot,
Laurent Lamy,
Philippe Zarka,
Corentin Louis,
Sébastien L. G. Hess,
Doc R. Evans,
Danièle Boucon
Abstract:
Voyager/PRA (Planetary Radio Astronomy) data from digitized tapes archived at CNES have been reprocessed and recalibrated. The data cover the Jupiter and Saturn flybys of both Voyager probes. We have also reconstructed goniopolarimetric datasets (flux and polarization) at full resolution. These datasets are currently not available to the scientific community, but they are of primary interest for t…
▽ More
Voyager/PRA (Planetary Radio Astronomy) data from digitized tapes archived at CNES have been reprocessed and recalibrated. The data cover the Jupiter and Saturn flybys of both Voyager probes. We have also reconstructed goniopolarimetric datasets (flux and polarization) at full resolution. These datasets are currently not available to the scientific community, but they are of primary interest for the analysis of the Cassini data at Saturn, and the Juno data at Jupiter, as well as for the preparation of the JUICE mission. We present the first results derived from the re-analysis of this dataset.
△ Less
Submitted 28 October, 2017;
originally announced October 2017.
-
Imaging Jupiter's radiation belts down to 127 MHz with LOFAR
Authors:
J. N. Girard,
P. Zarka,
C. Tasse,
S. Hess,
I. de Pater,
D. Santos-Costa,
Q. Nenon,
A. Sicard,
S. Bourdarie,
J. Anderson,
A. Asgekar,
M. E. Bell,
I. van Bemmel,
M. J. Bentum,
G. Bernardi,
P. Best,
A. Bonafede,
F. Breitling,
R. P. Breton,
J. W. Broderick,
W. N. Brouw,
M. Brüggen,
B. Ciardi,
S. Corbel,
A. Corstanje
, et al. (49 additional authors not shown)
Abstract:
Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127…
▽ More
Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 ° $\pm$ 25 °. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.
△ Less
Submitted 29 November, 2015;
originally announced November 2015.
-
Extremely Metal-Poor Galaxies: The Environment
Authors:
M. E. Filho,
J. Sánchez Almeida,
C. Muñoz-Tuñón,
S. E. Nuza,
F. Kitaura,
S. Heß
Abstract:
We have analyzed bibliographical observational data and theoretical predictions, in order to probe the environment in which extremely metal-poor dwarf galaxies (XMPs) reside. We have assessed the HI component and its relation to the optical galaxy, the cosmic web type (voids, sheets, filaments and knots), the overdensity parameter and analyzed the nearest galaxy neighbours. The aim is to understan…
▽ More
We have analyzed bibliographical observational data and theoretical predictions, in order to probe the environment in which extremely metal-poor dwarf galaxies (XMPs) reside. We have assessed the HI component and its relation to the optical galaxy, the cosmic web type (voids, sheets, filaments and knots), the overdensity parameter and analyzed the nearest galaxy neighbours. The aim is to understand the role of interactions and cosmological accretion flows in the XMP observational properties, particularly the triggering and feeding of the star formation. We find that XMPs behave similarly to Blue Compact Dwarfs; they preferably populate low-density environments in the local Universe: ~60% occupy underdense regions, and ~75% reside in voids and sheets. This is more extreme than the distribution of irregular galaxies, and in contrast to those regions preferred by elliptical galaxies (knots and filaments). We further find results consistent with previous observations; while the environment does determine the fraction of a certain galaxy type, it does not determine the overall observational properties. With the exception of five documented cases (four sources with companions and one recent merger), XMPs do not generally show signatures of major mergers and interactions; we find only one XMP with a companion galaxy within a distance of 100 kpc, and the HI gas in XMPs is typically well-behaved, demonstrating asymmetries mostly in the outskirts. We conclude that metal-poor accretion flows may be driving the XMP evolution. Such cosmological accretion could explain all the major XMP observational properties: isolation, lack of interaction/merger signatures, asymmetric optical morphology, large amounts of unsettled, metal-poor HI gas, metallicity inhomogeneities, and large specific star formation.
△ Less
Submitted 27 January, 2015;
originally announced January 2015.
-
Cosmic flows and the expansion of the Local Universe from nonlinear phase-space reconstructions
Authors:
Steffen Hess,
Francisco-Shu Kitaura
Abstract:
This work investigates the impact of cosmic flows and density perturbations on Hubble constant $H_0$ measurements using nonlinear phase-space reconstructions of the Local Universe (LU). In particular, we rely on 25 constrained N-body simulations of the LU using the 2MRS galaxy sample within distances of 90 Mpc/h. These have been randomly extended up to volumes enclosing 360 Mpc/h with augmented La…
▽ More
This work investigates the impact of cosmic flows and density perturbations on Hubble constant $H_0$ measurements using nonlinear phase-space reconstructions of the Local Universe (LU). In particular, we rely on 25 constrained N-body simulations of the LU using the 2MRS galaxy sample within distances of 90 Mpc/h. These have been randomly extended up to volumes enclosing 360 Mpc/h with augmented Lagrangian perturbation theory (750 simulations in total), accounting in this way for gravitational mode coupling from larger scales, correcting for periodic boundary effects, and estimating systematics of missing attractors ($σ_{\rm large}=134$ km/s). We report on Local Group speed reconstructions, which are compatible with CMB-dipole measurements: $|v_{\rm LG}|=685\pm137$ km/s. The direction $(l,b)=(260.5\pm 13.3,39.1\pm 10.4)^\circ$ is compatible with observations after considering the variance of large scales. Accounting for large scales, our local bulk flow estimations assuming a $Λ$CDM model are compatible with estimates based on velocity data derived from the Tully-Fisher relation.
We focus on low redshift supernova measurements ($0.01<z<0.025$), which have been found to disagree with probes at larger distances. Our analysis indicates that there are two effects contributing to this tension. First, the anisotropic distribution of supernovae aligns with the velocity dipole and induces a systematic boost in $H_0$. Second, a divergent region surrounding the Virgo Supercluster is responsible for an additional positive bias in $H_0$. Taking these effects into account yields a correction of $ΔH_0=-1.76 \pm 0.21 $ km/s/Mpc, thereby reducing the tension between local probes and more distant probes. Effectively $H_0$ is lower by about $2\%$.
△ Less
Submitted 22 December, 2015; v1 submitted 23 December, 2014;
originally announced December 2014.
-
MultiDark simulations: the story of dark matter halo concentrations and density profiles
Authors:
Anatoly Klypin,
Gustavo Yepes,
Stefan Gottlober,
Francisco Prada,
Steffen Hess
Abstract:
Accurately predicting structural properties of dark matter halos is one of the fundamental goals of modern cosmology. We use the new suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. The MultiDark simulations cover a large range of masses 1e10-1e15Msun and volumes upto 50Gpc**3. The total number of da…
▽ More
Accurately predicting structural properties of dark matter halos is one of the fundamental goals of modern cosmology. We use the new suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. The MultiDark simulations cover a large range of masses 1e10-1e15Msun and volumes upto 50Gpc**3. The total number of dark matter halos in all the simulations exceeds 60 billion. We find that in order to understand the structure of dark matter halos and to make ~1% accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the traditional ratio of the virial radius to the core radius in the NFW profile. For massive halos the averge density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter alpha in the Einasto approximation. Combining results from different redshifts, masses and cosmologies, we show that halos progress through three stages of evolution. (1) They start as rare density peaks that experience very fast and nearly radial infall. This radial infall brings mass closer to the center producing a high concentrated halo. Here, the halo concentration increases with the increasing halo mass and the concentration is defined by the alpha parameter with nearly constant core radius. Later halos slide into (2) the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. Now the concentration does not depend on halo mass. (3) Once the rate of accretion slows down, halos move into the domain of declining concentration-mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. We provide accurate analytical fits to the numerical results for halo density profiles and concentrations.
△ Less
Submitted 3 February, 2016; v1 submitted 14 November, 2014;
originally announced November 2014.
-
The cosmic web of the Local Universe: cosmic variance, matter content and its relation to galaxy morphology
Authors:
Sebastian E. Nuza,
Francisco-Shu Kitaura,
Steffen Hess,
Noam I. Libeskind,
Volker Mueller
Abstract:
We present, for the first time, a Local Universe (LU) characterization using high precision constrained $N$-body simulations based on self-consistent phase-space reconstructions of the large-scale structure in the Two-Micron All-Sky Galaxy Redshift Survey. We analyse whether we live in a special cosmic web environment by estimating cosmic variance from a set of unconstrained $Λ$CDM simulations as…
▽ More
We present, for the first time, a Local Universe (LU) characterization using high precision constrained $N$-body simulations based on self-consistent phase-space reconstructions of the large-scale structure in the Two-Micron All-Sky Galaxy Redshift Survey. We analyse whether we live in a special cosmic web environment by estimating cosmic variance from a set of unconstrained $Λ$CDM simulations as a function of distance to random observers. By computing volume and mass filling fractions for voids, sheets, filaments and knots, we find that the LU displays a typical scatter of about $1σ$ at scales $r>15\,h^{-1}\,{\rm Mpc}$, in agreement with $Λ$CDM, converging to a fair unbiased sample when considering spheres of about $60\,h^{-1}\,{\rm Mpc}$ radius. Additionally, we compute the matter density profile of the LU and found a reasonable agreement with the estimates of Karachentsev (2012) only when considering the contribution of dark haloes. This indicates that observational estimates may be biased towards low density values. As a first application of our reconstruction, we investigate the likelihood of different galaxy morphological types to inhabit certain cosmic web environments. In particular, we find that, irrespective of the method used to define the web, either based on the density or the peculiar velocity field, elliptical galaxies show a clear tendency to preferentially reside in clusters as opposed to voids (up to a level of $5.3σ$ and $9.8σ$ respectively) and conversely for spiral galaxies (up to a level of $5.6σ$ and $5.4σ$ respectively). These findings are compatible with previous works, albeit at higher confidence levels.
△ Less
Submitted 28 August, 2014; v1 submitted 4 June, 2014;
originally announced June 2014.
-
Cosmic variance of the local Hubble flow in large-scale cosmological simulations
Authors:
Radoslaw Wojtak,
Alexander Knebe,
William A. Watson,
Ilian T. Iliev,
Steffen Hess,
David Rapetti,
Gustavo Yepes,
Stefan Gottloeber
Abstract:
The increasing precision in the determination of the Hubble parameter has reached a per cent level at which large-scale cosmic flows induced by inhomogeneities of the matter distribution become non-negligible. Here we use large-scale cosmological N-body simulations to study statistical properties of the local Hubble parameter as measured by local observers. We show that the distribution of the loc…
▽ More
The increasing precision in the determination of the Hubble parameter has reached a per cent level at which large-scale cosmic flows induced by inhomogeneities of the matter distribution become non-negligible. Here we use large-scale cosmological N-body simulations to study statistical properties of the local Hubble parameter as measured by local observers. We show that the distribution of the local Hubble parameter depends not only on the scale of inhomogeneities, but also on how one defines the positions of observers in the cosmic web and what reference frame is used. Observers located in random dark matter haloes measure on average lower expansion rates than those at random positions in space or in the centres of cosmic voids, and this effect is stronger from the halo rest frames compared to the CMB rest frame. We compare the predictions for the local Hubble parameter with observational constraints based on type Ia supernovae (SNIa) and CMB observations. Due to cosmic variance, for observers located in random haloes we show that the Hubble constant determined from nearby SNIa may differ from that measured from the CMB by 0.8 per cent at 1sigma statistical significance. This scatter is too small to significantly alleviate a recently claimed discrepancy between current measurements assuming a flat LCDM model. However, for observers located in the centres of the largest voids permitted by the standard LCDM model, we find that Hubble constant measurements from SNIa would be biased high by 5 per cent, rendering this tension inexistent in this extreme case.
△ Less
Submitted 23 December, 2013; v1 submitted 1 December, 2013;
originally announced December 2013.
-
Halo Abundance Matching: accuracy and conditions for numerical convergence
Authors:
Anatoly Klypin,
Francisco Prada,
Gustavo Yepes,
Steffen Hess,
Stefan Gottlober
Abstract:
Accurate predictions of the abundance and clustering of dark matter haloes play a key role in testing the standard cosmological model. Here, we investigate the accuracy of one of the leading methods of connecting the simulated dark matter haloes with observed galaxies -- the Halo Abundance Matching (HAM) technique. We show how to choose the optimal values of the mass and force resolution in large-…
▽ More
Accurate predictions of the abundance and clustering of dark matter haloes play a key role in testing the standard cosmological model. Here, we investigate the accuracy of one of the leading methods of connecting the simulated dark matter haloes with observed galaxies -- the Halo Abundance Matching (HAM) technique. We show how to choose the optimal values of the mass and force resolution in large-volume N-body simulations so that they provide accurate estimates for correlation functions and circular velocities for haloes and their subhaloes -- crucial ingredients of the HAM method. At the 10% accuracy, results converge for 50 particles for haloes and 150 particles for progenitors subhaloes. In order to achieve this level of accuracy a number of conditions should be satisfied. The force resolution for the smallest resolved (sub)haloes should be in the range (0.1-0.3)rs, where rs is the scale radius of (sub)haloes. The number of particles for progenitors of subhaloes should be 150. We also demonstrate that the two-body scattering plays a minor role for the accuracy of N-body simulations thanks to the relatively small number of crossing-times of dark matter in haloes, and the limited force resolution of cosmological simulations.
△ Less
Submitted 14 October, 2013;
originally announced October 2013.
-
Simulating Structure Formation of the Local Universe
Authors:
Steffen Heß,
Francisco-Shu Kitaura,
Stefan Gottloeber
Abstract:
In this work we present cosmological N-body simulations of the Local Universe with initial conditions constrained by the Two-Micron Redshift Survey (2MRS) within a cubic volume of 180 Mpc/h side-length centred at the Local Group. We use a self-consistent Bayesian based approach to explore the joint parameter space of primordial density fluctuations and peculiar velocity fields, which are compatibl…
▽ More
In this work we present cosmological N-body simulations of the Local Universe with initial conditions constrained by the Two-Micron Redshift Survey (2MRS) within a cubic volume of 180 Mpc/h side-length centred at the Local Group. We use a self-consistent Bayesian based approach to explore the joint parameter space of primordial density fluctuations and peculiar velocity fields, which are compatible with the 2MRS galaxy distribution after cosmic evolution. This method (the KIGEN-code) includes the novel ALPT (Augmented Lagrangian Perturbation Theory) structure formation model which combines second order LPT (2LPT) on large scales with the spherical collapse model on small scales. Furthermore we describe coherent flows with 2LPT and include a dispersion term to model fingers-of-god (fogs). These implementations are crucial to avoid artificial filamentary structures, which appear when using a structure formation model with 2LPT and data with compressed fogs. We assume LCDM cosmology throughout our method. The recovered initial Gaussian fields are used to perform a set of 25 constrained simulations. Statistically this ensemble of simulations is in agreement with a reference set of 25 simulations based on randomly seeded Gaussian fluctuations in terms of matter statistics, power-spectra and mass functions. Considering the entire volume, we obtain correlation coefficients of about 98.3% for the comparison between the simulated density fields and the galaxy density field in log-space with Gaussian smoothing scales of r_S=3.5 Mpc/h (74% for r_S=1.4 Mpc/h). The cross power-spectra show correlations with the galaxy distribution which weaken towards smaller length scales until they vanish at scales of 2.2-3.0 Mpc/h. The simulations we present provide a fully nonlinear density and velocity field with a high level of correlation with the observed galaxy distribution at scales of a few Mpc.
△ Less
Submitted 2 August, 2013; v1 submitted 24 April, 2013;
originally announced April 2013.
-
Cosmological Structure Formation with Augmented Lagrangian Perturbation Theory
Authors:
Francisco-Shu Kitaura,
Steffen Heß
Abstract:
We present a new fast and efficient approach to model structure formation with Augmented Lagrangian Perturbation Theory (ALPT). Our method is based on splitting the displacement field into a long and a short-range component. The long-range component is computed by second order LPT (2LPT). This approximation contains a tidal nonlocal and nonlinear term. Unfortunately, 2LPT fails on small scales due…
▽ More
We present a new fast and efficient approach to model structure formation with Augmented Lagrangian Perturbation Theory (ALPT). Our method is based on splitting the displacement field into a long and a short-range component. The long-range component is computed by second order LPT (2LPT). This approximation contains a tidal nonlocal and nonlinear term. Unfortunately, 2LPT fails on small scales due to severe shell crossing and a crude quadratic behaviour in the low density regime. The spherical collapse (SC) approximation has been recently reported to correct for both effects by adding an ideal collapse truncation. However, this approach fails to reproduce the structures on large scales where it is significantly less correlated with the N-body result than 2LPT or linear LPT (the Zeldovich approximation). We propose to combine both approximations using for the short-range displacement field the SC solution. A Gaussian filter with a smoothing radius r_S is used to separate between both regimes. We use the result of 25 dark matter only N-body simulations to benchmark at z=0 the different approximations: 1st, 2nd, 3rd order LPT, SC and our novel combined ALPT model. This comparison demonstrates that our method improves previous approximations at all scales showing ~25% and ~75% higher correlation than 2LPT with the N-body solution at k = 1 and 2 h Mpc^-1, respectively. We conduct a parameter study to determine the optimal range of smoothing radii and find that the maximum correlation is achieved with r_S = 4 - 5 h^-1 Mpc. This structure formation approach could be used for various purposes, such as setting-up initial conditions for N-body simulations, generating mock galaxy catalogues, cosmic web analysis or for reconstructions of the primordial density fluctuations.
△ Less
Submitted 9 July, 2013; v1 submitted 14 December, 2012;
originally announced December 2012.
-
Cosmic vorticity and the origin of halo spins
Authors:
Noam I Libeskind,
Yehuda Hoffman,
Matthias Steinmetz,
Stefan Gottloeber,
Alexander Knebe,
Steffen Hess
Abstract:
In the standard model of cosmology, structure emerges out of non-rotational flow and the angular momentum of collapsing halos is induced by tidal torques. The growth of halo angular momentum in the linear and quasi-linear phases is associated with a shear, curl-free, flow and it is well described within the linear framework of tidal torque theory (TTT). However, TTT is rendered irrelevant as haloe…
▽ More
In the standard model of cosmology, structure emerges out of non-rotational flow and the angular momentum of collapsing halos is induced by tidal torques. The growth of halo angular momentum in the linear and quasi-linear phases is associated with a shear, curl-free, flow and it is well described within the linear framework of tidal torque theory (TTT). However, TTT is rendered irrelevant as haloes approach turn around and virialization. At that stage the flow field around halos has non-zero vorticity. Using a cosmological simulation, we have examined the importance of the curl of the velocity field (vorticity) in determining halo spin, finding a strong alignment between the two. We have also examined the alignment of vorticity with the principle axes of the shear tensor, finding that it tends to be perpendicular to the axis along which material is collapsing fastest (e1). This behavior is independent of halo masses and cosmic web environment. Our results agree with previous findings on the tendency of halo spin to be perpendicular to e1, and of the spin of (simulated) halos and (observed) galaxies to be aligned with the large-scale structure. Our results imply that angular momentum growth proceeds in two distinct phases. In the first phase angular momentum emerges out of a shear, curl-free, potential flow, as described by TTT. In the second phase, in which haloes approach virialization, the angular momentum emerges out of a vortical flow and halo spin becomes strongly aligned with the vorticity of the ambient flow field.
△ Less
Submitted 6 December, 2012;
originally announced December 2012.
-
Gas stripping and mixing in galaxy clusters: A numerical comparison study
Authors:
Steffen Heß,
Volker Springel
Abstract:
The ambient hot intra-halo gas in clusters of galaxies is constantly fed and stirred by in-falling galaxies, a process that can be studied in detail with cosmological hydro- dynamical simulations. However, different numerical methods yield discrepant predic- tions for crucial hydrodynamical processes, leading for example to different entropy profiles in clusters of galaxies. In particular, the wid…
▽ More
The ambient hot intra-halo gas in clusters of galaxies is constantly fed and stirred by in-falling galaxies, a process that can be studied in detail with cosmological hydro- dynamical simulations. However, different numerical methods yield discrepant predic- tions for crucial hydrodynamical processes, leading for example to different entropy profiles in clusters of galaxies. In particular, the widely used Lagrangian smoothed particle hydrodynamics (SPH) scheme is suspected to strongly damp fluid instabili- ties and turbulence, which are both crucial to establish the thermodynamic structure of clusters. In this study, we test to which extent our recently developed Voronoi particle hydrodynamics (VPH) scheme yields different results for the stripping of gas out of infalling galaxies, and for the bulk gas properties of cluster. We consider both the evolution of isolated galaxy models that are exposed to a stream of intra cluster medium or are dropped into cluster models, as well as non-radiative cosmological sim- ulations of cluster formation. We also compare our particle-based method with results obtained with a fundamentally different discretisation approach as implemented in the moving-mesh code AREPO. We find that VPH leads to noticeably faster stripping of gas out of galaxies than SPH, in better agreement with the mesh-code than with SPH. We show that despite the fact that VPH in its present form is not as accurate as the moving mesh code in our investigated cases, its improved accuracy of gradient estimates makes VPH an attractive alternative to SPH.
△ Less
Submitted 1 August, 2012;
originally announced August 2012.
-
Particle hydrodynamics with tessellation techniques
Authors:
S. Hess,
V. Springel
Abstract:
Lagrangian smoothed particle hydrodynamics (SPH) is a well-established approach to model fluids in astrophysical problems, thanks to its geometric flexibility and ability to automatically adjust the spatial resolution to the clumping of matter. However, a number of recent studies have emphasized inaccuracies of SPH in the treatment of fluid instabilities. The origin of these numerical problems ca…
▽ More
Lagrangian smoothed particle hydrodynamics (SPH) is a well-established approach to model fluids in astrophysical problems, thanks to its geometric flexibility and ability to automatically adjust the spatial resolution to the clumping of matter. However, a number of recent studies have emphasized inaccuracies of SPH in the treatment of fluid instabilities. The origin of these numerical problems can be traced back to spurious surface effects across contact discontinuities, and to SPH's inherent prevention of mixing at the particle level. We here investigate a new fluid particle model where the density estimate is carried out with the help of an auxiliary mesh constructed as the Voronoi tessellation of the simulation particles instead of an adaptive smoothing kernel. This Voronoi-based approach improves the ability of the scheme to represent sharp contact discontinuities. We show that this eliminates spurious surface tension effects present in SPH and that play a role in suppressing certain fluid instabilities. We find that the new `Voronoi Particle Hydrodynamics' described here produces comparable results than SPH in shocks, and better ones in turbulent regimes of pure hydrodynamical simulations. We also discuss formulations of the artificial viscosity needed in this scheme and how judiciously chosen correction forces can be derived in order to maintain a high degree of particle order and hence a regular Voronoi mesh. This is especially helpful in simulating self-gravitating fluids with existing gravity solvers used for N-body simulations.
△ Less
Submitted 26 April, 2010; v1 submitted 3 December, 2009;
originally announced December 2009.
-
The XMM-SSC survey of hard-spectrum XMM-Newton sources 1: optically bright sources
Authors:
M. J. Page,
I. Lehmann,
Th. Boller,
M. G. Watson,
T. Dwelly,
S. Hess,
I. Matute,
N. S. Loaring,
S. Rosen,
H Ziaeepour,
A. Schwope,
G. Lamer,
F. J. Carrera,
J. Tedds,
R. Della Ceca,
P. Severgnini,
R. G. McMahon,
W. Yuan
Abstract:
We present optical and X-ray data for a sample of serendipitous XMM-Newton sources that are selected to have 0.5-2 keV vs 2-4.5 keV X-ray hardness ratios which are harder than the X-ray background. The sources have 2-4.5 keV X-ray flux >= 10^-14 cgs, and in this paper we examine a subsample of 42 optically bright (r < 21) sources; this subsample is 100 per cent spectroscopically identified. All…
▽ More
We present optical and X-ray data for a sample of serendipitous XMM-Newton sources that are selected to have 0.5-2 keV vs 2-4.5 keV X-ray hardness ratios which are harder than the X-ray background. The sources have 2-4.5 keV X-ray flux >= 10^-14 cgs, and in this paper we examine a subsample of 42 optically bright (r < 21) sources; this subsample is 100 per cent spectroscopically identified. All but one of the optical counterparts are extragalactic, and we argue that the single exception, a Galactic M star, is probably a coincidental association. The X-ray spectra are consistent with heavily absorbed power laws (21.8 < log NH < 23.4), and all of them appear to be absorbed AGN. The majority of the sources show only narrow emission lines in their optical spectra, implying that they are type-2 AGN. Only a small fraction of the sources (7/42) show broad optical emission lines, and all of these have NH < 10^23 cm^-2. This implies that ratios of X-ray absorption to optical/UV extinction equivalent to > 100 times the Galactic gas-to-dust ratio are rare in AGN absorbers (at most a few percent of the population), and may be restricted to broad absorption-line QSOs. Seven objects appear to have an additional soft X-ray component in addition to the heavily absorbed power law. We consider the implications of our results in the light of the AGN unified scheme. We find that the soft components in narrow-line objects are consistent with the unified scheme provided that > 4 per cent of broad-line AGN have ionised absorbers that attenuate their soft X-ray flux by >50 per cent. In at least one of the X-ray absorbed, broad-line AGN in our sample the X-ray spectrum requires an ionised absorber, consistent with this picture.
△ Less
Submitted 23 April, 2007;
originally announced April 2007.
-
XMM-Newton spectral properties of the Ultraluminous IRAS Galaxy Mrk 273
Authors:
I. Balestra,
Th. Boller,
L. Gallo,
D. Lutz,
S. Hess
Abstract:
We present a 23 ks XMM-Newton observation of the ULIRG Mrk 273. The hard X-ray spectrum can be modeled by a highly absorbed power law plus an Fe Kalpha emission line. The iron line is broad, suggesting possible superposition of a neutral iron line at 6.4 keV, and a blend of ionized iron lines from Fe XXV and Fe XXVI. Given the relatively short exposure, the three line components can not be singu…
▽ More
We present a 23 ks XMM-Newton observation of the ULIRG Mrk 273. The hard X-ray spectrum can be modeled by a highly absorbed power law plus an Fe Kalpha emission line. The iron line is broad, suggesting possible superposition of a neutral iron line at 6.4 keV, and a blend of ionized iron lines from Fe XXV and Fe XXVI. Given the relatively short exposure, the three line components can not be singularly resolved with high statistical significance: the neutral component is detected at ~2.5sigma and the Fe XXV line at ~2sigma c.l., while for the Fe XXVI line we can only estimate an upper limit. The broad band spectrum requires, in addition to a highly absorbed power law, at least three collisionally ionized plasma components, which may be associated with star-forming regions. The temperatures of the three plasmas are about 0.3, 0.8 and 6 keV, where the highest of the three is sufficient to produce ionized iron emission lines. An alternative interpretation for the origin of the soft emission might also be given in terms of reflection off some photoionized gas, as has been observed in a number of nearby Compton-thick Seyfert 2 galaxies. A hot gas, photoionized by the primary, continuum can also produce ionized iron lines. Unfortunately, given the limited statistics and the lack of high resolution spectroscopy, it is not possible to distinguish between the two models investigated. We further compare the XMM-Newton findings with the Chandra data obtaining consistent spectral results. The absorption corrected hard X-ray (2-10 keV) luminosity of Mrk 273 corresponds to ~0.2% of the far-IR luminosity, similar to typical values found in pure starbursts. The thermal contribution to the soft X-ray luminosity is ~0.2-0.7x10^42 erg s^-1, comparable to those found in NGC 6240 and other starburst dominated ULIRGs.
△ Less
Submitted 13 July, 2005;
originally announced July 2005.