We analyse the correlations between continuum properties and emission line equivalent widths of star-forming and active galaxies from the Sloan Digital Sky Survey. Since upcoming large sky surveys will make broad-band observations only, including strong emission lines into theoretical modelling of spectra will be essential to estimate physical properties of photometric galaxies. We show that emission line equivalent widths can be fairly well reconstructed from the stellar continuum using local multiple linear regression in the continuum principal component analysis (PCA) space. Line reconstruction is good for star-forming galaxies and reasonable for galaxies with active nuclei. We propose a practical method to combine stellar population synthesis models with empirical modelling of emission lines. The technique will help generate more accurate model spectra and mock catalogues of galaxies to fit observations of the new surveys. More accurate modelling of emission lines is also expected to improve template-based photometric redshift estimation methods. We also show that, by combining PCA coefficients from the pure continuum and the emission lines, automatic distinction between hosts of weak active galactic nuclei (AGNs) and quiescent star-forming galaxies can be made. The classification method is based on a training set consisting of high-confidence starburst galaxies and AGNs, and allows for the similar separation of active and star-forming galaxies as the empirical curve found by Kauffmann et al. We demonstrate the use of three important machine learning algorithms in the paper: k-nearest neighbour finding, k-means clustering and support vector machines.

Understanding the diversity in spectra is the key to determining the physical parameters of galaxies. The optical spectra of galaxies are highly convoluted with continuum and lines which are potentially sensitive to different physical parameters. Defining the wavelength regions of interest is therefore an important question. In this work, we identify informative wavelength regions in a single-burst stellar populations model by using the CUR Matrix Decomposition. Simulating the Lick/IDS spectrograph configuration, we recover the widely used Dn(4000), Hbeta, and HdeltaA to be most informative. Simulating the SDSS spectrograph configuration with a wavelength range 3450-8350 Angstrom and a model-limited spectral resolution of 3 Angstrom, the most informative regions are: first region-the 4000 Angstrom break and the Hdelta line; second region-the Fe-like indices; third region-the Hbeta line; fourth region-the G band and the Hgamma line. A Principal Component Analysis on the first region shows that the first eigenspectrum tells primarily the stellar age, the second eigenspectrum is related to the age-metallicity degeneracy, and the third eigenspectrum shows an anti-correlation between the strengths of the Balmer and the Ca K and H absorptions. The regions can be used to determine the stellar age and metallicity in early-type galaxies which have solar abundance ratios, no dust, and a single-burst star formation history. The region identification method can be applied to any set of spectra of the user's interest, so that we eliminate the need for a common, fixed-resolution index system. We discuss future directions in extending the current analysis to late-type galaxies.

We present a new catalog of spectroscopically-confirmed white dwarf stars from the Sloan Digital Sky Survey Data Release 7 spectroscopic catalog. We find 20,407 white dwarf spectra, representing 19,712 stars, and provide atmospheric model fits to 14,120 DA and 1011 DB white dwarf spectra from 12,843 and 923 stars, respectively. These numbers represent a more than factor of two increase in the total number of white dwarf stars from the previous SDSS white dwarf catalog based on DR4 data. Our distribution of subtypes varies from previous catalogs due to our more conservative, manual classifications of each star in our catalog, supplementing our automatic fits. In particular, we find a large number of magnetic white dwarf stars whose small Zeeman splittings mimic increased Stark broadening that would otherwise result in an overestimated log(g) if fit as a non-magnetic white dwarf. We calculate mean DA and DB masses for our clean, non-magnetic sample and find the DB mean mass is statistically larger than that for the DAs.

In this work we present an atlas of composite spectra of galaxies based on the data of the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Galaxies are classified by colour, nuclear activity and star-formation activity to calculate average spectra of high signal-to-noise ratio and resolution (S/N = 132 - 4760 at {Dlambda = 1 A), using an algorithm that is robust against outliers. Besides composite spectra, we also compute the first five principal components of the distributions in each galaxy class to characterize the nature of variations of individual spectra around the averages. The continua of the composite spectra are fitted with BC03 stellar population synthesis models to extend the wavelength coverage beyond the coverage of the SDSS spectrographs. Common derived parameters of the composites are also calculated: integrated colours in the most popular filter systems, line strength measurements, and continuum absorption indices (including Lick indices). These derived parameters are compared with the distributions of parameters of individual galaxies and it is shown on many examples that the composites of the atlas cover much of the parameter space spanned by SDSS galaxies. By co-adding thousands of spectra, a total integration time of several months can be reached, which results in extremely low noise composites. The variations in redshift not only allow for extending the spectral coverage bluewards to the original wavelength limit of the SDSS spectrographs, but also make higher spectral resolution achievable. The composite spectrum atlas is available online at http://www.vo.elte.hu/compositeatlas.

Star formation is one of the most important processes in galaxy formation. The luminosity of Halpha recombination and [OII] forbidden emissions remain to be most used in measuring formation rate of massive stars in galaxies. Here we report the inclination dependency of continuum-subtracted and aperture-corrected nebular luminosities, including Halpha, Hbeta, Hgamma, [OII], [NII], of disk-dominated galaxies in the local universe. Their luminosities decrease by a factor of three from face-on to edge-on (axis ratio limit = 0.17) orientations. This dependence is deduced to be caused by extinction due to diffuse dust within the disks with an amplitude of 1.2 mag. The line-luminosity--inclination relation provides a novel way to remove extinction in emission lines and present star formation rate of disk galaxies out to redshift of 1.6.

The inclination of galaxies induces both reddening and extinction to their observed spectral energy distribution, which in turn impact the derived properties of the galaxies. Here we report a significant dependence of the error in photometric redshift (photo-z) on the inclination of disk galaxies from the Sloan Digital Sky Survey. The bias in the photo-z based on the template-fitting approach increases from -0.015 in face-on to 0.021 in edge-on galaxies. A Principal Component Analysis on the full sample of photometry reveals the inclination of the galaxies to be represented by the 2nd mode. The corresponding eigenspectrum resembles an extinction curve. The isolation of the inclination effect in a low-order mode demonstrates the significant reddening induced on the observed colors, leading to the over-estimated photo-z in galaxies of high inclinations. We present approaches to correct the photo-z and the other properties of the disk galaxies against the inclination effect.

The Void Galaxy Survey (VGS) is a multi-wavelength program to study $\sim$60 void galaxies. Each has been selected from the deepest interior regions of identified voids in the SDSS redshift survey on the basis of a unique geometric technique, with no a prior selection of intrinsic properties of the void galaxies. The project intends to study in detail the gas content, star formation history and stellar content, as well as kinematics and dynamics of void galaxies and their companions in a broad sample of void environments. It involves the HI imaging of the gas distribution in each of the VGS galaxies. Amongst its most tantalizing findings is the possible evidence for cold gas accretion in some of the most interesting objects, amongst which are a polar ring galaxy and a filamentary configuration of void galaxies. Here we shortly describe the scope of the VGS and the results of the full analysis of the pilot sample of 15 void galaxies.

Void galaxies, residing within the deepest underdensities of the Cosmic Web, present an ideal population for the study of galaxy formation and evolution in an environment undisturbed by the complex processes modifying galaxies in clusters and groups, as well as provide an observational test for theories of cosmological structure formation. We have completed a pilot survey for the HI imaging aspects of a new Void Galaxy Survey (VGS), imaging 15 void galaxies in HI in local (d < 100 Mpc) voids. HI masses range from 3.5 x 10^8 to 3.8 x 10^9 M_sun, with one nondetection with an upper limit of 2.1 x 10^8 M_sun. Our galaxies were selected using a structural and geometric technique to produce a sample that is purely environmentally selected and uniformly represents the void galaxy population. In addition, we use a powerful new backend of the Westerbork Synthesis Radio Telescope that allows us to probe a large volume around each targeted galaxy, simultaneously providing an environmentally constrained sample of fore- and background control sample of galaxies while still resolving individual galaxy kinematics and detecting faint companions in HI. This small sample makes up a surprisingly interesting collection of perturbed and interacting galaxies, all with small stellar disks. Four galaxies have significantly perturbed HI disks, five have previously unidentified companions at distances ranging from 50 to 200 kpc, two are in interacting systems, and one was found to have a polar HI disk. Our initial findings suggest void galaxies are a gas-rich, dynamic population which present evidence of ongoing gas accretion, major and minor interactions, and filamentary alignment despite the surrounding underdense environment.

Extinction in galaxies affects their observed properties. In scenarios describing the distribution of dust and stars in individual disk galaxies, the amplitude of the extinction can be modulated by the inclination of the galaxies. In this work we investigate the inclination dependency in composite spectra of star-forming disk galaxies from the Sloan Digital Sky Survey Data Release 5. In a volume-limited sample within a redshift range 0.065-0.075 and a r-band Petrosian absolute magnitude range -19.5 to -22 which exhibits a flat distribution of inclination, the inclined relative to face-on extinction in the stellar continuum is found empirically to increase with inclination in the g, r, and i bands. Within the central 0.5 intrinsic half-light radius of the galaxies, the g-band relative extinction in the stellar continuum for the highly-inclined objects (axis ratio b/a = 0.1) is 1.2 mag, agreeing with previous studies. The extinction curve of the disk galaxies is given in the restframe wavelengths 3700-8000 angstrom, identified with major optical emission and absorption lines in diagnostics. The Balmer decrement remains constant with inclination, suggesting a different kind of dust configuration and/or reddening mechanism in the HII region from that in the stellar continuum. One factor is shown to be the presence of spatially non-uniform interstellar extinction, presumably caused by clumped dust in the vicinity of the HII region.

Many analyses and parameter estimations undertaken in astronomy require a large set (> 10^5) of non-analytical, theoretical spectra, each of these defined by multiple parameters. We describe the construction of an N-dimensional grid which is suitable for generating such spectra. The theoretical spectra are designed to correspond to a targeted parameter grid but otherwise to random positions in the parameter space, and they are interpolated on-the-fly through a pre-calculated grid of spectra. The initial grid is designed to be relatively low in parameter resolution and small in occupied hard disk space and therefore can be updated efficiently when a new model is desired. In a pilot study of stellar population synthesis of galaxies, the mean square errors on the estimated parameters are found to decrease with the targeted grid resolution. This scheme of generating a large model grid is general for other areas of studies, particularly if they are based on multi-dimensional parameter space and are focused on contrasting model differences.

We have completed a pilot survey imaging 15 SDSS selected void galaxies in HI in local (d=50 to 100 Mpc) voids. This small sample makes up a surprisingly interesting collection of galaxies, consisting of galaxies with asymmetric and perturbed HI disks, previously unidentified companions, and ongoing interactions. One was found to have a polar HI disk with no stellar counterpart. While our small number statistics so far are limiting, results support past findings that most void galaxies are typically late type galaxies with gas rich disks and small scale clustering similar to field galaxies despite their large scale underdense environment.

We present a novel technique to overcome the limitations of the applicability of Principal Component Analysis to typical real-life data sets, especially astronomical spectra. Our new approach addresses the issues of outliers, missing information, large number of dimensions and the vast amount of data by combining elements of robust statistics and recursive algorithms that provide improved eigensystem estimates step-by-step. We develop a generic mechanism for deriving reliable eigenspectra without manual data censoring, while utilising all the information contained in the observations. We demonstrate the power of the methodology on the attractive collection of the VIMOS VLT Deep Survey spectra that manifest most of the challenges today, and highlight the improvements over previous workarounds, as well as the scalability of our approach to collections with sizes of the Sloan Digital Sky Survey and beyond.

Under the unified model for active galactic nuclei (AGNs), narrow-line (Type 2) AGNs are, in fact, broad-line (Type 1) AGNs but each with a heavily obscured accretion disk. We would therefore expect the optical continuum emission from Type 2 AGN to be composed mainly of stellar light and non-variable on the time-scales of months to years. In this work we probe the spectroscopic variability of galaxies and narrow-line AGNs using the multi-epoch data in the Sloan Digital Sky Survey (SDSS) Data Release 6. The sample contains 18,435 sources for which there exist pairs of spectroscopic observations (with a maximum separation in time of ~700 days) covering a wavelength range of 3900-8900 angstrom. To obtain a reliable repeatability measurement between each spectral pair, we consider a number of techniques for spectrophotometric calibration resulting in an improved spectrophotometric calibration of a factor of two. From these data we find no obvious continuum and emission-line variability in the narrow-line AGNs on average -- the spectroscopic variability of the continuum is 0.07+/-0.26 mag in the g band and, for the emission-line ratios log10([NII]/Halpha) and log10([OIII]/Hbeta), the variability is 0.02+/-0.03 dex and 0.06+/-0.08 dex, respectively. From the continuum variability measurement we set an upper limit on the ratio between the flux of varying spectral component, presumably related to AGN activities, and that of host galaxy to be ~30%. We provide the corresponding upper limits for other spectral classes, including those from the BPT diagram, eClass galaxy classification, stars and quasars.

I summarize the current spectroscopic sky surveys and some of the scientific results, emphasizing the largest sky survey to-date, the Sloan Digital Sky Survey (SDSS). Techniques used commonly in spectral analyses are discussed, followed by the present needs and challenges for solving some of the unknown problems. I discuss how the Virtual Observatory (VO) can help astronomers in carrying out related research.

The evolution of stellar disks is of great importance for understanding many aspects of galaxy formation. In this work we perform stellar population synthesis on radially resolved photometry of 564 disk galaxies from the SDSS DR5, selected to have both spectra of the central regions and photometry. To explore fully the multi-dimensional likelihood space defined by the output parameters of the spectral synthesis, we use Markov Chain Monte Carlo to quantify the expectation values, the uncertainties and the degeneracies of the parameters. We find good agreement between the parameter values obtained using the SDSS broad-band colors and the spectra respectively. In general the derived mean stellar age and the best-fit stellar metallicity decline in value from the galaxy center to the outer regions (around 1.5 half-light radii), based on sub-samples defined by concentration index. We also find that the radial dependency of the stellar population parameters exhibits a significant variation, and this diversity is likely related to morphology and the physics of star formation.

Using subspace methods, we study the distribution of physical components of galaxies in wavelength space. We find that it is valid to assume that the stellar and the gaseous components of galaxies span complementary subspaces. To first order, stellar and gaseous spectral features can be extracted from galaxy spectra through a simple matrix multiplication. By comparing the stellar continua obtained respectively using the model-based and the empirical approach through a commonality measure, we conclude that the latter may lose higher-order spectral features.

We have studied a sample of 809 Mg II absorption systems with 1.0 < z_abs < 1.86 in the spectra of SDSS QSOs, with the aim of understanding the nature and abundance of the dust and the chemical abundances in the intervening absorbers. Normalized, composite spectra were derived, for abundance measurements, for the full sample and several sub-samples, chosen on the basis of the line strengths and other absorber and QSO properties. Average extinction curves were obtained for the sub-samples by comparing their geometric mean spectra with those of matching samples of QSOs without absorbers in their spectra. There is clear evidence for the presence of dust in the intervening absorbers. The 2175 A feature is not present in the extinction curves, for any of the sub-samples. The extinction curves are similar to the SMC extinction curve with a rising UV extinction below 2200 A. The absorber rest frame colour excess, E(B-V), derived from the extinction curves, depends on the absorber properties and ranges from < 0.001 to 0.085 for various sub-samples. The column densities of several ions do not show such a correspondingly large variation. The depletion pattern is similar to halo clouds in the Galaxy. Assuming an SMC gas-to-dust ratio we find a trend of increasing abundance with decreasing extinction; systems with N_H I \~ 10^{20} cm^{-2} show solar abundance of Zn. The large velocity spread of strong Mg II systems seems to be mimicked by weak lines of other elements. The ionization of the absorbers, in general appears to be low. QSOs with absorbers are, in general, at least three times as likely to have highly reddened spectra as compared to QSOs without any absorption systems in their spectra.

Using an eigenspectrum decomposition technique, we separate the host galaxy from the broad line active galactic nucleus (AGN) in a set of 4666 spectra from the Sloan Digital Sky Survey (SDSS), from redshifts near zero up to about 0.75. The decomposition technique uses separate sets of galaxy and quasar eigenspectra to efficiently and reliably separate the AGN and host spectroscopic components. The technique accurately reproduces the host galaxy spectrum, its contributing fraction, and its classification. We show how the accuracy of the decomposition depends upon S/N, host galaxy fraction, and the galaxy class. Based on the eigencoefficients, the sample of SDSS broad-line AGN host galaxies spans a wide range of spectral types, but the distribution differs significantly from inactive galaxies. In particular, post-starburst activity appears to be much more common among AGN host galaxies. The luminosities of the hosts are much higher than expected for normal early-type galaxies, and their colors become increasingly bluer than early-type galaxies with increasing host luminosity. Most of the AGNs with detected hosts are emitting at between 1% and 10% of their estimated Eddington luminosities, but the sensitivity of the technique usually does not extend to the Eddington limit. There are mild correlations among the AGN and host galaxy eigencoefficients, possibly indicating a link between recent star formation and the onset of AGN activity. The catalog of spectral reconstruction parameters is available as an electronic table.

We find evidence for dust in the intervening QSO absorbers from the spectra of QSOs in the Sloan Digital Sky Survey Data Release 1. No evidence is found for the 2175 A feature which is present in the Milky Way dust extinction curve. The extinction curve resembles the SMC extinction curve. The observed Delta(g-i) excess for QSOs with strong absorption systems appears to be a result of the reddening due to dust in the intervening absorbers.

We study 16,707 quasar spectra from the SDSS using the Karhunen-Lo\`eve (KL) transform (or Principal Component Analysis, PCA). The quasar eigenspectra of the full catalog reveal the following: 1st order - the mean spectrum; 2nd order - a host-galaxy component; 3rd order - the UV-optical continuum slope; 4th order - the correlations of Balmer emission lines. We find that the spectral classification of quasars is redshift and luminosity dependent, as such there does not exist a compact set (< 10 modes) of eigenspectra (covering 900 to 8000 angstrom) which can describe most variations (> 95%) of the entire catalog. We therefore construct several sets of eigenspectra in different redshift and luminosity bins. From these eigenspectra we find that quasar spectra can be classified (by the first two eigenspectra) into a sequence that is defined by a simple progression in the steepness of the slope of the continuum. We also find a dependence on redshift and luminosity in the eigencoefficients. The dominant redshift effect is a result of the evolution of the blended FeII emission (optical) and the Balmer continuum (the ``small bump''). A luminosity dependence is also present in the eigencoefficients and is related to the Baldwin effect. The eigenspectra constructed from the wavelength-selected SDSS spectra are found to agree with the principal components by Francis et al. (1992) and the well-known ``Eigenvector-1'' (Boroson and Green 1992). (Abridged)

We perform an objective classification of 170,000 galaxy spectra in the Sloan Digital Sky Survey (SDSS) using the Karhunen-Lo\`eve (KL) transform. With about one-sixth of the total set of galaxy spectra which will be obtained by the survey, we are able to carry out the most extensive analysis of its kind to date. The formalism proposed by Connolly and Szalay is adopted to correct for gappy regions in the spectra, and to derive eigenspectra and eigencoefficients. From this analysis, we show that this gap-correction formalism leads to a converging set of eigenspectra and KL-repaired spectra. Furthermore, KL eigenspectra of galaxies are found to be convergent not only as a function of iteration, but also as a function of the number of randomly selected galaxy spectra used in the analysis. From these data a set of ten eigenspectra of galaxy spectra are constructed, with rest-wavelength coverage 3450-8350 angstrom. The eigencoefficients describing these galaxies naturally place the spectra into several classes defined by the plane formed by the first three eigencoefficients of each spectrum. Spectral types, corresponding to different Hubble-types and galaxies with extreme emission lines, are identified for the 170,000 spectra and are shown to be complementary to existing spectral classifications. From a non-parametric classification technique, we find that the population of galaxies can be divided into three classes which correspond to early late- through to intermediate late-types galaxies. (Abridged)

We construct eigenspectra from the DR1 quasars in the SDSS using the Karhunen-Lo\`eve (KL) transform (or Principal Component Analysis, PCA) in different redshift and luminosity bins. We find that the quasar spectra can be classified, by the first two eigenspectra, into a continuous sequence in the variation of the spectral slope. We also find a dependence on redshift and luminosity in the eigencoefficients. The dominant redshift effect is the evolution of the blended Fe~II emission (optical) and the Balmer continuum (the ``small-bump'', $\lambda_{rest} \approx 2000-4000$\AA), while the luminosity effect is related to the Baldwin effect. Correlations among several major broad emission lines are found, including the well-known ``Eigenvector-1''.

Using a large sample of quasar spectra from the SDSS, we examine the composite spectral trends of quasars as functions of both redshift and luminosity, independently of one another. Aside from the well known Baldwin effect (BE) -- the decrease of line equivalent width with luminosity -- the average spectral properties are remarkably similar. Host galaxy contamination and the BE are the primary causes for apparent changes in the average spectral slope of the quasars. The BE is detected for most emission lines, including the Balmer lines, but with several exceptions including NV1240A. Emission line shifts of several lines are associated with the BE. The BE is mainly a function of luminosity, but also partly a function of redshift in that line equivalent widths become stronger with redshift. Some of the complex iron features change with redshift, particularly near the small blue bump region.