Molecular line and continuum studies of the early stages of star formation

New stars form in dense interstellar clouds of gas and dust called molecular clouds. The actual sites where the process of star formation takes place are the dense clumps and cores deeply embedded in molecular clouds. The details of the star formation process are complex and not completely understood. Thus, determining the physical and chemical properties of molecular cloud cores is necessary for a better understanding of how stars are formed. Some of the main features of the origin of low-mass stars, like the Sun, are already relatively well-known, though many details of the process are still under debate. The mechanism through which high-mass stars form, on the other hand, is poorly understood. Although it is likely that the formation of high-mass stars shares many properties similar to those of low-mass stars, the very first steps of the evolutionary sequence are unclear. Observational studies of star formation are carried out particularly at infrared, submillimetre, millimetre, and radio wavelengths. Much of our knowledge about the early stages of star formation in our Milky Way galaxy is obtained through molecular spectral line and dust continuum observations. The continuum emission of cold dust is one of the best tracers of the column density of molecular hydrogen, the main constituent of molecular clouds. Consequently, dust continuum observations provide a powerful tool to map large portions across molecular clouds, and to identify the dense star-forming sites within them. Molecular line observations, on the other hand, provide information on the gas kinematics and temperature. Together, these two observational tools provide an efficient way to study the dense interstellar gas and the associated dust that form new stars. The properties of highly obscured young stars can be further examined through radio continuum observations at centimetre wavelengths. For example, radio continuum emission carries useful information on conditions in the protostar+disk interaction region where protostellar jets are launched. In this PhD thesis, we study the physical and chemical properties of dense clumps and cores in both low- and high-mass star-forming regions. The sources are mainly studied in a statistical sense, but also in more detail. In this way, we are able to examine the general characteristics of the early stages of star formation, cloud properties on large scales (such as fragmentation), and some of the initial conditions of the collapse process that leads to the formation of a star. The studies presented in this thesis are mainly based on molecular line and dust continuum observations. These are combined with archival observations at infrared wavelengths in order to study the protostellar content of the cloud cores. In addition, centimetre radio continuum emission from young stellar objects (YSOs; i.e., protostars and pre-main sequence stars) is studied in this thesis to determine their evolutionary stages. The main results of this thesis are as follows: i) filamentary and sheet-like molecular cloud structures, such as infrared dark clouds (IRDCs), are likely to be caused by supersonic turbulence but their fragmentation at the scale of cores could be due to gravo-thermal instability; ii) the core evolution in the Orion B9 star-forming region appears to be dynamic and the role played by slow ambipolar diffusion in the formation and collapse of the cores may not be significant; iii) the study of the R CrA star-forming region suggests that the centimetre radio emission properties of a YSO are likely to change with its evolutionary stage; iv) the IRDC G304.74+01.32 contains candidate high-mass starless cores which may represent the very first steps of high-mass star and star cluster formation; v) SiO outflow signatures are seen in several high-mass star-forming regions which suggest that high-mass stars form in a similar way as their low-mass counterparts, i.e., via disk accretion. The results presented in this thesis provide constraints on the initial conditions and early stages of both low- and high-mass star formation. In particular, this thesis presents several observational results on the early stages of clustered star formation, which is the dominant mode of star formation in our Galaxy.

A Triggering Mechanism for Enhanced Star-Formation in Colliding Galaxies

We propose a physical mechanism to explain the origin of the intense burst of massive-star formation seen in colliding/merging, gas-rich, field spiral galaxies. We explicitly take account of the different parameters for the two main mass components, H-2 and H I, of the interstellar medium within a galaxy and follow their consequent different evolution during a collision between two galaxies. We also note that, in a typical spiral galaxy-like our galaxy, the Giant Molecular Clouds (GMCs) are in a near-virial equilibrium and form the current sites of massive-star formation, but have a low star formation rate. We show that this star formation rate is increased following a collision between galaxies. During a typical collision between two field spiral galaxies, the H I clouds from the two galaxies undergo collisions at a relative velocity of approximately 300 km s-1. However, the GMCs, with their smaller volume filling factor, do not collide. The collisions among the H I clouds from the two galaxies lead to the formation of a hot, ionized, high-pressure remnant gas. The over-pressure due to this hot gas causes a radiative shock compression of the outer layers of a preexisting GMC in the overlapping wedge region. This makes these layers gravitationally unstable, thus triggering a burst of massive-star formation in the initially barely stable GMCs.The resulting value of the typical IR luminosity from the young, massive stars from a pair of colliding galaxies is estimated to be approximately 2 x 10(11) L., in agreement with the observed values. In our model, the massive-star formation occurs in situ in the overlapping regions of a pair of colliding galaxies. We can thus explain the origin of enhanced star formation over an extended, central area approximately several kiloparsecs in size, as seen in typical colliding galaxies, and also the origin of starbursts in extranuclear regions of disk overlap as seen in Arp 299 (NGC 3690/IC 694) and in Arp 244 (NGC 4038/39). Whether the IR emission from the central region or that from the surrounding extranuclear galactic disk dominates depends on the geometry and the epoch of the collision and on the initial radial gas distribution in the two galaxies. In general, the central starburst would be stronger than that in the disks, due to the higher preexisting gas densities in the central region. The burst of star formation is expected to last over a galactic gas disk crossing time approximately 4 x 10(7) yr. We can also explain the simultaneous existence of nearly normal CO galaxy luminosities and shocked H-2 gas, as seen in colliding field galaxies.This is a minimal model, in that the only necessary condition for it to work is that there should be a sufficient overlap between the spatial gas distributions of the colliding galaxy pair.

Star formation in giant extragalactic H II regions

Star formation properties in Giant Extragalactic H II Regions (GEHRs) are investigated using optical photometry and evolutionary population synthesis models. Photometric data in $BVR$ bands and in the emission line of H-alpha are obtained by CCD imaging at Vainu Bappu Observatory, Kavalur. Aperture photometry is performed for 180 GEHRs in galaxies NGC 1365, 1566, 2366, 2903, 2997, 3351, 4303, 4449, 4656 and 5253. Thirty six of these GEHRs having published spectroscopic data are studied for star formation properties. The population synthesis model is constructed based on Maeder's stellar evolutionary and Kurucz stellar atmosphere models, to synthesize observational quantities of embedded clusters in GEHRs. The observed H-alpha luminosity is a measure of the number of massive stars while the contribution to BVR bands is from intermediate mass (5-15 solar mass) stars when the cluster is young and from evolving supergiants when the cluster is old (age >/= 6~Myr). Differential reddening between gas and embedded stars is essential to constrain the dereddened cluster colors within the range of youngest clusters. Obscuring dust closely associated with gas, which is distributed in filaments and clumps, as in the case of 30 Doradus, is the most likely configuration giving rise to net reduction of extinction towards stars. The fraction of the stellar photons escaping the nebula unattenuated is estimated to be 50%. GEHRs are rarely found to be simple systems containing stars from single generation. In the present sample such regions in addition to being older than 3~Myr, have their Lyman continuum luminosity reduced by as much as 60%, compared to the observed $B$ band luminosity for a normal IMF. The missing ionizing photons may be escaping the nebula, leading to the ionization of extra-H II region ionized medium. Co-existence of young (age

A search for star formation around the Galactic halo B-type star PHL 346 with the 2dF spectrograph

A preliminary search for stars that may have formed coevally with the apparently young halo B-type star PHL 346 has been performed with the 2dF multifibre spectrograph on the Anglo- Australian Telescope (AAT). Candidates were selected for spectroscopy from APM scans of B and R Schmidt plates centred on PHL 346. A total of 476 stars of spectral type A or F were found; radial velocity estimates and more accurate spectral type assignments narrowed the number of possible coeval candidates to 6 A-type and 14 F-type stars. A statistical analysis of these results using a comparison with a control field suggests that the number of A-type or F-type candidate stars around PHL 346 is not unexpected, and that they need not be associated with PHL 346. A number of ways to improve the project are suggested.

SCUBA imaging of high mass star formation regions

We present 450 and 850 mu m images of five ultracompact HII regions (G10.47, G12.21, G13.87, G31.41 and G43.89) taken at 9

A comparison between star formation rate diagnostics and rate of core collapse supernovae within 11 Mpc

Aims. The core collapse supernova rate provides a strong lower limit for the star formation rate (SFR). Progress in using it as a cosmic SFR tracer requires some confidence that it is consistent with more conventional SFR diagnostics in the nearby Universe. This paper compares standard SFR measurements based on H alpha, far ultraviolet (FUV) and total infrared (TIR) galaxy luminosities with the observed core collapse supernova rate in the same galaxy sample. The comparison can be viewed from two perspectives. Firstly, by adopting an estimate of the minimum stellar mass to produce a core collapse supernova one can determine a SFR from supernova numbers. Secondly, the radiative SFR can be assumed to be robust and then the supernova statistics provide a constrain on the minimum stellar mass for core collapse supernova progenitors.

YSOVAR: Mid-infrared Variability Among YSOs in the Star Formation Region GGD12-15

We present an IR-monitoring survey with the Spitzer Space Telescope of the star-forming region GGD 12-15. More than 1000 objects were monitored, including about 350 objects within the central 5′, which is found to be especially dense in cluster members. The monitoring took place over 38 days and is part of the Young Stellar Object VARiability project. The region was also the subject of a contemporaneous 67 ks Chandra observation. The field includes 119 previously identified pre-main sequence star candidates. X-rays are detected from 164 objects, 90 of which are identified with cluster members. Overall, we find that about half the objects in the central 5′ are young stellar objects (YSOs) based on a combination of their spectral energy distribution, IR variability, and X-ray emission. Most of the stars with IR excess relative to a photosphere show large amplitude (>0.1 mag) mid-infrared (mid-IR) variability. There are 39 periodic sources, and all but one is found to be a cluster member. Almost half of the periodic sources do not show IR excesses. Overall, more than 85% of the Class I, flat spectrum, and Class II sources are found to vary. The amplitude of the variability is larger in more embedded YSOs. Most of the Class I/ II objects exhibit redder colors in a fainter state, which is compatible with time-variable extinction. A few become bluer when fainter, which can be explained with significant changes in the structure of the inner disk. A search for changes in the IR due to X-ray events is carried out, but the low number of flares prevented an analysis of the direct impact of X-ray flares on the IR light curves. However, we find that X-ray detected Class II sources have longer timescales for change in the MIR than a similar set of non-X-ray detected Class IIs.

Local star formation triggered by supernova shocks in magnetized diffuse neutral clouds

In this work, considering the impact of a supernova remnant (SNR) with a neutral magnetized cloud we derived analytically a set of conditions that are favourable for driving gravitational instability in the cloud and thus star formation. Using these conditions, we have built diagrams of the SNR radius, R(SNR), versus the initial cloud density, n(c), that constrain a domain in the parameter space where star formation is allowed. This work is an extension to previous study performed without considering magnetic fields (Melioli et al. 2006, hereafter Paper I). The diagrams are also tested with fully three-dimensional MHD radiative cooling simulations involving a SNR and a self-gravitating cloud and we find that the numerical analysis is consistent with the results predicted by the diagrams. While the inclusion of a homogeneous magnetic field approximately perpendicular to the impact velocity of the SNR with an intensity similar to 1 mu G within the cloud results only a small shrinking of the star formation zone in the diagram relative to that without magnetic field, a larger magnetic field (similar to 10 mu G) causes a significant shrinking, as expected. Though derived from simple analytical considerations these diagrams provide a useful tool for identifying sites where star formation could be triggered by the impact of a supernova blast wave. Applications of them to a few regions of our own Galaxy (e.g. the large CO shell in the direction of Cassiopeia, and the Edge Cloud 2 in the direction of the Scorpious constellation) have revealed that star formation in those sites could have been triggered by shock waves from SNRs for specific values of the initial neutral cloud density and the SNR radius. Finally, we have evaluated the effective star formation efficiency for this sort of interaction and found that it is generally smaller than the observed values in our own Galaxy (SFE similar to 0.01-0.3). This result is consistent with previous work in the literature and also suggests that the mechanism presently investigated, though very powerful to drive structure formation, supersonic turbulence and eventually, local star formation, does not seem to be sufficient to drive global star formation in normal star-forming galaxies, not even when the magnetic field in the neutral clouds is neglected.

Star formation in H I tails: HCG 92, HCG 100 and six interacting systems

We present new Gemini spectra of 14 new objects found within the H?i tails of Hickson Compact Groups (HCGs) 92 and 100. Nine of them are Galaxy Evolution Explorer (GALEX) far-ultraviolet (FUV) and near-ultraviolet (NUV) sources. The spectra confirm that these objects are members of the compact groups and have metallicities close to solar, with an average value of 12+log(O/H) similar to 8.5. They have average FUV luminosities 7 x 10(40)?erg?s-1 and very young ages (<100?Myr), and two of them resemble tidal dwarf galaxy (TDG) candidates. We suggest that they were created within gas clouds that were ejected during galaxygalaxy interactions into the intergalactic medium, which would explain the high metallicities of the objects, inherited from the parent galaxies from which the gas originated. We conduct a search for similar objects in six interacting systems with extended H?i tails: NGC 2623, NGC 3079, NGC 3359, NGC 3627, NGC 3718 and NGC 4656. We found 35 ultraviolet (UV) sources with ages < 100?Myr; however, most of them are on average less luminous/massive than the UV sources found around HCG 92 and HCG 100. We speculate that this might be an environmental effect and that compact groups of galaxies are more favourable to TDG formation than other interacting systems.

The star formation rate in the inner Milky Way Galaxy

The present star formation rate (SFR) in the inner Galaxy is puzzling for the chemical evolution models (CEM). No static CEM is able to reproduce the peak of the SFR in the 4 kpc ring. The main reason is probably a shortage of gas, which could be due to the dynamical effects produced by the galactic bar, not considered by these models. We developed a CEM that includes radial gas flows in order to mimic the effects of the galactic bar in the first 5 kpc of the galactic disk. In this model, the star formation (SF) is a two-step process: first, the diffuse gas forms molecular clouds. Then, stars form from cloud-cloud collisions or by the interaction between massive stars and the molecular gas. The former is called spontaneous and the latter induced SF. The mass in the different phases of each region changes by the processes associated with the stellar formation and death by: the SF due to spontaneous fragmentation of gas in the halo; formation of gas clouds in the disk from the diffuse gas; induced SF in the disk due to the interaction between massive stars and gas clouds; and finally, the restitution of the diffuse gas associated to these process of cloud and star formation. In the halo, the star formation rate for the diffuse gas follows a Schmidt law with a power n = 1.5. In the disk, the stars form in two steps: first, molecular clouds are formed from the diffuse gas also following a Schmidt law with n=1.5 and a proportionality factor. Including a specific pattern of radial gas flows, the CEM is able to reproduce with success the peak in the SFR at 4 kpc (fig. 1).

Molecular gas in the galaxy M83 : Its distribution, kinematics, and relation to star formation

The barred spiral galaxy M83 (NGC5236) has been observed in the 12CO J=1–0 and J=2–1 millimetre lines with the Swedish-ESO Submillimetre Telescope (SEST). The sizes of the CO maps are 100×100, and they cover the entire optical disk. The CO emission is strongly peaked toward the nucleus. The molecular spiral arms are clearly resolved and can be traced for about 360º. The total molecular gas mass is comparable to the total Hi mass, but H2 dominates in the optical disk. Iso-velocity maps show the signature of an inclined, rotating disk, but also the effects of streaming motions along the spiral arms. The dynamical mass is determined and compared to the gas mass. The pattern speed is determined from the residual velocity pattern, and the locations of various resonances are discussed. The molecular gas velocity dispersion is determined, and a trend of decreasing dispersion with increasing galactocentric radius is found. A total gas (H2+Hi+He) mass surface density map is presented, and compared to the critical density for star formation of an isothermal gaseous disk. The star formation rate (SFR) in the disk is estimated using data from various star formation tracers. The different SFR estimates agree well when corrections for extinctions, based on the total gas mass map, are made. The radial SFR distribution shows features that can be associated with kinematic resonances. We also find an increased star formation efficiency in the spiral arms. Different Schmidt laws are fitted to the data. The star formation properties of the nuclear region, based on high angular resolution HST data, are also discussed.

The Universe 10 Billion years ago: Morphologies, Star-formation rates and Galactic-scale winds in z ∼ 2 Galaxies

This thesis is devoted to the study of the properties of high-redsfhit galaxies in the epoch 1 < z < 3, when a substantial fraction of galaxy mass was assembled, and when the evolution of the star-formation rate density peaked. Following a multi-perspective approach and using the most recent and high-quality data available (spectra, photometry and imaging), the morphologies and the star-formation properties of high-redsfhit galaxies were investigated. Through an accurate morphological analyses, the built up of the Hubble sequence was placed around z ~ 2.5. High-redshift galaxies appear, in general, much more irregular and asymmetric than local ones. Moreover, the occurrence of morphological k-­correction is less pronounced than in the local Universe. Different star-formation rate indicators were also studied. The comparison of ultra-violet and optical based estimates, with the values derived from infra-red luminosity showed that the traditional way of addressing the dust obscuration is problematic, at high-redshifts, and new models of dust geometry and composition are required. Finally, by means of stacking techniques applied to rest-frame ultra-violet spectra of star-forming galaxies at z~2, the warm phase of galactic-scale outflows was studied. Evidence was found of escaping gas at velocities of ~ 100 km/s. Studying the correlation of inter-­stellar absorption lines equivalent widths with galaxy physical properties, the intensity of the outflow-related spectral features was proven to depend strongly on a combination of the velocity dispersion of the gas and its geometry.

Tracing AGN accretion and star formation in Herschel galaxies

In this Thesis, we study the physical properties and the cosmic evolution of AGN and their host galaxies since z∼3. Our analysis exploits samples of star forming galaxies detected with Herschel at far-IR wavelengths (from 70 up to 500 micron) in different extragalactic surveys, such as COSMOS and the deep GOODS (South and North) fields. The broad-band ancillary data available in COSMOS and the GOODS fields, allows us to implement Herschel and Spitzer photometry with multi-wavelength ancillary data. We perform a multicomponent SED-fitting decomposition to decouple the emission due to star formation from that due to AGN accretion, and to estimate both host-galaxy parameters (such as stellar mass, M* and star formation rate, SFR), and nuclear intrinsic bolometric luminosities. We use the individual estimates of AGN bolometric luminosity obtained through SED-fitting decomposition to reconstruct the redshit evolution of the AGN bolometric luminosity function since z∼3. The resulting trends are used to estimate the overall AGN accretion rate density at different cosmic epochs and to trace the first ever estimate of the AGN accretion history from an IR survey. Later on, we focus our study on the connection between AGN accretion and integrated galaxy properties. We analyse the relationships of AGN accretion with galaxy properties in the SFR-M* plane and at different cosmic epochs. Finally, we infer what is the parameter that best correlates with AGN accretion, comparing our results with previous studies and discussing their physical implications in the context of current scenarios of AGN/galaxy evolution.