Abstract visibility
Abstract text size

Papers for Friday, Sep 03 2021

Papers with local authors

Alina Sabyr, Zoltán Haiman, José Manuel Zorrilla Matilla, Tianhuan Lu

15 pages, 15 figures, submitted to PRD

0 votes
Paper 1 — arXiv:2109.00547
0 votes
Paper 1 — arXiv:2109.00547

In order to extract full cosmological information from next-generation large and high-precision weak lensing (WL) surveys (e.g. Euclid, Roman, LSST), higher-order statistics that probe the small-scale, non-linear regime of large scale structure (LSS) need to be utilized. WL peak counts, which trace overdensities in the cosmic web, are one promising and simple statistic for constraining cosmological parameters. The physical origin of WL peaks have previously been linked to dark matter halos along the line of sight and this peak-halo connection has been used to develop various semi-analytic halo-based models for predicting peak counts. Here, we study the origin of WL peaks and the effectiveness of halo-based models for WL peak counts using a suite of ray-tracing N-body simulations. We compare WL peaks in convergence maps from the full simulations to those in maps created from only particles associated with halos -- the latter playing the role of a "perfect" halo model. We find that while halo-only contributions are able to replicate peak counts qualitatively well, halos do not explain all WL peaks. Halos particularly underpredict negative peaks, which are associated with local overdensities in large-scale underdense regions along the line of sight. In addition, neglecting non-halo contributions to peaks counts leads to a significant bias on the parameters ($\Omega_{\rm m}$, $\sigma_{8}$) for surveys larger than $\geq$ 100 deg$^{2}$. We conclude that other elements of the cosmic web, outside and far away from dark matter halos, need to be incorporated into models of WL peaks in order to infer unbiased cosmological constraints.

Jason Glenn, Charles M. Bradford, Erik Rosolowsky, Rashied Amini, Katherine Alatalo, Lee Armus, Andrew J. Benson, Tzu-Ching Chang, Jeremy Darling, Peter K. Day, Jeanette Domber, Duncan Farrah, Brandon Hensley, Sarah Lipscy, Bradley Moore, Seb Oliver, Joanna Perido, David Redding, Michael Rodgers, Raphael Shirley, Howard A. Smith, John B. Steeves, Carole Tucker, Jonas Zmuidzinas

84 pages, 24 figures, published in JATIS

0 votes
Paper 11 — arXiv:2109.00614
0 votes
Paper 11 — arXiv:2109.00614

The Galaxy Evolution Probe (GEP) is a concept for a mid- and far-infrared space observatory to measure key properties of large samples of galaxies with large and unbiased surveys. GEP will attempt to achieve zodiacal light and Galactic dust emission photon background-limited observations by utilizing a 6 Kelvin, 2.0 meter primary mirror and sensitive arrays of kinetic inductance detectors. It will have two instrument modules: a 10 - 400 micron hyperspectral imager with spectral resolution R = 8 (GEP-I) and a 24 - 193 micron, R = 200 grating spectrometer (GEP-S). GEP-I surveys will identify star-forming galaxies via their thermal dust emission and simultaneously measure redshifts using polycyclic aromatic hydrocarbon emission lines. Galaxy luminosities derived from star formation and nuclear supermassive black hole accretion will be measured for each source, enabling the cosmic star formation history to be measured to much greater precision than previously possible. Using optically thin far-infrared fine-structure lines, surveys with GEP-S will measure the growth of metallicity in the hearts of galaxies over cosmic time and extraplanar gas will be mapped in spiral galaxies in the local universe to investigate feedback processes. The science case and mission architecture designed to meet the science requirements are described, and the kinetic inductance detector and readout electronics state of the art and needed developments are described. This paper supersedes the GEP concept study report cited in it by providing new content, including: a summary of recent mid-infrared KID development, a discussion of microlens array fabrication for mid-infrared KIDs, and additional context for galaxy surveys. The reader interested in more technical details may want to consult the concept study report.

Meng Zhang, Mao-Sheng Xiang, Hua-Wei Zhang, Yuan-Sen Ting, Hans-Walter Rix, Ya-qian Wu, Yang Huang, Wei-Xiang Sun, Zhi-Jia Tian, Chun Wang, Xiao-Wei Liu

17 pages, 14 figures, 2 tables, accepted for publication in ApJ

0 votes
Paper 20 — arXiv:2109.00746
0 votes
Paper 20 — arXiv:2109.00746

Recent observations have revealed a population of $\alpha$-element abundances enhanced giant stars with unexpected high masses ($\gtrsim$1 $M_\odot$) from asteroseismic analysis and spectroscopy. Assuming single-star evolution, their masses imply young ages ($\tau<6$Gyr) incompatible with the canonical Galactic chemical evolution scenario. Here we study the chemistry and kinematics of a large sample of such $\alpha$-rich, high-mass red giant branch (RGB) stars drawn from the LAMOST spectroscopic surveys. Using LAMOST and Gaia, we found these stars share the same kinematics as the canonical high-$\alpha$ old stellar population in the Galactic thick disk. The stellar abundances show that these high-$\alpha$ massive stars have $\alpha$- and iron-peak element abundances similar to those of the high-$\alpha$ old thick disk stars. However, a portion of them exhibit higher [(N+C)/Fe] and [Ba/Fe] ratios, which implies they have gained C- and Ba-rich materials from extra sources, presumably asymptotic giant branch (AGB) companions. The results support the previous suggestion that these RGB stars are products of binary evolution. Their high masses thus mimic "young" single stars, yet in fact they belong to an intrinsic old stellar population. To fully explain the stellar abundance patterns of our sample stars, a variety of binary evolution channels, such as, main-sequence (MS) + RGB, MS + AGB, RGB + RGB and RGB + AGB, are required, pointing to diverse formation mechanisms of these seemly rejuvenated cannibals. With this larger sample, our results confirm earlier findings that most, if not all, $\alpha$-rich stars in the Galactic disk seem to be old.

All other papers

Ingrid Pelisoli, T. R. Marsh, R. P. Ashley, Pasi Hakala, A. Aungwerojwit, K. Burdge, E. Breedt, A. J. Brown, K. Chanthorn, V. S. Dhillon, M. J. Dyer, M. J. Green, P. Kerry, S. P. Littlefair, S. G. Parsons, D. I. Sahman, J. F. Wild, S. Yotthanathong

8 pages, 6 figures. Accepted for publication in MNRAS

Accreting magnetic white dwarfs offer an opportunity to understand the interplay between spin-up and spin-down torques in binary systems. Monitoring of the white dwarf spin may reveal whether the white dwarf spin is currently in a state of near-equilibrium, or of uni-directional evolution towards longer or shorter periods, reflecting the recent history of the system and providing constraints for evolutionary models. This makes the monitoring of the spin history of magnetic white dwarfs of high interest. In this paper we report the results of a campaign of follow-up optical photometry to detect and track the 39 sec white dwarf spin pulses recently discovered in Hubble Space Telescope data of the cataclysmic variable V1460 Her. We find the spin pulsations to be present in g-band photometry at a typical amplitude of 0.4%. Under favourable observing conditions, the spin signal is detectable using 2-meter class telescopes. We measured pulse-arrival times for all our observations, which allowed us to derive a precise ephemeris for the white dwarf spin. We have also derived an orbital modulation correction that can be applied to the measurements. With our limited baseline of just over four years, we detect no evidence yet for spin-up or spin-down of the white dwarf, obtaining a lower limit of |P/Pdot|> 4e7 years, which is already 4 to 8 times longer than the timescales measured in two other cataclysmic variable systems containing rapidly rotating white dwarfs, AE Aqr and AR Sco.

Subhajit Sarkar, Nikku Madhusudhan

Accepted for publication in MNRAS, 21 pages, 9 figures, 4 tables

The recently developed JWST Exoplanet Observation Simulator (JexoSim) simulates transit spectroscopic observations of exoplanets by JWST with each of its four instruments using a time-domain approach. Previously we reported the validation of JexoSim against Pandexo and instrument team simulators. In the present study, we report a substantially enhanced version, JexoSim 2.0, which improves on the original version through incorporation of new noise sources, enhanced treatment of stellar and planetary signals and instrumental effects, as well as improved user-operability and optimisations for increased speed and efficiency. A near complete set of instrument modes for exoplanet time-series observations is now included. In this paper we report the implementation of JexoSim 2.0 and assess performance metrics for JWST in end-member scenarios using the hot Jupiter HD 209458 b and the mini-Neptune K2-18 b. We show how JexoSim can be used to compare performance across the different JWST instruments, selecting an optimal combination of instrument and subarray modes, producing synthetic transmission spectra for each planet. These studies indicate that the 1.4 {\mu}m water feature detected in the atmosphere of K2-18 b using the Hubble WFC3 might be observable in just one transit observation with JWST with either NIRISS or NIRSpec. JexoSim 2.0 can be used to investigate the impact of complex noise and systematic effects on the final spectrum, plan observations and test the feasibility of novel science cases for JWST. It can also be customised for other astrophysical applications beyond exoplanet spectroscopy. JexoSim 2.0 is now available for use by the scientific community.

Francois Foucart, Philipp Moesta, Teresita Ramirez, Alex James Wright, Siva Darbha, Daniel Kasen

The determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae, and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound and what the asymptotic velocity of ejected material will be. We then show that properly accounting for both heating and cooling during r-process nucleosynthesis is important to accurately predict the mass and kinetic energy of the outflows. We also derive a model accounting for both of these effects that can be easily implemented in merger simulations. We show, however, that the detailed velocity distribution and geometry of the outflows can currently only be captured by full 3D fluid simulations of the outflows, as non-local effect ignored by the simple criteria used in merger simulations cannot be safely neglected when modeling these effects. Finally, we propose the introduction of simple source terms in the fluid equations to approximately account for heating/cooling from r-process nucleosynthesis in future seconds-long 3D simulations of merger remnants, without the explicit inclusion of out-of-nuclear statistical equilibrium reactions in the simulations.

Alberto Traina, Stefano Marchesi, Cristian Vignali, Núria Torres-Albà, Marco Ajello, Andrealuna Pizzetti, Ross Silver, Xiurui Zhao, Tonima Tasnim Ananna, Mislav Baloković, Peter Boorman, Poshak Gandhi, Roberto Gilli, Giorgio Lanzuisi

15 pages, 22 figures, Accepted for publication on ApJ

We present the joint Chandra, XMM-Newton and NuSTAR analysis of two nearby Seyfert galaxies, NGC 3081 and ESO 565-G019. These are the only two having Chandra data in a larger sample of ten low redshift ($z \le 0.05$), candidates Compton-thick Active Galactic Nuclei (AGN) selected in the 15-150 keV band with Swift-BAT that were still lacking NuSTAR data. Our spectral analysis, performed using physically-motivated models, provides an estimate of both the line-of-sight (l.o.s.) and average (N$_{H,S}$) column densities of the two torii. NGC 3081 has a Compton-thin l.o.s. column density N$_{H,z}$=[0.58-0.62] $\times 10^{24}$cm$^{-2}$, but the N$_{H,S}$, beyond the Compton-thick threshold (N$_{H,S}$=[1.41-1.78] $\times 10^{24}$cm$^{-2}$), suggests a "patchy" scenario for the distribution of the circumnuclear matter. ESO 565-G019 has both Compton-thick l.o.s. and N$_{H,S}$ column densities (N$_{H,z}>$2.31 $\times 10^{24}$cm$^{-2}$ and N$_{H,S} >$2.57 $\times 10^{24}$cm$^{-2}$, respectively). The use of physically-motivated models, coupled with the broad energy range covered by the data (0.6-70 keV and 0.6-40 keV, for NGC 3081 and ESO 565-G019, respectively) allows us to constrain the covering factor of the obscuring material, which is C$_{TOR}$=[0.63-0.82] for NGC 3081, and C$_{TOR}$=[0.39-0.65] for ESO 565-G019.

Gitika Shukla, Raghunathan Srianand, Neeraj Gupta, Patrick Petitjean, Andrew J. Baker, Jens-Kristian Krogager, Pasquier Noterdaeme

21 pages + Appendix, 18 figures, submitted to MNRAS

We use Southern African Large Telescope (SALT) to perform long-slit spectroscopic observations of 23 newly discovered radio-loud quasars (RLQs) at $2.7<z<3.3$. The sample consists of powerful AGN brighter than 200 mJy at 1.4 GHz and is selected on the basis of mid-infrared colors i.e., unbiased to the presence of dust. We report 7 confirmed and 5 tentative detections of diffuse Ly$\alpha$ emission in the sample. We present the properties of diffuse Ly$\alpha$ emission and discuss in detail its relationship to different quasar properties. We find strong dependence of Ly$\alpha$ halo detection rate on the extent of radio source, spectral luminosity of RLQ at 420 MHz ($L_{\rm 420MHz}$), presence of associated C IV absorption and nuclear He II emission line equivalent width. As seen in previous surveys, the FWHM of diffuse Ly$\alpha$ emission in the case of confirmed detections are much higher (i.e $>$1000 km/s). Using the samples of high-$z$ radio-loud quasars and galaxies from literature, we confirm the correlation between the Ly$\alpha$ halo luminosity and its size with radio power ($L_{\rm 420MHz}$). The same quantities are found to be correlating weakly with the projected linear size of the radio emission. Our sample is the second largest sample of RLQs being studied for the presence of diffuse Ly$\alpha$ emission and fills in a redshift gap between previous such studies. Integral Field Spectroscopy is required to fully understand the relationship between the large scale radio emission and the overall distribution, kinematics and over density of Ly$\alpha$ emission in the field of these RLQs.

Wonki Lee, M. James Jee, Kim HyeongHan, Ruta Kale, Hyein Yoon, William Forman, Ralph Kraft, Christine Jones, Aeree Chung

12 pages, 8 figures. Submitted to ApJ. Comments are welcome

We report a discovery of a double radio relic in the cluster merger ZwCl1447.2+2619 ($z=0.376$) with uGMRT observations at $420\rm~MHz$ and $700\rm~MHz$. The linear sizes of the northern and southern relics are $\sim0.3~$Mpc and $\sim1.2~$Mpc, respectively, which is consistent with the theoretical expectation that a larger relic is produced in the less massive subcluster side. However, ZwCl1447.2+2619 is unlike other known double radio relic systems, where the larger relics are much more luminous by several factors. In this merger the higher surface brightness of the smaller northern relic makes its total radio luminosity comparable to that of the much larger southern relic. The surface brightness ratio $\sim0.1$ between the two radio relics differs significantly from the relation observed in other double radio relic systems. From our radio spectral analysis, we find that both relics signify similar weak shocks with Mach numbers of $2.9\pm0.8$ and $2.0\pm0.7$ for the northern and southern relics, respectively. Moreover, the northern relic is connected to a discrete radio source with an optical counterpart, which indicates the possible presence of cosmic ray injection and re-acceleration. Therefore, we propose that this atypical surface brightness ratio can be explained with the particle acceleration efficiency precipitously dropping in the weak shock regime and/or with re-acceleration of fossil cosmic rays. Our multi-wavelength analysis and numerical simulation suggest that ZwCl1447.2+2619 is a post-merger, which has experienced a near head-on collision $\sim0.7\rm~Gyr$ ago.

N. Torres-Albà, S. Marchesi, X. Zhao, M. Ajello, R. Silver, T. T. Ananna, M. Baloković, P. B. Boorman, A. Comastri, R. Gilli, G. Lanzuisi, K. Murphy, C. M. Urry, C. Vignali

31 pages, 13 figures, accepted for publication in ApJ

We present the analysis of simultaneous NuSTAR and XMM-Newton data of 8 Compton-thick (CT-) active galactic nuclei (AGN) candidates selected in the Swift-Burst Alert Telescope (BAT) 100 month survey. This work is part of an ongoing effort to find and characterize all CT-AGN in the local ($z\leq$0.05) Universe. We used two physically motivated models, MYTorus and borus02, to characterize the sources in the sample, finding 5 of them to be confirmed CT-AGN. These results represent an increase of $\sim19$% over the previous NuSTAR-confirmed, BAT-selected CT-AGN at $z\leq0.05$, bringing the total number to 32. This corresponds to an observed fraction of $\sim 8$\% of all AGN within this volume-limited sample, although it increases to $20\pm5$% when limiting the sample to $z\leq0.01$. Out of a sample of 48 CT-AGN candidates, selected using BAT and soft (0.3$-$10 keV) X-ray data, only 24 are confirmed as CT-AGN with the addition of the NuSTAR data. This highlights the importance of NuSTAR when classifying local obscured AGN. We also note that most of the sources in our full sample of 48 Seyfert 2 galaxies with NuSTAR data have significantly different line-of-sight and average torus column densities, favouring a patchy torus scenario.

Ryoma Ouchi, Keiichi Maeda, Joseph P. Anderson, Ryo Sawada

19 pages, 17 figures, accepted for publication in The Astrophysical Journal

Recent works have indicated that the $^{56}$Ni masses estimated for Stripped Envelope SNe (SESNe) are systematically higher than those estimated for SNe II. Although this may suggest a distinct progenitor structure between these types of SNe, the possibility remains that this may be caused by observational bias. One important possible bias is that SESNe with low $^{56}$Ni mass are dim, and therefore they are more likely to escape detection. By investigating the distributions of the $^{56}$Ni mass and distance for the samples collected from the literature, we find that the current literature SESN sample indeed suffers from a significant observational bias, i.e., objects with low $^{56}$Ni mass - if they exist - will be missed, especially at larger distances. Note, however, that those distant objects in our sample are mostly SNe Ic-BL. We also conducted mock observations assuming that the $^{56}$Ni mass distribution for SESNe is intrinsically the same with that for SNe II. We find that the $^{56}$Ni mass distribution of the detected SESNe samples moves toward higher mass than the assumed intrinsic distribution, because of the difficulty in detecting the low-$^{56}$Ni mass SESNe. These results could explain the general trend of the higher $^{56}$Ni mass distribution (than SNe II) of SESNe found thus far in the literature. However, further finding clear examples of low-$^{56}$Ni mass SESNe ($\leq 0.01M_{\odot}$) is required to add weight to this hypothesis. Also, the objects with high $^{56}$Ni mass ($\gtrsim 0.2 M_{\odot}$) are not explained by our model, which may require an additional explanation.

Jialin Li, Andrew Skemer

7 pages, 6 figures, SPIE Conference Proceedings

Santa Cruz Array of Lenslets for Exoplanet Spectroscopy (SCALES) is an instrument being designed for direct imaging of exoplanets in the infrared with the Adaptive Optics System of the W.M. Keck Observatory. The performance of SCALES will be largely affected by thermal transmission and emission from various sources, including the adaptive optics and instrument structures. The placement of a cold stop and a Lyot stop can preserve maximal and stable throughput while limiting the emission of instrument structures such as primary mirror segment gaps, secondary structures, and spider arms. Here we propose and compare three cold stops, a circular inner mask paired with circular, hexagonal, and serrated outer masks, as well as one Lyot stop design. Taking into account the pupil nutation and mirror emissivity, we model the throughput and the background emission for all designs to optimize the dimensions of the cold stop and the Lyot stop.

M. K. Rodriguez Wimberly (1), M. C. Cooper (1), D. C. Baxter (1), M. Boylan-Kolchin (2), J. S. Bullock (1), S. P. Fillingham (3), A. P. Ji (4 and 5), L. V. Sales (6), J. D. Simon (4) ((1) University of California Irvine, (2) The University of Texas at Austin, (3) University of Washington, (4) Observatories of the Carnegie Institution for Science, (5) The University of Chicago, and (6) University of California Riverside)

Key Figures: 4 and 5; 15 pages; 8 figures; submitted to MNRAS; Comments Welcomed!

As the Milky Way and its satellite system become more entrenched in near field cosmology efforts, the need for an accurate estimate of the Milky Way dark matter halos mass is increasingly critical. With the second and early third data releases of stellar proper motions from $Gaia$, several groups calculated full 6D phase-space information for the population of Milky Way satellite galaxies. Utilizing these data in comparison to subhalo properties drawn from the Phat ELVIS simulations, we constrain the Milky Way dark matter halo mass to be $\sim 1-1.2 \times 10^{12}$ M$_{\odot}$. We find that the kinematics of subhalos drawn from more- or less-massive hosts (i.e. $> 1.2 \times 10^{12}$ M$_{\odot}$ or $< 10^{12}$ M$_{\odot}$) are inconsistent with the observed velocities of the Milky Way satellites. The preferred host halo mass for the Milky Way is largely insensitive to the exclusion of systems associated the Large Magellanic Cloud, changes in galaxy formation thresholds, and variations in observational completeness. As more Milky Way satellites are discovered, their velocities (radial, tangential, and total) plus Galactocentric distances will provide further insight into the mass of the Milky Way dark matter halo.

Celia Escamilla-Rivera, Maryi Carvajal, Cristian Zamora, Martin Hendry

20 pages, 4 figures. Tables 1, 2 and 3 contains newly calibrated high redshift GRB

Gamma-ray bursts (GRBs) detected at high redshift can be used to trace the cosmic expansion history. However, the calibration of their luminosity distances is not an easy task in comparison to Type Ia Supernovae (SNeIa). To calibrate these data, correlations between their luminosity and other observed properties of GRBs need to be identified, and we must consider the validity of our assumptions about these correlations over their entire observed redshift range. In this work, we propose a new method to calibrate GRBs as cosmological distance indicators using SNeIa observations with a completely model-independent deep learning architecture. An overview of this machine learning technique was developed in [1] to study the evolution of dark energy models at high redshift. The aim of the method developed in this work is to combine two networks: a Recurrent Neural Network (RNN) and a Bayesian Neural Network (BNN). Using this computational approach, denoted RNN+BNN, we extend the network's efficacy by adding the computation of covariance matrices to the Bayesian process. Once this is done, the SNeIa distance-redshift relation can be tested on the full GRB sample and therefore used to implement a cosmographic reconstruction of the distance-redshift relation in different regimes. Thus, our newly-trained neural network is used to constrain the parameters describing the kinematical state of the Universe via a cosmographic approach at high redshifts (up to $z\approx 10$), wherein we require a very minimal set of assumptions that do not rely on dynamical equations for any specific theory of gravity.

Ziteng Wang, David L. Kaplan, Tara Murphy, Emil Lenc, Shi Dai, Ewan Barr, Dougal Dobie, B. M. Gaensler, George Heald, James K. Leung, Andrew O'Brien, Sergio Pintaldi, Joshua Pritchard, Nanda Rea, Gregory R. Sivakoff, B. W. Stappers, Adam Stewart, E. Tremou, Yuanming Wang, Patrick A. Woudt, Andrew Zic

22 pages, 15 figures, Accepted by ApJ

We report the discovery of a highly-polarized, highly-variable, steep-spectrum radio source, ASKAP J173608.2-321635, located $\sim$4\,deg from the Galactic center in the Galactic plane. The source was detected six times between 2020 January and 2020 September as part of the Australian Square Kilometre Array Pathfinder Variables and Slow Transients (ASKAP VAST) survey at 888\,MHz. It exhibited a high degree ($\sim 25$\%) of circular polarization when it was visible. We monitored the source with the MeerKAT telescope from 2020 November to 2021 February on a 2--4 week cadence. The source was not detected with MeerKAT before 2021 February 07 when it appeared and reached a peak flux density of 5.6\,mJy. The source was still highly circularly polarized, but also showed up to 80\% linear polarization, and then faded rapidly with a timescale of one day. The rotation measure of the source varied significantly, from $-11.8\pm0.8$\,rad\,m$^{-2}$ to $-64.0\pm1.5$\,rad\,m$^{-2}$, over three days. No X-ray counterpart was found in follow-up \textit{Swift} or \textit{Chandra} observations about a week after the first MeerKAT detection, with upper limits of $\sim 5.0\times10^{31}$\,erg\,s$^{-1}$ (0.3--8\,keV, assuming a distance $\sim10$ kpc). No counterpart is seen in new or archival near-infrared observations down to $J=20.8$\,mag. We discuss possible identifications for ASKAP J173608.2-321635 including a low-mass star/substellar object with extremely low infrared luminosity, a pulsar with scatter-broadened pulses, a transient magnetar, or a Galactic Center Radio Transient: none of these fully explains the observations, which suggests that ASKAP J173608.2-321635 may represent part of a new class of objects being discovered through radio imaging surveys.

Huaizhong Gao, Dongxin Yang, Jiaxing Wen, Xutao Zheng, Ming Zeng, Jirong Cang, Weihe Zeng, Xiaofan Pan, Qimin Zhou, Yihui Liu, Hua Feng, Binbin Zhang, Zhi Zeng, Yang Tian, GRID Collaboration

13 pages, 12 figures, submitted to Experimental Astronomy

The Gamma-Ray Integrated Detectors (GRID) are a space project to monitor the transient gamma-ray sky in the multi-messenger astronomy era using multiple detectors on-board CubeSats. The second GRID detector, GRID-02, was launched in 2020. The performance of the detector, including the energy response, effective area, angular response, and temperature-bias dependence, is calibrated in the laboratory and presented here. These measurements are compared with particle tracing simulations and validate the Geant4 model that will be used for generating detector responses.

J. Yana Galarza, R. López-Valdivia, J. Meléndez, D. Lorenzo-Oliveira

9 pages, 3 figures and 3 tables

Binary stars are supposed to be chemically homogeneous, as they are born from the same molecular cloud. However, high precision chemical abundances show that some binary systems display chemical differences between the components, which could be due to planet engulfment. In this work, we determine precise fundamental parameters and chemical abundances for the binary system HIP 71726/ HIP 71737. Our results show that the pair is truly conatal, coeval and comoving. We also find that the component HIP 71726 is more metal-rich than HIP 71737 in the refractory elements such as iron, with $\Delta$[Fe/H] $= 0.11 \pm 0.01$ dex. Moreover, HIP 71726 has a lithium abundance 1.03 dex higher than HIP 71737, which is the largest difference in Li detected in twin-star binary systems with $\Delta$ $T_{\rm eff}$ $\leq$ 50 K. The ingestion of $9.8^{+2.0}_{-1.6}$ M$_{\oplus}$ of rocky material fully explains both the enhancement in refractory elements and the high Li content observed in HIP 71726, thereby reinforcing the planet engulfment scenario in some binary systems.

Jielai Zhang

This is a a thesis submitted in conformity with the requirements for the degree of Doctor of Philosophy Department of Astronomy and Astrophysics University of Toronto, 2018

The low surface brightness visible wavelength Universe below 29 mag arcsec$^{-2}$ is teeming with unexplored astrophysical phenomena. Structures fainter than this surface brightness are extremely difficult to image due to systematic errors of sky subtraction and scattered light in the atmosphere and in the telescope. In Chapter 1, I show how The Dragonfly Telephoto Array (Dragonfly for short) addresses these systematics via a combination of hardware and software and is able to image at a level of 30 mag arcsec$^{-2}$ or fainter. In Chapter 2, I describe the Dragonfly Pipeline and how it is optimized for low surface brightness imaging, how it automatically rejects problematic exposures, and its cloud-orchestration. In Chapter 3, I present a study of the outer disk of the nearby spiral galaxy NGC 2841 using Dragonfly as well as archival data in UV from the Galaxy Evolution Explorer Satellite and rest frame 21 cm data using the Very Large Array. While it is commonly accepted that gas dominates over stars in galaxy outer disks, I find that in NGC 2841, this is not the case. The stellar disk extends to five times R25, and there is more stellar than gas mass at all radii. Surprisingly there is a constant ratio of stellar to gas mass beyond 30 kpc, where the disk is also warped. I propose the most likely formation mechanism for this outer disk is co-planar satellite accretion. In Chapter 4, I present a study of thermally emitted and scattered light from dust in the optically thin regions of the Spider HI Cloud, using Dragonfly and Herschel Space Observatory data. In closing the thesis (Chapter 5), I look forward to further improvements in the Dragonfly Pipeline, a population study of the formation mechanisms of galaxy disks and to carrying out tests of dust models.

Pei Wang, Di Li, Colin J. Clark, Pablo Saz Parkinson, Xian Hou, Weiwei Zhu, Lei Qian, Youling Yue, Zhichen Pan, Zhijie Liu, Xuhong Yu, Xiaoyao Xie, Qijun Zhi, Hui Zhang, Jumei Yao, Jun Yan, Chengmin Zhang, Paul S. Ray, Matthew Kerr, David A. Smith, Peter F. Michelson, Elizabeth C. Ferrara, David J. Thompson, Zhiqiang Shen, Na Wang, FAST-FermiLAT Collaboration

7 pages, 2 figures

High sensitivity radio searches of unassociated $\gamma$-ray sources have proven to be an effective way of finding new pulsars. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during its commissioning phase, we have carried out a number of targeted deep searches of \textit{Fermi} Large Area Telescope (LAT) $\gamma$-ray sources. On Feb. 27$^{th}$, 2018 we discovered an isolated millisecond pulsar (MSP), PSR J0318+0253, coincident with the unassociated $\gamma$-ray source 3FGL J0318.1+0252. PSR J0318+0253 has a spin period of $5.19$ milliseconds, a dispersion measure (DM) of $26$ pc cm$^{-3}$ corresponding to a DM distance of about $1.3$ kpc, and a period-averaged flux density of $\sim$11 $\pm$ 2 $\mu$Jy at L-band (1.05-1.45 GHz). Among all high energy MSPs, PSR J0318+0253 is the faintest ever detected in radio bands, by a factor of at least $\sim$4 in terms of L-band fluxes. With the aid of the radio ephemeris, an analysis of 9.6 years of \textit{Fermi}-LAT data revealed that PSR J0318+0253 also displays strong $\gamma$-ray pulsations. Follow-up observations carried out by both Arecibo and FAST suggest a likely spectral turn-over around 350 MHz. This is the first result from the collaboration between FAST and the \textit{Fermi}-LAT teams as well as the first confirmed new MSP discovery by FAST, raising hopes for the detection of many more MSPs. Such discoveries will make a significant contribution to our understanding of the neutron star zoo while potentially contributing to the future detection of gravitational waves, via pulsar timing array (PTA) experiments.

Lu Wen, Ke-Yao Wu, Huan Yu, Jun Fang

8 pages, 5 figures, RAA Accepted

Investigating the energy distribution of the high-energy particles in the Crab nebula

Chun-qian Li, Jian-rong Shi, Hong-liang Yan, Jian-Ning Fu, Jia-dong Li, Yong-Hui Hou

15 pages, 13figures, accepted for publication in ApJS

The LAMOST Medium-Resolution Spectroscopic Survey (LAMOST-MRS) provides an unprecedented opportunity for detecting multi-line spectroscopic systems. Based on the method of Cross-Correlation Function (CCF) and successive derivatives, we search for spectroscopic binaries and triples and derive their radial velocities (RVs) from the LAMOST-MRS spectra. A Monte-Carlo simulation is adopted to estimate the RV uncertainties. After examining over 1.3 million LAMOST DR7 MRS blue arm spectra, we obtain 3,133 spectroscopic binary (SB) and 132 spectroscopic triple (ST) candidates, which account for 1.2% of the LAMOST-MRS stars. Over 95% of the candidates are newly discovered. It is found that all of the ST candidates are on the main sequence, while around 10% of the SB candidates may have one or two components on the red giant branch.

Guang-Xuan Lan, Jun-Jie Wei, Hou-Dun Zeng, Ye Li, Xue-Feng Wu

17 pages, 6 figures, 2 tables, accepted by MNRAS

In this work, we update and enlarge the long gamma-ray burst (GRB) sample detected by the {\it Swift} satellite. Given the incomplete sampling of the faint bursts and the low completeness in redshift measurement, we carefully select a subsample of bright {\it Swift} bursts to revisit the GRB luminosity function (LF) and redshift distribution by taking into account the probability of redshift measurement. Here we also explore two general expressions for the GRB LF, i.e., a broken power-law LF and a triple power-law LF. Our results suggest that a strong redshift evolution in luminosity (with an evolution index of $\delta=1.92^{+0.25}_{-0.37}$) or in density ($\delta=1.26^{+0.33}_{-0.34}$) is required in order to well account for the observations, independent of the assumed expression of the GRB LF. However, in a one-on-one comparison using the Akaike information criterion, the best-fitting evolution model involving the triple power-law LF is statistically preferred over the best-fitting one involving the broken power-law LF with a relative probability of $\sim94.3$\% versus $\sim5.7$\%. Extrapolating our fitting results to the flux limit of the whole {\it Swift} sample, and considering the trigger probability of {\it Swift}/Burst Alert Telescope in detail, we find that the expectations from our evolution models provide a good representation of the observed distributions of the whole sample without the need for any adjustment of the model free parameters. This further confirms the reliability of our analysis results.

Naoya Kitajima, Yuichiro Tada, Shuichiro Yokoyama, Chul-Moon Yoo

25 pages, 9 figures

In this paper, we update the peak theory for the estimation of the primordial black hole (PBH) abundance, particularly by implementing the critical behavior in the estimation of the PBH mass and employing the averaged compaction function for the PBH formation criterion to relax the profile dependence. We apply our peak theory to a specific non-Gaussian feature called the exponential tail, which is characteristic in ultra slow-roll models of inflation. With this type of non-Gaussianity, the probability of a large perturbation is not suppressed by the Gaussian factor but decays only exponentially, so the PBH abundance is expected to be much enhanced. Not only do we confirm this enhancement even compared to the case of the corresponding nonlinearity parameter $f_\mathrm{NL}=5/2$, but also we find that the resultant PBH mass spectrum has a characteristic maximal mass which is not seen in the simple Press--Schechter approach.

Florian Peißker, Andreas Eckart, Basel Ali

12 pages, 8 figures, published at ApJ

Given the increased attention towards the detection of faint stars in the Galactic center, we would like to address the detectability of S62 in its apoapsis with SINFONI (VLT) and NIRC2 (KECK) in this work. Because of the nearby stars and the chance of confusion, we are using Lucy-Richardson deconvolved images to detect S62 on its Keplerian orbit around Sgr~A* with a period of less than 10 years. We use the same dataset as for S62 to trace additionally the S-cluster star S29 at the expected position based on the orbital elements presented in this work. To verify the results of the filtering technique, we are analysing K-band continuum data of the same epoch independently observed with NIRC2/KECK. Based on the well-derived orbit of S62, we find the star in projection at the expected position in 2019.4 and 2019.5. By analyzing the SINFONI data of 2019.5, we confirm the $16.1\,\pm\,0.2$ mag for S62 that was formerly derived with NACO (VLT). We base our NACO imaging analysis on the robust data set that was previously used to investigate the Schwarzschild precision of S2. We also present a critical discussion of the elsewhere proposed linear trajectory of S62 and its disputed identification with a 19 mag star found with GRAVITY mounted at the VLT Interferometer.

Christoph Lhotka, Lei Zhou

Accepted manuscript in Communications in Nonlinear Science and Numerical Simulation Available online 1 September 2021, 106024, 32 pages, 14 figures

We investigate the dynamics of charged dust interacting with the interplanetary magnetic field in a Parker spiral type model and subject to the solar wind and Poynting-Robertson effect in the vicinity of the 1:1 mean motion resonance with planet Jupiter. We estimate the shifts of the location of the minimum libration amplitude solutions close to the location of the L4 and L5 points of the classical - gravitational - problem and provide the extension of the 'librational regimes of motion' and the width of the resonance in dependency of the nongravitational parameters related to the dust grain size and surface potential of the particles. Our study is based on numerical simulations in the framework of the spatial, elliptic restricted three-body problem and semi-analytical estimates obtained by averaging of Gauss' planetary equations of motion.

Lorenzo Nativi (1), Gavin P. Lamb (2), Stephan Rosswog (1), Christoffer Lundman (1), Grzegorz Kowal (3) ((1) Department of Astronomy and Oskar Klein Centre Stockholm University, (2) School of Physics and Astronomy University of Leicester, (3) Escola de Artes Ciências e Humanidades Universidade de São Paulo)

12 pages, 6 figures, 1 table

Jets can become collimated as they propagate through dense environments and understanding such interactions is crucial for linking physical models of the environments to observations. In this work, we use 3D special-relativistic simulations to study how jets propagate through the environment created around a neutron star merger remnant by neutrino-driven winds. We simulate four jets with two different initial structures, top-hat and Gaussian, and two luminosities. After jet breakout, we study the angular jet structures and the resulting afterglow light curves. We find that the initial angular structures are efficiently washed out during the propagation, despite the small wind mass of only $\sim 10^{-3}$ M$_\odot$. The final structures depend, however, on the jet luminosity, as less energetic jets are more strongly collimated. Although entrainment of baryons leads to only moderate outflow Lorentz factors ($\approx 40$), all simulated jets can well reproduce the afterglow observed in the aftermath of GW170817. The inferred physical parameters (e.g. inclination angle, ambient particle number density), however, vary substantially between the fits and appear to be sensitive to smaller details of the angular jet shape, indicating that observationally inferred parameters may depend sensitively on the employed jet models.

Joachim Saur, Clarissa Willmes, Christian Fischer, Alexandre Wennmacher, Lorenz Roth, Allison Youngblood, Darrell F. Strobel, Ansgar Reiners

Astronomy & Astrophysic (in press)

Context: Observations of auroral emissions are powerful means to remotely sense the space plasma environment around planetary bodies and ultracool dwarfs. Therefore successful searches and characterization of aurorae outside the Solar System will open new avenues in the area of extrasolar space physics. Aims: We aim to demonstrate that brown dwarfs are ideal objects to search for UV aurora outside the Solar System. We specifically search for UV aurora on the late-type T6.5 brown dwarf 2MASS J12373919+6526148 (in the following 2MASS J1237+6526). Methods: Introducing a parameter referred to as auroral power potential, we derive scaling models for auroral powers for rotationally driven aurora applicable to a broad range of wavelengths. We also analyze Hubble Space Telescope observations obtained with the STIS camera at near-UV, far-UV, and Ly-$\alpha$ wavelengths of 2MASS J1237+6526. Results: We show that brown dwarfs, due to their typically strong surface magnetic fields and fast rotation, can produce auroral UV powers on the order of 10$^{19}$ watt or more. Considering their negligible thermal UV emission, their potentially powerful auroral emissions make brown dwarfs ideal candidates for detecting extrasolar aurorae. We find possible emission from 2MASS J1237+6526, but cannot conclusively attribute it to the brown dwarf due to low signal-to-noise values in combination with nonsystematic trends in the background fluxes. The observations provide upper limits for the emission at various UV wavelength bands. The upper limits for the emission correspond to a UV luminosity of $\sim$1 $\times$ 10$^{19}$ watt, which lies in the range of the theoretically expected values. Conclusions: The possible auroral emission from the dwarf could be produced by a close-in companion and/or magnetospheric transport processes.

We study accretion environments of active galactic nuclei when a super-massive black hole wanders in a circum-nuclear region and passes through an interstellar medium there. It is expected that a Bondi-Hoyle-Lyttleton type accretion of the interstellar matter takes place and an accretion stream of matter trapped by the black hole gravitational field appears from a tail shock region. Since the trapped matter is likely to have a certain amount of specific angular momentum, the accretion stream eventually forms an accretion ring around the black hole. According to the recent study, the accretion ring consists of a thick envelope and a thin core, and angular momenta are transfered from the inner side facing to the black hole to the opposite side respectively in the envelope and the core. As a result, a thick accretion flow and a thick excretion flow extend from the envelope, and a thin accretion disk and a thin excretion disk do from the core. The thin excretion disk is predicted to terminate at some distance forming an excretion ring, while the thick excretion flow is considered to become a super-sonic wind flowing to the infinity. The thick excretion flow from the accretion ring is expected to interact with the accretion stream toward the accretion ring and to be collimated to bi-polar cones. These pictures provide a likely guide line to interpret the overall accretion environments suggested from observations.

Gareth D. Smith, Edward Gillen, Didier Queloz, Lynne A. Hillenbrand, Jack S. Acton, Douglas R. Alves, David R. Anderson, Daniel Bayliss, Joshua T. Briegal, Matthew R. Burleigh, Sarah L. Casewell, Laetitia Delrez, Georgina Dransfield, Elsa Ducrot, Samuel Gill, Michaël Gillon, Michael R. Goad, Maximilian N. Günther, Beth A. Henderson, James S. Jenkins, Emmanuël Jehin, Maximiliano Moyano, Catriona A. Murray, Peter P. Pedersen, Daniel Sebastian, Samantha Thompson, Rosanna H. Tilbrook, Amaury H.M.J. Triaud, Jose I. Vines, Peter J. Wheatley

21 pages, 9 figures. Accepted for publication in MNRAS

We present the discovery and characterisation of an eclipsing binary identified by the Next Generation Transit Survey in the $\sim$115 Myr old Blanco 1 open cluster. NGTS J0002-29 comprises three M dwarfs: a short-period binary and a companion in a wider orbit. This system is the first well-characterised, low-mass eclipsing binary in Blanco 1. With a low mass ratio, a tertiary companion and binary components that straddle the fully convective boundary, it is an important benchmark system, and one of only two well-characterised, low-mass eclipsing binaries at this age. We simultaneously model light curves from NGTS, TESS, SPECULOOS and SAAO, radial velocities from VLT/UVES and Keck/HIRES, and the system's spectral energy distribution. We find that the binary components travel on circular orbits around their common centre of mass in $P_{\rm orb} = 1.09800524 \pm 0.00000038$ days, and have masses $M_{\rm pri}=0.3978\pm 0.0033$ M$_{\odot}$ and $M_{\rm sec}=0.2245\pm 0.0018$ M$_{\odot}$, radii $R_{\rm pri}=0.4037\pm 0.0048$ R$_{\odot}$ and $R_{\rm sec}=0.2759\pm 0.0055$ R$_{\odot}$, and effective temperatures $T_{\rm pri}=3372\,^{+44}_{-37}$ K and $T_{\rm sec}=3231\,^{+38}_{-31}$ K. We compare these properties to the predictions of seven stellar evolution models, which typically imply an inflated primary. The system joins a list of 19 well-characterised, low-mass, sub-Gyr, stellar-mass eclipsing binaries, which constitute some of the strongest observational tests of stellar evolution theory at low masses and young ages.

Markus Janson, Thomas Henning, Sascha P. Quanz, Ruben Asensio-Torres, Lars Buchhave, Oliver Krause, Enric Palle

Published in special issue of Experimental Astronomy together with other White Papers submitted to ESA in 2019 in the context of the Voyage 2050 long-term planning of the Science Programme. This is the manuscript version; the published journal version is available with open access at: this http URL

Direct detection and characterization of Earth-like planets around Sun-like stars is a core task for evaluating the prevalence of habitability and life in the Universe. Here, we discuss a promising option for achieving this goal, which is based on placing an occulter in orbit and having it project its shadow onto the E- ELT at the surface of Earth, thus providing a sufficient contrast for imaging and taking spectra of Earth-like planets in the habitable zones of Sun-like stars. Doing so at a sensible fuel budget will require tailored orbits, an occulter with a high area-to-mass ratio, and appropriate instrumentation at the E-ELT. In this White Paper, submitted in response to the ESA Voyage 2050 Call, we outline the fundamental aspects of the concept, and the most important technical developments that will be required to develop a full mission.

Jin-Wei Wang, Xiao-Jun Bi, Peng-Fei Yin

22 pages, 8 figs, and 2 tables

As a well-motivated dark matter candidate, axions can be detected through the axion-photon resonant conversion in the magnetospheres of magnetic white dwarf stars or neutron stars. In this work, we utilize Omega Centauri, which is the largest globular cluster in the Milky Way and is suggested to be the remnant core of a dwarf galaxy, to probe the axion dark matter through radio signals that originate from all the neutron stars and magnetic white dwarf stars in it. With 100 hours of observation, the combination of SKA phase 1 and LOFAR can effectively probe the parameter space of the axion-photon coupling $g_{a\gamma}$ up to $10^{-14}\sim 10^{-15}~\text{GeV}^{-1}$ for the axion mass range of $0.1\sim 30 ~\mu\text{eV}$. Depending on the choice of neutron star evolution model, this limitation is two or three and a half orders of magnitude higher than that of the single neutron star or magnetic white dwarf.

Gerald Cecil, Alexander Y. Wagner, Joss Bland-Hawthorn, Geoffrey V. Bicknell, Dipanjan Mukherjee

Accepted to The Astrophysical Journal. 30 pages, 32 figures in main text, 17 interactive figure and control Javascript in /anc subdirectory. Please retrieve paper with full resolution figures and interactive figure at this http URL

MeerKAT radio continuum and XMM-Newton X-ray images have recently revealed a spectacular bipolar channel at the Galactic Center that spans several degrees ($\sim$0.5 kpc). An intermittent jet likely formed this channel and is consistent with earlier evidence of a sustained, Seyfert-level outburst fueled by black-hole accretion onto SgA* several Myr ago. Therefore, to trace an intermittent jet that perhaps penetrated, deflected, and percolated along multiple paths through the interstellar medium, relevant interactions are identified and quantified in archival X-ray images, Hubble Space Telescope Paschen-$\alpha$ images and ALMA mm-wave spectra, and new SOAR telescope IR spectra. Hydrodynamical simulations are used to show how a currently weak jet can explain these structures and inflate the ROSAT/eROSITA X-ray and Fermi $\gamma$-ray bubbles that extend $\pm$60 deg from the Galactic plane. Thus, our Galactic outflow has features in common with energetic, jet-driven structures in the prototypical Seyfert galaxy NGC 1068.

I. Baraffe, J. Pratt, D. G. Vlaykov, T. Guillet, T. Goffrey, A. Le Saux, T. Constantino

13 pages, 14 figures; accepted for publication in A&A

We perform 2D, fully compressible, time-implicit simulations of convection in a solar-like model with the MUSIC code. Our main motivation is to explore the impact of a common tactic adopted in numerical simulations of convection that use realistic stellar conditions. This tactic is to artificially increase the luminosity and to modify the thermal diffusivity of the reference stellar model. This work focuses on the impact of these modifications on convective penetration (or overshooting) at the base of the convective envelope of a solar-like model. We explore a range of enhancement factors for the energy input and confirm the increase in the characteristic overshooting depth with the increase in the energy input. Our results highlight the importance of the impact of penetrative downflows on the thermal background below the convective boundary. This is a result of compression and shear which induce local heating and thermal mixing. The artificial increase in the energy flux intensifies the heating process by increasing the velocities in the convective zone and at the convective boundary, revealing a subtle connection between the local heating of the thermal background and the plume dynamics. This heating also increases the efficiency of heat transport by radiation which may counterbalance further heating and helps to establish a steady state. The modification of the thermal background by penetrative plumes impacts the width of the overshooting layer. Our results suggest that an artificial modification of the radiative diffusivity in the overshooting layer, rather than only accelerating the thermal relaxation, could also alter the dynamics of the penetrating plumes and thus the width of the overshooting layer. Results from simulations with an artificial modification of the energy flux and of the thermal diffusivity should be regarded with caution if used to determine an overshooting width.

M. Benedettini, A. Traficante, L. Olmi, S. Pezzuto, A. Baldeschi, S. Molinari, D. Elia, E. Schisano, M. Merello, F. Fontani, K. L. J. Rygl, J. Brand, M. T. Beltran, R. Cesaroni, S. J. Liu, L. Testi

14 pages, 11 figures, accepted by Astronomy and Astrophysics

New-generation spectroscopic surveys of the Milky Way plane have been revealing the structure of the interstellar medium, allowing the simultaneous study of dense structures from single star-forming objects or systems to entire spiral arms. We present the catalogue of molecular clouds extracted from the $^{13}$CO(1-0) data cubes of the Forgotten Quadrant Survey, which mapped the Galactic plane in the range 220\deg<l<240\deg, and -2.5\deg<b<0\deg in $^{12}$CO(1-0) and $^{13}$CO(1-0).The catalogue contains 87 molecular clouds for which the main physical parameters such as area, mass, distance, velocity dispersion, and virial parameter were derived. These structures are overall less extended and less massive than the molecular clouds identified in the $^{12}$CO(1-0) data-set because they trace the brightest and densest part of the $^{12}$CO(1-0) clouds. Conversely, the distribution of aspect ratio, equivalent spherical radius, velocity dispersion, and virial parameter in the two catalogues are similar. The mean value of the mass surface density of molecular clouds is 87$\pm$55 M$_{\odot}$ pc$^{-2}$ and is almost constant across the galactocentric radius, indicating that this parameter, which is a proxy of star formation, is mostly affected by local conditions.In data of the Forgotten Quadrant Survey, we find a good agreement between the total mass and velocity dispersion of the clouds derived from $^{12}$CO(1-0) and $^{13}$CO(1-0). This is likely because in the surveyed portion of the Galactic plane, the H$_2$ column density is not particularly high, leading to a CO emission with a not very high optical depth. This mitigates the effects of the different line opacities between the two tracers on the derived physical parameters. This is a common feature in the outer Galaxy, but our result cannot be readily generalised to the entire Milky Way.

F. P. Pijpers, M. P. Di Mauro, R. Ventura

11 pages, 8 figures

During the past 10 years the unprecedented quality and frequency resolution of asteroseismic data provided by space photometry has revolutionized the study of red-giant stars providing us with the possibility to probe the interior of thousands of these targets. Our aim is to present an asteroseismic tool which allows to determine the total angular momentum of stars, without a priori inference of their internal rotational profile. We adopt and adapt to red giants the asteroseismic inversion technique developed by Pijpers (2003} for the case of the Sun. The method has been tested assuming different artificial sets of data, including also modes with harmonic degree l> 1. We estimate with an accuracy of 14.5% the total angular momentum of the red-giant star KIC 4448777 observed by Kepler during the first four consecutive years of operation.} Our results indicate that the measurement of the total angular momentum of red-giant stars can be determined with a fairly high precision by means of asteroseismology by using a small set of rotational splittings of only dipolar modes and that our method, based on observations of stellar pulsations, provides a powerful mean for testing and modeling transport of angular momentum in stars.

C. L. Hale, D. McConnell, A. J. M. Thomson, E. Lenc, G. H. Heald, A. W. Hotan, J. K. Leung, V. A. Moss, T. Murphy, J. Pritchard, E. M. Sadler, A. J. Stewart, M. T. Whiting

Accepted for Publication in PASA; 29 pages, 20 figures, 6 tables

The Rapid ASKAP Continuum Survey (RACS) is the first large sky survey using the Australian Square Kilometre Array Pathfinder (ASKAP), covering the sky south of +41$^\circ$ declination. With ASKAP's large, instantaneous field of view, $\sim 31$ deg$^2$, RACS observed the entire sky at a central frequency of 887.5 MHz using 903 individual pointings with 15 minute observations. This has resulted in the deepest radio survey of the full Southern sky to date at these frequencies. In this paper, we present the first Stokes I catalogue derived from the RACS survey. This catalogue was assembled from 799 tiles that could be convolved to a common resolution of 25$^{\prime\prime}$, covering a large contiguous region in the declination range $\delta=-$80$^\circ$ to +30$^\circ$. The catalogue provides an important tool for both the preparation of future ASKAP surveys and for scientific research. It consists of $\sim$2.1 million sources and excludes the $|b|<5^\circ$ region around the Galactic plane. This provides a first extragalactic catalogue with ASKAP covering the majority of the sky ($\delta<+30^\circ$). We describe the methods to obtain this catalogue from the initial RACS observations and discuss the verification of the data, to highlight its quality. Using simulations, we find this catalogue detects 95% of point sources at an integrated flux density of $\sim$5 mJy. Assuming a typical sky source distribution model, this suggests an overall 95% point source completeness at an integrated flux density $\sim$3 mJy. The catalogue will be available through the CSIRO ASKAP Science Data Archive (CASDA).

A. M. Cherepashchuk, T. S. Khruzina, A. I. Bogomazov

18 pages, 17 figures, MNRAS accepted

We modelled optical light curves of Sco~X-1 obtained by the Kepler space telescope during K2 mission. Modelling was performed for the case of the strong heating of the optical star and accretion disc by X-rays. In the considered model the optical star fully filled its Roche lobe. We investigated the inverse problem in wide ranges of values of model parameters and estimated following parameters of Sco X-1: the mass ratio of components $q=M_x/M_v=3.6$ ($3.5-3.8$), where $M_x$ and $M_v$ were masses of the neutron and optical stars correspondingly, the orbital inclination was $i=30^{\circ}$ ($25^{\circ}-34^{\circ}$). In the brackets uncertainties of parameters $q$ and $i$ were shown, they originated due to uncertainties of characteristics of the physical model of Sco X-1. The temperature of non-heated optical star was $T_2 = 2500-3050$ K, its radius was $R_2=1.25R_{\odot}=8.7\times 10^{10}$ cm, and its bolometric luminosity was $L_{bol}=(2.1-4.6)\times 10^{32}$ erg s$^{-1}$. The mass of the star was $M_v\simeq 0.4M_{\odot}$. The contribution of the X-ray heated accretion disc dominated in the total optical emission of Sco~X-1. The transition between low and high states occurred due to the increase of X-ray luminosity by a factor $2-3$.

Optical and near-infrared observations are compiled for the three gamma-ray binaries hosting Be stars: PSR B1259-63, LSI+61 303, and HESS J0632+057. The emissions from the Be disk are considered to vary according to the changes in its structure, some of which are caused by interactions with the compact object (e.g., tidal forces). Due to the high eccentricity and large orbit of these systems, the interactions -- and, hence the resultant observables -- depend on the orbital phase. To explore such variations, multi-band photometry and linear polarization were monitored for the three considered systems, using two 1.5 m-class telescopes: IRSF at the South African Astronomical Observatory and Kanata at the Higashi-Hiroshima Observatory.

Eloy Peña-Asensio, Josep M. Trigo-Rodríguez, Marco Langbroek, Albert Rimola, Antonio J. Robles

Astrodynamics Journal (in press)

On February 16, 2021, an artificial object was recorded by the Spanish Meteor Network (SPMN) moving slowly over the Mediterranean. From the astrometric measurements, we identify this event as the reentry engine burn of a SpaceX Falcon 9 launch vehicle's upper stage. To study this event in detail, we adapted the plane intersection method for near-straight meteoroid trajectories to analyze slow and curved orbits associated with artificial objects. To corroborate our results, we approximated the orbital elements for the upper stage using four pieces of "debris" cataloged by the U.S. Government Combined Space Operations Center (CSpOC). Based on these calculations, we also estimated the possible deorbit hazard zone using the MSISE90 model atmosphere. We warn of the interference that these artificial bolides might have in fireball studies. In addition, given that artificial bolides will be probably more frequent in the future, we point out the new role that ground-based detection networks can play in the monitoring of potentially hazardous artificial objects in near-Earth space and determining the strewn field of artificial space debris.

Renate Mauland, Øystein Elgarøy

15 pages, 15 figures

In Berezhiani & Khoury (2015) a Superfluid Dark Matter (SFDM) model is introduced, where dark matter condenses and forms a superfluid on galactic scales. In the superfluid state phonons interact with baryons, resulting in a behavior similar to that of Modified Newtonian Dynamics (MOND). If one assumes that the DM condensate rotates along with the galaxy, a grid of vortices should form throughout the superfluid component if the rotation is fast enough. We aim to investigate the size and impact of the vortices on surrounding baryons, and to further investigate the parameter space of the model. We also look for a possible vortex solution of the Lagrangian presented for the SFDM theory. We first take a simple approach and investigate vortex properties in a constant density DM halo, applying knowledge from condensed matter physics. We then use the zero-temperature condensate density profile as a template to vary the DM particle mass and the energy scale, $\Lambda$, of the SFDM model. Further, we attempt to find a vortex solution of the theory by extracting the Euler-Lagrange equation with respect to the modulus of the condensate wavefunction from the full relativistic SFDM Lagrangian. For the constant density approach we find that the vortices are on millimeter scale, and separated by distances $\sim0.002\,\rm{AU}$. The parameter space of the model is found to be substantial and a reduction in the DM particle mass leads to larger vortices with a higher energy. However, none of the parameter combinations explored here give both realistic values of $\Lambda$ and vortices energetic enough to have an observational impact on the galaxy as a whole. The vortex equation extracted from the Lagrangian of the model is unstable, and no solution exhibiting the standard properties of a vortex solution is found.

Theodore Kareta, Vishnu Reddy, Neil Pearson, Juan A. Sanchez, Walter M. Harris

32 pages, 5 body figures, 2 appendix figures. Accepted in the Planetary Science Journal on August 5, 2021

The relationship between the Near-Earth Objects (3200) Phaethon and (155140) 2005 UD is unclear. While both are parents to Meteor Showers, (the Geminids and Daytime Sextantids, respectively), have similar visible-wavelength reflectance spectra and orbits, dynamical investigations have failed to find any likely method to link the two objects in the recent past. Here we present the first near-infrared reflectance spectrum of 2005 UD, which shows it to be consistently linear and red-sloped unlike Phaethon's very blue and concave spectrum. Searching for a process that could alter some common starting material to both of these end states, we hypothesized that the two objects had been heated to different extents, motivated by their near-Sun orbits, the composition of Geminid meteoroids, and previous models of Phaethon's surface. We thus set about building a new laboratory apparatus to acquire reflectance spectra of meteoritic samples after heating to higher temperatures than available in the literature to test this hypothesis and were loaned a sample of the CI Chondrite Orgueil from the Vatican Meteorite Collection for testing. We find that while Phaethon's spectrum shares many similarities with different CI Chondrites, 2005 UD's does not. We thus conclude that the most likely relationship between the two objects is that their similar properties are only by coincidence as opposed to a parent-fragment scenario, though the ultimate test will be when JAXA's DESTINY+ mission visits one or both objects later this decade. We also discuss possible paths forward to understanding Phaethon's properties from dynamical and compositional grounds.

A. Anitra, T. Di Salvo, R. Iaria, L. Burderi, A.F. Gambino, S.M. Mazzola, A. Marino, A. Sanna, A. Riggio

The X-ray source 4U 1822-371 is an eclipsing low-mass X-ray binary and X-ray pulsar, hosting a NS that shows periodic pulsations in the X-ray band. The inclination angle of the system is so high that in principle, it should be hard to observe both the direct thermal emission of the central object and the reflection component of the spectrum because they are hidden by the outer edge of the accretion disc. Assuming that the source accretes at the Eddington limit, we analysed non-simultaneous XMM-Newton and NuSTAR observations and studied the average broadband spectrum, with the aim to investigate the presence of a reflection component. No such component has been observed before in a high-inclination source such as 4U 1822-371. We modelled the spectral emission of the source using two different reflection models, Diskline plus Pexriv and the self-consistent model RfxConv. In our analysis, we find significant evidence of a reflection component in the spectrum, in addition to two lines associated with neutral or mildly ionised iron. The continuum spectrum is well fitted by a saturated Comptonisation model and a thermal black-body component emitted by the accretion disc at a lower temperature. We updated the ephemeris, adding two new eclipse times to the most recent ephemeris reported in literature. In our proposed scenario, the source is accreting at the Eddington limit with an intrinsic luminosity of $10^{38}$ erg/s, while the observed luminosity is two orders of magnitude lower. Despite the high inclination, we find that a reflection component is required to fit residuals at the Fe line range and the hard excess observed in the spectrum. The best-fit value of the inner disc radius is still uncertain and model dependent. More observations are therefore needed to confirm these results, which can give important information on this enigmatic and peculiar source.

S. Alan Stern, Brian Keeney, Rachael Hoover, Silvia Protopapa, Oliver White, Will Grundy, Dale P. Cruikshank, the New Horizons Team

36 pages, 14 figures

Sputnik Planitia, Pluto's gigantic ice glacier, hosts numerous scientific mysteries, including the presence of thousands of elongated pit structures. We examine various attributes of these pit structures in New Horizons data sets, revealing their length, aspect ratios, and orientation properties; we also study their interior reflectivities, colors, and compositions, and compare these attributes to some other relevant regions on Pluto. We then comment on origin mechanisms of the pits and also the fate of the missing volatiles represented by the pits on Sputnik Planitia.

David Benisty, Eduardo I. Guendelman, Armin van de Venn, David Vasak, Jürgen Struckmeier, Horst Stoecker

6 pages, 1 figure

An extension to the Einstein-Cartan (EC) action is discussed in terms of cosmological solutions. The torsion incorporated in the EC Lagrangian is assumed to be totally anti-symmetric, and written by of a vector $S^\mu$. Then this torsion model, compliant with the Cosmological Principle, is made dynamical by introducing its quadratic, totally anti-symmetric derivative. The EC Lagrangian then splits up into the Einstein-Hilbert portion and a (mass) term $\sim s_0^2$. While for the quintessence model, dark energy arises from the potential, here the kinetic term, $\frac{1}{\mu^2} \dot{s}_0^2$, plays the role of dark energy. The quadratic torsion term, on the other hand, gives rise to a stiff fluid that leads to a bouncing solution. A bound on the bouncing solution is calculated.

S. D. Bale, T. S. Horbury, M. Velli, M. I. Desai, J. S. Halekas, M. D. McManus, O. Panasenco, S. T. Badman, T. A. Bowen, B. D. G. Chandran, J. F. Drake, J. C. Kasper, R. Laker, A. Mallet, L Matteini, T. D. Phan, N. E. Raouafi, J. Squire, L. D. Woodham, T. Wooley

18 pages, 9 figures, submitted to Astrophys. J

One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfv\'enic magnetic field reversals termed 'switchbacks'. These $\delta B_R/B \sim \mathcal{O}(1$) fluctuations occur on a range of timescales and in {\em patches} separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate that patches of switchbacks are localized within the extensions of plasma structures originating at the base of the corona. These structures are characterized by an increase in alpha particle abundance, Mach number, plasma $\beta$ and pressure, and by depletions in the magnetic field magnitude and electron temperature. These intervals are in pressure-balance, implying stationary spatial structure, and the field depressions are consistent with overexpanded flux tubes. The structures are asymmetric in Carrington longitude with a steeper leading edge and a small ($\sim$1$^\circ$) edge of hotter plasma and enhanced magnetic field fluctuations. Some structures contain suprathermal ions to $\sim$85 keV that we argue are the energetic tail of the solar wind alpha population. The structures are separated in longitude by angular scales associated with supergranulation. This suggests that these switchbacks originate near the leading edge of the diverging magnetic field funnels associated with the network magnetic field - the primary wind sources. We propose an origin of the magnetic field switchbacks, hot plasma and suprathermals, alpha particles in interchange reconnection events just above the solar transition region and our measurements represent the extended regions of a turbulent outflow exhaust.

Anne Inkenhaag, Peter G. Jonker, Giacomo Cannizzaro, Daniel Mata Sánchez, Richard D. Saxton

7 pages, 6 figures, 1 table, 1 appendix (2 additional figures, 1 additional table). Accepted for publication in MNRAS

We present results of our analysis of spectra of the host galaxies of the candidate Tidal Disruption Events (TDEs) XMMSL1 J111527.3+180638 and PTF09axc to determine the nature of these transients. We subtract the starlight component from the host galaxy spectra to determine the origin of the nuclear emission lines. Using a Baldwin-Phillips-Terlevich (BPT) diagram we conclude that the host galaxy of XMMSL1 J111527.3+180638 is classified as a Seyfert galaxy, suggesting this transient is likely to be caused by (extreme) variability in the active galactic nucleus. We find that the host of PTF09axc falls in the 'star-forming' region of the BPT-diagram, implying that the transient is a strong TDE candidate. For both galaxies we find a WISE-colour difference of $W1-W2<0.8$, which means there is no indication of a dusty torus and therefore an active galactic nucleus, seemingly contradicting our BPT finding for the host of XMMSL1 J111527.3+180638. We discuss possible reasons for the discrepant results obtained through the two methods.

Giorgio Orlando, P. Daniel Meerburg, Subodh P. Patil

36 pages, 6 figures

Cross-correlations between Cosmic Microwave Background (CMB) temperature and polarization anisotropies and $\mu$-spectral distortions have been considered to measure (squeezed) primordial scalar bispectra in a range of scales inaccessible to primary CMB bispectra. In this work we address whether it is possible to constrain tensor non-Gaussianities with these cross-correlations. We find that only primordial tensor bispectra with anisotropies leave distinct signatures, while isotropic tensor bispectra leave either vanishing or highly suppressed signatures. We discuss how the kind of anisotropies and the parity state in primordial bispectra determine the non-zero cross-correlations. By employing the so-called BipoSH formalism to capture the observational effects of these anisotropies, we make Fisher-forecasts to assess the detection prospects from $\mu T$, $\mu E$ and $\mu B$ cross-correlations. Observing anisotropies in squeezed $\langle \gamma \gamma \gamma\rangle$ and $\langle \gamma \gamma \zeta\rangle$ bispectra is going to be challenging as the imprint of tensor perturbations on $\mu$-distortions is subdominant to scalar perturbations, therefore requiring a large, independent amplification of the effect of tensor perturbations in the $\mu$-epoch. In absence of such a mechanism, anisotropies in squeezed $\langle \zeta \zeta \gamma\rangle$ bispectrum are the most relevant sources of $\mu T$, $\mu E$ and $\mu B$ cross-correlations. In particular, we point out that in models where anisotropies in $\langle \zeta \zeta \zeta \rangle$ leave potentially observable signatures in $\mu T$ and $\mu E$, the detection prospects of anisotropies in $\langle \zeta \zeta \gamma\rangle$ from $\mu B$ are enhanced.

A pictorial atlas of the spectroscopic evolution at optical wavelengths is presented for the first 18 days of the 2021 outburst of the recurrent nova RS Oph, prior to the emergence of high ionization emission lines. The spectra presented here have been obtained at daily cadence with the Asiago 1.22m + B&C (3200-7900 Ang, 2.3 \AA/pix) and Varese 0.84m + Echelle telescopes (4250-8900 \AA, resolving power 18,000). The spectra have been fully calibrated in IRAF, absolutely fluxed, and heliocentric corrected. The Echelle spectra have been also corrected for telluric absorptions.

Sergei D. Odintsov, Tanmoy Paul, Indrani Banerjee, Ratbay Myrzakulov, Soumitra SenGupta

"Physics of the Dark Universe" accepted. Abstract is reduced due to arxiv limitation

We propose a cosmological scenario in which the universe undergoes through a non-singular bounce, and after the bounce, it decelerates having a matter-like dominated evolution during some regime of the deceleration era, and finally at the present epoch it evolves through an accelerating stage. Our aim is to study such evolution in the context of Chern-Simons corrected F(R) gravity theory and confront the model with various observational data. Using the reconstruction technique, and in addition by employing suitable boundary conditions, we determine the form of F(R) for the entire possible range of the cosmic time. The form of F(R) seems to unify a non-singular bounce with a dark energy epoch, in particular, from a non-singular bounce to a deceleration epoch and from a deceleration epoch to a late time acceleration era. It is important to mention that the bouncing scenario in the present context is an asymmetric bounce, in particular, the Hubble radius monotonically increases and asymptotically diverges at the late contracting era, while it seems to decrease with time at the present epoch. Such evolution of the Hubble radius leads to the primordial perturbation modes generate at the deep contracting era when all the perturbation modes lie within the horizon. We calculate the scalar and tensor power spectra, and as a result, the primordial observables are found to be in agreement with the latest Planck 2018 constraints. In this regard, the Chern-Simons term seems to have considerable effects on the tensor perturbation evolution, however keeping intact the scalar part of the perturbation with that of in the case of a vacuum F(R) model, and as a result, the Chern-Simons term proves to play an important role in making the observable quantities consistent with the Planck results. Furthermore the theoretical expectation of the dark energy observables are confronted with the Planck+SNe+BAO data.

Rotating black holes without equatorial reflection symmetry can naturally arise in effective low-energy theories of fundamental quantum gravity, in particular, when parity-violating interactions are introduced. Adopting a theory-agnostic approach and considering a recently proposed Kerr-like black hole model, we investigate the structure and properties of accretion disk around a rotating black hole without reflection symmetry. In the absence of reflection symmetry, the accretion disk is in general a curved surface in shape, rather than a flat disk lying on the equatorial plane. Furthermore, the parameter $\epsilon$ that controls the reflection asymmetry would shrink the size of the innermost stable circular orbits, and enhance the efficiency of the black hole in converting rest-mass energy to radiation during accretion. In addition, we find that spin measurements based on the gravitational redshift observations of the disk, assuming a Kerr geometry, may overestimate the true spin values if the central object is actually a Kerr-like black hole with conspicuous equatorial reflection asymmetry.

R. Munini, E. Vannuccini, M. Boezio, P. von Doetinchem, C. Gerrity, A. Lenni, N. Marcelli, S. Quinn, F. Rogers, J.L. Ryan, A. Stoessl, M. Xiao, N. Saffold, A. Tiberio, M. Yamatani

The General AntiParticle Spectrometer (GAPS) is an Antarctic balloon-borne detector designed to measure low-energy cosmic antinuclei (< 0.25 GeV/n), with a specific focus on antideuterons, as a distinctive signal from dark matter annihilation or decay in the Galactic halo. The instrument consists of a tracker, made up of ten planes of lithium-drifted Silicon Si(Li) detectors, surrounded by a plastic scintillator Time-of-Flight system. GAPS uses a novel particle identification method based on exotic atom capture and decay with the emission of pions, protons, and atomic X-rays from a common annihilation vertex. An important ingredient for the antinuclei identification is the reconstruction of the "annihilation star" topology. A custom antinucleus annihilation reconstruction algorithm, called the "star-finding" algorithm, was developed to reconstruct the annihilation star fully, determining the annihilation vertex position and reconstructing the tracks of the primary and secondary charged particles. The reconstruction algorithm and its performances were studied on simulated data obtained with the Geant4-based GAPS simulation software, which fully reproduced the detector geometry. This custom algorithm was found to have better performance in the vertex resolution and reconstruction efficiency compared with a standard Hough-3D algorithm.

We discuss a field transformation from fields $\psi_a$ to other fields $\phi_i$ that involves derivatives, $\phi_i = \bar \phi_i(\psi_a, \partial_\alpha \psi_a, \ldots ;x^\mu)$, and derive conditions for this transformation to be invertible, primarily focusing on the simplest case that the transformation maps between a pair of two fields and involves up to their first derivatives. General field transformation of this type changes number of degrees of freedom, hence for the transformation to be invertible, it must satisfy certain degeneracy conditions so that additional degrees of freedom do not appear. Our derivation of necessary and sufficient conditions for invertible transformation is based on the method of characteristics, which is used to count the number of independent solutions of a given differential equation. As applications of the invertibility conditions, we show some non-trivial examples of the invertible field transformations with derivatives, and also give a rigorous proof that a simple extension of the disformal transformation involving a second derivative of the scalar field is not invertible.