Search

We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec- IFU spectroscopy of the young supernova remnant Cassiopeia A (Cas A). We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO emission is stronger at the outer layers than the Ar ejecta, which indicates the reformation of CO molecules behind the reverse shock. NIRSpec-IFU spectra (3 - 5.5 microns) were obtained toward two representative knots in the NE and S fields. Both regions are dominated by the bright fundamental rovibrational band of CO in the two R and P branches, with strong [Ar VI] and relatively weaker, variable strength ejecta lines of [Si IX], [Ca IV], [Ca V] and [Mg IV]. The NIRSpec-IFU data resolve individual ejecta knots and filaments spatially and in velocity space. The fundamental CO band in the JWST spectra reveals unique shapes of CO, showing a few tens of sinusoidal patterns of rovibrational lines with pseudo-continuum underneath, which is attributed to the high-velocity widths of CO lines. The CO also shows high J lines at different vibrational transitions. Our results with LTE modeling of CO emission indicate a temperature of 1080 K and provide unique insight into the correlations between dust, molecules, and highly ionized ejecta in supernovae, and have strong ramifications for modeling dust formation that is led by CO cooling in the early Universe., Comment: accepted for the ApJ letter (17 pages and 10 figures)

JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the equatorial ring. The emission of the remnant in the NIRCam 1-2.3 micron images is mostly due to line emission, which is mostly emitted in the ejecta and in the hot spots within the equatorial ring. In contrast, the NIRCam 3-5 micron images are dominated by continuum emission. In the ejecta, the continuum is due to dust, obscuring the centre of the ejecta. In contrast, in the ring and exterior to the ring, synchrotron emission contributes a substantial fraction to the continuum. Dust emission contributes to the continuum at outer spots and diffuse emission exterior to the ring, but little within the ring. This shows that dust cooling and destruction time scales are shorter than the synchrotron cooling time scale, and the time scale of hydrogen recombination in the ring is even longer than the synchrotron cooling time scale. With the advent of high sensitivity and high angular resolution images provided by JWST/NIRCam, our observations of SN 1987A demonstrate that NIRCam opens up a window to study particle-acceleration and shock physics in unprecedented details, probed by near-infrared synchrotron emission, building a precise picture of how a SN evolves., Comment: Accepted for publication in MNRAS; 18 pages

Upon its release the Gaia DR3 catalogue has led to tremendous progress in multiple fields of astronomy by providing the complete astrometric solution for nearly 1.5 billion sources. We analysed the photometric and astrometric results for Hydrogen-deficient Carbon (HdC), Extreme Helium (EHe), and DYPer type stars to identify any potential biases and thus to select stars suitable for kinematic and spatial distribution studies. We investigated the impact of photometric declines in R Coronae Borealis (RCB) stars on Gaia astrometry using information from the Gaia IPD process cross-matched with light curves. We have reached the conclusion that the astrometric fits for numerous RCB stars are not valid due to the Gaia PSF chromaticity effect in both shape and centroid. The astrometric results of all stars with a significant time-dependent colour variation should be similarly affected. RCB stars might thus be promising sources to correct this effect in future Gaia releases. Furthermore, after validating the Gaia astrometric results for 92 stars, we observed that the majority of HdC and EHe stars are distributed across the three old stellar structures, the thick disk, the bulge and the halo. However, we have also uncovered evidence indicating that some of them exhibit orbits characteristic of the thin disk. This is also particularly true for all DYPer type stars studied. Finally, we have produced a list of star memberships for each Galactic substructure, and provided a list of heliocentric radial velocities and associated errors for targets not observed by Gaia DR3. We are beginning to observe a relationship between kinematics, stellar population, and metallicity in RCB and EHe stars. That relation can be explained, within the double degenerate scenario, by the large range in the delay time distribution expected from population synthesis simulations, particularly through the HybCO merger channel., Comment: 22 pages, 20 figures, accepted for publication in A&A v2: RV section from accompanying paper ArXiv2309.10139 moved into that one +New section added on the comparison between our classification of HdC and EHe stars within a Galactic substructure and the one published lately by Philip Monai et al. (2023), and between both results on dynamical orbital parameters as determined with Galpy simulations

The Gaia DR3 release includes heliocentric radial velocity measurements and velocity variability indices for tens of millions of stars observed over 34 months.In this study, we utilise these indices to investigate the intrinsic radial velocity variations of Hydrogen-deficient Carbon (HdC) stars and Extreme Helium (EHe) stars across their large ranges of temperature and brightness. Taking advantage of the newly defined HdC temperature classes, we examine the evolution of the total velocity amplitude with effective temperature. Additionally, we analyse the variation in the dust production rate of R Coronae Borealis (RCB) stars with temperature using two different proxies for the photometric state of RCB stars: one from Gaia and another from the 2MASS survey. Our observations revealed a trend in the evolution of the maximum radial velocity amplitude across each HdC temperature class. Similarly, we also observed a correlation between stellar temperature and the dust production rate. Interestingly, we possibly observed for the first time some variations of the intrinsic radial velocity amplitude and the dust production rate with HdC temperature class. If confirmed, these variations would indicate that the helium shell-burning giant stage starts with strong atmospheric motions that decrease in strength, up to $\sim$6000 K, before picking up again as the HdC star atmosphere shrinks further in size and reaches warmer temperatures. Moreover, the dust formation rate appears to be much higher in colder RCB stars compared to warmer ones., Comment: 7 pages, 4 figures, accepted in A&A V2: version accepted, radial velocity section of V1 moved into the accompanying paper ArXiv2309.10148

We present a study of the dust associated with the core-collapse supernova SN~1995N. Infrared emission detected 14--15 years after the explosion was previously attributed to thermally echoing circumstellar material associated with the SN progenitor. We argue that this late-time emission is unlikely to be an echo, and is more plausibly explained by newly formed dust in the supernova ejecta, indirectly heated by the interaction between the ejecta and the CSM. Further evidence in support of this scenario comes from emission line profiles in spectra obtained 22 years after the explosion; these are asymmetric, showing greater attenuation on the red wing, consistent with absorption by dust within the expanding ejecta. The spectral energy distribution and emission line profiles at epochs later than $\sim$5000 days are both consistent with the presence of about 0.4~M$_\odot$ of amorphous carbon dust. The onset of dust formation is apparent in archival optical spectra, taken between 700 and 1700 days after the assumed explosion date. As this is considerably later than most other instances where the onset of dust formation has been detected, we argue that the explosion date must be later than previously assumed., Comment: 14 pages, 16 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

This effort is focused on examining the behavior of reinforcement learning systems in personalization environments and detailing the differences in policy entropy associated with the type of learning algorithm utilized. We demonstrate that Policy Optimization agents often possess low-entropy policies during training, which in practice results in agents prioritizing certain actions and avoiding others. Conversely, we also show that Q-Learning agents are far less susceptible to such behavior and generally maintain high-entropy policies throughout training, which is often preferable in real-world applications. We provide a wide range of numerical experiments as well as theoretical justification to show that these differences in entropy are due to the type of learning being employed.

Core-collapse supernovae (CCSNe) are potentially capable of producing large quantities of dust, with strong evidence that ejecta dust masses can grow significantly over extended periods of time. Red-blue asymmetries in the broad emission lines of CCSNe can be modelled using the Monte Carlo radiative transfer code DAMOCLES, to determine ejecta dust masses. To facilitate easier use of DAMOCLES, we present a Tkinter graphical user interface (GUI) running DAMOCLES. The GUI was tested by high school students as part of the Original Research By Young Twinkle Students (ORBYTS) programme, who used it to measure the dust masses formed at two epochs in two Type IIP CCSNe: SN 2012aw and iPTF14hls, demonstrating that a wide range of people can contribute significantly to scientific advancement. Bayesian methods were used to quantify uncertainties on our model parameters. From the presence of a red scattering wing in the day 1863 H$\alpha$ profile of SN 2012aw, we were able to constrain the dust composition to large (radius $>0.1 \mu$m) silicate grains, with a dust mass of $6.0^{+21.9}_{-3.6}\times10^{-4} M_\odot$. From the day 1158 H$\alpha$ profile of SN 2012aw, we found a dust mass of $3.0^{+14}_{-2.5}\times10^{-4}$ M$_\odot$. For iPTF14hls, we found a day 1170 dust mass of 8.1 $^{+81}_{-7.6}\times10^{-5}$ M$_{\odot}$ for a dust composition consisting of 50% amorphous carbon and 50% astronomical silicate. At 1000 days post explosion, SN 2012aw and iPTF14hls have formed less dust than SN 1987A, suggesting that SN 1987A could form larger dust masses than other Type IIP's., Comment: Accepted by MNRAS on 28/11/2022, 10 pages, 6 figures. Author accepted manuscript

Decades after their discovery, only four hydrogen-deficient carbon (HdC) stars were known to have no circumstellar dust shell. This is in complete contrast to the $\sim$130 known Galactic HdC stars that are notorious for being heavy dust producers, i.e. the R Coronae Borealis (RCB) stars. Together they form a rare class of supergiant stars that are thought to originate from the merger of CO/He white dwarf binary systems, otherwise known as the double-degenerate scenario. We searched for new dustless HdC (dLHdC) stars. We used primarily the 2MASS and GAIA eDR3 catalogues to select candidates that were followed-up spectroscopically. We discovered 27 new dLHdC stars, one new RCB star and two new EHe stars. Surprisingly, 20 of the new dLHdC stars share a characteristic of the known dLHdC star HD 148839, having lower atmospheric hydrogen deficiencies. The uncovered population of dLHdC stars exhibit a Bulge-like distribution, like the RCB stars, but show multiple differences from those that indicate they are a different population of HdC stars following its own evolutionary sequence with a fainter luminosity and also a narrow range of effective temperature. We found indication of a current low dust production activity for four of the new dLHdC stars which could be typical RCB stars passing through a transition time. We have evidence for the first time of a large range of absolute magnitudes in the overall population of HdC stars, spanning over 3 mag. In the favoured formation framework, this is explained by a large range in the initial total WD binary mass which leads to a series of evolutionary sequences with distinct maximum brightness and initial temperature. The cold Galactic RCB stars are also noticeably fainter than the Magellanic ones, possibly due to a difference in metallicity resulting in different WD mass ratio. In our Galaxy, there could be as many dLHdC stars as RCB stars., Comment: 22 pages, 18 figures, accepted in A&A. V3: minor corrections on Table 2 and 3 related to HD stars names and (V-I)0 colours (Fig.10 changed accordingly). V4: changed the colour to light blue for the two new EHe stars indicated in Fig.10, page 13

Luminous Red Variables (LRVs) are most likely eruptions that are the outcome of stellar mergers. V838 Mon is one of the best-studied members of this class, representing an archetype for stellar mergers resulting from B-type stars. As result of the merger event, nova-like eruptions occur driving mass-loss from the system. As the gas cools considerable circumstellar dust is formed. V838 Mon erupted in 2002 and is undergoing very dynamic changes in its dust composition, geometry, and infrared luminosity providing a real-time laboratory to validate mineralogical condensation sequences in stellar mergers and evolutionary scenarios. We discuss recent NASA Stratospheric Observatory for Infrared Astronomy SOFIA 5 to 38 micron observations combined with archival NASA Spitzer spectra that document the temporal evolution of the freshly formed (within the last 20 yrs) circumstellar material in the environs of V838 Mon. Changes in the 10 micron spectral region are strong evidence that we are witnessing a classical dust condensation sequence expected to occur in oxygen-rich environments where alumina formation is followed by that of silicates at the temperature cools., Comment: 12 Pages, 6 figures, 2 tables, To be published in the Astronomical Journal

In this work, we propose a novel adaptive stochastic gradient-free (ASGF) approach for solving high-dimensional nonconvex optimization problems based on function evaluations. We employ a directional Gaussian smoothing of the target function that generates a surrogate of the gradient and assists in avoiding bad local optima by utilizing nonlocal information of the loss landscape. Applying a deterministic quadrature scheme results in a massively scalable technique that is sample-efficient and achieves spectral accuracy. At each step we randomly generate the search directions while primarily following the surrogate of the smoothed gradient. This enables exploitation of the gradient direction while maintaining sufficient space exploration, and accelerates convergence towards the global extrema. In addition, we make use of a local approximation of the Lipschitz constant in order to adaptively adjust the values of all hyperparameters, thus removing the careful fine-tuning of current algorithms that is often necessary to be successful when applied to a large class of learning tasks. As such, the ASGF strategy offers significant improvements when solving high-dimensional nonconvex optimization problems when compared to other gradient-free methods (including the so called "evolutionary strategies") as well as iterative approaches that rely on the gradient information of the objective function. We illustrate the improved performance of this method by providing several comparative numerical studies on benchmark global optimization problems and reinforcement learning tasks.

We first propose a novel criterion that guarantees that an $s$-sparse signal is the local minimizer of the $\ell_1/\ell_2$ objective; our criterion is interpretable and useful in practice. We also give the first uniform recovery condition using a geometric characterization of the null space of the measurement matrix, and show that this condition is easily satisfied for a class of random matrices. We also present analysis on the robustness of the procedure when noise pollutes data. Numerical experiments are provided that compare $\ell_1/\ell_2$ with some other popular non-convex methods in compressed sensing. Finally, we propose a novel initialization approach to accelerate the numerical optimization procedure. We call this initialization approach \emph{support selection}, and we demonstrate that it empirically improves the performance of existing $\ell_1/\ell_2$ algorithms., Comment: 25 pages, 7 figures

The luminous Type IIn SN 2010jl shows strong signs of interaction between the SN ejecta and dense circumstellar material. Dust may be present in the unshocked ejecta, the cool, dense shell between the shocks in the interaction region, or in the circumstellar medium. We present and model new optical and infrared photometry and spectroscopy of SN 2010jl from 82 to 1367 days since explosion. We evaluate the photometric and spectroscopic evolution using the radiative transfer codes MOCASSIN and DAMOCLES, respectively. We propose an interaction scenario and investigate the resulting dust formation scenarios and dust masses. We find that SN 2010jl has been continuously forming dust based on the evolution of its infrared emission and optical spectra. There is evidence for pre-existing dust in the circumstellar medium as well as new dust formation in the cool, dense shell and/or ejecta. We estimate that 0.005--0.01M$_{\odot}$ of predominantly carbon dust grains has formed in SN 2010jl by ~1400d post-outburst., Comment: Accepted by ApJ. Author accepted manuscript. Accepted on 26/03/2020. Deposited on 03/04/20. 21 pages, 10 figures

We present an analysis of the evolution of circumstellar dust and molecules in the environment of the very late thermal pulse object V4334 Sgr (Sakurai's Object) over a $\sim20$-year period, drawing on ground-, airborne- and space-based infrared photometry and spectroscopy. The dust emission, which started in 1997, resembles a blackbody that cooled from $\sim1200$K in 1998 August to $\sim180$K in 2016 July. The dust mass, assuming amorphous carbon, was $\sim5\times10^{-10}$M$_\odot$ in 1998 August, and we estimate that the total dust mass was $\sim2\times10^{-5}$M$_\odot$ by $\sim2016$. The appearance of a near infrared excess in 2008 suggests a new episode of (or renewed) mass loss began then. We infer lower limits on the bolometric luminosity of the embedded star from that of the dust shell, which rose to $\sim16000$L$_\odot$ before declining to $\sim3000$L$_\odot$. There is evidence for weak 6-7$\mu$m absorption, which we attribute to hydrogenated amorphous carbon formed in material ejected by Sakurai's Object during a mass ejection phase that preceded the 1997 event. We detect small hydrocarbon and other molecules in the spectra, and trace the column densities in hydrogen cyanide (HCN) and acetylene (C$_2$H$_2$). We use the former to determine the $^{12}$C/$^{13}$C ratio to be $6.4\pm0.7$, 14 times smaller than the Solar System value.

Developing accurate, efficient, and robust closure models is essential in the construction of reduced order models (ROMs) for realistic nonlinear systems, which generally require drastic ROM mode truncations. We propose a deep residual neural network (ResNet) closure learning framework for ROMs of nonlinear systems. The novel ResNet-ROM framework consists of two steps: (i) In the first step, we use ROM projection to filter the given nonlinear PDE and construct a filtered ROM. This filtered ROM is low-dimensional, but is not closed (because of the PDE nonlinearity). (ii) In the second step, we use ResNet to close the filtered ROM, i.e., to model the interaction between the resolved and unresolved ROM modes. We emphasize that in the new ResNet-ROM framework, data is used only to complement classical physical modeling (i.e., only in the closure modeling component), not to completely replace it. We also note that the new ResNet-ROM is built on general ideas of spatial filtering and deep learning and is independent of (restrictive) phenomenological arguments, e.g., of eddy viscosity type. The numerical experiments for the 1D Burgers equation show that the ResNet-ROM is significantly more accurate than the standard projection ROM. The new ResNet-ROM is also more accurate and significantly more efficient than other modern ROM closure models.

We present a nonlocal variational image completion technique which admits simultaneous inpainting of multiple structures and textures in a unified framework. The recovery of geometric structures is achieved by using general convolution operators as a measure of behavior within an image. These are combined with a nonlocal exemplar-based approach to exploit the self-similarity of an image in the selected feature domains and to ensure the inpainting of textures. We also introduce an anisotropic patch distance metric to allow for better control of the feature selection within an image and present a nonlocal energy functional based on this metric. Finally, we derive an optimization algorithm for the proposed variational model and examine its validity experimentally with various test images.

In this effort, we propose a convex optimization approach based on weighted $\ell_1$-regularization for reconstructing objects of interest, such as signals or images, that are sparse or compressible in a wavelet basis. We recover the wavelet coefficients associated to the functional representation of the object of interest by solving our proposed optimization problem. We give a specific choice of weights and show numerically that the chosen weights admit efficient recovery of objects of interest from either a set of sub-samples or a noisy version. Our method not only exploits sparsity but also helps promote a particular kind of structured sparsity often exhibited by many signals and images. Furthermore, we illustrate the effectiveness of the proposed convex optimization problem by providing numerical examples using both orthonormal wavelets and a frame of wavelets. We also provide an adaptive choice of weights which is a modification of the iteratively reweighted $\ell_1$-minimization method., Comment: 16 pages and 20 figures

In this effort we propose a novel approach for reconstructing multivariate functions from training data, by identifying both a suitable network architecture and an initialization using polynomial-based approximations. Training deep neural networks using gradient descent can be interpreted as moving the set of network parameters along the loss landscape in order to minimize the loss functional. The initialization of parameters is important for iterative training methods based on descent. Our procedure produces a network whose initial state is a polynomial representation of the training data. The major advantage of this technique is from this initialized state the network may be improved using standard training procedures. Since the network already approximates the data, training is more likely to produce a set of parameters associated with a desirable local minimum. We provide the details of the theory necessary for constructing such networks and also consider several numerical examples that reveal our approach ultimately produces networks which can be effectively trained from our initialized state to achieve an improved approximation for a large class of target functions., Comment: 11 pages, 6 figures, submitted to NeurIPS 2019, corrected several typos and included new examples

We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses., Comment: 25 pages; 28 figures (Accepted to MNRAS)

In order to understand the contribution of core-collapse supernovae to the dust budget of the early universe, it is important to understand not only the mass of dust that can form in core-collapse supernovae but also the location and rate of dust formation. SN 2005ip is of particular interest since dust has been inferred to have formed in both the ejecta and the post-shock region behind the radiative reverse shock. We have collated eight optical archival spectra that span the lifetime of SN 2005ip and we additionally present a new X-shooter optical-near-IR spectrum of SN 2005ip at 4075d post-discovery. Using the Monte Carlo line transfer code DAMOCLES, we have modelled the blueshifted broad and intermediate width H$\alpha$, H$\beta$ and He I lines from 48d to 4075d post-discovery using an ejecta dust model. We find that dust in the ejecta can account for the asymmetries observed in the broad and intermediate width H$\alpha$, H$\beta$ and He I line profiles at all epochs and that it is not necessary to invoke post-shock dust formation to explain the blueshifting observed in the intermediate width post-shock lines. Using a Bayesian approach, we have determined the evolution of the ejecta dust mass in SN 2005ip over 10 years presuming an ejecta dust model, with an increasing dust mass from ~10$^{-8}$ M$_{\odot}$ at 48d to a current dust mass of $\sim$0.1 M$_{\odot}$., Comment: Accepted by MNRAS, 17 pages, 11 figures. Author accepted manuscript. Accepted on 04/03/19. Deposited on 07/03/19

In this effort we propose a data-driven learning framework for reduced order modeling of fluid dynamics. Designing accurate and efficient reduced order models for nonlinear fluid dynamic problems is challenging for many practical engineering applications. Classical projection-based model reduction methods generate reduced systems by projecting full-order differential operators into low-dimensional subspaces. However, these techniques usually lead to severe instabilities in the presence of highly nonlinear dynamics, which dramatically deteriorates the accuracy of the reduced-order models. In contrast, our new framework exploits linear multistep networks, based on implicit Adams-Moulton schemes, to construct the reduced system. The advantage is that the method optimally approximates the full order model in the low-dimensional space with a given supervised learning task. Moreover, our approach is non-intrusive, such that it can be applied to other complex nonlinear dynamical systems with sophisticated legacy codes. We demonstrate the performance of our method through the numerical simulation of a two-dimensional flow past a circular cylinder with Reynolds number Re = 100. The results reveal that the new data-driven model is significantly more accurate than standard projection-based approaches.

It is increasingly suspected that the rare R Coronae Borealis (RCB) stars - hydrogen-deficient and carbon-rich supergiant stars - are the products of mergers of CO/He white-dwarf binary systems in the intermediate mass regime ($0.615000$ K). Forty of these belong to the Milky Way and five are located in the Magellanic Clouds. We also confirmed that the long lasting candidate KDM 5651 is indeed a new Magellanic RCB star, increasing the total number of Magellanic Cloud RCB stars to 30. We have increased by $\sim$50\% the total number of RCB stars known, now reaching 147. We also include a list of 14 strong RCB candidates, most certainly observed during a dust obscuration phase. From the detection efficiency and success rate so far, we estimate that there should be no more than 500 RCB stars/HdC stars in the Milky Way., Comment: 34 pages, 20 figures

R Coronae Borealis stars (RCBs) are the best candidates to be the products of white dwarf mergers in the intermediate mass range ($0.4

The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21 Treasury program on the Hubble Space Telescope. It studied 50 nearby star-forming galaxies in five bands from the near UV to the I-band, combining new Wide Field Camera 3 observations with archival Advanced Camera for Surveys data. LEGUS was designed to investigate how star formation occurs and develops on both small and large scales, and how it relates to the galactic environments. In this paper we present the photometric catalogs for all the apparently single stars identified in the 50 LEGUS galaxies. Photometric catalogs and mosaicked images for all filters are available for download. We present optical and near UV color-magnitude diagrams for all the galaxies. For each galaxy we derived the distance from the tip of the red giant branch. We then used the NUV color-magnitude diagrams to identify stars more massive than 14 Mo, and compared their number with the number of massive stars expected from the GALEX FUV luminosity. Our analysis shows that the fraction of massive stars forming in star clusters and stellar associations is about constant with the star formation rate. This lack of a relation suggests that the time scale for evaporation of unbound structures is comparable or longer than 10 Myr. At low star formation rates this translates to an excess of mass in clustered environments as compared to model predictions of cluster evolution, suggesting that a significant fraction of stars form in unbound systems., Comment: 71 pages, 15 figures, 4 tables, Accpeted for publication on ApJS

JWST was conceived and built to answer one of the most fundamental questions that humans can address empirically: "How did the Universe make its first stars?". Our First Lights At REionization (FLARE) project transforms the quest for the epoch of reionization from the static to the time domain. It targets the complementary question: "What happened to those first stars?". It will be answered by observations of the most luminous events: supernovae and accretion on to black holes formed by direct collapse from the primordial gas clouds. These transients provide direct constraints on star-formation rates (SFRs) and the truly initial Initial Mass Function (IMF), and they may identify possible stellar seeds of supermassive black holes (SMBHs). Furthermore, our knowledge of the physics of these events at ultra-low metallicity will be much expanded. JWST's unique capabilities will detect these most luminous and earliest cosmic messengers easily in fairly shallow observations. However, these events are very rare at the dawn of cosmic structure formation and so require large area coverage. Time domain astronomy can be advanced to an unprecedented depth by means of a shallow field of JWST reaching 27 mag AB in 2 and 4.4 microns over a field as large as 0.1 square degree visited multiple times each year. Such a survey may set strong constraints or detect massive Pop III SNe at redshifts beyond 10, pinpointing the redshift of the first stars, or at least their death. Based on our current knowledge of superluminous supernovae (SLSNe), such a survey will find one or more SLSNe at redshifts above 6 in five years and possibly several direct collapse black holes. Although JWST is not designed as a wide field survey telescope, we show that such a wide field survey is possible with JWST and is critical in addressing several of its key scientific goals.

The Infrared Spectrograph (IRS) on the {\em Spitzer Space Telescope} observed nearly 800 point sources in the Large Magellanic Cloud (LMC), taking over 1,000 spectra. 197 of these targets were observed as part of the Sage-Spec Spitzer Legacy program; the remainder are from a variety of different calibration, guaranteed time and open time projects. We classify these point sources into types according to their infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, using a decision-tree classification method. We then refine the classification using supplementary information from the astrophysical literature. We find that our IRS sample is comprised substantially of YSO and H\,{\sc ii} regions, post-Main Sequence low-mass stars: (post-)AGB stars and planetary nebulae and massive stars including several rare evolutionary types. Two supernova remnants, a nova and several background galaxies were also observed. We use these classifications to improve our understanding of the stellar populations in the Large Magellanic Cloud, study the composition and characteristics of dust species in a variety of LMC objects, and to verify the photometric classification methods used by mid-IR surveys. We discover that some widely-used catalogues of objects contain considerable contamination and others are missing sources in our sample., Comment: Accepted for publication in MNRAS

We report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628., Comment: accepted for publication in ApJ; 27 pages; complete stellar cluster catalogues will be released in June 2017 at https://legus.stsci.edu

In recent years, the use of sparse recovery techniques in the approximation of high-dimensional functions has garnered increasing interest. In this work we present a survey of recent progress in this emerging topic. Our main focus is on the computation of polynomial approximations of high-dimensional functions on $d$-dimensional hypercubes. We show that smooth, multivariate functions possess expansions in orthogonal polynomial bases that are not only approximately sparse, but possess a particular type of structured sparsity defined by so-called lower sets. This structure can be exploited via the use of weighted $\ell^1$ minimization techniques, and, as we demonstrate, doing so leads to sample complexity estimates that are at most logarithmically dependent on the dimension $d$. Hence the curse of dimensionality - the bane of high-dimensional approximation - is mitigated to a significant extent. We also discuss several practical issues, including unknown noise (due to truncation or numerical error), and highlight a number of open problems and challenges.

This work focuses on weighted Lagrange interpolation on an unbounded domain, and analyzes the Lebesgue constant for a sequence of weighted Leja points. The standard Leja points are a nested sequence of points defined on a compact subset of the real line, and can be extended to unbounded domains with the introduction of a weight function $w:\mathbb{R}\rightarrow [0,1]$. Due to a simple recursive formulation in one dimension, such abscissas provide a foundation for high-dimensional approximation methods such as sparse grid collocation, deterministic least squares, and compressed sensing. Just as in the unweighted case of interpolation on a compact domain, we use results from potential theory to prove that the Lebesgue constant for the Leja points grows subexponentially with the number of interpolation nodes.

This work proposes and analyzes a compressed sensing approach to polynomial approximation of complex-valued functions in high dimensions. Of particular interest is the setting where the target function is smooth, characterized by a rapidly decaying orthonormal expansion, whose most important terms are captured by a lower (or downward closed) set. By exploiting this fact, we present an innovative weighted $\ell_1$ minimization procedure with a precise choice of weights, and a new iterative hard thresholding method, for imposing the downward closed preference. Theoretical results reveal that our computational approaches possess a provably reduced sample complexity compared to existing compressed sensing techniques presented in the literature. In addition, the recovery of the corresponding best approximation using these methods is established through an improved bound for the restricted isometry property. Our analysis represents an extension of the approach for Hadamard matrices in [5] to the general case of continuous bounded orthonormal systems, quantifies the dependence of sample complexity on the successful recovery probability, and provides an estimate on the number of measurements with explicit constants. Numerical examples are provided to support the theoretical results and demonstrate the computational efficiency of the novel weighted $\ell_1$ minimization strategy., Comment: 33 pages, 3 figures