Abstract visibility
Abstract text size

Papers for Thursday, May 05 2022

Papers with local authors

Tea Temim, Patrick Slane, John C. Raymond, Daniel Patnaude, Emily Murray, Parviz Ghavamian, Mathieu Renzo, Taylor Jacovich

19 pages, 4 tables, 11 figures, accepted for publication in ApJ

0 votes
Paper 19 — arXiv:2205.01798
0 votes
Paper 19 — arXiv:2205.01798

We present a study of the Galactic supernova remnant (SNR) G292.0+1.8, a classic example of a core-collapse SNR that contains oxygen-rich ejecta, circumstellar material, a rapidly moving pulsar, and a pulsar wind nebula (PWN). We use hydrodynamic simulations of the remnant evolution to show that the SNR reverse shock is interacting with the PWN and has most likely shocked the majority of supernova ejecta. In our models, such a scenario requires a total ejecta mass of $\lesssim 3\: \rm M_{\odot}$ and implies that there is no significant quantity of cold ejecta in the interior of the reverse shock. In light of these results, we compare the estimated elemental masses and abundance ratios in the reverse-shocked ejecta to nucleosynthesis models and find that they are consistent with a progenitor star with an initial mass of 12-16 $\: \rm M_{\odot}$. We conclude that the progenitor of G292.0+1.8 was likely a relatively low mass star that experienced significant mass loss through a binary interaction and would have produced a stripped-envelope supernova explosion. We also argue that the region known as the "spur" in G292.0+1.8 arises as a result of the pulsar's motion through the supernova ejecta and that its dynamical properties may suggest a line-of-sight component to the pulsar's velocity, leading to a total space velocity of $\sim 600\: \rm km\:s^{-1}$ and implying a significant natal kick. Finally, we discuss binary mass loss scenarios relevant to G292.0+1.8 and their implications for the binary companion properties and future searches.

All other papers

Filipp Romanov

5 pages, 6 figures, 2 tables. Published in the bimonthly Journal of the International Meteor Organization, Vol. 49, No. 6, p. 158

I present the results of my observations (visual and photographic) of the Geminid meteor shower in 2016, 2018, 2019 and 2020. I observed meteors from different regions (Moscow and Primorsky Krai) of Russia, under different observation conditions: light pollution, Moon phases and weather. I used a DSLR camera with a lens to photograph meteor tracks. I compare the results of my visual observations in different years and determine the coordinates of the meteors from the photographs to graphically demonstrate the radiant.

C. Donoso-Oliva, I. Becker, P. Protopapas, G. Cabrera-Vives, Vishnu M., Harsh Vardhan

Taking inspiration from natural language embeddings, we present ASTROMER, a transformer-based model to create representations of light curves. ASTROMER was trained on millions of MACHO R-band samples, and it can be easily fine-tuned to match specific domains associated with downstream tasks. As an example, this paper shows the benefits of using pre-trained representations to classify variable stars. In addition, we provide a python library including all functionalities employed in this work. Our library includes the pre-trained models that can be used to enhance the performance of deep learning models, decreasing computational resources while achieving state-of-the-art results.

Sophia Lilleengen, Michael S. Petersen, Denis Erkal, Jorge Peñarrubia, Sergey E. Koposov, Ting S. Li, Lara R. Cullinane, Alexander P. Ji, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Andrew B. Pace, Nora Shipp, Daniel B. Zucker, Joss Bland-Hawthorn, Tariq Hilmi

16 pages, 10 figures. Submitted to MNRAS. Comments welcome

It has recently been shown that the Large Magellanic Cloud (LMC) has a substantial effect on the Milky Way's stellar halo and stellar streams. Here, we explore how deformations of the Milky Way and LMC's dark matter haloes affect stellar streams, and whether these effects are observable. In particular, we focus on the Orphan-Chenab (OC) stream which passes particularly close to the LMC, and spans a large portion of the Milky Way's halo. We represent the Milky Way--LMC system using basis function expansions that capture their evolution in an $N$-body simulation. We present the properties of this system, such as the evolution of the densities and force fields of each galaxy. The OC stream is evolved in this time-dependent, deforming potential, and we investigate the effects of the various moments of the Milky Way and the LMC. We find that the simulated OC stream is strongly influenced by the deformations of both the Milky Way and the LMC, and that this effect is much larger than current observational errors. In particular, the Milky Way dipole has the biggest impact on the stream, followed by the evolution of the LMC's monopole, and the LMC's quadrupole. Detecting these effects would confirm a key prediction of collisionless, cold dark matter, and would be a powerful test of alternative dark matter and alternative gravity models.

Savvas Constantinou, Nikku Madhusudhan

Accepted for publication in MNRAS, 19 pages, 12 figures

The upcoming James Webb Space Telescope (JWST) promises a generational shift in the study of temperate mini-Neptune atmospheres using transit spectroscopy. High-altitude clouds however threaten to impede their atmospheric characterisation by muting spectral features. In this study, we systematically investigate JWST instrument configurations for characterising cloudy mini-Neptune atmospheres, assessing the importance of instrument choice and wavelength coverage, focusing on NIRISS and NIRSpec. We consider two temperate mini-Neptunes orbiting nearby M dwarfs, K2-18 b and TOI-732 c, with equilibrium temperatures below 400 K, as case studies and assess observations using different instrument configurations with one transit per instrument. We find that their JWST transmission spectra with modest observing time and adequate wavelength coverage can provide precise abundance constraints of key molecules H2O, CH4, and NH3 even in the presence of clouds at significantly high altitudes. The best constraints are obtained by combining all three high-resolution NIRSpec gratings (G140H+G235H+G395H) that together span the ~1-5 $\mu$m range. Single-transit observations with this three-instrument configuration allow precise abundance constraints for cloud-top pressures as low as 3 mbar and 0.1 mbar for K2-18 b and TOI-732 c, respectively, assuming a nominal 10x solar metallicity. The constraints vary with instrument combinations. We find that NIRSpec G235H+G395H is the optimal two-instrument configuration, while NIRISS or NIRSpec G235H is optimal for single-instrument observations. Absent high-altitude clouds, even single-instrument observations can provide good abundance constraints for these planets. Our findings underscore the promise of JWST transmission spectroscopy for characterising temperate mini-Neptunes orbiting nearby M dwarfs.

Ore Gottlieb, Serena Moseley, Teresita Ramirez-Aguilar, Ariadna Murguia-Berthier, Matthew Liska, Alexander Tchekhovskoy

For movies of the simulation, see this http URL

Short $\gamma$-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over 5 orders of magnitude in distance before reaching the photosphere. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if: (i) the bound ejecta isotropic equivalent mass along the pole at the time of the BH formation is $ M_p\lesssim10^{-4}~M_\odot $, setting a limit on the delay time between the merger and BH formation, otherwise the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) post-merger remnant disk contains strong large-scale vertical magnetic field, $\gtrsim10^{15}$ G; (iii) if the jet are weak ($\lesssim10^{50}$ erg), $ M_p $ must be small ($\lesssim10^{-2}~M_\odot$). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they reach the photosphere at $\sim10^{10.5}$ cm while retaining significant magnetization ($\lesssim1$), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope $\gtrsim3$), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope $\sim1.5$).

Floor S. Broekgaarden, Simon Stevenson, Eric Thrane

All code and data publicly available at this https URL

The spins of merging binary black holes offer insights into their formation history. Recently it has been argued that in isolated binary evolution of two massive stars the firstborn black hole is slowly rotating, whilst the progenitor of the second-born black hole can be tidally spun up if the binary is tight enough. Naively, one might therefore expect that only the less massive black hole in merging binaries exhibits non-negligible spin. However, if the mass ratio of the binary is "reversed" (typically during the first mass transfer episode), it is possible for the tidally spun up second-born to become the more massive black hole. We study the properties of such mass-ratio reversed (MRR) binary black hole mergers using a large set of 560 population synthesis models. We find that the more massive black hole is formed second in $\gtrsim 70\%$ of binary black holes observable by LIGO, Virgo, and KAGRA for most model variations we consider, with typical total masses $\gtrsim 20$ M$_{\odot}$ and mass ratios $q = m_2 / m_1 \sim 0.7$ (where $m_1 > m_2$). The formation history of these systems typically involves only stable mass transfer episodes. The second-born black hole has non-negligible spin ($\chi > 0.05$) in up to $25\%$ of binary black holes, with among those the more (less) massive black hole spinning in $0\%$--$80\%$ ($20\%$--$100\%$) of cases, varying greatly in our models. We discuss our models in the context of several observed gravitational-wave events and the observed mass ratio - effective spin correlation.

Michael G. Jones, David J. Sand, Michele Bellazzini, Kristine Spekkens, Ananthan Karunakaran, Elizabeth A. K. Adams, Giuseppina Battaglia, Giacomo Beccari, Paul Bennet, John M. Cannon, Giovanni Cresci, Denija Crnojevic, Nelson Caldwell, Jackson Fuson, Puragra Guhathakurta, Martha P. Haynes, John L. Inoue, Laura Magrini, Ricardo R. Munoz, Burcin Mutlu-Pakdil, Anil Seth, Jay Strader, Elisa Toloba, Dennis Zaritsky

Submitted to ApJ

We discuss five blue stellar systems in the direction of the Virgo cluster, analogous to the enigmatic object SECCO 1 (AGC 226067). These objects were identified based on their optical and UV morphology and followed up with HI observations with the VLA (and GBT), MUSE/VLT optical spectroscopy, and HST imaging. These new data indicate that one system is a distant group of galaxies. The remaining four are extremely low mass ($M_\ast \sim 10^5 \; \mathrm{M_\odot}$), are dominated by young, blue stars, have highly irregular and clumpy morphologies, are only a few kpc across, yet host an abundance of metal-rich, $12 + \log (\mathrm{O/H}) > 8.2$, HII regions. These high metallicities indicate that these stellar systems formed from gas stripped from much more massive galaxies. Despite the young age of their stellar populations, only one system is detected in HI, while the remaining three have minimal (if any) gas reservoirs. Furthermore, two systems are surprisingly isolated and have no plausible parent galaxy within $\sim$30' ($\sim$140 kpc). Although tidal stripping cannot be conclusively excluded as the formation mechanism of these objects, ram pressure stripping more naturally explains their properties, in particular their isolation, owing to the higher velocities, relative to the parent system, that can be achieved. Therefore, we posit that most of these systems formed from ram pressure stripped gas removed from new infalling cluster members, and survived in the intracluster medium long enough to become separated from their parent galaxies by hundreds of kiloparsecs, and that they thus represent a new type of stellar system.

Shyam Balaji, Guillem Domenech, Joseph Silk

25 pages, 6 figures, 2 appendices

The primordial spectrum of fluctuations may present a large peak as a result of enhancing features during inflation. This may include, but is not limited to, bumps in the inflaton's potential, phases of ultra-slow-roll or turns in multi-field space. However, in many models, inflation does not end immediately after the enhancing feature and it is likely to continue with a second phase of slow-roll. We show that the resulting induced gravitational waves may probe the primordial spectrum from the second inflationary phase, even if its amplitude is too small to directly induce detectable gravitational waves. This is because, if there are sharp peaks in the primordial spectrum, the total gravitational wave spectrum is not simply the sum of gravitational waves induced by a peaked and scale-invariant primordial spectra separately, but cross terms from interaction between these modes also become important. We also find that such cross terms always have a characteristic slope. We discuss the parameter space that may be probed by future gravitational waves detectors in the presence of these signals.

Daniel Farid, Han Aung, Daisuke Nagai, Arya Farahi, Eduardo Rozo

11 pages, 11 figures, submitted to MNRAS

We present a classification algorithm based on the Random Forest machine learning model that differentiates galaxies into orbiting, infalling and background (interloper) populations, using phase space information as input. We train and test our model with the galaxies from UniverseMachine mock catalog based on Multi-Dark Planck 2 N-body simulations. We show that we can recover the distribution of orbiting and infalling galaxies in 3D and 2D position and velocity distribution with $<1\%$ error when using probabilistic predictions in the presence of interlopers in the projected phase space. Our machine learning-based galaxy classification method enables a percent level measurements of the dynamics of cluster galaxies. We discuss potential applications of this technique to improve cluster cosmology and galaxy quenching.

Billy K. K. Lee, William R. Coulton, Leander Thiele, Shirley Ho

18 pages, 7 figures

With the advent of high-resolution, low-noise CMB measurements, the ability to extract cosmological information from thermal Sunyaev-Zel'dovich effect and kinetic Sunyaev-Zel'dovich effect will be limited not by statistical uncertainties but rather by systematic and theoretical uncertainties. The theoretical uncertainty is driven by the lack of knowledge about the electron pressure and density. Thus we explore the electron pressure and density distributions in the IllustrisTNG hydrodynamical simulations, and we demonstrate that the cluster properties exhibit a strong dependence on the halo concentration -- providing some of the first evidence of cluster assembly bias in the electron pressure and density. Further, our work shows evidence for a broken power-law mass dependence, with lower pressure in lower mass halos than previous work and a strong evolution with mass of the radial correlations in the electron density and pressure. Both of these effects highlight the differing impact of active galactic nuclei and supernova feedback on the gas in galaxy groups compared to massive clusters. We verified that we see qualitatively similar features in the SIMBA hydro-dynamical simulations, suggesting these effects could be generic features. Finally, we provide a parametric formula for the electron pressure and density profile as a function of dark matter halo mass, halo concentration, and redshift. These fitting formulae can reproduce the distribution of density and pressure of clusters and will be useful in extracting cosmological information from upcoming CMB surveys.

Chaoran Wang, S. Peng Oh, M. Ruszkowski

There is considerable evidence for widespread subsonic turbulence in galaxy clusters, most notably from {\it Hitomi}. Turbulence is often invoked to offset radiative losses in cluster cores, both by direct dissipation and by enabling turbulent heat diffusion. However, in a stratified medium, buoyancy forces oppose radial motions, making turbulence anisotropic. This can be quantified via the Froude number ${\rm Fr}$, which decreases inward in clusters as stratification increases. We exploit analogies with MHD turbulence to show that wave-turbulence interactions increase cascade times and reduces dissipation rates $\epsilon \propto {\rm Fr}$. Equivalently, for a given energy injection/dissipation rate $\epsilon$, turbulent velocities $u$ must be higher compared to Kolmogorov scalings. High resolution hydrodynamic simulations show excellent agreement with the $\epsilon \propto {\rm Fr}$ scaling, which sets in for ${\rm Fr} < 0.1$. We also compare previously predicted scalings for the turbulent diffusion coefficient $D \propto {\rm Fr}^2$ and find excellent agreement, for ${\rm Fr} < 1$. However, we find a different normalization, corresponding to stronger diffusive suppression by more than an order of magnitude. Our results imply that turbulent diffusion is more heavily suppressed by stratification, over a much wider radial range, than turbulent dissipation. Thus, the latter potentially dominates. Furthermore, this shift implies significantly higher turbulent velocities required to offset cooling, compared to previous models. These results are potentially relevant to turbulent metal diffusion (which is likewise suppressed), and to planetary atmospheres.

The improved astrometry precision of Gaia-eDR3 allows us to perform a detailed study of the Upper Scorpius OB association and revisit its kinematic, spatial, and age substructure. We achieve this by combining the Convergent Point method with clustering techniques, complementing with age estimations based on Gaia photometry. We obtain a set of 4028 candidate members for Upper Scorpius with completeness $\sim$99\% and contamination $\sim$11\%, extending the current census with Gaia by at least $\sim$9\%. We extract an astrometrically clean sample of 3114 sources with contamination $\sim$6\%. We show that Upper Scorpius can be divided into at least 3 main kinematic groups. We systematically investigate and characterize the Upper Scorpius' internal structure, revealing that at least $\sim 34\%$ of its stellar populations are contained in 7 spatial substructures, with well defined boundaries, kinematics and relative ages, with suggested names: $\pi$ Scorpii (20 $^{\pm2}_{\pm1}$ Myr), $\alpha$ Scorpii (14$^{\pm2}_{\pm1}$ Myr), $\delta$ Scorpii (9$^{\pm2}_{\pm1}$ Myr), $\beta$ Scorpii (8$^{\pm1}_{\pm1}$ Myr), $\omega$ Scorpii (8$^{\pm1}_{\pm1}$ Myr), $\nu$ Scorpii (7$^{\pm1}_{\pm1}$ Myr), after their brightest member, and the well known $\rho$ Ophiuchi (4$^{\pm1}_{\pm1}$ Myr). We find a clear correlation in (1) density-age suggesting that these substructures are expanding at a measurable rate, and (2) tangential velocity-age, providing constrains on the dynamics of these substructures, and the position of potential past triggering events. We estimate the time at which 4 potential supernovae events occurred in Upper Scorpius. Based on these findings, we suggest a star formation scenario for Upper Scorpius with unprecedented temporal resolution.

Anirudh Chiti, Anna Frebel, Alexander P. Ji, Mohammad K. Mardini, Xiaowei Ou, Joshua D. Simon, Kaitlin C. Rasmussen, Helmut Jerjen, Dongwon Kim, John E. Norris

24 pages, 7 figures, 6 tables; submitted to ApJ. Full version of Table 5 in source. Comments are welcome

We present chemical abundances and velocities of five stars between 0.3 kpc to 1.1 kpc from the center of the Tucana II ultra-faint dwarf galaxy (UFD) from high-resolution Magellan/MIKE spectroscopy. We find that every star is deficient in metals (-3.6 < [Fe/H] < -1.9) and in neutron-capture elements as is characteristic of UFD stars, unambiguously confirming their association with Tucana II. Other chemical abundances (e.g., C, iron-peak) largely follow UFD trends and suggest that faint core-collapse supernovae (SNe) dominated the early evolution of Tucana II. We see a downturn in [$\alpha$/Fe] at [Fe/H] $\approx -2.8$, indicating the onset of Type Ia SN enrichment and somewhat extended chemical evolution. The most metal-rich star has strikingly low [Sc/Fe] = $-1.29 \pm 0.48$ and [Mn/Fe] = $-1.33 \pm 0.33$, implying significant enrichment by a sub-Chandrasekhar mass Type Ia SN. We do not detect a radial velocity gradient in Tucana II ($\text{d}v_{\text{helio}}/\text{d}\theta_1=-2.6^{+3.0}_{-2.9}$ km s$^{-1}$ kpc$^{-1}$) reflecting a lack of evidence for tidal disruption, and derive a dynamical mass of $M_{1/2} (r_h) = 1.6^{+1.1}_{-0.7}\times 10^6$ M$_{\odot}$. We revisit formation scenarios of the extended component of Tucana II in light of its stellar chemical abundances. We find no evidence that Tucana II had abnormally energetic SNe, suggesting that if SNe drove in-situ stellar halo formation then other UFDs should show similar such features. Although not a unique explanation, the decline in [$\alpha$/Fe] is consistent with an early galactic merger triggering later star formation. Future observations may disentangle such formation channels of UFD outskirts.

Jesse A. Miller, Brian D. Fields

18 pages, 11 figures. Comments welcome

The widespread detection of 60Fe in geological and lunar archives provides compelling evidence for recent nearby supernova explosions within $\sim 100$ pc around 3 Myr and 7 Myr ago. The blasts from these explosions had a profound effect on the heliosphere. We perform new calculations to study the compression of the heliosphere due to a supernova blast. Assuming a steady but non-isotropic solar wind, we explore a range of properties appropriate for supernova distances inspired by recent 60Fe data, and for a 20 pc supernova proposed to account for mass extinctions at the end-Devonian period. We examine the locations of the termination shock decelerating the solar wind and the heliopause that marks the boundary between the solar wind and supernova material. Pressure balance scaling holds, consistent with studies of other astrospheres. Solar wind anisotropy does not have an appreciable effect on shock geometry. We find that supernova explosions at 50 pc (95 pc) lead to heliopause locations at 16 au (23 au) when the forward shock arrives. Thus, the outer solar system was directly exposed to the blast, but the inner planets -- including the Earth -- were not. This finding reaffirms that the delivery of supernova material to the Earth is not from the blast plasma itself, but likely is from supernova dust grains. After the arrival of the forward shock, the weakening supernova blast will lead to a gradual rebound of the heliosphere, taking $\sim100$s of kyr to expand beyond 100 au. Prospects for future work are discussed.

We introduce a novel evolutionary method that takes leverage from the MCMC method that can be used for constraining the parameters and theoretical models of Cosmology. Unlike the MCMC technique, which is essentially a non-parallel algorithm by design, the newly proposed algorithm is able to obtain the full potential of multi-core machines. With this algorithm, we could obtain the best-fit parameters of the Lambda CDM cosmological model and identify the discrepancy in the Hubble parameter $H_0$. In the present work we discuss the design principle of this novel approach and also the results from the analysis of Pantheon, OHD and Planck datasets are reported here. The estimation of parameters shows significant consistency with the previously reported values as well as a higher computational performance compared to the other similar exercises.

Celine Guervilly

submitted for publication to JGR: Planets

Stably-stratified layers may be present at the top of the electrically-conducting fluid layers of many planets either because the temperature gradient is locally subadiabatic or because a stable composition gradient is maintained by the segregation of chemical elements. Here we study the double-diffusive processes taking place in such a stable layer, considering the case of Mercury's core where the temperature gradient is stable but the composition gradient is unstable over a 800km-thick layer. The large difference in the molecular diffusivities leads to the development of buoyancy-driven instabilities that drive radial flows known as fingering convection. We model fingering convection using hydrodynamical simulations in a rotating spherical shell and varying the rotation rate and the stratification strength. For small Rayleigh numbers (i.e. weak background temperature and composition gradients), fingering convection takes the form of columnar flows aligned with the rotation axis and with an azimuthal size comparable with the layer thickness. For larger Rayleigh numbers, the flows retain a columnar structure but the azimuthal size is drastically reduced leading to thin sheet-like structures that are elongated in the meridional direction. The azimuthal length decreases when the thermal stratification increases, following closely the scaling law expected from the linear non-rotating planar theory (Stern, 1960). We find that the radial flows always remain laminar with local Reynolds number of order 1-10. Equatorially-symmetric zonal flows form due to latitudinal variations of the axisymmetric composition. The zonal velocity exceeds the non-axisymmetric velocities at the largest Rayleigh numbers. We discuss plausible implications for planetary magnetic fields.

Adam R. Kobelski, Lucas A. Tarr, Sarah A. Jaeggli, Nicholas Luber, Harry P. Warren, Sabrina L. Savage

22 pages, 15 figures, accepted by ApJS

New instruments sensitive to chromospheric radiation at X-ray, UV, Visible, IR, and sub-mm wavelengths have become available that significantly enhance our ability to understand the bi-directional flow of energy through the chromosphere. We describe the calibration, co-alignment, initial results, and public release of a new data set combining a large number of these instruments to obtain multi-wavelength photospheric, chromospheric, and coronal observations capable of improving our understanding of the connectivity between the photosphere and the corona via transient brightenings and wave signatures. The observations center on a bipolar region of enhanced network magnetic flux near disk center on SOL2017-03-17T14:00-17:00. The comprehensive data set provides one of the most complete views of chromospheric activity related to small scale brightenings in the corona and chromosphere to date. Our initial analysis shows strong spatial correspondence between the areas of broadest width of the Hydrogen-$\alpha$ spectral line and the hottest temperatures observed in ALMA Band 3 radio data, with a linear coefficient of $6.12\times 10^{-5}$\AA{}/K. The correspondence persists for the duration of co-temporal observations ($\approx 60$ minutes). Numerous transient brightenings were observed in multiple data series. We highlight a single, well observed transient brightening along a set of thin filamentary features with a duration of 20 minutes. The timing of the peak intensity transitions from the cooler (ALMA, 7000 K) to hotter (XRT, 3 MK) data series.

Charles J. Law, Sage Crystian, Richard Teague, Karin I. Öberg, Evan A. Rich, Sean M. Andrews, Jaehan Bae, Kevin Flaherty, Viviana V. Guzmán, Jane Huang, John D. Ilee, Joel H. Kastner, Ryan A. Loomis, Feng Long, Laura M. Pérez, Sebastián Pérez, Chunhua Qi, Giovanni P. Rosotti, Dary Ruíz-Rodríguez, Takashi Tsukagoshi, David J. Wilner

31 pages, 17 figures, accepted for publication in ApJ. Image cubes available at this https URL

High spatial resolution CO observations of mid-inclination (30-75{\deg}) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$.

Matthew Fong, Jiaxin Han, Jun Zhang, Xiaohu Yang, Hongyu Gao, Jiaqi Wang, Hekun Li, Antonios Katsianis, Pedro Alonso

17 pages, 13 figures, accepted: MNRAS May 03, 2022

We use weak lensing observations to make the first measurement of the characteristic depletion radius, one of the three radii that characterize the region where matter is being depleted by growing haloes. The lenses are taken from the halo catalog produced by the extended halo-based group/cluster finder applied to DESI Legacy Imaging Surveys DR9, while the sources are extracted from the DECaLS DR8 imaging data with the Fourier_Quad pipeline. We study halo masses $12 < \log ( M_{\rm grp} ~[{\rm M_{\odot}}/h] ) \leq 15.3$ within redshifts $0.2 \leq z \leq 0.3$. The virial and splashback radii are also measured and used to test the original findings on the depletion region. When binning haloes by mass, we find consistency between most of our measurements and predictions from the CosmicGrowth simulation, with exceptions to the lowest mass bins. The characteristic depletion radius is found to be roughly $2.5$ times the virial radius and $1.7 - 3$ times the splashback radius, in line with an approximately universal outer density profile, and the average enclosed density within the characteristic depletion radius is found to be roughly $29$ times the mean matter density of the Universe in our sample. When binning haloes by both mass and a proxy for halo concentration, we do not detect a significant variation of the depletion radius with concentration, on which the simulation prediction is also sensitive to the choice of concentration proxy. We also confirm that the measured splashback radius varies with concentration differently from simulation predictions.

Zhen-Yi Cai (USTC/China), Mattia Negrello (Cardiff/UK), Gianfranco De Zotti (INAF-Padova/Italy)

accepted by ApJ

The results of the recently published spectroscopically complete survey of dusty star-forming galaxies detected by the South Pole Telescope (SPT) over 2500 deg^2 proved to be challenging for galaxy formation models that generally underpredict the observed abundance of high-z galaxies. In this paper we interpret these results in the light of a physically grounded model for the evolution of spheroidal galaxies. The model accurately reproduces the measured redshift distribution of galaxies without any adjustment of the parameters. The data do not support the indications of an excess of z > 4 dusty galaxies reported by some analyses of Herschel surveys.

Petr Kabáth, Priyanka Chaturvedi, Phillip J. MacQueen, Marek Skarka, Ján Šubjak, Massimilliano Esposito, William D. Cochran, Salvatore E. Bellomo, Raine Karjalainen, Eike W. Guenther, Michael Endl, Szilárd Csizmadia, Marie Karjalainen, Artie Hatzes, Jiří Žák, Davide Gandolfi, Henri M.J. Boffin, Jose I. Vines, John H. Livingston, Rafael A. García, Savita Mathur, Lucía González-Cuesta, Martin Blažek, Douglas A. Caldwell, Knicole D. Colón, Hans Deeg, Anders Erikson, Vincent Van Eylen, William Fong, Malcolm Fridlund, Akihiko Fukui, Gábor Fűrész, Robert F. Goeke, Elisa Goffo, Steve Howell, Jon M. Jenkins, Peter Klagyivik, Judith Korth, David W. Latham, Rafael Luque, Dan Moldovan, Felipe Murgas, Norio Narita, Jaume Orell-Miquel, et al. (13 additional authors not shown)

MNRAS accepted

We present the confirmation and characterization of three hot Jupiters, TOI-1181b, TOI-1516b, and TOI-2046b, discovered by the TESS space mission. The reported hot Jupiters have orbital periods between 1.4 and 2.05 days. The masses of the three planets are $1.18\pm0.14$ M$_{\mathrm{J}}$, $3.16\pm0.12$\, M$_{\mathrm{J}}$, and 2.30 $\pm 0.28$ M$_{\mathrm{J}}$, for TOI-1181b, TOI-1516b, and TOI-2046b, respectively. The stellar host of TOI-1181b is a F9IV star, whereas TOI-1516b and TOI-2046b orbit F main sequence host stars. The ages of the first two systems are in the range of 2-5 Gyrs. However, TOI-2046 is among the few youngest known planetary systems hosting a hot Jupiter, with an age estimate of 100-400 Myrs. The main instruments used for the radial velocity follow-up of these three planets are located at Ond\v{r}ejov, Tautenburg and McDonald Observatory, and all three are mounted on 2-3 meter aperture telescopes, demonstrating that mid-aperture telescope networks can play a substantial role in the follow-up of gas giants discovered by \textit{TESS} and in the future by \textit{PLATO}.

Guillaume Payeur, Étienne Artigau, Laurence Perreault-Levasseur, René Doyon

17 pages, 16 figures, To be published in Astronomical Journal (Accepted 2022-04-21)

We present a new procedure rooted in deep learning to construct science images from data cubes collected by astronomical instruments using HxRG detectors in low-flux regimes. It improves on the drawbacks of the conventional algorithms to construct 2D images from multiple readouts by using the readout scheme of the detectors to reduce the impact of correlated readout noise. We train a convolutional recurrent neural network on simulated astrophysical scenes added to laboratory darks to estimate the flux on each pixel of science images. This method achieves a reduction of the noise on constructed science images when compared to standard flux-measurement schemes (correlated double sampling, up-the-ramp sampling), which results in a reduction of the error on the spectrum extracted from these science images. Over simulated data cubes created in a low signal-to-noise ratio regime where this method could have the largest impact, we find that the error on our constructed science images falls faster than a $1/\sqrt{N}$ decay, and that the spectrum extracted from the images has, averaged over a test set of three images, a standard error reduced by a factor of 1.85 in comparison to the standard up-the-ramp pixel sampling scheme. The code used in this project is publicly available on GitHub

Agnibha De Sarkar, Nayantara Gupta

15 pages, 4 figures, 1 table; accepted for publication in The Astrophysical Journal (ApJ)

Recent observations by the Large High Altitude Air Shower Observatory (LHAASO) have paved the way for the observational detection of PeVatrons in the Milky Way Galaxy, thus revolutionizing the field of $\gamma$-ray astrophysics. In this paper, we study one such detected source, LHAASO J1908+0621, and explore the origin of multi-TeV $\gamma$-ray emission from this source. A middle-aged radio supernova remnant SNR G40.5-0.5 and a GeV pulsar PSR J1907+0602 are co-spatial with LHAASO J1908+0621. Dense molecular clouds are also found to be associated with SNR G40.5-0.5. We explain the multi-TeV $\gamma$-ray emission observed from the direction of LHAASO J1908+0621, by the hadronic interaction between accelerated protons that escaped from the SNR shock front and cold protons present inside the dense molecular clouds, and the leptonic emission from the pulsar wind nebula (PWN) associated with the pulsar J1907+0602. Moreover, we explain lower energy $\gamma$-ray emission by considering the radiative cooling of the electrons that escaped from SNR G40.5-0.5. Finally, the combined lepto-hadronic scenario was used to explain the multi-wavelength spectral energy distribution (SED) of LHAASO J1908+0621. Although not yet significant, an ICECUBE hotspot of neutrino emission is spatially associated with LHAASO J1908+0621, indicating a possible hadronic contribution. In this paper, we show that if a hadronic component is present in LHAASO J1908+0621, then the second generation ICECUBE observatory will detect neutrino from this source.

I. Ceyhun Andaç, Benoît Cerutti, Guillaume Dubus, K. Yavuz Ekşi

7 pages, 6 figures, accepted for publication in Astronomy & Astrophysics

Pulsars show irregularities in their pulsed radio emission that originate from propagation effects and the intrinsic activity of the source. In this work, we investigate the role played by magnetic reconnection and the formation of plasmoids in the pulsar wind current sheet as a possible source of intrinsic pulse-to-pulse variability in the incoherent, high-energy emission pattern. We used a two-dimensional particle-in-cell simulation of an orthogonal pulsar magnetosphere restricted to the plane perpendicular to the star spin axis. We evolved the solution for several tens of pulsar periods to gather a statistically significant sample of synthetic pulse profiles. The formation of plasmoids leads to strong pulse-to-pulse variability in the form of multiple short, bright subpulses, which appear only on the leading edge of each main pulse. These secondary peaks of emission are dominated by the dozen plasmoids that can grow up to macroscopic scales. They emerge from the high end of the hierarchical merging process occurring along the wind current layer. The flux of the subpulses is correlated with their width in phase. Although the full-scale separation is not realistic, we argue that the simulation correctly captures the demographics and the properties of the largest plasmoids, and therefore of the brightest subpulses. The prediction of subpulses at specific pulse phases provides a new observational test of the magnetic reconnection scenario as the origin of the pulsed incoherent emission. High-time-resolution observations of the Crab pulsar in the optical range may be the most promising source to target for this purpose.

S.A. Tyul'bashev, V.S. Tyul'bashev, M.A. Kitaeva, A.I. Chernyshova, V.M. Malofeev, I.V. Chashey, V.I. Shishov, R.D. Dagkesamanskii, S.V. Klimenko, I.N. Nikitin, L.D. Nikitina

published in Astronomy Report, translated by Yandex translator with correction of scientific lexis, 9 pages, 5 figures, 1 table

Pulsars with periods more than 0.4 seconds in the declination range -9o < decj < 42o and in the right ascension range 0h < r.a.< 24h were searched in parallel with the program of interplanetary scintillations monitoring of a large number of sources with the radiotelescope LPA LPI. Four-year observations carried out at the central frequency 110.25 MHz in the band 2.5 MHz and in six frequency channels with the time resolution 100 ms were used for the search. Initial detections of candidates for pulsars were based on the Fourier power spectra averaged over the whole observations period. The candidates for pulsars were checked in more details from observations analysis with higher time-frequency resolution: 32 frequency channels with time resolution of 12.5 ms. 18 new pulsars was discovered, their main characteristics are presented.

Devika Tharakkal, Andrew P. Snodin, Graeme R. Sarson, Anvar Shukurov

19 pages, 18 figures

The spatial distribution of cosmic ray (CR) particles in the interstellar medium (ISM) is of major importance in radio astronomy, where its knowledge is essential for the interpretation of observations, and in theoretical astrophysics, where CR contribute to the structure and dynamics of the ISM. Local inhomogeneities in interstellar magnetic field strength and structure can affect the local diffusivity and ensemble dynamics of the cosmic ray particles. Magnetic traps (regions between magnetic mirrors located on the same magnetic line) can lead to especially strong and persistent features in the CR spatial distribution. Using test particle simulations, we study the spatial distribution of an ensemble of CR particles (both protons and electrons) in various magnetic field configurations, from an idealized axisymmetric trap to those that emerge in intermittent (dynamo-generated) random magnetic fields. We demonstrate that both the inhomogeneity in the CR sources and the energy losses by the CR particles can lead to persistent local inhomogeneities in the CR distribution and that the protons and electrons have different spatial distributions. Our results can have profound implications for the interpretation of the synchrotron emission from astronomical objects, and in particular its random fluctuations.

Beibei Liu, Sean N. Raymond, Seth A. Jacobson

Authors' version, 12 pages, 4 figures

The Solar System's orbital structure is thought to have been sculpted by an episode of dynamical instability among the giant planets. However, the instability trigger and timing have not been clearly established. Hydrodynamical modeling has shown that while the Sun's gaseous protoplanetary disk was present the giant planets migrated into a compact orbital configuration in a chain of resonances. Here we use dynamical simulations to show that the giant planets' instability was likely triggered by the dispersal of the gaseous disk. As the disk evaporated from the inside-out, its inner edge swept successively across and dynamically perturbed each planet's orbit in turn. The associated orbital shift caused a dynamical compression of the exterior part of the system, ultimately triggering instability. The final orbits of our simulated systems match those of the Solar System for a viable range of astrophysical parameters. The giant planet instability therefore took place as the gaseous disk dissipated, constrained by astronomical observations to be a few to ten million years after the birth of the Solar System. Terrestrial planet formation would not complete until after such an early giant planet instability; the growing terrestrial planets may even have been sculpted by its perturbations, explaining the small mass of Mars relative to Earth.

Denis Bernard, Tanmoy Chattopadhyay, Fabian Kislat, Nicolas Produit

44 pages

While the scientific potential of high-energy X-ray and gamma-ray polarimetry has long been recognized, measuring the polarization of high-energy photons is challenging. To date, there has been very few significant detections from an astrophysical source. However, recent technological developments raise the possibility that this may change in the not-too-distant future. Significant progress has been made in the development of Gamma-ray Burst (GRB) polarimeters and polarization sensitive Compton telescopes. A second-generation dedicated GRB polarimeter, POLAR-2, is under development for launch in 2024, and COSI a second-generation polarization sensitive Compton Telescope has been selected by NASA for launch in 2025. This chapter reviews basic concepts and experimental approaches of scattering polarimetry of hard X-rays to MeV {\gamma}-rays, and pair production polarimetry of higher-energy photons

Yoshiaki Sofue

Accepted for PASJ Letters. 4 pages 3 figures

We report the discovery of an expanding cylinder of HI gas of radius 1 kpc and vertical extent 800 pc by analyzing the 21-cm line survey data from the literature. The cylinder is expanding at 150 km s$^{-1}$ and rotating at 100 km s$^{-1}$, and is interpreted as due to a high-velocity conical wind at $\sim 180$ km s$^{-1}$ from the Galactic Center. The total mass of the cylinder is estimated to be $\sim 8.5\times 10^5 M_\odot$ and kinetic energy $\sim 3\times 10^{53}$ ergs.

We investigate the prospects of electromagnetic follow-ups for binary neutron star (BNS) mergers, with the help of early warnings from decihertz gravitational-wave (GW) observatories, B-DECIGO and DO-Optimal. Extending our previous work in Liu et al. (2022), we not only give quick assessments of joint short $\gamma$-ray burst (sGRB) detection rates for different $\gamma$-ray satellites and BNS population models, but also elaborate on the analyses and results on multi-band kilonova detections for survey telescopes with different limiting magnitudes. During an assumed 4-year mission time for decihertz GW observatories, we find that for the goals of electromagnetic follow-ups, DO-Optimal performs better than B-DECIGO as a whole on the detection rate, and has a larger detectable distance for joint sGRB/kilonova searches. Taking the log-normal population model for BNS mergers and a $\text{one-day}$ early-warning time as an example, we discuss the accuracy in localization and timing, as well as the redshift distributions for various synergy observations with electromagnetic facilities and decihertz GW detectors. Based on our analyses, we propose a feasible "wait-for" pattern as a novel detecting mode for future multi-messenger astrophysics.

Étienne Martel, Geoffroy Lesur

25 pages, 30 figures, accepted to A&A

Protoplanetary discs (PPDs) are cold, dense and weakly ionised environments that witness the planetary formation. Among these discs, transition discs (TDs) are characterised by a wide cavity in the dust and gas distribution. Despite this lack of material, many TDs strongly accrete onto their central star, possibly indicating that a mechanism is driving fast accretion in TDs cavities. The presence of radially extended 'dead zones' in PPDs recently revived the interest in magnetised disc winds (MDWs), where accretion is driven by a large magnetic field. We propose that TDs could be subject to similar winds, explaining their fast-accreting and long-lived cavities. We present the results of the first 2.5D global numerical simulations of TDs harbouring MDWs using the PLUTO code. We impose a cavity in the gas distribution and consider a power law distribution for the large-scale magnetic field strength. We assume weakly ionised discs subject to ambipolar diffusion, as expected in this range of densities and temperatures. Our simulated TDs reach a steady state with an inner cavity and an outer 'standard' disc. The accretion rates remain approximately constant through the entire discs, reaching $10^{-7}~M_\odot/\text{yrs}^{-1}$ for typical surface density values. The MDW launched from the cavity is more magnetised with a much larger lever arm than the MDW launched from the outer disc. The material in the cavity is accreted at sonic velocities, and the cavity itself is rotating at 70% of the Keplerian velocity due to the efficient magnetic braking imposed by the MDW. Overall, our cavity matches the dynamical properties of an inner jet emitting disc (JED) and of magnetically arrested discs (MADs) in black hole physics. Finally, kinematic diagnostics (wind speeds, orbital and accretion velocities) could disentangle classical photo-evaporation from MDW models.

Benjamin Lynch, Nada Al-Haddad, Wenyuan Yu, Erika Palmerio, Noé Lugaz

25 pages, 13 figures, 3 tables. Accepted for publication in Adv. Space Res

We present a comprehensive analysis of the three-dimensional magnetic flux rope structure generated during the Lynch et al. (2019) magnetohydrodynamic (MHD) simulation of a global-scale, 360 degree-wide streamer blowout coronal mass ejection (CME) eruption. We create both fixed and moving synthetic spacecraft to generate time series of the MHD variables through different regions of the flux rope CME. Our moving spacecraft trajectories are derived from the spatial coordinates of Parker Solar Probe's past encounters 7 and 9 and future encounter 23. Each synthetic time series through the simulation flux rope ejecta is fit with three different in-situ flux rope models commonly used to characterize the large-scale, coherent magnetic field rotations observed in a significant fraction of interplanetary CMEs (ICMEs). We present each of the in-situ flux rope model fits to the simulation data and discuss the similarities and differences between the model fits and the MHD simulation's flux rope spatial orientations, field strengths and rotations, expansion profiles, and magnetic flux content. We compare in-situ model properties to those calculated with the MHD data for both classic bipolar and unipolar ICME flux rope configurations as well as more problematic profiles such as those with a significant radial component to the flux rope axis orientation or profiles obtained with large impact parameters. We find general agreement among the in-situ flux rope fitting results for the classic profiles and much more variation among results for the problematic profiles. We also examine the force-free assumption for a subset of the flux rope models and quantify properties of the Lorentz force within MHD ejecta intervals. We conclude that the in-situ flux rope models are generally a decent approximation to the field structure, but all the caveats associated with in-situ flux rope models will still apply...

Francesca Bisconti, Andrea Chiavassa

Prepared for submission to JINST, 16 pages

In the framework of the development of the SWGO experiment we have performed a detailed study of the single unit of an extensive air shower observatory based on an array of water Cherenkov detectors. Indeed, one of the possible water Cherenkov detector unit configurations for SWGO consists of tanks, and to reach a high detection efficiency and discrimination capability between gamma-ray and hadronic air showers, different tank designs are under investigation. In this study, we considered double-layer tanks with several sizes, shapes and number of photo-multiplier tubes (PMTs). Muons, electrons, and gamma-rays with energies typical of secondary particles in extensive air showers have been simulated entering the tanks with zenith angles from 0 to 60 degrees. The tank response was evaluated considering the number of photoelectrons produced by the PMTs, the detection efficiency, and the time resolution of the measurement of the first photon. This analysis allowed to compare the performance of tanks with different size, configuration of PMTs, and with circular, hexagonal and square geometry. The method used and the results will be discussed in this paper.

Gabriele Franciolini, Anshuman Maharana, Francesco Muia

50 pages, 8 figures

Light primordial black holes may comprise a dominant fraction of the dark matter in our Universe. This paper critically assesses whether planned and future gravitational wave detectors in the ultra-high-frequency band could constrain the fraction of dark matter composed of sub-solar primordial black holes. Adopting the state-of-the-art description of primordial black hole merger rates, we compare various signals with currently operating and planned detectors. As already noted in the literature, our findings confirm that detecting individual primordial black hole mergers with currently existing and operating proposals remains difficult. Current proposals involving gravitational wave to electromagnetic wave conversion in a static magnetic field and microwave cavities feature a technology gap with respect to the loudest gravitational wave signals from primordial black holes of various orders of magnitude. However, we point out that one recent proposal involving resonant LC circuits represents the best option in terms of individual merger detection prospects in the range $(1\div 100) \, \text{MHz}$. In the same frequency range, we note that alternative setups involving resonant cavities, whose concept is currently under development, might represent a promising technology to detect individual merger events. We also show that a detection of the stochastic gravitational wave background produced by unresolved binaries is possible only if the theoretical sensitivity of the proposed Gaussian beam detector is achieved. Such a detector, whose feasibility is subject to various caveats, may be able to rule-out some scenarios for asteroidal mass primordial black hole dark matter. We conclude that pursuing dedicated studies and developments of gravitational wave detectors in the ultra-high-frequency band remains motivated and may lead to novel probes on the existence of light primordial black holes.

Andrés E. Piatti

9 pages, 4 figures. Accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2111.02496

We study the outer regions of the Milky Way globular cluster NGC 7089 based on new Dark Energy Camera (DECam) observations. The resulting background cleaned stellar density profile reveals the existence of an extended envelope. We confirm previous results that cluster stars are found out to ~ 1deg from the cluster's centre, which is nearly three times the value of the most robust tidal radii estimations. We also used results from direct N-body simulations in order to compare with the observations. We found a fairly good agreement between the observed and numerically generated stellar density profiles. Because of the existence of gaps and substructures along globular cluster tidal tails, we closely examined the structure of the outer cluster region beyond the Jacobi radius. We extended the analysis to a sample of 35 globular clusters, 20 of them with observed tidal tails. We found that if the stellar density profile follows a power law ~ r**(-alpha), the alpha slope correlates with the globular cluster present mass, in the sense that, the more massive the globular cluster the smaller the alpha value. This trend is not found in globular clusters without observed tidal tails. The origin of such a phenomenon could be related, among other reasons, to the proposed so-called potential escapers or to the formation of globular clusters within dark matter mini-haloes.

A. Agarwal, Ashwani Pandey, Aykut Özdönmez, Ergün Ege, Avik Kumar Das, Volkan Karakulak

17 pages, 8 figures, 7 tables, accepted for publication in The Astrophysical Journal (ApJ)

We report the results from our study of the blazar S5 1803+784 carried out using the quasi-simultaneous $B$, $V$, $R$, and $I$ observations from May 2020 to July 2021 on 122 nights. Our observing campaign detected the historically bright optical flare during MJD 59063.5$-$MJD 59120.5. We also found the source in its brightest ($R_{mag}$= 13.617) and faintest ($R_{mag}$= 15.888) states till date. On 13 nights, covering both flaring and non-flaring periods, we searched for the intraday variability using the power-enhanced $F-$test and the nested ANOVA test. We found significant variability in 2 out of these 13 nights. However, no such variability was detected during the flaring period. From the correlation analysis, we observed that the emission in all optical bands were strongly correlated with a time lag of $\sim$ 0 days. To get insights into its dominant emission mechanisms, we generated the optical spectral energy distributions of the source on 79 nights and estimated the spectral indices by fitting the simple power law. Spectral index varied from 1.392 to 1.911 and showed significant variations with time and $R-$band magnitude. We have detected a mild bluer-when-brighter trend (BWB) during the whole monitoring period while a much stronger BWB trend during the flare. We also carried out a periodicity search using four different methods and found no significant periodicity during our observation duration. Based on the analysis during the flaring state of the source one can say that the emissions most likely originate from the jet rather than the accretion disk.

Christiana Erba, Richard Ignace

15 pages, 3 figures, accepted for publication in The Astrophysical Journal

The mass-loss rates from single massive stars are high enough to form radio photospheres at large distances from the stellar surface where the wind is optically thick to (thermal) free-free opacity. Here we calculate the far-infrared, millimeter, and radio band spectral energy distributions (SEDs) that can result from the combination of free-free processes and synchrotron emission, to explore the conditions for non-thermal SEDs. Simplifying assumptions are adopted in terms of scaling relations for the magnetic field strength and the spatial distribution of relativistic electrons. The wind is assumed to be spherically symmetric, and we consider the effect of Razin suppression on the synchrotron emission. Under these conditions, long-wavelength SEDs with synchrotron emission can be either more steep or more shallow than the canonical asymptotic power-law SED from a non-magnetic wind. When non-thermal emission is present, the resultant SED shape is generally not a power-law; however, the variation in slope can change slowly with wavelength. Consequently, over a limited range of wavelengths, the SED can masquerade as approximately a power law. While most observed non-thermal long-wavelength spectra are associated with binarity, synchroton emission can have only a mild influence on single-star SEDs, requiring finer levels of wavelength sampling for detection of the effect.

Shuo Xiao, Yan-Qiu Zhang, Zi-Pei Zhu, Shao-Lin Xiong, He Gao, Dong Xu, Shuang-Nan Zhang, Wen-Xi Peng, Xiao-Bo Li, Peng Zhang, Fang-Jun Lu, Lin Lin, Liang-Duan Liu, Zhen Zhang, Ming-Yu Ge, You-Li Tuo, Wang-Chen Xue, Shao-Yu Fu, Xing Liu, An Li, Tian-Cong Wang, Chao Zheng, Yue Wang, Shuai-Qing Jiang, Jin-Da Li, Jia-Cong Liu, Zhou-Jian Cao, Ce Cai, Qi-Bin Yi, Yi Zhao, Sheng-Lun Xie, Cheng-Kui Li, Qi Luo, Jin-Yuan Liao, Li-Ming Song, Shu Zhang, Jin-Lu Qu, Cong-Zhan Liu, Xu-Fang Li, Yu-Peng Xu, Ti-Pei Li

37 pages, 14 figures

The milestone of GW 170817-GRB 170817A-AT 2017gfo has shown that gravitational wave (GW) is produced during the merger of neutron star-neutron star/black hole and that in electromagnetic (EM) wave a gamma-ray burst (GRB) and a kilonovae (KN) are generated in sequence after the merger. Observationally, however, EM property during a merger is still unclear. Here we report a peculiar precursor in a KN-associated long GRB 211211A. The duration of the precursor is $\sim$ 0.2 s, and the waiting time between the precursor and the main emission (ME) of the burst is $\sim$ 1 s, which is about the same as the time interval between GW 170817 and GRB 170817A. Quasi-Periodic Oscillations (QPO) with frequency $\sim$22 Hz (at $>5\sigma$ significance) are found throughout the precursor, the first detection of periodic signals from any {\it bona fide} GRBs. This indicates most likely that a magnetar participated in the merger, and the precursor might be produced due to a catastrophic flare accompanying with torsional or crustal oscillations of the magnetar. The strong seed magnetic field of $\sim 10^{14-15}$ G at the surface of the magnetar may also account for the prolonged duration of GRB 211211A. However, it is a challenge to reconcile the rather short lifetime of a magnetar \cite{kaspi2017magnetars} with the rather long spiraling time of a binary neutron star system only by gravitational wave radiation before merger.

Dennis Zaritsky, Richard Donnerstein, Ananthan Karunakaran, C.E. Barbosa, Arjun Dey, Jennifer Kadowaki, Kristine Spekkens, Huanian Zhang

22 pages, 21 figures, accepted for publication in ApJS

We present a catalog of 5598 ultra-diffuse galaxy (UDG) candidates with effective radius $r_e > 5.3$ arcsec distributed throughout the southern portion of the DESI Legacy Imaging Survey covering $\sim$ 15000 deg$^2$. The catalog is most complete for physically large ($r_e > 2.5$ kpc) UDGs lying in the redshift range $1800 \lesssim cz/{\rm km\ s}^{-1} \lesssim 7000$, where the lower bound is defined by where incompleteness becomes significant for large objects on the sky and the upper bound by our minimum angular size selection criterion. Because physical size is integral to the definition of a UDG, we develop a method {of distance estimation} using existing redshift surveys. With three different galaxy samples, two of which contain UDGs with spectroscopic redshifts, we estimate that the method has a redshift accuracy of $\sim$ 75% when the method converges, although larger, more representative spectroscopic UDG samples are needed to fully understand the behavior of the method. We are able to estimate distances for 1079 of our UDG candidates (19%). Finally, to illustrate uses of the catalog, we present distance independent and dependent results. In the latter category we establish that the red sequence of UDGs lies on the extrapolation of the red sequence relation for bright ellipticals and that the environment-color relation is at least qualitatively similar to that of high surface brightness galaxies. Both of these results challenge some of the models proposed for UDG evolution.

Emma R. Beasor, Nathan Smith

17 pages, 10 figures. Accepted with minor revision to ApJ

Quiescent mass-loss during the red supergiant (RSG) phase has been shown to be far lower than prescriptions typically employed in single-star evolutionary models. Importantly, RSG winds are too weak to drive the production of Wolf-Rayets (WRs) and stripped-envelope supernovae (SE-SNe) at initial masses of roughly 20--40$M_{\odot}$. If single-stars are to make WRs and SE-SNe, this shifts the burden of mass-loss to rare dust-enshrouded RSGs (DE-RSGs), objects claimed to represent a short-lived high mass-loss phase. Here, we take a fresh look at the purported DE-RSGs. By modeling the mid-IR excesses of the full sample of RSGs in the LMC, we find that only one RSG has both a high mass-loss rate (\mdot $\ge$ 10$^{-4}$ $M_{\odot}$ yr$^{-1}$) and a high optical circumstellar dust extinction (7.92 mag). This one RSG is WOH G64, and it is the only one of the 14 originally proposed DE-RSGs that is actually dust enshrouded. The rest appear to be either normal RSGs without strong infrared-excess, or lower-mass asymptotic giant branch (AGB) stars. Only one additional object in the full catalog of RSGs (not previously identified as a DE-RSG) shows strong mid-IR excess. We conclude that if DE-RSGs do represent a pre-SN phase of enhanced \mdot\ in single-stars, it is extremely short-lived, only capable of removing $\leq$2\msun\ of material. This rules out the single-star post-RSG pathway for the production of WRs, LBVs, and SE-SN. Single-star models should not employ \mdot-prescriptions based on these extreme objects for any significant fraction of the RSG phase.

M. Bonavita, C. Fontanive, R. Gratton, K. Muzic, S. Desidera, B. Biller, A. Scholz, A. Sozzetti, V. Squicciarini

Accepter for publication in MNRAS. 20 pages, 8 figures, 8 tables. Full version of Table 2, 3 and 4 available with the source files

The last decade of direct imaging (DI) searches for sub-stellar companions has uncovered a widely diverse sample that challenges the current formation models, while highlighting the intrinsically low occurrence rate of wide companions, especially at the lower end of the mass distribution. These results clearly show how blind surveys, crucial to constrain the underlying planet and sub-stellar companion population, are not an efficient way to increase the sample of DI companions. It is therefore becoming clear that efficient target selection methods are essential to ensure a larger number of detections. We present the results of the COPAINS Survey conducted with SPHERE/VLT, searching for sub-stellar companions to stars showing significant proper motion differences (Delta mu) between different astrometric catalogues. We observed twenty-five stars and detected ten companions, including four new brown dwarfs: HIP 21152 B, HIP 29724 B, HD 60584 B and HIP 63734 B. Our results clearly demonstrates how astrometric signatures, in the past only giving access to stellar companions, can now thanks to Gaia reveal companions well in the sub-stellar regime. We also introduce FORECAST (Finley Optimised REtrieval of Companions of Accelerating STars), a tool which allows to check the agreement between position and mass of the detected companions with the measured Delta mu. Given the agreement between the values of the masses of the new sub-stellar companions from the photometry with the model-independent ones obtained with FORECAST, the results of COPAINS represent a significant increase of the number of potential benchmarks for brown dwarf and planet formation and evolution theories.

We use numerical relativity simulations to explore the conditions for a canonical scalar field $\phi$ minimally coupled to Einstein gravity to generate an extended phase of slow contraction that robustly smooths the universe for a wide range of initial conditions and then sets the conditions for a graceful exit stage. We show that to achieve robustness it suffices that the potential $V(\phi)$ is negative and $M_{\rm Pl}|V_{,\phi}/V|\gtrsim5$ during the smoothing phase. We also show that, to exit slow contraction, the potential must have a minimum. Beyond the minimum, we find no constraint on the uphill slope including the possibility of ending on a positive potential plateau or a local minimum with $V_{\rm min}>0$. Our study establishes ultralocality for a wide range of potentials as a key both to robust smoothing and to graceful exit.

We investigate the detectability of gravitational waves (GWs) lensed by a system that consists of binary black holes as lenses using time-domain numerical simulations. The gravitational lensing potential of this system is no longer static but evolves with time. When GWs from the source pass through the binary lens, their amplitudes can be modulated, which is similar to the phenomenon of amplitude modulation (AM) in radio communication. We find that even the frequency of the binary lens itself is too low to be detected by the LISA detection band, the sidebands in the spectrum of the lensed GWs due to AM can still be within the sensitive range of the detection band. Moreover, we also calculate the relative differences of SNR (mismatch) between the lensed and unlensed GWs. We find that the {\it mismatch} can be as significant as 9.18%. Since mismatch does not depend on the amplitude of wavefrom, the differences between the binary lensed and unlensed waveforms are substantial. This provides a robust way to identify the lensing event for the LISA project in the future.

Raymond T. Co, Yann Mambrini, Keith A. Olive

6 pages, 3 figures

We consider the generation of the baryon asymmetry in models with right-handed neutrinos produced through gravitational scattering of the inflaton during reheating. The right-handed neutrinos later decay and generate a lepton asymmetry, which is partially converted to a baryon asymmetry by Standard Model sphaleron processes. We find that a sufficient asymmetry can be generated for a wide range of right-handed neutrino masses and reheating temperatures. We also show that the same type of gravitational scattering produces Standard Model Higgs bosons, which can achieve inflationary reheating consistent with the production of a baryon asymmetry.

Massive particles produced during the cosmic inflation can imprint in the primordial non-Gaussianities as characteristic oscillating functions of various momentum ratios, known as cosmological collider signals. We initiate a study of the phase of the oscillating signals which can be unambiguously defined and measured. The phase can provide useful new information about the spin and the couplings of the intermediate heavy particles that cannot be obtained from the signal frequency and angular dependences alone. We also present new analytical results for full nonlocal signals from two typical 1-loop processes, enabling precise determination of the signal phase away from the squeezed limit.

The excitation of $f$-mode in a neutron star member of coalescing binaries accelerates the merger course, and thereby introduces a phase shift in the gravitational waveform. Emphasising on the tidal phase shift by rotating stars, we provide an accurate, yet economical, method to generate $f$-mode-involved, pre-merger waveforms using realistic spin-modulated $f$-mode frequencies for some viable equations of state. We find for slow-rotating stars that the dephasing effects of the dynamical tides can be uniquely, EOS-independently determined by the direct observables (chirp mass ${\cal M}$, symmetric ratio $\eta$ and the mutual tidal deformability ${\tilde \Lambda}$), while this universality is gradually lost for increasing spin. Although a high cutoff waveform frequency combined with large signal-to-noise ratio (SNR) is needed to trace the tidal dephasing if binary members rotate slowly, for binaries with fast rotating members ($\lesssim800\text{ Hz}$) the phase shift due to $f$-mode will exceed the uncertainty in the waveform phase at reasonable SNR ($\rho=25$) and cutoff frequency of $\gtrsim400\text{ Hz}$. In addition, a significant phase shift of $\gtrsim100$ rads can be found for a high cutoff frequency of $10^3\text{ Hz}$. For systems involving a rapidly-spinning star (potentially the secondary of GW190814), neglecting $f$-mode effect in the waveform templates can therefore lead to considerable systemic errors in the relevant analysis.

İzzet Sakallı, Sara Kanzi

Published in Turkish Journal of Physics. 69 pages. For original article, please follow the link: this https URL

We give a pedagogical introduction to black holes (BHs) greybody factors (GFs) and quasinormal modes (QNMs) and share the recent developments on those subjects. In this study, our primary focus will be on the bosonic and fermionic GFs and QNMs of various BH and brane geometries and reveal the fingerprints of the invisibles with the radiation spectra to be obtained by the WKB approximation and bounding the Bogoliubov coefficients (together with the Miller-Good transformation) methods. (*Due to the notification of arXiv "The Abstract field cannot be longer than 1,920 characters", the appeared Abstract is shortened. For the full Abstract, please download the Article.)

In this paper, we review the nonminimal derivative coupling (NDC) between gravity and a scalar field, in face of the constraint imposed by the detection of gravitational waves. We study in detail the relation between the asymptotic inflationary solution obtained from NDC and its relation with the constraint, by means of a dynamical system analysis. We thus show that such solutions are incompatible with the given constraints.

R.D. Strauss, Nigussie M. Giday, Ephrem B. Seba, Daniel A. Chekole, Gemechu F. Garuma, Bereket H. Kassa, Teshome Dugassa, C. Diedericks

Submitted to Advances in Space Research spacial issue: Space and Geophysical Observations and Recent Results related to the African Continent

We discuss a newly established neutron monitor station installed at the ENTOTO Observatory Research Center outside of Addis Ababa, Ethiopia. This is a version of a mini-neutron monitor, recently upgraded to detect individual neutrons and able to calculate the waiting time distribution between neutron pulses down to $\sim 1 $ $\mu$s. From the waiting time distribution we define and calculate a new quantity, the correlation ratio, as the ratio of correlated to uncorrelated neutrons measured inside the monitor. We propose that this quantity can, in future, be used as a proxy for spectral index of atmospheric particles incident on the monitor and show that this quantity has a weak pressure dependence. We believe that future measurements from the ENTOTO mini-neutron monitor will contribute towards the understanding of cosmic ray acceleration and transport in the heliosphere.

Stefano Gariazzo, Martina Gerbino, Thejs Brinckmann, Massimiliano Lattanzi, Olga Mena, Thomas Schwetz, Shouvik Roy Choudhury, Katherine Freese, Steen Hannestad, Christoph A. Ternes, Mariam Tórtola

10 pages, 4 figures

The extraction of the neutrino mass ordering is one of the major challenges in particle physics and cosmology, not only for its implications for a fundamental theory of mass generation in nature, but also for its decisive role in the scale of future neutrinoless double beta decay experimental searches. It has been recently claimed that current oscillation, beta decay and cosmological limits on the different observables describing the neutrino mass parameter space provide robust decisive Bayesian evidence in favor of the normal ordering of the neutrino mass spectrum [arXiv:2203.14247]. We further investigate these strong claims using a rich and wide phenomenology, with different sampling techniques of the neutrino parameter space. Contrary to the findings of Jimenez et al [arXiv:2203.14247], no decisive evidence for the normal mass ordering is found. Neutrino mass ordering analyses must rely on priors and parameterizations that are ordering-agnostic: robust results should be regarded as those in which the preference for the normal neutrino mass ordering is driven exclusively by the data, while we find a difference of up to a factor of 33 in the Bayes factors among the different priors and parameterizations exploited here. An ordering-agnostic prior would be represented by the case of parameterizations sampling over the two mass splittings and a mass scale, or those sampling over the individual neutrino masses via normal prior distributions only. In this regard, we show that the current significance in favor of the normal mass ordering should be taken as $2.7\sigma$ (i.e. moderate evidence), mostly driven by neutrino oscillation data.