Nuclear and extended spectra of NGC 1068-II. Near-infrared stellar population synthesis

We performed stellar population synthesis on the nuclear and extended regions of NGC 1068 by means of near-infrared spectroscopy to disentangle their spectral energy distribution components. This is the first time that such a technique is applied to the whole 0.8-2.4 mu m wavelength interval in this galaxy. NGC 1068 is one of the nearest and probably the most studied Seyfert 2 galaxy, becoming an excellent laboratory to study the interaction between black holes, the jets that they can produce and the medium in which they propagate. Our main result is that traces of young stellar population are found at similar to 100 pc south of the nucleus. The contribution of a power-law continuum in the centre is about 25 per cent, which is expected if the light is scattered from a Seyfert 1 nucleus. We find peaks in the contribution of the featureless continuum about 100-150 pc from the nucleus on both sides. They might be associated with regions where the jet encounters dense clouds. Further support to this scenario is given by the peaks of hot dust distribution found around these same regions and the H(2) emission-line profile, leading us to propose that the peaks might be associated to regions where stars are being formed. Hot dust also has an important contribution to the nuclear region, reinforcing the idea of the presence of a dense, circumnuclear torus in this galaxy. Cold dust appears mostly in the south direction, which supports the view that the south-west emission is behind the plane of the galaxy and is extinguished very likely by dust in the plane. Intermediate-age stellar population contributes significantly to the continuum, especially in the inner 200 pc.

An analysis of the composite stellar population in M32

We obtained long-slit spectra of high signal-to-noise ratio of the galaxy M32 with the Gemini Multi-Object Spectrograph at the Gemini-North telescope. We analysed the integrated spectra by means of full spectral fitting in order to extract the mixture of stellar populations that best represents its composite nature. Three different galactic radii were analysed, from the nuclear region out to 2 arcmin from the centre. This allows us to compare, for the first time, the results of integrated light spectroscopy with those of resolved colour-magnitude diagrams from the literature. As a main result we propose that an ancient and an intermediate-age population co-exist in M32, and that the balance between these two populations change between the nucleus and outside one effective radius (1r(eff)) in the sense that the contribution from the intermediate population is larger at the nuclear region. We retrieve a smaller signal of a young population at all radii whose origin is unclear and may be a contamination from horizontal branch stars, such as the ones identified by Brown et al. in the nuclear region. We compare our metallicity distribution function for a region 1 to 2 arcmin from the centre to the one obtained with photometric data by Grillmair et al. Both distributions are broad, but our spectroscopically derived distribution has a significant component with [Z/Z(circle dot)] <= -1, which is not found by Grillmair et al.

The 6dF Galaxy Survey: stellar population trends across and through the Fundamental Plane

We present results from an analysis of stellar population parameters for 7132 galaxies in the 6dF Galaxy Survey Fundamental Plane (FP) sample. We bin the galaxies along the axes, v1, v2 and v3, of the tri-variate Gaussian to which we have fitted the galaxy distribution in effective radius, surface brightness and central velocity dispersion (FP space), and compute median values of stellar age, [Fe/H], [Z/H] and [a/Fe]. We determine the directions of the vectors in FP space along which each of the binned stellar population parameters vary most strongly. In contrast to previous work, we find stellar population trends not just with velocity dispersion and FP residual, but with radius and surface brightness as well. The most remarkable finding is that the stellar population parameters vary through the plane (v1 direction) and across the plane (v3 direction), but show no variation at all along the plane (v2 direction). The v2 direction in FP space roughly corresponds to luminosity density. We interpret a galaxys position along this vector as being closely tied to its merger history, such that early-type galaxies with lower luminosity density are more likely to have undergone major mergers. This conclusion is reinforced by an examination of the simulations of Kobayashi, which show clear trends of merger history with v2.

Irradiated disks and planet population synthesis

Recent planet population synthesis models (Alibert et al. 2010, submitted) have emphasized the key role played by the proto-planetary disk properties in determining the overall planet population characteristics. We present a disk model that takes into account viscous heating and irradiation by a central star. We consider the case of an equilibrium flaring angle. We illustrate the consequences of the resulting changes in the disk structure on the planet population by the synthetic populations corresponding to each of the different structures.

Planet Population Synthesis

With the increasing number of exoplanets discovered, statistical properties of the population as a whole become unique constraints on planet-formation models, provided a link between the description of the detailed processes playing a role in this formation and the observed population can be established. Planet population synthesis provides such a link. This approach allows us to study how different physical models of individual processes (e.g., protoplanetary disk structure and evolution, planetesimal formation, gas accretion, migration, etc.) affect the overall properties of the population of emerging planets. By necessity, planet population synthesis relies on simplified descriptions of complex processes. These descriptions can be obtained from more detailed specialized simulations of these processes. The objective of this chapter is twofold: (1) provide an overview of the physics entering in the two main approaches to planet population synthesis, and (2) present some of the results achieved as well as illustrate how it can be used to extract constraints on the models and to help interpret observations.

Population synthesis of isolated neutron stars with magneto-rotational evolution

We revisit the population synthesis of isolated radio-pulsars incorporating recent advances on the evolution of the magnetic field and the angle between the magnetic and rotational axes from new simulations of the magneto-thermal evolution and magnetosphere models, respectively. An interesting novelty in our approach is that we do not assume the existence of a death line. We discuss regions in parameter space that are more consistent with the observational data. In particular, we find that any broad distribution of birth spin periods with P0 ≲ 0.5 s can fit the data, and that if the alignment angle is allowed to vary consistently with the torque model, realistic magnetospheric models are favoured compared to models with classical magneto-dipolar radiation losses. Assuming that the initial magnetic field is given by a lognormal distribution, our optimal model has mean strength 〈log B0[G]〉 ≈ 13.0–13.2 with width σ(log B0) = 0.6–0.7. However, there are strong correlations between parameters. This degeneracy in the parameter space can be broken by an independent estimate of the pulsar birth rate or by future studies correlating this information with the population in other observational bands (X-rays and γ-rays).

Population synthesis of isolated neutron stars with magneto-rotational evolution – II. From radio-pulsars to magnetars

Population synthesis studies constitute a powerful method to reconstruct the birth distribution of periods and magnetic fields of the pulsar population. When this method is applied to populations in different wavelengths, it can break the degeneracy in the inferred properties of initial distributions that arises from single-band studies. In this context, we extend previous works to include X-ray thermal emitting pulsars within the same evolutionary model as radio-pulsars. We find that the cumulative distribution of the number of X-ray pulsars can be well reproduced by several models that, simultaneously, reproduce the characteristics of the radio-pulsar distribution. However, even considering the most favourable magneto-thermal evolution models with fast field decay, lognormal distributions of the initial magnetic field overpredict the number of visible sources with periods longer than 12 s. We then show that the problem can be solved with different distributions of magnetic field, such as a truncated lognormal distribution, or a binormal distribution with two distinct populations. We use the observational lack of isolated neutron stars (NSs) with spin periods P > 12 s to establish an upper limit to the fraction of magnetars born with B > 1015 G (less than 1 per cent). As future detections keep increasing the magnetar and high-B pulsar statistics, our approach can be used to establish a severe constraint on the maximum magnetic field at birth of NSs.

The VLT-FLAMES Tarantula Survey XXI. Stellar spin rates of O-type spectroscopic binaries

Context: The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are, however, found in multiple systems. 

Aims: By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary subpopulations with one another and with that of presumed-single stars in the same region, we aim to constrain the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the spin rates of massive stars. 

Methods: We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (ν_esini) for components of 114 spectroscopic binaries in 30 Doradus. The ν_esini values are derived from the full width at half maximum (FWHM) of a set of spectral lines, using a FWHM vs. ν_esini calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. 

Results: The overall ν_esini distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at ν_esini<200 kms^-1) and a shoulder at intermediate velocities (200 <ν_esini<300 kms^-1). The distributions of binaries and single stars, however, differ in two ways. First, the main peak at ν_esini ~ 100kms^-1 is broader and slightly shifted towards higher spin rates in the binary distribution than that of the presumed-single stars. This shift is mostly due to short-period binaries (P_orb~<10 d). Second, the ν_esini distribution of primaries lacks a significant population of stars spinning faster than 300 kms^-1, while such a population is clearly present in the single-star sample. The ν_esini distribution of binaries with amplitudes of radial velocity variation in the range of 20 to 200 kms^-1 (mostly binaries with P_orb ~ 10-1000 d and/or with q<0.5) is similar to that of single O stars below ν_esini_~<170kms^-1. 

Conclusions: Our results are compatible with the assumption that binary components formed with the same spin distribution as single stars, and that this distribution contains few or no fast-spinning stars. The higher average spin rate of stars in short-period binaries may either be explained by spin-up through tides in such tight binary systems, or by spin-down of a fraction of the presumed-single stars and long-period binaries through magnetic braking (or by a combination of both mechanisms). Most primaries and secondaries of SB2 systems with P_orb~<10 d appear to have similar rotational velocities. This is in agreement with tidal locking in close binaries where the components have similar radii. The lack of very rapidly spinning stars among binary systems supports the idea that most stars with ν_esini_~> 300kms^-1 in the single-star sample are actually spun-up post-binary interaction products. Finally, the overall similarities (low-velocity peak and intermediate-velocity shoulder) of the spin distribution of binary and single stars argue for a massive star formation process in which the initial spin is set independently of whether stars are formed as single stars or as components of a binary system.

GALAXY EVOLUTION THROUGH THE COSMIC TIME

Elliptical galaxies are the best systems to study the early star formation activity in the universe. This work aims to understand the formation and evolution of these objects through the study of the integrated properties of their stellar populations. Here an evolutionary model is developed and their predicted spectrophotometric properties are presented.