Source author record

Jenny G. Sorce

Jenny G. Sorce 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

25works
2topics
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

25 published item(s)

preprint2025arXiv

Velocity fields and turbulence from cosmic filaments to galaxy clusters

Galaxy clusters are currently the endpoint of the hierarchical structure formation; they form via the accretion of dark matter and cosmic gas from their local environment. In particular, filaments contribute grandly by accreting gas from cosmic matter sheets and underdense regions and feeding it to the galaxy clusters. Along the way, the gas in filaments is shocked and heated, which, together with the velocity structure within the filament, induces swirling and, thus, turbulence. In this work, we study a constrained hydrodynamical simulation replica of the Virgo cluster to characterise the velocity field in the two cosmic filaments connected to the cluster with unprecedented high resolution. First, we conduct a qualitative examination of slices extracted from the simulation. We study the temperature, the velocity field, and derived quantities in longitudinal cuts to study the general structure of the filaments and in transverse cuts to study their inner organisation and connection to cosmic matter sheets and underdense regions. Then, we conduct a quantitative study of velocities in Virgo's filaments by computing the 2D energy spectrum from 1 and 5~Mpc square maps extracted from the slices and centred on the core of the filaments. We show that the velocity field goes from mostly compressive far in the filaments to mostly solenoidal in Virgo's core. Moreover, we observe that the total energy spectrum in the filaments gains in amplitude and steepens towards Virgo.

preprint2022arXiv

Anisotropic satellite accretion onto the Local Group with HESTIA

How the cosmic web feeds halos, and fuels galaxy formation is an open question with wide implications. This study explores the mass assembly in the Local Group within the context of the local cosmography by employing simulations whose initial conditions have been constrained to reproduce the local environment. The goal of this study is to inspect whether the direction of accretion of satellites on to the Milky Way and Andromeda galaxies, is related to the cosmic web. The analysis considers the three high-resolution simulations available in the HESTIA simulation suite, as well as the derived velocity shear and tidal tensors. We notice two eras in the Local Group accretion history, delimited by an epoch around $z \approx 0.7$. We also find that satellites can travel up to $\sim 4$ Mpc, relative to their parent halo before crossing its viral radius $R_{200}$. Finally, we observe a strong alignment of the infall direction with the axis of slowest collapse $\vec{e_3}$ of both tidal and shear tensors, implying satellites of the Local Group originated from one particular region of the cosmic web and were channeled towards us via the process of accretion.This alignment is dominated by the satellites that enter during the early infall era, i.e $z>0.7$.

preprint2022arXiv

Cold and Hot gas distribution around the Milky-Way-M31 system in the HESTIA simulations

Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the HESTIA (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The HESTIA suite comprise of a set of three high-resolution {\sc arepo}-based simulations of the LG, run using the Auriga galaxy formation model. For this paper, we focus only on the $z = 0$ simulation datasets and generate mock skymaps along with a power spectrum analysis to show that the distributions of ions tracing low-temperature gas (HI and SiIII) are more clumpy in comparison to warmer gas tracers (OVI, OVII and OVIII). We compare to the spectroscopic CGM observations of M31 and low-redshift galaxies. HESTIA under-produces the column densities of the M31 observations, but the simulations are consistent with the observations of low-redshift galaxies. A possible explanation for these findings is that the spectroscopic observations of M31 are contaminated by gas residing in the CGM of the Milky Way.

preprint2022arXiv

Hermeian haloes: Field haloes that interacted with both the Milky Way and M31

The Local Group is a unique environment in which to study the astrophysics of galaxy formation. The proximity of the Milky Way and M31 enhances the frequency of interactions of the low-mass halo population with more massive dark matter haloes, which increases their concentrations and strips them of gas and other material. Some low-mass haloes pass through the haloes of the Milky Way or M31 and are either ejected into the field or exchanged between the two primary hosts. We use high resolution gas-dynamical simulations to describe a new class of field haloes that passed through the haloes of both the Milky Way and M31 at early times and are almost twice as concentrated as field haloes that do not interact with the primary pair. These 'Hermeian' haloes are distributed anisotropically at larger distances from the Local Group barycentre than the primary haloes and appear to cluster along the line connecting the Milky Way and M31. Hermeian haloes facilitate the exchange of dark matter, gas, and stars between the Milky Way and M31 and can enhance the star formation rate of the gas in the primary haloes during their interactions with them. We also show that some Hermeian haloes can host galaxies that, because they are embedded in haloes that are more concentrated than regular field haloes, are promising targets for indirect dark matter searches beyond the Milky Way virial radius and can produce signals that are competitive with those of some dwarf galaxies. Hermeian galaxies in the Local Group should be detectable by forthcoming wide-field imaging surveys.

preprint2022arXiv

NIHAO-LG: The uniqueness of Local Group dwarf galaxies

Recent observational and theoretical studies of the Local Group (LG) dwarf galaxies have highlighted their unique star-formation history, stellar metallicity, gas content, and kinematics. We investigate the commonality of these features by comparing constrained LG and field central dwarf halo simulations in the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project. Our simulations, performed with NIHAO-like hydrodynamics which track the evolution of the Milky Way (MW) and M31 along with ~100 dwarfs in the LG, reveal the total gas mass and stellar properties (velocity dispersion, evolution history, etc.) of present-day LG dwarfs to be similar to field systems. However, relative to field galaxies, LG dwarfs have more cold gas in their central parts and more metal-rich gas in the halo stemming from interactions with other dwarfs living in a high-density environment like the LG. Interestingly, the direct impact of massive MW/M31 analogues on the metallicity evolution of LG dwarfs is minimal; LG dwarfs accrete high-metallicity gas mostly from other dwarfs at late times. We have also tested for the impact of metal diffusion on the chemical evolution of LG dwarfs, and found that it does not affect the stellar or gaseous content of LG dwarfs. Our simulations suggest that the stellar components of LG dwarfs offer a unique and unbiased local laboratory for galaxy-formation tests and comparisons, especially against the overall dwarf population in the Universe.

preprint2022arXiv

Reionization time of the Local Group and Local-Group-like halo pairs

Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired halos within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully-coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation and isolation criteria but in other environments, alongside isolated halos within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out given our reionization model. Within this model, halos in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ~625 Myr against 660+/-4 Myr (z~8.25 against 7.87+/-0.02) on average.

preprint2022arXiv

Scattering of Lyα Photons through the Reionizing Intergalactic Medium: I. Spectral Energy Distribution

During reionization, a fraction of galactic Ly$α$ emission is scattered in the intergalactic medium (IGM) and appears as a diffuse light extending megaparsecs from the source. We investigate how to probe the properties of the early galaxies and their surrounding IGM using this scattered light. We create a Monte Carlo algorithm to track individual photons and reproduce several test cases from previous literature. Then, we run our code on the simulated IGM of the CoDaII simulation. We find that the scattered light can leave an observable imprint on the emergent spectrum if collected over several square arcminutes. Scattering can redden the emission by increasing the path lengths of photons, but it can also make the photons bluer by upscattering them according to the peculiar motion of the scatterer. The photons emitted on the far blue side of the resonance appear more extended in both frequency and space compared to those emitted near the resonance. This provides a discriminating feature for the blueward emission, which cannot be constrained from the unscattered light coming directly from the source. The ionization state of the IGM also affects the scattered light spectrum. When the source is in a small HII region, the emission goes through more scatterings in the surrounding HI region regardless of the initial frequency and ends up more redshifted and spatially extended. This can result in a weakening of the scattered light toward high $z$ during reionization. Our results provide a framework for interpreting the scattered light to be measured by high-$z$ integral-field-unit surveys.

preprint2022arXiv

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$).

preprint2020arXiv

Efficiently estimating mean, uncertainty and unconstrained large scale fraction of local Universe simulations with paired fixed fields

Provided a random realization of the cosmological model, observations of our cosmic neighborhood now allow us to build simulations of the latter down to the non-linear threshold. The resulting local Universe models are thus accurate up to a given residual cosmic variance. Namely some regions and scales are apparently not constrained by the data and seem purely random. Drawing conclusions together with their uncertainties involves then statistics implying a considerable amount of computing time. By applying the constraining algorithm to paired fixed fields, this paper diverts the original techniques from their first use to efficiently disentangle and estimate uncertainties on local Universe simulations obtained with random fields. Paired fixed fields differ from random realizations in the sense that their Fourier mode amplitudes are fixed and they are exactly out of phase. Constrained paired fixed fields show that only 20% of the power spectrum on large scales (> tens of megaparsecs) is purely random. Namely 80% of it is partly constrained by the large scale / small scale data correlations. Additionally, two realizations of our local environment obtained with paired fixed fields of the same pair constitute an excellent non-biased average or quasi-linear realization of the latter, namely the equivalent of hundreds of constrained simulations. The variance between these two realizations gives the uncertainty on the achievable local Universe simulations. These two simulations will permit enhancing faster our local cosmic web understanding thanks to a drastically reduced required computational time to appreciate its modeling limits and uncertainties.

preprint2020arXiv

Galactic ionising photon budget during the Epoch of Reionisation in the Cosmic Dawn II simulation

Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmission inside high-mass haloes is bi-modal -- ionized cells are transparent, while neutral cells absorb the photons their stars produce - and the halo escape fraction $f_{esc}$ reflects the balance of star formation rate ("SFR") between these modes. The latter is increasingly prevalent at higher halo mass, driving down $f_{esc}$ (we provide analytical fits to our results), whereas halo escape luminosity, proportional to $f_{esc} \times$SFR, increases with mass. Haloes with dark matter masses within $6.10^{8} M_\odot < M_h < 3.10^{10} M_\odot$ produce $\sim 80$% of the escaping photons at z=7, when the Universe is 50% ionized, making them the main drivers of cosmic reionization. Less massive haloes, though more numerous, have low SFRs and contribute less than 10% of the photon budget then, despite their high $f_{esc}$. High mass haloes are too few and too opaque, contributing $<10$% despite their high SFRs. The dominant mass range is lower (higher) at higher (lower) redshift, as mass function and reionization advance together (e.g. at z$=8.5$, x$_{\rm HI}=0.9$, $M_h < 5.10^9 M_\odot$ haloes contributed $\sim$80%). Galaxies with UV magnitudes $M_{AB1600}$ between $-12$ and $-19$ dominated reionization between z$=6$ and 8.

preprint2020arXiv

The Hestia project: simulations of the Local Group

We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter only simulations, and selected for higher resolution magneto hydrodynamic simulation using the Arepo code. Baryonic physics follows the Auriga model of galaxy formation. The simulations contain a high resolution region of 3-5 Mpc in radius from the Local Group midpoint embedded in the correct cosmographic landscape. Within this region a simulated Local Group consisting of a Milky Way and Andromeda like galaxy forms, whose description is in excellent agreement with observations. The simulated Local Group galaxies resemble the Milky Way and Andromeda in terms of their halo mass, mass ratio, stellar disc mass, morphology separation, relative velocity, rotation curves, bulge-disc morphology, satellite galaxy stellar mass function, satellite radial distribution and in some cases, the presence of a Magellanic cloud like object. Because these simulations properly model the Local Group in their cosmographic context, they provide a testing ground for questions where environment is thought to play an important role.

preprint2020arXiv

The Past History of Galaxy Clusters told by their present neighbors

Galaxy clusters can play a key role in modern cosmology provided their evolution is properly understood. However, observed clusters give us only a single timeframe of their dynamical state. Therefore, finding present observable data of clusters that are well correlated to their assembly history constitutes an inestimable tool for cosmology. Former studies correlating environmental descriptors of clusters to their formation history are dominated by halo mass - environment relations. This paper presents a mass-free correlation between the present neighbor distribution of cluster-size halos and the latter mass assembly history. From the Big Multidark simulation, we extract two large samples of random halos with masses ranging from Virgo to Coma cluster sizes. Additionally, to find the main environmental culprit for the formation history of the Virgo cluster, we compare the Virgo-size halos to 200 Virgo-like halos extracted from simulations that resemble the local Universe. The number of neighbors at different cluster-centric distances permits discriminating between clusters with different mass accretion histories. Similarly to Virgo-like halos, clusters with numerous neighbors within a distance of about 2 times their virial radius experience a transition at z~1 between an active period of mass accretion, relative to the mean, and a quiet history. On the contrary, clusters with few neighbors share an opposite trend: from passive to active assembly histories. Additionally, clusters with massive companions within about 4 times their virial radius tend to have recent active merging histories. Therefore, the radial distribution of cluster neighbors provides invaluable insights into the past history of these objects.

preprint2016arXiv

Constrained Local UniversE Simulations: A Local Group Factory

Near field cosmology is practiced by studying the Local Group (LG) and its neighbourhood. The present paper describes a framework for simulating the near field on the computer. Assuming the LCDM model as a prior and applying the Bayesian tools of the Wiener filter (WF) and constrained realizations of Gaussian fields to the Cosmicflows-2 (CF2) survey of peculiar velocities, constrained simulations of our cosmic environment are performed. The aim of these simulations is to reproduce the LG and its local environment. Our main result is that the LG is likely a robust outcome of the LCDM scenario when subjected to the constraint derived from CF2 data, emerging in an environment akin to the observed one. Three levels of criteria are used to define the simulated LGs. At the base level, pairs of halos must obey specific isolation, mass and separation criteria. At the second level the orbital angular momentum and energy are constrained and on the third one the phase of the orbit is constrained. Out of the 300 constrained simulations 146 LGs obey the first set of criteria, 51 the second and 6 the third. The robustness of our LG factory enables the construction of a large ensemble of simulated LGs. Suitable candidates for high resolution hydrodynamical simulations of the LG can be drawn from this ensemble, which can be used to perform comprehensive studies of the formation of the LG

preprint2016arXiv

Cosmicflows-3

The Cosmicflows database of galaxy distances that in the 2nd edition contained 8,188 entries is now expanded to 17,669 entries. The major additions are 2,257 distances that we have derived from the correlation between galaxy rotation and luminosity with photometry at 3.6 microns obtained with Spitzer Space Telescope and 8,885 distances based on the Fundamental Plane methodology from the 6dFGS collaboration. There are minor augmentations to the Tip of the Red Giant Branch and Type Ia supernova compilations. A zero point calibration of the supernova luminosities give a value for the Hubble Constant of 76.2 +-3.4(r) +-2.7(s) km/s/Mpc. Alternatively, a restriction on the peculiar velocity monopole term representing global infall/outflow implies H_0 = 75 +-2 km/s/Mpc.

preprint2016arXiv

How did the Virgo cluster form?

While the Virgo cluster is the nearest galaxy cluster and therefore the best observed one, little is known about its formation history. In this paper, a set of cosmological simulations that resemble the Local Universe is used to shed the first light on this mystery. The initial conditions for these simulations are constrained with galaxy peculiar velocities of the second catalog of the Cosmicflows project using algorithms developed within the Constrained Local UniversE Simulation project. Boxes of 500 Mpc/h on a side are set to run a series of dark matter only constrained simulations. In each simulation, a unique dark matter halo can be reliably identified as Virgo's counterpart. The properties of these Virgo halos are in agreement at a 10-20% level with the global properties of the observed Virgo cluster. Their zero-velocity masses agree at one-sigma with the observational mass estimate. In all the simulations, the matter falls onto the Virgo objects along a preferential direction that corresponds to the observational filament and the slowest direction of collapse. A study of the mass accretion history of the Virgo candidates reveals the most likely formation history of the Virgo cluster, namely a quiet accretion over the last 7 Gigayears.

preprint2015arXiv

Cosmicflows Constrained Local UniversE Simulations

This paper combines observational datasets and cosmological simulations to generate realistic numerical replicas of the nearby Universe. These latter are excellent laboratories for studies of the non-linear process of structure formation in our neighborhood. With measurements of radial peculiar velocities in the Local Universe (cosmicflows-2) and a newly developed technique, we produce Constrained Local UniversE Simulations (CLUES). To assess the quality of these constrained simulations, we compare them with random simulations as well as with local observations. The cosmic variance, defined as the mean one-sigma scatter of cell-to-cell comparison between two fields, is significantly smaller for the constrained simulations than for the random simulations. Within the inner part of the box where most of the constraints are, the scatter is smaller by a factor 2 to 3 on a 5 Mpc/h scale with respect to that found for random simulations. This one-sigma scatter obtained when comparing the simulated and the observation-reconstructed velocity fields is only 104 +/- 4 km/s i.e. the linear theory threshold. These two results demonstrate that these simulations are in agreement with each other and with the observations of our neighborhood. For the first time, simulations constrained with observational radial peculiar velocities resemble the Local Universe up to a distance of 150 Mpc/h on a scale of a few tens of megaparsecs. When focusing on the inner part of the box, the resemblance with our cosmic neighborhood extends to a few megaparsecs (< 5 Mpc/h). The simulations provide a proper Large Scale environment for studies of the formation of nearby objects.

preprint2015arXiv

Minimization of Biases in Galaxy Peculiar Velocity Catalogs

Galaxy distances and derived radial peculiar velocity catalogs constitute valuable datasets to study the dynamics of the Local Universe. However, such catalogs suffer from biases whose effects increase with the distance. Malmquist biases and lognormal error distribution affect the catalogs. Velocity fields of the Local Universe reconstructed with these catalogs present a spurious overall infall onto the Local Volume if they are not corrected for biases. Such an infall is observed in the reconstructed velocity field obtained when applying the BayesianWiener-Filter technique to the raw second radial peculiar velocity catalog of the Cosmicflows project. In this paper, an iterative method to reduce spurious non-Gaussianities in the radial peculiar velocity distribution, to retroactively derive overall better distance estimates resulting in a minimization of the effects of biases, is presented. This method is tested with mock catalogs. To control the cosmic variance, mocks are built out of different cosmological constrained simulations which resemble the Local Universe. To realistically reproduce the effects of biases, the mocks are constructed to be look-alikes of the second data release of the Cosmicflows project, with respect to the size, distribution of data and distribution of errors. Using a suite of mock catalogs, the outcome of the correction is verified to be affected neither by the added error realization, nor by the datapoint selection, nor by the constrained simulation. Results are similar for the different tested mocks. After correction, the general infall is satisfactorily suppressed. The method allows us to obtained catalogs which together with the Wiener-Filter technique give reconstructions approximating non biased velocity fields at 100-150 km/s (2-3 Mpc/h in terms of linear displacement), the linear theory threshold.

preprint2015arXiv

Properties of Galaxies in the Disc Central Surface Brightness Gap

Intermediate surface brightness (ISB) galaxies are less numerous than their counterparts at high and low surface brightness (HSB and LSB). Investigating ISB characteristics from a sample from the S4G survey, complete down to M_B=-16, we find that they have intermediate stellar, gas and baryonic masses and on average as much gas as stars. They lie on the (baryonic) Tully-Fisher relation between HSBs and LSBs, although they present a higher scatter than the latter. Their stellar to baryonic mass ratios have intermediate values unlike their condensed baryonic fractions. By comparing their environments, as classified by the eigenvalues of the velocity shear tensor of local constrained simulations, ISBs have a 5-10% probability higher (smaller) to be in sheets (filaments) with respect to HSBs and LSBs. Additionally, for galaxies in filaments (with close neighbors), the mass and mu_0 are correlated at 2.5 (2) sigma more than for those in sheets. ISBs live in regions where the divergence of the velocity field is smaller than where HSBs and LSBs live, a result at more than 50% significance. ISBs may exist as an unstable transition state between LSBs and HSBs, the low flow activity environment maximally encouraging their formation. Interaction events altering the central baryon fraction could happen at a lower rate in these less dense environment, whilst in the higher density environments the LSBs are primarily satellite galaxies, whose accretion is sufficiently constrained that it fails to promote them to HSBs.

preprint2015arXiv

The alignment of galaxy spin with the shear field in observations

Tidal torque theory suggests that galaxies gain angular momentum in the linear stage of structure formation. Such a theory predicts alignments between the spin of haloes and tidal shear field. However, non-linear evolution and angular momentum acquisition may alter this prediction significantly. In this paper, we use a reconstruction of the cosmic shear field from observed peculiar velocities combined with spin axes extracted from galaxies within $115\, \mathrm{Mpc} $ ($\sim8000 \, {\mathrm {km}}{\mathrm s}^{-1}$) from 2MRS catalog, to test whether or not galaxies appear aligned with principal axes of shear field. Although linear reconstructions of the tidal field have looked at similar issues, this is the first such study to examine galaxy alignments with velocity-shear field. Ellipticals in the 2MRS sample, show a statistically significant alignment with two of the principal axes of the shear field. In general, elliptical galaxies have their short axis aligned with the axis of greatest compression and perpendicular to the axis of slowest compression. Spiral galaxies show no signal. Such an alignment is significantly strengthened when considering only those galaxies that are used in velocity field reconstruction. When examining such a subsample, a weak alignment with the axis of greatest compression emerges for spiral galaxies as well. This result indicates that although velocity field reconstructions still rely on fairly noisy and sparse data, the underlying alignment with shear field is strong enough to be visible even when small numbers of galaxies are considered - especially if those galaxies are used as constraints in the reconstruction.

preprint2013arXiv

Bimodality of Galaxy Disk Central Surface Brightness Distribution in the Spitzer 3.6 micron band

We report on measurements of the disk central surface brightnesses (mu0) at 3.6 microns for 438 galaxies selected by distance and absolute magnitude cutoffs from the 2350+ galaxies in the Spitzer Survey of Stellar Structure in Galaxies (S4G), one of the largest and deepest homogeneous mid-infrared datasets of nearby galaxies. Our sample contains nearly 3 times more galaxies than the most recent study of the mu0 distribution. We demonstrate that there is a bimodality in the distribution of mu0. Between the low and high surface brightness galaxy regimes there is a lack of intermediate surface brightness galaxies. Caveats invoked in the literature from small number statistics to the knowledge of the environmental influences, and possible biases from low signal to noise data or corrections for galaxy inclination are investigated. Analyses show that the bimodal distribution of mu0 cannot be due to any of these biases or statistical fluctuations. It is highly probable that galaxies settle in two stable modes: a dark matter dominated mode where the dark matter dominates at all radii - this gives birth to low surface brightness galaxies - and a baryonic matter dominated mode where the baryons dominate the dark matter in the central parts - this gives rise to the high surface brightness disks. The lack of intermediate surface brightness objects suggests that galaxies avoid (staying in) a mode where dark matter and baryons are co-dominant in the central parts of galaxies.

preprint2013arXiv

Calibration of the Mid-Infrared Tully-Fisher Relation

Distance measures on a coherent scale around the sky are required to address the outstanding cosmological problems of the Hubble Constant and of departures from the mean cosmic flow. The correlation between galaxy luminosities and rotation rates can be used to determine distances to many thousands of galaxies in a wide range of environments potentially out to 200 Mpc. Mid-infrared (3.6 microns) photometry with the Spitzer Space Telescope is particularly valuable as the source of the luminosities because it provides products of uniform quality across the sky. From a perch above the atmosphere, essentially the total magnitude of targets can be registered in exposures of a few minutes. Extinction is minimal and the flux is dominated by the light from old stars which is expected to correlate with the mass of the targets. In spite of the superior photometry, the correlation between mid-infrared luminosities and rotation rates extracted from neutral hydrogen profiles is slightly degraded from the correlation found with I band luminosities. A color correction recovers a correlation that provides comparable accuracy to that available at I band (~20% 1sigma in an individual distance) while retaining the advantages identified above. Without the color correction the relation between linewidth and [3.6] magnitudes is M^{b,i,k,a}_{[3.6]} = -20.34 - 9.74 (log W_{mx}^{i} -2.5). This description is found with a sample of 213 galaxies in 13 clusters that define the slope and 26 galaxies with Cepheid or tip of the red giant branch distances that define the zero point. A color corrected parameter M_{C_{[3.6]}} is constructed that has reduced scatter: M_{C_{[3.6]}} = -20.34 - 9.13 (log W_{mx}^{i} -2.5). Consideration of the 7 calibration clusters beyond 50 Mpc, outside the domain of obvious peculiar velocities, provides a preliminary Hubble Constant estimate of H_0=74+/-5 km/s/Mpc.

preprint2013arXiv

Cosmicflows-2: The Data

Cosmicflows-2 is a compilation of distances and peculiar velocities for over 8000 galaxies. Numerically the largest contributions come from the luminosity-linewidth correlation for spirals, the TFR, and the related Fundamental Plane relation for E/S0 systems, but over 1000 distances are contributed by methods that provide more accurate individual distances: Cepheid, Tip of the Red Giant Branch, Surface Brightness Fluctuation, SNIa, and several miscellaneous but accurate procedures. Our collaboration is making important contributions to two of these inputs: Tip of the Red Giant Branch and TFR. A large body of new distance material is presented. In addition, an effort is made to assure that all the contributions, our own and those from the literature, are on the same scale. Overall, the distances are found to be compatible with a Hubble Constant H_0 = 74.4 +-3.0 km/s/Mpc. The great interest going forward with this data set will be with velocity field studies. Cosmicflows-2 is characterized by a great density and high accuracy of distance measures locally, falling to sparse and coarse sampling extending to z=0.1.

preprint2013arXiv

Simulations of the Local Universe Constrained by Observational Peculiar Velocities

Peculiar velocities, obtained from direct distance measurements, are data of choice to achieve constrained simulations of the Local Universe reliable down to a scale of a few Megaparsecs. Unlike redshift surveys, peculiar velocities are direct tracers of the underlying gravitational field as they trace both baryonic and dark matter. This paper presents the first attempt to use solely observational peculiar velocities to constrain cosmological simulations of the nearby universe. In order to set up Initial Conditions, a Reverse Zel'dovich Approximation (RZA) is used to displace constraints from their positions at z=0 to their precursors' locations at higher redshifts. An additional new feature replaces original observed radial peculiar velocity vectors by their full 3D reconstructions provided by the Wiener-Filter (WF) estimator. Subsequently, the Constrained Realization of Gaussian fields technique (CR) is applied to build various realizations of the Initial Conditions. The WF/RZA/CR method is first tested on realistic mock catalogs built from a reference simulation similar to the Local Universe. These mocks include errors on peculiar velocities, on data-point positions and a large continuous zone devoid of data in order to mimic galactic extinction. Large scale structures are recovered with a typical accuracy of 5 Megaparsecs/h in position, the best realizations reaching a 2-3 Mpc/h precision, the limit imposed by the RZA linear theory. Then, the method is applied to the first observational radial peculiar velocity catalog of the project Cosmicflows. This paper is a proof of concept that the WF/RZA/CR method can be applied to observational peculiar velocities to successfully build constrained Initial Conditions.

preprint2012arXiv

The Mid-Infrared Tully-Fisher Relation: Calibration of the SNIa Scale and Ho

This paper builds on a calibration of the SNIa absolute distance scale begun with a core of distances based on the correlation between galaxy rotation rates and optical Ic band photometry. This new work extends the calibration through the use of mid-infrared photometry acquired at 3.6 microns with Spitzer Space Telescope. The great virtue of the satellite observations is constancy of the photometry at a level better than 1% across the sky. The new calibration is based on 39 individual galaxies and 8 clusters that have been the sites of well observed SNIa. The new 3.6 micron calibration is not yet as extensively based as the Ic band calibration but is already sufficient to justify a preliminary report. Distances based on the mid-infrared photometry are 2% greater in the mean than reported at Ic band. This difference is only marginally significant. The Ic band result is confirmed with only a small adjustment. Incorporating a 1% decrease in the LMC distance, the present study indicates Ho = 75.2 +/- 3.0 km/s/Mpc.

preprint2012arXiv

The mid-infrared Tully-Fisher relation: Spitzer Surface Photometry

The availability of photometric imaging of several thousand galaxies with the Spitzer Space Telescope enables a mid-infrared calibration of the correlation between luminosity and rotation in spiral galaxies. The most important advantage of the new calibration in the 3.6 micron band, IRAC ch.1, is photometric consistency across the entire sky. Additional advantages are minimal obscuration, observations of flux dominated by old stars, and sensitivity to low surface brightness levels due to favorable backgrounds. Through Spitzer cycle 7 roughly 3000 galaxies had been observed and images of these are available at the Spitzer archive. In cycle 8 a program called Cosmic Flows with Spitzer has been initiated that will increase by 1274 the available sample of spiral galaxies with inclinations greater than 45 degrees from face-on suitable for distance measurements. This paper describes procedures based on the photometry package Archangel that are being employed to analyze both the archival and the new data in a uniform way. We give results for 235 galaxies, our calibrator sample for the Tully-Fisher relation. Galaxy magnitudes are determined with uncertainties held below 0.05 mag for normal spiral systems. A subsequent paper will describe the calibration of the [3.6] luminosity-rotation relation.