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

We present a scenario for the production of dynamically cold rings around the Milky Way via a high-eccentricity, flyby encounter. These initial conditions are more cosmologically motivated than those considered in previous works. We find that the encounters we examine generically produce a series of nearly dynamically cold ring-like features on low-eccentricity orbits that persist over timescales of ~2-4 Gyr via the tidal response of the primary galaxy to the close passage of the satellite. Moreover, they are both qualitatively and quantitatively similar to the distribution, kinematics, and stellar population of the Monoceros ring. Therefore, we find that a high eccentricity flyby by a satellite galaxy represents a cosmologically appealing scenario for forming kinematically distinct ring-like features around the Milky Way., Comment: accepted to ApJL; 4 Figures

We construct a physically motivated model for predicting the properties of the remnants of gaseous galaxy mergers, given the properties of the progenitors and the orbit. The model is calibrated using a large suite of SPH merger simulations. It implements generalized energy conservation while accounting for dissipative energy losses and star formation. The dissipative effects are evaluated from the initial gas fractions and from the orbital parameters via an "impulse" parameter, which characterizes the strength of the encounter. Given the progenitor properties, the model predicts the remnant stellar mass, half-mass radius, and velocity dispersion to an accuracy of 25%. The model is valid for both major and minor mergers. We provide an explicit recipe for semi-analytic models of galaxy formation., Comment: 14 pages, 9 figures. MNRAS accepted

We assess the impact of starburst and AGN feedback-driven winds on the CO emission from galaxy mergers, and, in particular, search for signatures of these winds in the simulated CO morphologies and emission line profiles. We do so by combining a 3D non-LTE molecular line radiative transfer code with smoothed particle hydrodynamics (SPH) simulations of galaxy mergers that include prescriptions for star formation, black hole growth, a multiphase interstellar medium (ISM), and the winds associated with star formation and black hole growth. Our main results are: (1) Galactic winds can drive outflows of masses ~10^8-10^9 Msun which may be imaged via CO emission line mapping. (2) AGN feedback-driven winds are able to drive imageable CO outflows for longer periods of time than starburst-driven winds owing to the greater amount of energy imparted to the ISM by AGN feedback compared to star formation. (3) Galactic winds can control the spatial extent of the CO emission in post-merger galaxies, and may serve as a physical motivation for the sub-kiloparsec scale CO emission radii observed in local advanced mergers. (4) Secondary emission peaks at velocities greater than the circular velocity are seen in the CO emission lines in all models. In models with winds, these high velocity peaks are seen to preferentially correspond to outflowing gas entrained in winds, which is not the case in the model without winds. The high velocity peaks seen in models without winds are typically confined to velocity offsets (from the systemic) < 1.7 times the circular velocity, whereas the models with AGN feedback-driven winds can drive high velocity peaks to ~2.5 times the circular velocity., Comment: Accepted by ApJ; Minor revisions; Resolution tests included

We employ numerical simulations of galaxy mergers to explore the effect of galaxy mass ratio on merger--driven starbursts. Our numerical simulations include radiative cooling of gas, star formation, and stellar feedback to follow the interaction and merger of four disk galaxies. The galaxy models span a factor of 23 in total mass and are designed to be representative of typical galaxies in the local Universe. We find that the merger--driven star formation is a strong function of merger mass ratio, with very little, if any, induced star formation for large mass ratio mergers. We define a burst efficiency that is useful to characterize the merger--driven star formation and test that it is insensitive to uncertainties in the feedback parameterization. In accord with previous work we find that the burst efficiency depends on the structure of the primary galaxy. In particular, the presence of a massive stellar bulge stabilizes the disk and suppresses merger--driven star formation for large mass ratio mergers. Direct, co--planar merging orbits produce the largest tidal disturbance and yield that most intense burst of star formation. Contrary to naive expectations, a more compact distribution of gas or an increased gas fraction both decrease the burst efficiency. Owing to the efficient feedback model and the newer version of SPH employed here, the burst efficiencies of the mergers presented here are smaller than in previous studies., Comment: 26 pages, 21 figures, submitted to MNRAS

We use hydrodynamic simulations of minor mergers of galaxies to investigate the nature of surface brightness excesses at large radii observed in some spiral galaxies: antitruncated stellar disks. We find that this process can produce the antitruncation via two competing effects: (1) merger-driven gas inflows that concentrate mass in the center of the primary galaxy and contract its inner density profile; and (2) angular momentum transferred outwards by the interaction, causing the outer disk to expand. In our experiments, this requires both a significant supply of gas in the primary disk, and that the encounter be prograde with moderate orbital angular momentum. The stellar surface mass density profiles of our remnants both qualitatively and quantitatively resemble the broken exponentials observed in local face--on spirals that display antitruncations. Moreover, the observed trend towards more frequent antitruncation relative to classical truncation in earlier Hubble types is consistent with a merger-driven scenario., Comment: 10 pages, 11 figures, accepted to ApJ

We use a N--body/hydrodynamic simulation to forecast the future encounter between the Milky Way and the Andromeda galaxies, given current observational constraints on their relative distance, relative velocity, and masses. Allowing for a comparable amount of diffuse mass to fill the volume of the Local Group, we find that the two galaxies are likely to collide in a few billion years - within the Sun's lifetime. During the the interaction, there is a chance that the Sun will be pulled away from its present orbital radius and reside in an extended tidal tail. The likelihood for this outcome increases as the merger progresses, and there is a remote possibility that our Sun will be more tightly bound to Andromeda than to the Milky Way before the final merger. Eventually, after the merger has completed, the Sun is most likely to be scattered to the outer halo and reside at much larger radii (>30 kpc). The density profiles of the stars, gas and dark matter in the merger product resemble those of elliptical galaxies. Our Local Group model therefore provides a prototype progenitor of late--forming elliptical galaxies., Comment: accepted to MNRAS

Recent proper motion measurements of the Large and Small Magellanic Clouds (LMC and SMC, respectively) by Kallivayalil et al (2006a,b) suggest that the 3D velocities of the Clouds are substantially higher (~100 km/s) than previously estimated and now approach the escape velocity of the Milky Way (MW). Previous studies have also assumed that the Milky Way can be adequately modeled as an isothermal sphere to large distances. Here we re-examine the orbital history of the Clouds using the new velocities and a LCDM-motivated MW model with virial mass Mvir = 1e12 Msun (e.g. Klypin et al 2002). We conclude that the Clouds are either currently on their first passage about the MW or, if the MW can be accurately modeled by an isothermal sphere to distances >200 kpc (i.e., Mvir > 2e12 Msun), that their orbital period and apogalacticon distance must be a factor of two larger than previously estimated, increasing to 3 Gyr and 200 kpc, respectively. A first passage scenario is consistent with the fact that the LMC and SMC appear to be outliers when compared to other satellite galaxies of the MW: they are irregular in appearance and are moving faster. We discuss the implications of this orbital analysis for our understanding of the star formation history, the nature of the warp in the MW disk and the origin of the Magellanic Stream (MS), a band of HI gas trailing the LMC and SMC that extends roughly 100 degrees across the sky. Specifically, as a consequence of the new orbital history of the Clouds, the origin of the MS may not be explainable by current tidal and ram pressure stripping models.

We study the effects of dust in hydrodynamic simulations of spiral galaxies when different radial metallicity gradients are assumed. SUNRISE, a Monte-Carlo radiative-transfer code, is used to make detailed calculations of the internal extinction of disk galaxies caused by their dust content. SUNRISE is used on eight different Smooth Particle Hydrodynamics (SPH) simulations of isolated spiral galaxies. These galaxies vary mainly in mass and hence luminosity, spanning a range in luminosities from -16 to -22 magnitudes in the B band. We focus on the attenuation in different wavelength bands as a function of the disk inclination and the luminosity of the models, and compare this to observations. Observations suggest different metallicity gradients for galaxies of different luminosities. These metallicity gradients were explored in our different models, finding that the resulting dust attenuation matches observations for edge-on galaxies, but do not show a linear behaviour in log axis ratio as some observations have suggested. A quadratic law describing the dependence of attenuation on inclination, as proposed by more recent observations, reconciles the attenuation of the simulations at intermediate inclinations with observations. We also compare the total infrared-to-ultraviolet flux ratios for the simulated galaxies with those of the SINGS sample and find general agreement. Finally we compare our results with those from simpler models that do not take into account structure such as spiral arms, finding that the inclusion of sub-structure on the size scale of spiral arms does not change conclusions about the attenuation dependence on inclination or wavelength., Comment: 12 pages, 8 figures, all figures require colors, accepted for publication in MNRAS. Major changes from v1 include: Comparison of the calculated infrared-to-ultraviolet ratios in our models to those observed by SIGNS, clearer fonts in our plots and clearer justification of the rationale behind the adopted gradients. Changes from v2: Added journal-ref

[abridged]We calculate multi-wavelength spectral energy distributions (SEDs) (spanning optical to millimeter wavelengths) from simulations of major galaxy mergers with black hole feedback which produce submillimeter bright galaxies (SMGs), using a self-consistent three-dimensional radiative transfer code. We reproduce correlations for local AGN observed in Spitzer Space Telescope's IRAC bands, and make definitive predictions for infrared X-ray correlations that should be testable by combining observations by Spitzer and the upcoming Herschel mission with X-ray surveys. Our dynamical approach allows us to directly correlate observed clustering in the data as seen in IRAC color-color plots with the relative amount of time the system spends in a region of color-color space. We also find that this clustering is positively correlated with the stars dominating in their contribution to the total bolometric luminosity. We present photo albums spanning the lifetime of SMGs, from their infancy in the pre-merger phase to the final stage as an elliptical galaxy, as seen in the observed $3.6 \micron$ and $450 \micron$ band to visually illustrate some of the morphological differences between mergers of differing orbital inclination and progenitor redshift. We find that SMGs are a broader class of systems than starbursts or quasars. We introduce a simple, heuristic classification scheme on the basis of the $L_{\rm IR}/L_{\rm x}$ ratios of these galaxies, which may be interpreted qualitatively as an evolutionary scheme, as these galaxies evolve in $L_{\rm IR}/L_{\rm x}$ while transiting from the pre-merger stage, through the quasar phase, to a merger remnant., Comment: 19 pages, 27 figures, Submitted to ApJ, higher resolution figures can be found at http://www.cfa.harvard.edu/~schakrab/smg.ps

We use numerical simulations to study the kinematic structure of remnants formed from mergers of equal-mass disk galaxies. In particular, we show that remnants of dissipational mergers, which include the radiative cooling of gas, star formation, feedback from supernovae, and the growth of supermassive black holes, are smaller, rounder, have, on average, a larger central velocity dispersion, and show significant rotation compared to remnants of dissipationless mergers. The increased rotation speed of dissipational remnants owes its origin to star formation that occurs in the central regions during the galaxy merger. We have further quantified the anisotropy, three-dimensional shape, minor axis rotation, and isophotal shape of each merger remnant, finding that dissipational remnants are more isotropic, closer to oblate, have the majority of their rotation along their major axis, and are more disky than dissipationless remnants. Individual remnants display a wide variety of kinematic properties. A large fraction of the dissipational remnants are oblate isotropic rotators. Many dissipational, and all of the dissipationless, are slowly rotating and anisotropic. The remnants of gas-rich major mergers can well-reproduce the observed distribution of projected ellipticities, rotation parameter (V/\sigma)*, kinematic misalignments, Psi, and isophotal shapes. The dissipationless remnants are a poor match to this data. Our results support the merger hypothesis for the origin of low-luminosity elliptical galaxies provided that the progenitor disks are sufficiently gas-rich, however our remnants are a poor match to the bright ellipticals that are slowly rotating and uniformly boxy., Comment: 22 pages, 17 figures, accepted to ApJ

We calculate infrared spectral energy distributions (SEDs) from simulations of major galaxy mergers and study the effect of AGN and starburst driven feedback on the evolution of the SED as a function of time. We use a self-consistent three-dimensional radiative equilibrium code to calculate the emergent SEDs and to make images. To facilitate a simple description of our findings, we describe our results in reference to an approximate analytic solution for the far-IR SED. We focus mainly on the luminous infrared galaxy (LIRG) and ultraluminous infrared galaxy (ULIRG) phases of evolution. We contrast the SEDs of simulations performed with AGN feedback to simulations performed with starburst driven wind feedback. We find that the feedback processes critically determine the evolution of the SED. Changing the source of illumination (whether stellar or AGN) has virtually no impact on the reprocessed far-infrared SED. We find that AGN feedback is particularly effective at dispersing gas and rapidly injecting energy into the ISM. The observational signature of such powerful feedback is a warm SED. In general, simulations performed with starburst driven winds have colder spectra and reprocess more of their emission into the infrared, resulting in higher infrared to bolometric luminosities compared to (otherwise equivalent) simulations performed with AGN feedback. We depict our results in IRAS bands, as well as in Spitzer's MIPS bands, and in Herschel's PACS bands., Comment: Accepted to ApJ, 11 figures, 26 pages, added 2 figures

Elements of kinematical and dynamical modeling of elliptical galaxies (Es) are presented. In projection, NFW models resemble Sersic models, but with a very narrow range of shapes (m=3+/-1). The total density profile of Es cannot be NFW-like because the predicted local M/L and aperture velocity dispersion within an effective radius (Re) are much lower than observed. Stars must then dominate Es out to a few Re. Fitting an NFW model to the total density profile of Sersic+NFW (stars+dark matter [DM]) Es results in very high concentration parameters, as found by X-ray observers. Kinematical modeling of Es assuming an isotropic NFW DM model underestimates M/L at the virial radius by a factor of 1.6 to 2.4, because dissipationless LCDM halos have slightly different density profiles and slightly radial velocity anisotropy. In N-body+gas simulations of Es as merger remnants of spirals embedded in DM halos, the slope of the DM density profile is steeper when the initial spiral galaxies are gas-rich. The Hansen & Moore (2006) relation between anisotropy and the slope of the density profile breaks down for gas and DM, but the stars follow an analogous relation with slightly less radial anisotropies for a given density slope. Using kurtosis (h_4) to infer anisotropy in Es is dangerous, as h_4 is also sensitive to small levels of rotation. The stationary Jeans equation provides accurate masses out to 8 Re. The discrepancy between the modeling of Romanowsky et al. (2003), indicating a dearth of DM in Es, and the simulations analyzed by Dekel et al. (2005), which match the spectroscopic observations of Es, is partly due to radial anisotropy and to observing oblate Es face-on. However, one of the 15 solutions to the orbit modeling of Romanowsky et al. is found to have an amount and concentration of DM consistent with LCDM predictions., Comment: Invited talk, 21st IAP meeting, Mass Profiles and Shapes of Cosmological Structures, ed. G. A. Mamon, F. Combes, C. Deffayet & B. Fort (Paris: EDP). 10 pages, 7 figures (8 plots)

(ABRIDGED) We examine the fundamental scaling relations of elliptical galaxies formed through mergers. Using hundreds of simulations to judge the impact of progenitor galaxy properties on merger remnants, we find that gas dissipation provides an important contribution to tilt in the Fundamental Plane relation. Dissipationless mergers of disks produce remnants that occupy the virial plane. As the gas content of disk galaxies is increased, the tilt of the Fundamental Plane relation increases and the slope of the Re-M_* relation steepens. For gas fractions fgas > 30%, the simulated Fundamental Plane scalings approach those observed in the K-band. In our simulations, feedback from supermassive black hole growth has only a minor influence on the stellar-mass scaling relations of spheroidal galaxies, but may play a role in maintaining the observed Fundamental Plane tilt at optical wavelengths by suppressing residual star formation in merger remnants. We estimate that \approx 40-100% of the Fundamental Plane tilt induced by structural properties owes to trends in the central total-to-stellar mass ratio M_total/M_* produced by dissipation. Lower mass systems obtain greater phase- space densities than higher mass systems, producing a galaxy mass-dependent central M_total/M_* and a corresponding tilt in the Fundamental Plane., Comment: Version accepted by ApJ, 20 pages, 18 figures, resolution reduced for size

We examine X-ray emission produced from hot gas during collisions and mergers of disk galaxies. To study this process, we employ simulations that incorporate cosmologically motivated disk-galaxy models and include the effects of radiative cooling, star formation, supernova feedback, and accreting supermassive black holes. We find that during a merger, the colliding gas in the disks is shock-heated to X-ray-emitting temperatures. The X-ray luminosity is spatially extended, rises during the initial stages of the merger, and peaks when the galactic centers coalesce. When a physical model for accreting black holes is included, the resulting feedback can drive powerful winds that contribute significantly to the amount and metallicity of hot gas, both of which increase the X-ray luminosity. In terms of their stellar kinematics and structural properties, the merger remnants in our simulations resemble elliptical galaxies. We find that the X-ray luminosities of the remnants with B-band luminosities in the range L_B ~ 10^10 - 10^11 Lsun are consistent with observations, while remnants with smaller or larger masses are underluminous in X-rays. Moreover, because the majority of the merger remnants are broadly consistent with the observed scaling relations between temperature, B-band luminosity and X-ray luminosity we conclude that major mergers are a viable mechanism for producing the X-ray halos of large, luminous elliptical galaxies., Comment: 16 pages, 15 figures, matches version accepted to ApJ

Using hydrodynamic simulations of disc-galaxy major mergers, we investigate the star formation history and remnant properties when various parametrizations of a simple stellar feedback model are implemented. The simulations include radiative cooling, a density-dependent star formation recipe and a model for feedback from massive stars. The feedback model stores supernova feedback energy within individual gas particles and dissipates this energy on a time-scale specified by two free parameters; tau_fb, which sets the dissipative time-scale, and n, which sets the effective equation of state in star-forming regions. Using a self-consistent disc galaxy, modelled after a local Sbc spiral, in both isolated and major-merger simulations, we investigate parametrizations of the feedback model that are selected with respect to the quiescent disc stability. These models produce a range of star formation histories that are consistent with the star formation relation found by Kennicutt. All major mergers produce a population of new stars that is highly centrally concentrated, demonstrating a distinct break in the r1/4 surface density profile, consistent with previous findings. The half-mass radius and one-dimensional velocity dispersion are affected by the feedback model used. Finally, we compare our results to those of previous simulations of star formation in disc-galaxy major mergers, addressing the effects of star formation normalization, the version of smoothed particle hydrodynamics (SPH) employed and assumptions about the interstellar medium., Comment: update to match published version

A new Monte-Carlo radiative-transfer code, Sunrise, is used in conjunction with hydrodynamic simulations of major galaxy mergers to calculate the effects of dust in such systems. The simulations are in good agreement with observations of dust absorption in starburst galaxies, and the dust has a profound effect on their appearance. The dust attenuation increases with luminosity such that at peak luminosities ~90% of the bolometric luminosity is absorbed by dust. In general, the detailed appearance of the merging event depends on the stage of the merger and the geometry of the encounter. The fraction of bolometric energy absorbed by the dust, however, is a robust quantity that can be predicted from the intrinsic properties bolometric luminosity, baryonic mass, star-formation rate, and metallicity of the system. This paper presents fitting formulae, valid over a wide range of masses and metallicities, from which the absorbed fraction of luminosity (and consequently also the infrared dust luminosity) can be predicted. The attenuation of the luminosity at specific wavelengths can also be predicted, albeit with a larger scatter due to the variation with viewing angle. These formulae for dust attenuation appear to be valid for both isolated and interacting galaxies, are consistent with earlier studies, and would be suitable for inclusion in theoretical models, e.g. semi-analytic models of galaxy formation., Comment: 12 pages, 10 figures, submitted to ApJ

There is strong evidence that the mass in the Universe is dominated by dark matter, which exerts gravitational attraction but whose exact nature is unknown. In particular, all galaxies are believed to be embedded in massive haloes of dark matter. This view has recently been challenged by surprisingly low random stellar velocities in the outskirts of ordinary elliptical galaxies, which were interpreted as indicating a lack of dark matter (Mendez et al. 2001; Romanowsky et al. 2003). Here we show that the low velocities are in fact compatible with galaxy formation in dark-matter haloes. Using numerical simulations of disc-galaxy mergers, we find that the stellar orbits in the outer regions of the resulting ellipticals are very elongated. These stars were torn by tidal forces from their original galaxies during the first close passage and put on outgoing trajectories. The elongated orbits, combined with the steeply falling density profile of the observed tracers, explain the observed low velocities even in the presence of large amounts of dark matter. Projection effects when viewing a triaxial elliptical can lead to even lower observed velocities along certain lines of sight., Comment: Letter to Nature, 13+15 pages, 2+11 figures, improved text, extended Supplementary Information added

A new Monte-Carlo radiative-transfer code, Sunrise, is used to study the effects of dust in N-body/hydrodynamic simulations of interacting galaxies. Dust has a profound effect on the appearance of the simulated galaxies. At peak luminosities, about 90% of the bolometric luminosity is absorbed, and the dust obscuration scales with luminosity in such a way that the brightness at UV/visual wavelengths remains roughly constant. A general relationship between the fraction of energy absorbed and the ratio of bolometric luminosity to baryonic mass is found. Comparing to observations, the simulations are found to follow a relation similar to the observed IRX-Beta relation found by Meurer et al (1999) when similar luminosity objects are considered. The highest-luminosity simulated galaxies depart from this relation and occupy the region where local (U)LIRGs are found. This agreement is contingent on the presence of Milky-Way-like dust, while SMC-like dust results in far too red a UV continuum slope to match observations. The simulations are used to study the performance of star-formation indicators in the presence of dust. The far-infrared luminosity is found to be reliable. In contrast, the H-alpha and far-UV luminosity suffer severely from dust attenuation, and dust corrections can only partially remedy the situation., Comment: 4 pages, 5 figures, to appear in the proceedings of the conference "The Spectral Energy Distribution of Gas-Rich Galaxies", eds. C.C. Popescu & R.J. Tuffs (Heidelberg, October 2004)

We report on the merger-induced generation of a shock-heated gas wind and formation of a remnant gas halo in simulations of colliding disk galaxies. The simulations use cosmologically motivated initial conditions and include the effects of radiative cooling, star formation, stellar feedback and the non-adiabatic heating of gas. The non-adiabatic heating, i.e. shocks, generated in the final merger forces gas out of the central region of the merger remnant and into the dark-matter halo. We demonstrate that the amount of heating depends on the size of the progenitor disk galaxy as well as the initial orbit the galaxies are placed on. Based upon these dependencies, we motivate a possible recipe for including this effect in semi-analytic models of galaxy formation., Comment: 4 pages, 2 figures, uses emulateapj.cls, submitted to ApJ Letters

Chemically active climate compounds are either primary compounds such as methane (CH4), removed by oxidation in the atmosphere, or secondary compounds such as ozone (O3), sulfate and organic aerosols, formed and removed in the atmosphere. Man-induced climate-chemistry interaction is a two-way process: Emissions of pollutants change the atmospheric composition contributing to climate change through the aforementioned climate components, and climate change, through changes in temperature, dynamics, the hydrological cycle, atmospheric stability, and biosphere-atmosphere interactions, affects the atmospheric composition and oxidation processes in the troposphere. Here we present progress in our understanding of processes of importance for climate-chemistry interactions, and their contributions to changes in atmospheric composition and climate forcing. A key factor is the oxidation potential involving compounds such as O3 and the hydroxyl radical (OH). Reported studies represent both current and future changes. Reported results include new estimates of radiative forcing based on extensive model studies of chemically active climate compounds such as O3, and of particles inducing both direct and indirect effects. Through EU projects such as ACCENT, QUANTIFY, and the AEROCOM project, extensive studies on regional and sector-wise differences in the impact on atmospheric distribution are performed. Studies have shown that land-based emissions have a different effect on climate than ship and aircraft emissions, and different measures are needed to reduce the climate impact. Several areas where climate change can affect the tropospheric oxidation process and the chemical composition are identified. This can take place through enhanced stratospheric-tropospheric exchange of ozone, more frequent periods with stable conditions favouring pollution build up over industrial areas, enhanced temperature-induced biogenic emissions, methane releases from permafrost thawing, and enhanced concentration through reduced biospheric uptake. During the last 510 years, new observational data have been made available and used for model validation and the study of atmospheric processes. Although there are significant uncertainties in the modelling of composition changes, access to new observational data has improved modelling capability. Emission scenarios for the coming decades have a large uncertainty range, in particular with respect to regional trends, leading to a significant uncertainty range in estimated regional composition changes and climate impact.

The NASA Dryden Flight Research Center and the Technical University of Munich are cooperating in a research program to assess the impact of unstable long-term dynamics on the flying qualities of planes in hypersonic flight. These flying qualities issues are being investigated with a dedicated flight simulator for hypersonic vehicles located at NASA Dryden. Several NASA research pilots have flown the simulator through well-defined steady-level turns with varying phugoid and height mode instabilities. Th data collected include pilot ratings and comments, performance measurements, and pilot workload measurements. The results presented in this paper include design guidelines for height and phugoid mode instabilities, an evaluation of the tapping method used to measure pilot workload, a discussion of techniques developed by the pilots to control large instabilities, and a discussion of how flying qualities of unstable long-term dynamics influence control power design requirements.

Over a 4-year period, compared scholastic attainment of children from disadvantaged homes with that of children from working class homes. Results supported the hypothesis that scholastic learning in culturally disadvantaged children shows a cumulative deficit over a time, although not all differences between the groups were statistically significant. (GC)