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.

Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and extragalactic surveys will be carried out simultaneously, so the combined target density will strongly vary. In this paper, we describe a new tiling algorithm that can naturally deal with the large target density variations on the sky and which automatically handles the different exposure times of targets. The tiling pattern is modelled as a marked point process, which is characterised by a probability density that integrates the requirements imposed by the 4MOST survey. The optimal tilling pattern with respect to the defined model is estimated by the tiles configuration that maximises the proposed probability density. In order to achieve this maximisation a simulated annealing algorithm is implemented. The algorithm automatically finds an optimal tiling pattern and assigns a tentative sky brightness condition and exposure time for each tile, while minimising the total execution time that is needed to observe the list of targets in the combined input catalogue of all surveys. Hence, the algorithm maximises the long-term observing efficiency and provides an optimal tiling solution for the survey. While designed for the 4MOST survey, the algorithm is flexible and can with simple modifications be applied to any other multi-object spectroscopic survey.

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.

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.

The decline in abundance of Lyman-$\alpha$ (Ly$\alpha$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$\alpha$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$. Our results mainly confirm previous studies, i.e., we find a declining Ly$\alpha$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$\alpha$ peaks at high redshift, we also analyze the IGM transmission on the blue side, which shows a rapid decline at $z\gtrsim 6$ of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\,\mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 5\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$\alpha$ peak. We discuss our results in the context of the interpretation of observations.

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RAMSES-CUDATON as CoDa I in Ocvirk et al. (2016), CoDa II modified and re-calibrated the subgrid star-formation algorithm, making reionization end earlier, at $z \gtrsim 6$, thereby better matching the observations of intergalactic Lyman-alpha opacity from quasar spectra and electron-scattering optical depth from cosmic microwave background fluctuations. CoDa II predicts a UV continuum luminosity function in good agreement with observations of high-z galaxies, especially at $z = 6$. As in CoDa I, reionization feedback suppresses star formation in haloes below $\sim 2 \times 10^9$ M$_{\odot}$, though suppression here is less severe, a possible consequence of modifying the star-formation algorithm. Suppression is environment-dependent, occurring earlier (later) in overdense (underdense) regions, in response to their local reionization times. Using a constrained realization of $\Lambda$CDM constructed from galaxy survey data to reproduce the large-scale structure and major objects of the present-day Local Universe, CoDa II serves to model both global and local reionization. In CoDa II, the Milky Way and M31 appear as individual islands of reionization, i.e. they were not reionized by the progenitor of the Virgo cluster, nor by nearby groups, nor by each other.

To date at least 10 highly flattened planes of dwarf galaxies are claimed to have been discovered in the Local Universe. The origin of these planes of galaxies remains unknown. One suggestion is that they are related to the large-scale structure of the cosmic web. A recent study found that the normal of a number of these dwarf galaxy planes are very closely aligned with the eigenvector of the shear tensor corresponding to the direction of greatest collapse obtained by reconstructing the full velocity field in the linear regime. Here we extend that work by both considering an additional 5 planes beyond the 5 examined previously and by examining the alignment with respect to the quasi-linear field, a more sophisticated reconstruction, which is a better approximation on smaller (quasi-linear) scales. Our analysis recovers the previous result while not finding a significantly tight alignment with the additional 5 planes. However, the additional 5 planes normals also do not appear to be randomly oriented. We conclude that this could either be due to the normals of the new planes being poorly defined and described; the quasi-linear field at those locations being poorly constrained; or different formation mechanisms for the orientation of planes of dwarf galaxies.

It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the $\Lambda$CDM model. One type of anisotropy recently detected in the SDSS (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called "lopsided". Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set up when comparing pairs of galaxies in $\Lambda$CDM simulations indicating that such a set up is not inconsistent with $\Lambda$CDM, the origin has yet to be explained. Here we examine the origin of such lopsided setups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at $z=0$ display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards $z=0$. Second, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the $z=0$ signal is driven primarily by satellites that are on first approach, who have yet to experience a "flyby". This suggests that the signal seen in the observations is also dominated by dynamically young accretion events.

Galaxy clusters constitute powerful cosmological probes thanks to comparisons between observed and simulated clusters. As such Virgo constitutes a formidable source for detailed observations facilitated by its proximity. However, the diversity of clusters complicates the comparisons on a one-to-one basis. Simulated clusters must be carefully selected, a daunting task since most properties are unknown. Alternatively, lookalikes produced in the proper large scale environment can be used. Additionally, their statistical study give access to the mean properties of the observed cluster including its most probable history as well as its deviation from an average cluster. This paper presents such a statistical study with 200 Virgo-like and 400+ cluster-size random dark matter halos. Only 18%(0.5%) of these random halos comply within 3(2)-sigma with the mean values (radius, velocity dispersion, number of substructures, spin, velocity, concentration, center of mass offset with respect to the spherical center) of Virgo halos at z=0 and abide by a similar merging history up to redshift 4. None are within 1-sigma because of environmentally induced properties (number of substructures and velocity). For further comparisons, random halos are selected to reproduce the mass distribution of the lookalikes to cancel mass bias effects. Redshift 1 appears then as a turning point: random to Virgo-like property ratios are alternatively smaller/larger than 1. This highlights the importance of studying clusters within their proper large scale environment: simulated galaxy population, grandly affected by the cluster history, can then be compared with the observed one in details. Direct lookalikes simplify grandly the challenge.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission.

We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We investigate the impact of different magneto-genesis scenarios, both, primordial and astrophysical, on the propagation of cosmic rays. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements. The distribution of nearby sources causes anisotropy at very high energies, independent of the magnetic field model. We compare our results to the dipole signal measured by the Pierre Auger Observatory. All our models could reproduce the observed dipole amplitude with a pure iron injection composition. This is due to clustering of secondary nuclei in direction of nearby sources of heavy nuclei. A light injection composition is disfavoured by the non-observation of anisotropy at energies of 4 - 8 EeV.

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

Galaxy groups and clusters are formidable cosmological probes. They permit the studying of the environmental effects on galaxy formation. A reliable detection of galaxy groups is an open problem and is important for ongoing and future cosmological surveys. We propose a probabilistic galaxy group detection algorithm based on marked point processes with interactions. The pattern of galaxy groups in a catalogue is seen as a random set of interacting objects. The positions and the interactions of these objects are governed by a probability density. The estimator of the unknown cluster pattern is given by the configuration of objects maximising the proposed probability density. Adopting the Bayesian framework, the proposed probability density is maximised using a simulated annealing algorithm. The method provides "for free" additional information such as the probabilities that a point or two points in the observation domain belong to the cluster pattern, respectively. These supplementary tools allow the construction of tests and techniques to validate and to refine the detection result. To test the feasibility of the proposed methodology, we applied it to the well-studied 2MRS data set. Compared to previously published Friends-of-Friends (FoF) group finders, the proposed Bayesian group finder gives overall similar results. The proposed Bayesian group finder is tested on a galaxy redshift survey, but more detailed analyses are needed to understand the actual capabilities of the algorithm regarding upcoming cosmological surveys. The presented mathematical framework permits adapting it easily for other data sets (in astronomy and in other fields of sciences). In cosmology, one promising application is the detection of galaxy groups in photometric galaxy redshift surveys, while taking into account the full photometric redshift posteriors. (abridged)

The local Universe provides a unique opportunity for testing cosmology and theories of structure formation. To facilitate this opportunity we present a new method for the reconstruction of the quasi-linear matter density and velocity fields from galaxy peculiar velocities and apply it to the Cosmicflows-2 data. The method consists of constructing an ensemble of cosmological simulations, constrained by the standard cosmological model and the observational data. The quasi-linear density field is the geometric mean and variance of the fully non-linear density fields of the simulations. The main nearby clusters (Virgo, Centaurus, Coma), superclusters (Shapley, Perseus-Pisces) and voids (Dipole Repeller) are robustly reconstructed. Galaxies are born biased with respect to the underlying dark matter distribution. Using our quasi-linear framework we demonstrate that the luminosity-weighted density field derived from the 2M++ redshift compilations is non-linearly biased with respect to the matter density field. The bias diminishes in the linear regime.

Galaxy clusters are excellent cosmological probes provided that their formation and evolution within the large scale environment are precisely understood. Therefore studies with simulated galaxy clusters have flourished. However detailed comparisons between simulated and observed clusters and their population - the galaxies - are complicated by the diversity of clusters and their surrounding environment. An original way initiated by Bertschinger as early as 1987, to legitimize the one-to-one comparison exercise down to the details, is to produce simulations constrained to resemble the cluster under study within its large scale environment. Subsequently several methods have emerged to produce simulations that look like the local Universe. This paper highlights one of these methods and its essential steps to get simulations that not only resemble the local Large Scale Structure but also that host the local clusters. It includes a new modeling of the radial peculiar velocity uncertainties to remove the observed correlation between the decreases of the simulated cluster masses and of the amount of data used as constraints with the distance from us. This method has the particularity to use solely radial peculiar velocities as constraints: no additional density constraints are required to get local cluster simulacra. The new resulting simulations host dark matter halos that match the most prominent local clusters such as Coma. Zoom-in simulations of the latter and of a volume larger than the 30 Mpc/h radius inner sphere become now possible to study local clusters and their effects. Mapping the local Sunyaev-Zel'dovich and Sachs-Wolfe effects can follow.

To study the full formation and evolution history of galaxy clusters and their population, high resolution simulations of the latter are flourishing. However comparing observed clusters to the simulated ones on a one-to-one basis to refine the models and theories down to the details is non trivial. The large variety of clusters limits the comparisons between observed and numerical clusters. Simulations resembling the local Universe down to the cluster scales permit pushing the limit. Simulated and observed clusters can be matched on a one-to-one basis for direct comparisons provided that clusters are well reproduced besides being in the proper large scale environment. Comparing random and local-Universe like simulations obtained with differently grouped observational catalogs of peculiar velocities, this paper shows that the grouping scheme used to remove non-linear motions in the catalogs that constrain the simulations affects the quality of the numerical clusters. With a less aggressive grouping scheme - galaxies still falling onto clusters are preserved - combined with a bias minimization scheme, the mass of the dark matter halos, simulacra for 5 local clusters - Virgo, Centaurus, Coma, Hydra and Perseus - is increased by 39% closing the gap with observational mass estimates. Simulacra are found on average in 89% of the simulations, an increase of 5% with respect to the previous grouping scheme. The only exception is Perseus. Since the Perseus-Pisces region is not well covered by the used peculiar velocity catalog, the latest release let us foresee a better simulacrum for Perseus in a near future.

As our closest cluster-neighbor, the Virgo cluster of galaxies is intensely studied by observers to unravel the mysteries of galaxy evolution within clusters. At this stage, cosmological numerical simulations of the cluster are useful to efficiently test theories and calibrate model. However, it is not trivial to select the perfect simulacrum of the Virgo cluster to fairly compare in detail its observed and simulated galaxy populations that are affected by the type and history of the cluster. Determining precisely the properties of Virgo for a later selection of simulated clusters becomes essential. It is still not clear how to access some of these properties such as the past history of the Virgo cluster from current observations. Therefore, directly producing effective simulacra of the Virgo cluster is inevitable. Efficient simulacra of the Virgo cluster can be obtained via simulations that resemble the local Universe down to the cluster scale. In such simulations, Virgo-like halos form in the proper local environment and permit assessing the most probable formation history of the cluster. Studies based on these simulations have already revealed that the Virgo cluster has had a quiet merging history over the last seven gigayears and that the cluster accretes matter along a preferential direction. This paper reveals that in addition such Virgo halos have had on average only one merger larger than about a tenth of their mass at redshift zero within the last four gigayears. This second branch (by opposition to main branch) formed in a given sub-region and merged recently (within the last gigayear). These properties are not shared with a set of random halos within the same mass range. This study extends [truncated]

We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We use constrained initial conditions with the cosmological magnetohydrodynamics code {\sc ENZO}. The resulting models of the distribution of magnetic fields in the local Universe are used in the \crpropa code to simulate the propagation of ultra-high energy cosmic rays. We investigate the impact of six different magneto-genesis scenarios, both primordial and astrophysical, on the propagation of cosmic rays over cosmological distances. Moreover, we study the influence of different source distributions around the Milky Way. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements of ultra-high energy cosmic rays. However, at high energies above the GZK-limit, there is anisotropy caused by the distribution of nearby sources, independent of the magnetic field model. This provides a chance to identify cosmic ray sources with future full-sky measurements and high number statistics at the highest energies. Finally, we compare our results to the dipole signal measured by the Pierre Auger Observatory. All our source models and magnetic field models could reproduce the observed dipole amplitude with a pure iron injection composition. Our results indicate that the dipole is observed due to clustering of secondary nuclei in direction of nearby sources of heavy nuclei. A light injection composition is disfavoured by the non-observation of anisotropy at energies of $4-8 \rm\ EeV$.

The Zone of Avoidance (ZOA), whose emptiness is an artifact of our Galaxy dust, has been challenging observers as well as theorists for many years. Multiple attempts have been made on the observational side to map this region in order to better understand the local flows. On the theoretical side, however, this region is often simply statistically populated with structures but no real attempt has been made to confront theoretical and observed matter distributions. This paper takes a step forward using constrained realizations of the local Universe shown to be perfect substitutes of local Universe-like simulations for smoothed high density peak studies. Far from generating completely `random' structures in the ZOA, the reconstruction technique arranges matter according to the surrounding environment of this region. More precisely, the mean distributions of structures in a series of constrained and random realizations differ: while densities annihilate each other when averaging over 200 random realizations, structures persist when summing 200 constrained realizations. The probability distribution function of ZOA grid cells to be highly overdense is a Gaussian with a 15% mean in the random case, while that of the constrained case exhibits large tails. This implies that areas with the largest probabilities host most likely a structure. Comparisons between these predictions and observations, like those of the Puppis 3 cluster, show a remarkable agreement and allow us to assert the presence of the, recently highlighted by observations, Vela supercluster at about 180 Mpc/h, right behind the thickest dust layers of our Galaxy.

High quality reconstructions of the three dimensional velocity and density fields of the local Universe are essential to study the local Large Scale Structure. In this paper, the Wiener Filter reconstruction technique is applied to galaxy radial peculiar velocity catalogs to understand how the Hubble constant (H0) value and the grouping scheme affect the reconstructions. While H0 is used to derive radial peculiar velocities from galaxy distance measurements and total velocities, the grouping scheme serves the purpose of removing non linear motions. Two different grouping schemes (based on the literature and a systematic algorithm) as well as five H0 values ranging from 72 to 76 km/s/Mpc are selected. The Wiener Filter is applied to the resulting catalogs. Whatever grouping scheme is used, the larger H0 is, the larger the infall onto the local Volume is. However, this conclusion has to be strongly mitigated: a bias minimization scheme applied to the catalogs after grouping suppresses this effect. At fixed H0, reconstructions obtained with catalogs grouped with the different schemes exhibit structures at the proper location in both cases but the latter are more contrasted in the less aggressive scheme case: having more constraints permits an infall from both sides onto the structures to reinforce their overdensity. Such findings highlight the importance of a balance between grouping to suppress non linear motions and preserving constraints to produce an infall onto structures expected to be large overdensities. Such an observation is promising to perform constrained simulations of the local Universe including its massive clusters.

The Wiener Filter (WF) technique enables the reconstruction of density and velocity fields from observed radial peculiar velocities. This paper aims at identifying the optimal design of peculiar velocity surveys within the WF framework. The prime goal is to test the dependence of the quality of the reconstruction on the distribution and nature of data points. Mock datasets, extending to 250 Mpc/h, are drawn from a constrained simulation that mimics the local Universe to produce realistic mock catalogs. Reconstructed fields obtained with these mocks are compared to the reference simulation. Comparisons, including residual distributions, cell-to-cell and bulk velocities, imply that the presence of field data points is essential to properly measure the flows. The fields reconstructed from mocks that consist only of galaxy cluster data points exhibit poor quality bulk velocities. In addition, the quality of the reconstruction depends strongly on the grouping of individual data points into single points to suppress virial motions in high density regions. Conversely, the presence of a Zone of Avoidance hardly affects the reconstruction. For a given number of data points, a uniform sample does not score any better than a sample with decreasing number of data points with the distance. The best reconstructions are obtained with a grouped survey containing field galaxies: Assuming no error, they differ from the simulated field by less than 100 km/s up to the extreme edge of the catalogs or up to a distance of three times the mean distance of data points for non-uniform catalogs. The overall conclusions hold when errors are added.

The mass of the Local Group (LG) is a crucial parameter for galaxy formation theories. However, its observational determination is challenging - its mass budget is dominated by dark matter which cannot be directly observed. To meet this end the posterior distributions of the LG and its massive constituents have been constructed by means of constrained and random cosmological simulations. Two priors are assumed - the LCDM model that is used to set up the simulations and a LG model,which encodes the observational knowledge of the LG and is used to select LG-like objects from the simulations. The constrained simulations are designed to reproduce the local cosmography as it is imprinted onto the Cosmicflows-2 database of velocities. Several prescriptions are used to define the LG model, focusing in particular on different recent estimates of the tangential velocity of M31. It is found that (a) different $v_{tan}$ choices affect the peak mass values up to a factor of 2, and change mass ratios of $M_{M31}$ to $M_{MW}$ by up to 20%; (b) constrained simulations yield more sharply peaked posterior distributions compared with the random ones; (c) LG mass estimates are found to be smaller than those found using the timing argument; (d) preferred MW masses lie in the range of $(0.6 - 0.8)\times10^{12} M_{\odot}$ whereas (e) $M_{M31}$ is found to vary between $(1.0 - 2.0)\times10^{12} M_{\odot}$, with a strong dependence on the $v_{tan}$ values used.

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.

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.

Determining the precise value of the tangential component of the velocity of M31 is a non trivial astrophysical issue, that relies on complicated modeling. This has recently lead to con- flicting estimates, obtained by several groups that used different methodologies and assump- tions. This letter addresses the issue by computing a Bayesian posterior distribution function of this quantity, in order to measure the compatibility of those estimates with LambdaCDM. This is achieved using an ensemble of local group (LG) look-alikes collected from a set of Con- strained Simulations (CSs) of the local Universe, and a standard unconstrained LambdaCDM. The latter allows us to build a control sample of LG-like pairs and to single out the influence of the environment in our results. We find that neither estimate is at odds with LambdaCDM; how- ever, whereas CSs favour higher values of vtan , the reverse is true for estimates based on LG samples gathered from unconstrained simulations, overlooking the environmental element

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

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.

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.

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.

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.

We identify some of the most HI massive and fastest rotating disk galaxies in the local universe with the aim of probing the processes that drive the formation of these extreme disk galaxies. By combining data from the Cosmic Flows project, which has consistently reanalyzed archival galaxy HI profiles, and 3.6$\mu$m photometry obtained with the Spitzer Space Telescope, with which we can measure stellar mass, we use the baryonic Tully-Fisher (BTF) relationship to explore whether these massive galaxies are distinct. We discuss several results, but the most striking is the systematic offset of the HI-massive sample above the BTF. These galaxies have both more gas and more stars in their disks than the typical disk galaxy of similar rotational velocity. The "condensed" baryon fraction, $f_C$, the fraction of the baryons in a dark matter halo that settle either as cold gas or stars into the disk, is twice as high in the HI-massive sample than typical, and almost reaches the universal baryon fraction in some cases, suggesting that the most extreme of these galaxies have little in the way of a hot baryonic component or cold baryons distributed well outside the disk. In contrast, the star formation efficiency, measured as the ratio of the mass in stars to that in both stars and gas, shows no difference between the HI-massive sample and the typical disk galaxies. We conclude that the star formation efficiency is driven by an internal, self-regulating process, while $f_C$ is affected by external factors. We also found that the most massive HI detected galaxies are located preferentially in filaments. We present the first evidence of an environmental effect on galaxy evolution using a dynamical definition of a filament.

This paper involves a data release of the observational campaign: Cosmicflows with Spitzer (CFS). Surface photometry of the 1270 galaxies constituting the survey is presented. An additional ~ 400 galaxies from various other Spitzer surveys are also analyzed. CFS complements the Spitzer Survey of Stellar Structure in Galaxies, that provides photometry for an additional 2352 galaxies, by extending observations to low galactic latitudes (|b|<30 degrees). Among these galaxies are calibrators, selected in K band, of the Tully-Fisher relation. The addition of new calibrators demonstrate the robustness of the previously released calibration. Our estimate of the Hubble constant using supernova host galaxies is unchanged, H0 = 75.2 +/- 3.3 km/s/Mpc. Distance-derived radial peculiar velocities, for the 1935 galaxies with all the available parameters, will be incorporated into a new data release of the Cosmicflows project. The size of the previous catalog will be increased by 20%, including spatial regions close to the Zone of Avoidance.

In order to explore local large-scale structures and velocity fields, accurate galaxy distance measures are needed. We now extend the well-tested recipe for calibrating the correlation between galaxy rotation rates and luminosities -- capable of providing such distance measures -- to the all-sky, space-based imaging data from the Wide-field Infrared Survey Explorer (WISE) W1 ($3.4\mu$m) and W2 ($4.6\mu$m) filters. We find a linewidth to absolute magnitude correlation (known as the Tully-Fisher Relation, TFR) of $\mathcal{M}^{b,i,k,a}_{W1} = -20.35 - 9.56 (\log W^i_{mx} - 2.5)$ (0.54 magnitudes rms) and $\mathcal{M}^{b,i,k,a}_{W2} = -19.76 - 9.74 (\log W^i_{mx} - 2.5)$ (0.56 magnitudes rms) from 310 galaxies in 13 clusters. We update the I-band TFR using a sample 9% larger than in Tully & Courtois (2012). We derive $\mathcal{M}^{b,i,k}_I = -21.34 - 8.95 (\log W^i_{mx} - 2.5)$ (0.46 magnitudes rms). The WISE TFRs show evidence of curvature. Quadratic fits give $\mathcal{M}^{b,i,k,a}_{W1} = -20.48 - 8.36 (\log W^i_{mx} - 2.5) + 3.60 (\log W^i_{mx} - 2.5)^2$ (0.52 magnitudes rms) and $\mathcal{M}^{b,i,k,a}_{W2} = -19.91 - 8.40 (\log W^i_{mx} - 2.5) + 4.32 (\log W^i_{mx} - 2.5)^2$ (0.55 magnitudes rms). We apply an I-band -- WISE color correction to lower the scatter and derive $\mathcal{M}_{C_{W1}} = -20.22 - 9.12 (\log W^i_{mx} - 2.5)$ and $\mathcal{M}_{C_{W2}} = -19.63 - 9.11 (\log W^i_{mx} - 2.5)$ (both 0.46 magnitudes rms). Using our three independent TFRs (W1 curved, W2 curved and I-band), we calibrate the UNION2 supernova Type Ia sample distance scale and derive $H_0 = 74.4 \pm 1.4$(stat) $\pm\ 2.4$(sys) kms$^{-1}$ Mpc$^{-1}$ with 4% total error.

We combine data from the Spitzer Survey for Stellar Structure in Galaxies (S$^4$G), a recently calibrated empirical stellar mass estimator from Eskew et al., and an extensive database of HI spectral line profiles to examine the baryonic Tully-Fisher (BTF) relation. We find 1) that the BTF has lower scatter than the classic Tully-Fisher (TF) relation and is better described as a linear relationship, confirming similar previous results, 2) that the inclusion of a radial scale in the BTF decreases the scatter but only modestly, as seen previously for the TF relation, and 3) that the slope of the BTF, which we find to be $3.5\pm 0.2$ ($\Delta$ log $M_{baryon}/\Delta$ log $v_c$), implies that on average a nearly constant fraction ($\sim 0.4$) of all baryons expected to be in a halo are "condensed" onto the central region of rotationally supported galaxies. The condensed baryon fraction, $M_{baryon}/M_{total}$, is, to our measurement precision, nearly independent of galaxy circular velocity (our sample spans circular velocities, $v_c$, between 60 and 250 km s$^{-1}$, but is extended to $v_c\sim 10$ km s$^{-1}$ using data from the literature). The observed galaxy-to-galaxy scatter in this fraction is generally $\le$ a factor of 2 despite fairly liberal selection criteria. These results imply that cooling and heating processes, such as cold vs. hot accretion, mass loss due to stellar winds, and AGN driven feedback, to the degree that they affect the global galactic properties involved in the BTF, are independent of halo mass for galaxies with $10 < v_c < 250$ km/s and typically introduce no more than a factor of two range in the resulting $M_{baryon}/M_{total}$. Recent simulations by Aumer et al. of a small sample of disk galaxies are in excellent agreement with our data, suggesting that current simulations are capable of reproducing the global properties of individual disk galaxies.

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.

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.

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.

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.

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.

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.