We present high-spectral resolution observations from 12.96 - 13.33 microns towards Orion IRc2 using the mid-infrared spectrograph, EXES, on SOFIA. These observations probe the physical and chemical conditions of the Orion hot core, which is sampled by a bright, compact, mid-infrared background continuum source in the region, IRc2. All ten of the rovibrational C2H2 transitions expected in our spectral coverage, are detected with high S/N, yielding continuous coverage of the R-branch lines from J=9-8 to J=18-17, including both ortho and para species. Eight of these rovibrational transitions are newly reported detections. The isotopologue, 13CCH2, is clearly detected with high signal-to-noise. This enabled a direct measurement of the 12C/13C isotopic ratio for the Orion hot core of 14 +/- 1 and an estimated maximum value of 21. We also detected several HCN rovibrational lines. The ortho and para C2H2 ladders are clearly separate and tracing two different temperatures, 226 K and 164 K, respectively, with a non-equilibrium ortho to para ratio (OPR) of 1.7 +\- 0.1. Additionally, the ortho and para V_LSR values differ by about 1.8 +/- 0.2 km/s, while, the mean line widths differ by 0.7 +/- 0.2 km/s, suggesting that these species are not uniformly mixed along the line of sight to IRc2. We propose that the abnormally low C2H2 OPR could be a remnant from an earlier, colder phase, before the density enhancement (now the hot core) was impacted by shocks generated from an explosive event 500 yrs ago.

Accurate modeling of physical and chemical processes in the interstellar medium requires detailed knowledge of how atoms and molecule adsorb on dust grains. However, the sticking coefficient, a number between 0 and 1 that measures the first step in the interaction of a particle with a surface, is usually assumed in simulations of ISM environments to be either 0.5 or 1. Here we report on the determination of the sticking coefficient of H$_2$, D$_2$, N$_2$, O$_2$, CO, CH$_4$, and CO$_2$ on non-porous amorphous solid water (np-ASW). The sticking coefficient was measured over a wide range of surface temperatures using a highly collimated molecular beam. We showed that the standard way of measuring the sticking coefficient --- the King-Wells method --- leads to the underestimation of trapping events in which there is incomplete energy accommodation of the molecule on the surface. Surface scattering experiments with the use of a pulsed molecular beam are used instead to measure the sticking coefficient. Based on the values of the measured sticking coefficient we suggest a useful general formula of the sticking coefficient as a function of grain temperature and molecule-surface binding energy. We use this formula in a simulation of ISM gas-grain chemistry to find the effect of sticking on the abundance of key molecules both on grains and in the gas-phase.

We measured the binding energy of N$_2$, CO, O$_2$, CH$_4$, and CO$_2$ on non-porous (compact) amorphous solid water (np-ASW), of N$_2$ and CO on porous amorphous solid water (p-ASW), and of NH$_3$ on crystalline water ice. We were able to measure binding energies down to a fraction of 1\% of a layer, thus making these measurements more appropriate for astrochemistry than the existing values. We found that CO$_2$ forms clusters on np-ASW surface even at very low coverages. The binding energies of N$_2$, CO, O$_2$, and CH$_4$ decrease with coverage in the submonolayer regime. Their values at the low coverage limit are much higher than what is commonly used in gas-grain models. An empirical formula was used to describe the coverage dependence of the binding energies. We used the newly determined binding energy distributions in a simulation of gas-grain chemistry for cold cloud and hot core models. We found that owing to the higher value of desorption energy in the sub-monlayer regime a fraction of all these ices stays much longer and up to higher temperature on the grain surface compared to the single value energies currently used in the astrochemical models.

According to traditional gas-phase chemical models, O2 should be abundant in molecular clouds, but until recently, attempts to detect interstellar O2 line emission with ground- and space-based observatories have failed. Following the multi-line detections of O2 with low abundances in the Orion and rho Oph A molecular clouds with Herschel, it is important to investigate other environments, and we here quantify the O2 abundance near a solar-mass protostar. Observations of O2, at 487 GHz toward a deeply embedded low-mass Class 0 protostar, NGC 1333-IRAS 4A, are presented, using the HIFI instrument on the Herschel Space Observatory. Complementary data of the chemically related NO and CO molecules are obtained as well. The high spectral resolution data are analysed using radiative transfer models to infer column densities and abundances, and are tested directly against full gas-grain chemical models. The deep HIFI spectrum fails to show O2 at the velocity of the dense protostellar envelope, implying one of the lowest abundance upper limits of O2/H2 at <6x10^-9 (3 sigma). However, a tentative (4.5 sigma) detection of O2 is seen at the velocity of the surrounding NGC 1333 molecular cloud, shifted by 1 km/s relative to the protostar. For the protostellar envelope, pure gas-phase models and gas-grain chemical models require a long pre-collapse phase (~0.7-1x10^6 years), during which atomic and molecular oxygen are frozen out onto dust grains and fully converted to H2O, to avoid overproduction of O2 in the dense envelope. The same model also reproduces the limits on the chemically related NO molecule. The tentative detection of O2 in the surrounding cloud is consistent with a low-density PDR model with small changes in reaction rates. The low O2 abundance in the collapsing envelope around a low-mass protostar suggests that the gas and ice entering protoplanetary disks is very poor in O2.

Evolution of grain mantles in various interstellar environment is studied. We concentrate mainly on water, methanol, carbon di-oxide, which constitute nearly 90% of the grain mantle. We investigate how the production rates of these molecules depend on the relative gas phase abundances of oxygen and carbon monoxide and constrain the relevant parameter space which reproduces these molecules closed to the observed abundances. Allowing to accrete only H, O and CO on the grains and using the Monte-Carlo method we follow the chemical processes for a few million years. We allow formation of multi-layers on the grains and incorporate the freeze-out effects of accreting O and CO. We find that the formation of these molecules depends on the initial conditions as well as the average cloud density. Specifically, when the number density of accreting O is less than 3 times more than that of CO, methanol is always over-produced. Using available reaction pathways it appears to be difficult to match the exact observed abundances of all the three molecules simultaneously. Only in a narrow region of parameter space all these three molecules are produced within the observed limit. In addition to this, we found that the incorporation of the freeze-outs of O and CO leads to almost steady state on the grain surface. The mantle thickness grows anywhere between 60 to 500 layers in a period of two million years. In addition, we consider a case where the gas number density changes with time due to gradual collapse of the molecular cloud and present the evolution of composition of different species as a function of radius of the collapsing cloud.

We study the formation of water and methanol in the dense cloud conditions to find the dependence of its production rate on the binding energies, reaction mechanisms, temperatures, and grain site number. We wish to find the effective grain surface area available for chemical reaction and the effective recombination timescales as functions of grain and gas parameters. We used a Monte Carlo simulation to follow the chemical processes occurring on the grain surface. We find that the formation rate of various molecules is strongly dependent on the binding energies. When the binding energies are high, it is very difficult to produce significant amounts of the molecular species. Instead, the grain is found to be full of atomic species. The production rates are found to depend on the number density in the gas phase. We show that the concept of the effective grain surface area, which we introduced in our earlier work, plays a significant role in grain chemistry. We compute the abundance of water and methanol and show that the results strongly depend on the density and composition in the gas phase, as well as various grain parameters. In the rate equation, it is generally assumed that the recombination efficiencies are independent of the grain parameters, and the surface coverage. Presently, our computed parameter $\alpha$ for each product is found to depend on the accretion rate, the grain parameters and the surface coverage of the grain. We compare our results obtained from the rate equation and the one from the effective rate equation, which includes $\alpha$. At the end we compare our results with the observed abundances.

In the interstellar clouds, molecular hydrogens are formed from atomic hydrogen on grain surfaces. An atomic hydrogen hops around till it finds another one with which it combines. This necessarily implies that the average recombination time, or equivalently, the effective grain surface area depends on the relative numbers of atomic hydrogen influx rate and the number of sites on the grain. Our aim is to discover this dependency. We perform a numerical simulation to study the recombination of hydrogen on grain surfaces in a variety of cloud conditions. We use a square lattice (with a periodic boundary condition) of various sizes on two types of grains, namely, amorphous carbon and olivine. We find that the steady state results of our simulation match very well with those obtained from a simpler analytical consideration provided the `effective' grain surface area is written as $\sim S^{\alpha}$, where, $S$ is the actual physical grain area and $\alpha$ is a function of the flux of atomic hydrogen which is determined from our simulation. We carry out the simulation for various astrophysically relevant accretion rates. For high accretion rates, small grains tend to become partly saturated with $H$ and $H_2$ and the subsequent accretion will be partly inhibited. For very low accretion rates, the number of sites to be swept before a molecular hydrogen can form is too large compared to the actual number of sites on the grain, implying that $\alpha$ is greater than unity.

We study the formation and evolution of several molecules in a collapsing interstellar cloud using a reasonably large reaction network containing more then four hundred atomic and molecular species. We employ a time dependent, spherically symmetric, hydrodynamics code to follow the hydrodynamic and chemical evolution of the collapsing cloud. The flow is assumed to be self-gravitating. We use two models to study the hydrodynamic evolution: in the first model, we inject matter into an initially low density region and in the second model, we start with a constant density cloud and let it collapse due to self-gravity. We study the evolution of the central core for both the cases. We include the grain chemistry to compute the formation of molecular hydrogen and carried out the effect of gas and grain chemistry at each time step. We follow the collapse for more than $10^{14}$s (about 3 million years) and present the time evolution of the globally averaged abundances of various simple but biologically important molecules, such as glycine, alanine etc. We compare our results with those obtained from observations found that for lighter molecules the agreement is generally very good. For complex molecules we tend to under predict the abundances. This indicates that other pathways could be present to form these molecules or more accurate reaction rates were needed. Keywords: hydrodynamics; star formation; ISM; chemical evolution

A knowledge of the recombination time on the grain surfaces has been a major obstacle in deciding the production rate of molecular hydrogen and other molecules in the interstellar medium. We present a numerical study to compute this time for molecular hydrogen for various cloud and grain parameters. We also find the time dependence, particularly when a grain is freshly injected into the system. Apart from the fact that the recombination times seem to be functions of the grain parameters such as the activation barrier energy, temperature etc, our result also shows the dependence on the number of sites in the grain $S$ and the effective accretion rate per site $a_s$ of atomic hydrogen. Simply put, the average time that a pair of atomic hydrogens will take to produce one molecular hydrogen depends on how heavily the grain is already populated by atomic and molecular hydrogens and how fast the hopping and desorption times are. We show that if we write the average recombination time as $T_r \sim S^\alpha/A_H$, where, $A_H$ is the hopping rate, then $\alpha$ could be much greater than 1 for all astrophysically relevant accretion rates. Thus the average formation rate of $H_2$ is also dependent on the grain parameters, temperature and the accretion rate. We believe that our result will affect the overall rate of the formation of complex molecules such as methanol which require successive hydrogenation on the grain surfaces in the interstellar medium.

Solid O2 has been proposed as a possible reservoir for oxygen in dense clouds through freeze-out processes. The aim of this work is to characterize quantitatively the physical processes that are involved in the desorption kinetics of CO-O2 ices by interpreting laboratory temperature programmed desorption (TPD) data. This information is used to simulate the behavior of CO-O2 ices under astrophysical conditions. The TPD spectra have been recorded under ultra high vacuum conditions for pure, layered and mixed morphologies for different thicknesses, temperatures and mixing ratios. An empirical kinetic model is used to interpret the results and to provide input parameters for astrophysical models. Binding energies are determined for different ice morphologies. Independent of the ice morphology, the desorption of O2 is found to follow 0th-order kinetics. Binding energies and temperature-dependent sticking probabilities for CO-CO, O2-O2 and CO-O2 are determined. O2 is slightly less volatile than CO, with binding energies of 912+-15 versus 858+-15 K for pure ices. In mixed and layered ices, CO does not co-desorb with O2 but its binding energies are slightly increased compared with pure ice whereas those for O2 are slightly decreased. Lower limits to the sticking probabilities of CO and O2 are 0.9 and 0.85, respectively, at temperatures below 20K. The balance between accretion and desorption is studied for O2 and CO in astrophysically relevant scenarios. Only minor differences are found between the two species, i.e., both desorb between 16 and 18K in typical environments around young stars. Thus, clouds with significant abundances of gaseous CO are unlikely to have large amounts of solid O2.

We report our first results of monitoring sudden ionospheric disturbances (SID). We present data taken continuously for two weeks during 20th Sept. 2002 and 4th Oct. 2002. We compare the effects of solar flares of the VLF signal with those obtained by GOES Satellite of NASA monitoring in X-rays and found excellent agreement.

Using a Gyrator-II Loop antenna tuned at 19.0Khz, we monitored the meteor shower during 17-24th November, 2002. We observe the primary peak at 3h58m (UT) on the 19th of November, 2002. We distinctly observed several `beadlike' and `exponentially dropping' signals. The `beadlike' signals were more in abundance on the 18th of November, 2002, one day prior to the actual encounter.

We present results of several numerical simulations of two dimensional advective flows which include cooling processes. We show that the computed light curve is similar to the $\chi$ state in GRS 1915+105. The power density spectrum (PDS) also shows presence of QPOs near the break frequency.

We study the evolution of molecular hydrogen on the grain surfaces and in the gas phase using both the rate equation (which tracks the average number of molecules) and the master equation (which tracks the expectation values of molecules). We show that above a certain critical accretion rate of $H$ on the grains, the results from these two methods become identical. We used this result to follow the collapse of a dense interstellar cloud and studied the formation of molecular hydrogen for two different temperatures (T=10K and 12K) and two different masses ($1M_\odot$ and $10M_\odot$) of the cloud when olivine grains were used. Since at higher temperatures, the recombination is very small for these grains, we also studied a similar hydrodynamic processes at higher temperatures (T=20K and 25K) when amorphous carbon grains were used. We find that generally, for olivine grains, more than 90% $H$ is converted to $H_2$ within $\sim 10^{5-7}$yr whereas for amorphous grains it takes $\sim 10^{6-7}$yr. $H_2$ formed in this manner can be adequate to produce the observed complex molecules.

We present the results of several numerical simulations of two dimensional axi-symmetric accretion flows around black holes using Smoothed Particle Hydrodynamics (SPH) in the presence of cooling effects. We consider both stellar black holes and super-massive black holes. We observe that due to {\it both radial and vertical oscillation} of shock waves in the accretion flow, the luminosity and average thermal energy content of the inner disk exhibit very interesting behaviour. When power density spectra are taken, quasi-periodic variabilities are seen at a few Hz and also occasionally at hundreds of Hz for stellar black holes. For super-massive black holes, the time scale of the oscillations ranges from hours to weeks. The power density spectra have a flat top behavior with average {\it rms} amplitude of a few percent and a broken power-law behavior. The break frequency is generally found to be close to the frequency where the shock oscillates.

We present results of several numerical simulations of two dimensional axi-symmetric accretion flows around black holes using the Smoothed Particle Hydrodynamics (SPH). We consider both stellar black holes and as well as super-massive black holes. We assume bremsstrahlung to be the only source of cooling as it is simpler to implement in numerical simulations. We observe that due to both radial and vertical oscillation of shock waves in the accretion flow, the luminosity and average thermal energy content of the inner disk exhibit very interesting behaviors. When power density spectra are taken, quasi-periodic oscillations are seen at a few Hz and also occasionally at hundreds of Hz for stellar black holes. For super-massive black holes, the time-scale of the oscillations ranges from hours to weeks. The power density spectra have usual flat top behavior with average {\it rms} amplitude a few percent and a broken power-law behavior. The break frequency is generally found to be close to the QPO frequency as seen in the observed power spectra of black holes.