The clouds of gas and dust that constitute the Interstellar Medium (ISM) within the Milky Way and other galaxies can be studied through the spectral lines of the atoms and molecules. The ISM follows a lifecycle in which each of its phases can be traced through spectral lines in the Terahertz (THz) portion of the electromagnetic spectrum, loosely defined as 0.3 - 3 THz. Using the high spectral resolution afforded by heterodyne instruments, astronomers can potentially disentangle the large-scale structure and kinematics within these clouds. In order to study the ISM over large size scales, large format THz heterodyne arrays are needed. The research presented in this dissertation focuses on the development of two heterodyne array receiver systems for ISM studies, SuperCam and a Super-THz (>3 THz) receiver. SuperCam is a 64-pixel heterodyne imaging array designed for use on ground-based submillimeter telescopes to observe the astrophysically important CO J=3-2 emission line at 345 GHz. The SuperCam focal plane stacks eight, 1x8 mixer subarrays. Each pixel in the array has its own integrated superconductor-insulator-superconductor (SIS) mixer and Low Noise Amplifier (LNA). In spring 2012, SuperCam was installed on the University of Arizona Submillimeter Telescope (SMT) for its first engineering run with 32 active pixels. A second observing run in May 2013 had 52 active pixels. With the outliers removed, the median double sideband receiver temperature was 104 K. The Super-THz receiver is designed to observe the astrophysically important neutral atomic oxygen line at 4.7448 THz. The local oscillator is a third-order distributed feedback Quantum Cascade Laser operating in continuous wave mode at 4.741 THz. A quasi-optical hot electron bolometer is used as the mixer. We record a double sideband receiver noise temperature of 815 K, which is ~7 times the quantum noise limit and an Allan variance time of 15 seconds at an effective noise fluctuation bandwidth of 18 MHz. Heterodyne performance is confirmed by measuring a methanol line spectrum. By combining knowledge of large array formats from SuperCam and quasi-optical mixers, initial tests and designs are presented to expand the single pixel 4.7 THz receiver into a quasi-optical 16-pixel array.

I investigate the emission properties of the molecular interstellar medium in protoplanetary disks and galaxy mergers, though focus largely on the latter topic. I utilize both numerical models as well as observations to relate the emission characteristics to physical models for the formation and evolution of gas giant planets and galaxies. The main results of this thesis follow. (1) Gas giant protoplanets may be detectable via self-absorption signatures in molecular emission lines with sufficiently high critical density. Given the spatial resolution of e.g. ALMA, gas giant planets in formation may be directly imageable. (2) Starburst and AGN feedback-driven winds in galaxies can leave imprints on the molecular line emission properties via morphological outflows and high velocity peaks in the emission line spectra. Methods for distinguishing between high velocity peaks driven by dynamics versus those driven by winds are discussed. (3) CO line widths on average trace the virial velocity of z ∼ 6 quasar host halos. Thus, if the earliest quasars formed in ∼1013 M ⊙ halos, they are predicted to have broad molecular line widths. Selection effects may exist which tend quasars selected for optical luminosity toward molecular line widths narrower than the slightline-dependent mean. (4) Using the SMT, I observe a roughly linear relation between infrared luminosity and CO (J=3-2) luminosity in local galaxies confirming the results of recently observed L(IR)-HCN (J=1-0) relations. Subsequent modeling shows that observed SFR-molecular line luminosity relations owe to the average fraction of subthermally excited gas in galaxies, and are simply reflective of the assumed Schmidt law governing the SFR.

Moving groups are associations of stars which originated from the same star forming region. These groups are typically young (< 200 Myr) since they have not dissipated into the galactic field population. Over the last 15 years, roughly 10 such moving groups have been found with distances < 150 pc (7 with distances < 100 pc), each with a unique velocity and position. This work first investigates the likelihood to resolve star from two moving groups (AB Doradus and Beta Pictoris) using high spacial resolution optical interferrometry and found 5 AB Doradus stars and 1 Beta Pictoris star with declinations > -30 could be spacially resolved. To more deeply characterize individual groups, we used the 2.7m telescope at the McDonald Observatory to observe 10 proposed AB Doradus stars and 5 proposed Octans-Near stars (3 probable members, 2 possible) with high resolution (R ~60,000) optical spectroscopy. Each group is characterized in three ways: (1) Chemical analysis to determine the homogeneity among members, (2) Kinematic traceback to determine the origin, and (3) Isochrone fitting to determine the age. We find the 8 stars in our AB Doradus sample are chemically homogeneous with [M/H] = -0.03 ± 0.06 dex, traceback to an age of 125 Myr, and the stars in this mass range are on the main sequence. The two deviants are a metal rich, potentially younger member and a metal poor, young star likely not associated with AB Doradus. In our Octans-Near sample, we find the 3 probable members have [M/H] = -0.06 ± 0.11, the stars do not trace back to a common origin, and the probable members are on the main sequence. In addition to these tests, we found that the probable members are slightly more lithium depleted than the Pleiades, implying an age between 125 and 200 Myr. Finally, we investigate systematic trends in fundamental stellar parameters from the use of different techniques. Preliminary results find differences in temperatures between interferrometric and spectroscopic techniques to be a function of temperature with a interferrometric temperatures being cooler by an average of 36 ± 115 K. We also calculated the chemical abundances as a function of condensation temperature for our moving group sample and predict 2 stars in AB Doradus could represent the initial star forming environment and discuss the implications for planet hosting stars in nearby moving groups. This updated characterization technique allows for a deeper understanding of the moving group environment. As future, high precision instruments emerge in astronomy (Jame Webb Space Telescope, GAIA, 30m class telescopes), moving groups are ideal targets since these associations will help us understand star forming regions, stellar evolution at young ages, constrain stellar evolutionary models, and identify planetary formation and evolution mechanisms.

The construction of the James Webb Space Telescope has brought attention to infrared astronomy and cosmology. The potential information about our universe to be gained by this mission and future infrared telescopes is staggering, but infrared observation faces many obstacles. These telescopes face large amounts of noise by many phenomena, from emission off of the mirrors to the cosmic infrared background. Infrared telescopes need to be designed in such a way that noise is minimized to achieve sufficient signal to noise ratio on high redshift objects. We will investigate current and planned space and ground based telescopes, model the noise they encounter, and discover their limitations. The ultimate goal of our investigation is to compare the sensitivity of these missions in the near and mid IR and to propose new missions. Our investigation is broken down into four major sections: current missions, noise, signal, and proposed missions. In the proposed missions section we investigate historical and current infrared telescopes with attention given to their location and properties. The noise section discusses the noise that an infrared telescope will encounter and set the background limit. The signal section will look at the spectral energy distributions (SED) of a few significant objects in our universe. We will calculate the intensity of the objects at various points on Earth and in orbit. In the final section we use our findings in the signal and noise sections to model integration times (observation time) for a variety of missions to achieve a given signal to noise ratio (SNR).

A two dimensional hydrodynamic model describing the solar wind interaction with the local interstellar medium, which surrounds the solar system, is used to study the heliosphere both as a steady-state- and dynamic structure. The finite volume method used to solve the associated system of hydrodynamic equations numerically is discussed in detail. Subsequently the steady state heliosphere is studied for both the case where the solar wind and the interstellar medium are assumed to consist of protons only, as well as the case where the neutral hydrogen population in the interstellar medium is taken into account. It is shown that the heliosphere forms as three waves, propagating away from the initial point of contact between the solar wind and interstellar matter, become stationary. Two of these waves become stationary at sonic points, forming the termination shock and bow shock respectively. The third wave becomes stationary as a contact discontinuity, called the heliopause. It is shown that the position and geometry of the termination shock, heliopause and bow shock as well as the plasma flow characteristics of the heliosphere largely depend on the dynamic pressure of either the solar wind or interstellar matter. The heliosphere is modelled as a dynamic structure, including both the effects of the solar cycle and short term variations in the solar wind observed by a range of spacecraft over the past ~ 30 years. The dynamic model allows the calculation of an accurate record of the heliosphere state over the past ~ 30 years. This record is used to predict the time at which the Voyager 2 spacecraft will cross the termination shock. Voyager 1 observations of 10 MeV cosmic ray electrons are then used in conjunction with a cosmic ray modulation model to constrain the record of the heliosphere further. It is shown that the dynamic hydrodynamic model describes the heliosphere accurately within a margin of error of ±0.7 years and ±3 AU. The model predicts that Voyager 2 crossed the termination shock in 2007, corresponding to preliminary results from observations indicating that the crossing occurred in August 2007. / Thesis (M.Sc. (Physics))--North-West University, Potchefstroom Campus, 2008.

This document describes studies of two independent regions of the interstellar medium (ISM). These studies have the common element that both pertain to regions in our Galaxy that are known to be associated with present-day star formation. These studies aim to help us understand the ISM, star formation, and ultimately where we came from, since, after all, our star, the Sun, is itself the product of star formation 4.5 billion years ago. The first project measured the Zeeman Effect on the 21 cm H I absorption line in order to create a map of the line of sight magnetic fields near a star forming region called W3. From the map of the field, it was possible to create a three dimensional model of the magnetic field morphology. Also, calculating the various energies associated with W3 revealed that it is most likely in virial equilibrium, not expanding or contracting. The second project used an instrument on the Hubble Space Telescope (HST) to measure the abundance of iron in a neutral region near the Orion Nebula called Orion’s Veil. One of the goals of this project is to investigate whether solid dust grains can be destroyed by ionizing radiation by comparing the amount of solid iron in Orion’s Veil to the amount in the nearby ionized regions. By measuring the depletion of iron in the neutral Veil and comparing it to the depletion of iron the H+ regions of the Orion Nebula, it was possible to conclude that iron was not being released into the gas phase by ultraviolet photons. In addition, oscillator strengths for two Fe II transitions were measured.

RR Lyrae stars are old, iron-poor, Helium-burning variable stars. RR Lyraes are extremely useful for tracing phase-space structures and metallicities within the galaxy because they are easy to identify, have consistent luminosities, and are found in large numbers in the galactic disk, bulge, and halo. Here we present a new set of spectroscopic metallicity calibrations that use the equivalent widths of the Ca II K, Hγ, and Hδ lines to calculate metallicity values. Applied to spectroscopic survey data, these calibrations will help shed light on the evolution of the Milky Way and other galaxies.

A two dimensional hydrodynamic model describing the solar wind interaction with the local interstellar medium, which surrounds the solar system, is used to study the heliosphere both as a steady-state- and dynamic structure. The finite volume method used to solve the associated system of hydrodynamic equations numerically is discussed in detail. Subsequently the steady state heliosphere is studied for both the case where the solar wind and the interstellar medium are assumed to consist of protons only, as well as the case where the neutral hydrogen population in the interstellar medium is taken into account. It is shown that the heliosphere forms as three waves, propagating away from the initial point of contact between the solar wind and interstellar matter, become stationary. Two of these waves become stationary at sonic points, forming the termination shock and bow shock respectively. The third wave becomes stationary as a contact discontinuity, called the heliopause. It is shown that the position and geometry of the termination shock, heliopause and bow shock as well as the plasma flow characteristics of the heliosphere largely depend on the dynamic pressure of either the solar wind or interstellar matter. The heliosphere is modelled as a dynamic structure, including both the effects of the solar cycle and short term variations in the solar wind observed by a range of spacecraft over the past ~ 30 years. The dynamic model allows the calculation of an accurate record of the heliosphere state over the past ~ 30 years. This record is used to predict the time at which the Voyager 2 spacecraft will cross the termination shock. Voyager 1 observations of 10 MeV cosmic ray electrons are then used in conjunction with a cosmic ray modulation model to constrain the record of the heliosphere further. It is shown that the dynamic hydrodynamic model describes the heliosphere accurately within a margin of error of ±0.7 years and ±3 AU. The model predicts that Voyager 2 crossed the termination shock in 2007, corresponding to preliminary results from observations indicating that the crossing occurred in August 2007. / Thesis (M.Sc. (Physics))--North-West University, Potchefstroom Campus, 2008.

26Al(p,p)26Al and 26Al(p,p’)26Al* scattering reactions were performed at the Holifield Radioactive Ion Beam Facility (HRIBF) at the Oak Ridge National Laboratory (ORNL). The purpose of the elastic scattering study was to determine properties of previously uncharacterized 27Si levels above the proton threshold in the energy range E(c.m.) ~ 0.5 - 1.5 MeV and to calculate reaction rates for the 26Al(p,γ[gamma])27Si reaction that destroys 26Al. The inelastic scattering reaction was also evaluated to investigate the reaction that produces the metastable state of 26Al at E(c.m.) = 228 keV, which would in turn destroy 26Al in the stellar environment. Pure 26Al beams (E(beam) = 13 - 41 MeV) with intensities of ~2*106 26Al/s bombarded a thin polypropylene target of 46 μ[micro]g/cm2 thickness for 7 days. Scattered protons were detected in the Silicon Detector Array (SIDAR), covering laboratory angles 18 to 41 degrees. Background events were rejected by detecting these protons in coincidence with recoiled 26Al particles in an ionization chamber, and proton yields were measured at 45 energies from E(c.m.) = 0.49 - 1.53 MeV. A thick 2.4 mg/cm2 polypropylene target was also bombarded with a 32 MeV 26Al beam for 1.5 days for comparison with thin-target excitation functions. Little evidence for the inelastic scattering reaction was observed, indicating that this is not a significant destruction pathway. For the first time, however, an upper limit for the cross section of this reaction was estimated, and it has been set at 5*10-2 barns. The first upper limits were also established for possible resonances of the elastic scattering reaction with angular momentum transfers up to L = 3 that were not directly observed by this study. Thin-target elastic scattering data suggested a potential resonance at E(r) = 544 keV, which had not been previously observed, with (9/2, 11/2)+ spin and proton width Γp[Gamma_p] ≤ 1 keV. Thick-target analysis appeared to confirm this result. An upper limit for the strength of this resonance was estimated to be 1.4*10-5 keV or 1.6*10-5 keV for a 9/2+ or 11/2+ state, respectively, moderately increasing the total 26Al(p,γ[gamma])27Si resonant reaction rate at supernova temperatures.

A new methodology is developed to calculate density of states of interstellar and atmospheric clusters that takes account of their loosely bound nature and incorporates kinetically important angular momentum constraints explicitly. The method is based on classical phase space integration for the intermonomer modes of the cluster with imposition of the constraints of selected total energy and total angular momentum. It achieves considerable efficiency via essentially analytic evaluation of the momentum space integrals coupled with efficient Monte Carlo sampling of configurations. The derivation for the equation for the density of states is outlined and all steps in the simplification of the accessible momentum space volume are detailed. The method is tested rigorously against an entirely analytic result obtained for the ideal case of a dimer with spherical top fragments and no interaction potential. Interstellar applications of the new approach are presented for (HCN)2 and (CO)2. The new intermononmer density of states has been integrated over metastable states to obtain the intermonomer partition function, which in turn is used to calculate the metastable equilibrium constants for interstellar clusters, which in turn is used tocalculate the second order rate constant of overall dimer formation in the interstellar environment. Atmospheric applications of the new approach are presented for (H2O)2. The new intermonomer density of states is convoluted with the intramonomer density of states to obtain the convoluted density of states. This convoluted density of states is then integrated over total energy and angular momentum to obtain the anharmonic partition function, which in turn is used to calculate the equilibrium constant for atmospheric clusters, which in turn is used to calculate the third order rate constant for overall dimer formation in the atmospheric environment. Kinetic quantities are also calculated with the intermonomer and convoluted density of states for interstellar and atmospheric clusters, respectively. These densities of states are combined with RRKM theory to compute unimolecular dissociation rate constants, which are then averaged with respect to the thermal capture flux distribution to compute average lifetimes as a function of temperature.