WISE Bands Research Articles

Exploring Be phenomena in OBA stars: A mid-infrared search

As early-type stars with a rotation speed close to their critical velocity, Be stars experience an event called the Be phenomenon. The material in their equator is ejected into outside space during the Be phenomenon and forms a circumstellar disk. The mechanism triggering these events remains poorly understood, and observations of these events are limited because the duration of these events ranges from months to years. Long-term epoch photometry in the infrared bands is expected to be ideal for detecting Be phenomena because the brightness variation is larger than that in the optical, and the effect of interstellar extinction is weaker as well. We conducted a systematic search for Be phenomena among Milky Way OBA stars in the mid-infrared. We examined the brightness and colour variations of known classical Be stars using the WISE W1 and W2 photometry bands to quantify their characteristics. Subsequently, we established a set of criteria to identify similar photometric variations in a large sample of OBA stars. We found 916 OBA stars that show Be phenomena in the past 13 yr, 736 of which are newly discovered. The peak-to-peak variations in magnitude and colour were found to be correlated, indicating that a decretion disk is common. The increase in colour was observed to be strongly correlated with the emission of the Hα line, providing further evidence of the association with circumstellar disks. The brightness variation of a star with Be phenomena can be up to 1.5 mag, and the colour variations can be up to 0.4 mag. The median durations for the disk build-up and decay phases are 474 and 524 days, respectively (durations shorter than 180 days are not sampled). The search for Be phenomena in the WISE bands greatly enlarges the number of stars showing disk variation, and it enables multi-band photometry analysis of these events with the help of current and future optical photometry surveys.

A versatile classification tool for galactic activity using optical and infrared colors

Context.The overwhelming majority of diagnostic tools for galactic activity are focused mainly on the classes of active galaxies. Passive or dormant galaxies are often excluded from these diagnostics, which usually employ emission-line features (e.g., forbidden emission lines). Thus, most of them focus on specific types of activity or only on one activity class, for example active galactic nucleus (AGN) galaxiesAims.In this work we used infrared and optical colors to build an all-inclusive galactic activity diagnostic tool that can discriminate between star-forming, AGN, low-ionization nuclear emission-line region, composite, and passive galaxies, and which can be used in local and low-redshift galaxies.Methods.We used the random forest algorithm to define a new activity diagnostic tool. As the ground truth for the training of the algorithm, we considered galaxies that have been classified based on their optical spectral lines. We explored classification criteria based on infrared colors from the first three WISE bands (bands 1, 2, and 3) supplemented with optical colors from theu, g, andrSDSS bands. From them, we sought the combination with the minimum number of colors that provides optimal results. Furthermore, to mitigate biases related to aperture effects, we introduced a new WISE photometric scheme that combines apertures of different sizes.Results.Using machine learning methods, we developed a diagnostic tool that accommodates both active and passive galaxies under one unified classification scheme using just three colors. We find that the combination ofW1-W2, W2-W3, and g-r colors offers a good performance, while the broad availability of these colors for a large number of galaxies ensures it can be applied to large galaxy samples. The overall accuracy is ~81%, and the achieved completeness for each class is ~81% for star-forming, ~56% for AGN, ~68% for LINER, ~65% for composite, and ~85% for passive galaxies.Conclusions.Our diagnostic represents a significant improvement over existing infrared diagnostics because it includes all types of active galaxies, as well as passive galaxies, extending their application to the local Universe. The inclusion of the optical colors improves its ability to identify low-luminosity AGN galaxies, which are generally confused with star-forming galaxies, and helps us identify cases of starbursts with extreme mid-infrared colors that mimic obscured AGN galaxies, a well-known problem for most infrared diagnostics.

Long-period maser-bearing Miras in the Galactic center

We establish a sample of 370 Mira variables that are likely near the Galactic center (GC). The sources have been selected from the OGLE and BAaDE surveys based on their sky coordinates, OGLE classifications, and BAaDE maser-derived line-of-sight velocities. As the distance to the GC is known to a high accuracy, this sample is a test bed for reddening and extinction studies toward the GC and in Mira envelopes. We calculated separate interstellar- and circumstellar-extinction values for individual sources, showing that there is a wide range of circumstellar extinction values (up to four magnitudes in the Ks band) in the sample, and that circumstellar reddening is statistically different from interstellar reddening laws. Further, the reddening laws in the circumstellar environments of our sample and the circumstellar environments of Large Magellanic Cloud (LMC) Miras are strikingly similar despite the different metallicities of the samples. Period-magnitude relations for the mid-infrared (MIR) WISE and MSX bands are also explored, and in the WISE bands we compare these to period-magnitude relationships derived from Miras in the LMC as it is important to compare these LMC relations to those in a higher metallicity environment. Emission from the envelope itself may contaminate MIR magnitudes altering the relations, especially for sources with thick envelopes.

Star formation rate and stellar mass calibrations based on infrared photometry and their dependence on stellar population age and extinction

Context. The stellar mass (M⋆) and the star formation rate (SFR) are among the most important features that characterize galaxies. Measuring these fundamental properties accurately is critical for understanding the present state of galaxies, their history, and future evolution. Infrared (IR) photometry is widely used to measure the M⋆ and SFR of galaxies because the near-IR traces the continuum emission of the majority of their stellar populations (SPs), and the mid/far-IR traces the dust emission powered by star-forming activity. Aims. This work explores the dependence of the IR emission of galaxies on their extinction, and the age of their SPs. It aims to provide accurate and precise IR-photometry SFR and M⋆ calibrations that account for SP age and extinction while providing quantification of their scatter. Methods. We used the CIGALE spectral energy distribution (SED) fitting code to create model SEDs of galaxies with a wide range of star formation histories, dust content, and interstellar medium properties. We fit the relations between M⋆ and SFR with IR and optical photometry of the model-galaxy SEDs with the Markov chain Monte Carlo (MCMC) method. As an independent confirmation of the MCMC fitting method, we performed a machine-learning random forest (RF) analysis on the same data set. The RF model yields similar results to the MCMC fits, thus validating the latter. Results. This work provides calibrations for the SFR using a combination of the WISE bands 1 and 3, or the JWST NIR-F200W and MIRI-F2100W. It also provides mass-to-light ratio calibrations based on the WISE band-1, the JWST NIR-F200W, and the optical u − r or g − r colors. These calibrations account for the biases attributed to the SP age, while they are given in the form of extinction-dependent and extinction-independent relations. Conclusions. The proposed calibrations show robust estimations while minimizing the scatter and biases throughout a wide range of SFRs and stellar masses. The SFR calibration offers better results, especially in dust-free or passive galaxies where the contributions of old SPs or biases from the lack of dust are significant. Similarly, the M⋆ calibration yields significantly better results for dusty and high-SFR galaxies where dust emission can otherwise bias the estimations.

Mid-infrared and Maser Flux Variability Correlation in Massive Young Stellar Object G036.70+00.09

We present the discovery of the simultaneous flux variation of a massive young stellar object (MYSO) G036.70+00.09 (G036.70) both in the maser emission and mid-infrared (MIR; λ = 3–5 μm) bands. Using the ALLWISE and NEOWISE archival databases that cover a long time span of approximately 10 yr with a cadence of 6 months, we confirm that G036.70 indicates a stochastic year-long MIR variability with no signs of a WISE band color change of W1 (3.4 μm) −W2 (4.6 μm). Cross-matching the MIR data set with the high-cadence 6.7 GHz class II methanol maser flux using the Hitachi 32 m radio telescope that discovered its periodicity in the methanol maser of 53.0–53.2 days, we also determine the flux correlations between the two bands at two different timescales, year-long and day-long, both of which have never been reported in MYSOs, except when they are in the accretion burst phase. The results of our study support the scenario that a class II methanol maser is pumped up by infrared emission from accreting disks of MYSOs. We also discuss the possible origins of MIR and maser variability. To explain the two observed phenomena, a stochastic year-long MIR variability with no signs of significant color change and maser-MIR variability correlation or a change in mass accretion rate and line-of-sight extinction because of the nonaxisymmetric dust density distribution in a rotating accretion disk are possible origins. Observations through spectroscopic monitoring of accretion-related emission lines are essential for determining the origin of the observed variability in G036.70.

A Resolved Search for AGN in the Centers of Nearby Galaxies with WISE

We investigate the mid-infrared detection of active galactic nucleus (AGN) in nearby galaxies and the effects of contamination by starlight and dust emission from the host galaxy. We used a sample of ∼14,000 nearby galaxies from the z = 0 Multiwavelength Galaxy Synthesis to study the effect of isolating the galaxy center on AGN detection. We identify 39 known AGN and 9 additional previously selected AGN candidates using color selection in WISE Bands 1 and 2. Most AGN identified using the smallest apertures around the nucleus are lost from the WISE color selection at larger aperture radii. This suggests that when low spatial resolution or integrated galaxy measurements are used to identify AGN in the mid-IR, lower luminosity AGN may go undetected due to contamination from the galaxy’s dust emission and starlight.

High-contrast, high-angular resolution view of the GJ 367 exoplanet system

ABSTRACT We search for additional companions in the GJ 367 exoplanet system and aim to better constrain its age and evolutionary status. We analyse high-contrast direct imaging observations obtained with HST/NICMOS, VLT/NACO, and VLT/SPHERE. We investigate and critically discuss conflicting age indicators based on theoretical isochrones and models for Galactic dynamics. A comparison of GAIA EDR3 parallax and photometric measurements with theoretical isochrones suggests a young age ≤60 Myr for GJ 367. The star’s Galactic kinematics exclude membership to any nearby young moving group or stellar stream. Its highly eccentric Galactic orbit, however, is atypical for a young star. Age estimates considering Galactic dynamical evolution are most consistent with an age of 1–8 Gyr. We find no evidence for a significant mid-infrared excess in the WISE bands, suggesting the absence of warm dust in the GJ 367 system. The direct imaging data provide significantly improved detection limits compared to previous studies. At 530 mas (5 au) separation, the SPHERE data achieve a 5σ contrast of 2.6 × 10−6. The data exclude the presence of a stellar companion at projected separations ≥0.4 au. At projected separations ≥5 au we can exclude substellar companions with a mass ≥1.5 MJup for an age of 50 Myr, and ≥20 MJup for an age of 5 Gyr. By applying the stellar parameters corresponding to the 50 Myr isochrone, we derive a bulk density of ρplanet = 6.2 g cm−3 for GJ 367 b, which is 25 per cent smaller than a previous estimate.

Herbig Ae/Be Stars toward the Dark Cloud LDN 1667*

Abstract We report the discovery of a new emission-line object, named SPH 4-South = (GAIA EDR3 5616553300192230272), toward the dark cloud LDN 1667. This object came to our attention after inspecting public images that show a faint diffuse nebula a few arcseconds south of SPH 4, an emission-line object previously classified as a T Tauri star. We present high-resolution spectra and analyzed JHK photometry of SPH 4 and SPH 4-South and new narrowband and archival broadband images of these objects. A comparison of the spectra of SPH 4 and SPH 4-South with high-resolution ones of DG Cir and R Mon strongly suggests that SPH 4 and SPH 4-South are Herbig Ae/Be stars. The classification of SPH 4-South is further supported by using a k-NN algorithm to its position in an H−K versus J−H color–color diagram. Both stars are detected in the four WISE bands and the WISE colors allow us to classify SPH 4 as a Class I and SPH 4-South as a Class II source. We also show that the faint nebula is most probably associated with SPH 4-South. Using published results on LDN 1667 and the Gaia Early Data Release 3, we conclude that SPH 4 is a member of LDN 1667. The case of SPH 4-South is not clear because the determination of its distance and proper motion could be affected by the nebulosity around the star, although membership of SPH 4-South to LDN 1667 cannot be ruled out.

A possible far-ultraviolet flux-dependent core mass function in NGC 6357

Context. NGC 6357 is a galactic star-forming complex (d ~ 1.7 kpc) composed of several HII regions, a few young stellar clusters, and giant molecular clouds. In particular, the HII regions G353.2+0.9, G353.1+0.6, and G353.2+0.7 are associated with three young clusters; the most prominent of these, Pismis 24, contains some of the most massive stars known. Aims. We aim to derive the properties of the densest compact gas structures (cores) in the region as well as the effects of an intense far-ultraviolet (FUV) radiation field on their global properties. Methods. We mapped the NGC 6357 region at 450 and 850 μm with SCUBA-2 and in the CO(3–2) line with HARP at the JCMT. We also made use of the Herschel Hi-GAL data at 70 and 160 μm. We used the algorithm Gaussclumps to retrieve the compact cores embedded in the diffuse sub-millimetre emission and constructed their spectral energy distribution from 70 to 850 μm, from which we derived mass and temperature. We divided the observed area into an ‘active’ region (i.e. the eastern half, which is exposed to the FUV radiation from the more massive members of the three clusters) and a ‘quiescent’ region (i.e. the western half, which is less affected by FUV radiation). We compared the core mass functions and the temperature distributions in the two areas to look for any differences that could be due to the different levels of FUV radiation. Results. We retrieved 686 dense cores, 411 in the active region and 275 in the quiescent region, with an estimated mass completeness limit of ~5 M⊙. We also attempted to select a sample of pre-stellar cores based on cross-correlation with 70 μm emission and red WISE point sources, which unfortunately is biased due to distance, emission at 70 μm from the dust on the surface of the cores that is heated by the FUV radiation, and saturation in the WISE bands. Most of the cores above the mass completeness limit are likely to be gravitationally bound. The fraction of gas in dense cores is very low, 1.4%. We found a mass-size relation log(M∕M⊙) ~ a × log(D∕arcsec), with a in the range 2.0–2.4, depending on the precise selection of the sample. The temperature distributions in the two sub-regions are clearly different, peaking at ~25 K in the quiescent region and at ~35 K in the active region. The core mass functions are different as well, at a 2σ level, consistent with a Salpeter initial mass function in the quiescent region and flatter than that in the active region. The dense cores lying close to the HII regions are consistent with pre-existing cores being gradually engulfed by a photon dominated region and photoevaporating. A comparison of the obtained distribution of core masses with those derived from simulations of cloud-cloud collisions yields no conclusive evidence of ongoing cloud-cloud collisions. Conclusions. We attribute the different global properties of dense cores in the two sub-regions to the influence of the FUV radiation field.

Forward modelling the multiwavelength properties of active galactic nuclei: application to X-ray and WISE mid-infrared samples

ABSTRACT An empirical forward-modelling framework is developed to interpret the multiwavelength properties of active galactic nuclei (AGNs) and provide insights into the overlap and incompleteness of samples selected at different parts of the electromagnetic spectrum. The core of the model are observationally derived probabilites on the occupation of galaxies by X-ray-selected AGN. These are used to seed mock galaxies drawn from stellar-mass functions with accretion events and then associate them with spectral energy distributions that describe both the stellar and AGN emission components. This approach is used to study the complementarity between X-ray and WISE mid-infrared AGN selection methods. We first show that the basic observational properties of the X-ray and WISE AGN (magnitude and redshift distributions) are adequately reproduced by the model. We then infer the level of contamination of the WISE selection and show that this is dominated by non-AGN at redshifts z < 0.5. These are star-forming galaxies that scatter into the WISE AGN selection wedge because of photometric uncertainties affecting their colours. Our baseline model shows a sharp drop in the number density of heavily obscured AGN above the Compton-thick limit in the WISE bands. The model also overpredicts by a factor of 1.5 the fraction of X-ray associations in the WISE AGN selection box compared to observations. This suggests a population of X-ray faint sources that is not reproduced by the model. This discrepancy is discussed in the context of either heavily obscured or intrinsically X-ray weak AGN. Evidence is found in favour of the latter.

The XXL Survey

The combination of optical and mid-infrared (MIR) photometry has been extensively used to select red active galactic nuclei (AGNs). Our aim is to explore the obscuration properties of these red AGNs with both X-ray spectroscopy and spectral energy distributions (SEDs). In this study, we re-visit the relation between optical/MIR extinction and X-ray absorption. We use IR selection criteria, specifically the W1 and W2 WISE bands, to identify 4798 AGNs in the XMM-XXL area (∼25 deg2). Application of optical/MIR colours (r−W2 > 6) reveals 561 red AGNs (14%). Of these, 47 have available X-ray spectra with at least 50 net (background-subtracted) counts per detector. For these sources, we construct SEDs from the optical to the MIR using the CIGALE code. The SED fitting shows that 44 of these latter 47 sources present clear signs of obscuration based on the AGN emission and the estimated inclination angle. Fitting the SED also reveals ten systems (∼20%) which are dominated by the galaxy. In these cases, the red colours are attributed to the host galaxy rather than AGN absorption. Excluding these ten systems from our sample and applying X-ray spectral fitting analysis shows that up to 76% (28/37) of the IR red AGNs present signs of X-ray absorption. Thus, there are nine sources (∼20% of the sample) that although optically red, are not substantially X-ray absorbed. Approximately 50% of these sources present broad emission lines in their optical spectra. We suggest that the reason for this apparent discrepancy is that the r−W2 criterion is sensitive to smaller amounts of obscuration relative to the X-ray spectroscopy. In conclusion, it appears that the majority of red AGNs present considerable obscuration levels as shown by their SEDs. Their X-ray absorption is moderate with a mean of NH ∼ 1022 cm−2.

2MASS J15460752−6258042: a mid-M dwarf hosting a prolonged accretion disc

ABSTRACT We report the discovery of the oldest (∼55 Myr) mid-M type star known to host ongoing accretion. 2MASS J15460752–6258042 (2M1546, spectral type M5, 59.2 pc) shows spectroscopic signs of accretion such as strong H α, He i, and [O i] emission lines, from which we estimate an accretion rate of ∼10−10 M⊙ yr−1. Considering the clearly detected infrared excess in all WISE bands, the shape of its spectral energy distribution (SED) and its age, we believe that the star is surrounded by a transitional disc, clearly with some gas still present at inner radii. The position and kinematics of the star from Gaia DR2 and our own radial-velocity measurements suggest membership in the nearby ∼55 Myr-old Argus moving group. At only 59 pc from Earth, 2M1546 is one of the nearest accreting mid-M dwarfs, making it an ideal target for studying the upper limit on the lifetimes of gas-rich discs around low-mass stars.

Sporadic and intense accretion in a 1 Myr-old brown dwarf candidate

Context. Studying the accretion process in very low-mass objects has important implications for understanding their formation mechanism. Many nearby late-M dwarfs that have previously been identified in the field are in fact young brown dwarf members of nearby young associations. Some of them are still accreting. They are therefore excellent targets for further studies of the accretion process in the very low-mass regime at different stages. Aims. We aim to search for accreting young brown dwarf candidates in a sample of 85 nearby late-M dwarfs. Methods. Using photometric data from DENIS, 2MASS, and WISE, we constructed the spectral energy distribution of the late- M dwarfs based on BT-Settl models to detect infrared excesses. We then searched for lithium and Hα emission in candidates that exhibit infrared excesses to confirm their youth and the presence of accretion. Results. Among the 85 late-M dwarfs, only DENIS-P J1538317−103850 (M5.5) shows strong infrared excesses in WISE bands. The detection of lithium absorption in the M5.5 dwarf and its Gaia trigonometric parallax indicate an age of ~1 Myr and a mass of 47 MJ. The Hα emission line in the brown dwarf shows significant variability that indicates sporadic accretion. This 1 Myr-old brown dwarf also exhibits intense accretion bursts with accretion rates of up to 10−7.9 M⊙ yr−1. Conclusions. Our detection of sporadic accretion in one of the youngest brown dwarfs might imply that sporadic accretion at early stages could play an important role in the formation of brown dwarfs. Very low-mass cores would not be able to accrete enough material to become stars, and thus they end up as brown dwarfs.

Infrared properties of planetary nebulae with [WR] and wels central stars

ABSTRACT We report the infrared (IR) properties of planetary nebulae (PNe) with Wolf–Rayet (WR) type and wels central stars known to date and compare them with the IR properties of a sample of PNe with H-rich central stars. We use near-, mid-, and far-IR photometric data from archives to derive the IR properties of PNe. We have constructed IR colour–colour diagrams of PNe using measurements from 2MASS, IRAS, WISE, and Akari bands. [WR] PNe have a larger near-IR emission from the hot dust component and also show a tendency for stronger 12 μm emission as compared to the other two groups. Cool asymptotic giant branch dust properties of all PNe are found to be similar. We derived the dust colour temperatures, dust masses, dust-to-gas mass ratios, IR luminosities, and IR excess (IRE) of PNe for these three groups. [WR] PNe and wels-PNe tend to have larger mean values for dust mass when compared to the third group. The average dust-to-gas mass ratio is found to be similar for the three groups of PNe. While there is a strong correlation of dust temperature and IR luminosity with the age for the three groups of PNe, the dust mass, dust-to-gas mass ratios, and IRE are found to be non-varying as the PNe evolve. [WR] PNe and wels-PNe show very similar distribution of excitation classes and also show similar distribution with Galactic latitude.

ATM: An open-source tool for asteroid thermal modeling and its application to NEOWISE data

We publicly release ATM, a Python package designed to model asteroid flux measurements to estimate an asteroid's size, surface temperature distribution, and emissivity. A number of the most popular static asteroid thermal models (NEATM, STM, FRM) are implemented with the reflected solar light contribution and Kirchhoff's law accounted for. Priors for fitted parameters can be easily specified and the solution, including the full multi-dimensional posterior probability density function, is found using Markov Chain Monte Carlo (MCMC). We describe the package's architecture and discuss model and fitting validation. Data files with ~10 million WISE flux measurements for ~150,000 unique asteroids and additional Minor Planet Center data are also included with the package, as well as Python Jupyter Notebooks with examples of how to select subsamples of objects, filter and process data, and use ATM to fit for the desired model parameters. The entirety of the analysis presented here, including all the figures, tables, and catalogs, can be easily reproduced with these publicly released Notebooks. We show that ATM can match the best-fit size estimates for well-observed asteroids published in 2016 by the NEOWISE team (Mainzer et al., 2016) with a sub-percent bias and a scatter of only 6%. Our analysis of various sources of random and systematic size uncertainties shows that for the majority of over 100,000 objects with WISE-based size estimates random uncertainties (precision) are about 10%; systematic uncertainties within the adopted model framework, such as NEATM, and with assumed emissivity for WISE W3 and W4 bands, are likely in the range of 10–20%. Hence, the accuracy of WISE-based asteroid size estimates is approximately in the range of 15–20% for most objects, except for unknown errors due to a possibly over-simplified modeling framework (e.g., spherical asteroid approximation). We compare model families to data in WISE color-color diagrams and derive a simple method to estimate size that only uses WISE W3 band flux; we show that it matches estimates based on all four WISE bands to within 10%. We also study optical data collected by the Sloan Digital Sky Survey (SDSS) and show that correlations of optical colors and WISE-based best-fit model parameters indicate robustness of the latter. Our analysis gives support to the claim by Harris & Drube (2014) that candidate metallic asteroids can be selected using the best-fit temperature parameter and infrared albedo. We investigate a correlation between SDSS colors and optical albedo derived using WISE-based size estimates and show that this correlation can be used to estimate asteroid sizes with optical data alone, with a precision of about 17% relative to WISE-based size estimates. After accounting for systematic errors, the difference in accuracy between infrared and optical color-based size estimates becomes less than a factor of two.

An empirical examination of WISE/NEOWISE asteroid analysis and results

Asteroid observations by the WISE space telescope and the analysis of those observations by the NEOWISE project have provided more information about the diameter, albedo, and other properties of approximately 164,000 asteroids, more than all other sources combined. The raw data set from this mission will likely be the largest and most important such data on asteroids available for many years. To put this trove of data to productive use, we must understand its strengths and weaknesses, and we need clear and reproducible methods for analyzing the data set. This study critically examines the WISE observational data and the NEOWISE results published in both the original papers and the NASA Planetary Data System (PDS). There seem to be multiple areas where the analysis might benefit from improvement or independent verification. The NEOWISE results were obtained by the application of 10 different modeling methods, many of which are not adequately explained or even defined, to 12 different combinations of WISE band data. More than half of NEOWISE results are based on a single band of data. The majority of curve fits to the data in the NEOWISE results are of poor quality, frequently missing most or all of the data points on which they are based. Complete misses occur for about 30% of single-band results, and among the results derived from the most common multiple-band combinations, about 43% miss all data points in at least one band. The NEOWISE data analysis relies on assumptions that are in many cases inconsistent with each other. A substantial fraction of WISE data was systematically excluded from the NEOWISE analysis. Building on methods developed by Hanuš et al. (2015), I show that error estimates for the WISE observational data were not well characterized, and all observations have true uncertainty at least a factor of 1.3–2.5 times larger than previously described, depending on the band. I also show that the error distribution is not well fit by a normal distribution. These findings are important because the Monte Carlo error-analysis method used by the NEOWISE project depends on both the observational errors and the normal distribution. An empirical comparison of published NEOWISE diameters to those in the literature that were estimated by using radar, occultation, or spacecraft (ROS) measurements shows that, for 129 results involving 105 asteroids, the NEOWISE diameters presented in tables of thermal-modeling results exactly match prior ROS results from the literature. While these are only a tiny fraction (0.06%) of the asteroids analyzed, they are important because they represent the only independent check on NEOWISE diameter accuracy. After removing the exact matches and adding additional ROS results, I find that the accuracy of diameter estimates for NEOWISE results depends strongly on the choice of data bands and on which of the 10 models was used. I show that systematic errors in the diameter estimates are much larger than previously described and range from − 5% to + 23%. In addition, random errors range from − 15% to + 19% when all four WISE bands were used, and from − 39% to + 57% in cases employing only the W2 band. The empirical results presented here show that much work remains to be done in analyzing data from the WISE/NEOWISE mission and interpreting it for asteroid science.

Asteroid thermal modeling in the presence of reflected sunlight

A new derivation of simple asteroid thermal models is presented, investigating the need to account correctly for Kirchhoff's law of thermal radiation when IR observations contain substantial reflected sunlight. The framework applies to both the NEATM and related thermal models. A new parameterization of these models eliminates the dependence of thermal modeling on visible absolute magnitude H, which is not always available. Monte Carlo simulations are used to assess the potential impact of violating Kirchhoff's law on estimates of physical parameters such as diameter and IR albedo, with an emphasis on NEOWISE results. The NEOWISE papers use ten different models, applied to 12 different combinations of WISE data bands, in 47 different combinations. The most prevalent combinations are simulated and the accuracy of diameter estimates is found to be depend critically on the model and data band combination. In the best case of full thermal modeling of all four band the errors in an idealized model the 1σ (68.27%) confidence interval is −5% to +6%, but this combination is just 1.9% of NEOWISE results. Other combinations representing 42% of the NEOWISE results have about twice the CI at −10% to +12%, before accounting for errors due to irregular shape or other real world effects that are not simulated. The model and data band combinations found for the majority of NEOWISE results have much larger systematic and random errors. Kirchhoff's law violation by NEOWISE models leads to errors in estimation accuracy that are strongest for asteroids with W1, W2 band emissivity ϵ12 in both the lowest (0.605≤ϵ12≤0.780), and highest decile (0.969≤ϵ12≤0.988), corresponding to the highest and lowest deciles of near-IR albedo pIR. Systematic accuracy error between deciles ranges from a low of 5% to as much as 45%, and there are also differences in the random errors. Kirchhoff's law effects also produce large errors in NEOWISE estimates of pIR, particularly for high values. IR observations of asteroids in bands that have substantial reflected sunlight can largely avoid these problems by adopting the Kirchhoff law compliant modeling framework presented here, which is conceptually straightforward and comes without computational cost.

The Optical–Mid-infrared Extinction Law of the l = 165° Sightline in the Galactic Plane: Diversity of the Extinction Law in the Diffuse Interstellar Medium

Abstract Understanding the effects of dust extinction is important to properly interpret observations. The optical total-to-selective extinction ratio, , is widely used to describe extinction variations in ultraviolet and optical bands. Since the extinction curve adequately represents the average extinction law of diffuse regions in the Milky Way, it is commonly used to correct observational measurements along sightlines toward diffuse regions in the interstellar medium. However, the value may vary even along different diffuse interstellar medium sightlines. In this paper, we investigate the optical–mid-infrared (mid-IR) extinction law toward a very diffuse region at in the Galactic plane, which was selected based on a CO emission map. Adopting red clump stars as extinction tracers, we determine the optical–mid-IR extinction law for our diffuse region in two APASS bands ( ), three XSTPS-GAC bands ( ), three 2MASS bands ( ), and two WISE bands ( ). Specifically, 18 red clump stars were selected from the APOGEE–RC catalog based on spectroscopic data in order to explore the diversity of the extinction law. We find that the optical extinction curves exhibit appreciable diversity. The corresponding ranges from 1.7 to 3.8, while the mean value of 2.8 is consistent with the widely adopted average value of 3.1 for Galactic diffuse clouds. There is no apparent correlation between value and color excess in the range of interest, from 0.2 to 0.6 mag, or with specific visual extinction per kiloparsec, .

Quantifying the (X/peanut)-shaped structure of the Milky Way – new constraints on the bar geometry

The nature, size and orientation of the dominant structural components in the Milky Way's inner ~4 kpc - specifically the bulge and bar - have been the subject of conflicting interpretations in the literature. We present a different approach to inferring the properties of the long bar which extends beyond the inner bulge, via the information encoded in the Galaxy's X/peanut (X/P)-shaped structure. We perform a quantitative analysis of the X/P feature seen in WISE wide-field imaging at 3.4 $\mu$m and 4.6 $\mu$m. We measure the deviations of the isophotes from pure ellipses, and quantify the X/P structure via the radial profile of the Fourier n=6 harmonic (cosine term $B_6$). In addition to the vertical height and integrated `strength' of the X/P instability, we report an intrinsic radius of $R_{\Pi;{\rm int}} = 1.67\pm0.27$ kpc, and an orientation angle of $\alpha = {37^\circ}^{+7^\circ}_{-10^\circ}$ with respect to our line-of-sight to the Galactic Centre. Based on X/P structures observed in other galaxies, we make three assumptions: (i) the peanut is intrinsically symmetric, (ii) the peanut is aligned with the long Galactic bar, and (iii) their sizes are correlated. Thus the implication for the Galactic bar is that it is oriented at the same $37^\circ$ angle and has an expected radius of ~4.2 kpc, but possibly as low as ~3.2 kpc. We further investigate how the Milky Way's X/P structure compares with other analogues, and find that the Galaxy is broadly consistent with our recently established scaling relations, though with a moderately stronger peanut instability than expected. We additionally perform a photometric decomposition of the Milky Way's major axis surface brightness profile, accounting for spiral structure, and determine an average disc scale length of $h = 2.54\pm0.16$ kpc in the WISE bands, in good agreement with the literature.

Red clump stars and Gaia: calibration of the standard candle using a hierarchical probabilistic model

Distances to individual stars in our own Galaxy are critical in order to piece together the nature of its velocity and spatial structure. Core helium burning red clump (RC) stars have similar luminosities, are abundant throughout the Galaxy, and thus constitute good standard candles. We build a hierarchical probabilistic model to quantify the quality of RC stars as standard candles using parallax measurements from the first Gaia data release. A unique aspect of our methodology is to fully account for (and marginalize over) parallax, photometry, and dust corrections uncertainties, which leads to more robust results than standard approaches. We determine the absolute magnitude and intrinsic dispersion of the RC in 2MASS bands J, H, Ks, Gaia G band, and WISE bands W1, W2, W3, and W4. We find that the absolute magnitude of the RC is $-1.61 \pm$ 0.01 (in Ks), $+0.44 \pm$ 0.01 (in G) , $-0.93 \pm$ 0.01 (in J), $-1.46 \pm$ 0.01 (in H), $-1.68 \pm$ 0.02 (in W1), $-1.69\pm$ 0.02 (in W2), $-1.67 \pm$ 0.02 (in W3), $1.76 \pm$ 0.01 mag (in W4). The mean intrinsic dispersion is $\sim 0.17 \pm$ 0.03 mag across all bands (yielding a typical distance precision of $\sim$ 8%). Thus RC stars are reliable and precise standard candles. In addition, we have also re-calibrated the zero point of the absolute magnitude of the RC in each band, which provide a benchmark for future studies to estimate distances to RC stars. Finally, the parallax error shrinkage in the hierarchical model outlined in this work can be used to obtain more precise parallaxes than Gaia for the most distant RC stars across the Galaxy.