We derive the UV-optical stellar dust attenuation curve of galaxies at z=1.4-2.6 as a function of gas-phase metallicity. We use a sample of 218 star-forming galaxies, excluding those with very young or heavily obscured star formation, from the MOSFIRE Deep Evolution Field (MOSDEF) survey with H$\alpha$, H$\beta$, and [NII]$\lambda 6585$ spectroscopic measurements. We constrain the shape of the attenuation curve by comparing the average flux densities of galaxies sorted into bins of dust obscuration using Balmer decrements, i.e., H$\alpha$-to-H$\beta$ luminosities. The average attenuation curve for the high-metallicity sample (12+log(O/H)>8.5, corresponding to $M_*\gtrsim10^{10.4}\,M_{\odot}$) has a shallow slope, identical to that of the Calzetti local starburst curve, and a significant UV 2175A extinction bump that is $\sim 0.5\times$ the strength of the Milky Way bump. On the other hand, the average attenuation curve of the low-metallicity sample (12+log(O/H) $\sim 8.2-8.5$) has a steeper slope similar to that of the SMC curve, only consistent with the Calzetti slope at the $3\sigma$ level. The UV bump is not detected in the low-metallicity curve, indicating the relative lack of the small dust grains causing the bump at low metallicities. Furthermore, we find that on average the nebular reddening (E(B-V)) is a factor of 2 times larger than that of the stellar continuum for galaxies with low metallicities, while the nebular and stellar reddening are similar for galaxies with higher metallicities. The latter is likely due to a high surface density of dusty clouds embedding the star forming regions but also reddening the continuum in the high-metallicity galaxies.

We provide a coherent, uniform measurement of the evolution of the logarithmic star formation rate (SFR) - stellar mass ($M_*$) relation, called the main sequence of star-forming galaxies (MS), for galaxies out to $z\sim5$. We measure the MS using mean stacks of 3 GHz radio continuum images to derive average SFRs for $\sim$200,000 mass-selected galaxies at $z>0.3$ in the COSMOS field. We describe the MS relation adopting a new model that incorporates a linear relation at low stellar mass (log($M_*$/M$_\odot$)$<$10) and a flattening at high stellar mass that becomes more prominent at low redshift ($z<1.5$). We find that the SFR density peaks at $1.5<z<2$ and at each epoch there is a characteristic stellar mass ($M_* = 1 - 4 \times 10^{10}\mathrm{M}_\odot$) that contributes the most to the overall SFR density. This characteristic mass increases with redshift, at least to $z\sim2.5$. We find no significant evidence for variations in the MS relation for galaxies in different environments traced by the galaxy number density at $0.3<z<3$, nor for galaxies in X-ray groups at $z\sim0.75$. We confirm that massive bulge-dominated galaxies have lower SFRs than disk-dominated galaxies at a fixed stellar mass at $z<1.2$. As a consequence, the increase in bulge-dominated galaxies in the local star-forming population leads to a flattening of the MS at high stellar masses. This indicates that "mass-quenching" is linked with changes in the morphological composition of galaxies at a fixed stellar mass.

We update the spectral modeling code MAGPHYS to include a 2175\AA\ absorption feature in its UV-to-near-IR dust attenuation prescription. This allows us to determine the strength of this feature and the shape of the dust attenuation curve in ~5000 star-forming galaxies at 0.1<z<3 in the COSMOS field. We find that a 2175\AA\ absorption feature of ~1/3 the strength of that in the Milky Way is required for models to minimize residuals. We characterize the total effective dust attenuation curves as a function of several galaxy properties and find that the UV slopes of the attenuation curve for COSMOS galaxies show a strong dependence with star formation rate (SFR) and total dust attenuation ($A_V$), such that galaxies with higher SFR and $A_V$ have shallower curves and vice versa. These results are consistent with expectations from radiative transfer that attenuation curves become shallower as the effective dust optical depth increases. We do not find significant trends in the strength of the 2175\AA\ absorption feature as a function of galaxy properties, but this may result from the high uncertainties associated with this measurement. The updated code is publicly available online.

The rich information on (sub)millimeter dust continuum emission from distant galaxies in the public Atacama Large Millimeter/submillimeter Array (ALMA) archive is contained in thousands of inhomogeneous observations from individual PI-led programs. To increase the usability of these data for studies deepening our understanding of galaxy evolution, we have developed automated mining pipelines for the ALMA archive in the COSMOS field (A3COSMOS) that efficiently exploit the available information for large numbers of galaxies across cosmic time and keep the data products in sync with the increasing public ALMA archive: (a) a dedicated ALMA continuum imaging pipeline, (b) two complementary photometry pipelines for both blind source extraction and prior source fitting, (c) a counterpart association pipeline utilizing the multiwavelength data available (including quality assessment based on machine-learning techniques), (d) an assessment of potential (sub)millimeter line contribution to the measured ALMA continuum, and (e) extensive simulations to provide statistical corrections to biases and uncertainties in the ALMA continuum measurements. Application of these tools yields photometry catalogs with ~1000 (sub)millimeter detections (spurious fraction ~8%-12%) from over 1500 individual ALMA continuum images. Combined with ancillary photometric and redshift catalogs and the above quality assessments, we provide robust information on redshift, stellar mass, and star formation rate for ~700 galaxies at redshifts 0.5-6 in the COSMOS field (with undetermined selection function). The ALMA photometric measurements and galaxy properties are released publicly within our blind extraction, prior fitting, and galaxy property catalogs, plus the images. These products will be updated on a regular basis in the future.

The largest Herschel extragalactic surveys, H-ATLAS and HerMES, have selected a sample of "ultrared" dusty, star-forming galaxies (DSFGs) with rising SPIRE flux densities ($S_{500} > S_{350} > S_{250}$; so-called "500 $\mu$m-risers") as an efficient way for identifying DSFGs at higher redshift ($z > 4$). In this paper, we present a large Spitzer follow-up program of 300 Herschel ultrared DSFGs. We have obtained high-resolution ALMA, NOEMA, and SMA data for 63 of them, which allow us to securely identify the Spitzer/IRAC counterparts and classify them as gravitationally lensed or unlensed. Within the 63 ultrared sources with high-resolution data, $\sim$65% appear to be unlensed, and $\sim$27% are resolved into multiple components. We focus on analyzing the unlensed sample by directly performing multi-wavelength spectral energy distribution (SED) modeling to derive their physical properties and compare with the more numerous $z \sim 2$ DSFG population. The ultrared sample has a median redshift of 3.3, stellar mass of 3.7 $\times$ 10$^{11}$ $M_{\odot}$, star formation rate (SFR) of 730 $M_{\odot}$yr$^{-1}$, total dust luminosity of 9.0 $\times$ 10$^{12}$ $L_{\odot}$, dust mass of 2.8 $\times$ 10$^9$ $M_{\odot}$, and V-band extinction of 4.0, which are all higher than those of the ALESS DSFGs. Based on the space density, SFR density, and stellar mass density estimates, we conclude that our ultrared sample cannot account for the majority of the star-forming progenitors of the massive, quiescent galaxies found in infrared surveys. Our sample contains the rarer, intrinsically most dusty, luminous and massive galaxies in the early universe that will help us understand the physical drivers of extreme star formation.

We present an enhanced version of the multiwavelength spectral modeling code MAGPHYS that allows the estimation of galaxy photometric redshift and physical properties (e.g., stellar mass, star formation rate, dust attenuation) simultaneously, together with robust characterization of their uncertainties. The self-consistent modeling over ultraviolet to radio wavelengths in MAGPHYS+photo-z is unique compared to standard photometric redshift codes. The broader wavelength consideration is particularly useful for breaking certain degeneracies in color vs. redshift for dusty galaxies with limited observer-frame ultraviolet and optical data (or upper limits). We demonstrate the success of the code in estimating redshifts and physical properties for over 4,000 infrared-detected galaxies at 0.4<z<6.0 in the COSMOS field with robust spectroscopic redshifts. We achieve high photo-z precision ($\sigma_{\Delta z/(1+z_{spec})}\lesssim0.04$), high accuracy (i.e., minimal offset biases; median$(\Delta z/(1+z_{spec}))\lesssim0.02$), and low catastrophic failure rates ($\eta\simeq4\%$) over all redshifts. Interestingly, we find that a weak 2175A absorption feature in the attenuation curve models is required to remove a subtle systematic photo-z offset ($z_{phot}-z_{spec}\simeq-0.03$) that occurs when this feature is not included. As expected, the accuracy of derived physical properties in MAGPHYS+photo-z decreases strongly as redshift uncertainty increases. The all-in-one treatment of uncertainties afforded with this code is beneficial for accurately interpreting physical properties of galaxies in large photometric datasets. Finally, we emphasize that MAGPHYS+photo-z is not intended to replace existing photo-z codes, but rather offer flexibility to robustly interpret physical properties when spectroscopic redshifts are unavailable. The MAGPHYS+photo-z code is publicly available online.

We investigate the relation between gas and star formation in sub-galactic regions, ~360 pc to ~1.5 kpc in size, within the nearby starburst dwarf NGC4449, in order to separate the underlying relation from the effects of sampling at varying spatial scales. Dust and gas mass surface densities are derived by combining new observations at 1.1 mm, obtained with the AzTEC instrument on the Large Millimeter Telescope, with archival infrared images in the range 8-500 micron from the Spitzer Space Telescope and the Herschel Space Observatory. We extend the dynamic range of our mm (and dust) maps at the faint end, using a correlation between the far-infrared/millimeter colors F(70)/F(1100) [and F(160)/F(1100)] and the mid-infrared color F(8)/F(24) that we establish for the first time for this and other galaxies. Supplementing our data with maps of the extinction-corrected star formation rate (SFR) surface density, we measure both the SFR-molecular gas and the SFR-total gas relations in NGC4449. We find that the SFR-molecular gas relation is described by a power law with exponent that decreases from ~1.5 to ~1.2 for increasing region size, while the exponent of the SFR-total gas relation remains constant with value ~1.5 independent of region size. We attribute the molecular law behavior to the increasingly better sampling of the molecular cloud mass function at larger region sizes; conversely, the total gas law behavior likely results from the balance between the atomic and molecular gas phases achieved in regions of active star formation. Our results indicate a non-linear relation between SFR and gas surface density in NGC4449, similar to what is observed for galaxy samples.

We characterize the influence that inclination has on the shape and normalization in average dust attenuation curves derived from a sample of ~10,000 local star-forming galaxies. To do this we utilize aperture-matched multi-wavelength data from the Galaxy Evolution Explorer, the Sloan Digital Sky Survey, the United Kingdom Infrared Telescope, and the Two Micron All-Sky Survey. We separate our sample into groups according to axial ratio (b/a) and obtain an estimate of their average total-to-selective attenuation $k(\lambda)$. The attenuation curves are found to be shallower at UV wavelengths with increasing inclination, whereas the shape at longer wavelengths remains unchanged. The highest inclination subpopulation (b/a<0.42) exhibits a NUV excess in its average selective attenuation, which, if interpreted as a 2175\AA\ feature, is best fit with a bump strength of 17-26% of the MW value. No excess is apparent in the average attenuation curve of lower inclination galaxies. The differential reddening between the stellar continuum and ionized gas is found to decrease with increasing inclination. We find that higher inclination galaxies have slightly higher values of $R_V$, although this is poorly constrained given the uncertainties. We outline possible explanations for these trends within a two component dust model (dense cloud+diffuse dust) and find that they can be naturally explained if carriers of the 2175\AA\ feature are preferentially destroyed in star-forming regions (UV-bright regions).

We characterize the near-infrared (NIR) dust attenuation for a sample of ~5500 local (z<0.1) star-forming galaxies and obtain an estimate of their average total-to-selective attenuation $k(\lambda)$. We utilize data from the United Kingdom Infrared Telescope (UKIRT) and the Two Micron All-Sky Survey (2MASS), which is combined with previously measured UV-optical data for these galaxies. The average attenuation curve is slightly lower in the far-UV than local starburst galaxies, by roughly 15%, but appears similar at longer wavelengths with a total-to-selective normalization at V-band of $R_V=3.67\substack{+0.44 \\ -0.35}$. Under the assumption of energy balance, the total attenuated energy inferred from this curve is found to be broadly consistent with the observed infrared dust emission ($L_{\rm{TIR}}$) in a small sample of local galaxies for which far-IR measurements are available. However, the significant scatter in this quantity among the sample may reflect large variations in the attenuation properties of individual galaxies. We also derive the attenuation curve for sub-populations of the main sample, separated according to mean stellar population age (via $D_n4000$), specific star formation rate, stellar mass, and metallicity, and find that they show only tentative trends with low significance, at least over the range which is probed by our sample. These results indicate that a single curve is reasonable for applications seeking to broadly characterize large samples of galaxies in the local Universe, while applications to individual galaxies would yield large uncertainties and is not recommended.

We directly detect dust emission in an optically-detected, multiply-imaged galaxy lensed by the Frontier Fields cluster MACSJ0717.5+3745. We detect two images of the same galaxy at 1.1mm with the AzTEC camera on the Large Millimeter Telescope leaving no ambiguity in the counterpart identification. This galaxy, MACS071_Az9, is at z>4 and the strong lensing model (mu=7.5) allows us to calculate an intrinsic IR luminosity of 9.7e10 Lsun and an obscured star formation rate of 14.6 +/- 4.5 Msun/yr. The unobscured star formation rate from the UV is only 4.1 +/- 0.3 Msun/yr which means the total star formation rate (18.7 +/- 4.5 Msun/yr) is dominated (75-80%) by the obscured component. With an intrinsic stellar mass of only 6.9e9Msun, MACS0717_Az9 is one of only a handful of z>4 galaxies at these lower masses that is detected in dust emission. This galaxy lies close to the estimated star formation sequence at this epoch. However, it does not lie on the dust obscuration relation (IRX-beta) for local starburst galaxies and is instead consistent with the Small Magellanic Cloud (SMC) attenuation law. This remarkable lower mass galaxy showing signs of both low metallicity and high dust content may challenge our picture of dust production in the early Universe.

The tidal disruption of a star by a massive black hole (MBH) is thought to produce a transient luminous event. Such tidal disruption events (TDEs) may play an important role in the detection and characterization of MBHs and probe the properties and dynamics of their nuclear stellar clusters (NSCs) hosts. Previous studies estimated the recent rates of TDEs in the local universe. However, the long-term evolution of the rates throughout the history of the universe has been hardly explored. Here we consider the TDE history, using evolutionary models for the evolution of galactic nuclei. We use a 1D Fokker-Planck approach to explore the evolution of MBH-hosting NSCs, and obtain the disruption rates of stars during their evolution. We complement these with an analysis of TDEs history based on N-body simulation data, and find them to be comparable. We consider NSCs that are built-up from close-in star-formation (SF) or from SF/clusters-dispersal far-out, a few pc from the MBH. We also explore cases where primordial NSCs exist and later evolve through additional star-formation/cluster-dispersal processes. We study the dependence of the TDE history on the type of galaxy, as well as the dependence on the MBH mass. These provide several scenarios, with a continuous increase of the TDE rates over time for cases of far-out SF and a more complex behavior for the close-in SF cases. Finally, we integrate the TDE histories of the various scenarios to provide a total TDE history of the universe, which can be potentially probed with future large surveys (e.g. LSST).

The dust attenuation for a sample of $\sim$10000 local ($z\lesssim0.1$) star forming galaxies is constrained as a function of their physical properties. We utilize aperture-matched multi-wavelength data available from the Galaxy Evolution Explorer (GALEX) and the Sloan Digital Sky Survey (SDSS) to ensure that regions of comparable size in each galaxy are being analyzed. We follow the method of Calzetti et al. (1994) and characterize the dust attenuation through the UV power-law index, $\beta$, and the dust optical depth, which is quantified using the difference in Balmer emission line optical depth, $\tau_B^l=\tau_{\mathrm{H}\beta} - \tau_{\mathrm{H}\alpha}$. The observed linear relationship between $\beta$ and $\tau_B^l$ is similar to the local starburst relation, but the large scatter ($\sigma_{\mathrm{int}}=0.44$) suggests there is significant variation in the local Universe. We derive a selective attenuation curve over the range $1250\mathrm{\AA}<\lambda<8320\mathrm{\AA}$ and find that a single attenuation curve is effective for characterizing the majority of galaxies in our sample. This curve has a slightly lower selective attenuation in the UV compared to previously determined curves. We do not see evidence to suggest that a 2175\AA\ feature is significant in the average attenuation curve. Significant positive correlations are seen between the amount of UV and optical reddening and galaxy metallicity, mass, star formation rate (SFR), and SFR surface density. This provides a potential tool for gauging attenuation where the stellar population is unresolved, such as at high-$z$.

We present continuous, monochromatic star formation rate (SFR) indicators over the mid-infrared wavelength range of 6-70 micron. We use a sample of 58 star forming galaxies (SFGs) in the Spitzer-SDSS-GALEX Spectroscopic Survey (SSGSS) at z<0.2, for which there is a rich suite of multi-wavelength photometry and spectroscopy from the ultraviolet through to the infrared. The data from the Spitzer infrared spectrograph (IRS) of these galaxies, which spans 5-40 micron, is anchored to their photometric counterparts. The spectral region between 40-70 micron is interpolated using dust model fits to the IRS spectrum and Spitzer 70 and 160 micron photometry. Since there are no sharp spectral features in this region, we expect these interpolations to be robust. This spectral range is calibrated as a SFR diagnostic using several reference SFR indicators to mitigate potential bias. Our band-specific continuous SFR indicators are found to be consistent with monochromatic calibrations in the local universe, as derived from Spitzer, WISE, and Herschel photometry. Our local composite template and continuous SFR diagnostics are made available for public use through the NASA/IPAC Infrared Science Archive (IRSA) and have typical dispersions of 30% or less. We discuss the validity and range of applicability for our SFR indicators in the context of unveiling the formation and evolution of galaxies. Additionally, in the era of the James Webb Space Telescope this will become a flexible tool, applicable to any SFG up to z~3.

We present a new implementation of the Monte Carlo method to simulate the evolution of star clusters. The major improvement with respect to the previously developed codes is the treatment of the external tidal field taking into account for both the loss of stars from the cluster boundary and the disk/bulge shocks. We provide recipes to handle with eccentric orbits in complex galactic potentials. The first calculations for stellar systems containing 21000 and 42000 equal-mass particles show good agreement with direct N-body simulations in terms of the evolution of both the enclosed mass and the Lagrangian radii provided that the mass-loss rate does not exceed a critical value.

In this paper, we consider how gas damping affects the dynamical evolution of gas-embedded star clusters. Using a simple three-component (i.e. one gas and two stellar components) model, we compare the rates of mass segregation due to two-body relaxation, accretion from the interstellar medium, and gas dynamical friction in both the supersonic and subsonic regimes. Using observational data in the literature, we apply our analytic predictions to two different astrophysical environments, namely galactic nuclei and young open star clusters. Our analytic results are then tested using numerical simulations performed with the NBSymple code, modified by an additional deceleration term to model the damping effects of the gas. The results of our simulations are in reasonable agreement with our analytic predictions, and demonstrate that gas damping can significantly accelerate the rate of mass segregation. A stable state of approximate energy equilibrium cannot be achieved in our model if gas damping is present, even if Spitzer's Criterion is satisfied. This instability drives the continued dynamical decoupling and subsequent ejection (and/or collisions) of the more massive population. Unlike two-body relaxation, gas damping causes overall cluster contraction, reducing both the core and half-mass radii. If the cluster is mass segregated (and/or the gas density is highest at the cluster centre), the latter contracts faster than the former, accelerating the rate of core collapse.

The mass fraction of dense gas within giant molecular clouds (GMCs) of the Milky Way is investigated using 13CO data from the FCRAO Galactic Plane Surveys and the Bolocam Galactic Plane Survey (BGPS) of 1.1 mm dust continuum emission. A sample of 860 compact dust sources are selected from the BGPS catalog and kinematically linked to 344 clouds of extended (>3') 13CO J=1-0 emission. Gas masses are tabulated for the full dust source and subregions within the dust sources with mass surface densities greater than 200 M_solar/pc^2, which are assumed to be regions of enhanced volume density. Masses of the parent GMCs are calculated assuming optically thin 13CO J=1-0 emission and LTE conditions. The mean fractional mass of dust sources to host GMC mass is 0.11 (-0.06, +0.12). The high column density subregions comprise 0.07 (-0.05, +0.13) of the mass of the cloud. Owing to our assumptions, these values are upper limits to the true mass fractions. The fractional mass of dense gas is independent of GMC mass and gas surface density. The low dense gas mass fraction suggests that the formation of dense structures within GMCs is the primary bottleneck for star formation. The distribution of velocity differences between the dense gas and the low density material along the line of sight is also examined. We find a strong, centrally peaked distribution centered on zero velocity displacement. This distribution of velocity differences is modeled with radially converging flows towards the dense gas position that are randomly oriented with respect to the observed line of sight. These models constrain the infall velocities to be 2-4 km/s for various flow configurations.

We introduce the Stony Brook / SMARTS Atlas of (mostly) Southern Novae. This atlas contains both spectra and photometry obtained since 2003. The data archived in this atlas will facilitate systematic studies of the nova phenomenon and correlative studies with other comprehensive data sets. It will also enable detailed investigations of individual objects. In making the data public we hope to engender more interest on the part of the community in the physics of novae. The atlas is on-line at \url{http://www.astro.sunysb.edu/fwalter/SMARTS/NovaAtlas/} .

The application of high end computing to astrophysical problems, mainly in the galactic environment, is under development since many years at the Dep. of Physics of Sapienza Univ. of Roma. The main scientific topic is the physics of self gravitating systems, whose specific subtopics are: i) celestial mechanics and interplanetary probe transfers in the solar system; ii) dynamics of globular clusters and of globular cluster systems in their parent galaxies; iii) nuclear clusters formation and evolution; iv) massive black hole formation and evolution; v) young star cluster early evolution. In this poster we describe the software and hardware computational resources available in our group and how we are developing both software and hardware to reach the scientific aims above itemized.

We present Cosmic Origins Spectrograph (COS) measurements of metal abundances in eight 0.083<z<0.321 damped Lyman-\alpha (DLA) and sub-damped Ly-\alpha\ absorption systems serendipitously discovered in the COS-Halos survey. We find that these systems show a large range in metallicities, with -1.10<[Z/H]<0.31, similar to the spread found at higher redshifts. These low-redshift systems on average have subsolar metallicities, but do show a rise in metallicity over cosmic time when compared to higher-redshift systems. We find the average sub-DLA metallicity is higher than the average DLA metallicity at all redshifts. Nitrogen is underabundant with respect to \alpha-group elements in all but perhaps one of the absorbers. In some cases, [N/\alpha] is significantly below the lowest nitrogen measurements in nearby galaxies. Systems for which depletion patterns can be studied show little, if any, depletion, which is characteristic of Milky Way halo-type gas. We also identify affiliated galaxies for 3 of the sub-DLAs using spectra obtained from Keck/LRIS. None of these sub-DLAs arise in the stellar disks of luminous galaxies; instead, these absorbers may exist in galaxy halos at impact parameters ranging from 38 to 92 kpc. Multiple galaxies are present near two of the sub-DLAs, and galaxy interactions may play a role in the dispersal of the gas. Many of these low-redshift absorbers exhibit simple kinematics, but one sub-DLA has a complicated mix of at least 13 components spread over 150 km/s. We find three galaxies near this sub-DLA, which also suggests that galaxy interactions roil the gas. This study reinforces the view that DLAs have a variety of origins, and low-redshift studies are crucial for understanding absorber-galaxy connections.

We present the results of a detailed analysis of the projected velocity dispersion of the globular cluster Palomar 14 performed using recent high-resolution spectroscopic data and extensive Monte Carlo simulations. The comparison between the data and a set of dynamical models (differing in fraction of binaries, degree of anisotropy, mass-to-light ratio M/L, cluster orbit and theory of gravity) shows that the observed velocity dispersion of this stellar system is well reproduced by Newtonian models with a fraction of binaries f_b<30% and a M/L compatible with the predictions of stellar evolution models. Instead, models computed with a large fraction of binaries systematically overestimate the cluster velocity dispersion. We also show that, across the parameter space sampled by our simulations, models based on the Modified Newtonian Dynamics theory can be reconciled with observations only assuming values of M/L lower than those predicted by stellar evolution models under standard assumptions.