Search

We present results from a high-resolution, cosmological, $\Lambda$CDM simulation of a group of field dwarf galaxies with the "superbubble" model for clustered SN feedback, accounting for thermal conduction and cold gas evaporation. The initial conditions and the galaxy formation physics, other than SN feedback, are the same as in Shen et al. (2014). The simulated luminous galaxies have blue colors, low star formation efficiencies and metallicities, and high cold gas content, reproducing the observed scaling relations of dwarfs in the Local Volume. Bursty star formation histories and superbubble-driven outflows lead to the formation of kpc-size DM cores when stellar masses reaches $M_{*} > 10^6$ $M_{\odot}$, similar to previous findings. However, the superbubble model appears more effective in destroying DM cusps than the previously adopted "blastwave" model, reflecting a higher coupling efficiency of SN energy with the ISM. On larger scale, superbubble-driven outflows have a more moderate impact: galaxies have higher gas content, more extended stellar disks, and a smaller metal-enriched region in the CGM. The two halos with $M_{vir} \sim 10^9$ $M_{\odot}$, which formed ultra-faint dwarf galaxies in Shen et al. (2014), remain dark due to the different impact of metal-enriched galactic winds from two nearby luminous galaxies. The column density distributions of H I, Si II, C IV and O VI are in agreement with recent observations of CGM around isolated dwarfs. While H I is ubiquitous with a covering fraction of unity within the CGM, Si II and C IV are less extended. O VI is more extended, but its mass is only 11% of the total CGM oxygen budget, as the diffuse CGM is highly ionised by the UVB. Superbubble feedback produces C IV and O VI an order of magnitude higher column densities than those with blastwave feedback. The CGM and DM cores are most sensitive probes of feedback mechanisms., Comment: 16 pages, 15 figures

Stars in globular clusters (GCs) lose a non negligible amount of mass during their post-main sequence evolution. This material is then expected to build up a substantial intra-cluster medium (ICM) within the GC. However, the observed gas content in GCs is a couple of orders of magnitude below these expectations. Here we follow the evolution of this stellar wind material through hydrodynamical simulations to attempt to reconcile theoretical predictions with observations. We test different mechanisms proposed in the literature to clear out the gas such as ram-pressure stripping by the motion of the GC in the Galactic halo medium and ionisation by UV sources. We use the code ramses to run 3D hydrodynamical simulations to study for the first time the ICM evolution within discretised multi-mass GC models including stellar winds and full radiative transfer. We find that the inclusion of both ram-pressure and ionisation is mandatory to explain why only a very low amount of ionised gas is observed in the core of GCs. The same mechanisms operating in ancient GCs that clear the gas could also be efficient at younger ages, meaning that young GCs would not be able to retain gas and form multiple generations of stars as assumed in many models to explain "multiple populations". However, this rapid clearing of gas is consistent with observations of young massive clusters., Comment: 11 pages, 5 figures, accepted for publication in MNRAS

Exoplanet characterization is one of the main foci of current exoplanetary science. For super-Earths and sub-Neptunes, we mostly rely on mass and radius measurements, which allow to derive the body's mean density and give a rough estimate of the planet's bulk composition. However, the determination of planetary interiors is a very challenging task. In addition to the uncertainty in the observed fundamental parameters, theoretical models are limited due to the degeneracy in determining the planetary composition. We aim to study several aspects that affect internal characterization of super-Earths and sub-Neptunes: observational uncertainties, location on the M-R diagram, impact of additional constraints as bulk abundances or irradiation, and model assumptions. We use a full probabilistic Bayesian inference analysis that accounts for observational and model uncertainties. We employ a Nested Sampling scheme to efficiently produce the posterior probability distributions for all the planetary structural parameter of interest. We include a structural model based on self-consistent thermodynamics of core, mantle, high-pressure ice, liquid water, and H-He envelope. Regarding the effect of mass and radius uncertainties on the determination of the internal structure, we find three different regimes: below the Earth-like composition line and above the pure-water composition line smaller observational uncertainties lead to better determination of the core and atmosphere mass respectively, and between them structure characterization only weakly depends on the observational uncertainties. We show that small variations in the temperature or entropy profiles lead to radius variations that are comparable to the observational uncertainty, suggesting that uncertainties linked to model assumptions can become more relevant to determine the internal structure than observational uncertainties., Comment: 12 pages, 12 figures

Using 324 numerically modelled galaxy clusters we investigate the radial and galaxy-halo alignment of dark matter subhaloes and satellite galaxies orbiting within and around them. We find that radial alignment depends on distance to the centre of the galaxy cluster but appears independent of the dynamical state of the central host cluster. Furthermore, we cannot find a relation between radial alignment of the halo or galaxy shape with its own mass. We report that backsplash galaxies, i.e. objects that have already passed through the cluster radius but are now located in the outskirts, show a stronger radial alignment than infalling objects. We further find that there exists a population of well radially aligned objects passing very close to the central cluster's centre which were found to be on highly radial orbit., 14 pages, accepted for publication in MNRAS

A promising route for revealing the existence of dark matter structures on mass scales smaller than the faintest galaxies is through their effect on strong gravitational lenses. We examine the role of local, lens-proximate clustering in boosting the lensing probability relative to contributions from substructure and unclustered line-of-sight (LOS) halos. Using two cosmological simulations that can resolve halo masses of $M_{\rm halo} \simeq 10^{9}\ M_{\odot}$ (in a simulation box of length $L_{\rm box}{\sim}100\,{\rm Mpc}$) and $10^{7}\ M_{\odot}$ ($L_{\rm box}\sim20\,{\rm Mpc}$), we demonstrate that clustering in the vicinity of the lens host produces a clear enhancement relative to an assumption of unclustered halos that persists to $> 20\,R_{\rm vir}$. This enhancement exceeds estimates that use a two-halo term to account for clustering, particularly within $2-5\,R_{\rm vir}$. We provide an analytic expression for this excess, clustered contribution. We find that local clustering boosts the expected count of $10^9 \ M_\odot$ perturbing halos by ${\sim}35\%$ compared to substructure alone, a result that will significantly enhance expected signals for low-redshift ($z_l \simeq 0.2$) lenses, where substructure contributes substantially compared to LOS halos. We also find that the orientation of the lens with respect to the line of sight (e.g., whether the line of sight passes through the major axis of the lens) can also have a significant effect on the lensing signal, boosting counts by an additional $\sim 50\%$ compared to a random orientations. This could be important if discovered lenses are biased to be oriented along their principal axis., Accepted to MNRAS

We present a numerical weak-lensing analysis that is fully relativistic and non-perturbative for the scalar part of the gravitational potential and first-order in the vector part, frame dragging. Integrating the photon geodesics backwards from the observer to the emitters, we solve the Sachs optical equations and study in detail the weak-lensing convergence, ellipticity and rotation. For the first time, we apply such an analysis to a high-resolution relativistic N-body simulation, which consistently includes the leading-order corrections due to general relativity on both large and small scales. These are related to the question of gauge choice and to post-Newtonian corrections, respectively. We present the angular power spectra and one-point probability distribution functions for the weak-lensing variables, which we find are broadly in agreement with comparable Newtonian simulations. Our geometric approach, however, is more robust and flexible, and can therefore be applied consistently to non-standard cosmologies and modified theories of gravity., Comment: 20 pages, 17 figures, corresponds to published version

Venus’ mass and radius are similar to those of Earth. However, dissimilarities in atmospheric properties, geophysical activity, and magnetic field generation could hint toward significant differences in the chemical composition and interior evolution of the two planets. Although various explanations for the differences between Venus and Earth have been proposed, the currently available data are insufficient to discriminate among the different solutions. Here we investigate the possible range of models for Venus’ structure. We assume that core segregation happened as a single-stage event. The mantle composition is inferred from the core composition using a prescription for metal-silicate partitioning. We consider three different cases for the composition of Venus defined via the bulk Si and Mg content, and the core’s S content. Permissible ranges for the core size, mantle, and core composition as well as the normalized moment of inertia (MoI) are presented for these compositions. A solid inner core could exist for all compositions. We estimate that Venus’ MoI is 0.317–0.351 and its core size 2930–4350 km for all assumed compositions. Higher MoI values correspond to more oxidizing conditions during core segregation. A determination of the abundance of FeO in Venus’ mantle by future missions could further constrain its composition and internal structure. This can reveal important information on Venus’ formation and evolution, and, possibly, the reasons for the differences between Venus and our home planet.

We hypothesise that differences in the temperatures at which the rocky material condensed out of the nebula gas can lead to differences in the composition of key rocky species (e.g., Fe, Mg, Si, Ca, Al, Na) and thus planet bulk density. Such differences in the observed bulk density of planets may occur as a function of radial location and time of planet formation. In this work we show that the predicted differences are on the cusp of being detectable with current instrumentation. In fact, for HD 219134, the 10 % lower bulk density of planet b compared to planet c could be explained by enhancements in Ca, Al rich minerals. However, we also show that the 11 % uncertainties on the individual bulk densities are not sufficiently accurate to exclude the absence of a density difference as well as differences in volatile layers. Besides HD 219134 b, we demonstrate that 55 Cnc e and WASP-47 e are similar candidates of a new Super-Earth class that have no core and are rich in Ca and Al minerals which are among the first solids that condense from a cooling proto-planetary disc. Planets of this class have densities 10-20% lower than Earth-like compositions and may have very different interior dynamics, outgassing histories and magnetic fields compared to the majority of Super-Earths., Comment: accepted for publication in MNRAS, 17 pages, 16 figures

Molecular clouds are turbulent structures whose star formation efficiency (SFE) is strongly affected by internal stellar feedback processes. In this paper we determine how sensitive the SFE of molecular clouds is to randomised inputs in the star formation feedback loop, and to what extent relationships between emergent cloud properties and the SFE can be recovered. We introduce the yule suite of 26 radiative magnetohydrodynamic (RMHD) simulations of a 10,000 solar mass cloud similar to those in the solar neighbourhood. We use the same initial global properties in every simulation but vary the initial mass function (IMF) sampling and initial cloud velocity structure. The final SFE lies between 6 and 23 percent when either of these parameters are changed. We use Bayesian mixed-effects models to uncover trends in the SFE. The number of photons emitted early in the cluster's life and the length of the cloud provide are the strongest predictors of the SFE. The HII regions evolve following an analytic model of expansion into a roughly isothermal density field. The more efficient feedback is at evaporating the cloud, the less the star cluster is dispersed. We argue that this is because if the gas is evaporated slowly, the stars are dragged outwards towards surviving gas clumps due to the gravitational attraction between the stars and gas. While star formation and feedback efficiencies are dependent on nonlinear processes, statistical models describing cloud-scale processes can be constructed., 24 pages, 16 figures, 6 tables. Accepted to MNRAS, version updated with published title

The first visibility of the lunar crescent signals the start of a new month in the Islamic calendar. The eminent astronomer Ḥabash al-Ḥāsib sib developed a method of uncompromising complexity for predicting the visibility of the lunar crescent. He derived his threshold value from a moonwatch carried out at different places in Iraq on November 17th, 860 CE. We will allude to a few modern visibility criteria as well and highlight the uncertainties of today’s calculations when converting historical Arabic into Julian or Gregorian dates. Tables of first visibility of the lunar crescent for different locations are provided for the purpose of date conversions. Since the Islamic calendar was based on the observation of the lunar crescent, historical dates imply information on positive and negative sightings of the lunar crescent. From such information estimations of cloudiness in different regions of the Islamic world can be extracted.

It has been proposed that close interactions with satellite galaxies can significantly perturb the morphology of the main galaxy. However, the dynamics of an already formed bar following the interaction with the external environment has not been studied in detail in a fully cosmological context. In this work, analysing the cosmological zoom-in simulation Eris2k, we study the effects that a very unequal-mass flyby crossing the stellar disc has on the stability of the pre-existing bar. We characterize the evolution of the bar strength and length showing that the perturbation exerted by the flyby shuffles the orbits of stars for less than one Gyr. After this time, the bar shows a remarkable resilience, reforming with properties comparable to those it had before the interaction. Our work shows that close unequal-mass encounters, the most frequent interactions occurring during the evolution of cosmic structures, have (i) an overall minor impact on the global evolution of the bar in the long term, still (ii) the effect is destructive and (iii) a very weak interaction is sufficient to dismantle a strong bar leading to its "apparent death". As a consequence, due to the non-negligible duration of the bar-less period, a fraction of observed spiral galaxies classified as non-barred could be prone to bar formation., 6 pages, 3 Figure; Accepted for publication in MNRAS

Journal of Geophysical Research: Earth Surface, 123 (6), ISSN:0148-0227, ISSN:2169-9003, ISSN:2169-9011

The dearth of dwarf galaxies in the local universe is hard to reconcile with the large number of low mass haloes expected within the concordance $\Lambda$CDM paradigm. In this paper we perform a systematic evaluation of the uncertainties affecting the measurement of DM halo abundance using galaxy kinematics. Using a large sample of dwarf galaxies with spatially resolved kinematic data we derive a correction to obtain the observed abundance of galaxies as a function of their halo maximum circular velocity from the line-of-sight velocity function in the Local Volume. This estimate provides a direct means of comparing the predictions of theoretical models and simulations (including nonstandard cosmologies and novel galaxy formation physics) to the observational constraints. The new "galactic $V_{max}$" function is steeper than the line-of-sight velocity function but still shallower than the theoretical CDM expectation, showing that some unaccounted physical process is necessary to reduce the abundance of galaxies and/or drastically modify their density profiles compared to CDM haloes. Using this new galactic $V_{max}$ function, we investigate the viability of baryonic solutions such as feedback-powered outflows and photoevaporation of gas from an ionising radiation background. At the 3-$\sigma$ confidence level neither energetic feedback nor photoevaporation are effective enough to reconcile the disagreement. In the case of maximum baryonic effects, the theoretical estimate still deviates significantly from the observations for $V_{max} < 60$ km/s. CDM predicts at least 1.8 times more galaxies with $V_{max} = 50$ km/s and 2.5 times more than observed at $30$ km/s. Recent hydrodynamic simulations seem to resolve the discrepancy but disagree with the properties of observed galaxies with resolved kinematics. (abridged), Comment: 17 pages, 22 figures; major revisions include clarification of the method, expanded comparison with simulations with a new figure, analysis of uncertainties in model as well as pressure support corrections, and a new table with nomenclature

Cosmic far-infrared background (CFIRB) is a powerful probe of the history of star formation rate (SFR) and the connection between baryons and dark matter across cosmic time. In this work, we explore to which extent the CFIRB anisotropies can be reproduced by a simple physical framework for galaxy evolution, the gas regulator (bathtub) model. This model is based on continuity equations for gas, stars, and metals, taking into account cosmic gas accretion, star formation, and gas ejection. We model the large-scale galaxy bias and small-scale shot noise self-consistently, and we constrain our model using the CFIRB power spectra measured by Planck. Because of the simplicity of the physical model, the goodness of fit is limited. We compare our model predictions with the observed correlation between CFIRB and gravitational lensing, bolometric infrared luminosity functions, and submillimetre source counts. The strong clustering of CFIRB indicates a large galaxy bias, which corresponds to haloes of mass $10^{12.5} M_\odot$ at z=2, higher than the mass associated with the peak of the star formation efficiency. We also find that the far-infrared luminosities of haloes above $10^{12} M_\odot$ are higher than the expectation from the SFR observed in ultraviolet and optical surveys., Comment: 23 pages, 21 figures, 4 tables; replaced to match published version

Context. The formation of the most massive quasars observed at high redshifts requires extreme inflows of gas down to the length scales of the central compact object. Aims. Here we estimate the maximum inflow rate allowed by gravity down to the surface of supermassive stars, the possible progenitors of these supermassive black holes. Methods. We use the continuity equation and the assumption of spherical symmetry and free fall to derive the maximum allowed inflow rates for various density profiles. We apply our approach to the mass–radius relation of rapidly accreting supermassive stars to estimate an upper limit to the accretion rates allowed during the formation of these objects. Results. We find that, as long as the density of the accreted gas is smaller than or equal to the average density of the accretor, the maximum allowed rate, Ṁmax, is given uniquely by the compactness of the accretor. We argue that a density inversion between accreting matter and the accretor is inconsistent with gravitational collapse. For the compactness of rapidly accreting supermassive stars, Ṁmax is related to the stellar mass, M, by a power law, Ṁmax ∝ M3/4. The rates of atomically cooled halos (0.1−10 M⊙ yr−1) are allowed as soon as M ≳ 1 M⊙. The largest rates expected in galaxy mergers (104 − 105 M⊙ yr−1) become accessible once the accretor is supermassive (M ≳ 104 M⊙). Conclusions. These results suggest that supermassive stars can accrete up to masses > 106 M⊙ before they collapse via the general-relativistic instability. At such masses, the collapse is expected to lead to the direct formation of a supermassive black hole, even within metal-rich gas, resulting in a black hole seed that is significantly heavier than in conventional direct collapse models for atomic cooling halos.

Although skinfold-derived equations seem to be practical for field application in estimating body fat percentage (BF%) and minimum body mass in Olympic wrestlers, prediction equations applied first need to be cross-validated in Olympic wrestlers to define the best prediction equation. This study aimed to evaluate the most accurate field method to predict BF% in Olympic wrestlers compared to BF% estimated by air displacement plethysmography (ADP). Sixty-one male (body mass 72.4 ± 13.5 kg, height 170.3 ± 7.0 cm, body mass index (BMI) 24.9 ± 3.5 kg.m−2, BF% 8.5 ± 4.9%) and twenty-five female wrestlers (body mass 60.3 ± 9.9 kg, height 161.3 ± 7.1 cm, BMI 23.1 ± 2.5 kg.m−2, BF% 18.7 ± 4.7%) undertook body composition assessments including ADP and nine-site skinfold measurements. Correlations, bias, limits of agreement, and standardized differences between alterations in BF% measured by ADP and other prediction equations were evaluated to validate measures, and multiple regression analyses to develop an Olympic wrestlers-specific prediction formula. The Stewart and Hannan equation for male wrestlers and the Durnin and Womersley equation for female wrestlers provided the most accurate BF% compared to the measured BF% by ADP, with the lowest bias and presented no significant differences between the measured and predicted BF%. A new prediction equation was developed using only abdominal skinfold and sex as variables, predicting 83.2% of the variance. The findings suggest the use of the new wrestler-specific prediction equation proposed in the study as a valid and accurate alternative to ADP to quantify BF% among Olympic wrestlers.

We present mock observations of forming planets with ALMA. The possible detections of circumplanetary disks (CPDs) were investigated around planets of Saturn, 1, 3, 5, and 10 Jupiter-masses that are placed at 5.2 AU from their star. The radiative, three dimensional hydrodynamic simulations were then post-processed with RADMC3D and the ALMA Observation Simulator. We found that even though the CPDs are too small to be resolved, they are hot due to the accreting planet in the optically thick limit, therefore the best chance to detect them with continuum observations in this case is at the shortest ALMA wavelengths, such as Band 9 (440 microns). Similar fluxes were found in the case of Saturn and Jupiter-mass planets, as for the 10 $\mathrm{M_{Jup}}$ gas-giant, due to temperature weighted optical depth effects: when no deep gap is carved, the planet region is blanketed by the optically thick circumstellar disk leading to a less efficient cooling there. A test was made for a 52 AU orbital separation, showed that optically thin CPDs are also detectable in band 7 but they need longer integration times ($>$5hrs). Comparing the gap profiles of the same simulation at various ALMA bands and the hydro simulation confirmed that they change significantly, first because the gap is wider at longer wavelengths due to decreasing optical depth; second, the beam convolution makes the gap shallower and at least 25% narrower. Therefore, caution has to be made when estimating planet masses based on ALMA continuum observations of gaps., Accepted for publication at MNRAS. Typos are corrected since previous version. 11 pages, 5 tables, 4 figures

Most stars in the Galaxy are believed to be formed within star clusters from collapsing molecular clouds. However, the complete process of star formation, from the parent cloud to a gas-free star cluster, is still poorly understood. We perform radiation-hydrodynamical simulations of the collapse of a turbulent molecular cloud using the RAMSES-RT code. Stars are modelled using sink particles, from which we self-consistently follow the propagation of the ionising radiation. We study how different feedback models affect the gas expulsion from the cloud and how they shape the final properties of the emerging star cluster. We find that the star formation efficiency is lower for stronger feedback models. Feedback also changes the high mass end of the stellar mass function. Stronger feedback also allows the establishment of a lower density star cluster, which can maintain a virial or sub-virial state. In the absence of feedback, the star formation efficiency is very high, as well as the final stellar density. As a result, high energy close encounters make the cluster evaporate quickly. Other indicators, such as mass segregation, statistics of multiple systems and escaping stars confirm this picture. Observations of young star clusters are in best agreement with our strong feedback simulation., Comment: 20 pages, 14 figures. Accepted by MNRAS

Chondrules are crystallised droplets of silicates formed by rapid heating to high temperatures (> 1800 K) of solid precursors followed by hours or days of cooling. Dating of chondrules is consistent with the formation timescale of Jupiter in the core-accretion model (1-4 Myrs). Here we investigate if the shocks generated by a massive planet could generate flash heating episodes necessary to form chondrules using high resolution 2D simulations with the multi-fluid code ROSSBI. We use different radiative cooling prescriptions, i.e. different cooling rates and models, and vary planet mass, orbit and disk models. For long disk cooling rates (> 1000 orbits) and a massive protoplanet (> 0.75 jupiter masses), we obtain hot enough temperatures for chondrule formation, while using more realistic thermodynamics is not successful in the Minimum Mass Solar Nebula (MMSN) model. However, sufficient flash heating is achieved by a Jupiter mass planet in a 5 times more massive disk, which is a conceivable scenario for the young solar nebula and exoplanetary systems. Interestingly, a gap-forming massive planet triggers vortices, which can trap dust, i.e. chondrule precursors, and generates a high pressure environment that is consistent with cosmochemical evidence from chondrules. A massive gas giant can thus, in principle, both stimulate the concentration of chondrule precursors in the vicinity of the shocking regions, and flash-heat them via the shocks., 12 pages, 7 figures, Submitted to ApJ

Adequate sleep quantity and quality is required to maintain vigilance, cognitive and learning processes. A decrease of sleep quantity preflight and on the International Space Station (ISS) has been reported. Recent counter-measures have been implemented to better regulate sleep opportunities on ISS. In our study, astronauts were allocated enough time for sleep the night before the recordings. However, for proper sleep recovery, the quality of sleep is also critical. Unfortunately, data on sleep quality have yet to be acquired from the ISS. Here, we investigate sleep pressure markers during wakefulness in five astronauts throughout their 6-month space mission by the mean of electroencephalographic recordings. We show a global increase of theta oscillations (5–7 Hz) on the ISS compared to on Earth before the mission. We also show that local sleep-like events, another marker of sleep pressure, are more global in space (p < 0.001). By analysing the performances of the astronauts during a docking simulation, we found that local sleep-like events are more global when reaction times are slower (R 2 = 0.03, p = 0.006) and there is an increase of reaction times above 244 ms after 2 months in space (p = 0.012). Our analyses provide first evidence for increased sleep pressure in space and raise awareness on possible impacts on visuomotor performances in space., SCOPUS: ar.j, info:eu-repo/semantics/published

Rock glaciers are creeping periglacial landforms experiencing strong acceleration during recent atmospheric warming and raising concerns with regard to their future behaviour and stability. High resolution kinematic observations show strong seasonal and multi-annual variations in rock glacier creep, but the linking mechanisms to environmental forcing remain poorly understood and lack quantitative models. Here we investigate the interaction between rock glacier creep and climatic forcing - temperature and precipitation - by developing a novel conceptual and numerical modelling approach. The model is constrained and the results are compared with data from the Dirru Rock Glacier (Vallis - CH). We are able to reproduce the observed velocity variations both in magnitude and phase on seasonal and inter-annual time scales. We find that water from liquid precipitation and snow melt, rather than air temperature, is the main driver of variations in rock glacier creep. Our results imply that the influence of water on rock glacier creep is fundamental and must be considered when investigating the historic and future evolution of rock glaciers.

We present a revised version of Peters' (1964) time-scale for the gravitational-wave (GW) induced decay of two point masses. The new formula includes the effects of the first-order post-Newtonian perturbation and additionally provides a simple fit to account for the Newtonian self-consistent evolution of the eccentricity. The revised time-scale is found by multiplying Peters' estimate by two factors, $R(e_0)= 8^{1-\sqrt{1-e_0}}$ and $Q_{\rm f}(p_0) = \exp \left(2.5 (r_{\rm S}/p_0) \right)$, where $e_0$ and $p_0$ are the initial eccentricity and periapsis, respectively, and $r_{\rm S}$ the Schwarzschild radius of the system. Their use can correct errors of a factor of 1-10 that arise from using the original Peters' formula. We apply the revised time-scales to a set of typical sources for existing ground-based laser interferometers and for the future Laser Interferometer Space Antenna (LISA), at the onset of their GW driven decay. We argue that our more accurate model for the orbital evolution will affect current event- and detection-rate estimates for mergers of compact object binaries, with stronger deviations for eccentric LISA sources, such as extreme and intermediate mass-ratio inspirals. We propose the correction factors $R$ and $Q_{\rm f}$ as a simple prescription to quantify decay time-scales more accurately in future population synthesis models. We also suggest that the corrected time-scale may be used as a computationally efficient alternative to numerical integration in other applications that include the modelling of radiation reaction for eccentric sources., Accepted for publication in MNRAS

We study star formation histories (SFHs) of $\simeq500$ dwarf galaxies (stellar mass $M_\ast = 10^5 - 10^9\,M_\odot$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher-mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of "satellite vs. central" and "LG vs. individual MWvs. isolated dwarf central." Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the "near field" have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated ("true field") dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for $M_\ast = 10^6 - 10^7\,M_\odot$ but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG., Main text: 11 pages, 8 figures; appendices: 4 pages, 4 figures. Submitted to MNRAS; comments welcome

We simulate the flux emitted from galaxy halos in order to quantify the brightness of the circumgalactic medium (CGM). We use dedicated zoom-in cosmological simulations with the hydrodynamical Adaptive Mesh Refinement code RAMSES, which are evolved down to z=0 and reach a maximum spatial resolution of 380 $h^{-1}$pc and a gas mass resolution up to 1.8$\times 10^{5} h^{-1} \rm{M}_{\odot}$ in the densest regions. We compute the expected emission from the gas in the CGM using CLOUDY emissivity models for different lines (e.g. Ly$\alpha$, CIV, OVI, CVI, OVIII) considering UV background fluorescence, gravitational cooling and continuum emission. In the case of Ly$\alpha$ we additionally consider the scattering of continuum photons. We compare our predictions to current observations and find them to be in good agreement at any redshift after adjusting the Ly$\alpha$ escape fraction. We combine our mock observations with instrument models for FIREBall-2 (UV balloon spectrograph) and HARMONI (visible and NIR IFU on the ELT) to predict CGM observations with either instrument and optimise target selections and observing strategies. Our results show that Ly$\alpha$ emission from the CGM at a redshift of 0.7 will be observable with FIREBall-2 for bright galaxies (NUV$\sim$18 mag), while metal lines like OVI and CIV will remain challenging to detect. HARMONI is found to be well suited to study the CGM at different redshifts with various tracers., Comment: accepted for publication in MNRAS

Soon after the recent first ever detection of gravitational waves from merging black holes it has been suggested that their origin is primordial. Appealingly, a sufficient number of primordial black holes (PBHs) could also partially or entirely constitute the dark matter (DM) in our Universe. However, recent studies on PBHs in ultra-faint dwarf galaxies (UFDGs) suggest that they would dynamically heat up the stellar component due to two-body relaxation processes. From the comparison with the observed stellar velocity dispersions and the stellar half-light radii, it was claimed that only PBHs with masses $\lesssim 10\, {\rm M}_\odot$ can significantly contribute to the DM. In this work, we improve the latter constraints by considering the largest observational sample of UFDGs and by allowing the PBH masses to follow an extended (lognormal) distribution. By means of collisional Fokker–Planck simulations, we explore a wide parameter space of UFDGs containing PBHs. The analysis of the half-light radii and velocity dispersions resulting from the simulations leads to three general findings that exclude PBHs with masses $\sim \mathcal {O}(1\operatorname{-}100)\, {\rm M}_\odot {}$ from constituting all of the DM: (i) we identify a critical sub-sample of UFDGs that only allows for $\sim \mathcal {O}(1)\, {\rm M}_\odot$ PBH masses; (ii) for any PBH mass, there is an UFDG in our sample that disfavours it; (iii) the spatial extensions of a majority of simulated UFDGs containing PBHs are too large to match the observed.

The harvest of exoplanet discoveries has opened the area of exoplanet characterisation. But this cannot be achieved without a careful analysis of the host star parameters. The system of HD219134 hosts two transiting exoplanets and at least two additional non-transiting exoplanets. We used the VEGA/CHARA interferometer to measure the angular diameter of HD219134, leading to a stellar radius of $R_{\star}=0.726\pm0.014 R_{\odot}$. We also derived the stellar density from the transits light curves ($\rho_{\star}=1.82\pm0.19 \rho_{\odot}$), which finally gives a direct estimate of the mass ($M_{\star}=0.696\pm0.078 M_{\odot}$) with a correlation of 0.46 between $R_{\star}$ and $M_{\star}$. This new mass is smaller than that derived from the C2kSMO stellar evolutionary model, which provides a mass range of 0.755$-$0.810 ($\pm 0.040$) $M_{\odot}$. This allows us to infer the mass, radius and density of the two transiting exoplanets of the system. We then use an inference model to obtain the internal parameters of these two transiting exoplanets. Moreover, we find that planet $b$ and $c$ have smaller radii than previously estimated ($1.500\pm0.057$ and $1.415\pm0.049 R_{\oplus}$, respectively); this clearly puts these planets out of the gap in the exoplanetary radii distribution and validates their super-Earth nature. Planet $b$ is more massive than planet $c$, but possibly less dense. We investigate whether this could be caused by partial melting of the mantle and find that tidal heating due to non-zero eccentricity of planet $b$ may be powerful enough. The system of HD219134 constitutes a very valuable benchmark for both stellar physics and exoplanetary science. The direct determination of the stellar density, radius and mass should be more extensively applied to provide accurate exoplanets properties and calibrate stellar models., Comment: 12 pages, 6 figures, 7 tables; published in A$\&$A. (This version includes language editing corrections.)

We present the first detailed study of the spatially-resolved dust continuum emission of simulated galaxies at 1, Comment: 12 pages, 8 figures; accepted for publication in MNRAS

Galactic halos are of great importance for our understanding of both the dark matter nature and primordial non-Gaussianity in the perturbation spectrum, a powerful discriminant of the physical mechanisms that generated the cosmological fluctuations observed today. In this paper we analyze {\it Planck} data towards the galaxy M104 (Sombrero) and find an asymmetry in the microwave temperature which extends up to about $1 \degr$ from the galactic center. This frequency-independent asymmetry is consistent with that induced by the Doppler effect due to the galactic rotation and we find a probability of less than about $0.2\%$ that it is due to a random fluctuation of the microwave background. In addition, {\it Planck} data indicate the relatively complex dynamics of the M104 galactic halo, and this appears to be in agreement with previous studies. In view of our previous analysis of the dark halos of nearby galaxies, this finding confirms the efficiency of the method used in revealing and mapping the dark halos around relatively nearby edge-on galaxies., Comment: 6 pages, 5 Figures, in press on Astronomy and Astrophyiscs

Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disk galaxies, in which gas rotation accurately traces the gravitational potential, the interstellar medium (ISM) of z>1 galaxies is typically more dynamic and exhibits elevated turbulence. If not properly modeled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyze the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a range of galaxy properties from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M_* > 10^11 M_sun at z=1). Selecting only snapshots where a disk is present, we calculate the rotational profile v_phi(r) of the cool (10^3.5 K < T < 10^4.5 K) gas and compare it to the circular velocity v_c=sqrt(GM/r). In the simulated galaxies, the gas rotation traces the circular velocity at intermediate radii, but the two quantities diverge significantly in the center and in the outer disk. Our simulations appear to over-predict observed rotational velocities in the centers of massive galaxies (likely from a lack of black hole feedback), so we focus on larger radii. Gradients in the turbulent pressure at these radii can provide additional radial support and bias dynamical mass measurements low by up to 40%. In both the interior and exterior, the gas' motion can be significantly non-circular due to e.g. bars, satellites, and inflows/outflows. We discuss the accuracy of commonly-used analytic models for pressure gradients (or "asymmetric drift") in the ISM of high-redshift galaxies., 16 pages, 9 figures, MNRAS accepted version

Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimeter broadband flux to total infrared luminosity (L_IR), and hence star formation rate, in high-z galaxies. However, different definitions of dust temperatures have been used in the literature, leading to different physical interpretations of how ISM conditions change with, e.g., redshift and star formation rate. In this paper, we analyse the dust temperatures of massive (M* > 10^10 Msun) z=2-6 galaxies with the help of high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. At z~2, our simulations successfully predict dust temperatures in good agreement with observations. We find that dust temperatures based on the peak emission wavelength increase with redshift, in line with the higher star formation activity at higher redshift, and are strongly correlated with the specific star formation rate. In contrast, the mass-weighted dust temperature does not strongly evolve with redshift over z=2-6 at fixed IR luminosity but is tightly correlated with L_IR at fixed z. The mass-weighted temperature is important for accurately estimating the total dust mass. We also analyse an 'equivalent' dust temperature for converting (sub)millimeter flux density to total IR luminosity, and provide a fitting formula as a function of redshift and dust-to-metal ratio. We find that galaxies of higher equivalent (or higher peak) dust temperature ('warmer dust') do not necessarily have higher mass-weighted temperatures. A 'two-phase' picture for interstellar dust can explain the different scaling relations of the various dust temperatures., Comment: 26 pages, 15 figures, accepted for publication in MNRAS

Bars are a key factor in the long-term evolution of spiral galaxies, in their unique role in redistributing angular momentum and transporting gas and stars on large scales. The Eris-suite simulations are cosmological zoom-in, N-body, smoothed-particle hydrodynamic simulations built to follow the formation and evolution of a Milky Way-sized galaxy across the build-up of the large scale structure. Here we analyse and describe the outcome of two particular simulations taken from the Eris suite - ErisBH and Eris2k - which mainly differ in the prescriptions employed for gas cooling, star formation, and feedback from supernovae and black holes. Our study shows that the enhanced effective feedback in Eris2k, due to the collective effect of the different micro-physics implementations, results in a galaxy which is less massive than its ErisBH counterpart till z~1. However, when the stellar content is large enough so that global dynamical instabilities can be triggered, the galaxy in Eris2k develops a stronger and more extended bar with respect to ErisBH. We demonstrate that he structural properties and time evolution of the two bars are very different. Our results highlight the importance of accurate sub-grid prescriptions in cosmological zoom-in simulations of the process of galaxy formation and evolution, and the possible use of a statistical sample of barred galaxies to assess the strength of the stellar feedback., Comment: 16 pages, 14 Figure; accepted for publication in MNRAS; added references and figures, more details, results unchanged

Recently, seven Earth-sized planets were discovered around the M-dwarf star TRAPPIST-1. Thanks to transit-timing variations, the masses and therefore the bulk densities of the planets have been constrained, suggesting that all TRAPPIST-1 planets are consistent with water mass fractions on the order of 10%. These water fractions, as well as the similar planet masses within the system, constitute strong constraints on the origins of the TRAPPIST-1 system. In a previous work, we outlined a pebble-driven formation scenario. In this paper we investigate this formation scenario in more detail. We used a Lagrangian smooth-particle method to model the growth and drift of pebbles and the conversion of pebbles to planetesimals through the streaming instability. We used the N-body code \texttt{MERCURY} to follow the composition of planetesimals as they grow into protoplanets by merging and accreting pebbles. This code is adapted to account for pebble accretion, type-I migration, and gas drag. In this way, we modelled the entire planet formation process (pertaining to planet masses and compositions, not dynamical configuration). We find that planetesimals form in a single, early phase of streaming instability. The initially narrow annulus of planetesimals outside the snowline quickly broadens due to scattering. Our simulation results confirm that this formation pathway indeed leads to similarly-sized planets and is highly efficient in turning pebbles into planets ($\sim$50% solids-to-planets conversion efficiency). [...] The water content of planets resulting from our simulations is on the order of 10%, and our results predict a `V-shaped' trend in the planet water fraction with orbital distance: from relatively high (innermost planets) to relatively low (intermediate planets) to relatively high (outermost planets)., Comment: 16 pages, accepted for publication in A&A, abstract abridged

Natural floodplains are characterized by a complex habitat mosaic. However, damming, water storage, and hydropower production affect many floodplains by altering their natural habitat diversity. Field sampling data and imaging spectroscopy are used in combination with statistical models to assess resource allocation strategies of willow stands in perialpine floodplains. Three contrasting floodplain reaches located along two rivers in Switzerland serve as test beds: The Sarine River is partitioned into an upstream and downstream segment under the influence of a dam and a hydropower plant, while the Sense River represents an undisturbed, natural floodplain. Airborne imaging spectrometer data allow mapping of spatially distributed Competitor/Stress tolerator/Ruderal (CSR) strategies using a partial least square modeling approach. Using cross validation, we demonstrate that a statistical modeling approach can reveal variations in CSR scores based on the StrateFy model. Such intraspecific variation of CSR scores cannot be captured by a strategy categorization based solely on the species. Results reveal that willows shifted toward more competition and less stress tolerance along hydrologically altered reaches compared to the willows strategy along the natural control. Moreover, the overall distribution of strategies indicates that stress factors (i.e., limiting growth factors), rather than disturbance (i.e., events leading to partial or total destruction), shape the plant traits of alluvial willow trees. Detailed assessments of resource allocation strategies contribute to a more complete understanding of the continuous and reciprocal shaping between flow regimes, landforms, and alluvial vegetation.

Regulating the available gas mass inside galaxies proceeds through a delicate balance between inflows and outflows, but also through the internal depletion of gas due to star formation. At the same time, stellar feedback is the internal engine that powers the strong outflows. Since star formation and stellar feedback are both small scale phenomena, we need a realistic and predictive subgrid model for both. We describe the implementation of supernova momentum feedback and star formation based on the turbulence of the gas in the RAMSES code. For star formation, we adopt the so-called multi-freefall model. The resulting star formation efficiencies can be significantly smaller or bigger than the traditionally chosen value of $1\%$. We apply these new numerical models to a prototype cosmological simulation of a massive halo that features a major merger which results in the formation of an early-type galaxy without using AGN feedback. We find that the feedback model provides the first order mechanism for regulating the stellar and baryonic content in our simulated galaxy. At high redshift, the merger event pushes gas to large densities and large turbulent velocity dispersions, such that efficiencies come close to $10\%$, resulting in large SFR. We find small molecular gas depletion time during the starburst, in perfect agreement with observations. Furthermore, at late times, the galaxy becomes quiescent with efficiencies significantly smaller than $1\%$, resulting in small SFR and long molecular gas depletion time., 15 pages, 11 figures

It has been a common assumption of interior models that the outer planets of our solar system are convective, and that the internal temperature distributions are therefore adiabatic. This assumption is also often applied to exoplanets. However, if a large portion of the thermal flux can be transferred by conduction, or if convection is inhibited, the thermal profile could be substantially different and would therefore affect the inferred planetary composition. Here we investigate how the assumption of non-adiabatic temperature profiles in Uranus and Neptune affects their internal structures and compositions. We use a set of plausible temperature profiles together with density profiles that match the measured gravitational fields to derive the planets' compositions. We find that the inferred compositions of both Uranus and Neptune are quite sensitive to the assumed thermal profile in the outer layers, but relatively insensitive to the thermal profile in the central, high pressure region. The overall value of the heavy element mass fraction, $Z$, for these planets is between 0.8 and 0.9. Finally, we suggest that large parts of Uranus' interior might be conductive, a conclusion that is consistent with Uranus dynamo models and a hot central inner region., 40 pages, 5 tables, 5 figures, to appear in MNRAS

HESS J1825-137 is a bright very-high-energy (VHE) gamma-ray source that has been firmly established as a pulsar wind nebula (PWN), and one of the most extended gamma-ray objects within this category. The progenitor supernova remnant (SNR) for this PWN has not been firmly established. We carried out an analysis of gamma-ray observations in the region of HESS J1825-137 with the Fermi-LAT which reveal emission in the direction away from the Galactic plane. The region lies beyond the PWN and reaches a distance from the pulsar compatible with the supposed location of the SNR shock front. The spectrum of the gamma-rays is hard with a photon index of $\sim 1.9$ in the 10-250 GeV range. Several scenarios for the origin of the emission are discussed, including the SNR as a source of high-energy particles and the "leakage"' of leptons from the PWN., Comment: To be published in MNRAS