Search

The origin of magnetic fields in white dwarfs (WDs) remains mysterious. Magnetic WDs are traditionally associated with field strengths $\gtrsim1\,\mathrm{MG}$, set by the sensitivity of typical spectroscopic magnetic field measurements. Informed by recent developments in red giant magnetoasteroseismology, we revisit the use of WD pulsations as a seismic magnetometer. WD pulsations primarily probe near-surface magnetic fields, whose effect on oscillation mode frequencies is to asymmetrize rotational multiplets and, if strong enough, suppress gravity-mode propagation altogether. The sensitivity of seismology to magnetic fields increases strongly with mode period and decreases quickly with the depth of the partial ionization-driven surface convective zone. We place upper limits for magnetic fields in $24$ pulsating WDs: $20$ hydrogen-atmosphere (DAV) and three helium-atmosphere (DBV) carbon-oxygen WDs, and one extremely low-mass (helium-core) pulsator. These bounds are typically $\sim1$-$10\,\mathrm{kG}$, although they can reach down to $\sim10$-$100\,\mathrm{G}$ for DAVs and helium-core WDs in which lower-frequency modes are excited. Seismic magnetometry may enable new insights into the formation and evolution of WD magnetism., Comment: 21 pages, 9 figures, submitted to ApJ

Mergers between helium white dwarfs and main-sequence stars are likely common, producing red giant-like remnants making up roughly a few percent of all low-mass ($\lesssim2M_\odot$) red giants. Through detailed modeling, we show that these merger remnants possess distinctive photometric, asteroseismic, and surface abundance signatures through which they may be identified. During hydrogen shell burning, merger remnants reach higher luminosities and possess pulsations which depart from the usual degenerate sequence on the asteroseismic $\Delta\nu$--$\Delta\Pi$ diagram for red giant branch stars. For sufficiently massive helium white dwarfs, merger remnants undergo especially violent helium flashes which can dredge up a large amount of core material (up to $\sim0.1M_\odot$) into the envelope. Such post-dredge-up remnants are more luminous than normal red clump stars, are surface carbon-, helium-, and possibly lithium-rich, and possess a wider range of asteroseismic g-mode period spacings and mixed-mode couplings. Recent asteroseismically determined low-mass ($\lesssim0.8M_\odot$) red clump stars may be core helium-burning remnants of mergers involving lower-mass helium white dwarfs., Comment: 24 pages, 16 figures, accepted to the Open Journal of Astrophysics!

Strong magnetic fields are expected to significantly modify the pulsation frequencies of waves propagating in the cores of red giants or in the radiative envelopes of intermediate- and high-mass main-sequence stars. We calculate the g-mode frequencies of stars with magnetic dipole fields which are aligned with their rotational axes, treating both the Lorentz and Coriolis forces non-perturbatively. We provide a compact asymptotic formula for the g-mode period spacing, and universally find that strong magnetism decreases this period spacing substantially more than is predicted by perturbation theory. These results are validated with explicit numerical mode calculations for realistic stellar models. The approach we present is highly versatile: once the eigenvalues $\lambda$ of a certain differential operator are precomputed as a function of the magnetogravity and rotational frequencies (in units of the mode frequency), the non-perturbative impact of the Coriolis and Lorentz forces is understood under a broad domain of validity, and is readily incorporated into asteroseismic modeling., Comment: 18 pages, 8 figures; accepted to MNRAS

When main-sequence stars expand into red giants, they are expected to engulf close-in planets. Until now, the absence of planets with short orbital periods around post-expansion, core-helium-burning red giants has been interpreted as evidence that short-period planets around Sun-like stars do not survive the giant expansion phase of their host stars. Here we present the discovery that the giant planet 8 Ursae Minoris b orbits a core-helium-burning red giant. At a distance of only 0.5 au from its host star, the planet would have been engulfed by its host star, which is predicted by standard single-star evolution to have previously expanded to a radius of 0.7 au. Given the brief lifetime of helium-burning giants, the nearly circular orbit of the planet is challenging to reconcile with scenarios in which the planet survives by having a distant orbit initially. Instead, the planet may have avoided engulfment through a stellar merger that either altered the evolution of the host star or produced 8 Ursae Minoris b as a second-generation planet. This system shows that core-helium-burning red giants can harbour close planets and provides evidence for the role of non-canonical stellar evolution in the extended survival of late-stage exoplanetary systems., Comment: Published in Nature on 28 June 2023. In press

Strong magnetic fields in the cores of stars are expected to significantly modify the behavior of gravity waves: this is likely the origin of suppressed dipole modes observed in many red giants. However, a detailed understanding of how such fields alter the spectrum and spatial structure of magnetogravity waves has been elusive. For a dipole field, we analytically characterize the horizontal eigenfunctions of magnetogravity modes, assuming that the wavevector is primarily radial. For axisymmetric modes ($m=0$), the magnetogravity wave eigenfunctions become Hough functions, and they have a radial turning point for sufficiently strong magnetic fields. For non-axisymmetric modes ($m\neq0$), the interaction between the discrete $g$ mode spectrum and a continuum of Alfv\'en waves produces nearly discontinuous features in the fluid displacements at critical latitudes associated with a singularity in the fluid equations. effects. We find that magnetogravity modes cannot propagate in regions with sufficiently strong magnetic fields, instead becoming evanescent. When encountering strong magnetic fields, ingoing gravity waves are likely refracted into outgoing slow magnetic waves. These outgoing waves approach infinite radial wavenumbers, which are likely to be damped efficiently. However, it may be possible for a small fraction of the wave power to escape the stellar core as pure Alfv\'en waves or magnetogravity waves confined to a very narrow equatorial band. The artificially sharp features in the WKB-separated solutions suggest the need for global mode solutions which include small terms neglected in our analysis., Comment: 22 pages, 10 figures, 1 table, comments welcome! (submitted to MNRAS)

We report CO(5$\rightarrow$4) and CO(6$\rightarrow$5) line observations in the dusty starbursting galaxy CRLE ($z = 5.667$) and the main-sequence (MS) galaxy HZ10 ($z = 5.654$) with the Northern Extended Millimeter Array (NOEMA). CRLE is the most luminous $z>5$ starburst in the COSMOS field and HZ10 is the most gas-rich "normal" galaxy currently known at $z>5$. We find line luminosities for CO(5$\rightarrow$4) and CO(6$\rightarrow$5) of (4.9 $\pm$ 0.5) and (3.8 $\pm$ 0.4) $\times$ 10$^{10}$ K km s$^{-1}$ pc$^{2}$ for CRLE and upper limits of $< 0.76$ and $< 0.60$ $\times$ 10$^{10}$ K km s$^{-1}$ pc$^{2}$ for HZ10, respectively. The CO excitation of CRLE appears comparable to other $z>5$ dusty star-forming galaxies (DSFGs). For HZ10, these line luminosity limits provide the first significant constraints of this kind for a MS galaxy at $z > 5$. We find the upper limit of $L'_{5\rightarrow4}/L'_{2\rightarrow1}$ in HZ10 could be similar to the average value for MS galaxies around $z\approx 1.5$, suggesting that MS galaxies with comparable gas excitation may already have existed one billion years after the Big Bang. For CRLE we determine the most likely values for the H$_2$ density, kinetic temperature and dust temperature based on excitation modeling of the CO line ladder. We also derive a total gas mass of $(7.1 \pm 1.3) \times 10^{10} M_\odot$. Our findings provide some of the currently most detailed constraints on the gas excitation that sets the conditions for star formation in a galaxy protocluster environment at $z > 5$., Comment: 14 pages, 6 figures, accepted to ApJ

The globular cluster 47~Tucanae (47~Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars ($N\gtrsim 10^6$) and high density. Here we present detailed and self-consistent 47~Tuc models computed with the \texttt{Cluster Monte Carlo} code (\texttt{CMC}). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47~Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions \citep[as in, e.g.,][]{Kroupa2001}, and an initial Elson profile \citep{Elson1987} that is overfilling the cluster's tidal radius. We include new prescriptions in \texttt{CMC} for the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47~Tuc; our best-fit model contains $\sim 50$ millisecond pulsars, $70\%$ of which are formed through giant collisions and tidal captures. Our models also suggest that 47~Tuc presently contains up to $\sim 200$ stellar-mass black holes, $\sim 5$ binary black holes, $\sim 15$ low-mass X-ray binaries, and $\sim 300$ cataclysmic variables., Comment: 23 pages, 15 figures, 4 tables, published at ApJ

Stellar mergers are important processes in stellar evolution, dynamics, and transient science. However, it is difficult to identify merger remnant stars because they cannot easily be distinguished from single stars based on their surface properties. We demonstrate that merger remnants can potentially be identified through asteroseismology of red giant stars using measurements of the gravity mode period spacing together with the asteroseismic mass. For mergers that occur after the formation of a degenerate core, remnant stars have over-massive envelopes relative to their cores, which is manifested asteroseismically by a g~mode period spacing smaller than expected for the star's mass. Remnants of mergers which occur when the primary is still on the main sequence or whose total mass is less than $\approx\! 2 \, M_\odot$ are much harder to distinguish from single stars. Using the red giant asteroseismic catalogs of Vrard et al. 2016 and Yu et al. 2018, we identify $24$ promising candidates for merger remnant stars. In some cases, merger remnants could also be detectable using only their temperature, luminosity, and asteroseismic mass, a technique that could be applied to a larger population of red giants without a reliable period spacing measurement., Comment: 13 pages, 6 figures, 1 table, submitted to MNRAS

We describe the public release of the Cluster Monte Carlo Code (CMC) a parallel, star-by-star $N$-body code for modeling dense star clusters. CMC treats collisional stellar dynamics using H\'enon's method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release of CMC is pinned directly to the COSMIC population synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest ($N = 10^8$) star-by-star $N$-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, a recently reported merger with component masses in the pulsational pair-instability mass gap)., Comment: Code is available at https://clustermontecarlo.github.io/ 25 pages, 8 Figures, Matches version accepted by ApJS

Numerical and observational evidence suggests that massive white dwarfs dominate the innermost regions of core-collapsed globular clusters by both number and total mass. Using NGC 6397 as a test case, we constrain the features of white dwarf populations in core-collapsed clusters, both at present day and throughout their lifetimes. The dynamics of these white dwarf subsystems have a number of astrophysical implications. We demonstrate that the collapse of globular cluster cores is ultimately halted by the dynamical burning of white dwarf binaries. We predict core-collapsed clusters in the local universe yield a white dwarf merger rate of $\mathcal{O}(10\rm{)\,Gpc}^{-3}\,\rm{yr}^{-1}$, roughly $0.1-1\%$ of the observed Type Ia supernova rate. We show that prior to merger, inspiraling white dwarf binaries will be observable as gravitational wave sources at milli- and decihertz frequencies. Over $90\%$ of these mergers have a total mass greater than the Chandrasekhar limit. If the merger/collision remnants are not destroyed completely in an explosive transient, we argue the remnants may be observed in core-collapsed clusters as either young neutron stars/pulsars/magnetars (in the event of accretion-induced collapse) or as young massive white dwarfs offset from the standard white dwarf cooling sequence. Finally, we show collisions between white dwarfs and main sequence stars, which may be detectable as bright transients, occur at a rate of $\mathcal{O}(100\rm{)\,Gpc}^{-3}\,\rm{yr}^{-1}$ in the local universe. We find that these collisions lead to depletion of blue straggler stars and main sequence star binaries in the centers of core-collapsed clusters., Comment: Submitted to ApJ, 26 pages, 12 figures, 3 tables. Comments welcome

Recently, \citet{vitral2021does} detected a central concentration of dark objects in the core-collapsed globular cluster NGC 6397, which could be interpreted as a subcluster of stellar-mass black holes. However, it is well established theoretically that any significant number of black holes in the cluster would provide strong dynamical heating and is fundamentally inconsistent with this cluster's core-collapsed profile. Claims of intermediate-mass black holes in core-collapsed clusters should similarly be treated with suspicion, for reasons that have been understood theoretically for many decades. Instead, the central dark population in NGC 6397 is exactly accounted for by a compact subsystem of white dwarfs, as we demonstrate here by inspection of a previously published model that provides a good fit to this cluster. These central subclusters of heavy white dwarfs are in fact a generic feature of core-collapsed clusters, while central black hole subclusters are present in all {\em non\/}-collapsed clusters., Comment: 3 pages, 1 figure, to submit to AAS journals

As ancient, gravitationally bound stellar populations, globular clusters are abundant, vibrant laboratories characterized by high frequencies of dynamical interactions coupled to complex stellar evolution. Using surface brightness and velocity dispersion profiles from the literature, we fit $59$ Milky Way globular clusters to dynamical models from the \texttt{CMC Cluster Catalog}. Without doing any interpolation, and without any directed effort to fit any particular cluster, $26$ globular clusters are well-matched by at least one of our models. We discuss in particular the core-collapsed clusters NGC 6293, NGC 6397, NGC 6681, and NGC 6624, and the non-core-collapsed clusters NGC 288, NGC 4372, and NGC 5897. As NGC 6624 lacks well-fitting snapshots on the main \texttt{CMC Cluster Catalog}, we run six additional models in order to refine the fit. We calculate metrics for mass segregation, explore the production of compact object sources such as millisecond pulsars, cataclysmic variables, low-mass X-ray binaries, and stellar-mass black holes, finding reasonable agreement with observations. Additionally, closely mimicking observational cuts, we extract the binary fraction from our models, finding good agreement except in the dense core regions of core-collapsed clusters. Accompanying this paper are a number of \textsf{python} methods for examining the publicly accessible \texttt{CMC Cluster Catalog}, as well as any other models generated using \texttt{CMC}., Comment: 22 pages, 11 figures, 5 tables; accepted to ApJ

In this study, we develop a methodology for model reduction and selection informed by global sensitivity analysis (GSA) methods. We apply these techniques to a control model that takes systolic blood pressure and thoracic tissue pressure data as inputs and predicts heart rate in response to the Valsalva maneuver (VM). The study compares four GSA methods based on Sobol' indices (SIs) quantifying the parameter influence on the difference between the model output and the heart rate data. The GSA methods include standard scalar SIs determining the average parameter influence over the time interval studied and three time-varying methods analyzing how parameter influence changes over time. The time-varying methods include a new technique, termed limited-memory SIs, predicting parameter influence using a moving window approach. Using the limited-memory SIs, we perform model reduction and selection to analyze the necessity of modeling both the aortic and carotid baroreceptor regions in response to the VM. We compare the original model to three systematically reduced models including (i) the aortic and carotid regions, (ii) the aortic region only, and (iii) the carotid region only. Model selection is done quantitatively using the Akaike and Bayesian Information Criteria and qualitatively by comparing the neurological predictions. Results show that it is necessary to incorporate both the aortic and carotid regions to model the VM.

We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models to investigate the role that stellar-mass black holes play in the process of core collapse. Furthermore, we study how dynamical interactions affect the formation and evolution of several important types of sources in GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers, cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging compact binaries. While our focus here is on old, low-metallicity GCs, our $\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time, and they include a wide range of metallicities (up to solar), so that our results can also be used to study younger and higher-metallicity star clusters., Comment: 50 pages, 15 figures, 10 tables. Accepted for publication in ApJS

Delay differential equations (DDEs) are widely used in mathematical modeling to describe physical and biological systems. Delays can impact model dynamics, resulting in oscillatory behavior. In physiological systems, this instability may signify (i) an attempt to return to homeostasis or (ii) system dysfunction. In this study, we analyze a nonlinear, nonautonomous, nonhomogeneous open-loop neurological control model describing the autonomic nervous system response to the Valsalva maneuver. Unstable modes have been identified as a result of parameter interactions between the sympathetic delay and time-scale. In a two-parameter bifurcation analysis, we examine both the homogeneous and nonhomogeneous systems. Discrepancies between solutions result from the presence of the forcing functions which stabilize the system. We use analytical methods to determine stability regions for the homogeneous system, identifying transcendental relationships between the parameters. We also use computational methods to determine stability regions for the nonhomogeneous system. The presence of a Hopf bifurcation within the system is discussed and solution types from the sink and stable focus regions are compared to two control patients and a patient with postural orthostatic tachycardia syndrome (POTS). The model and its analysis support the current clinical hypotheses that patients suffering from POTS experience altered nervous system activity., Comment: Mathematical Biosciences, 2019

The Quintuplet star cluster is one of only three known young ($<10$ Myr) massive (M $>10^4$ M$_\odot$) clusters within $\sim100$ pc of the Galactic Center. In order to explore star cluster formation and evolution in this extreme environment, we analyze the Quintuplet's dynamical structure. Using the HST WFC3-IR instrument, we take astrometric and photometric observations of the Quintuplet covering a $120''\times120''$ field-of-view, which is $19$ times larger than those of previous proper motion studies of the Quintuplet. We generate a catalog of the Quintuplet region with multi-band, near-infrared photometry, proper motions, and cluster membership probabilities for $10,543$ stars. We present the radial density profile of $715$ candidate Quintuplet cluster members with $M\gtrsim4.7$ M$_\odot$ out to $3.2$ pc from the cluster center. A $3\sigma$ lower limit of $3$ pc is placed on the tidal radius, indicating the lack of a tidal truncation within this radius range. Only weak evidence for mass segregation is found, in contrast to the strong mass segregation found in the Arches cluster, a second and slightly younger massive cluster near the Galactic Center. It is possible that tidal stripping hampers a mass segregation signature, though we find no evidence of spatial asymmetry. Assuming that the Arches and Quintuplet formed with comparable extent, our measurement of the Quintuplet's comparatively large core radius of $0.62^{+0.10}_{-0.10}$ pc provides strong empirical evidence that young massive clusters in the Galactic Center dissolve on a several Myr timescale., Comment: 25 pages (21-page main text, 4-page appendix), 18 figures, submitted to ApJ

We report the first detailed measurement of the shape of the CO luminosity function at high redshift, based on $>$320 hr of the NSF's Karl G. Jansky Very Large Array (VLA) observations over an area of $\sim$60 arcmin$^2$ taken as part of the CO Luminosity Density at High Redshift (COLDz) survey. COLDz "blindly" selects galaxies based on their cold gas content through CO($J$=1$\to$0) emission at $z$$\sim$2-3 and CO($J$=2$\to$1) at $z$$\sim$5-7 down to a CO luminosity limit of log($L'_{\rm CO}$/K km s$^{-1}$ pc$^2$)$\simeq$9.5. We find that the characteristic luminosity and bright end of the CO luminosity function are substantially higher than predicted by semi-analytical models, but consistent with empirical estimates based on the infrared luminosity function at $z$$\sim$2. We also present the currently most reliable measurement of the cosmic density of cold gas in galaxies at early epochs, i.e., the cold gas history of the universe, as determined over a large cosmic volume of $\sim$375,000 Mpc$^3$. Our measurements are in agreement with an increase of the cold gas density from $z$$\sim$0 to $z$$\sim$2-3, followed by a possible decline towards $z$$\sim$5-7. These findings are consistent with recent surveys based on higher-$J$ CO line measurements, upon which COLDz improves in terms of statistical uncertainties by probing $\sim$50-100 times larger areas and in the reliability of total gas mass estimates by probing the low-$J$ CO lines accessible to the VLA. Our results thus appear to suggest that the cosmic star-formation rate density follows an increased cold molecular gas content in galaxies towards its peak about 10 billion years ago, and that its decline towards the earliest epochs is likely related to a lower overall amount of cold molecular gas (as traced by CO) bound in galaxies towards the first billion years after the Big Bang., Comment: 16 pages, 6 figures, 6 tables, to appear in ApJ (accepted January 04, 2019)

We describe the CO Luminosity Density at High-z (COLDz) survey, the first spectral line deep field targeting CO(1-0) emission from galaxies at $z=1.95-2.85$ and CO(2-1) at $z=4.91-6.70$. The main goal of COLDz is to constrain the cosmic density of molecular gas at the peak epoch of cosmic star formation. By targeting both a wide ($\sim$51 arcmin$^2$) and a deep area ($\sim$9 arcmin$^2$), the survey is designed to robustly constrain the bright end and the characteristic luminosity of the CO(1-0) luminosity function. An extensive analysis of the reliability of our line candidates, and new techniques provide detailed completeness and statistical corrections as necessary to determine the best constraints to date on the CO luminosity function. Our blind search for CO(1-0) uniformly selects starbursts and massive Main Sequence galaxies based on their cold molecular gas masses. Our search also detects CO(2-1) line emission from optically dark, dusty star-forming galaxies at $z>5$. We find a range of spatial sizes for the CO-traced gas reservoirs up to $\sim40$ kpc, suggesting that spatially extended cold molecular gas reservoirs may be common in massive, gas-rich galaxies at $z\sim2$. Through CO line stacking, we constrain the gas mass fraction in previously known typical star-forming galaxies at $z=2$-3. The stacked CO detection suggests lower molecular gas mass fractions than expected for massive Main Sequence galaxies by a factor of $\sim3-6$. We find total CO line brightness at $\sim34\,$GHz of $0.45\pm0.2\,\mu$K, which constrains future line intensity mapping and CMB experiments., Comment: 53 pages, 22 figures, 8 tables, ApJ in press. Full resolution images and publicly available data may be found on the survey webpage at coldz.astro.cornell.edu

We report interferometric measurements of [NII] 205 um fine-structure line emission from a representative sample of three galaxies at z=5-6 using the Atacama Large (sub)Millimeter Array (ALMA). These galaxies were previously detected in [CII] and far-infrared continuum emission and span almost two orders of magnitude in star formation rate (SFR). Our results show at least two different regimes of ionized inter-stellar medium properties for galaxies in the first billion years of cosmic time, separated by their L_[CII]/L_[NII] ratio. We find extremely low [NII] emission compared to [CII] (L_ [CII]/L_[NII]=68 [+200/-28]) from a "typical" L*_UV star-forming galaxy, likely directly or indirectly (by its effect on the radiation field) related to low dust abundance and low metallicity. The infrared-luminous modestly star-forming Lyman Break Galaxy (LBG) in our sample is characterized by an ionized-gas fraction (L_[CII]/L_[NII]<=20) typical of local star-forming galaxies and shows evidence for spatial variations in its ionized-gas fraction across an extended gas reservoir. The extreme SFR, warm and compact dusty starburst AzTEC-3 shows an ionized fraction higher than expected given its star-formation rate surface density (L_[CII]/L_[NII]=22+/-8) suggesting that [NII] dominantly traces a diffuse ionized medium rather than star-forming HII regions in this type of galaxy. This highest redshift sample of [NII] detections provides some of the first constraints on ionized and neutral gas modeling attempts and on the structure of the inter-stellar medium at z=5-6 in "normal" galaxies and starbursts., Comment: 12 pages, 5 figures, 2 tables, to appear in ApJ (accepted September 15, 2016)

We investigate the global cold dust properties of 85 nearby (z < 0.5) QSOs, chosen from the Palomar-Green sample of optically luminous quasars. We determine their infrared spectral energy distributions and estimate their rest-frame luminosities by combining Herschel data from 70 to 500 microns with near-infrared and mid-infrared measurements from the Two Micron All Sky Survey (2MASS) and the Wide-Field Infrared Survey Explorer (WISE). In most sources the far-infrared (FIR) emission can be attributed to thermally heated dust. Single temperature modified black body fits to the FIR photometry give an average dust temperature for the sample of 33~K, with a standard deviation of 8~K, and an average dust mass of 7E6 Solar Masses with a standard deviation of 9E6 Solar Masses. Estimates of star-formation that are based on the FIR continuum emission correlate with those based on the 11.3 microns PAH feature, however, the star-formation rates estimated from the FIR continuum are higher than those estimated from the 11.3 microns PAH emission. We attribute this result to a variety of factors including the possible destruction of the PAHs and that, in some sources, a fraction of the FIR originates from dust heated by the active galactic nucleus and by old stars., Comment: accepted for publication in ApJS

We report interferometric imaging of [CII] and OH emission toward the center of the galaxy protocluster associated with the z=5.3 submillimeter galaxy (SMG) AzTEC-3, using the Atacama Large (sub)Millimeter Array (ALMA). We detect strong [CII], OH, and rest-frame 157.7 um continuum emission toward the SMG. The [CII] emission is distributed over a scale of 3.9 kpc, implying a dynamical mass of 9.7 x 10^10 Msun, and a star formation rate (SFR) surface density of Sigma_SFR = 530 Msun/yr/kpc2. This suggests that AzTEC-3 forms stars at Sigma_SFR approaching the Eddington limit for radiation pressure supported disks. We find that the OH emission is slightly blueshifted relative to the [CII] line, which may indicate a molecular outflow associated with the peak phase of the starburst. We also detect and dynamically resolve [CII] emission over a scale of 7.5 kpc toward a triplet of Lyman-break galaxies with moderate UV-based SFRs in the protocluster at ~95kpc projected distance from the SMG. These galaxies are not detected in the continuum, suggesting far-infrared SFRs of <18-54 Msun/yr, consistent with a UV-based estimate of 22 Msun/yr. The spectral energy distribution of these galaxies is inconsistent with nearby spiral and starburst galaxies, but resembles those of dwarf galaxies. This is consistent with expectations for young starbursts without significant older stellar populations. This suggests that these galaxies are significantly metal-enriched, but not heavily dust-obscured, "normal" star-forming galaxies at z>5, showing that ALMA can detect the interstellar medium in "typical" galaxies in the very early universe., Comment: 15 pages, 12 figures, 4 tables, to appear in ApJ (accepted October 15, 2014)

The HII region luminosity function (LF) is an important tool for deriving the birthrates and mass distribution of OB associations, and is an excellent tracer of the newly formed massive stars and associations. To date, extensive work (predominantly in H{\alpha}) has been done from the ground, which is hindered by dust extinction and the severe blending of adjacent (spatially or in projection) HII regions. Reliably measuring the properties of HII regions requires a linear resolution <40 pc, but analyses satisfying this requirement have been done only in a handful of galaxies, so far. As the first space-based work using a galaxy sample, we have selected 12 galaxies from our HST NICMOS Pa{\alpha} survey and studied the luminosity function and size distribution of HII regions both in individual galaxies and cumulatively, using a virtually extinction-free tracer of the ionizing photon rate. The high angular resolution and low sensitivity to diffuse emission of NICMOS also offer an advantage over ground-based imaging by enabling a higher degree of de-blending of the HII regions. We do not confirm the broken power-law LFs found in ground-based studies. Instead, we find that the LFs, both individual and co-added, follow a single power law dN(L)/dlnL ~ L^-1, consistent with the mass function of star clusters in nearby galaxies, and in agreement with the results of the existing analyses with HST data. The individual and co-added size distribution of HII regions are both roughly consistent with dN(D)/dlnD ~ D^-3, but the power-law scaling is probably contaminated by blended regions or complexes., Comment: 17 pages, 11 figures, 2 tables, accepted for publication in ApJ. The full catalog of the 1907 HII regions is available at: http://www.pha.jhu.edu/~liu/papers/h2lf.tab2.dat