Testing hydrodynamics schemes in galaxy disc simulations

We examine how three fundamentally different numerical hydrodynamics codes follow the evolution of an isothermal galactic disc with an external spiral potential. We compare an adaptive mesh refinement code (RAMSES), a smoothed particle hydrodynamics code (SPHNG), and a volume-discretised meshless code (GIZMO). Using standard refinement criteria, we find that RAMSES produces a disc that is less vertically concentrated and does not reach such high densities as the SPHNG or gizmo runs. The gas surface density in the spiral arms increases at a lower rate for the RAMSES simulations compared to the other codes. There is also a greater degree of substructure in the SPHNG and GIZMOruns and secondary spiral arms are more pronounced. By resolving the Jeans’ length with a greater number of grid cells we achieve more similar results to the Lagrangian codes used in this study. Other alterations to the refinement scheme (adding extra levels of refinement and refining based on local density gradients) are less successful in reducing the disparity between RAMSES and SPHNG/GIZMO. Although more similar, SPHNG displays different density distributions and vertical mass profiles to all modes of gizmo (including the smoothed particle hydrodynamics version). This suggests differences also arise which are not intrinsic to the particular method but rather due to its implementation. The discrepancies between codes (in particular, the densities reached in the spiral arms) could potentially result in differences in the locations and timescales for gravitational collapse, and therefore impact star formation activity in more complex galaxy disc simulations

Feedback from massive stars and gas expulsion from proto-globular clusters

© 2015. The American Astronomical Society. All rights reserved. Globular clusters (GCs) are considerably more complex structures than previously thought, harboring at least two stellar generations that present clearly distinct chemical abundances. Scenarios explaining the abundance patterns in GCs mostly assume that originally the clusters had to be much more massive than today, and that the second generation of stars originates from the gas shed by stars of the first generation (FG). The lack of metallicity spread in most GCs further requires that the supernova-enriched gas ejected by the FG is completely lost within ∼30 Myr, a hypothesis never tested by means of three-dimensional hydrodynamic simulations. In this paper, we use 3D hydrodynamic simulations including stellar feedback from winds and supernovae, radiative cooling and self-gravity to study whether a realistic distribution of OB associations in a massive proto-GC of initial mass M tot ∼ 10 7 M o is sufficient to expel its entire gas content. Our numerical experiment shows that the coherence of different associations plays a fundamental role: as the bubbles interact, distort, and merge, they carve narrow tunnels that reach deeper and deeper toward the innermost cluster regions, and through which the gas is able to escape. Our results indicate that after 3 Myr, the feedback from stellar winds is responsible for the removal of ∼40% of the pristine gas, and that after 14 Myr, 99% of the initial gas mass has been removed

Properties of simulated Milky Way-mass galaxies in loose group and field environments

We test the validity of comparing simulated field disk galaxies with the empirical properties of systems situated within environments more comparable to loose groups, including the Milky Way's Local Group. Cosmological simulations of Milky Way-mass galaxies have been realised in two different environment samples: in the field and in environments with similar properties to the Local Group. Apart from the environments of the galaxies, the samples are kept as homogeneous as possible with equivalent ranges in last major merger time, halo mass and halo spin. Comparison of these two samples allow for systematic differences in the simulations to be identified. Metallicity gradients, disk scale lengths, colours, magnitudes and age-velocity dispersion relations are studied for each galaxy in the suite and the strength of the link between these and environment of the galaxies is studied. The bulge-to-disk ratio of the galaxies show that these galaxies are less spheroid dominated than many other simulated galaxies in literature with the majority of both samples being disk dominated. We find that secular evolution and mergers dominate the spread of morphologies and metallicity gradients with no visible differences between the two environment samples. In contrast with this consistency in the two samples there is tentative evidence for a systematic difference in the velocity dispersion-age relations of galaxies in the different environments. Loose group galaxies appear to have more discrete steps in their velocity dispersion-age relations. We conclude that at the current resolution of cosmological galaxy simulations field environment galaxies are sufficiently similar to those in loose groups to be acceptable proxies for comparison with the Milky Way provided that a similar assembly history is considered.Comment: 16 pages, 11 figures, abstract abridged for arXiv. Accepted for publication in Astronomy & Astrophysic

The role of stellar radial motions in shaping galaxy surface brightness profiles

Aims. The physics driving features such as breaks observed in galaxy surface brightness (SB) profiles remains contentious. Here, we assess the importance of stellar radial motions in shaping their characteristics. Methods. We use the simulated Milky Way-mass cosmological discs from the Ramses Disc Environment Study (RaDES) to characterise the radial redistribution of stars in galaxies displaying type-I (pure exponentials), II (downbending), and III (upbending) SB profiles. We compare radial profiles of the mass fractions and the velocity dispersions of different sub-populations of stars according to their birth and current location. Results. Radial redistribution of stars is important in all galaxies regardless of their light profiles. Type-II breaks seem to be a consequence of the combined effects of outward-moving and accreted stars. The former produce shallower inner profiles (lack of stars in the inner disc) and accumulate material around the break radius and beyond, strengthening the break; the latter can weaken or even convert the break into a pure exponential. Further accretion from satellites can concentrate material in the outermost parts, leading to type-III breaks that can coexist with type-II breaks, but situated further out. Type-III galaxies would be the result of an important radial redistribution of material throughout the entire disc, as well as a concentration of accreted material in the outskirts. In addition, type-III galaxies display the most efficient radial redistribution and the largest number of accreted stars, followed by type-I and II systems, suggesting that type-I galaxies may be an intermediate case between types-II and III. In general, the velocity dispersion profiles of all galaxies tend to flatten or even increase around the locations where the breaks are found. The age and metallicity profiles are also affected, exhibiting different inner gradients depending on their SB profile, being steeper in the case of type-II systems (as found observationally). The steep type-II profiles might be inherent to their formation rather than acquired via radial redistribution

Evolution of starburst galaxies in the Illustris simulation

© 2018 The Author(s). Published by Oxford University Press on behalf of The Royal Astronomical Society. There is a consensus in the literature that starburst galaxies are triggered by interaction events. However, it remains an open question as to what extent both merging and nonmerging interactions have in triggering starbursts. In this study, we make use of the Illustris simulation to test how different triggeringmechanisms can affect starburst events.We examine the star formation rate, colour, and environment of starburst galaxies to determine if this could be why we witness a bimodality in post-starburst populations within observational studies. Further, we briefly test the extent of quenching due to active galactic nuclei feedback. From Illustris, we select 196 starburst galaxies at z = 0.15 and split them into post-merger and pre-merger/harassment-driven starburst samples. We find that 55 % of this sample have not undergone a merger in the past 2 Gyr. Both of our samples are located in low-density environments within the filament regions of the cosmic web; however, we find that premerger/ harassment-driven starbursts are in higher-density environments than post-mergerdriven starbursts. We also find that pre-merger/harassment starbursts are redder than postmerger starbursts; this could be driven by environmental effects. Both, however, produce nuclear starbursts of comparable strengths

The stellar metallicity distribution of disc galaxies and bulges in cosmological simulations

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.Comment: 20 pages, 14 figures, MNRAS, accepte

Risk map as a library management information dashboard: a case study in adapting a configural display

In this paper, we report on our application of Cognitive Work Analysis to create an Abstraction Hierarchy model that helps librarians identify key functional relationships for managing the overall performance of a library. By themselves, functional relationships are not as useful in providing insights into the reasons for good or poor performance. However, when these functional relationships are set against the context of system invariants and constraints, they can provide library managers information useful for diagnosis and localization of problems. We propose the Risk Map visualization technique as an information dashboard to cognitively access these functional relationships. Furthermore, when these functional relationships are portrayed over time, trends and patterns can be detected with relative ease

Floods and health in Gambella region, Ethiopia: a qualitative assessment of the strengths and weaknesses of coping mechanisms

BACKGROUND: Floods are the most frequent and devastating type of natural disaster worldwide, causing unprecedented deaths, diseases, and destruction of property and crops. Flooding has a greater impact in developing countries due to lack of sufficient disaster management structures and a lack of economic resources. OBJECTIVE: This study was conducted with the aim of contributing to the knowledge base of development strategies that reduce flood-related health risks in developing countries. The study focused particularly on assessing the flood risks and health-related issues in the Gambella region of Ethiopia; with the intent of producing relevant information to assist with the improvements in the efficacy of the current flood coping strategies in the region. METHODS: Data were gathered through interviews with 14 officers from different government and non-governmental organizations and a questionnaire survey given to 35 flood victims in Itang woreda. A qualitative approach was applied and the data were analyzed using content analysis. RESULTS: It was found that flooding is a common problem in Gambella region. The findings also indicate that the flood frequency and magnitude has increased rapidly during the last decade. The increase in floods was driven mainly by climate change and changes in land use, specifically deforestation. The reported main impacts of flooding on human health in Gambella region were deaths, injuries, and diseases such as malaria and diarrhea. Another notable consequence of flooding was crop destruction and subsequent malnutrition. CONCLUSIONS: Three weaknesses that were identified in the current coping strategies for flood-related health impacts in Gambella region were a lack of flood-specific policy, absence of risk assessment, and weak institutional capacity. This study recommends new policy approaches that will increase the effectiveness of the current flood coping strategies to sustainably address the impact of flooding on human health

The properties, origin and evolution of stellar clusters in galaxy simulations and observations

We investigate the properties and evolution of star particles in two simulations of isolated spiral galaxies, and two galaxies from cosmological simulations. Unlike previous numerical work, where typically each star particle represents one ‘cluster’, for the isolated galaxies we are able to model features we term ‘clusters’ with groups of particles. We compute the spatial distribution of stars with different ages, and cluster mass distributions, comparing our findings with observations including the recent LEGUS survey. We find that spiral structure tends to be present in older (100s Myrs) stars and clusters in the simulations compared to the observations. This likely reflects differences in the numbers of stars or clusters, the strength of spiral arms, and whether the clusters are allowed to evolve. Where we model clusters with multiple particles, we are able to study their evolution. The evolution of simulated clusters tends to follow that of their natal gas clouds. Massive, dense, long-lived clouds host massive clusters, whilst short-lived clouds host smaller clusters which readily disperse. Most clusters appear to disperse fairly quickly, in basic agreement with observational findings. We note that embedded clusters may be less inclined to disperse in simulations in a galactic environment with continuous accretion of gas onto the clouds than isolated clouds and correspondingly, massive young clusters which are no longer associated with gas tend not to occur in the simulations. Caveats of our models include that the cluster densities are lower than realistic clusters, and the simplistic implementation of stellar feedback