The South African Astronomical Observatory 200 Virtual Symposium

No abstrac

The effect of the f-component of the pseudopotential on selected properties of 5d transition metal systems

Thesis (M.Sc.) - University of KwaZulu-Natal, Pietermaritzburg, 2008.Cohesive energies, bulk moduli, and equilibrium lattice constants have been calculated for the 5d transition atoms (Hf, Ta, W, Re, Os, Ir and Pt) in face–centred cubic crystal lattices. We have used the ab initio pseudopotential method for the total energy calculations within the local density approximation. Two calculations have been performed for each element, one using only the s, p and d angular momentum components and another including the s, p and d components as well as the unoccupied 5f orbital in the ionic pseudopotentials. The pseudo–wave functions and charge densities of the valence electrons have been represented by a basis of plane waves. For the 5d metals the changes in the electronic structure of the solid are small and they produce small changes in the bulk properties

Submitted in fulfillment of the requirements for the degree of Doctor of Philosophy in the School of Physics, University of the Western Cape

Philosophiae Doctor - PhDNumerical simulations play a crucial role in testing current cosmological models of the formation and evolution of the cosmic structure observed in the modern Universe. Simulations of the collapse of both baryonic and non-baryonic matter under the influence of gravity have yielded important results in our understanding of the large scale structure of the Universe. In addition to the underlying large scale structure, simulations which include gas dynamics can give us valuable insight into, and allow us to make testable predictions on, the nature and distribution of baryonic matter on a wide range of scales. In this work we give an overview of cosmological simulations and the methods employed in the solution of many body problems. We then present three projects focusing on scales ranging from individual galaxies to the cosmic web connecting clusters of galaxies thereby demonstrating the potential and diversity of numerical simulations in the fields of cosmology and astrophysics. We firstly investigate the environmental dependance of neutral hydrogen in the intergalactic medium by utilising high resolution cosmological hydrodynamic simulations in Chapter 3. We find that the extent of the neutral hydrogen radial profile is dependant on both the environment of the galaxy and its classification within the group ie. whether it is a central or satellite galaxy. We investigate whether this effect could arise from ram pressure forces exerted on the galaxies and find good agreement between galaxies experiencing high ram pressure forces and those with a low neutral hydrogen content. In Chapter 4 we investigate the velocity–shape alignment of clusters in a dark matter only simulation and the effect of such an alignment on measurements of the kinetic Sunyaev–Zeldovich (kSZ) effect. We find an alignment not only exists but can lead to an enhancement in the kSZ signal of up to 60% when the cluster is orientated along the line-of-sight. Finally we attempt to identify shocked gas in clusters and filaments using intermediate resolution cosmological hydrodynamic simulations in Chapter 5 with a view to predicting the synchrotron emission from these areas, something that may be detectable with the Square Kilometer Array

Quest Volume 19 Number 1

Department of Science and Innovation; Academy of Science of South Afric

Galaxy and Mass Assembly (GAMA): formation and growth of elliptical galaxies in the group environment

There are many proposed mechanisms driving the morphological transformation of disc galaxies to elliptical galaxies. In this paper, we determine if the observed transformation in low-mass groups can be explained by the merger histories of galaxies. We measured the group mass– morphology relation for groups from the Galaxy and Mass Assembly group catalogue with masses from 1011 to 1015 M. Contrary to previous studies, the fraction of elliptical galaxies in our more complete group sample increases significantly with group mass across the full range of group mass. The elliptical fraction increases at a rate of 0.163 ± 0.012 per dex of group mass for groups more massive than 1012.5 M. If we allow for uncertainties in the observed group masses, our results are consistent with a continuous increase in elliptical fraction from group masses as low as 1011 M. We tested if this observed relation is consistent with the merger activity using a GADGET-2 dark matter simulation of the galaxy groups. We specified that a simulated galaxy would be transformed to an elliptical morphology either if it experienced a major merger or if its cumulative mass gained from minor mergers exceeded 30 per cent of its final mass. We then calculated a group mass–morphology relation for the simulations. The position and slope of the simulated relation were consistent with the observational relation, with a gradient of 0.184 ± 0.010 per dex of group mass. These results demonstrate a strong correlation between the frequency of merger events and disc-to-elliptical galaxy transformation in galaxy group environments.This research was conducted by the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020. SB acknowledges funding support from the Australian Research Council through a Future Fellowship (FT140101166). GAMA is a joint European-Australasian project based around a spectroscopic campaign using the AAT. The GAMA input catalogue is based on data taken from the SDSS and the UKIRT Infrared Deep Sky Survey. Complementary imaging of the GAMA regions is being obtained by a number of independent survey programmes including GALEX MIS, VST KIDS, VISTA VIKING, WISE, Herschel-ATLAS, GMRT and ASKAP providing ultraviolet to radio coverage. GAMA is funded by the STFC (UK), the ARC (Australia), the AAO and the participating institutions. The GAMA web site is http://www.gama-survey.org/

nIFTy galaxy cluster simulations – II. Radiative models

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $\Lambda$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.Comment: 20 pages, 13 figures, submitted to MNRA

nIFTy galaxy cluster simulations - IV. Quantifying the influence of baryons on halo properties

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern smoothed particle hydrodynamics (SPH) approaches. For each code represented we have a dark-matteronly (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<≈ 20 per cent) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include active galactic nucleus feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers

Cosmic CARNage I: on the calibration of galaxy formation models

We present a comparison of nine galaxy formation models, eight semi-analytical, and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of comoving width 125h−1 Mpc, with a dark-matter particle mass of 1.24 × 109h−1M) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some ‘memory’ of any previous calibration that served as the starting point (especially for the manually tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3σ. The second calibration extended the observational data to include the z = 2 SMF alongside the z ∼ 0 star formation rate function, cold gas mass, and the black hole–bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models, we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.Instituto de Astrofísica de La PlataFacultad de Ciencias Astronómicas y Geofísica

Cosmic CARNage I: on the calibration of galaxy formation models

We present a comparison of nine galaxy formation models, eight semi-analytical, and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of comoving width 125h−1 Mpc, with a dark-matter particle mass of 1.24 × 109h−1M) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some ‘memory’ of any previous calibration that served as the starting point (especially for the manually tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3σ. The second calibration extended the observational data to include the z = 2 SMF alongside the z ∼ 0 star formation rate function, cold gas mass, and the black hole–bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models, we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.Instituto de Astrofísica de La PlataFacultad de Ciencias Astronómicas y Geofísica

Cosmic CARNage I: on the calibration of galaxy formation models

We present a comparison of nine galaxy formation models, eight semi-analytical, and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of comoving width 125h−1 Mpc, with a dark-matter particle mass of 1.24 × 109h−1M) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some ‘memory’ of any previous calibration that served as the starting point (especially for the manually tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3σ. The second calibration extended the observational data to include the z = 2 SMF alongside the z ∼ 0 star formation rate function, cold gas mass, and the black hole–bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models, we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.Instituto de Astrofísica de La PlataFacultad de Ciencias Astronómicas y Geofísica