Stellar Metallicity Distribution Research Articles

SN1991bg-like supernovae are associated with old stellar populations

Abstract SN1991bg-like supernovae are a distinct subclass of thermonuclear Type Ia supernovae (SNe Ia). Their spectral and photometric peculiarities indicate that their progenitors and explosion mechanisms differ from ‘normal’ SNe Ia. One method of determining information about supernova progenitors we cannot directly observe is to observe the stellar population adjacent to the apparent supernova explosion site to infer the distribution of stellar population ages and metallicities. We obtain integral field observations and analyse the spectra extracted from regions of projected radius $\sim\,\text{kpc}$ about the apparent SN explosion site for 11 91bg-like SNe in both early- and late-type galaxies. We utilise full-spectrum spectral fitting to determine the ages and metallicities of the stellar population within the aperture. We find that the majority of the stellar populations that hosted 91bg-like supernovae have little recent star formation. The ages of the stellar populations suggest that that 91bg-like SN progenitors explode after delay times of >6 Gyr, much longer than the typical delay time of normal SNe Ia, which peaks at $\sim$ 1 Gyr.

Barlenses in the CALIFA survey: Combining photometric and stellar population analyses

Aims: It is theoretically predicted that, at low galaxy inclinations, boxy/peanut bar components have a barlens appearance of a round central component embedded in the narrow bar. We investigate barlenses in the Calar Alto Legacy Integral Field Area (CALIFA) survey galaxies, studying their morphologies, stellar populations, and metallicities. We show that, when present, barlenses account for a significant portion of light of photometric bulges, i.e., the excess light on top of the disks, which highlights the importance of bars in accumulating central galaxy mass concentrations in the cosmic timescale. Methods: We made multi-component decompositions for a sample of 46 barlens galaxies drawn from the CALIFA survey, where M⋆/M⊙ = 109.7 − 1011.4 and z = 0.005 − 0.03. Unsharp masks of the Sloan Digital Sky Survey (SDSS) r′-band mosaics were used to identify the boxy/peanut or X-shaped features. Barlenses are identified in the images using our simulation snapshots as an additional guide. Our decompositions with GALFIT include bulges, disks, and bars as well as barlenses as a separate component. For 26 of the decomposed galaxies the CALIFA DR2 V500 grating data cubes were used to explore stellar ages and metallicities at the regions of various structure components. Results: We find that 25 ± 2% of the 1064 galaxies in the whole CALIFA sample show either X-shaped or barlens features. In the decomposed galaxies with barlenses, on average 13% ± 2% of the total galaxy light belongs to this component, leaving less than 10% for possible separate bulge components. Most importantly, bars and barlenses are found to have similar cumulative stellar age and metallicity distributions. The metallicities in barlenses are on average near solar, but exhibit a large range. In some of the galaxies barlenses and X-shaped features appear simultaneously, in which case the bar origin of the barlens is unambiguous. Conclusion: This is the first time that a combined morphological and stellar population analysis is used to study barlenses. We show that their stars are accumulated in a prolonged time period concurrently with the evolution of the narrow bar.

Pushing back the limits: detailed properties of dwarf galaxies in a ΛCDM universe

We present the results of a set of high-resolution chemo-dynamical simulations of dwarf galaxies in a ΛCDM cosmology. Out of an original (3.4 Mpc/h)3 cosmological box, a sample of 27 systems are re-simulated from z = 70 to z = 0 using a zoom-in technique. Gas and stellar properties are confronted to the observations in the greatest details: in addition to the galaxy global properties, we investigated the model galaxy velocity dispersion profiles, half-light radii, star formation histories, stellar metallicity distributions, and [Mg/Fe] abundance ratios. The formation and sustainability of the metallicity gradients and kinematically distinct stellar populations are also tackled. We show how the properties of six Local Group dwarf galaxies, NGC 6622, Andromeda II, Sculptor, Sextans, Ursa Minor and Draco are reproduced, and how they pertain to three main galaxy build-up modes. Our results indicate that the interaction with a massive central galaxy could be needed for a handful of Local Group dwarf spheroidal galaxies only, the vast majority of the systems and their variety of star formation histories arising naturally from a ΛCDM framework. We find that models fitting well the local Group dwarf galaxies are embedded in dark haloes of mass between 5 × 108 to a few 109 M⊙, without any missing satellite problem. We confirm the failure of the abundance matching approach at the mass scale of dwarf galaxies. Some of the observed faint however gas-rich galaxies with residual star formation, such as Leo T and Leo P, remain challenging. They point out the need of a better understanding of the UV-background heating.

The VMC survey – XXXI: The spatially resolved star formation history of the main body of the Small Magellanic Cloud

We recover the spatially resolved star formation history across the entire main body and Wing of the Small Magellanic Cloud (SMC), using fourteen deep tile images from the VISTA survey of the Magellanic Clouds (VMC), in the YJKs filters. The analysis is performed on 168 subregions of size 0.143 deg2, covering a total contiguous area of 23.57 deg2. We apply a colour-magnitude diagram (CMD) reconstruction method that returns the best-fitting star formation rate SFR(t), age--metallicity relation, distance and mean reddening, together with their confidence intervals, for each subregion. With respect to previous analyses, we use a far larger set of VMC data, updated stellar models, and fit the two available CMDs (Y-Ks versus Ks and J-Ks versus Ks) independently. The results allow us to derive a more complete and more reliable picture of how the mean distances, extinction values, star formation rate, and metallicities vary across the SMC, and provide a better description of the populations that form its Bar and Wing. We conclude that the SMC has formed a total mass of (5.31+-0.05)x10^8 Msun in stars over its lifetime. About two thirds of this mass is expected to be still locked in stars and stellar remnants. 50 per cent of the mass was formed prior to an age of 6.3 Gyr, and 80 per cent was formed between 8 and 3.5 Gyr ago. We also illustrate the likely distribution of stellar ages and metallicities in different parts of the CMD, to aid the interpretation of data from future astrometric and spectroscopic surveys of the SMC.

No Metallicity Correlation Associated with the Kepler Dichotomy

Abstract NASA’s Kepler mission has discovered thousands of planetary systems, ∼20% of which are found to host multiple transiting planets. This relative paucity (compared to the high fraction of single transiting systems) is postulated to result from a distinction in the architecture between multi-transiting systems and those hosting a single transiting planet: a phenomenon usually referred to as the Kepler dichotomy. In this paper, we investigate the hypothesis that external giant planets are the main cause behind the over-abundance of single- relative to multi-transiting systems, which would be signaled by higher metallicities in the former sample. To this end, we perform a statistical analysis on the stellar metallicity distribution with respect to planet multiplicity in the Kepler data. We perform our analysis on a variety of samples taken from a population of 1166 Kepler main-sequence planetary hosts, using precisely determined metallicities from the California-Kepler survey and Swift et al. Contrary to some predictions, we do not find a significant difference between the stellar metallicities of the single- and multiple-transiting planet systems. However, we do find a 55% upper bound for systems with a single non-giant planet that could also host a hidden giant planet, based on metallicity considerations. While the presence of external giant planets might be one factor behind the Kepler dichotomy, our results also favor alternative explanations. We suggest that additional radial velocity and direct imaging measurements are necessary to constrain the presence of gas giants in systems with a single transiting planet.

Modelling chemical abundance distributions for dwarf galaxies in the Local Group: the impact of turbulent metal diffusion

We investigate stellar metallicity distribution functions (MDFs), including Fe and ${\alpha}$-element abundances, in dwarf galaxies from the Feedback in Realistic Environments (FIRE) project. We examine both isolated dwarf galaxies and those that are satellites of a Milky Way-mass galaxy. In particular, we study the effects of including a sub-grid turbulent model for the diffusion of metals in gas. Simulations that include diffusion have narrower MDFs and abundance ratio distributions, because diffusion drives individual gas and star particles toward the average metallicity. This effect provides significantly better agreement with observed abundance distributions of dwarf galaxies in the Local Group, including the small intrinsic scatter in [${\alpha}$/Fe] vs. [Fe/H] (less than 0.1 dex). This small intrinsic scatter arises in our simulations because the interstellar medium (ISM) in dwarf galaxies is well-mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to within 0.1 dex. Thus, most of the scatter in abundances at z = 0 arises from redshift evolution and not from instantaneous scatter in the ISM. We find similar MDF widths and intrinsic scatter for satellite and isolated dwarf galaxies, which suggests that environmental effects play a minor role compared with internal chemical evolution in our simulations. Overall, with the inclusion of metal diffusion, our simulations reproduce abundance distribution widths of observed low-mass galaxies, enabling detailed studies of chemical evolution in galaxy formation.

Stellar Populations in a semi-analytic model I: Bulges of Milky Way-like galaxies

We study the stellar populations of bulges of Milky Way-like (MW-like) galaxies with the aim of identifying the physical processes involved in the formation of the bulge of our Galaxy. We use the semi-analytic model of galaxy formation and evolution SAG adapted to this aim; this kind of models can trace the properties of galaxies and their components like stellar discs, bulges and halos, but resolution limits prevent them from reaching the scale of stellar populations (SPs). Properties of groups of stars formed during single star formation events are stored and tracked in the model and results are compared with observations of stars in the galactic bulge. MW-like galaxies are selected using two different criteria. One of them considers intrinsic photo-metric properties and the second is focused on the cosmological context of the local group of galaxies (LG). We compare our model results with spectroscopic and photometric stellar metallicity distributions. We find that 87% of stars in bulges of MWtype galaxies in our model are accreted and formed in starbursts during disc instability events. Mergers contribute to 13% of the mass budget of the bulge and are responsible for the low metallicity tail of the distribution. Abundance ratios of {\alpha} elements with respect to iron, [{\alpha}/Fe], are measured in SPs of model galaxies. The patterns found in the model for SPs with different origins help to explain the lack of a gradient of [{\alpha}/Fe] ratios in observed stars along the minor axis of the bulge.

Observational Evidence for High Neutronization in Supernova Remnants: Implications for Type Ia Supernova Progenitors

Abstract The physical process whereby a carbon–oxygen white dwarf explodes as a Type Ia supernova (SN Ia) remains highly uncertain. The degree of neutronization in SN Ia ejecta holds clues to this process because it depends on the mass and the metallicity of the stellar progenitor, and on the thermodynamic history prior to the explosion. We report on a new method to determine ejecta neutronization using Ca and S lines in the X-ray spectra of Type Ia supernova remnants (SNRs). Applying this method to Suzaku data of Tycho, Kepler, 3C 397, and G337.2−0.7 in the Milky Way, and N103B in the Large Magellanic Cloud, we find that the neutronization of the ejecta in N103B is comparable to that of Tycho and Kepler, which suggests that progenitor metallicity is not the only source of neutronization in SNe Ia. We then use a grid of SN Ia explosion models to infer the metallicities of the stellar progenitors of our SNRs. The implied metallicities of 3C 397, G337.2−0.7, and N103B are major outliers compared to the local stellar metallicity distribution functions, indicating that progenitor metallicity can be ruled out as the origin of neutronization for these SNRs. Although the relationship between ejecta neutronization and equivalent progenitor metallicity is subject to uncertainties stemming from the 12C + 16O reaction rate, which affects the Ca/S mass ratio, our main results are not sensitive to these details.

The Stellar Metallicity Distribution of the Galactic Halo Based on SCUSS and SDSS Data

Abstract Based on the Sloan Digital Sky Survey and South Galactic Cap u-band Sky Survey (SCUSS), we simulate the photometric metallicity distribution functions (MDFs) of stars in the Galactic halo. The photometric metallicity of stars was estimated by a new Monte-Carlo method. Due to the use of a more reliable metallicity calibration method and more accurate u-band deep measurements from SCUSS, we can obtain more accurate MDFs of a large sample of distant stars in the Galactic halo. In this study, we select 78,092 F/G main-sequence turnoff stars (MSTO) in the south Galactic cap, with 0.2 < (g − r)0 < 0.4, as tracers of the stellar MDFs in the Galactic halo. The sample stars are divided into two height intervals above the Galactic plane: −8 < z < −4 kpc and −12 < z < −8 kpc. The MDFs of selected stars in each interval are well fit by a three-Gaussian model, with peaks at [Fe/H] ≈ −0.63, −1.45, and −2.0. The two metal-poor components correspond to the inner halo and outer halo, respectively. The fraction of the metal-rich component, which may be contributed by the substructure (such as Sagittarius stream or other streams) is about 10%. With limited kinematic estimation, we find the correlations between metallicity and kinematics. Our results provide additional supporting evidence of duality of the Galactic halo.

Chemical enrichment of stars due to accretion from the ISM during the Galaxy's assembly

Using the Eris zoom-in cosmological simulation of assembly of a Milky Way analog, we study chemical enrichment of stars due to accretion of metal-enriched gas from the interstellar medium during the Galaxy's development. We consider metal-poor and old stars in both Galactic halo and bulge and make use of stellar orbits, gas density and metallicity distributions in Eris. Assuming spherically symmetric Bondi-Hoyle accretion, we find that halo and bulge stars accrete metals at the rate of about 10^-24 solar mass per year and 10^-22 solar mass per year, respectively, at redshifts z < 3, but this accretion rate increases hundred-fold to about 10^-20 solar mass per year at higher redshifts due to increased gas density. Bulge and halo stars accrete similar amounts of metals at high redshifts as kinematically distinct bulge and halo are not yet developed at these redshifts and both sets of stars encounter similar metal distribution in the ISM on average. Accretion alone can enrich main-sequence stars up to [Fe/H] -2 in extreme cases.Median enrichment level due to accretion in these stars is about [Fe/H]~-6 to -5.Because accretion mostly takes place at high redshifts, it is alpha-enriched to [alpha/Fe]~0.5. We find that accretive metal enrichment is significant enough to affect the predicted metallicity distribution function of halo stars at [Fe/H] < -5.This suggests that attempts to infer the natal chemical environment of the most metal-poor stars from their observed enrichment today can be hindered due to metal accretion. Peculiar enrichment patterns such as those predicted to arise from pair-instability supernovae could help in disentangling natal and accreted metal content of stars.

Crater 2: An Extremely Cold Dark Matter Halo

Abstract We present results from MMT/Hectochelle spectroscopy of 390 red giant candidate stars along the line of sight to the recently discovered Galactic satellite Crater 2. Modeling the joint distribution of stellar positions, velocities, and metallicities as a mixture of Crater 2 and Galactic foreground populations, we identify ∼62 members of Crater 2, for which we resolve a line-of-sight velocity dispersion of = km s−1 and a mean velocity of = km s−1 (solar rest frame). We also resolve a metallicity dispersion of = dex and a mean of = dex that is 0.28 ± 0.14 dex poorer than estimated from photometry. Despite Crater 2's relatively large size (projected halflight radius R h ∼ 1 kpc) and intermediate luminosity (M V ∼ −8), its velocity dispersion is the coldest that has been resolved for any dwarf galaxy. These properties make Crater 2 the most extreme low-density outlier in dynamical as well as structural scaling relations among the Milky Way’s dwarf spheroidals. Even so, under assumptions of dynamical equilibrium and negligible contamination by unresolved binary stars, the observed velocity distribution implies a gravitationally dominant dark matter halo, with a dynamical mass of M ⊙ and a mass-to-light ratio of enclosed within a radius of ∼1 kpc, where the equivalent circular velocity is km s−1.

The Metallicity Distribution and Hot Jupiter Rate of the Kepler Field: Hectochelle High-resolution Spectroscopy for 776 Kepler Target Stars

Abstract The occurrence rate of hot Jupiters from the Kepler transit survey is roughly half that of radial velocity surveys targeting solar neighborhood stars. One hypothesis to explain this difference is that the two surveys target stars with different stellar metallicity distributions. To test this hypothesis, we measure the metallicity distribution of the Kepler targets using the Hectochelle multi-fiber, high-resolution spectrograph. Limiting our spectroscopic analysis to 610 dwarf stars in our sample with &gt; 3.5, we measure a metallicity distribution characterized by a mean of , in agreement with previous studies of the Kepler field target stars. In comparison, the metallicity distribution of the California Planet Search radial velocity sample has a mean of , and the samples come from different parent populations according to a Kolmogorov–Smirnov test. We refit the exponential relation between the fraction of stars hosting a close-in giant planet and the host star metallicity using a sample of dwarf stars from the California Planet Search with updated metallicities. The best-fit relation tells us that the difference in metallicity between the two samples is insufficient to explain the discrepant hot Jupiter occurrence rates; the metallicity difference would need to be ≃0.2–0.3 dex for perfect agreement. We also show that (sub)giant contamination in the Kepler sample cannot reconcile the two occurrence calculations. We conclude that other factors, such as binary contamination and imperfect stellar properties, must also be at play.

SDSS-IV MaNGA: environmental dependence of stellar age and metallicity gradients in nearby galaxies

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

A MONTE CARLO METHOD FOR MAKING THE SDSS u-BAND MAGNITUDE MORE ACCURATE

ABSTRACT We develop a new Monte Carlo-based method to convert the Sloan Digital Sky Survey (SDSS) u-band magnitude to the south Galactic Cap of the u-band Sky Survey (SCUSS) u-band magnitude. Due to the increased accuracy of SCUSS u-band measurements, the converted u-band magnitude becomes more accurate compared with the original SDSS u-band magnitude, in particular at the faint end. The average u-magnitude error (for both SDSS and SCUSS) of numerous main-sequence stars with increases as the g-band magnitude becomes fainter. When g = 19.5, the average magnitude error of the SDSS u is 0.11. When g = 20.5, the average SDSS u error rises to 0.22. However, at this magnitude, the average magnitude error of the SCUSS u is just half as much as that of the SDSS u. The SDSS u-band magnitudes of main-sequence stars with and are converted, therefore the maximum average error of the converted u-band magnitudes is 0.11. The potential application of this conversion is to derive a more accurate photometric metallicity calibration from SDSS observations, especially for the more distant stars. Thus, we can explore stellar metallicity distributions either in the Galactic halo or some stream stars.

The SLUGGS Survey: A New Mask Design to Reconstruct the Stellar Populations and Kinematics of Both Inner and Outer Galaxy Regions

AbstractIntegral 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.

KINETIC ENERGY FROM SUPERNOVA FEEDBACK IN HIGH-RESOLUTION GALAXY SIMULATIONS

We describe a new method for adding a prescribed amount of kinetic energy to simulated gas modeled on a cartesian grid by directly altering grid cells' mass and velocity in a distributed fashion. The method is explored in the context of supernova feedback in high-resolution ($\sim 10$ pc) hydrodynamic simulations of galaxy formation. Resolution-dependence is a primary consideration in our application of the method and simulations of isolated explosions (performed at different resolutions) motivate a resolution-dependent scaling for the injected fraction of kinetic energy that we apply in cosmological simulations of a $10^9$ Msun dwarf halo. We find that in high density media ($\gtrsim$ 50 cm$^{-3}$) with coarse resolution ($\gtrsim 4$ pc per cell), results are sensitive to the initial kinetic energy fraction due to early and rapid cooling. In our galaxy simulations, the deposition of small amounts of supernova energy in kinetic form (as little as 1%) has a dramatic impact on the evolution of the system, resulting in an order of magnitude suppression of stellar mass. The overall behavior of the galaxy in the two highest resolution simulations we perform appears to converge. We discuss the resulting distribution of stellar metallicities, an observable sensitive to galactic wind properties, and find that while the new method demonstrates increased agreement with observed systems, significant discrepancies remain, likely due to simplistic assumptions that neglect contributions from Type Ia supernovae and stellar winds.

Impact of radial migration on stellar and gas radial metallicity distribution

Radial migration is defined as the change in guiding centre radius of stars and gas caused by gains or losses of angular momentum that result from gravitational interaction with non-axisymmetric structure. This has been shown to have significant impact on the metallicity distribution in galactic discs, and therefore affects the interpretation of Galactic archeology. We use a simulation of a Milky Way-sized galaxy to examine the effect of radial migration on the star and gas radial metallicity distribution. We find that both the star and gas component show significant radial migration. The stellar radial metallicity gradient remains almost unchanged but the radial metallicity distribution of the stars is broadened to produce a greater dispersion at all radii. However, the metallicity dispersion of the gas remains narrow. We find that the main drivers of the gas metallicity distribution evolution are metal enrichment and mixing: more efficient metal enrichment in the inner region maintains a negative slope in the radial metallicity distribution, and the metal mixing ensures the tight relationship of the gas metallicity with the radius. The metallicity distribution function reproduces the trend in the age-metallicity relation found from observations for stars younger than 1.0 Gyr in the Milky Way.

Synthetic galaxy images and spectra from the Illustris simulation

We present our methods for generating a catalog of 7,000 synthetic images and 40,000 integrated spectra of redshift z = 0 galaxies from the Illustris Simulation. The mock data products are produced by using stellar population synthesis models to assign spectral energy distributions (SED) to each star particle in the galaxies. The resulting synthetic images and integrated SEDs therefore properly reflect the spatial distribution, stellar metallicity distribution, and star formation history of the galaxies. From the synthetic data products it is possible to produce monochromatic or color-composite images, perform SED fitting, classify morphology, determine galaxy structural properties, and evaluate the impacts of galaxy viewing angle. The main contribution of this paper is to describe the production, format, and composition of the image catalog that makes up the Illustris Simulation Obsevatory. As a demonstration of this resource, we derive galactic stellar mass estimates by applying the SED fitting code FAST to the synthetic galaxy products, and compare the derived stellar masses against the true stellar masses from the simulation. We find from this idealized experiment that systematic biases exist in the photometrically derived stellar mass values that can be reduced by using a fixed metallicity in conjunction with a minimum galaxy age restriction.

Stellar population gradients in galaxy discs from the CALIFA survey

While studies of gas-phase metallicity gradients in disc galaxies are common, very little has been done in the acquisition of stellar abundance gradients in the same regions. We present here a comparative study of the stellar metallicity and age distributions in a sample of 62 nearly face-on, spiral galaxies with and without bars, using data from the CALIFA survey. We measure the slopes of the gradients and study their relation with other properties of the galaxies. We find that the mean stellar age and metallicity gradients in the disc are shallow and negative. Furthermore, when normalized to the effective radius of the disc, the slope of the stellar population gradients does not correlate with the mass or with the morphological type of the galaxies. Contrary to this, the values of both age and metallicity at $\sim$2.5 scale-lengths correlate with the central velocity dispersion in a similar manner to the central values of the bulges, although bulges show, on average, older ages and higher metallicities than the discs. One of the goals of the present paper is to test the theoretical prediction that non-linear coupling between the bar and the spiral arms is an efficient mechanism for producing radial migrations across significant distances within discs. The process of radial migration should flatten the stellar metallicity gradient with time and, therefore, we would expect flatter stellar metallicity gradients in barred galaxies. However, we do not find any difference in the metallicity or age gradients in galaxies with without bars. We discuss possible scenarios that can lead to this absence of difference.

Deep SDSS optical spectroscopy of distant halo stars

We analyze a sample of tens of thousands of spectra of halo turnoff stars, obtained with the optical spectrographs of the Sloan Digital Sky Survey (SDSS), to characterize the stellar halo population "in situ" out to a distance of a few tens of kpc from the Sun. In this paper we describe the derivation of atmospheric parameters. We also derive the overall stellar metallicity distribution based on F-type stars observed as flux calibrators for the Baryonic Oscillations Spectroscopic Survey (BOSS). Our analysis is based on an automated method that determines the set of parameters of a model atmosphere that reproduces each observed spectrum best. We used an optimization algorithm and evaluate model fluxes by means of interpolation in a precomputed grid. In our analysis, we account for the spectrograph's varying resolution as a function of fiber and wavelength. Our results for early SDSS (pre-BOSS upgrade) data compare well with those from the SEGUE Stellar Parameter Pipeline (SSPP), except for stars with logg (cgs units) lower than 2.5. An analysis of stars in the globular cluster M13 reveals a dependence of the inferred metallicity on surface gravity for stars with logg < 2.5, confirming the systematics identified in the comparison with the SSPP. We find that our metallicity estimates are significantly more precise than the SSPP results. We obtain a halo metallicity distribution that is narrower and more asymmetric than in previous studies. The lowest gravity stars in our sample, at tens of kpc from the Sun, indicate a shift of the metallicity distribution to lower abundances, consistent with what is expected from a dual halo system in the Milky Way.