Search

Chemical abundances are an essential tool in untangling the Milky Way's enrichment history. However, the evolution of the interstellar medium abundance gradient with cosmic time is lost as a result of radial mixing processes. For the first time, we quantify the evolution of many observational abundances across the Galactic disk as a function of lookback time and birth radius, $R_\text{birth}$. Using an empirical approach, we derive $R_\text{birth}$ estimates for 145,447 APOGEE DR17 red giant disk stars, based solely on their ages and [Fe/H]. We explore the detailed evolution of 6 abundances (Mg, Ca ($\alpha$), Mn (iron-peak), Al, C (light), Ce (s-process)) across the Milky Way disk using 87,426 APOGEE DR17 red giant stars. We discover that the interstellar medium had three fluctuations in the metallicity gradient $\sim 9$, $\sim 6$, and $\sim4$ Gyr ago. The first coincides with the end of high-$\alpha$ sequence formation around the time of the Gaia-Sausage-Enceladus disruption, while the others are likely related to passages of the Sagittarius dwarf galaxy. A clear distinction is found between present-day observed radial gradients with age and the evolution with lookback time for both [X/Fe] and [X/H], resulting from the significant flattening and inversion in old populations due to radial migration. We find the [Fe/H]--[$\alpha$/Fe] bimodality is also seen as a separation in the $R_\text{birth}$--[X/Fe] plane for the light and $\alpha$-elements. Our results recover the chemical enrichment of the Galactic disk over the past 12 Gyr, providing tight constraints on Galactic disk chemical evolution models., Comment: accepted for publication in MNRAS

Recovering the birth radii of observed stars in the Milky Way is one of the ultimate goals of Galactic Archaeology. One method to infer the birth radius and the evolution of the interstellar medium (ISM) metallicity assumes a linear relation between the ISM metallicity with radius at any given look-back time. Here, we test the reliability of this assumption by using four zoom-in cosmological hydrodynamic simulations from the NIHAO-UHD project. We find that one can infer precise birth radii only when the stellar disc starts to form, which for our modelled galaxies happens ∼10 Gyr ago, in agreement with recent estimates for the Milky Way. With a current day measurement of ISM metallicity gradient of −0.05 dex and a dispersion of 0.03 dex, the intrinsic uncertainty in inferring Rbirth is ∼0.6 kpc. At later times, the linear correlation between the ISM metallicity and radius increases, as stellar motions become more ordered and the azimuthal variations of the ISM metallicity start to drop. The formation of a central bar and perturbations from mergers can increase this uncertainty in the inner and outer disc, respectively.

The present-day age–metallicity relation (AMR) is a record of the star formation history of galaxies, as this traces the chemical enrichment of the gas over time. We use a zoomed-in cosmological simulation that reproduces key signatures of the Milky Way (MW), g2.79e12 from the NIHAO-UHD project, to examine how stellar migration and satellite infall shape the AMR across the disc. We find in the simulation, similar to the MW, the AMR in small spatial regions (R, z) shows turning points that connect changes in the direction of the relations. The turning points in the AMR in the simulation are a signature of late satellite infall. This satellite infall has a mass radio similar as that of the Sagittarius dwarf to the MW (∼0.001). Stars in the apex of the turning points are young and have nearly not migrated. The late satellite infall creates the turning points via depositing metal-poor gas in the disc, triggering star formation of stars in a narrow metallicity range compared to the overall AMR. The main effect of radial migration on the AMR turning points is to widen the metallicity range of the apex. This can happen when radial migration brings stars born from the infallen gas in other spatial bins, with slightly different metallicities, into the spatial bin of interest. These results indicate that it is possible that the passage of the Sagittarius dwarf galaxy played a role in creating the turning points that we see in the AMR in the Milky Way.

Stars move away from their birth places over time via a process known as radial migration, which blurs chemo-kinematic relations used for reconstructing the Milky Way formation history. One of the ultimate goals of Galactic Archaeology, therefore, is to find stars’ birth aggregates in the disk via chemical tagging. Here we show that stellar birth radii can be derived directly from the data with minimum prior assumptions on the Galactic enrichment history. We recover the time evolution of the stellar birth metallicity gradient, d[Fe/H](R, τ)/dR, through its inverse relation to the metallicity range as a function of age today, allowing us to place any star with age and metallicity measurements back to its birthplace, Rb. Applying our method to a high-precision large data set of Milky Way disk subgiant stars, we find a steepening of the birth metallicity gradient from 11 to 8 Gyr ago, which coincides with the time of the last major merger, Gaia-Sausage-Enceladus (GSE). This transition appears to play a major role in shaping both the age-metallicity relation and the bimodality in the [α/Fe]-[Fe/H] plane. By dissecting the disk into mono-Rb populations, clumps in the low-[α/Fe] sequence appear, which are not seen in the total sample and coincide in time with known star-formation bursts. We estimated that the Sun was born at 4.5 +- 0.4 kpc from the Galactic center. Our Rb estimates provide the missing piece needed to recover the Milky Way formation history, while the by-product, [Fe/H](R, τ), can be used as the thus-far missing prior for chemical evolution modeling.

With the advent of large spectroscopic surveys the amount of high quality chemo-dynamical data in the Milky Way (MW) increased tremendously. Accurately and correctly capturing and explaining the detailed features in the high-quality observational data is notoriously difficult for state-of-the-art numerical models. In order to keep up with the quantity and quality of observational datasets, improved prescriptions for galactic chemical evolution need to be incorporated into the simulations. Here we present a new, flexible, time resolved chemical enrichment model for cosmological simulations. Our model allows to easily change a number of stellar physics parameters such as the shape of the initial mass function (IMF), stellar lifetimes, chemical yields or SN Ia delay times. We implement our model into the Gasoline2 code and perform a series of cosmological simulations varying a number of key parameters, foremost evaluating different stellar yield sets for massive stars from the literature. We find that total metallicity, total iron abundance and gas phase oxygen abundance are robust predictions from different yield sets and in agreement with observational relations. On the other hand, individual element abundances, especially $\alpha$-elements show significant differences across different yield sets and none of our models can simultaneously match constraints on the dwarf and MW mass scale. This offers a unique way of observationally constraining model parameters. For MW mass galaxies we find for most yield tables tested in this work a bimodality in the $[\alpha$/Fe] vs. [Fe/H] plane of rather low intrinsic scatter potentially in tension with the observed abundance scatter., Comment: main text 19 pages and 11 figures, 4 pages of appendix, 23 pages total, python code and data at https://github.com/TobiBu/chemical_enrichment.git. submitted to MNRAS, comments very welcome

We investigate the impact of cosmic rays (CR) and different modes of CR transport on the properties of Milky Way-mass galaxies in cosmological magneto-hydrodynamical simulations in the context of the AURIGA project. We systematically study how advection, anisotropic diffusion and additional Alfv\'en-wave cooling affect the galactic disc and the circum-galactic medium (CGM). Global properties such as stellar mass and star formation rate vary little between simulations with and without various CR transport physics, whereas structural properties such as disc sizes, CGM densities or temperatures can be strongly affected. In our simulations, CRs affect the accretion of gas onto galaxies by modifying the CGM flow structure. This alters the angular momentum distribution which manifests itself as a difference in stellar and gaseous disc size. The strength of this effect depends on the CR transport model: CR advection results in the most compact discs while the Alfv\'en-wave model resembles more the AURIGA model. The advection and diffusion models exhibit large ($r\sim50$ kpc) CR pressure-dominated gas haloes causing a smoother and partly cooler CGM. The additional CR pressure smoothes small-scale density peaks and compensates for the missing thermal pressure support at lower CGM temperatures. In contrast, the Alfv\'en-wave model is only CR pressure dominated at the disc-halo interface and only in this model the gamma-ray emission from hadronic interactions agrees with observations. In contrast to previous findings, we conclude that details of CR transport are critical for accurately predicting the impact of CR feedback on galaxy formation., Comment: accepted by MNRAS, 26 pages (5 appendix), 20 figures (5 appendix)

Cosmological galaxy formation simulations are still limited by their spatial/mass resolution and cannot model from first principles some of the processes, like star formation, that are key in driving galaxy evolution. As a consequence they still rely on a set of 'effective parameters' that try to capture the scales and the physical processes that cannot be directly resolved in the simulation. In this study we show that it is possible to use Machine Learning techniques applied to real and simulated images of galaxies to discriminate between different values of these parameters by making use of the full information content of an astronomical image instead of collapsing it into a limited set of values like size, or stellar/ gas masses. In this work we apply our method to the NIHAO simulations and the THINGS and VLA-ANGST observations of HI maps in nearby galaxies to test the ability of different values of the star formation density threshold $n$ to reproduce observed HI maps. We show that observations indicate the need for a high value of $n \gtrsim 80$ ,cm$^{-3}$ (although the exact numerical value is model-dependent), which has important consequences for the dark matter distribution in galaxies. Our study shows that with innovative methods it is possible to take full advantage of the information content of galaxy images and compare simulations and observations in an interpretable, non-parametric and quantitative manner., Comment: 8 pages, 5 figures, accepted for publication on MNRAS

Recent observational and theoretical studies of the Local Group (LG) dwarf galaxies have highlighted their unique star-formation history, stellar metallicity, gas content, and kinematics. We investigate the commonality of these features by comparing constrained LG and field central dwarf halo simulations in the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project. Our simulations, performed with NIHAO-like hydrodynamics which track the evolution of the Milky Way (MW) and M31 along with ~100 dwarfs in the LG, reveal the total gas mass and stellar properties (velocity dispersion, evolution history, etc.) of present-day LG dwarfs to be similar to field systems. However, relative to field galaxies, LG dwarfs have more cold gas in their central parts and more metal-rich gas in the halo stemming from interactions with other dwarfs living in a high-density environment like the LG. Interestingly, the direct impact of massive MW/M31 analogues on the metallicity evolution of LG dwarfs is minimal; LG dwarfs accrete high-metallicity gas mostly from other dwarfs at late times. We have also tested for the impact of metal diffusion on the chemical evolution of LG dwarfs, and found that it does not affect the stellar or gaseous content of LG dwarfs. Our simulations suggest that the stellar components of LG dwarfs offer a unique and unbiased local laboratory for galaxy-formation tests and comparisons, especially against the overall dwarf population in the Universe., 16 pages, 12 figures, Accepted for publication in MNRAS

The Milky Way's stellar disk exhibits a bimodality in the [Fe/H] vs. [$\alpha$/Fe] plane, showing a distinct high-$\alpha$ and low-$\alpha$ sequence whose origin is still under debate. We examine the [Fe/H]-[$\alpha$/Fe] abundance plane in cosmological hydrodynamical simulations of Milky Way like galaxies from the NIHAO-UHD project and show that the bimodal $\alpha$-sequence is a generic consequence of a gas-rich merger at some time in the Galaxy's evolution. The high-$\alpha$ sequence evolves first in the early galaxies, extending to high metallicities, while it is the low-$\alpha$ sequence that is formed after the gas-rich merger. The merger brings in fresh metal-poor gas diluting the interstellar medium's metallicity while keeping the [$\alpha$/Fe] abundance almost unchanged. The kinematic, structural and spatial properties of the bimodal $\alpha$-sequence in our simulations reproduces that of observations. In all simulations, the high-$\alpha$ disk is old, radially concentrated towards the galaxy's center and shows large scale heights. In contrast, the low-$\alpha$ disk is younger, more radially extended and concentrated to the disk mid-plane. Our results show that the abundance plane is well described by these two populations that have been distributed radially across the disk by migration: at present-day in the solar neighbourhood, low-$\alpha$ stars originate from both the inner and outer disk while most of the high-$\alpha$ stars have migrated from the inner disk. We show that age dating the stars in the [Fe/H]-[$\alpha$/Fe] plane can constrain the time of the low-$\alpha$ sequence forming merger and conclude that $\alpha$-bimodality is likely a not uncommon feature of disk galaxies., Comment: 13 pages, 9 main text, 7 figures plus 3 in appendix, accepted by MNRAS

The phenomenological basis for Modified Newtonian Dynamics (MOND) is the radial-acceleration-relation (RAR) between the observed acceleration, $a=V^2_{rot}(r)/r$, and the acceleration accounted for by the observed baryons (stars and cold gas), $a_{bar}=V_{bar}^2(r)/r$. We show that the RAR arises naturally in the NIHAO sample of 89 high-resolution LCDM cosmological galaxy formation simulations. The overall scatter from NIHAO is just 0.079 dex, consistent with observational constraints. However, we show that the scatter depends on stellar mass. At high masses ($10^9, Comment: 15 pages, 15 figures, accepted to MNRAS

The kinematics of the most metal-poor stars provide a window into the early formation and accretion history of the Milky Way. Here, we use 5~high-resolution cosmological zoom-in simulations ($\sim~5\times10^6$ star particles) of Milky Way-like galaxies taken from the NIHAO-UHD project, to investigate the origin of low-metallicity stars ([Fe/H]$\leq-2.5$). The simulations show a prominent population of low-metallicity stars confined to the disk plane, as recently discovered in the Milky Way. The ubiquity of this finding suggests that the Milky Way is not unique in this respect. Independently of the accretion history, we find that $\gtrsim~90$ per cent of the retrograde stars in this population are brought in during the initial build-up of the galaxies during the first few Gyrs after the Big Bang. Our results therefore highlight the great potential of the retrograde population as a tracer of the early build-up of the Milky Way. The prograde planar population, on the other hand, is accreted during the later assembly phase and samples the full galactic accretion history. In case of a quiet accretion history, this prograde population is mainly brought in during the first half of cosmic evolution ($t\lesssim7$~Gyr), while, in the case of an on-going active accretion history, later mergers on prograde orbits are also able to contribute to this population. Finally, we note that the Milky Way shows a rather large population of eccentric, very metal-poor planar stars. This is a feature not seen in most of our simulations, with the exception of one simulation with an exceptionally active early building phase., Comment: Submitted to MNRAS, 8 figures, 13 pages

To understand the formation and evolution of the Milky Way disk, we must connect its current properties to its past. We explore hydrodynamical cosmological simulations to investigate how the chemical abundances of stars might be linked to their origins. Using hierarchical clustering of abundance measurements in two Milky Way-like simulations with distributed and steady star formation histories, we find that abundance clusters of stars comprise different groups in birth place ($R_\text{birth}$) and time (age). Simulating observational abundance errors (0.05 dex), we find that to trace discrete groups of ($R_\text{birth}$, age) requires a large vector of abundances. Using 15-element abundances (Fe, O, Mg, S, Si, C, P, Mn, Ne, Al, N, V, Ba, Cr, Co), up to $\approx$ 10 clusters can be defined with $\approx$ 25% overlap in ($R_\text{birth}$, age). We build a simple model to show that it is possible to infer a star's age and $R_\text{birth}$ from abundances with precisions of $\pm$0.06 Gyr and $\pm$1.17 kpc respectively. We find that abundance clustering is ineffective for a third simulation, where low-$\alpha$ stars form distributed in the disc and early high-$\alpha$ stars form more rapidly in clumps that sink towards the galactic center as their constituent stars evolve to enrich the interstellar medium. However, this formation path leads to large age-dispersions across the [$\alpha$/Fe]-[Fe/H] plane, which is inconsistent with the Milky Way's observed properties. We conclude that abundance clustering is a promising approach toward charting the history of our Galaxy., Comment: Submitted to ApJ

We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the response of cold dark matter (CDM) haloes to baryonic processes. Previous work has shown that the halo response is primarily a function of the ratio between galaxy stellar mass and total virial mass, and the density threshold above which gas is eligible to form stars, $n [{\rm cm}^{-3}]$. At low $n$ all simulations in the literature agree that dwarf galaxy haloes are cuspy, but at high $n\ge 100$ there is no consensus. We trace halo contraction in dwarf galaxies with $n\ge 100$ reported in some previous simulations to insufficient spatial resolution. Provided the adopted star formation threshold is appropriate for the resolution of the simulation, we show that the halo response is remarkably stable for $n\ge 5$, up to the highest star formation threshold that we test, $n=500$. This free parameter can be calibrated using the observed clustering of young stars. Simulations with low thresholds $n\le 1$ predict clustering that is too weak, while simulations with high star formation thresholds $n\ge 5$, are consistent with the observed clustering. Finally, we test the CDM predictions against the circular velocities of nearby dwarf galaxies. Low thresholds predict velocities that are too high, while simulations with $n\sim 10$ provide a good match to the observations. We thus conclude that the CDM model provides a good description of the structure of galaxies on kpc scales provided the effects of baryons are properly captured., 15 pages, 10 figures, published in MNRAS

We use hydrodynamical cosmological simulations to show that it is possible to create, via tidal interactions, galaxies lacking dark matter in a dark matter dominated universe. We select dwarf galaxies from the NIHAO project, obtained in the standard Cold Dark Matter model and use them as initial conditions for simulations of satellite-central interactions. After just one pericentric passage on an orbit with a strong radial component, NIHAO dwarf galaxies can lose up to 80 per~cent of their dark matter content, but, most interestingly, their central ($\approx 8$~kpc) dark matter to stellar ratio changes from a value of ${\sim}25$, as expected from numerical simulations and abundance matching techniques, to roughly unity as reported for NGC1052-DF2 and NGC1054-DF4. The stellar velocity dispersion drops from ${\sim}30$ ${\rm km\,s^{-1}}$ before infall to values as low as $6\pm 2$~ ${\rm km\,s^{-1}}$. These, and the half light radius around 3 kpc, are in good agreement with observations from van Dokkum and collaborators. Our study shows that it is possible to create a galaxy "without" dark matter starting from typical dwarf galaxies formed in a dark matter dominated universe, provided they live in a dense environment., 8 pages, 11 figures. Two extra figures added, accepted for publication in MNRAS

The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,\varphi)$ of luminous red giant stars within the mid-plane of the Milky Way reveal a spiral signature, which could plausibly reflect the response to a non-axisymmetric perturbation of the gravitational potential in the Galactic disk. We apply a simple steady-state toy model of a logarithmic spiral to interpret these observations, and find a good qualitative and quantitative match. Presuming that the amplitude of the gravitational potential perturbation is proportionate to that in the disk's surface mass density, we estimate the surface mass density amplitude to be $\Sigma_{\rm max} (R_{\odot})\approx 5.5\,\rm M_{\odot}\,pc^{-2}$ at the solar radius when choosing a fixed pattern speed of $\Omega_{\mathrm p}=12\,\rm km\,s^{-1}\,kpc^{-1}$. Combined with the local disk density, this implies a surface mass density contrast between the arm and inter-arm regions of approximately $\pm 10\%$ at the solar radius, with an increases towards larger radii. Our model constrains the pitch angle of the dynamical spiral arms to be approximately $12^{\circ}$., Comment: Accepted for publication in ApJ

The spin, $\lambda$, of dark matter (DM) halos in cosmological simulations follows a log normal distribution and has little correlation with galaxy observables. As such, there is currently no way to infer the $\lambda$ parameter of individual halos hosting observed galaxies. We present here a first attempt to measure $\lambda$ starting from the dynamically distinct stellar components identified in high-resolution cosmological simulations with Galactic Structure Finder. In a subsample of NIHAO galaxies, we find tight correlations between the total angular momentum (AM) of the DM halos, $J_h$, and the azimuthal AM, $J_z$, of the stellar components of the form: log($J_h$)=$\alpha$+$\beta\cdot$log($J_z$). The stellar halos have the tightest relation with $\alpha=9.50\pm0.42$ and $\beta=0.46\pm0.04$. The other tight relation is with the disks: $\alpha=6.15\pm0.92$ and $\beta=0.68\pm0.07$. We used Gaia DR2 and APOGEE to generate a combined kinematics-abundance space, where the Galaxy's thin and thick stellar disks stars can be neatly separated and their rotational velocity profiles, $v_{\phi}(R)$, can be computed. For both disks, $v_{\phi}(R)$ decreases with radius with $\sim$2 km s$^{-1}$ kpc$^{-1}$ for $R\gtrsim5$ kpc, resulting in $v_{\phi,thin}\backsimeq221$ km s$^{-1}$ and $v_{\phi,thick}\backsimeq188$ km s$^{-1}$ at $R_{\odot}$. These velocity profiles together with the Galaxy mass model of Cautun et al. (2020) result in the AM for the two disks: $J_{z,thin}=(3.26\pm0.43)\times10^{13}$ and $J_{z,thick}=(1.20\pm0.30)\times10^{13}$ M$_{\odot}$ kpc km s$^{-1}$, where the DM halo is assumed to have a contracted NFW profile. Adopting the correlation found in simulations, the spin estimate of the Galaxy's DM halo is $\lambda_{MW}=0.061^{+0.022}_{-0.016}$. If the DM halo has a NFW profile instead, the spin becomes $\lambda_{MW}=0.088^{+0.024}_{-0.020}$, making the Galaxy a more extreme outlier., Comment: 15 pages, 7 figures, accepted for publication in A&A

We use a sample of 83 high-resolution cosmological zoom-in simulations and a semi-analytic model (SAM) to study the stochasticity of galaxy formation in haloes ranging from dwarf to Milky Way masses. Our simulated galaxies reproduce the observed inefficiency of galaxy formation as expressed through the stellar, gas and baryonic Tully-Fisher relations. For HI velocities in the range (70 less than or similar to V less than or similar to 220 km s(-1)), the scatter is just 0.08 to 0.14 dex, consistent with the observed intrinsic scatter at these scales. At low velocities (20 less than or similar to V less than or similar to 70 km s(-1)), the simulated scatter is 0.2-0.25 dex, which could be tested with future observations. The scatter in the stellar mass versus dark halo velocity relation is constant for 30 less than or similar to V less than or similar to 180 km s(-1), and smaller (similar or equal to 0.17 dex) when using the maximum circular velocity of the dark-matter-only simulation, V-max(DMO), compared to the virial velocity (V-200 or V-200(DMO)). The scatter in stellar mass is correlated with halo concentration, and is minimized when using a circular velocity at a fixed fraction of the virial radius similar or equal to 0.4R(200) or with V-alpha = V-200(DMO)(V-max(DMO) / V-200(DMO))(alpha) with alpha similar or equal to 0.7, consistent with constraints from halo clustering. Using the SAM we show the correlation between halo formation time and concentration is essential in order to reproduce this result. This uniformity in galaxy formation efficiency we see in our hydrodynamical simulations and an SAM proves the simplicity and self-regulating nature of galaxy formation in a Lambda cold dark matter (Lambda CDM) universe.

In the paper 'NIHAO IV: Core creation and destruction in dark matter density profiles across cosmic time' published in MNRAS main journal, Volume 456, Issue 4, p. 3542-3552, an error in table 2 came to our attention. It appears that there is a mismatch between the values on the second line of table 2 and the corresponding curve in fig. 5. Indeed the parameters given for the fit of the relation between α and Mh are incorrect. We realised that the published parameters correspond to a different fitting formula that was used internally during thiswork and that does not have the right asymptotic behavior when the mass goes to zero. We apologize for this mistake. (Table presented) The correct fitting parameters to reproduce the fitting curve on fig. 5 are given in Table Errata. These parameters were provided only to allow the reader to reproduce our work and therefore this error does not affect our conclusions.

The direction of the spin vectors of disk galaxies change over time. We present the tilting rate of a sample of galaxies in the NIHAO suite of cosmological hydrodynamical simulations. Galaxies in our sample have been selected to be isolated and to have well determined spins. We compare the tilting rates to the predicted observing limit of Gaia, finding that our entire sample lies above the limit, in agreement with previous work. To test the role of dark matter and of gas we compare the weighted Pearson's correlation coefficients between the tilting rates and various properties. We find no correlation between the dark halo's tilting rate, shape, or misalignment with respect to the disc, and the tilting rate of the stellar disc. Therefore, we argue that, in the presence of gas, the dark halo plays a negligible role in the tilting of the stellar disc. On the other hand, we find a strong correlation between the tilting rate of the stellar disc and the misalignment of the cold gas warp. Adding the stellar mass fraction improves the correlation, while none of the dark matter's properties together with the cold gas misalignment improves the correlation to any significant extent. This implies that the gas cooling onto the disc is the principal driver of disc tilting., 12 pages, 13 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

We use 25 simulated galaxies from the NIHAO project to define and characterize a variety of kinematic stellar structures: thin and thick discs, large scale single discs, classical and pseudo bulges, spheroids, inner discs, and stellar haloes. These structures have masses, spins, shapes and rotational support in good agreement with theoretical expectations and observational data. Above a dark matter halo mass of $2.5\times10^{\rm~11}M_{\rm\odot}$, all galaxies have a classical bulge and 70\% have a thin and thick disc. The kinematic (thin) discs follow a power-law relation between angular momentum and stellar mass $J_{\rm *}=3.4M_{\rm *}^{\rm1.26\pm0.06}$, in very good agreement with the prediction based on the empirical stellar-to-halo mass relation in the same mass range, and show a strong correlation between maximum `observed' rotation velocity and dark matter halo circular velocity $v_{\rm c}=6.4v_{\rm max}^{0.64\pm0.04}$. Tracing back in time these structures' progenitors, we find all to lose a fraction $1-f_j$ of their maximum angular momentum. Thin discs are significantly better at retaining their high-redshift spins ($f_j\sim0.70$) than thick ones ($f_j\sim0.40$). Stellar haloes have their progenitor baryons assembled the latest ($z_{\rm~1/2}\sim1.1$) and over the longest timescales ($\tau\sim6.2$~Gyr), and have the smallest fraction of stars born in-situ ($f_{\rm in-situ}=0.35\pm0.14$). All other structures have $1.5\lesssim z_{\rm1/2}\lesssim3$, $\tau=4\pm2$~Gyr and $f_{\rm in-situ}\gtrsim0.9$., Comment: Accepted by MNRAS. First application of the methods described in arXiv:1804.05576

The stellar populations in the inner kiloparsecs of the Milky Way (MW) show complex kinematical and chemical structures. The origin and evolution of these structures is still under debate. Here we study the central region of a fully cosmological hydrodynamical simulation of a disk galaxy that reproduces key properties of the inner kiloparsecs of the MW: it has a boxy morphology and shows an overall rotation and dispersion profile in agreement with observations. We use a clustering algorithm on stellar kinematics to identify a number of discrete kinematic components: a high- and low-spin disk, a stellar halo and two bulge components; one fast rotating and one slow-rotating. We focus on the two bulge components and show that the slow rotating one is spherically symmetric while the fast rotating component shows a boxy/peanut morphology. Although the two bulge components are kinematically discrete populations at present-day, they are both mostly formed over similar time scales, from disk material. We find that stellar particles with lower initial birth angular momentum (most likely thick disc stars) end up in the slow-rotating low-spin bulge, while stars with higher birth angular momentum (most likely thin disc stars) are found in the high-spin bulge. This has the important consequence that a bulge population with a spheroidal morphology does not necessarily indicate a merger origin. In fact, we do find that only $\sim2.3$\% of the stars in the bulge components are ex-situ stars brought in by accreted dwarf galaxies early on. We identify these ex-situ stars as the oldest and most metal-poor stars on highly radial orbits with large vertical excursions from the disk., accepted by ApJ, 14 pages (plus 3 pages appendix), 13 + 8 figures

In this third paper of the series, we investigate the effects of warm dark matter with a particle mass of $m_\mathrm{WDM}=3\,\mathrm{keV}$ on the smallest galaxies in our Universe. We present a sample of 21 hydrodynamical cosmological simulations of dwarf galaxies and 20 simulations of satellite-host galaxy interaction that we performed both in a Cold Dark Matter (CDM) and Warm Dark Matter (WDM) scenario. In the WDM simulations, we observe a higher critical mass for the onset of star formation. Structure growth is delayed in WDM, as a result WDM haloes have a stellar population on average two Gyrs younger than their CDM counterparts. Nevertheless, despite this delayed star formation, CDM and WDM galaxies are both able to reproduce the observed scaling relations for velocity dispersion, stellar mass, size, and metallicity at $z=0$. WDM satellite haloes in a Milky Way mass host are more susceptible to tidal stripping due to their lower concentrations, but their galaxies can even survive longer than the CDM counterparts if they live in a dark matter halo with a steeper central slope. In agreement with our previous CDM satellite study we observe a steepening of the WDM satellites' central dark matter density slope due to stripping. The difference in the average stellar age for satellite galaxies, between CDM and WDM, could be used in the future for disentangling these two models., Comment: 10 pages, 11 figures, accepted for publication on MNRAS

A large fraction of the dwarf satellites orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and apparently kinematically coherent planar structure. Such a structure is not easily found in simulations based on the Cold Dark Matter model ({\Lambda}CDM). Using 21 high resolution cosmological simulations we analyse the kinematics of planes of satellites similar to the one around Andromeda. We find good agreement when co-rotation is characterized by the line-of-sight velocity. At the same time, when co-rotation is inferred by the angular momenta of the satellites, the planes are in agreement with the plane around our Galaxy. We find such planes to be common in our high concentration haloes. The number of co-rotating satellites obtained from the sign of the line-of-sight velocity shows large variations depending on the viewing angle and is consistent with that obtained from a sample with random velocities. We find that the clustering of angular momentum vectors of the satellites in the plane is a better measure of the kinematic coherence. Thus we conclude that the line-of- sight velocity is not well suited as a proxy for the kinematical coherence of the plane. Analysis of the kinematics of our planes shows a fraction of $\sim$30% chance aligned satellites. Tracking the satellites in the plane back in time reveals that these planes are a transient feature and not kinematically coherent as would appear at first sight. Thus we expect some of the satellites in the plane around Andromeda to have high velocities perpendicular to the plane., Comment: 19 pages, 15 figures, MNRAS accepted, Analysis strongly improved, paper shortened but main conclusions unchanged

We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo, Comment: 11 pages, 13 figures. Corrected typo in table 2 (middle row) with respect to the version published in MNRAS

We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the impact of the threshold for star formation on the response of the dark matter (DM) halo to baryonic processes. The fiducial NIHAO threshold, $n=10\, {\rm cm}^{-3}$, results in strong expansion of the DM halo in galaxies with stellar masses in the range $10^{7.5} < M_{star} < 10^{9.5} M_{\odot}$. We find that lower thresholds such as $n=0.1$ (as employed by the EAGLE/APOSTLE and Illustris/AURIGA projects) do not result in significant halo expansion at any mass scale. Halo expansion driven by supernova feedback requires significant fluctuations in the local gas fraction on sub-dynamical times (i.e., < 50 Myr at galaxy half-light radii), which are themselves caused by variability in the star formation rate. At one per cent of the virial radius, simulations with $n=10$ have gas fractions of $\simeq 0.2$ and variations of $\simeq 0.1$, while $n=0.1$ simulations have order of magnitude lower gas fractions and hence do not expand the halo. The observed DM circular velocities of nearby dwarf galaxies are inconsistent with CDM simulations with $n=0.1$ and $n=1$, but in reasonable agreement with $n=10$. Star formation rates are more variable for higher $n$, lower galaxy masses, and when star formation is measured on shorter time scales. For example, simulations with $n=10$ have up to 0.4 dex higher scatter in specific star formation rates than simulations with $n=0.1$. Thus observationally constraining the sub-grid model for star formation, and hence the nature of DM, should be possible in the near future., 18 pages, 13 figures, accepted to MNRAS

We used high resolution dark matter only cosmological simulations to investigate the structural properties of Lambda Cold Dark Matter ($\Lambda$CDM) haloes over cosmic time. The haloes in our study range in mass from $\sim 10^{10}$ to $\sim 10^{12} \mathrm{M}_\odot$, and are resolved with $10^5$ to $10^7$ particles. We fit the spherically averaged density profiles of DM haloes with the three parameter Einasto function. For our sample of haloes, the Einasto shape parameter, $\alpha$, is uncorrelated with the concentration, $c$, at fixed halo mass, and at all redshifts. Previous reports of an anti-correlation are traced to fitting degeneracies, which our fits are less sensitive to due to our higher spatial resolution. However, for individual haloes the evolution in $\alpha$ and $c$ is anti-correlated: at redshift $z=7$, $\alpha \simeq 0.4$ and decreases with time, while $c\simeq 3$ and increases with time. The evolution in structure is primarily due to accretion of mass at larger radii. We suggest that $\alpha$ traces the evolutionary state of the halo, with dynamically young haloes having high $\alpha$ (closer to a top-hat: $\alpha^{-1}=0$), and dynamically relaxed haloes having low $\alpha$ (closer to isothermal: $\alpha=0$). Such an evolutionary dependence reconciles the increase of $\alpha$ vs peak height, $\nu$, with the dependence on the slope of the power spectrum of initial density fluctuations found by previous studies., Comment: accepted for publication in MNRAS

We study the impact of the host on dwarf galaxy properties using four new Milky Way-like, ultra high-resolution simulations, ($N_{\rm part} >10^7$) from the NIHAO project. We split our sample in satellite ($R2.5 R_{\rm 200}$) galaxies. Simulated galaxies from all three groups are in excellent agreement with Local Group dwarf galaxies in terms of: stellar mass-velocity dispersion, stellar mass-metallicity relation, star formation histories, and stellar mass functions. Satellites and nearby galaxies show lower velocity dispersions and gas fractions compared to field galaxies. While field galaxies follow global abundance matching relations, satellites and nearby galaxies deviate from them, showing lower dark matter masses for given stellar mass. The reason for this deficit in dark matter mass is substantial mass loss experienced by satellites and $\sim80$\% of the nearby galaxies, while orbiting inside $R_{200}$ at earlier times. However, both satellites and nearby objects fall back onto the relation for field galaxies if we use the maximum of their virial mass instead of the present-day value. This allows us to provide estimates for the peak masses of observed Local Group galaxies. Finally, using radial velocities, distances, and the velocity dispersion-stellar mass relation from our simulations, we derive a metric to distinguish between galaxies harassed by the central object and unaffected ones. Applying this metric to observed objects we find that even far away dwarf galaxies like Eri II ($D\approx$ 370 kpc) have a strong probability ($\approx 83$\%) of having been affected by the Milky Way in the past. This naturally explains the lack of gas and recent star formation seen in Eri II., Comment: accepted by MNRAS, 31 pages, 17 pages main text, 14 pages appendix with tables of dwarf galaxy properties and convergence tests, 17 main figures

High resolution cosmological and hydrodynamical simulations have reached a resolution able to resolve in a self consistent way the disc of our galaxy, the galaxy center and the satellites orbiting around it. We present first results from the NIHAO-UHD project, a set of very high-resolution baryonic zoom-in simulations of Milky Way mass disc galaxies. These simulations model the full cosmological assembly history of the galaxies and their satellite system using the same, well tested physics as the NIHAO project. We show that these simulations can self-consistently reproduce the observed kinematical and morphological features of the X-shaped bulge observed in our own Milky Way.

In an effort to search for Ly α emission from circum- and intergalactic gas on scales of hundreds of kpc around z ̃ 3 quasars, and thus characterize the physical properties of the gas in emission, we have initiated an extensive fast survey with the Multi-Unit Spectroscopic Explorer (MUSE): Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission (QSO MUSEUM). In this work, we report the discovery of an enormous Ly α nebula (ELAN) around the quasar SDSS J102009.99+104002.7 at z = 3.164, which we followed-up with deeper MUSE observations. This ELAN spans ̃297 projected kpc, has an average Ly α surface brightness SBLy α ̃ 6.04 × 10-18 erg s-1 cm-2 arcsec-2(within the 2σ isophote) and is associated with an additional four previously unknown embedded sources: two Ly α emitters and two faint active galactic nuclei (one type-1 and one type-2 quasar). By mapping at high significance, the line-of-sight velocity in the entirety of the observed structure, we unveiled a large-scale coherent rotation-like pattern spanning ̃300 km s-1 with a velocity dispersion of -1, which we interpret as a signature of the inspiraling accretion of substructures within the quasar's host halo. Future multiwavelength data will complement our MUSE observations and are definitely needed to fully characterize such a complex system. None the less, our observations reveal the potential of new sensitive integral-field spectrographs to characterize the dynamical state of diffuse gas on large scales in the young Universe, and thereby witness the assembly of galaxies.

We show for the first time, that a fully cosmological hydrodynamical simulation can reproduce key properties of the innermost region of the Milky Way. Our high resolution simulation matches the profile and kinematics of the Milky Way's boxy/peanut-shaped bulge, and hence we can use it to reconstruct and understand the bulge assembly. In particular, the age dependence of the X-shape morphology of the simulated bulge parallels the observed metallicity dependent split in the red clump stars of the inner Galaxy. We use this feature to derive an observational metric that allows us to quantify when the bulge formed from the disk. The metric we propose can be employed with upcoming survey data to constrain the age of the Milky Way bar. From the split in stellar counts we estimate the formation of the 4~kpc scale bar in the simulation to have happened $t^{\rm bar}_{\rm form}\sim8^{+2}_{-2}$ Gyr ago, in good agreement with conventional methods to measure bar formation in simulations. We test the prospects for observationally differentiating the stars that belong to the bulge/bar compared to the surrounding disk, and find that the inner disk and bulge are practically indistinguishable in both chemistry and ages., Comment: 14+2 pages, 12+2 figures, accepted by ApJ

We use cosmological hydrodynamical zoom-in simulations with the SPH code gasoline of four haloes of mass M_{200} \sim 10^{13}\Msun to study the response of the dark matter to elliptical galaxy formation. Our simulations include metallicity dependent gas cooling, star formation, and feedback from massive stars and supernovae, but not active galactic nuclei (AGN). At z=2 the progenitor galaxies have stellar to halo mass ratios consistent with halo abundance matching, assuming a Salpeter initial mass function. However by z=0 the standard runs suffer from the well known overcooling problem, overpredicting the stellar masses by a factor of > 4. To mimic a suppressive halo quenching scenario, in our forced quenching (FQ) simulations, cooling and star formation are switched off at z=2. The resulting z=0 galaxies have stellar masses, sizes and circular velocities close to what is observed. Relative to the control simulations, the dark matter haloes in the FQ simulations have contracted, with central dark matter density slopes d\log��/d\log r \sim -1.5, showing that dry merging alone is unable to fully reverse the contraction that occurs at z>2. Simulations in the literature with AGN feedback however, have found expansion or no net change in the dark matter halo. Thus the response of the dark matter halo to galaxy formation may provide a new test to distinguish between ejective and suppressive quenching mechanisms., 21 pages, 18 figures, accepted to MNRAS

The satellites of the Milky Way and Andromeda represent the smallest galaxies we can observe in our Universe. In this series of papers we aim to shed light on their formation and evolution using cosmological hydrodynamical simulations. In this first paper we focus on the galaxy properties before accretion, by simulating twenty seven haloes with masses between $5\times 10^8$ and $10^{10} M_\odot$. Out of this set nineteen haloes successfully form stars, while eight remain dark. The simulated galaxies match quite well present day observed scaling relations between stellar mass, size and metallicity, showing that such relations are in place before accretion. Our galaxies show a large variety of star formation histories, from extended star formation periods to single bursts. As in more massive galaxies, large star formation bursts are connected with major mergers events, which greatly contribute to the overall stellar mass build up. The intrinsic stochasticity of mergers induces a large scatter in the stellar mass halo mass relation, up to two orders of magnitude. Despite the bursty star formation history, on these mass scales baryons are very ineffective in modifying the dark matter profiles, and galaxies with a stellar mass below $\approx 10^6 M_\odot$ retain their cuspy central dark matter distribution, very similar to results from pure N-body simulations., Comment: 12 pages, 15 figures. Two new figures added, discussion improved, several new references, conclusions unchanged. Accepted to be published in MNRAS. the second paper can found at arXiv:1707.01102

In the first paper we presented 27 hydrodynamical cosmological simulations of galaxies with total masses between $5 \times 10^8$ and $10^{10}\,\mathrm{M}_\odot$. In this second paper we use a subset of these cosmological simulations as initial conditions (ICs) for more than forty hydrodynamical simulations of satellite and host galaxy interaction. Our cosmological ICs seem to suggest that galaxies on these mass scales have very little rotational support and are velocity dispersion ($\sigma$) dominated. Accretion and environmental effects increase the scatter in the galaxy scaling relations (e.g. size - velocity dispersion) in very good agreement with observations. Star formation is substantially quenched after accretion. Mass removal due to tidal forces has several effects: it creates a very flat stellar velocity dispersion profiles, and it reduces the dark matter content at all scales (even in the centre), which in turn lowers the stellar velocity on scales around 0.5 kpc even when the galaxy does not lose stellar mass. Satellites that start with a cored dark matter profile are more prone to either be destroyed or to end up in a very dark matter poor galaxy. Finally, we found that tidal effects always increase the "cuspyness" of the dark matter profile, even for haloes that infall with a core., Comment: 13 pages, 15 figures, accepted to be published in MNRAS (with minor changes to the published version). This is the second paper of a series of two. The first paper can be found here arXiv:1707.01106

We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(\Omega_b/\Omega_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes., Comment: 20 pages, 15 figures, accepted to MNRAS

Many massive star forming disc galaxies in the redshift range 3 to 0.5 are observed to have a clumpy morphology showing giant clumps of size $\sim$1 kpc and masses of about $10^7M_{\odot}$ to $10^{10} M_{\odot}$. The nature and fate of these giant clumps is still under debate. In this work we use 19 high-resolution simulations of disc galaxies from the NIHAO sample to study the formation and the evolution of clumps in the discs of high redshift galaxies. We use mock HST - CANDELS observations created with the radiative transfer code GRASIL-3D to carry out, for the first time, a quantitative comparison of the observed fraction of clumpy galaxies and its evolution with redshift with simulations. We find a good agreement between the observed clumpy fraction and the one of the NIHAO galaxies. We find that dust attenuation can suppress intrinsically bright clumps and enhance less luminous ones. In our galaxy sample we only find clumps in light (u-band) from young stars but not in stellar mass surface density maps. This means that the NIHAO sample does not show clumpy stellar discs but rather a clumpy light distribution originating from clumpy star formation events. The clumps found in the NIHAO sample match observed age/color gradients as a function of distance from the galaxy center but they show no sign of inward migration. Clumps in our simulations disperse on timescales of a about a hundred Myr and their contribution to bulge growth is negligible., Comment: 17 pages, 18 figures, accepted by MNRAS

Recently it has been shown that a large fraction of the dwarf satellite galaxies orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and kinematically coherent planar structure. The presence of such a structure seems to challenge the current Cold Dark Matter paradigm of structure formation, which predicts a more uniform distribution of satellites around central objects. We show that it is possible to obtain a thin, extended, rotating plane of satellites resembling the one in Andromeda in cosmological collisionless simulations based on the Cold Dark Matter model. Our new high resolution simulations show a correlation between the formation time of the dark matter halo and the thickness of the plane of satellites. Our simulations have a high incidence of satellite planes as thin, extended, and as rich as the one in Andromeda and with a very coherent kinematic structure when we select high concentration/early forming halos. By tracking the formation of the satellites in the plane we show that they have been mainly accreted onto the main object along thin dark matter filaments at high redshift. Our results show that the presence of a thin, extended, rotating plane of satellites is not a challenge for the Cold Dark Matter paradigm, but actually supports one of the predictions of this paradigm related to the presence of filaments of dark matter around galaxies at high redshift., 6 pages, 5 figures, accepted for publication in ApJ Letters, minor changes to match the accepted version

Recently it has been shown by Ibata et al. (2013) that a large fraction of the dwarf satellite galaxies found in the PAndAS survey (McConnachie et al. 2009) and orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended, and kinematically coherent planar structure. The presence of such a structure seems to challenge the current Cold Dark Matter paradigm of structure formation (Ibata et al. 2014, Pawlowski et al. 2014), which predicts a more uniform distribution of satellites around central objects. We show that it is possible to obtain a thin, extended, rotating plane of satellites resembling the one in Andromeda in cosmological collisionless simulations based on this model. Our new 21 high-resolution simulations (see Buck et al. 2015) show a correlation between the formation time of the dark matter halo and the thickness of the plane of satellites. Our simulations have a high incidence of satellite planes as thin, extended, and as rich as the one in Andromeda and with a very coherent kinematic structure when we select early forming haloes. By tracking the formation of the satellites in the plane we show that they have mainly been accreted onto the main object along thin dark matter filaments at high redshift (Dekel et al. 2009, Libeskind et al. 2009, 2011). Our results show that the presence of a thin, extended, rotating plane of satellites is not a challenge for the Cold Dark Matter paradigm, but actually supports one of the predictions of this paradigm related to the presence of filaments of dark matter around galaxies at high redshift.

Fairchild und Dornier sind uber die langste Zeit ihrer Unternehmensexistenz getrennte Wege ohne gegenseitige Beruhrung gegangen. Erst gegen Ende ihres Firmenlebens im Jahre 1996 finden beide Unternehmen als Fairchild Dornier zusammen. Ziel ist es, mit einem neuen Produkt in dem sehr zyklischen, technologie- und damit auch kapitalintensiven Flugzeugmarkt erfolgreich zu sein.

Das Aufgabengebiet der Kontrolle hat zwischenzeitlich nicht nur in der Betriebswirtschaftslehre einen hohen Stellenwert erlangt.1 Bis zum heutigen Zeitpunkt findet man zahlreiche kontrollrelevante Konzeptionen, die sich uber eine breite Palette wissenschaftlicher Disziplinen erstrecken. Entsprechend vielgestaltig sind auch die Losungsansatze zur Kontrollproblematik. Je nach Aufgabengebiet werden sogar innerhalb einzelner Disziplinen unterschiedliche Begriffsverstandnisse von Kontrolle verwendet. Fur die Entwicklung eines zu erarbeitenden funktionsspezifischen Kontrollkonzepts ist es aus diesem Grunde zweckmasig, zunachst einen theoriegeleiteten Bezugsrahmen zu erarbeiten, der auch interdisziplinare Gesichtspunkte in die Betrachtung integriert. Dieser soll anschliesend eine Ableitung des hier zugrunde gelegten Kontrollbegriffs erlauben.

Bevor das Konzept einer integrierten Beschaffungskontrolle weiter prazisiert wird, soll zunachst der Entwicklungsstand der Beschaffungskontrolle in der deutschsprachigen Literatur untersucht werden.1 Zur Darstellung des heute als’ state of the art’ zu bezeichnenden Kontroll- bzw. Controlling-Know-how in der Beschaffungswirtschaft werden einige Forschungsarbeiten, die konzeptionelle Bedeutung beanspruchen konnen, in exemplarischer Weise ausgewertet.2

Die steigende Ambiguitat des Wettbewerbsumfeldes erfordert von der Betriebswirtschaftslehre eine zunehmend genauere Ausdifferenzierung der umfassenden Fuhrungsaufgaben. Um diese Ausdifferenzierung zu vollziehen, findet man in der Literatur diesbezuglich zwei grundsatzliche Systematisierungen: Die Unterscheidung nach Realfunktionen strukturiert die Prozesse der Leistungserstellung (Beschaffung, Produktion, Absatz), wahrend die Managementfunktionen auf den Prozes der Unternehmensfuhrung (Planung, Realisation und Kontrolle) rekurrieren.1 Die vorliegende Arbeit verbindet mit Beschaffung und Kontrolle jeweils ein Element der genannten Teilfunktionen. Sie zielt auf eine sachlich angemessene Betrachtung der Kontrolle im Beschaffungsmanagement unter integrativen Gesichtspunkten ab. Folgende Uberlegungen sollen die Wahl dieses Themenschwerpunktes kurz erlautern: (1) In den letzten zehn Jahren hat das Beschaffungsmanagement sowohl in der Literatur als auch in der Praxis eine starke Aufwertung erfahren. Im Zusammenhang mit veranderten marktlichen Rahmenbedingungen, der Notwendigkeit einer Konzentration auf Kernkompetenzen und der damit einhergehenden Reduzierung der Fertigungstiefe mus die Beschaffung zunehmend in die strategische Unternehmensfuhrung eingebunden werden. Die Beschaffung ist nunmehr nicht langer eine operativ ausgerichtete Ausfuhrungsfunktion bereits vorgedachter strategischer Uberlegungen. Strategische beschaffungsbezogene Handlungspotentiale fuhren zu einer Verpflichtung des Management, die Beschaffung als Quelle zur Erzielung strategischer Wettbewerbsvorteile fur das gesamte Unternehmen anzuerkennen und entsprechend auszugestalten. (2) Ahnlich wie die Beschaffung innerhalb der Realfunktionen der Betriebswirtschaftslehre reprasentiert auch die Kontrolle eine — insbesondere im Vergleich zur betriebswirtschaftlichen Planung — traditionell vernachlassigte Fuhrungsfunktion innerhalb der strategischen Managementlehre. Die Kontrolle wurde uber Jahrzehnte hinweg lediglich als Teilfunktion von uberwiegend untergeordnetem Stellenwert angesehen, der es als letzter Phase im Managementprozes oblag, die Erreichung von Handlungszielen durch Soll-lst-Vergleiche ex post festzustellen.2 In der deutschen Betriebswirtschaftslehre erfuhr die Kontrollfunktion erst seit Mitte der 80er Jahre insbesondere durch die Leistungen der Forschungsgruppe um Steinmann/Schreyogg einen bemerkenswerten Bedeutungsgewinn. Inhalte dieser Ansatze werden auch in dieser Arbeit als wichtige Bausteine einer integrierten Beschaffungskontrolle wiederzufinden sein. Kontrolle wird diesem Verstandnis folgend zu einer Aufgabe, die gleichgewichtig gegenuber anderen Managementfunktionen — vor allem gegenuber der Planung — einzuordnen ist. So besehen wandelt sich das Verstandnis von Kontrolle von einer reinen Komplementar- bzw. Korrekturfunktion zu einer praventiven Kompensationsfunktion mit Umsteuerungspotential.1