Search

We show that the hierarchical assembly of cold dark matter haloes preserves the memory of the initial conditions. Using N-body cosmological simulations, we demonstrate that the present-day spatial distribution and kinematics of objects that formed within early(z≳ 10) protogalactic systems (old stars, satellite galaxies, globular clusters, massive black holes, etc.) depends mostly on the rarity of the peak of the primordial density field to which they originally belonged. Only for objects forming at lower redshifts does the exact formation site within the progenitor halo (e.g. whether near the centre or in an extended disc) become important. In present-day haloes, material from the rarer early peaks is more centrally concentrated and falls off more steeply with radius compared to the overall mass distribution, has a lower velocity dispersion, moves on more radial orbits, and has a more elongated shape. Population II stars that formed within protogalactic haloes collapsing from ≥2.5σ fluctuations would follow today an r−3.5 density profile with a half-light radius of 17 kpc and a velocity anisotropy that increases from isotropic in the inner regions to nearly radial at the halo edge. This agrees well with the radial velocity dispersion profile of Galaxy halo stars from the recent work of Battaglia et al. and with the anisotropic orbits of nearby halo stars

We study the pairing of massive black holes embedded in a massive circumnuclear, rotationally supported disc, until they form a close binary. Using high resolution SPH simulations, we follow the black hole dynamics, and in particular the eccentricity evolution, as a function of the composition in stars and gas of the disc. Binarydisc interaction always leads to orbital decay and, in case of corotating black holes, to orbit circularization. We present also a higher resolution simulation performed using the particlesplitting technique showing that the binary orbital decay is efficient down to a separation of ~ 0.1 pc, comparable to our new resolution limit. We detail the gaseous mass profile bound to each black hole. Double nuclear activity is expected to occur on an estimated timescale of 10 Myrs. © 2008 Copyright International Astronomical Union 2008.

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge., 28 pages, 16 figures, MNRAS accepted. High resolution figures and additional images can be found at http://www.aip.de/People/cscannapieco/aquila/

There are different methods for finding exoplanets such as radial spectral shifts, astrometrical measurements, transits, timing, etc. Gravitational microlensing (including pixel-lensing) is among the most promising techniques with the potential of detecting Earth-like planets at distances about a few astronomical units from their host stars. Here we emphasize the importance of polarization measurements which can help to resolve degeneracies in theoretical models. In particular, the polarization angle could give additional information about the relative position of the lens with respect to the source.

Monthly Notices of the Royal Astronomical Society, 378 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

Modified Newtonian dynamics (MOND) has been shown to be able to fit spiral galaxy rotation curves as well as giving a theoretical foundation for empirically determined scaling relations, such as the Tully-Fisher law, without the need for a dark matter halo. As a complementary analysis, one should investigate whether MOND can also reproduce the dynamics of early-type galaxies (ETGs) without dark matter. As a first step, we here show that MOND can indeed fit the observed central velocity dispersion σ0 of a large sample of ETGs assuming a simple MOND interpolating functions and constant anisotropy. We also show that, under some assumptions on the luminosity dependence of the Sérsic n parameter and the stellar mass-to-light ratio (M/L), MOND predicts a Fundamental Plane for ETGs: a loglinear relation among the effective radius Reff, σ0 and the mean effective intensity 〈Ie〉. However, we predict a tilt between the observed and the MOND Fundamental Planes

Observations of turbulent velocity dispersions in the H i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud-cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H i velocity dispersion profiles and the characteristic value of ∼10 km s−1 observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area ≳10−3 M⊙ yr−1 kpc−2

Warm dark matter (WDM) has been invoked to resolve apparent conflicts of cold dark matter (CDM) models with observations on subgalactic scales. In this work, we provide a new and independent lower limit for the WDM particle mass (e.g. sterile neutrino) through the analysis of image fluxes in gravitationally lensed quasi-stellar objects (QSOs). Starting from a theoretical unperturbed cusp configuration, we analyse the effects of intergalactic haloes in modifying the fluxes of QSO multiple images, giving rise to the so-called anomalous flux ratio. We found that the global effect of such haloes strongly depends on their mass/abundance ratio and it is maximized for haloes in the mass range 106-108 M⊙. This result opens up a new possibility to constrain CDM predictions on small scales and test different warm candidates, since free streaming of WDM particles can considerably dampen the matter power spectrum in this mass range. As a consequence, while a (Λ)CDM model is able to produce flux anomalies at a level similar to those observed, a WDM model, with an insufficiently massive particle, fails to reproduce the observational evidences. Our analysis suggests a lower limit of a few keV (mν∼ 10) for the mass of WDM candidates in the form of a sterile neutrino. This result makes sterile neutrino WDM less attractive as an alternative to CDM, in good agreement with previous findings from Lyman α forest and cosmic microwave background analysis

The phase space coordinates of individual halo stars obtained by Galactic surveys enable the computation of their full 3-dimensional orbits. Spectral analysis of halo orbits can be used to construct "frequency maps" which provide a compact representation of the 6-dimensional phase space distribution function. Frequency maps identify important major orbit families, and the orbital abundances reflect the shape and orientation of the dark matter halo relative to the disk. We apply spectral analysis to halo orbits in a series of controlled simulations of disk galaxies. Although the shape of the simulated halo varies with radius, frequency maps of local samples of halo orbits confined to the inner halo contain most of the information about the global shape of the halo and its major orbit families. Quiescent or adiabatic disk formation results in significant trapping of halo orbits in resonant orbit families (i.e. orbits with commensurable frequencies). If a good estimate of the Galactic potential in the inner halo (within ~50kpc) is available, the appearance of strong, stable resonances in frequency maps of halo orbits will allow us to determine the degree of resonant trapping induced by the disk potential. The locations and strengths of these resonant families are determined both by the global shape of the halo and its distribution function. Identification of such resonances in the Milky Way's stellar halo would therefore provide evidence of an extended period of adiabatic disk growth. If the Galactic potential is not known exactly, a measure of the diffusion rate of large sample of 10^4 halo orbits can help distinguish between the true potential and an incorrect potential. The orbital spectral analysis methods described in this paper provide a strong complementarity to existing methods for constraining the potential of the Milky Way halo and its stellar distribution function (ABRIDGED)., 23 pages, 12 figures, 2 tables. Submitted to MNRAS

We use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas ($T\gtrsim 10^6$ K) filling $> 90$% of the volume. This gas reaches high pressures ($10^4 < P/k_\mathrm{B} < 10^7$ K cm$^{-3}$) due to the combination of SN explosions in the hot, low density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass ($>50$%), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ($\ll 1$%). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity dispersion in HI remains $\lesssim 10$ km s$^{-1}$ in all cases. For peak driving the velocity dispersion in H$_��$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants., 19 pages, 12 figures, 2 tables. Accepted for publication in MNRAS. Minor revisions to match published version

We present a statistical study of substructure within a sample of LCDM clusters and galaxies simulated with up to 25 million particles. With thousands of subhalos per object we can accurately measure their spatial clustering and velocity distribution functions and compare these with observational data. The substructure properties of galactic halos closely resembles those of galaxy clusters with a small scatter in the mass and circular velocity functions. The velocity distribution function is non-Maxwellian and flat topped with a negative kurtosis of about -0.7. Within the virial radius the velocity bias $b=\sigma_{\rm sub}/\sigma_{\rm DM}\sim 1.12 \pm 0.04$, increasing to b > 1.3 within the halo centers. Slow subhalos are much less common, due to physical disruption by gravitational tides early in the merging history. This leads to a spatially anti-biased subhalo distribution that is well fitted by a cored isothermal. Observations of cluster galaxies do not show such biases which we interpret as a limitation of pure dark matter simulations - we estimate that we are missing half of the halo population which has been destroyed by physical overmerging. High resolution hydrodynamical simulations are required to study these issues further. If CDM is correct then the cluster galaxies must survive the tidal field, perhaps due to baryonic inflow during elliptical galaxy formation. Spirals can never exist near the cluster centers and the elliptical galaxies there will have little remaining dark matter. This implies that the morphology-density relation is set {\it before} the cluster forms, rather than a subsequent transformation of disks to S0's by virtue of the cluster environment.

Theoretical studies of structure formation find an inverse proportionality between the concentration of dark matter haloes and virial mass. This trend has been recently confirmed for virial masses Mvir > ~6e12 Msun by the observation of the X-ray emission from the hot halo gas. We present an alternative approach to this problem, exploring the concentration of dark matter haloes over galaxy scales on a sample of 18 early-type systems. Our c-Mvir relation is consistent with the X-ray analysis, extending towards lower virial masses, covering the range from ~4e11 Msun up to 5e12 Msun. A combination of the lensing analysis along with photometric data allows us to constrain the baryon fraction within a few effective radii, which is compared with prescriptions for adiabatic contraction (AC) of the dark matter haloes. We find that the standard methods for AC are strongly disfavored, requiring additional mechanisms -- such as mass loss during the contraction process -- to play a role during the phases following the collapse of the haloes., 12 pages, 9 figures, 2 tables. Accepted for publication in MNRAS

We study the formation of discs via the cooling flow of gas within galactic haloes using smoothed particle hydrodynamic simulations. These simulations resolve mass scales of a few thousand solar masses in the gas component for the first time. Thermal instabilities result in the formation of numerous warm clouds that are pressure confined by the hot ambient halo gas. The clouds fall slowly on to the disc through non-spherical accretion from material flowing preferentially down the angular momentum axis. The rotational velocity of the infalling cold gas decreases as a function of height above the disc, closely resembling that of the extra-planar gas recently observed around the spiral galaxy, NGC 891

Globular clusters (GCs) are rich of millisecond pulsars (MSPs) and might also host single or binary intermediatemass black holes (IMBHs). We simulate 3- and 4-body encounters in order to test the possibility that an IMBH captures a MSP. The newly formed system could be revealed from the timing signal of the MSP, providing an unambiguous measure of the BH mass. In current surveys, the number of expected [IMBH,MSP] binaries in the Milky Way is ~0.1. If next-generation radio telescopes (e.g. SKA) will detect ~10 times more MSPs in GCs, we expect to observe at least one [IMBH,MSP] binary. © 2008 Copyright International Astronomical Union 2008.

Monthly Notices of the Royal Astronomical Society, 409 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We use 1-kpc resolution cosmological Adaptive Mesh Refinement (AMR) simulations of a Virgo-like galaxy cluster to investigate the effect of feedback from supermassive black holes on the mass distribution of dark matter, gas and stars. We compared three different models: (i) a standard galaxy formation model featuring gas cooling, star formation and supernovae feedback, (ii) a ‘quenching' model for which star formation is artificially suppressed in massive haloes and finally (iii) the recently proposed active galactic nucleus (AGN) feedback model of Booth and Schaye. Without AGN feedback (even in the quenching case), our simulated cluster suffers from a strong overcooling problem, with a stellar mass fraction significantly above observed values in M87. The baryon distribution is highly concentrated, resulting in a strong adiabatic contraction (AC) of dark matter. With AGN feedback, on the contrary, the stellar mass in the brightest cluster galaxy (BCG) lies below observational estimates and the overcooling problem disappears. The stellar mass of the BCG is seen to increase with increasing mass resolution, suggesting that our stellar masses converge to the correct value from below. The gas and total mass distributions are in better agreement with observations. We also find a slight deficit (∼10 per cent) of baryons at the virial radius, due to the combined effect of AGN-driven convective motions in the inner parts and shock waves in the outer regions, pushing gas to Mpc scales and beyond. This baryon deficit results in a slight adiabatic expansion of the dark matter distribution that can be explained quantitatively by AC theory

We study the relation between the halo and matter density fields -- commonly termed bias -- in the LCDM framework. In particular, we examine the local model of biasing at quadratic order in the matter density. This model is characterized by parameters b_1 and b_2. Using an ensemble of N-body simulations, we apply several statistical methods to estimate the parameters. We measure halo and matter fluctuations smoothed on various scales and find that the parameters vary with smoothing scale. We argue that, for real-space measurements, owing to the mixing of wavemodes, no scale can be found for which the parameters are independent of smoothing. However, this is not the case in Fourier space. We measure halo power spectra and construct estimates for an effective large-scale bias. We measure the configuration dependence of the halo bispectra B_hhh and reduced bispectra Q_hhh for very large-scale k-space triangles. From this we constrain b_1 and b_2. Using the lowest-order perturbation theory, we find that for B_hhh the best-fit parameters are in reasonable agreement with one another as the triangle scale is varied, but that the fits become poor as smaller scales are included. The same is true for Q_hhh. The best-fit parameters depend on the discreteness correction. This led us to consider halo-mass cross-bispectra. The results from these statistics support our earlier findings. We develop a test to explore the importance of missing higher-order terms in the models. We prove that low-order expansions are not able to correctly model the data, even on scales k_1~0.04 h/Mpc. If robust inferences are to be drawn from galaxy surveys, then accurate models for the full nonlinear matter bispectrum and trispectrum will be essential., Comment: 23 pages, 7 figures; accepted for publication in MNRAS

The search for evidence of past or present life on Mars will require the detection of markers that indicate the presence of life. Because deoxyribonucleic acid (DNA) is found in all known living organisms, it is considered to be a ‘biosignature’ of life. The main function of DNA is the long-term storage of genetic information, which is passed on from generation to generation as hereditary material. The Polymerase Chain Reaction (PCR) is a revolutionary technique which allows a single fragment or a small number of fragments of a DNA molecule to be amplified millions of times, making it possible to detect minimal traces of DNA. The compactness of the contemporary PCR instruments makes routine sample analysis possible with a minimum amount of laboratory space. Furthermore the technique is effective, robust and straightforward. Our goal was to establish a routine for the detection of DNA from micro-organisms using the PCR technique during the EuroGeoMars simulation campaign. This took place at the Mars Society's Mars Desert Research Station (MDRS) in Utah in February 2009 (organized with the support of the International Lunar Exploration Working Group (ILEWG), NASA Ames and the European Space Research and Technology Centre (ESTEC)). During the MDRS simulation, we showed that it is possible to establish a minimal molecular biology lab in the habitat for the immediate on-site analysis of samples by PCR after sample collection. Soil and water samples were taken at different locations and soil depths. The sample analysis was started immediately after the crew returned to the habitat laboratory. DNA was isolated from micro-organisms and used as a template for PCR analysis of the highly conserved ribosomal DNA to identify representatives of the different groups of micro-organisms (bacteria, archaea and eukarya). The PCR products were visualized by agarose gel electrophoresis and documented by transillumination and digital imaging. The microbial diversity in the collected samples was analysed with respect to sampling depth and the presence or absence of vegetation. For the first time, we have demonstrated that it is possible to perform direct on-site DNA analysis by PCR at MDRS, a simulated planetary habitat in an extreme environment that serves as a model for preparation and optimization of techniques to be used for future Mars exploration.

Spiral, fast-rotating galaxies like the Milky Way are the most common type in the Universe. One of the most pressing challenges faced by current models of galaxy formation is the origin of their angular momentum and disk. According to the standard tidal-torque theory the galactic spin is originated by tidal interactions between dark halos around galaxies and neighboring structures in the expanding Universe. We use a large cosmological N-body simulation to study the origin of possible correlations between the merging history and spin of cold dark matter halos. In particular, we examine claims that remnants of major mergers tend to have higher-than-average spins, and find that the effect is driven largely by unrelaxed systems: equilibrium dark matter halos show no significant correlation between spin and merging history. Out-of-equilibrium halos have, on average, higher spin than relaxed systems, suggesting that the virialization process leads to a net decrease in the value of the spin parameter. We present also high-resolution N-body/SPH cosmological simulations including cold gas and dark matter to investigate the processes by which gas loses its angular momentum during the protogalactic collapse phase, leading to simulated disk galaxies that are too compact with respect to the observations. We show that the gas and the dark matter have similar specific angular momenta until a merger event occurs at redshift 2. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and falls quickly into the center. Dynamical friction by small infalling substructures plays a minor role, in contrast to previous claims.

Radiative feedback and magnetic field are understood to have a strong impact on the protostellar collapse. We present high resolution numerical calculations of the collapse of a 1 solar mass dense core in solid body rotation, including both radiative transfer and magnetic field. Using typical parameters for low-mass cores, we study thoroughly the effect of radiative transfer and magnetic field on the first core formation and fragmentation. We show that including the two aforementioned physical processes does not correspond to the simple picture of adding them separately. The interplay between the two is extremely strong, via the magnetic braking and the radiation from the accretion shock., 4 pages, 2 figures ; to appear in "IAU Symposium 270: Computational Star formation", Eds. J. Alves, B. Elmegreen, J. Girart, V. Trimble

Current and future galaxy surveys cover a large fraction of the entire sky with a significant redshift range, and the recent theoretical development shows that general relativistic effects are present in galaxy clustering on very large scales. This trend has renewed interest in the wide angle effect in galaxy clustering measurements, in which the distant-observer approximation is often adopted. Using the full wide-angle formula for computing the redshift-space correlation function, we show that compared to the sample variance, the deviation in the redshift-space correlation function from the simple Kaiser formula with the distant-observer approximation is negligible in galaxy surveys such as the SDSS, Euclid and the BigBOSS, if the theoretical prediction from the Kaiser formula is properly averaged over the survey volume. We also find corrections to the wide-angle formula and clarify the confusion in literature between the wide angle effect and the velocity contribution in galaxy clustering. However, when the FKP method is applied, substantial deviations can be present in the power spectrum analysis in future surveys, due to the non-uniform distribution of galaxy pairs., 17 pages, 11 figures, accepted for publication in MNRAS

The existence of blue straggler stars (BSSs) in dwarf spheroidal galaxies (dSphs) is still an open question. In fact, many BSS candidates have been observed in the Local Group dSphs, but it is unclear whether they are real BSSs or young stars. Shedding light on the nature of these BSS candidates is crucial, in order to understand the star formation history of dSphs. In this paper, we consider BSS candidates in Sculptor and Fornax. In Fornax there are strong hints that the BSS population is contaminated by young stars, whereas in Sculptor there is no clear evidence of recent star formation. We derive the radial and luminosity distribution of BSS candidates from wide field imaging data extending beyond the nominal tidal radius of these galaxies. The observations are compared with the radial distribution of BSSs expected from dynamical simulations. In Sculptor the radial distribution of BSS candidates is consistent with that of red horizontal branch (RHB) stars and is in agreement with theoretical expectations for BSSs generated via mass transfer in binaries. On the contrary, in Fornax the radial distribution of BSS candidates is more concentrated than that of all the considered stellar populations. This result supports the hypothesis that most of BSS candidates in Fornax are young stars and is consistent with previous studies., 14 pages, 10 figures, accepted for publication in MNRAS

We show that a global relation between baryonic mass and virial velocity can be constructed from the scale of dwarf galaxies up to that of rich galaxy clusters. The slope of this relation is close to that expected if dark matter halos form in the standard hierarchical cosmogony and capture a universal baryon fraction, once the details of halo structure and the adiabatic contraction of halos due to cooling gas are taken into account. The scatter and deficiency of baryons within low mass halos ($V_{vir} < 50$ km/s) is consistent with the expected suppression of gas accretion by photo-evaporation due to the cosmic UV background at high redshift. The data are not consistent with significant gas removal from strong supernovae winds unless the velocities of galaxies measured from their gas kinematics are significantly lower than the true halo velocities for objects with $V_{vir} < 100$ km/s. Thus models such as $��$CDM with a steep mass function of halos may find it difficult to reproduce both the baryonic mass-velocity relation presented here whilst at the same time reproducing the flat luminosity/HI function of galaxies. Galaxies hold about 10% of the baryons in the Universe, which is close to the collapsed mass fraction expected within hierarchical models on these scales, suggesting a high efficiency for galaxy formation. Most of the baryons are expected to be evenly distributed between diffuse intergalactic gas in low density environments and the intra-galactic medium within galaxy groups., 8 pages, 4 figures, submitted to MNRAS

We present results from a high resolution cosmological galaxy formation simulation called Mare Nostrum and a ultra-high resimulation of the first 500 million years of a single, Milky Way (MW) sized galaxy. Using the cosmological run, we measure UV luminosity functions and assess their sensitivity to both cosmological parameters and dust extinction. We find remarkably good agreement with the existing data over the redshift range 4 < z < 7 provided we adopt the favoured cosmology (WMAP 5 year parameters) and a self-consistent treatment of the dust. Cranking up the resolution, we then study in detail a z = 9 protogalaxy sitting at the intersection of cold gas filaments. This high-z MW progenitor grows a dense, rapidly spinning, thin disk which undergoes gravitational fragmention. Star formation in the resulting gas clumps rapidly turns them into globular clusters. A far reaching galactic wind develops, co-powered by the protogalaxy and its cohort of smaller companions populating the filaments. Despite such an impressive blow out, the smooth filamentary material is hardly affected at these redshifts.

The coalescence of black hole binaries is a significant source of gravitational wave radiation. The typically asymmetric nature of this emission, which carries linear momentum, can result in the recoil of the black hole remnant with velocities in the range 100 < Vrecoil < 3750 km s−1. The detectability of recoiling massive black holes (MBH) as off-nuclear QSOs is tightly connected with the properties of the host galaxy, which determine the MBH's orbit and fuel reservoir. We present the results of N-body simulations of recoiling MBHs in high-resolution, non-axisymmetric potentials. We find that if the recoil velocities are high enough to reach regions of the galaxy dominated by the generally triaxial dark matter distribution, the return time is significantly extended when compared to a spherical distribution. We also perform simulations of recoiling MBHs traveling in gas merger remnants, where large amounts of gas have been funneled to the central regions, In this case, the MBHs remain within RVrecoil = 1200 km s−1) due to the compactness of the remnant galaxy's nuclear disk. We discuss the implications of both scenarios for detectability.

We investigate the effect of primordial non-Gaussianity of the local f_NL type on the auto- and cross-power spectrum of dark matter haloes using simulations of the LCDM cosmology. We perform a series of large N-body simulations of both positive and negative f_NL, spanning the range between 10 and 100. Theoretical models predict a scale-dependent bias correction ��b(k,f_NL) that depends on the linear halo bias b(M). We measure the power spectra for a range of halo mass and redshifts covering the relevant range of existing galaxy and quasar populations. We show that auto and cross-correlation analyses of bias are consistent with each other. We find that for low wavenumbers with k1.5. We show that a scale-independent bias correction improves the comparison between theory and simulations on smaller scales, where the scale-dependent effect rapidly becomes negligible. The current limits on f_NL from Slosar et al. (2008) come mostly from very large scales k, 13 pages, 10 figures. (v2): substantial changes. added a physically motivated scale-independent bias correction which improves significantly the agreement with the simulations (v3): matches published version

In a previous work (Pichardo et al. 2005), we studied stable configurations for circumstellar discs in eccentric binary systems. We searched for "invariant loops": closed curves (analogous to stable periodic orbits in time-independent potentials) that change shape with the binary orbital phase, as test particles in them move under the influence of the binary potential. This approach allows us to identify stable configurations when pressure forces are unimportant, and dissipation acts only to prevent gas clouds from colliding with one another. We now extend this work to study the main geometrical properties of circumbinary discs. We have studied more than 100 cases with a range in eccentricity 0 .le. e .le. 0.9, and mass ratio 0.1 .le. q .le. 0.9. Although gas dynamics may impose further restrictions, our study sets lower stable bounds for the size of the central hole in a simple and computationally cheap way, with a relation that depends on the eccentricity and mass ratio of the central binary. We extend our previous studies and focus on an important component of these systems: circumbinary discs. The radii for stable orbits that can host gas in circumbinary discs are sharply constrained as a function of the binary's eccentricity. The circumbinary disc configurations are almost circular, with eccentricity e_d < 0.15, but if the mass ratio is unequal the disk is offset from the center of mass of the system. We compare our results with other models, and with observations of specific systems like GG Tauri A, UY Aurigae, HD 98800 B, and Fomalhaut, restricting the plausible parameters for the binary., 12 pages, 10 figures and 6 tables. MNRAS, accepted

An automated search is carried out for microlensing events using a catalogue of 44554 variable superpixel lightcurves derived from our three-year monitoring program of M31. Each step of our candidate selection is objective and reproducible by a computer. Our search is unrestricted, in the sense that it has no explicit timescale cut. So, it must overcome the awkward problem of distinguishing long-timescale microlensing events from long-period stellar variables. The basis of the selection algorithm is the fitting of the superpixel lightcurves to two different theoretical models, using variable star and blended microlensing templates. Only if microlensing is preferred is an event retained as a possible candidate. Further cuts are made with regard to (i) sampling, (ii) goodness of fit of the peak to a Paczynski curve, (iii) consistency of the microlensing hypothesis with the absence of a resolved source, (iv) achromaticity, (v) position in the colour-magnitude diagram and (vi) signal-to-noise ratio. Our results are reported in terms of first-level candidates, which are the most trustworthy, and second-level candidates, which are possible microlensing but have lower signal-to-noise and are more questionable. The pipeline leaves just 3 first-level candidates, all of which have very short full-width half-maximum timescale (, 22 pages, 18 figures, MNRAS, to appear

We develop a subgrid model for the growth of supermassive black holes (BHs) and their associated active galactic nucleus (AGN) feedback in hydrodynamical cosmological simulations. This model transposes previous attempts to describe BH accretion and AGN feedback with the smoothed particle hydrodynamics (SPH) technique to the adaptive mesh refinement framework. It also furthers their development by implementing a new jet-like outflow treatment of the AGN feedback which we combine with the heating mode traditionally used in the SPH approach. Thus, our approach allows one to test the robustness of the conclusions derived from simulating the impact of self-regulated AGN feedback on galaxy formation vis-à-vis the numerical method. Assuming that BHs are created in the early stages of galaxy formation, they grow by mergers and accretion of gas at a Eddington-limited Bondi accretion rate. However this growth is regulated by AGN feedback which we model using two different modes: a quasar-heating mode when accretion rates on to the BHs are comparable to the Eddington rate, and a radio-jet mode at lower accretion rates which not only deposits energy, but also deposits mass and momentum on the grid. In other words, our feedback model deposits energy as a succession of thermal bursts and jet outflows depending on the properties of the gas surrounding the BHs. We assess the plausibility of such a model by comparing our results to observational measurements of the co-evolution of BHs and their host galaxy properties, and check their robustness with respect to numerical resolution. We show that AGN feedback must be a crucial physical ingredient for the formation of massive galaxies as it appears to be able to efficiently prevent the accumulation of and/or expel cold gas out of haloes/galaxies and significantly suppress star formation. Our model predicts that the relationship between BHs and their host galaxy mass evolves as a function of redshift, because of the vigorous accretion of cold material in the early Universe that drives Eddington-limited accretion on to BHs. Quasar activity is also enhanced at high redshift. However, as structures grow in mass and lose their cold material through star formation and efficient BH feedback ejection, the AGN activity in the low-redshift Universe becomes more and more dominated by the radio mode, which powers jets through the hot circumgalactic medium. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

The spatial distribution of matter in clusters of galaxies is mainly determined by the dominant dark matter component; however, physical processes involving baryonic matter are able to modify it significantly. We analyse a set of 500 pc resolution cosmological simulations of a cluster of galaxies with mass comparable to Virgo, performed with the AMR code ramses. We compare the mass density profiles of the dark, stellar and gaseous matter components of the cluster that result from different assumptions for the subgrid baryonic physics and galaxy formation processes. First, the prediction of a gravity-only N-body simulation is compared to that of a hydrodynamical simulation with standard galaxy formation recipes, and then all results are compared to a hydrodynamical simulation which includes thermal active galactic nucleus (AGN) feedback from supermassive black holes (SMBHs). We find the usual effects of overcooling and adiabatic contraction in the run with standard galaxy formation physics, but very different results are found when implementing SMBHs and AGN feedback. Star formation is strongly quenched, producing lower stellar densities throughout the cluster, and much less cold gas is available for star formation at low redshifts. At redshift z= 0 we find a flat density core of radius 10 kpc in both the dark and stellar matter density profiles. We speculate on the possible formation mechanisms able to produce such cores and we conclude that they can be produced through the coupling of different processes: (I) dynamical friction from the decay of black hole orbits during galaxy mergers; (II) AGN-driven gas outflows producing fluctuations of the gravitational potential causing the removal of collisionless matter from the central region of the cluster; (III) adiabatic expansion in response to the slow expulsion of gas from the central region of the cluster during the quiescent mode of AGN activity

Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major merger and the present day properties of the final remnant are typical of early type spiral galaxies, with an i band B/D ~0.65, a disk scale length of 7.2 kpc, g-r = 0.5 mag, an HI line width (W_{20}/2) of 238 km/sec and total magnitude i = -22.4. The key ingredients for the formation of a dominant stellar disk component after a major merger are: i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation during mergers and iii) relative fading of the spheroidal component. The gas fraction of the progenitors' disks does not exceed 25% at z, Comment: MNRAS in press. Minor additions to text. Movie at: http://www.astro.washington.edu/fabio/movies/Merger.mpg

Monthly Notices of the Royal Astronomical Society, 364 (2), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

Correlations between intrinsic shear and the density field on large scales, a potentially important contaminant for cosmic shear surveys, have been robustly detected at low redshifts with bright galaxies in SDSS data. Here we present a more detailed characterization of this effect, which can cause anti-correlations between gravitational lensing shear and intrinsic ellipticity (GI correlations). This measurement uses 36278 Luminous Red Galaxies (LRGs) from the SDSS spectroscopic sample with 0.153sigma detections of the effect for all galaxy subsamples within the SDSS LRG sample; for the 2SLAQ sample, we find a 2sigma detection for a bright subsample, and no detection for a fainter subsample. Fitting formulae are provided for the scaling of the GI correlations with luminosity, transverse separation, and redshift. We estimate contamination in the measurement of sigma_8 for future cosmic shear surveys on the basis of the fitted dependence of GI correlations on galaxy properties. We find contamination to the power spectrum ranging from -1.5 (optimistic) to -33 per cent (pessimistic) for a toy cosmic shear survey using all galaxies to a depth of R=24 using scales l~500. This corresponds to a bias in sigma_8 of Delta sigma_8=-0.004 (optimistic), -0.02 (central), or -0.10 (pessimistic). We provide a prescription for inclusion of this error in cosmological parameter estimation codes. The principal uncertainty is in the treatment of the L, Comment: Matches MNRAS accepted version. Extended/updated data tables and expanded comparison to theory

We present the RHAPSODY-G suite of cosmological hydrodynamic zoom simulations of 10 massive galaxy clusters at the Mvir ∼ 1015 M⊙ scale. These simulations include cooling and subresolution models for star formation and stellar and supermassive black hole feedback. The sample is selected to capture the whole gamut of assembly histories that produce clusters of similar final mass. We present an overview of the successes and shortcomings of such simulations in reproducing both the stellar properties of galaxies as well as properties of the hot plasma in clusters. In our simulations, a long-lived cool-core/non-cool-core dichotomy arises naturally, and the emergence of non-cool cores is related to low angular momentum major mergers. Nevertheless, the cool-core clusters exhibit a low central entropy compared to observations, which cannot be alleviated by thermal active galactic nuclei feedback. For cluster scaling relations, we find that the simulations match well the M500–Y500 scaling of Planck Sunyaev–Zeldovich clusters but deviate somewhat from the observed X-ray luminosity and temperature scaling relations in the sense of being slightly too bright and too cool at fixed mass, respectively. Stars are produced at an efficiency consistent with abundance-matching constraints and central galaxies have star formation rates consistent with recent observations. While our simulations thus match various key properties remarkably well, we conclude that the shortcomings strongly suggest an important role for non-thermal processes (through feedback or otherwise) or thermal conduction in shaping the intracluster medium., Monthly Notices of the Royal Astronomical Society, 470 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

In the last decades, a plethora of in vitro studies with living human cells contributed a vast amount of knowledge about cellular and molecular effects of microgravity. Previous studies focused mostly on the identification of gravity-responsive genes, whereas a multi-platform analysis at an integrative level, which specifically evaluates the extent and robustness of transcriptional response to an altered gravity environment was not performed so far. Therefore, we investigated the stability of gene expression response in non-activated human Jurkat T lymphocytic cells in different gravity environments through the combination of parabolic flights with a suborbital ballistic rocket and 2D clinostat and centrifuge experiments, using strict controls for excluding all possible other factors of influence. We revealed an overall high stability of gene expression in microgravity and identified olfactory gene expression in the chromosomal region 11p15.4 as particularly robust to altered gravity. We identified that classical reference genes ABCA5, GAPDH, HPRT1, PLA2G4A, and RPL13A were stably expressed in all tested gravity conditions and platforms, while ABCA5 and GAPDH were also known to be stably expressed in U937 cells in all gravity conditions. In summary, 10–20% of all transcripts remained totally unchanged in any gravitational environment tested (between 10−4 and 9 g), 20–40% remained unchanged in microgravity (between 10−4 and 10−2 g) and 97–99% were not significantly altered in microgravity if strict exclusion criteria were applied. Therefore, we suppose a high stability of gene expression in microgravity. Comparison with other stressors suggests that microgravity alters gene expression homeostasis not stronger than other environmental factors., Gene expression: potentially low-risk for long-term space missions Approximately 99% of genes in human lymphocytic cells have the same transcription activity on Earth as they do in microgravity environments. An international team led by Oliver Ullrich and Cora Thiel at the University of Zurich in Switzerland used a combination of parabolic aircraft and suborbital rocket flights and ground-based experiments to show how gene expression is stable through a range of gravity conditions. Screening of microarrays identified several reference genes that may serve as standardized controls for future immune cell trials. Their experiments also revealed that the chromosomal region associated with olfactory and taste receptor genes is particularly robust to altered gravity. Because altered transcripts were associated with fast cellular adaptation, the researchers predict that the risk of unusual gene behavior will be quite low during extended space deployments.

Accurate measurements of galaxy masses and sizes are key to tracing galaxy evolution over time. Cosmological zoom-in simulations provide an ideal test bed for assessing the recovery of galaxy properties from observations. Here, we utilize galaxies with M_* ~ 10^10 - 10^(11.5) M_☉ at z ~ 1.7–2 from the MassiveFIRE cosmological simulation suite, part of the Feedback in Realistic Environments (FIRE) project. Using mock multi-band images, we compare intrinsic galaxy masses and sizes to observational estimates. We find that observations accurately recover stellar masses, with a slight average underestimate of ~ 0.06 dex and a ~ 0.15 dex scatter. Recovered half-light radii agree well with intrinsic half-mass radii when averaged over all viewing angles, with a systematic offset of ~ 0.1 dex (with the half-light radii being larger) and a scatter of ~ 0.2 dex. When using color gradients to account for mass-to-light variations, recovered half-mass radii also exceed the intrinsic half-mass radii by ~ 1 dex. However, if not properly accounted for, aperture effects can bias size estimates by ~ 0.1 dex. No differences are found between the mass and size offsets for star-forming and quiescent galaxies. Variations in viewing angle are responsible for ~25% of the scatter in the recovered masses and sizes. Our results thus suggest that the intrinsic scatter in the mass–size relation may have previously been overestimated by ~25%. Moreover, orientation-driven scatter causes the number density of very massive galaxies to be overestimated by ~ 0.5 dex at M_* ~ 10^(11.5) M_☉_.

A conclusive model for the formation of dwarf spheroidal (dSph) galaxies still remains elusive. Owing to their proximity to the massive spirals Milky Way (MW) and M31, various environmental processes have been invoked to explain their origin. In this context, the tidal stirring model postulates that interactions with MW-sized hosts can transform rotationally supported dwarfs, resembling present-day dwarf irregular (dIrr) galaxies, into systems with the kinematic and structural properties of dSphs. Using N-body+SPH simulations, we investigate the dependence of this transformation mechanism on the gas fraction, fgas, in the disk of the progenitor dwarf. Our numerical experiments incorporate for the first time the combined effects of radiative cooling, ram-pressure stripping, star formation, supernova (SN) winds, and a cosmic UV background. For a given orbit inside the primary galaxy, rotationally supported dwarfs with gas fractions akin to those of observed dIrrs (fgas >= 0.5), demonstrate a substantially enhanced likelihood and efficiency of transformation into dSphs relative to their collisionless (fgas = 0) counterparts. We argue that the combination of ram-pressure stripping and SN winds causes the gas-rich dwarfs to respond more impulsively to tides, augmenting their transformation. When fgas >= 0.5, disky dwarfs on previously unfavorable low-eccentricity or large-pericenter orbits are still able to transform. On the widest orbits, the transformation is incomplete; the dwarfs retain significant rotational support, a relatively flat shape, and some gas, naturally resembling transition-type systems. We conclude that tidal stirring constitutes a prevalent evolutionary mechanism for shaping the structure of dwarf galaxies within the currently favored CDM cosmological paradigm., Accepted for publication in ApJ Letters, 8 pages, 2 figures, LaTeX (uses emulateapj.cls)

In this paper we present solutions to three short comings of Smoothed Particles Hydrodynamics (SPH) encountered in previous work when applying it to Giant Impacts. First we introduce a novel method to obtain accurate SPH representations of a planet's equilibrium initial conditions based on equal area tessellations of the sphere. This allows one to imprint an arbitrary density and internal energy profile with very low noise which substantially reduces computation because these models require no relaxation prior to use. As a consequence one can significantly increase the resolution and more flexibly change the initial bodies to explore larger parts of the impact parameter space in simulations. The second issue addressed is the proper treatment of the matter/vacuum boundary at a planet's surface with a modified SPH density estimator that properly calculates the density stabilizing the models and avoiding an artificially low density atmosphere prior to impact. Further we present a novel SPH scheme that simultaneously conserves both energy and entropy for an arbitrary equation of state. This prevents loss of entropy during the simulation and further assures that the material does not evolve into unphysical states. Application of these modifications to impact simulations for different resolutions up to $6.4 \cdot 10^6$ particles show a general agreement with prior result. However, we observe resolution dependent differences in the evolution and composition of post collision ejecta. This strongly suggests that the use of more sophisticated equations of state also demands a large number of particles in such simulations.

We present results from the "Ponos" simulation suite on the early evolution of a massive, $M_{\rm vir}(z=0)=1.2\times 10^{13}$ M$_{\odot}$ galaxy. At $z\gtrsim6$, before feedback from a central supermassive black hole becomes dominant, the main galaxy has a stellar mass $\sim 2\times 10^{9}$ M$_{\odot}$ and a star formation rate $\sim 20$ M$_{\odot}$ yr$^{-1}$. The galaxy sits near the expected main sequence of star-forming galaxies at those redshifts, and resembles moderately star-forming systems observed at $z>5$. The high specific star formation rate results in vigorous heating and stirring of the gas by supernovae feedback, and the galaxy develops a thick and turbulent disc, with gas velocity dispersion $\sim 40$ km s$^{-1}$, rotation to dispersion ratio $\sim 2$, and with a significant amount of gas at $\sim 10^5$ K. The Toomre parameter always exceeds the critical value for gravito-turbulence, $Q\sim 1.5-2$, mainly due to the contribution of warm/hot gas inside the disc. Without feedback, a nearly gravito-turbulent regime establishes with similar gas velocity dispersion and lower $Q$. We propose that the "hot and turbulent" disc regime seen in our simulations, unlike the "cold and turbulent" gravito-turbulent regime of massive clumpy disc galaxies at $z\sim 1-2$, is a fundamental characterisation of main sequence galaxies at $z\gtrsim 6$, as they can sustain star formation rates comparable to those of low-mass starbursts at $z=0$. This results in no sustained coherent gas inflows through the disc, and in fluctuating and anisotropic mass transport, possibly postponing the assembly of the bulge and causing the initial feeding of the central black hole to be highly intermittent., 22 pages, 18 figures, 1 table, MNRAS accepted, mistakes in figures corrected, results unchanged

We present two-dimensional hydrodynamical simulations of pairs of planets migrating simultaneously in the Type I regime in a protoplanetary disc. Convergent migration naturally leads to the trapping of these planets in mean-motion resonances. Once in resonance the planets' eccentricity grows rapidly, and disc-planet torques cause the planets to escape resonance on a time-scale of a few hundred orbits. The effect is more pronounced in highly viscous discs, but operates efficiently even in inviscid discs. We attribute this resonance-breaking to overstable librations driven by moderate eccentricity damping, but find that this mechanism operates differently in hydrodynamic simulations than in previous analytic calculations. Planets escaping resonance in this manner can potentially explain the observed paucity of resonances in Kepler multi-transiting systems, and we suggest that simultaneous disc-driven migration remains the most plausible means of assembling tightly-packed planetary systems., Comment: MNRAS; Accepted 16/10/2017, 12 pages, 9 figures, 1 table

Along with air temperatures, the freezing level height (FLH) has risen over the last decades. The mass balance of tropical glaciers in Peru is highly sensitive to a rise in the FLH, mainly due to a decrease in accumulation and increase of energy for ablation caused by reduced albedo. Knowledge of future changes in the FLH is thus crucial to estimating changes in glacier extents. Since in situ data are scarce at altitudes where glaciers exist (above ~4800 m above sea level (asl)), reliable FLH estimates must be derived from multiple data types. Here we assessed the FLHs and their spatiotemporal variability, as well as the related snow/rain transition in the two largest glacier-covered regions in Peru by combining data from two climate reanalysis products, Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar Bright Band data, Micro Rain Radar data, and meteorological ground station measurements. The mean annual FLH lies at 4900 and 5010 m asl, for the Cordillera Blanca and Vilcanota, respectively. During the wet season, the FLH in the Cordillera Vilcanota lies ~150 m higher compared to the Cordillera Blanca, which is in line with the higher glacier terminus elevations. Coupled Model Intercomparison Project version 5 (CMIP5) climate model projections reveal that by the end of the 21st century, the FLH will rise by 230 m (±190 m) for Representative Concentration Pathway (RCP) 2.6 and 850 m (±390 m) for RCP8.5. Even under the most optimistic scenario, glaciers may continue shrinking considerably, assuming a close relation between FLH and glacier extents. Under the most pessimistic scenario, glaciers may only remain at the highest summits above approximately 5800 m asl. © 2017. American Geophysical Union. Por pares

$\Lambda$-Warm Dark Matter (WDM) has been proposed as alternative scenario to $\Lambda$ cold dark matter (CDM), motivated by discrepancies at the scale of dwarf galaxies, with less small-scale power and realized by collisionless particles with energies in the range $1-3$ keV. We present a new approach to constrain the viability of such WDM models using star formation histories of the dwarf spheroidal galaxies (dSphs) in the Local Group. We compare their high time-resolution star formation histories (SFHs) obtained with HST-based color magnitude diagrams with the range of possible collapse redshifts of their dark matter halos expected in CDM and in different WDM scenarios. The collapse redshift is inferred after determining a plausible infall mass of the subhalo. This is based on the current mass of individual dwarf inferred from stellar kinematics combined with results of cosmological simulations providing information on the subhalo evolution. Since WDM subhalos close to the filtering mass scale form significantly later than CDM, we show that they are in the first place difficult to reconcile with a truncation of star formation occurring as early as $z\geq 3$. The Ultra-Faint Dwarfs (UFDs) provide the most stringent constraints. Using 6 UFDs with the best determination of the SFHs, we show that we can exclude a 1 keV warm particle to a 2-$\sigma$ confidence interval consistently with other methods reported in the literature. For some objects the $2$ keV model is also excluded. We discuss the various caveats of the method, most notably the low number of dwarfs with accurately determined star formation histories and the uncertainties in the determination of the infall mass of the subhalos. Our preliminary analysis serves as a pathfinder for future investigations that will combine upcoming accurate SFHs for more local dSphs with direct analysis of WDM cosmological simulations with baryons., Comment: 13 pages, 9 figures, 5 tables, matches the version published in ApJ (discussion extended, 2 figures added to account for the caveats in the mass estimations, overall conclusions not changed)

We analyse stellar masses of clumps drawn from a compilation of star-forming galaxies at 1.1, Comment: 8 pages, 3 figures, ApJL accepted

We quantify the gas-phase abundance of deuterium and fractional contribution of stellar mass loss to the gas in cosmological zoom-in simulations from the Feedback In Realistic Environments project. At low metallicity, our simulations confirm that the deuterium abundance is very close to the primordial value. The chemical evolution of the deuterium abundance that we derive here agrees quantitatively with analytical chemical evolution models. We furthermore find that the relation between the deuterium and oxygen abundance exhibits very little scatter. We compare our simulations to existing high-redshift observations in order to determine a primordial deuterium fraction of 2.549 +/- 0.033 x 10^-5 and stress that future observations at higher metallicity can also be used to constrain this value. At fixed metallicity, the deuterium fraction decreases slightly with decreasing redshift, due to the increased importance of mass loss from intermediate-mass stars. We find that the evolution of the average deuterium fraction in a galaxy correlates with its star formation history. Our simulations are consistent with observations of the Milky Way's interstellar medium: the deuterium fraction at the solar circle is 85-92 per cent of the primordial deuterium fraction. We use our simulations to make predictions for future observations. In particular, the deuterium abundance is lower at smaller galactocentric radii and in higher mass galaxies, showing that stellar mass loss is more important for fuelling star formation in these regimes (and can even dominate). Gas accreting onto galaxies has a deuterium fraction above that of the galaxies' interstellar medium, but below the primordial fraction, because it is a mix of gas accreting from the intergalactic medium and gas previously ejected or stripped from galaxies., Comment: Accepted for publication in MNRAS. Revised version: expanded discussion and added Figure 2 (residual dependence on iron abundance)