The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

LDN 1157, is one of the several clouds situated in the cloud complex, LDN 1147/1158, represents a coma-shaped morphology with a well-collimated bipolar outflow emanating from a Class 0 protostar, LDN 1157-mm. The main goals of this work are (a) to map the inter-cloud magnetic field (ICMF) geometry of the region surrounding LDN 1157 to investigate its relationship with the cloud morphology, with the outflow direction and with the core magnetic field (CMF) geometry inferred from the mm- and sub-mm polarization results from the literature, and (b) to investigate the kinematic structure of the cloud. We carried out R-band polarization observations of the stars projected on the cloud to map the pc-scale magnetic field geometry and made spectroscopic observations of the entire cloud in 12CO, C18O and N2H+ (J=1-0) lines to investigate its kinematic structure. We obtained a distance of 340$\pm$3 pc to the LDN 1147/1158, complex based on the Gaia DR2 parallaxes and proper motion values of the three YSOs associated with the complex. A single filament of $\sim1.2$ pc in length and $\sim0.09$ pc in width is found to run all along the coma-shaped cloud. Based on the relationships between the ICMF, CMF, filament orientations, outflow direction, and the presence of an hour-glass morphology of the magnetic field, it is likely that the magnetic field had played an important role in the star formation process in LDN 1157. Combining the proper motions of the YSOs and the radial velocity of LDN 1147/1158 and another complex LDN 1172/1174 which is situated $\sim2$\dgr~east of it, we found that both the complexes are moving collectively toward the Galactic plane. The filamentary morphology of the east-west segment of LDN 1157 may have formed as a result of mass lost by ablation due to the interaction of the moving cloud with the ambient interstellar medium.

We present a statistical study on the orientation of outflows with respect to large-scale filaments and the magnetic fields. Although filaments are widely observed toward Galactic star-forming regions, the exact role of filaments in star formation is unclear. Studies toward low-mass star-forming regions revealed both preferred and random orientation of outflows respective to the filament long-axes, while outflows in massive star-forming regions mostly oriented perpendicular to the host filaments, and parallel to the magnetic fields at similar physical scales. Here, we explore outflows in a sample of 11 protoclusters in HII regions, a more evolved stage compared to IRDCs, using ALMA CO (3-2) line observations. We identify a total of 105 outflow lobes in these protoclusters. Among the 11 targets, 7 are embedded within parsec-scale filamentary structures detected in $^{13}$CO line and 870 $\mu m$ continuum emissions. The angles between outflow axes and corresponding filaments ($\gamma_\mathrm{Fil}$) do not show any hint of preferred orientations (i.e., orthogonal or parallel as inferred in numerical models) with respect to the position angle of the filaments. Identified outflow lobes are also not correlated with the magnetic fields and Galactic plane position angles. Outflows associated with filaments aligned along the large-scale magnetic fields are also randomly orientated. Our study presents the first statistical results of outflow orientation respective to large-scale filaments and magnetic fields in evolved massive star-forming regions. The random distribution suggests a lack of alignment of outflows with filaments, which may be a result of the evolutionary stage of the clusters.

The supersonic winds produced by massive stars carry a large amount of kinetic power. In numerous scenarios such winds have been proven to produce shocks in which relativistic particles are accelerated emitting non-thermal radiation. Here, we report the first detection of non-thermal emission from a single stellar bubble, G2.4+1.4, associated with a WO star. We observed this source with the uGMRT in Band 4 ($550-850$ MHz) and Band 5 ($1050-1450$ MHz). We present intensity and spectral index maps for this source that are consistent with synchrotron emission (average spectral index, $\alpha = -0.83 \pm 0.10$). The fraction of the available kinetic wind power that is converted into cosmic ray acceleration is estimated to be of the order of a few per cent. This finding constitutes an observational breakthrough and gives new insight on the non-thermal physical processes taking place in the environments of isolated massive stars. In particular, our results show that non-runaway isolated massive stars are capable of accelerating relativistic particles and are therefore confirmed as sources of Galactic cosmic rays.

We present a multiwavelength study of the extended green object, G12.42+0.50 in this paper. The associated ionized, dust, and molecular components of this source are studied in detail employing various observations at near-, mid- and far-infrared, submillimeter and radio wavelengths. Radio continuum emission mapped at 610 and 1390 MHz, using the Giant Meterwave Radio Telescope, India, advocates for a scenario of coexistence of an UC H II region and an ionized thermal jet possibly powered by the massive young stellar object, IRAS 18079-1756 with an estimated spectral type of B1 - B0.5. Shock-excited lines of H2 and [FeII], as seen in the near-infrared spectra obtained with UKIRT-UIST, lend support to this picture. Cold dust emission shows a massive clump of mass 1375 M{\sun} enveloping G12.42+0.50. Study of the molecular gas kinematics using the MALT90 and JCMT archival data unravels the presence of both infall activity and large-scale outflow suggesting an early stage of massive star formation in G12.42+0.50. A network of filamentary features are also revealed merging with the massive clump mimicking a hub-filament layout. Velocity structure along these indicate bulk inflow motion.

We present low frequency wide band observations (300-500 MHz) of the star forming complex G351.7-1.2 using upgraded Giant Metrewave Radio Telescope (uGMRT), India. Combining this with the optical, infrared and submillimeter data, we analyse the large scale diffuse radio emission associated with the region that exhibits a broken shell morphology. The spectral index of the emission in the shell is -0.8, indicating non-thermal emission. H-alpha emission that mimics the morphology of the radio shell on a smaller scale is also detected here. Based on the non-thermal emission from the radio shell and the presence of its optical counterpart, we classify G351.7-1.2 as a candidate SNR. A gamma-ray source detected by Fermi LAT (1FGLJ1729.1-3641c) is located towards the south-west of the radio shell and could have a possible origin in the interaction between high velocity particles from the SNR and the ambient molecular cloud.

Detecting stellar clusters have always been an important research problem in Astronomy. Although images do not convey very detailed information in detecting stellar density enhancements, we attempt to understand if new machine learning techniques can reveal patterns that would assist in drawing better inferences from the available image data. This paper describes an unsupervised approach in detecting star clusters using Deep Variational Autoencoder combined with a Gaussian Mixture Model. We show that our method works significantly well in comparison with state-of-the-art detection algorithm in recognizing a variety of star clusters even in the presence of noise and distortion.

Low radio frequencies are favourable for the identification of emission from non-thermal processes such as synchrotron emission. The massive protostellar jet associated with IRAS 18162-2048 (also known as the HH80-81 system) has been imaged at low radio frequencies: 325, 610 and 1300 MHz, using the Giant Metrewave Radio Telescope, India. This is the first instance of detection of non-thermal emission from a massive protostellar jet at such low radio frequencies. The central region displays an elongated structure characteristic of the jet. In addition, the associated Herbig-Haro objects such as HH80, HH81, HH80N, and other condensations along the inner regions of the jet exhibit negative spectral indices. The spectral indices of most condensations are ~-0.7, higher than the value of -0.3 determined earlier using high frequency measurements. The magnetic field values derived using radio flux densities in the present work, under the assumption of equipartition near minimum energy condition, lie in the range 116-180 microgauss. We probe into the hard X-ray nature of a source that has been attributed to HH80, in an attempt to reconcile the non-thermal characteristics of radio and X-ray measurements. The flux densities of condensations at 610 MHz, measured nearly 11 yrs apart, display variability that could be ascribed to the cooling of condensations, and emphasize the importance of coeval or nearly-coeval measurements for estimation of spectral indices.

We present a multiwavelength study of the filamentary infrared dark cloud (IRDC) G333.73+0.37. The region contains two distinct mid-infrared sources S1 and S2 connected by dark lanes of gas and dust. Cold dust emission from the IRDC is detected at seven wavelength bands and we have identified 10 high density clumps in the region. The physical properties of the clumps such as temperature: 14.3-22.3 K and mass: 87-1530 M_sun are determined by fitting a modified blackbody to the spectral energy distribution of each clump between 160 micron and 1.2 mm. The total mass of the IRDC is estimated to be $~4700 M_sun. The molecular line emission towards S1 reveals signatures of protostellar activity. Low frequency radio emission at 1300 and 610 MHz is detected towards S1 (shell-like) and S2 (compact morphology), confirming the presence of newly formed massive stars in the IRDC. Photometric analysis of near and mid-infrared point sources unveil the young stellar object population associated with the cloud. Fragmentation analysis indicates that the filament is supercritical. We observe a velocity gradient along the filament, that is likely to be associated with accretion flows within the filament rather than rotation. Based on various age estimates obtained for objects in different evolutionary stages, we attempt to set a limit to the current age of this cloud.

We present a spectroscopic follow-up of photometrically-selected young stellar object (YSO) candidates in the Central Molecular Zone of the Galactic center. Our goal is to quantify the contamination of this YSO sample by reddened giant stars with circumstellar envelopes and to determine the star formation rate in the CMZ. We obtained KMOS low-resolution near-infrared spectra (R ~4000) between 2.0 and 2.5 um of sources, many of them previously identified, by mid-infrared photometric criteria, as massive YSOs in the Galactic center. Our final sample consists of 91 stars with good signal-to-noise ratio. We separate YSOs from cool late-type stars based on spectral features of CO and Br_gamma at 2.3 um and 2.16 um respectively. We make use of SED model fits to the observed photometric data points from 1.25 to 24 um in order to estimate approximate masses for the YSOs. Using the spectroscopically identified YSOs in our sample, we confirm that existing colour-colour diagrams and colour-magnitude diagrams are unable to efficiently separate YSOs and cool late-type stars. In addition, we define a new colour-colour criterion that separates YSOs from cool late-type stars in the H-Ks vs H-[8.0] diagram. We use this new criterion to identify YSO candidates in the |l| < 1.5, |b|<0.5 degree region and use model SED fits to estimate their approximate masses. By assuming an appropriate initial mass function (IMF) and extrapolating the stellar IMF down to lower masses, we determine a star formation rate (SFR) of ~0.046 +/- 0.026 Msun/yr assuming an average age of 0.75 +/- 0.25 Myr for the YSOs. This value is lower than estimates found using the YSO counting method in the literature. Our SFR estimate in the CMZ agrees with the previous estimates from different methods and reaffirms that star formation in the CMZ is proceeding at a lower rate than predicted by various star forming models.

Aims. We present a multiwavelength study of two southern Galactic H II regions G346.056-0.021 and G346.077-0.056 which are located at a distance of 10.9 kpc. The distribution of ionized gas, cold and warm dust and the stellar population associated with the two H II regions are studied in detail using measurements at near-infrared, mid-infrared, far-infrared, submillimeter and radio wavelengths. Methods. The radio continuum maps at 1280 and 610 MHz were obtained using the Giant Metrewave Radio Telescope to probe the ionized gas. The dust temperature, column density and dust emissivity maps were generated by using modified blackbody fits in the far-infrared wavelength range 160 - 500 {\mu}m. Various near- and mid-infrared colour and magnitude criteria were adopted to identify candidate ionizing star(s) and the population of young stellar objects in the associated field. Results. The radio maps reveal the presence diffuse ionized emission displaying distinct cometary morphologies. The 1280 MHz flux densities translate to ZAMS spectral types in the range O7.5V - O7V and O8.5V - O8V for the ionizing stars of G346.056-0.021 and G346.077-0.056, respectively. A few promising candidate ionizing star(s) are identified using near-infrared photometric data. The column density map shows the presence of a large, dense dust clump enveloping G346.077-0.056. The dust temperature map shows peaks towards the two H II regions. The submillimetre image shows the presence of two additional clumps one being associated with G346.056-0.021. The masses of the clumps are estimated to range between {\sim} 1400 to 15250 M{\sun}. Based on simple analytic calculations and the correlation seen between the ionized gas distribution and the local density structure, the observed cometary morphology in the radio maps is better explained invoking the champagne-flow model.

A multiwavelength investigation of the southern infrared dust bubble CS51 is presented in this paper. We probe the associated ionized, cold dust, molecular and stellar components. Radio continuum emission mapped at 610 and 1300 MHz, using the Giant Metrewave Radio Telescope, India, reveal the presence of three compact emission components (A, B, and C) apart from large-scale diffuse emission within the bubble interior. Radio spectral index map show the coexistence of thermal and non-thermal emission components. Modified blackbody fits to the thermal dust emission using Herschel PACS and SPIRE data is performed to generate dust temperature and column density maps. We identify five dust clumps associated with CS51 with masses and radius in the range 810 - 4600 M{\sun} and 1.0 - 1.9 pc, respectively. We further construct the column density probability distribution functions of the surrounding cold dust which display the impact of ionization feedback from high-mass stars. The estimated dynamical and fragmentation timescales indicate the possibility of collect and collapse mechanism in play at the bubble border. Molecular line emission from the MALT90 survey is used to understand the nature of two clumps which show signatures of expansion of CS51.

We probe the structure and kinematics of two neighbouring H II regions identified as cometary and bipolar, using radio recombination lines (RRL). The H172{\alpha} RRLs from these H II regions: G351.6-1.15 and G351.6-1.25, are mapped using GMRT, India. We also detect carbon RRLs C172{\alpha} towards both these regions. The hydrogen RRLs display the effects of pressure and dynamical broadening in the line profiles, with the dynamical broadening (15 km/s) playing a major role in the observed profile of G351.6-1.15. We investigate the kinematics of molecular gas species towards this H II region from the MALT90 pilot survey. The molecular gas is mostly distributed towards the north and north-west of the cometary head. The molecular line profiles indicate signatures of turbulence and outflow in this region. The ionized gas at the cometary tail is blue shifted by 8 km/s with respect to the ambient molecular cloud, consistent with the earlier proposed champagne flow scenario. The relative velocity of 5 km/s between the northern and southern lobes of the bipolar H II region G351.6-1.25 is consistent with the premise that the bipolar morphology is a result of the expanding ionized lobes within a flat molecular cloud.

A multi-wavelength investigation of the star forming complex IRAS 20286+4105, located in the Cygnus-X region, is presented here. Near-infrared K-band data is used to revisit the cluster / stellar group identified in previous studies. The radio continuum observations, at 610 and 1280 MHz show the presence of a HII region possibly powered by a star of spectral type B0 - B0.5. The cometary morphology of the ionized region is explained by invoking the bow-shock model where the likely association with a nearby supernova remnant is also explored. A compact radio knot with non-thermal spectral index is detected towards the centre of the cloud. Mid-infrared data from the Spitzer Legacy Survey of the Cygnus-X region show the presence of six Class I YSOs inside the cloud. Thermal dust emission in this complex is modelled using Herschel far-infrared data to generate dust temperature and column density maps. Herschel images also show the presence of two clumps in this region, the masses of which are estimated to be {\sim} 175 M{\sun} and 30 M{\sun}. The mass-radius relation and the surface density of the clumps do not qualify them as massive star forming sites. An overall picture of a runaway star ionizing the cloud and a triggered population of intermediate-mass, Class I sources located toward the cloud centre emerges from this multiwavelength study. Variation in the dust emissivity spectral index is shown to exist in this region and is seen to have an inverse relation with the dust temperature.

Stars and planetary systems are formed out of molecular clouds in the interstellar medium. Although the sequence of steps involved in star formation are generally known, a comprehensive theory which describes the details of the processes that drive formation of stars is still missing. The Square Kilometre Array (SKA), with its unprecedented sensitivity and angular resolution, will play a major role in filling these gaps in our understanding. In this article, we present a few science cases that the Indian star formation community is interested in pursuing with SKA, which include investigation of AU-sized structures in the neutral ISM, the origin of thermal and non-thermal radio jets from protostars and the accretion history of protostars, and formation of massive stars and their effect on the surrounding medium.

We present a multiwavelength investigation of star formation activity towards the southern H II regions associated with IRAS 17160-3707, located at a distance of 6.2 kpc with a bolometric luminosity of 830000 Lsun.The ionised gas distribution and dust clumps in the parental molecular cloud are examined in detail using measurements at infrared, submillimeter and radio wavelengths.The radio continuum images at 1280 and 610 MHz obtained using Giant Metrewave Radio Telescope reveal the presence of multiple compact sources as well as nebulous emission.At submillimeter wavelengths, we identify seven dust clumps and estimate their physical properties like temperature: 24 - 30 K, mass: 300 - 4800 Msun and luminosity: 900 - 31700 Lsun using modified blackbody fits to the spectral energy distributions between 70 and 870 um.We find 24 young stellar objects in the mid-infrared, with few of them coincident with the compact radio sources.The spectral energy distributions of young stellar objects have been fitted by the Robitaille models and the results indicate that those having radio compact sources as counterparts host massive objects in early evolutionary stages with best fit age <= 0.2 Myr.We compare the relative evolutionary stages of clumps using various signposts such as masers, ionised gas, presence of young stellar objects and infrared nebulosity and find six massive star forming clumps and one quiescent clump.Of the former, five are in a relatively advanced stage and one in an earlier stage.

An investigation in radio and infrared wavelengths of two high-mass star forming regions toward the southern Galactic bubble S10 is presented here. The two regions under study are associated with the broken bubble S10 and Extended Green Object, G345.99-0.02, respectively. Radio continuum emission mapped at 610 and 1280 MHz using the Giant Metrewave Radio Telescope, India is detected towards both the regions. These regions are estimated to be ionized by early B to late O type stars. Spitzer GLIMPSE mid-infrared data is used to identify young stellar objects associated with these regions. A Class I/II type source, with an estimated mass of 6.2 M{\sun} , lies {\sim} 7{\arcsec} from the radio peak. Pixel-wise, modified blackbody fits to the thermal dust emission using Herschel far-infrared data is performed to construct dust temperature and column density maps. Eight clumps are detected in the two regions using the 250 {\mu}m image. The masses and linear diameter of these range between {\sim} 300 - 1600 M{\sun} and 0.2 - 1.1 pc, respectively which qualifies them as high-mass star forming clumps. Modelling of the spectral energy distribution of these clumps indicates the presence of high luminosity, high accretion rate, massive young stellar objects possibly in the accelerating accretion phase. Further, based on the radio and MIR morphology, the occurrence of a possible bow-wave towards the likely ionizing star is explored.

We carried out optical polarimetry of an isolated cloud, Gal 110-13, to map the plane-of-the-sky magnetic field geometry. The main aim of the study is to understand the most plausible mechanism responsible for the unusual cometary shape of the cloud in the context of its magnetic field geometry. When unpolarized starlight passes through the intervening interstellar dust grains that are aligned with their short axes parallel to the local magnetic field, it gets linearly polarized. The plane-of-the-sky magnetic field component can therefore be traced by doing polarization measurements of background stars projected on clouds. Because the light in the optical wavelength range is most efficiently polarized by the dust grains typically found in the outer layers of the molecular clouds, optical polarimetry enables us to trace the magnetic field geometry of the outer layers of the clouds. We made R-band polarization measurements of 207 stars in the direction of Gal 110-13. The distance of Gal 110-13 was determined as $\sim450\pm80$ pc using our polarization and 2MASS near-infrared data. The foreground interstellar contribution was removed from the observed polarization values by observing a number of stars located in the vicinity of Gal 110-13 which has Hipparcos parallax measurements. The plane-of-the-sky magnetic field lines are found to be well ordered and aligned with the elongated structure of Gal 110-13. Using structure function analysis, we estimated the strength of the plane-of-the-sky component of the magnetic field as $\sim25\mu$G. Based on our results and comparing them with those from simulations, we conclude that compression by the ionization fronts from 10 Lac is the most plausible cause of the comet-like morphology of Gal 110-13 and of the initiation of subsequent star formation.

IRAS 17256-3631 is a southern Galactic massive star forming region located at a distance of 2 kpc. In this paper, we present a multiwavelength investigation of the embedded cluster, the HII region, as well as the parent cloud. Radio images at 325, 610 and 1372 MHz were obtained using GMRT, India while the near-infrared imaging and spectroscopy were carried out using UKIRT and Mt. Abu Infrared Telescope, India. The near-infrared K-band image reveals the presence of a partially embedded infrared cluster. The spectral features of the brightest star in the cluster, IRS-1, spectroscopically agrees with a late O or early B star and could be the driving source of this region. Filamentary H_2 emission detected towards the outer envelope indicates presence of highly excited gas. The parent cloud is investigated at far-infrared to millimeter wavelengths and eighteen dust clumps have been identified. The spectral energy distributions (SEDs) of these clumps have been fitted as modified blackbodies and the best-fit peak temperatures are found to range from 14-33 K, while the column densities vary from 0.7-8.5x10^22 cm^-2. The radio maps show a cometary morphology for the distribution of ionized gas that is density bounded towards the north-west and ionization bounded towards the south-east. This morphology is better explained with the champagne flow model as compared to the bow shock model. Using observations at near, mid and far-infrared, submillimeter and radio wavelengths, we examine the evolutionary stages of various clumps.

The FU Orionis (FUor) or EX Orionis (EXor) phenomenon has attracted increasing attention in recent years and is now accepted as a crucial element in the early evolution of low-mass stars. FUor and EXor eruptions of young stellar objects (YSOs) are caused by strongly enhanced accretion from the surrounding disk. FUors display optical outbursts of $\sim$ 4 mag or more and last for several decades, whereas EXors show smaller outbursts ($\Delta$m $\sim$ 2 - 3 mag) that last from a few months to a few years and may occur repeatedly. Therefore, FUor/EXor eruptions represent a rare but very important phenomenon in early stellar evolution, during which a young low-mass YSO brightens by up to several optical magnitudes. Hence, long-term observations of this class of eruptive variable are important to design theoretical models of low-mass star formation. In this paper, we present recent results from our long-term monitoring observations of three rare types of eruptive young variables with the 2-m Himalayan {\it Chandra} Telescope (HCT) and the 2-m IUCAA Girawali Observatory (IGO) telescope.

We present a study of MSX point sources in the Galactic bulge (|l| < 3 deg, 1 deg < |b| < 5 deg), observed at A, C, D and E-band (8 to 21 micron), with a total area ~ 48 deg^2. We discuss the nature of the MSX sources (mostly AGB stars), their luminosities, the interstellar extinction, the mass-loss rate distribution and the total mass-loss rate in the bulge. The MIR data of MSX point sources have been combined with the NIR (J, H and Ks) data of 2MASS survey. The cross-identification was restricted to Ks-band detected sources with Ks <= 11 mag. However, for those bright MSX D-band sources ([D] < 4.0 mag), which do not satisfy this criteria, we have set no Ks-band magnitude cut off. The relation between Mdot and (Ks-[15])0 was used to derive the mass-loss rate of each MSX source in the bulge fields. Except for very few post-AGB stars, PNe and OH/IR stars, a large fraction of the detected sources at 15 micron (MSX D-band) are AGB stars well above the RGB tip. A number of them show an excess in ([A]-[D])0 and (Ks-[D])0 colours, characteristic of mass-loss. These colours, especially (Ks-[D])0, enable estimation of the mass-loss rates (Mdot) of the sources in the bulge fields which range from 10^{-7} to 10^{-4} Msun/yr. Taking into consideration the completeness of the mass-loss rate bins, we find that the contribution to the integrated mass-loss is probably dominated by mass-loss rates larger than 3x10^{-7} Msun/yr and is about 1.96 x 10^{-4} Msun/yr/deg^2 in the "intermediate" and "outer" bulge fields of sources with mass-loss rates, Mdot > 3x10^{-7} Msun/yr. The corresponding integrated mass-loss rate per unit stellar mass is 0.48x10^{-11} yr^{-1}. Apart from this, the various MIR and NIR CC and CM diagrams are discussed in the paper to study the nature of the stellar population in the MSX bulge fields.

We present a multiwavelength study of the ultra compact HII region associated with IRAS 20178+4046. This enables us to probe the different components associated with this massive star forming region. The radio emission from the ionized gas was mapped at 610 and 1280 MHz using the Giant Metrewave Radio Telescope (GMRT), India. We have used 2MASS $J H K_{s}$ data to study the nature of the embedded sources associated with IRAS 20178+4046. Submillimetre emission from the cold dust at 450 and 850 $\mu$m was studied using JCMT-SCUBA. The high-resolution radio continuum maps at 610 and 1280 MHz display compact spherical morphology. The spectral type of the exciting source is estimated to be $\sim$ B0.5 from the radio flux densities. However, the near-infrared (NIR) data suggest the presence of several massive stars (spectral type earlier than O9) within the compact ionized region. Submillimetre emission shows the presence of two dense cloud cores which are probably at different evolutionary stages. The total mass of the cloud is estimated to be $\sim$ 700 -- 1500 $\rm M_{\odot}$ from the submillimetre emission at 450 and 850 $\mu$m. The multiwavelength study of this star forming complex reveals an interesting scenario where we see the presence of different evolutionary stages in star formation. The ultra compact HII region coinciding with the southern cloud core is at a later stage of evolution compared to the northern core which is likely to be a candidate protocluster.

We present a detailed study of the post-outburst phase of McNeil's nebula (V1647 Ori) using optical B,V,R,I and NIR J,H,K photometric and low resolution optical spectroscopic observations. The observations were carried out with the HFOSC, NIRCAM, TIRCAM and NICMOS cameras on the 2m HCT and 1.2m PRL telescopes during the period 2004 Feb-2005 Dec. The optical/NIR observations show a general decline in brightness of the exciting source of McNeil's nebula (V1647 Ori). Our recent optical images show that V1647 Ori has faded by more than 3 mags since Feb 2004. The optical/NIR photometric data also show a significant variation in the mags (Delta V = 0.78 mag, Delta R = 0.44 mag, Delta I = 0.21 mag, Delta J = 0.24 mag and Delta H = 0.20 mag) of V1647 Ori within a period of one month, which is possibly undergoing a phase similar to eruptive variables, like EXors or FUors. The optical spectra show a few features such as strong Halpha emission with blue-shifted absorption and the CaII IR triplet (8498A, 8542A and 8662A) in emission. As compared to the period just after outburst, there is a decrease in the depth and extent of the blue-shifted absorption component, indicating a weakening in the powerful stellar wind. The presence of the CaII IR triplet in emission confirms that V1647 Ori is a PMS star. The long-term, post-outburst photometric observations of V1647 Ori suggest an EXor, rather than an FUor event. An optical/IR comparison of the region surrounding McNeil's nebula shows that the optical nebula is more widely and predominantly extended to the north, whereas the IR nebula is relatively confined (dia ~ 60 arcsec), but definitely extended, to the south, too.

We present a multiwavelength study of the massive star forming region associated with IRAS 06055+2039 which reveals an interesting scenario of this complex where regions are at different stages of evolution of star formation. Narrow band near-infrared (NIR) observations were carried out with UKIRT-UFTI in molecular hydrogen and Br$\gamma$ lines to trace the shocked and ionized gases respectively. We have used 2MASS $J H K_{s}$ data to study the nature of the embedded cluster associated with IRAS 06055+2039. We obtain a power-law slope of 0.43$\pm$0.09 for the $K_{s}$-band Luminosity Function (KLF) which is in good agreement with other young embedded clusters. We estimate an age of 2 -- 3 Myr for this cluster. The radio emission from the ionized gas has been mapped at 610 and 1280 MHz using the Giant Metrewave Radio Telescope (GMRT), India. Apart from the diffuse emission, the high resolution 1280 MHz map also shows the presence of several discrete sources which possibly represent high density clumps. The morphology of shocked molecular hydrogen forms an arc towards the N-E of the central IRAS point source and envelopes the radio emission. Submillimetre emission using JCMT-SCUBA show the presence of a dense cloud core which is probably at an earlier evolutionary stage compared to the ionized region with shocked molecular gas lying in between the two. Emission from warm dust and the Unidentified Infrared Bands (UIBs) have been estimated using the mid-infrared (8 -- 21 $\mu$m) data from the MSX survey. From the submillimetre emission at 450 and 850 $\mu$m the total mass of the cloud is estimated to be $\sim$ 7000 -- 9000 $\rm M_{\odot}$.

Pulsation models are crucial for the interpretation of spectrophotometric and interferometric observations of Mira variables. Comparing predicted and observed spectra is one way of establishing the validity of such models. In this paper, we focus on the models published between 1996 and 1998 by Bessell, Hofmann, Scholz and Wood. A few new model spectra are added, to improve available phase coverage. We compare the synthetic spectra with observed low resolution spectra of optically selected oxygen-rich Miras, over a range of optical and near-IR wavelengths that encompasses most of the stellar energy output. We investigate the overall energy distributions, and specific spectral features in the near-IR wavelength range. The agreement between the observed and the model-predicted properties is found to be reasonably good. However, there are discrepancies seen especially in various molecular bands. We find that different combinations of stellar parameters and pulsation phases often result in very similar model spectra. Therefore the problem of deriving parameters of a Mira variable from its spectrum has no unique solution. More advanced models than presently available, providing even better fits to the data and covering a wider range of parameters, would be needed to achieve better discrimination.

Dynamic models of M-type Mira variables predict the occurrence of water "shells", i.e. of zones of high H2O density and high H2O absorption inside the stellar atmosphere. The density, position and width of these shells is closely correlated with different types of two-component shapes of the intensity distribution on the disk in the H, K and L near-continuum bandpasses. We investigate these correlations and highlight the role of a spectrophotometric H2O index that warns against serious complications in diameter measurements in the case of substantial water contamination of the bandpass of observation. Simultaneous spectrophotometric and interferometric measurements may allow observers to estimate real continuum diameters more precisely.

We present a statistical method to derive the mass functions of open clusters using sky survey data such as the 2 Micron All Sky Survey and the Guide Star Catalogue. We have used this method to derive the mass functions in the stellar/substellar regime of three young, nearby open clusters, namely IC 348, Sigma-Orionis and Pleiades. The mass function in the low mass range (M < 0.5 solar mass) is appreciably flatter than the stellar Salpeter function for all three open clusters. The contribution of objects below 0.5 solar mass to the total mass of the cluster is ~40% and the contribution of objects below 0.08 solar mass to the total is ~4%.

In this paper we present the mass functions in the substellar regime of three young open clusters, IC 348, $\sigma$ Orionis and Pleiades, as derived using the data from the 2 Micron All Sky Survey (2MASS) catalogue which has a limiting magnitude of K_s ~ 15, and the latest version of the Guide Star Catalogue (GSC) which has a limiting magnitude of F ~ 21 where F refers to the POSS II IIIa-F passband. Based on recent evolutionary models for low mass stars, we have formulated the selection criteria for stars with masses below 0.5 solar mass. Using a statistical approach to correct for the background contamination, we derive the mass function of objects with masses ranging from 0.5 solar mass down to the substellar domain, well below the Hydrogen Burning Mass Limit. The lowest mass bins in our analysis are 0.025, 0.045 and 0.055 solar masses for IC 348, sigma Orionis and Pleiades, respectively. The resultant slopes of the mass functions are 0.8 +/- 0.2 for IC 348, 1.2 +/- 0.2 for sigma Orionis and 0.5 +/- 0.2 for Pleiades, which are consistent with the previous results. The contribution of objects below 0.5 solar mass to the total mass of the cluster is ~40%, and the contribution of objects below 0.08 solar mass to the total mass is ~4%.