Search

Cultural Astronomy sits at the intersection of our study of the universe and the human experience, exploring how we observe the stars, interpret their motions, and incorporate them into our perspectives of the nature of reality. "Astronomy in Appalachia" is a new Planetarium show, developed for the Planetarium at East Tennessee State University, that explores how the inherited experience of the night sky in Appalachia can be used as a part of Astronomy education, and how that can serve as a cross-cultural point of connection. Using Earth's seasons as an anchor, the show highlights several Appalachian folk science narratives, examining how they are rooted in Astronomy, and relating them to modern scientific interpretations of the cosmos. This article describes the lessons our team learned while researching, scripting, and developing this new program., Comment: four pages, pg. 14-16, 27

The 21 cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation (EoR). This has led to the construction of low-frequency radio interferometric arrays, such as the Hydrogen Epoch of Reionization Array (HERA), aimed at systematically mapping this emission for the first time. Precision calibration, however, is a requirement in 21 cm radio observations. Due to the spatial compactness of HERA, the array is prone to the effects of mutual coupling, which inevitably lead to non-smooth calibration errors that contaminate the data. When unsmooth gains are used in calibration, intrinsically spectrally-smooth foreground emission begins to contaminate the data in a way that can prohibit a clean detection of the cosmological EoR signal. In this paper, we show that the effects of mutual coupling on calibration quality can be reduced by applying custom time-domain filters to the data prior to calibration. We find that more robust calibration solutions are derived when filtering in this way, which reduces the observed foreground power leakage. Specifically, we find a reduction of foreground power leakage by 2 orders of magnitude at k=0.5.

Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategies for mitigating mutual coupling. In this paper, we analyse 12 nights of data from the Hydrogen Epoch of Reionization Array and compare the data against simulations that include a computationally efficient and physically motivated semi-analytic treatment of mutual coupling. We find that simulated coupling features qualitatively agree with coupling features in the data; however, coupling features in the data are brighter than the simulated features, indicating the presence of additional coupling mechanisms not captured by our model. We explore the use of fringe-rate filters as mutual coupling mitigation tools and use our simulations to investigate the effects of mutual coupling on a simulated cosmological 21-cm power spectrum in a "worst case" scenario where the foregrounds are particularly bright. We find that mutual coupling contaminates a large portion of the "EoR Window", and the contamination is several orders-of-magnitude larger than our simulated cosmic signal across a wide range of cosmological Fourier modes. While our fiducial fringe-rate filtering strategy reduces mutual coupling by roughly a factor of 100 in power, a non-negligible amount of coupling cannot be excised with fringe-rate filters, so more sophisticated mitigation strategies are required., Comment: 19 pages, 12 figures, submitted to MNRAS

The long-term brightness evolution of the great comet C/1995 O1 (Hale-Bopp) presented a remarkable opportunity to study the behavior of its coma over four years. We used approximately 2200 total visual magnitudes published in the International Comet Quarterly taken from 17 observers during the period of 1995 July - 1999 September to create a secular lightcurve. In order to account for observer differences, we present a novel algorithm to reduce scatter and increase precision in a lightcurve compiled from many sources. It is implemented in a publicly available code, ICQSPLITTER. This code addresses the differences among observers by using a self-consistent statistical approach, leading to a sharper lightcurve, and improving the precision of the measured slopes. To first order, the comet's lightcurve approximately follows a r$^{-4}$ response for both pre- and post-perihelion distances. Interestingly, the pre-perihelion data are better fit with a fifth-order polynomial with inflection points at 4.0, 2.6, 2.1 and 1.1 au. We analyze these specific regions and find that they are associated with physical phenomena in the comet's evolution. Contrary to other reports, the lightcurve shows no evidence for the comet having been in outburst at discovery. Afrho values derived from the visual lightcurve data are consistent with a r$^{-1.5}$ dependence on heliocentric distance, which is similar in shape to those derived from spectroscopy and narrow-band photometry. We present correlation equations for visual magnitudes and CO and H2O production rates, which are consistent with the pre-perihelion visual magnitudes increasing almost entirely due to CO outgassing until a heliocentric distance of about 2.6 - 3.0 au. We also present two correlation equations that should prove highly useful for observation planning and data analysis, and can be generalized to be applicable to other comets., Comment: 30 pages, 7 figures, 4 tables, submitted to the AAS Journals

High-redshift submillimetre galaxies (SMGs) are some of the most rapidly star-forming galaxies in the Universe. Historically, galaxy formation models have had difficulty explaining the observed number counts of SMGs. We combine a semi-empirical model with 3-D hydrodynamical simulations and 3-D dust radiative transfer to predict the number counts of unlensed SMGs. Because the stellar mass functions, gas and dust masses, and sizes of our galaxies are constrained to match observations, we can isolate uncertainties related to the dynamical evolution of galaxy mergers and the dust radiative transfer. The number counts and redshift distributions predicted by our model agree well with observations. Isolated disc galaxies dominate the faint (S_1.1 < ~1 mJy, or S_850 < ~2 mJy) population. The brighter sources are a mix of merger-induced starbursts and galaxy-pair SMGs; the latter subpopulation accounts for ~30 - 50 per cent of all SMGs at all S_1.1 > ~0.5 mJy (S_850 > ~1 mJy). The mean redshifts are ~3.0 - 3.5, depending on the flux cut, and the brightest sources tend to be at higher redshifts. Because the galaxy-pair SMGs will be resolved into multiple fainter sources by ALMA, the bright ALMA counts should be as much as 2 times less than those observed using single-dish telescopes. The agreement between our model, which uses a Kroupa IMF, and observations suggests that the IMF in high-redshifts starbursts need not be top-heavy; if the IMF were top-heavy, our model would over-predict the number counts. We conclude that the difficulty some models have reproducing the observed SMG counts is likely indicative of more general problems -- such as an under-prediction of the abundance of massive galaxies or a star formation rate--stellar mass relation normalisation lower than that observed -- rather than a problem specific to the SMG population., Comment: Accepted for publication in MNRAS. Insignificant changes from the previous version

We study the motion of a vortex under the influence of a harmonic force in an approximately two dimensional trapped Bose-condensed gas. The Hall-Vinen-Iordanskii equations, modified to include a fluctuational force and an inertial mass term, are solved for the vortex motion. The mass of the vortex has a strong influence on the time it takes the vortex to escape the trap. Since the vortex mass also depends on the trap size we have an additional dependence on the trap size in the escape time which we compare to the massless case., Comment: Submitted to J. Low. Temp. Phys

We study the future orbital evolution and merging of the MW-M31-M33 system, using a combination of collisionless N-body simulations and semi-analytic orbit integrations. Monte-Carlo simulations are used to explore the consequences of varying the initial phase-space and mass parameters within their observational uncertainties. The observed M31 transverse velocity implies that the MW and M31 will merge t = 5.86 (+1.61-0.72) Gyr from now, after a first pericenter at t = 3.87 (+0.42-0.32) Gyr. M31 may (probability p=41%) make a direct hit with the MW (defined here as a first-pericenter distance less than 25 kpc). Most likely, the MW and M31 will merge first, with M33 settling onto an orbit around them. Alternatively, M33 may make a direct hit with the MW first (p=9%), or M33 may get ejected from the Local Group (p=7%). The MW-M31 merger remnant will resemble an elliptical galaxy. The Sun will most likely (p=85%) end up at larger radius from the center of the MW-M31 merger remnant than its current distance from the MW center, possibly further than 50 kpc (p=10%). The Sun may (p=20%) at some time in the next 10 Gyr find itself moving through M33 (within 10 kpc), but while dynamically still bound to the MW-M31 merger remnant. The arrival and possible collision of M31 (and possibly M33) with the MW is the next major cosmic event affecting the environment of our Sun and solar system that can be predicted with some certainty. (Abridged), Comment: 58 pages, 16 figures, to be published in ApJ. Version with high resolution figures and N-body movies available at http://www.stsci.edu/~marel/M31 . Press materials, graphics, and visualizations available at http://hubblesite.org/newscenter/archive/releases/2012/20

(Abridged) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift., Comment: 18 pages, 9 figures, 3 tables

In recent work (arXiv:1101.0002) we have suggested that the high-redshift (z ~ 2-4) bright submillimetre galaxy (SMG) population is heterogeneous, with major mergers contributing both at early stages, where quiescently star-forming discs are blended into one submm source (`galaxy-pair SMGs'), and late stages, where mutual tidal torques drive gas inflows and cause strong starbursts. Here we combine hydrodynamic simulations of major mergers with 3-D dust radiative transfer calculations to determine observational diagnostics that can distinguish between quiescently star-forming SMGs and starburst SMGs via integrated data alone. We fit the far-IR SEDs of the simulated galaxies with the optically thin single-temperature modified blackbody, the full form of the single-temperature modified blackbody, and a power-law temperature-distribution model. The effective dust temperature, T_dust, and power-law index of the dust emissivity in the far-IR, \beta, derived can significantly depend on the fitting form used, and the intrinsic \beta\ of the dust is not recovered. However, for all forms used here, there is a T_dust above which almost all simulated galaxies are starbursts, so a T_dust cut is very effective at selecting starbursts. Simulated merger-induced starbursts also have higher L_IR/M_gas and L_IR/L_FUV than quiescently star-forming galaxies and lie above the star formation rate-stellar mass relation. These diagnostics can be used to test our claim that the SMG population is heterogeneous and to observationally determine what star formation mode dominates a given galaxy population. We comment on applicability of these diagnostics to ULIRGs that would not be selected as SMGs. These `hot-dust ULIRGs' are typically starburst galaxies lower in mass than SMGs, but they can also simply be SMGs observed from a different viewing angle., Comment: 21 pages, 11 figures. Accepted for publication in MNRAS. Minor changes to text but otherwise identical to v1

We present a novel pair of numerical models of the interaction history between the Large and Small Magellanic Clouds (LMC and SMC, respectively) and our Milky Way (MW) in light of recent high precision proper motions (Kallivayalil et al. 2006a,b). Given the new velocities, cosmological simulations of structure formation favor a scenario where the Magellanic Clouds (MCs) are currently on their first infall towards our Galaxy (Boylan-Kolchin et al. 2011, Busha et al. 2011). We illustrate here that the observed irregular morphology and internal kinematics of the MCs (in gas and stars) are naturally explained by interactions between the LMC and SMC, rather than gravitational interactions with the MW. This picture further supports a first infall scenario (Besla et a. 2007). In particular, we demonstrate that the Magellanic Stream, a band of HI gas trailing behind the MCs 150 degrees across the sky, can be accounted for by the action of LMC tides on the SMC before the system was accreted by the MW. We further demonstrate that the off-center, warped stellar bar of the LMC and its one-armed spiral, can be naturally explained by a recent direct collision with the SMC. Such structures are key morphological characteristics of a class of galaxies referred to as Magellanic Irregulars (de Vaucouleurs & Freeman 1972), the majority of which are not associated with massive spiral galaxies. We infer that dwarf-dwarf galaxy interactions are important drivers for the morphological evolution of Magellanic Irregulars and can dramatically affect the efficiency of baryon removal from dwarf galaxies via the formation of extended tidal bridges and tails. Such interactions are important not only for the evolution of dwarf galaxies but also have direct consequences for the buildup of baryons in our own MW, as LMC-mass systems are believed to be the dominant building blocks of MW-type halos., Comment: 33 pages, 21 figures, Accepted for publication in MNRAS, Dec 23 2011

Calculating the galaxy merger rate requires both a census of galaxies identified as merger candidates, and a cosmologically-averaged `observability' timescale T_obs(z) for identifying galaxy mergers. While many have counted galaxy mergers using a variety of techniques, T_obs(z) for these techniques have been poorly constrained. We address this problem by calibrating three merger rate estimators with a suite of hydrodynamic merger simulations and three galaxy formation models. We estimate T_obs(z) for (1) close galaxy pairs with a range of projected separations, (2) the morphology indicator G-M20, and (3) the morphology indicator asymmetry A. Then we apply these timescales to the observed merger fractions at z < 1.5 from the recent literature. When our physically-motivated timescales are adopted, the observed galaxy merger rates become largely consistent. The remaining differences between the galaxy merger rates are explained by the differences in the range of mass-ratio measured by different techniques and differing parent galaxy selection. The major merger rate per unit co-moving volume for samples selected with constant number density evolves much more strongly with redshift (~ (1+z)^(+3.0 \pm 1.1)) than samples selected with constant stellar mass or passively evolving luminosity (~ (1+z)^(+0.1 \pm 0.4)). We calculate the minor merger rate (1:4 < M_{sat}/M_{primary} <~ 1:10) by subtracting the major merger rate from close pairs from the `total' merger rate determined by G-M20. The implied minor merger rate is ~3 times the major merger rate at z ~ 0.7, and shows little evolution with redshift., Comment: accepted by ApJ; 25 pages; 13 figures; 4 tables

Major mergers of disk galaxies are thought to be a substantial driver in galaxy evolution. To trace the fraction and the rate galaxies are in mergers over cosmic times, several observational techniques, including morphological selection criteria, have been developed over the last decade. We apply this morphological selection of mergers to 21 cm radio emission line (HI) column density images of spiral galaxies in nearby surveys. In this paper, we investigate how long a 1:1 merger is visible in HI from N-body simulations. We evaluate the merger visibility times for selection criteria based on four parameters: Concentration, Asymmetry, M20, and the Gini parameter of second order moment of the flux distribution (GM). Of three selection criteria used in the literature, one based on Concentration and M20 works well for the HI perspective with a merger time scale of 0.4 Gyr. Of the three selection criteria defined in our previous paper, the GM performs well and cleanly selects mergers for 0.69 Gyr. The other two criteria (A-M20 and C-M20), select isolated disks as well, but perform best for face-on, gas-rich disks (T(merger) ~ 1 Gyr). The different visibility scales can be combined with the selected fractions of galaxies in any large HI survey to obtain merger rates in the nearby Universe. All-sky surveys such as WALLABY with ASKAP and the Medium Deep Survey with the APETIF instrument on Westerbork are set to revolutionize our perspective on neutral hydrogen and will provide an accurate measure of the merger fraction and rate of the present epoch., Comment: 12 pages, 6 figures, 4 tables, accepted by MNRAS, appendix not included

We perform 3-D dust radiative transfer (RT) calculations on hydrodynamic simulations of isolated and merging disk galaxies in order to quantitatively study the dependence of observed-frame submillimeter (submm) flux density on galaxy properties. We find that submm flux density and star formation rate (SFR) are related in dramatically different ways for quiescently star-forming galaxies and starbursts. Because the stars formed in the merger-induced starburst do not dominate the bolometric luminosity and the rapid drop in dust mass and more compact geometry cause a sharp increase in dust temperature during the burst, starbursts are very inefficient at boosting submm flux density (e.g., a $\ga16$x boost in SFR yields a $\la 2$x boost in submm flux density). Moreover, the ratio of submm flux density to SFR differs significantly between the two modes; thus one cannot assume that the galaxies with highest submm flux density are necessarily those with the highest bolometric luminosity or SFR. These results have important consequences for the bright submillimeter-selected galaxy (SMG) population. Among them are: 1. The SMG population is heterogeneous. In addition to merger-driven starbursts, there is a subpopulation of galaxy pairs, where two disks undergoing a major merger but not yet strongly interacting are blended into one submm source because of the large ($\ga 15$", or $\sim 130$ kpc at $z = 2$) beam of single-dish submm telescopes. 2. SMGs must be very massive ($M_{\star} \ga 6 \times 10^{10} \msun$). 3. The infall phase makes the SMG duty cycle a factor of a few greater than what is expected for a merger-driven starburst. (Abridged.), Comment: 15 pages, 3 figures. Accepted for publication in ApJ. Minor changes from v1

Theoretical models have had difficulty matching the observed number density of sub-millimeter galaxies (SMGs), causing some authors (e.g., Baugh et al. 2005) to suggest that SMGs provide evidence for a top-heavy initial mass function (IMF). To test this claim, we have, for the first time, combined high-resolution 3-D hydrodynamic simulations of isolated and merging massive, gas-rich galaxies, radiative transfer, and a semi-empirical merger rate model to predict the number density of SMGs. Our model can reproduce the observed SMG number density even when using a standard (Kroupa) IMF. The agreement is due to a combination of relatively long sub-mm duty cycles for mergers (a few times 10^8 years for our most massive models), which owe to our combination of high-resolution 3-D hydrodynamic simulations and dust radiative transfer; sufficient number densities of massive, gas-rich mergers; and the decrease in sub-mm counts observed by recent deep/wide surveys (e.g., Austermann et al. 2010) relative to previous surveys. Our results suggest that the observed SMG number counts do not provide evidence for a top-heavy IMF at high redshift., Comment: Minor typo in the abstract corrected

Recent high precision proper motions from the Hubble Space Telescope (HST) suggest that the Large and Small Magellanic Clouds (LMC and SMC, respectively) are either on their first passage or on an eccentric long period (>6 Gyr) orbit about the Milky Way (MW). This differs markedly from the canonical picture in which the Clouds travel on a quasi-periodic orbit about the MW (period of ~2 Gyr). Without a short period orbit about the MW, the origin of the Magellanic Stream, a young (1-2 Gyr old) coherent stream of HI gas that trails the Clouds ~150 degrees across the sky, can no longer be attributed to stripping by MW tides and/or ram pressure stripping by MW halo gas. We propose an alternative formation mechanism in which material is removed by LMC tides acting on the SMC before the system is accreted by the MW. We demonstrate the feasibility and generality of this scenario using an N-body/SPH simulation with cosmologically motivated initial conditions constrained by the observations. Under these conditions we demonstrate that it is possible to explain the origin of the Magellanic Stream in a first infall scenario. This picture is generically applicable to any gas-rich dwarf galaxy pair infalling towards a massive host or interacting in isolation., Comment: 9 pages, 2 figures,1 table, submitted to apjl

The outer regions of galactic disks have received increased attention since ultraviolet observations with GALEX demonstrated that nearly 30% of galaxies have UV emission beyond their optical extents, indicating star formation activity. These galaxies have been termed extended UV (XUV) disks. Here, we address whether these observations contradict the gas surface density threshold for star formation inferred from Halpha radial profiles of galaxies. We run smoothed particle hydrodynamics simulations of isolated disk galaxies with fiducial star formation prescriptions and show that over-densities owing to the presence of spiral structure can induce star formation in extended gas disks. For direct comparison with observations, we use the 3-D radiative transfer code Sunrise to create simulated FUV and K_s band images. We find that galaxies classified as Type I XUV disks are a natural consequence of spiral patterns, but we are unable to reproduce Type II XUV disks. We also compare our results to studies of the Kennicutt-Schmidt relation in outer disks., Comment: Published in ApJ

Gas-rich galaxy mergers are more easily identified by their disturbed morphologies than mergers with less gas. Because the typical gas fraction of galaxy mergers is expected to increase with redshift, the under-counting of low gas-fraction mergers may bias morphological estimates of the evolution of galaxy merger rate. To understand the magnitude of this bias, we explore the effect of gas fraction on the morphologies of a series of simulated disc galaxy mergers. With the resulting g-band images, we determine how the time-scale for identifying major and minor galaxy mergers via close projected pairs and quantitative morphology (the Gini coefficient G, the second-order moment of the brightest 20% of the light M20, and asymmetry A) depends on baryonic gas fraction f(gas). Strong asymmetries last significantly longer in high gas-fraction mergers of all mass ratios, with time-scales ranging from >= 300 Myr for f(gas) ~ 20% to >= 1 Gyr for f(gas) ~ 50%. Therefore the strong evolution with redshift observed in the fraction of asymmetric galaxies may reflect evolution in the gas properties of galaxies rather than the global galaxy merger rate. On the other hand, the time-scale for identifying a galaxy merger via G-M20 is weakly dependent on gas-fraction (~ 200-400 Myr), consistent with the weak evolution observed for G-M20 mergers., Comment: 15 pages; resubmitted to MNRAS

The majority of galaxy mergers are expected to be minor mergers. The observational signatures of minor mergers are not well understood, thus there exist few constraints on the minor merger rate. This paper seeks to address this gap in our understanding by determining if and when minor mergers exhibit disturbed morphologies and how they differ from the morphology of major mergers. We simulate a series of unequal-mass moderate gas-fraction disc galaxy mergers. With the resulting g-band images, we determine how the time-scale for identifying galaxy mergers via projected separation and quantitative morphology (the Gini coefficient G, asymmetry A, and the second-order moment of the brightest 20% of the light M20) depends on the merger mass ratio, relative orientations and orbital parameters. We find that G-M20 is as sensitive to 9:1 baryonic mass ratio mergers as 1:1 mergers, with observability time-scales ~ 0.2-0.4 Gyr. In contrast, asymmetry finds mergers with baryonic mass ratios between 4:1 and 1:1 (assuming local disc galaxy gas-fractions). Asymmetry time-scales for moderate gas-fraction major disc mergers are ~ 0.2-0.4 Gyr, and less than 0.06 Gyr for moderate gas-fraction minor mergers. The relative orientations and orbits have little effect on the time-scales for morphological disturbances. Observational studies of close pairs often select major mergers by choosing paired galaxies with similar luminosities and/or stellar masses. Therefore, the various ways of finding galaxy mergers (G-M20, A, close pairs) are sensitive to galaxy mergers of different mass ratios. By comparing the frequency of mergers selected by different techniques, one may place empirical constraints on the major and minor galaxy merger rates., Comment: 16 pages; resubmitted to MNRAS

An updated version of the dust radiation transfer code Sunrise, including models for star-forming regions and a self-consistent calculation of the spatially dependent dust and PAH emission, is presented. Given a hydrodynamic simulation of a galaxy, this model can calculate a realistic 2-dimensional ultraviolet--submillimeter spectral energy distribution of the galaxy, including emission lines from HII regions, from any viewpoint. To model the emission from star-forming regions, the MAPPINGSIII photoionization code is used. The high wavelength resolution (~ 1000 wavelengths) is made possible by the polychromatic Monte-Carlo algorithm employed by Sunrise. From the 2-D spectral energy distributions, images in any filter bands or integrated galaxy SEDs can be created. Using a suite of hydrodynamic simulations of disc galaxies, the output broad-band images and spectral energy distributions are compared with observed galaxies from the multiwavelength SINGS and SLUGS galaxy surveys. Overall, the output spectral energy distributions show a good match with observed galaxies in colours ranging from GALEX far-UV to SCUBA submillimeter wavelengths. The only possible exception is the 160 micron/850 micron colour, which the simulations underestimate by a factor "of order 5" compared to the SINGS sample. However, the simulations here agree with the SLUGS galaxies, which consistently have significantly larger amounts of cold dust than the SINGS galaxies. The Sunrise model can be used to generate simulated observations of arbitrary hydrodynamic galaxy simulations. In this way, predictions of galaxy formation theories can be directly tested against observations of galaxies., Comment: 32 pages, 22 figures, submitted to MNRAS

We present a merger-driven evolutionary model for the production of luminous (LIRGs) and ultraluminous infrared galaxies (ULIRGs) with warm IR colours. Our results show that simulations of gas--rich major mergers including star formation, black hole growth, and feedback can produce warm (U)LIRGs. We also find that while the warm evolutionary phase is associated with increased AGN activity, star formation alone may be sufficient to produce warm IR colours. However, the transition can be suppressed entirely - even when there is a significant AGN contribution - when we assume a single-phase ISM, which maximizes the attenuation. Finally, our evolutionary models are consistent with the 25-to-60 micron flux density ratio vs. L_HX/L_IR relation for local LIRGs and ULIRGs, and predict the observed scatter in IR color at fixed L_HX/L_IR. Therefore, our models suggest a cautionary note in the interpretation of warm IR colours: while associated with periods of active black hole growth, they are probably produced by a complex mix of star formation and AGN activity intermediate between the cold star formation dominated phase and the birth of a bright, unobscured quasar., Comment: 5 pages, accepted for publication in MNRAS Letters

There is a large observational scatter toward low velocities in the stellar mass Tully-Fisher relation if disturbed and compact objects are included. However, this scatter can be eliminated if one replaces rotation velocity with $\rm S_{\rm 0.5}$, a quantity that includes a velocity dispersion term added in quadrature with the rotation velocity. In this work we use a large suite of hydrodynamic N-body galaxy merger simulations to explore a possible mechanism for creating the observed relations. Using mock observations of the simulations, we test for the presence of observational effects and explore the relationship between $\rm S_{\rm 0.5}$ and intrinsic properties of the galaxies. We find that galaxy mergers can explain the scatter in the TF as well as the tight $\rm S_{\rm 0.5}$-stellar mass relation. Furthermore, $\rm S_{\rm 0.5}$ is correlated with the total central mass of a galaxy, including contributions due to dark matter., Comment: ApJ accepted

We present a detailed comparison of the CO(3-2) emitting molecular gas between a local sample of luminous infrared galaxies (U/LIRGs) and a high redshift sample that comprises submm selected galaxies (SMGs), quasars, and Lyman Break Galaxies (LBGs). The U/LIRG sample consists of our recent CO(3-2) survey using the Submillimeter Array while the CO(3-2) data for the high redshift population are obtained from the literature. We find that the L(CO(3-2)) and L(FIR) relation is correlated over five orders of magnitude, which suggests that the molecular gas traced in CO(3-2) emission is a robust tracer of dusty star formation activity. The near unity slope of 0.93 +/- 0.03 obtained from a fit to this relation suggests that the star formation efficiency is constant to within a factor of two across different types of galaxies residing in vastly different epochs. The CO(3-2) size measurements suggest that the molecular gas disks in local U/LIRGs (0.3 - 3.1 kpc) are much more compact than the SMGs (3 - 16 kpc), and that the size scales of SMGs are comparable to the nuclear separation (5 - 40 kpc) of the widely separated nuclei of U/LIRGs in our sample. We argue from these results that the SMGs studied here are predominantly intermediate stage mergers, and that the wider line-widths arise from the violent merger of two massive gas-rich galaxies taking place deep in a massive halo potential., Comment: 16 pages, 5 figures, ApJ Accepted

The HST proper motion (PM) measurements of the Clouds have severe implications for their interaction history with the Milky Way (MW) and with each other. The Clouds are likely on their first passage about the MW and the SMC's orbit about the LMC is better described as quasi-periodic rather than circular. Binary L/SMC orbits that satisfy observational constraints on their mutual interaction history (e.g. the formation of the Magellanic Bridge during a collision between the Clouds ~300 Myr ago) can be located within 1 sigma of the mean PMs. However, these binary orbits are not co-located with the Magellanic Stream (MS) when projected on the plane of the sky and the line-of-sight velocity gradient along the LMC's orbit is significantly steeper than that along the MS. These combined results ultimately rule out a purely tidal origin for the MS: tides are ineffective without multiple pericentric passages and can neither decrease the velocity gradient nor explain the offset stream in a polar orbit configuration. Alternatively, ram pressure stripping of an extended gaseous disk may naturally explain the deviation. The offset also suggests that observations of the little-explored region between RA 21h and 23h are crucial for characterizing the full extent of the MS., Comment: 6 pages, 5 figures, to be published in proceedings of IAU Symposium No. 256: The Magellanic System: Stars, Gas, and Galaxies, Jacco Th. van Loon & Joana M. Oliveira, eds. 2009 (Cambridge: Cambridge University Press)

The outer regions of disk galaxies show a drop-off in optical and Halpha emission, suggesting a star formation threshold radius, assumed to owe to a critical surface density below which star formation does not take place. Signs of filamentary star formation beyond this threshold radius have been observed in individual galaxies in the Halpha and recent GALEX surveys have discovered that 30% of disk galaxies show UV emission out to 2-3 times the optical radius of the galaxy. We run smooth particle hydrodynamics simulations of disk galaxies with constant density extended gas disks to test whether over-densities owing to spiral structure in the outer disk can reproduce the observed star formation. We indeed find that spiral density waves from the inner disk propagate into the outer gas disk and raise local gas regions above the star formation density threshold, yielding features similar to those observed. Because the amount of star formation is low, we expect to see little optical emission in outer disks, as observed. Our results indicate that XUV disks can be simulated simply by adding an extended gas disk with a surface density near the threshold density to an isolated galaxy and evolving it with fiducial star formation parameters., Comment: Accepted for publication in ApJL

We present new data obtained with the Submillimeter Array for a sample of fourteen nearby luminous and ultraluminous infrared galaxies. The galaxies were selected to have luminosity distances D < 200 Mpc and far-infrared luminosities log(L_FIR) > 11.4. The galaxies were observed with spatial resolutions of order 1 kpc in the CO J=3-2, CO J=2-1, 13CO J=2-1, and HCO+ J=4-3 lines as well as the continuum at 880 microns and 1.3 mm. We have combined our CO and continuum data to measure an average gas-to-dust mass ratio of 120 +/- 28 (rms deviation 109) in the central regions of these galaxies, very similar to the value of 150 determined for the Milky Way. This similarity is interesting given the more intense heating from the starburst and possibly accretion activity in the luminous infrared galaxies compared to the Milky Way. We find that the peak H_2 surface density correlates with the far-infrared luminosity, which suggests that galaxies with higher gas surface densities inside the central kiloparsec have a higher star formation rate. The lack of a significant correlation between total H_2 mass and far-infrared luminosity in our sample suggests that the increased star formation rate is due to the increased availability of molecular gas as fuel for star formation in the central regions. In contrast to previous analyses by other authors, we do not find a significant correlation between central gas surface density and the star formation efficiency, as trace by the ratio of far-infrared luminosity to nuclear gas mass. Our data show that it is the star formation rate, not the star formation efficiency, that increases with increasing central gas surface density in these galaxies., Comment: 66 pages, 39 figures, aastex preprint format; to be published in ApJ Supplements. Version of paper with full resolution figures available at http://www.physics.mcmaster.ca/~wilson/www_xfer/ULIRGS_public

A key obstacle to understanding the galaxy merger rate and its role in galaxy evolution is the difficulty in constraining the merger properties and time-scales from instantaneous snapshots of the real universe.The most common way to identify galaxy mergers is by morphology, yet current theoretical calculations of the time-scales for galaxy disturbances are quite crude. We present a morphological analysis of a large suite of GADGET N-Body/hydro-dynamical equal-mass gas-rich disc galaxy mergers which have been processed through the Monte-Carlo radiative transfer code SUNRISE. With the resulting images, we examine the dependence of quantitative morphology (G, M20, C, A) in the SDSS g-band on merger stage, dust, viewing angle, orbital parameters, gas properties, supernova feedback, and total mass. We find that mergers appear most disturbed in G-M20 and asymmetry at the first pass and at the final coalescence of their nuclei, but can have normal quantitative morphologies at other merger stages. The merger observability time-scales depend on the method used to identify the merger as well as the gas fraction, pericentric distance, and relative orientation of the merging galaxies. Enhanced star formation peaks after and lasts significantly longer than strong morphological disturbances. Despite their massive bulges, the majority of merger remnants appear disc-like and dusty in g-band light because of the presence of a low-mass star-forming disc., Comment: accepted to MNRAS; 29 pages, 19 color figures; go to http://www.noao.edu/staff/lotz/pdfs/lotz_arXiv.0805.1246.pdf for full resolution

We present a series of simulations of the self--regulated growth of supermassive black holes (SMBHs) in galaxies via three different fueling mechanisms: major mergers, minor mergers, and disk instabilities. The SMBHs in all three scenarios follow the same black hole fundamental plane (BHFP) and correlation with bulge binding energy seen in simulations of major mergers, and observed locally. Furthermore, provided that the total gas supply is significantly larger than the mass of the SMBH, its limiting mass is not influenced by the amount of gas available or the efficiency of black hole growth. This supports the assertion that SMBHs accrete until they reach a critical mass at which feedback is sufficient to unbind the gas locally, terminating the inflow and stalling further growth. At the same time, while minor and major mergers follow the same projected correlations (e.g., the $M_{BH}-\sigma$ and Magorrian relations), SMBHs grown via disk instabilities do not, owing to structural differences between the host bulges. This finding is supported by recent observations of SMBHs in pseudobulges and bulges in barred systems, as compared to those hosted by classical bulges. Taken together, this provides support for the BHFP and binding energy correlations as being more "fundamental" than other proposed correlations in that they reflect the physical mechanism driving the co-evolution of SMBHs and spheroids., Comment: 15 pages, 16 figures, accepted for publication in ApJ