Synthetic Spectral Energy Distributions Research Articles

The Universe SPHEREx Will See: Empirically Based Galaxy Simulations and Redshift Predictions

We simulate galaxy properties and redshift estimation for SPHEREx, the next NASA Medium Class Explorer. To make robust models of the galaxy population and test the spectrophotometric redshift performance for SPHEREx, we develop a set of synthetic spectral energy distributions based on detailed fits to COSMOS2020 photometry spanning 0.1–8 μm. Given that SPHEREx obtains low-resolution spectra, emission lines will be important for some fraction of galaxies. Here, we expand on previous work, using better photometry and photometric redshifts from COSMOS2020 and tight empirical relations to predict robust emission-line strengths and ratios. A second galaxy catalog derived from the GAMA survey is generated to ensure the bright (m AB < 18 in the i band) sample is representative over larger areas. Using template fitting to estimate photometric continuum redshifts, we forecast the recovery of 19 million galaxies over 30,000 deg2 with redshifts better than σ z < 0.003(1 + z), 445 million with σ z < 0.1(1 + z), and 810 million with σ z < 0.2(1 + z). We also find through idealized tests that emission-line information from spectrally dithered flux measurements can yield redshifts with accuracy beyond that implied by the naive SPHEREx channel resolution, motivating the development of a hybrid continuum–line redshift estimation approach.

Variation of the stellar initial mass function in semi-analytical models

Context. In our previous work, we derived the CR-IGIMF, which is a new scenario for a variable stellar initial mass function (IMF) that combines numerical results on the role played by cosmic rays (CR) in setting the thermal state of star-forming gas with the analytical approach of the integrated galaxy-wide IMF (IGIMF). Aims. In this work, we study the implications of this scenario for the properties of local early-type galaxies (ETG) as inferred from dynamical, photometric, and spectroscopic studies. Methods. We implemented a library of CR-IGIMF shapes in the framework of the galaxy evolution and assembly (GAEA) model. GAEA provides predictions for the physical and photometric properties of model galaxies and for their chemical composition. Our realization includes a self-consistent derivation of the synthetic spectral energy distribution for each model galaxy, which allows a direct derivation of the mass fraction in the mean IMF of low-mass stars (i.e., the dwarf-to-giant ratio, fdg) and a comparison with IMF-sensitive spectral features. Results. The predictions of the GAEA model implementing the CR-IGIMF confirm our previous findings: It correctly reproduces both the observed excess of z ∼ 0 dynamical mass (mass-to-light ratio) with respect to spectroscopic (photometric) estimates assuming a universal MW-like IMF, and the observed increase in [α/Fe] ratios with stellar mass in spheroidal galaxies. Moreover, this realization reproduces the increasing trends of fdg and IMF-sensitive line strengths with velocity dispersion, although the predicted relations are significantly shallower than the observed ones. Conclusions. Our results show that the CR-IGIMF is a promising scenario that reproduces at the same time dynamical, photometric, and spectroscopic indications of a varying IMF in local ETGs. The shallow relations found for spectral indices suggest that either a stronger variability as a function of galaxy properties or additional dependences (e.g., as a function of star forming gas metallicity) might be required to match the strength of the observed trends.

Probing the inner Galactic halo with blue horizontal-branch stars

Context. Stars that are found on the blue horizontal-branch (BHB) evolved from low-mass stars that have completed their core hydrogen-burning main sequence (MS) stage and undergone the helium flash at the end of their red giant phase. Hence, they are very old objects that can be used as markers in studying galactic structure and formation history. The fact that their luminosity is virtually constant at all effective temperatures also makes them good standard candles. Aims. We have compiled a catalogue of BHB stars with stellar parameters calculated from spectral energy distributions (SEDs) constructed using data from multiple large-scale photometric surveys. In addition, we update our previous Gaia Early Data Release 3 (EDR3) catalogue of BHB stars with parallax errors less than 20% by using the SED results to define the selection criteria. The purpose of these catalogues is to create a set of BHB star candidates with reliable stellar parameters. In addition, they provide a more complete full-sky catalogue with candidate objects found along the whole BHB from where RR-Lyrae are found on the instability strip to the extreme horizontal-branch (EHB). Methods. We selected a large dataset of Gaia Data Release 3 (DR3) objects based only on their position on the colour-magnitude diagram (CMD), along with the tangential velocity and parallax errors. The SEDs were then used to evaluate contamination levels in the dataset and derive optimised data quality acceptance constraints. This allowed us to extend the Gaia DR3 colour and absolute magnitude criteria further towards the EHB. The level of contamination found using SED analysis was confirmed by acquiring spectra using the Ondrejov Echelle spectrograph, attached to the Perek 2m telescope at the Astronomical Institute of the Czech Academy of Sciences. Results. We present a catalogue of 9172 Galactic halo BHB candidate stars with atmospheric and stellar parameters calculated from synthetic SEDs. We also present an extended Gaia DR3-based catalogue of 22 335 BHB candidate stars with a wider range of effective temperatures and Gaia DR3 parallax errors of less than 20%. This represents an increase of 33% compared to the our 2021 catalogue, with a contamination level of 10%.

The improved component masses and parallaxes for the two close binary stars: HD 80671 and HD 97038

We present the fundamental stellar parameters, including the individual component masses, as well as the orbital parameters, and dynamical parallaxes of the two close binary stars; HD80671, and HD97038. The stellar parameters are spectrophotometrically estimated via Al-Wardat’s method for analyzing binary and multiple stellar systems, which is having a combination of the spectroscopic analysis and the photometric analysis to build the combined and individual synthetic spectral energy distributions of the individual components of the systems and so to estimate their fundamental parameters, metallicities, and ages. It employs Kurucz’s model atmospheres of single stars, while the orbital parameters are estimated using Tokovinin’s method. The individual spectrophotometric component masses are inferred with good accuracy, and found to be MSphA = 1.47±0.10M⊙ and MSphB=1.29±0.09M⊙ with an age of 1.0±0.09 Gyr for HD80671, and MSphA=1.17±0.09M⊙ and MSphB=1.12±0.08M⊙ with an age of 3.981±0.35 Gyr for HD97038. The improved dynamical parallaxes are given as πdyn=28.305±0.45 mas for HD80671, and πdyn=16.26±0.30 mas for HD97038. The evolutionary status of two binaries is discussed depending on the positions of the components on the isochrones and evolutionary tracks.

Beware the recent past: a bias in spectral energy distribution modelling due to bursty star formation

ABSTRACT We investigate how the recovery of galaxy star formation rates (SFRs) using energy-balance spectral energy distribution (SED) fitting codes depends on their recent star formation histories (SFHs). We use the magphys and prospector codes to fit 6706 synthetic SEDs of simulated massive galaxies at 1 &amp;lt; z &amp;lt; 8 from the Feedback in Realistic Environments project. We identify a previously unknown systematic error in the magphys results due to bursty star formation: the derived SFRs can differ from the truth by as much as 1 dex, at large statistical significance (&amp;gt;5σ), depending on the details of their recent SFH. SFRs inferred using prospector with non-parametric SFHs do not exhibit this trend. We show that using parametric SFHs (pSFHs) causes SFR uncertainties to be underestimated by a factor of up to 5×. Although this undoubtedly contributes to the significance of the systematic, it cannot explain the largest biases in the SFRs of the starbursting galaxies, which could be caused by details of the stochastic prior sampling or the burst implementation in the magphys libraries. We advise against using pSFHs and urge careful consideration of starbursts when SED modelling galaxies where the SFR may have changed significantly over the last ∼100 Myr, such as recently quenched galaxies, or those experiencing a burst. This concern is especially relevant, e.g. when fitting JWST observations of very high redshift galaxies.

An Updated Modular Set of Synthetic Spectral Energy Distributions for Young Stellar Objects

Measured properties of young stellar objects (YSOs) are key tools for research into pre-main-sequence stellar evolution. YSO properties are commonly measured by comparing observed radiation to existing grids of template YSO spectral energy distributions (SEDs) calculated by radiative transfer. These grids are often sampled and constructed using simple models of mass assembly/accretion over time. However, because we do not yet have a complete theory of star formation, the choice of model sets the tracked parameters and range of allowed values. By construction, then, the assumed model limits the measurements that can be made using the grid. Radiative transfer models not constrained by specific accretion histories would enable assessment of a wider range of theories. We present an updated version of the Robitaille set (2017) of YSO SEDs, a collection of models with no assumed evolutionary theory. We outline our newly calculated properties: envelope mass, weighted-average dust temperature, disk stability, and circumstellar A V. We also convolve the SEDs with new filters, including JWST, and provide users the ability to perform additional convolutions. We find a correlation between the average temperature and millimeter-wavelength brightness of optically thin dust in our models and discuss its ramifications for mass measurements of pre- and protostellar cores. We also compare the positions of YSOs of different observational classes and evolutionary stages in IR color space and use our models to quantify the extent to which class and stage may be confused due to observational effects. Our updated models are released to the public.

Emulating radiative transfer with artificial neural networks

ABSTRACT Forward-modeling observables from galaxy simulations enables direct comparisons between theory and observations. To generate synthetic spectral energy distributions (SEDs) that include dust absorption, re-emission, and scattering, Monte Carlo radiative transfer is often used in post-processing on a galaxy-by-galaxy basis. However, this is computationally expensive, especially if one wants to make predictions for suites of many cosmological simulations. To alleviate this computational burden, we have developed a radiative transfer emulator using an artificial neural network (ANN), ANNgelina, that can reliably predict SEDs of simulated galaxies using a small number of integrated properties of the simulated galaxies: star formation rate, stellar and dust masses, and mass-weighted metallicities of all star particles and of only star particles with age &amp;lt;10 Myr. Here, we present the methodology and quantify the accuracy of the predictions. We train the ANN on SEDs computed for galaxies from the IllustrisTNG project’s TNG50 cosmological magnetohydrodynamical simulation. ANNgelina is able to predict the SEDs of TNG50 galaxies in the ultraviolet (UV) to millimetre regime with a typical median absolute error of ∼7 per cent. The prediction error is the greatest in the UV, possibly due to the viewing-angle dependence being greatest in this wavelength regime. Our results demonstrate that our ANN-based emulator is a promising computationally inexpensive alternative for forward-modeling galaxy SEDs from cosmological simulations.

Energy balance SED modelling can be effective at high redshifts regardless of UV-FIR offsets

ABSTRACT Recent works have suggested that energy balance spectral energy distribution (SED) fitting codes may be of limited use for studying high-redshift galaxies for which the observed ultraviolet and far-infrared emission are offset (spatially ‘decoupled’). It has been proposed that such offsets could lead energy balance codes to miscalculate the overall energetics, preventing them from recovering such galaxies’ true properties. In this work, we test how well the SED fitting code magphys can recover the stellar mass, star formation rate (SFR), specific SFR, dust mass, and luminosity by fitting 6706 synthetic SEDs generated from four zoom-in simulations of dusty, high-redshift galaxies from the FIRE project via dust continuum radiative transfer. Comparing our panchromatic results (using wavelengths 0.4–500 μm, and spanning 1 &amp;lt; z &amp;lt; 8) with fits based on either the starlight ($\lambda _\mathrm{eff} \le 2.2\, \mu$m) or dust ($\ge 100\, \mu$m) alone, we highlight the power of considering the full range of multiwavelength data alongside an energy balance criterion. Overall, we obtain acceptable fits for 83 per cent of the synthetic SEDs, though the success rate falls rapidly beyond z ≈ 4, in part due to the sparser sampling of the priors at earlier times since SFHs must be physically plausible (i.e. shorter than the age of the universe). We use the ground truth from the simulations to show that when the quality of fit is acceptable, the fidelity of magphys estimates is independent of the degree of UV/FIR offset, with performance very similar to that previously reported for local galaxies.

Synthetic red supergiant explosion model grid for systematic characterization of Type II supernovae

Abstract A new model grid containing 228016 synthetic red supergiant explosions (Type II supernovae) is introduced. Time evolution of spectral energy distributions from 1 to 50000 Å (100 frequency bins in a log scale) is computed at each time step up to 500 d after explosion in each model. We provide light curves for the filters of Vera C, Rubin Observatory’s Legacy Survey of Space and Time (LSST), the Zwicky Transient Facility, the Sloan Digital Sky Survey, and the Neil Gehrels Swift Observatory, but light curves for any photometric filters can be constructed by convolving any filter response functions to the synthetic spectral energy distributions. We also provide bolometric light curves and photosphere information such as photospheric velocity evolution. The parameter space covered by the model grid is five progenitor masses (10, 12, 14, 16, and 18 M$_{\odot}$ at the zero-age main sequence, solar metallicity), ten explosion energies (0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, and 5.0 × 1051 erg), nine 56Ni masses (0.001, 0.01, 0.02, 0.04, 0.06, 0.08, 0.1, 0.2, and 0.3 M$_{\odot}$), nine mass-loss rates (10−5.0, 10−4.5, 10−4.0, 10−3.5, 10−3.0, 10−2.5, 10−2.0, 10−1.5, and 10−1.0 M$_{\odot}$ yr−1 with a wind velocity of 10 km s−1), six circumstellar matter radii (1, 2, 4, 6, 8, and 10 × 1014 cm), and ten circumstellar structures (β = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, and 5.0). 56Ni is assumed to be uniformly mixed up to the half-mass of a hydrogen-rich envelope. This model grid can be a base for rapid characterizations of Type II supernovae with sparse photometric sampling expected in LSST through a Bayesian approach, for example. The model grid is available at https://doi.org/10.5061/dryad.pnvx0k6sj.

Is the star-formation rate in z ∼ 6 quasars overestimated?

ABSTRACT The large total infrared (TIR) luminosities (LTIR; ≳1012 L⊙) observed in z ∼ 6 quasars are generally converted into high star-formation rates (SFRs; $\gtrsim\!{10}^2~{\rm M}_{\odot }\, {\rm yr}^{-1}$) of their host galaxies. However, these estimates rely on the assumption that dust heating is dominated by stellar radiation, neglecting the contribution from the central active galactic nucleus (AGN). We test the validity of this assumption by combining cosmological hydrodynamic simulations with radiative transfer calculations. We find that, when AGN radiation is included in the simulations, the mass (luminosity)-weighted dust temperature in the host galaxies increases from T ≈ 50 K (T ≈ 70 K) to T ≈ 80 K (T ≈ 200 K), suggesting that AGN effectively heats the bulk of dust in the host galaxy. We compute the AGN-host galaxy SFR from the synthetic spectral energy distribution by using standard SFR – LTIR relations, and compare the results with the ‘true’ values in the simulations. We find that the SFR is overestimated by a factor of ≈3 (≳10) for AGN bolometric luminosities of Lbol ≈ 1012 L⊙ (≳1013 L⊙), implying that the SFRs of z ∼ 6 quasars can be overestimated by over an order of magnitude.

Machine learning synthetic spectra for probabilistic redshift estimation: SYTH-Z

ABSTRACT Photometric redshift estimation algorithms are often based on representative data from observational campaigns. Data-driven methods of this type are subject to a number of potential deficiencies, such as sample bias and incompleteness. Motivated by these considerations, we propose using physically motivated synthetic spectral energy distributions in redshift estimation. In addition, the synthetic data would have to span a domain in colour-redshift space concordant with that of the targeted observational surveys. With a matched distribution and realistically modelled synthetic data in hand, a suitable regression algorithm can be appropriately trained; we use a mixture density network for this purpose. We also perform a zero-point recalibration to reduce the systematic differences between noise-free synthetic data and the (unavoidably) noisy observational data sets. This new redshift estimation framework, syth-z, demonstrates superior accuracy over a wide range of redshifts compared to baseline models trained on observational data alone. Approaches using realistic synthetic data sets can therefore greatly mitigate the reliance on expensive spectroscopic follow-up for the next generation of photometric surveys.

Physical and geometrical parameters of CVBS XIV: the two nearby systems HIP 19206 and HIP 84425

Data release 2 (DR2) from the Gaia mission was of great help in precise determination of fundamental parameters of Close Visual Binary and Multiple Systems (CVBMSs), especially masses of their components, which are crucial parameters in understanding formation and evolution of stars and galaxies. This article presents the complete set of fundamental parameters for two nearby close visual binary systems (CVBSs), which are HIP 19206 and HIP 84425. We utilised a combination of two methods; the first one is Tokovinin’s dynamical method to solve the orbit of the system and to estimate orbital elements and the dynamical mass sum, and the second one is Al-Wardat’s method for analysing CVBMSs to estimate the physical parameters of the individual components. The latest method employs grids of Kurucz line-blanketed plane parallel model atmospheres to build synthetic Spectral Energy Distributions (SEDs) of the individual components. Trigonometric parallax measurements given by Gaia DR2 and Hipparcos catalogues are used to analyse the two systems. The difference in these measurements yielded slight discrepancies in the fundamental parameters of the individual components, especially masses. So, a new dynamical parallax is suggested in this work based on the most convenient mass sum as given by each of the two methods. The new dynamical parallax for the system HIP 19205 of 22.97 ± 0.95 mas coincides well with the trigonometric one reported recently (in December 2020) by Gaia EDR3 of 22.3689 ± 0.4056 mas. The positions of the components of the two systems on the evolutionary tracks and isochrones are plotted, which suggest that all components are solar-type main sequence stars. Their most probable formation and evolution scenarios are also discussed.

The dust attenuation law in z ∼ 6 quasars

ABSTRACT We investigate the attenuation law in z ∼ 6 quasars by combining cosmological zoom-in hydrodynamical simulations of quasar host galaxies, with multifrequency radiative transfer calculations. We consider several dust models differing in terms of grain-size distributions, dust mass, and chemical composition, and compare the resulting synthetic spectral energy distributions with data from bright, early quasars. We show that only dust models with grain-size distributions in which small grains ($a\lesssim 0.1~\mu {\rm m}$, corresponding to $\approx 60{{\ \rm per\ cent}}$ of the total dust mass) are selectively removed from the dusty medium provide a good fit to the data. Removal can occur if small grains are efficiently destroyed in quasar environments and/or early dust production preferentially results in large grains. Attenuation curves for these models are close to flat, and consistent with recent data; they correspond to an effective dust-to-metal ratio fd ≃ 0.38, i.e. close to the Milky Way value.

Dwarf Satellites of High-z Lyman Break Galaxies: A Free Lunch for JWST

Abstract We show that the James Webb Space Telescope (JWST) will be able to detect dwarf satellites of high-z Lyman break galaxies (LBGs). To this end, we use cosmological simulations following the evolution of a typical M ⋆ ≃ 1010 M ⊙ LBG up to z ≃ 6, and analyze the observational properties of its five satellite dwarf galaxies (107 M ⊙ &lt; M ⋆ &lt; 109 M ⊙). Modeling their stellar emission and dust attenuation, we reconstruct their rest-frame ultraviolet-optical spectra for 6 &lt; z &lt; 6.5. JWST/NIRCam synthetic images show that the satellites can be spatially resolved from their host, and their emission is detectable by planned deep surveys. Moreover, we build synthetic spectral energy distributions and color–magnitude diagrams for the satellites. We conclude that the color F200W–F356W is a powerful diagnostic tool for understanding their physical properties once they have been identified. For example, F200W–F356W ≲ −0.25 can be used to identify star-bursting (SFR ∼ 5 M ⊙ yr−1), low-mass (M ⋆ ≲ 5 × 108 M ⊙) systems, with ∼80% of their stars being young and metal poor ].

Numerical study of the Kelvin-Helmholtz instability and its effect on synthetic emission from magnetized jets

Context.Nonthermal emission from active galactic nucleus (AGN) jets extends up to large scales, even though they are prone to many magnetohydrodynamic instabilities.Aims.The main focus of this study is to understand the effect of magnetohydrodynamic instabilities on the nonthermal emission from large-scale AGN jets.Methods.We performed high-resolution three-dimensional numerical magnetohydrodynamic simulations of a plasma column to investigate the dynamical and emission properties of jet configurations at kiloparsec scales with different magnetic field profiles, jet speeds, and density contrast. We also obtained synthetic nonthermal emission signatures for different viewing angles using an approach that assumes static particle spectra and that is obtained by evolving the particle spectra using Lagrangian macroparticles incorporating the effects of shock acceleration and radiative losses.Results.We find that shocks due to the Kelvin-Helmholtz (KH) instability in the axial magnetic field configurations can strongly affect the jet dynamics. Additionally, we also find weak biconical shocks in the underdense jet columns. The helical magnetic field hinders the vortex growth at the shear surface and thereby stabilizes the jet column. With the evolving particle spectra approach, the synthetic spectral energy distributions obtained for cases with strong KH instability show multiple humps ranging from the radio to the TeV gamma-ray band.Conclusions.We conclude that high-energy electrons that are accelerated in the vicinity of freshly formed shocks due to KH instability result in high X-ray emission.

Understanding the evolution and dust formation of carbon stars in the Large Magellanic Cloud via the JWST

Context. Carbon stars have been, and still are, extensively studied. Given their complex internal structure and their peculiar chemical composition, they are living laboratories in which we can test stellar structure and evolution theories of evolved stars. Furthermore, they are the most relevant dust manufacturers, thus playing a crucial role in the evolution of galaxies. Aims. We aim to study the dust mineralogy of the circumstellar envelope of carbon stars in the Large Magellanic Cloud (LMC) to achieve a better understanding of the dust formation process in the outflow of these objects. We intend to investigate the expected distribution of carbon stars in the observational planes built with the filters of the Mid-Infrared Instrument (MIRI) mounted onboard the James Webb Space Telescope (JWST) to select the best planes allowing an exhaustive characterisation of the stars. Methods. We compared the synthetic spectral energy distributions, obtained by modelling asymptotic giant branch stars and the dust formation process in the wind, with the spectra of carbon stars in the LMC, taken with the Infrared Spectrograph onboard the Spitzer Space Telescope. From the detailed comparison between synthetic modelling and observation we characterise the individual sources and derive the detailed mineralogy of the dust in the circumstellar envelope. Results. The sample of stars considered here is composed of stars of diverse mass, formation epoch, degree of obscuration, and metallicity. We find that precipitation of MgS on SiC seeds is common to all non-metal-poor carbon stars. Solid carbon is the dominant dust component, with percentages above 80% in all cases; a percentage between 10% and 20% of carbon dust is under the form of graphite, the remaining being amorphous carbon. Regarding the observational planes based on the MIRI filters, the colour-magnitude ([F770W]–[F1800W], [F1800W]) plane allows the best understanding of the degree of obscuration of the stars, while the ([F1800W]–[F2550W], [F1800W]) diagram allows better discrimination among stars of different metallicities.

Infrared emission of z ∼ 6 galaxies: AGN imprints

ABSTRACT We investigate the infrared (IR) emission of high-redshift (z ∼ 6), highly star-forming (${{\rm SFR}\gt 100\,{\rm M}_{\odot }\, {\rm yr}^{-1}}$) galaxies, with/without active galactic nuclei (AGN), using a suite of cosmological simulations featuring dust radiative transfer. Synthetic spectral energy distributions (SEDs) are used to quantify the relative contribution of stars/AGN to dust heating. In dusty (Md ≳ 3 × 107 M⊙) galaxies, ≳50–90 per cent of the ultraviolet (UV) radiation is obscured by dust inhomogeneities on scales ≳100 pc. In runs with AGN, a clumpy, warm (≈250 K) dust component coexists with a colder (≈60 K) and more diffuse one, heated by stars. Warm dust provides up to ${50 {{\ \rm per\ cent}}}$ of the total infrared (IR) luminosity, but only ${\lesssim}0.1 {{\ \rm per\ cent}}$ of the total mass content. The AGN boosts the MIR flux by 10–100 times with respect to star-forming galaxies, without significantly affecting the far-IR. Our simulations successfully reproduce the observed SED of bright (MUV ∼ −26) z ∼ 6 quasars, and show that these objects are part of complex, dust-rich merging systems, containing multiple sources (accreting black holes and/or star-forming galaxies) in agreement with recent HST and ALMA observations. Our results show that the proposed ORIGINS missions will be able to investigate the mid-IR (MIR) properties of dusty star-forming galaxies and to obtain good-quality spectra of bright quasars at z ∼ 6. Finally, the MIR-to-FIR flux ratio of faint (MUV ∼ −24) AGN is &amp;gt;10 times higher than for normal star-forming galaxies. This implies that combined JWST/ORIGINS/ALMA observations will be crucial to identify faint and/or dust-obscured AGN in the distant Universe.

The GRAVITY young stellar object survey

Context. Studies of the dust distribution, composition, and evolution of protoplanetary disks provide clues for understanding planet formation. However, little is known about the innermost regions of disks where telluric planets are expected to form. Aims. We aim constrain the geometry of the inner disk of the T Tauri star RY Lup by combining spectro-photometric data and interferometric observations in the near-infrared (NIR) collected at the Very Large Telescope Interferometer. We use PIONIER data from the ESO archive and GRAVITY data that were obtained in June 2017 with the four 8m telescopes. Methods. We use a parametric disk model and the 3D radiative transfer code MCFOST to reproduce the spectral energy distribution (SED) and match the interferometric observations. MCFOST produces synthetic SEDs and intensity maps at different wavelengths from which we compute the modeled interferometric visibilities and closure phases through Fourier transform. Results. To match the SED from the blue to the millimetric range, our model requires a stellar luminosity of 2.5 L⊙, higher than any previously determined values. Such a high value is needed to accommodate the circumstellar extinction caused by the highly inclined disk, which has been neglected in previous studies. While using an effective temperature of 4800 K determined through high-resolution spectroscopy, we derive a stellar radius of 2.29 R⊙. These revised fundamental parameters, when combined with the mass estimates available (in the range 1.3–1.5 M⊙), lead to an age of 0.5–2.0 Ma for RY Lup, in better agreement with the age of the Lupus association than previous determinations. Our disk model (that has a transition disk geometry) nicely reproduces the interferometric GRAVITY data and is in good agreement with the PIONIER ones. We derive an inner rim location at 0.12 au from the central star. This model corresponds to an inclination of the inner disk of 50°, which is in mild tension with previous determinations of a more inclined outer disk from SPHERE (70° in NIR) and ALMA (67 ± 5°) images, but consistent with the inclination determination from the ALMA CO spectra (55 ± 5°). Increasing the inclination of the inner disk to 70° leads to a higher line-of-sight extinction and therefore requires a higher stellar luminosity of 4.65 L⊙ to match the observed flux levels. This luminosity would translate to a stellar radius of 3.13 R⊙, leading to an age of 2–3 Ma, and a stellarmass of about 2 M⊙, in disagreement with the observed dynamical mass estimate of 1.3–1.5 M⊙. Critically, this high-inclination inner disk model also fails to reproduce the visibilities observed with GRAVITY. Conclusions. The inner dust disk, as traced by the GRAVITY data, is located at a radius in agreement with the dust sublimation radius. An ambiguity remains regarding the respective orientations of the inner and outer disk, coplanar and mildly misaligned, respectively.As our datasets are not contemporary and the star is strongly variable, a deeper investigation will require a dedicated multi-technique observing campaign.

KELT-17: a chemically peculiar Am star and a hot-Jupiter planet

Context. There is very little information to be found in the literature regarding the detection of planets orbiting chemically peculiar stars. Aims. Our aim is to determine the detailed chemical composition of the remarkable planet host star KELT-17. This object hosts a hot-Jupiter planet with 1.31 MJup detected by transits, and it is one of the more massive and rapidly rotating planet hosts seen to date. We set out to derive a complete chemical pattern for this star, in order to compare it with those of chemically peculiar stars. Methods. We carried out a detailed abundance determination in the planet host star KELT-17 via spectral synthesis. Stellar parameters were estimated iteratively by fitting Balmer line profiles and imposing the Fe ionization balance using the SYNTHE program together with plane-parallel ATLAS12 model atmospheres. Specific opacities for an arbitrary composition and microturbulence velocity vmicro were calculated through the opacity sampling (OS) method. The abundances were determined iteratively by fitting synthetic spectra to metallic lines of 16 different chemical species using SYNTHE. The complete chemical pattern of KELT-17 was compared to the recently published average pattern of Am stars. We estimated the stellar radius using two methods: a) comparing the synthetic spectral energy distribution with the available photometric data and the Gaia parallax, and b) using a Bayesian estimation of stellar parameters using stellar isochrones. Results. We found over-abundances of Ti, Cr, Mn, Fe, Ni, Zn, Sr, Y, Zr, and Ba, together with subsolar values of Ca and Sc. Notably, the chemical pattern agrees with those recently published for Am stars, making KELT-17 the first exoplanet host whose complete chemical pattern is unambiguously identified with this class. The stellar radius derived by two different methods agrees to each other and with those previously obtained in the literature.

Mentari: A pipeline to model the galaxy SED using semi analytic models

AbstractWe build a theoretical picture of how the light from galaxies evolves across cosmic time. In particular, we predict the evolution of the galaxy spectral energy distribution (SED) by carefully integrating the star formation and metal enrichment histories of semi-analytic model (SAM) galaxies and combining these with stellar population synthesis models which we call mentari. Our SAM combines prescriptions to model the interplay between gas accretion, star formation, feedback process, and chemical enrichment in galaxy evolution. From this, the SED of any simulated galaxy at any point in its history can be constructed and compared with telescope data to reverse engineer the various physical processes that may have led to a particular set of observations. The synthetic SEDs of millions of simulated galaxies from mentari can cover wavelengths from the far UV to infrared, and thus can tell a near complete story of the history of galaxy evolution.