The Galactic bulge is the least studied component of our Galaxy. Yet, its formation and evolution are key to understand the formation of the Galaxy itself. Studies on the Galactic bulge have increased significantly in the last years, but still there are many points of controversy. This volume contains several contributions from experts in different aspects of the bulge. Issues discussed include the following: the presence of an old spheroidal bulge, or identification of its old stellar population with the thick disk or halo; fraction of stars younger than 10 Gyr is estimated to be of < 5 to 22% depending on method and authors; multiple populations or only a metal-poor and a metal-rich ones; spheroidal or ellipsoidal distribution of RR Lyrae; formation of the bulge from early mergers or from secular evolution of the bar; different methods of mapping extinction; selection and identification of bulge globular clusters.

The Galactic bulge of the Milky Way is made up of stars with a broad range of metallicity, –3.0 < [Fe/H] < 1 dex. The mean of the metallicity distribution function decreases as a function of height z from the plane and, more weakly, with galactic radius RGC. The most metal-rich stars in the inner Galaxy are concentrated to the plane and the more metal-poor stars are found predominantly further from the plane, with an overall vertical gradient in the mean of the metallicity distribution function of about − 0.45 dex kpc−1. This vertical gradient is believed to reflect the changing contribution with height of different populations in the innermost region of the Galaxy. The more metal-rich stars of the bulge are part of the boxy/peanut structure and comprise stars in orbits which trace out the underlying X-shape. There is still a lack of consensus on the origin of the metal-poor stars ([Fe/H] < −0.5) in the region of the bulge. Some studies attribute the more metal-poor stars of the bulge to the thick disk and stellar halo that are present in the inner region, and other studies propose that the metal-poor stars are a distinct ‘old spheroid’ bulge population. Understanding the origin of the populations that make up the metallicity distribution function of the bulge, and identifying if there is a unique bulge population which has formed separately from the disk and halo, has important consequences for identifying the relevant processes in the formation and evolution of the Milky Way.

The stellar population of the Milky Way bulge is thoroughly studied, with a plethora of measurements from virtually the full suite of instruments available to astronomers. It is thus perhaps surprising that alongside well-established results lies some substantial uncertainty in its star-formation history. Cosmological models predict the bulge to host the Galaxy's oldest stars for [Fe/H] ≲ −1, and this is demonstrated by RR Lyrae stars and globular cluster observations. There is consensus that bulge stars with [Fe/H] ≲ 0 are older than t ≈ 10 Gyr. However, at super-solar metallicity, there is a substantial unresolved discrepancy. Data from spectroscopic measurements of the main-sequence turnoff and subgiant branch, the abundances of asymptotic giant branch stars, the period distribution of Mira variables, the chemistry and central-star masses of planetary nebulae, all suggest a substantial intermediate-age population (t ≈ 3 Gyr). This is in conflict with predictions from cosmologically motivated chemical evolution models and photometric studies of the main-sequence turnoff region, which both suggest virtually no stars younger than t ≈ 8 Gyr. A possible resolution to this conflict is enhanced helium-enrichment, as this would shift nearly all of the age estimates in the direction of decreasing discrepancy.

I review the literature covering the issue of interstellar extinction towards the Milky Way bulge, with emphasis placed on findings from planetary nebulae, RR Lyrae, and red clump stars. I also report on observations from HI gas and globular clusters. I show that there has been substantial progress in this field in recent decades, most particularly from red clump stars. The spatial coverage of extinction maps has increased by a factor ~ 100 × in the past 20 yr, and the total-to-selective extinction ratios reported have shifted by ~ 20–25%, indicative of the improved accuracy and separately, of a steeper-than-standard extinction curve. Problems remain in modelling differential extinction, explaining anomalies involving the planetary nebulae, and understanding the difference between bulge extinction coefficients and ‘standard’ literature values.

We review the observational evidences concerning the three-dimensional structure of the Galactic bulge. Although the inner few kpc of our Galaxy are normally referred to as the bulge, all the observations demonstrate that this region is dominated by a bar, i.e., the bulge is a bar. The bar has a boxy/peanut (X-shaped) structure in its outer regions, while it seems to become less and less elongated in its innermost region. A thinner and longer structure departing from the main bar has also been found, although the observational evidences that support the scenario of two separate structures has been recently challenged. Metal-poor stars ([Fe/H] ≲ −0.5 dex) trace a different structure, and also have different kinematics.

Recent large-scale surveys of galactic bulge stars allowed to build a detailed map of the bulge kinematics. The bulge exhibits cylindrical rotation consistent with a disky origin which evolved through bar-driven secular evolution. However, correlations between metallicity and kinematics complicate this picture. In particular a metal-poor component with distinct kinematic signatures has been detected. Its origin, density profile and link with the other Milky Way stellar populations are currently still poorly constrained.

The Galactic bulge, that is the prominent out-of-plane over-density present in the inner few kiloparsecs of the Galaxy, is a complex structure, as the morphology, kinematics, chemistry, and ages of its stars indicate. To understand the nature of its main components—those at [Fe/H] ≳ −1 dex—it is necessary to make an inventory of the stellar populations of the Galactic disc(s), and of their borders: the chemistry of the disc at the solar vicinity, well known from detailed studies of stars over many years, is not representative of the whole disc. This finding, together with the recent revisions of the mass and sizes of the thin and thick discs, constitutes a major step in understanding the bulge complexity. N-body models of a boxy-/peanut-shaped bulge formed from a thin disc through the intermediary of a bar have been successful in interpreting a number of global properties of the Galactic bulge, but they fail in reproducing the detailed chemo-kinematic relations satisfied by its components and their morphology. It is only by adding the thick disc to the picture that we can understand the nature of the Galactic bulge.

A view of the Galactic bulge by means of their globular clusters is fundamental for a deep understanding of its formation and evolution. Connections between the globular cluster and field star properties in terms of kinematics, orbits, chemical abundances, and ages should shed light on different stellar population components. Based on spatial distribution and metallicity, we define a probable best list of bulge clusters, containing 43 entries. Future work on newly discovered objects, mostly from the VVV survey, is suggested. These candidates might alleviate the issue of missing clusters on the far side of the bulge. We discuss the reddening law affecting the cluster distances towards the centre of the Galaxy, and conclude that the most suitable total-to-selective absorption value appears to be RV=3.2, in agreement with recent analyses. An update of elemental abundances for bulge clusters is provided.

The relation between the star-formation Rate and stellar mass (M⋆) of galaxies represents a fundamental constraint on galaxy formation, and has been studied extensively both in observations and cosmological hydrodynamic simulations. However, the observed amplitude of the star-formation rate—stellar mass relation has not been successfully reproduced in simulations, indicating either that the halo accretion history and baryonic physics are poorly understood/modelled or that observations contain biases. In this paper, we examine the evolution of the SFR − M⋆ relation of z ~ 1–4 galaxies and display the inconsistency between observed relations that are obtained using different techniques. We employ cosmological hydrodynamic simulations from various groups which are tuned to reproduce a range of observables and compare these with a range of observed SFR − M⋆ relations. We find that numerical results are consistent with observations that use Spectral Energy Distribution techniques to estimate star-formation rates, dust corrections, and stellar masses. On the contrary, simulations are not able to reproduce results that were obtained by combining only UV and IR luminosities (UV+IR). These imply star-formation rates at a fixed stellar mass that are larger almost by a factor of 5 than those of Spectral Energy Distribution measurements for z ~ 1.5–4. For z < 1.5, the results from simulations, Spectral Energy Distribution fitting techniques and IR+UV conversion agree well. We find that surveys that preferably select star-forming galaxies (e.g. by adopting Lyman-break or blue selection) typically predict a larger median/average star-formation rate at a fixed stellar mass especially for high mass objects, with respect to mass selected samples and hydrodynamic simulations. Furthermore, we find remarkable agreement between the numerical results from various authors who have employed different cosmological codes and run simulations with different resolutions. This is interesting for two reasons. (A) simulations can produce realistic populations of galaxies within representative cosmological volumes even at relatively modest resolutions. (B) It is likely that current numerical codes that rely on similar subgrid multiphase interstellar medium models and are tuned to reproduce statistical properties of galaxies, produce similar results for the SFR − M⋆ relation by construction, regardless of resolution, box size and, to some extent, the adopted feedback prescriptions.

The Millimetre Astronomy Legacy Team 90 GHz survey aims to characterise the physical and chemical evolution of high-mass clumps. Recently completed, it mapped 90 GHz line emission towards 3 246 high-mass clumps identified from the ATLASGAL 870 μm Galactic plane survey. By utilising the broad frequency coverage of the Mopra telescope’s spectrometer, maps in 16 different emission lines were simultaneously obtained. Here, we describe the first catalogue of the detected line emission, generated by Gaussian profile fitting to spectra extracted towards each clumps’ 870 μm dust continuum peak. Synthetic spectra show that the catalogue has a completeness of > 95%, a probability of a false-positive detection of < 0.3%, and a relative uncertainty in the measured quantities of < 20% over the range of detection criteria. The detection rates are highest for the (1–0) transitions of HCO+, HNC, N2H+, and HCN (~77–89%). Almost all clumps (~95%) are detected in at least one of the molecular transitions, just over half of the clumps (~53%) are detected in four or more of the transitions, while only one clump is detected in 13 transitions. We find several striking trends in the ensemble of properties for the different molecular transitions when plotted as a function of the clumps’ evolutionary state as estimated from Spitzer mid-IR images, including (1) HNC is relatively brighter in colder, less evolved clumps than those that show active star formation, (2) N2H+ is relatively brighter in the earlier stages, (3) that the observed optical depth decreases as the clumps evolve, and (4) the optically thickest HCO+ emission shows a ‘blue-red asymmetry’ indicating overall collapse that monotonically decreases as the clumps evolve. This catalogue represents the largest compiled database of line emission towards high-mass clumps and is a valuable data set for detailed studies of these objects.

We estimate spatial gradients in the ionosphere using the Global Positioning System and GLONASS (Russian global navigation system) observations, utilising data from multiple Global Positioning System stations in the vicinity of Murchison Radio-astronomy Observatory. In previous work, the ionosphere was characterised using a single-station to model the ionosphere as a single layer of fixed height and this was compared with ionospheric data derived from radio astronomy observations obtained from the Murchison Widefield Array. Having made improvements to our data quality (via cycle slip detection and repair) and incorporating data from the GLONASS system, we now present a multi-station approach. These two developments significantly improve our modelling of the ionosphere. We also explore the effects of a variable-height model. We conclude that modelling the small-scale features in the ionosphere that have been observed with the MWA will require a much denser network of Global Navigation Satellite System stations than is currently available at the Murchison Radio-astronomy Observatory.

The formation, assembly history, and environmental impact of the massive black holes (BH) that are ubiquitous in the nuclei of luminous galaxies today remain some of the main unsolved problems in cosmic structure formation. In the last several years, it has become clear that quasars are not just tracers of early and recent structure formation, but that they seem to have actively influenced galaxies and clusters through feedback mechanisms that are still not well understood. The discovery of more and more numerous quasars at redshift above 6, powered by BHs with masses similar to that of their local counterparts, further complicates this scenario. This emphasises the urgent need to better understand how and when such massive objects form and grow, what is the strength and scale of their impact on the evolution of their host galaxies, and what are the main physical processes driving and regulating this co-evolution.

In this paper, the hypothesis of the existence of a massive dark body (Nemesis, Tyche, Planet Nine, or any other trans-Plutonian planet) at the Solar system periphery is analysed. Basic physical properties and orbital characteristics of such massive bodies are considered. The problem of the definition of a scattering angle of a photon in the gravitational field of a spherical lens is studied. It is shown that, the required value of the scattering angle can be measured for the cases of Nemesis and Tyche. The formation of gravitational lensing images is studied here for a point mass event. It is demonstrated that in most cases of the close rapprochement of a source and the lens (for Nemesis and Tyche), it is possible to resolve two images. The possibility of resolving these images is one of the main arguments favouring the gravitational lensing method as its efficiency in searching for dark massive objects at the edge of the Solar System is higher than the one corresponding to other methods such as stellar occultation. For the cases of Planet Nine and any other trans-Plutonian planet, the strong gravitational lensing is impossible because at least one of the images is always eclipsed.

Prestellar cores, by definition, are gravitationally bound but starless pockets of dense gas. Physical conditions that could render a core starless (in the local Universe) is the subject of investigation in this work. To this end, we studied the evolution of four starless cores, B68, L694-2, L1517B, L1689, and L1521F, a VeLLO. We demonstrate: (i) cores contracted in quasistatic manner over a timescale on the order of ~ 105 yr. Those that remained starless briefly acquired a centrally concentrated density configuration that mimicked the profile of a unstable BonnorEbert sphere before rebounding, (ii) three cores viz. L694-2, L1689-SMM16, and L1521F remained starless despite becoming thermally super-critical. By contrast, B68 and L1517B remained sub-critical; L1521F collapsed to become a VeLLO only when gas-cooling was enhanced by increasing the size of dust-grains. This result is robust, for other starless cores viz. B68, L694-2, L1517B, and L1689 could also be similarly induced to collapse. The temperature-profile of starless cores and those that collapsed was found to be radically different. While in the former type, only very close to the centre of a core was there any evidence of decline in gas temperature, by contrast, a core of the latter type developed a more uniformly cold interior. Our principle conclusions are: (a) thermal super-criticality of a core is insufficient to ensure it will become protostellar, (b) potential star-forming cores (the VeLLO L1521F here), could be experiencing dust-coagulation that must enhance gasdust coupling and in turn lower gas temperature, thereby assisting collapse. This also suggests, mere gravitational/virial boundedness of a core is insufficient to ensure it will form stars.

Integral field unit spectrographs allow the 2D exploration of the kinematics and stellar populations of galaxies, although they are generally restricted to small fields-of-view. Using the large field-of-view of the DEIMOS multislit spectrograph on Keck and our Stellar Kinematics using Multiple Slits technique, we are able to extract sky-subtracted stellar light spectra to large galactocentric radii. Here, we present a new DEIMOS mask design named SuperSKiMS that explores large spatial scales without sacrificing high spatial sampling. We simulate a set of observations with such a mask design on the nearby galaxy NGC 1023, measuring stellar kinematics and metallicities out to where the galaxy surface brightness is orders of magnitude fainter than the sky. With this technique we also reproduce the results from literature integral field spectroscopy in the innermost galaxy regions. In particular, we use the simulated NGC 1023 kinematics to model its total mass distribution to large radii, obtaining comparable results with those from published integral field unit observation. Finally, from new spectra of NGC 1023, we obtain stellar 2D kinematics and metallicity distributions that show good agreement with integral field spectroscopy results in the overlapping regions. In particular, we do not find a significant offset between our Stellar Kinematics using Multiple Slits and the ATLAS3D stellar velocity dispersion at the same spatial locations.

Making use of a new N-body model to describe the evolution of a moderate-size globular cluster, we investigate the characteristics of the population of black holes within such a cluster. This model reaches core-collapse and achieves a peak central density typical of the dense globular clusters of the Milky Way. Within this high-density environment, we see direct confirmation of the merging of two stellar remnant black holes in a dynamically formed binary, a gravitational wave source. We describe how the formation, evolution, and ultimate ejection/destruction of binary systems containing black holes impacts the evolution of the cluster core. Also, through comparison with previous models of lower density, we show that the period distribution of black hole binaries formed through dynamical interactions in this high-density model favours the production of gravitational wave sources. We confirm that the number of black holes remaining in a star cluster at late times and the characteristics of the binary black hole population depend on the nature of the star cluster, critically on the number density of stars and by extension the relaxation timescale.

Although the universe at redshifts greater than six represents only the first one billion years (< 10%) of cosmic time, the dense nature of the early universe led to vigorous galaxy formation and evolution activity which we are only now starting to piece together. Technological improvements have, over only the past decade, allowed large samples of galaxies at such high redshifts to be collected, providing a glimpse into the epoch of formation of the first stars and galaxies. A wide variety of observational techniques have led to the discovery of thousands of galaxy candidates at z > 6, with spectroscopically confirmed galaxies out to nearly z = 9. Using these large samples, we have begun to gain a physical insight into the processes inherent in galaxy evolution at early times. In this review, I will discuss (i) the selection techniques for finding distant galaxies, including a summary of previous and ongoing ground and space-based searches, and spectroscopic follow-up efforts, (ii) insights into galaxy evolution gleaned from measures such as the rest-frame ultraviolet luminosity function, the stellar mass function, and galaxy star-formation rates, and (iii) the effect of galaxies on their surrounding environment, including the chemical enrichment of the universe, and the reionisation of the intergalactic medium. Finally, I conclude with prospects for future observational study of the distant universe, using a bevy of new state-of-the-art facilities coming online over the next decade and beyond.

We identify gravitationally bound structures in the Ursa Major region using positions, velocities and photometry from the Sloan Digital Sky Survey (SDSS DR7) and the Third Reference Catalogue of Bright Galaxies (RC3). A friends-of-friends algorithm is extensively tested on mock galaxy lightcones and then implemented on the real data to determine galaxy groups whose members are likely to be physically and dynamically associated with one another. We find several galaxy groups within the region that are likely bound to one another and in the process of merging. We classify 6 galaxy groups as the Ursa Major ‘supergroup’, which are likely to merge and form a poor cluster with a mass of ~ 8 × 1013 M⊙. Furthermore, the Ursa Major supergroup as a whole is likely bound to the Virgo cluster, which will eventually form an even larger system in the context of hierarchical structure formation. We investigate the evolutionary state of the galaxy groups in the Ursa Major region and conclude that these groups are in an early evolutionary state and the properties of their member galaxies are similar to those in the field.

The traditional cultures of Aboriginal Australians include a significant astronomical component, perpetuated through oral tradition, ceremony, and art. This astronomical knowledge includes a deep understanding of the motion of objects in the sky, which was used for practical purposes such as constructing calendars and for navigation. There is also evidence that traditional Aboriginal Australians made careful records and measurements of cyclical phenomena, recorded unexpected phenomena such as eclipses and meteorite impacts, and could determine the cardinal points to an accuracy of a few degrees. Putative explanations of celestial phenomena appear throughout the oral record, suggesting traditional Aboriginal Australians sought to understand the natural world around them, in the same way as modern scientists, but within their own cultural context. There is also a growing body of evidence for sophisticated navigational skills, including the use of astronomically based songlines. Songlines are effectively oral maps of the landscape, and are an efficient way of transmitting oral navigational skills in cultures that do not have a written language. The study of Aboriginal astronomy has had an impact extending beyond mere academic curiosity, facilitating cross-cultural understanding, demonstrating the intimate links between science and culture, and helping students to engage with science.

At a bulge latitude of b = −4°, the average [Fe/H] and [Mg/H] values are +0.06 and +0.17 dex, roughly 0.2 and 0.7 dex higher than the local thin and thick disk values, respectively, suggesting a large bulge effective yield, perhaps due to efficient retention of supernova ejecta.

The bulge vertical [Fe/H] gradient, at ~ 0.5 dex/kpc, appears to be due to a changing mixture of sub-populations (near +0.3 dex and −0.3 dex and one possibly near −0.7 dex) with latitude. At solar [Fe/H], the bulge [Al/Fe] and [α/Fe] ratios are ~ +0.15 dex. Below [Fe/H] < ![CDATA[$ ~ −0.5 dex, the bulge and local thick disk compositions are very similar; but the measured [Mg/Fe], [⟨SiCaTi⟩/Fe], [La/Eu] and dramatic [Cu/Fe] ratios suggest higher SFR in the bulge. However, these composition differences with the thick disk could be due to measurement errors and non-LTE effects.

Unusual zig-zag trends of [Cu/Fe] and [Na/Fe] suggest metallicity-dependent nucleosynthesis by core-collapse supernovae in the Type Ia supernova time-delay scenario.

The bulge sub-population compositions resemble the local thin and thick disks, but at higher [Fe/H], suggesting a radial [Fe/H] gradient of − $0.04]] > to − 0.05 dex/kpc for both the thin and thick disks. If the bulge formed through accretion of inner thin and thick disk stars, it appears that these stars retained vertical scale heights characteristic of their kinematic origin, resulting in the vertical [Fe/H] gradient and [α/Fe] trends seen today.