91̽»¨

Abstract:

To ascertain whether photometric decompositions of galaxies into bulges and discs are astrophysically meaningful, we have developed a new technique to decompose spectral data cubes into separate bulge and disc components, subject only to the constraint that they reproduce the conventional photometric decomposition. These decompositions allow us to study the kinematic and stellar population properties of the individual components and how they vary with position, in order to assess their plausibility as discrete elements, and to start to reconstruct their distinct formation histories. An initial application of this method to Calar Alto Integral Field Area integral field unit observations of three isolated S0 galaxies confirms that in regions where both bulge and disc contribute significantly to the flux, they can be physically and robustly decomposed into a rotating dispersion-dominated bulge component and a rotating low-dispersion disc component. Analysis of the resulting stellar populations shows that the bulges of these galaxies have a range of ages relative to their discs, indicating that a variety of processes are necessary to describe their evolution. This simple test case indicates the broad potential for extracting from spectral data cubes the full spectral data of a wide variety of individual galaxy components, and for using such decompositions to understand the interplay between these various structures, and hence how such systems formed.

Abstract:

I start by providing an updated summary of the penalized pixel-fitting (ppxf) method, which is used to extract the stellar and gas kinematics, as well as the stellar population of galaxies, via full spectrum fitting. I then focus on the problem of extracting the kinematic when the velocity dispersion σ is smaller than the velocity sampling ΔV, which is generally, by design, close to the instrumental dispersion σinst. The standard approach consists of convolving templates with a discretized kernel, while fitting for its parameters. This is obviously very inaccurate when σ ≲ ΔV=2, due to undersampling. Oversampling can prevent this, but it has drawbacks. Here I present a more accurate and efficient alternative. It avoids the evaluation of the under-sampled kernel, and instead directly computes its well-sampled analytic Fourier transform, for use with the convolution theorem. A simple analytic transform exists when the kernel is described by the popular Gauss-Hermite parametrization (which includes the Gaussian as special case) for the line-of-sight velocity distribution. I describe how this idea was implemented in a significant upgrade to the publicly available ppxf software. The key advantage of the new approach is that it provides accurate velocities regardless of σ. This is important e.g. for spectroscopic surveys targeting galaxies with σ << σinst, for galaxy redshift determinations, or for measuring line-of-sight velocities of individual stars. The proposed method could also be used to fix Gaussian convolution algorithms used in today’s popular software packages.