Abstract visibility
Abstract text size

Papers for Monday, Jun 05 2023

Himanshu Tiwari, Benjamin McKinley, Cathryn M. Trott, Nithyanandan Thyagarajan

17 pages, 14 figures and 5 tables, submitted to PASA

Radio interferometers can potentially detect the sky-averaged signal from the Cosmic Dawn (CD) and the Epoch of Reionisation (EoR) by studying the Moon as a thermal block to the foreground sky. The first step is to mitigate the Earth-based RFI reflections (Earthshine) from the Moon, which significantly contaminate the FM band $\approx 88-110$ MHz, crucial to CD-EoR science. We analysed MWA phase-I data from $72-180$ MHz at $40$ kHz resolution to understand the nature of Earthshine over three observing nights. We took two approaches to correct the Earthshine component from the Moon. In the first method, we mitigated the Earthshine using the flux density of the two components from the data, while in the second method, we used simulated flux density based on an FM catalogue to mitigate the Earthshine. Using these methods, we were able to recover the expected Galactic foreground temperature of the patch of sky obscured by the Moon. We performed a joint analysis of the Galactic foregrounds and the Moon's intrinsic temperature $(T_{\rm Moon})$ while assuming that the Moon has a constant thermal temperature throughout three epochs. We found $T_{\rm Moon}$ to be at $184.40\pm{2.65}\rm ~K$ and $173.77\pm{2.48}\rm ~K$ using the first and the second methods, respectively, and the best-fit values of the Galactic spectral index $(\alpha)$ were found to be within the $5\%$ uncertainty level when compared with the global sky model. Compared with our previous work, these results improved constraints on the Galactic spectral index and the Moon's intrinsic temperature. We also simulated the Earthshine at the MWA between November-December 2023 to find suitable observing times less affected by the Earthshine. Such time windows can be used to schedule future observations of CD-EoR using the MWA.

We show that an extragalactic jet with a velocity shear gives rise to Fermi like acceleration process for photons scattering withing the shear layers of the jet. Such photons then gain energy to produce a high energy power law. These power law spectra at high energies are frequently observed in several extragalactic objects such as Gamma Ray Bursts (GRBs). We implement the model on GRBs to show that the obtained range of the photon indices are well within their observed values. The analytic results are confirmed with numerical simulations following Monte Carlo approach.

Gregory Walsh, Sarah Burke-Spolaor, T. Joseph W. Lazio

17 pages, 7 figures, 3 tables; accepted for publication in ApJ

Mrk 1018 is a nearby changing-look AGN that has oscillated between spectral Type 1.9 and Type 1 over a period of 40 years. Recently, a recoiling supermassive black hole (rSMBH) scenario has been proposed to explain the spectral and flux variability observed in this AGN. Detections of rSMBHs are important for understanding the processes by which SMBH binaries merge and how rSMBHs influence their galactic environment through feedback mechanisms. However, conclusive identification of any rSMBHs has remained elusive to date. In this paper, we present an analysis of 6.5 years of multi-frequency Very Long Baseline Array (VLBA) monitoring of Mrk 1018. We find that the radio emission is compact down to 2.4 pc, and displays flux density and spectral variability over the length of our campaign, typical of a flat spectrum radio core. We observe proper motion in RA of the radio core at -36.4 $\pm$ 8.6 $\mu$as yr$^{-1}$ (4.2$\sigma$), or $0.10c \pm 0.02c$ at the redshift of Mrk 1018. No significant proper motion is found in DEC (31.3 $\pm$ 25.1 $\mu$as yr$^{-1}$). We discuss possible physical mechanisms driving the proper motion, including a rSMBH. We conclude that the apparent velocity we observe of the VLBI radio core is too high to reconcile with theoretical predictions of rSMBH velocities and that the proper motion is most likely dominated by an unresolved, outflowing jet component. Future observations may yet reveal the true nature of Mrk 1018. However, our observations are not able to confirm it as a true rSMBH.

Andrea Afruni, Gabriele Pezzulli, Filippo Fraternali, Asger Grønnow

18 pages, 10 Figures. Submitted to MNRAS, comments welcome

Galactic halos accrete material from the intergalactic medium (IGM) and part of this accretion is expected to be in the form of cool ($T\sim10^4$ K) gas. A signature of this process could reside in the detection of a large amount of clouds in the circumgalactic medium (CGM) of low-redshift galaxies. However, whether this material is able to accrete onto the galaxies and feed their star formation or, instead, evaporates into the CGM hot phase (corona, $T\sim10^6$ K), is not yet understood. In this work, we investigate the evolution of cool CGM clouds accreted from the IGM and falling through a hot corona, similar to what is believed to surround low-redshift disc galaxies, through 3D high-resolution hydrodynamical simulations. We include the effects of gravity due to the dark matter halo, isotropic thermal conduction, radiative cooling and an ionizing UV background. We explored different values of parameters such as the halo mass, coronal mass, initial cloud velocity and strength of the thermal conduction. We find that, in all our simulations, the clouds lose the vast majority of their mass at distances larger than half of the galaxy virial radius and are completely dissolved in the corona before reaching the central galaxy. Our results indicate that cold accretion from the IGM can not feed star formation in $z \sim 0$ star-forming galaxies in halos with masses of $10^{11.9}\ M_{\odot}$ or above. This suggests that present-day massive star-forming galaxies can sustain their star formation only via the spontaneous or induced cooling of their hot corona.

Elias R. Udnæs, Tiago M.D. Pereira

10 pages, 9 figures, accepted for publication in A&A

Context. Three-dimensional non-local thermodynamical equilibrium (NLTE) radiative transfer calculations are a fundamental tool for a detailed spectral analysis in stellar atmospheres, but require vast amounts of computer power. This prevents their broader application. Aims. We undertake a first exploration of the use of 3D irregular grids in stellar atmospheres. In particular, we aim to test whether irregular grids can be used to speed up the 3D NLTE problem, in the same way as depth optimisation can lead to faster running times in 1D. Methods. We created irregular grids based on 3D Voronoi diagrams, sampling different distributions from a 3D radiation-magnetohydrodynamic Bifrost simulation. We developed a method for solving radiation on the 3D irregular grid and implemented a simple NLTE solver using $\Lambda$-iteration and statistical equilibrium. We applied this to a simplified hydrogen-like atom and studied the convergence properties and accuracy of the irregular grid methods. For reference, we compared them to a standard short-characteristics solver on a regular grid. Results. We find that our method for radiation in irregular grids gives similar results to those from regular grids, and that it is possible to obtain nearly the same results with about ten times fewer points in the irregular grid for the continuum intensity in local thermodynamical equilibrium. We find that the irregular grid can give good results for the NLTE problem, but it takes four times longer per iteration than the regular grid, and it converges in about the same number of iterations. This makes it particularly inefficient. Our formulation therefore does not lead to an improvement. We also find that the design of the irregular grid is crucial for accurate results, and find it non-trivial to design an irregular grid that can work well across a wide range of heights.

Martin Millon, Frédéric Courbin, Aymeric Galan, Dominique Sluse, Xuheng Ding, Malte Tewes, S. G. Djorgovski

Author's pre-review version, published in Nature Astronomy, 13 pages, 5 figures

The tight correlations found between the mass of the supermassive black holes (SMBH) and their host galaxy luminosity, stellar mass, and velocity dispersion are often interpreted as a sign of their co-evolution. Studying these correlations across redshift provides a powerful insight into the evolutionary path followed by the quasar and its host galaxy. While the mass of the black hole is accessible from single-epoch spectra, measuring the mass of its host galaxy is challenging as the quasar largely overshines its host. Here, we present a novel technique to probe quasar-host relations beyond the local universe with strong gravitational lensing, hence overcoming the use of stellar population models or velocity dispersion measurements, both prone to degeneracies. We study in detail one of the three known cases of strong lensing by a quasar to accurately measure the mass of its host and to infer a total lensing mass of $\log_{10}(M_{\rm Tot, h}/M_{\odot}) = 10.27^{+0.06}_{-0.07}~$ within the Einstein radius of 1.2 kpc. The lensing measurement is more precise than any other alternative techniques and compatible with the local $M_{BH}$-$M_{\star, h}$ scaling relation. The sample of such quasar-galaxy or quasar-quasar lensing systems should reach a few hundreds with Euclid and Rubin-LSST, thus enabling the application of such a method with statistically significant sample sizes.

Mark R. Lovell (University of Iceland)

9 pages, 6 figures, to be submitted to MNRAS. Contact: lovell@hi.is

The JWST mission is in the process of probing the galaxy mass function at $z>10$, when conceivably any delay in halo assembly due to the presence of a dwarf galaxy-scale power spectrum cutoff may drastically suppress the number of galaxies relative to the cold dark matter (CDM) expectation. We employ N-body simulations of CDM and warm dark matter (WDM) to explore how the difference in halo collapse time between these models scales with $z=0$ descendant halo mass. We demonstrate that collapse begins first for the most massive haloes, and the delay in collapse time between CDM and WDM haloes correlates inversely with descendant mass. We thus infer that only present-day dwarf galaxies exhibit any difference in their assembly history between CDM and WDM at $z=10$, and therefore support previous studies that have found JWST is unlikely to determine whether our Universe is better described by the CDM cosmology or the WDM cosmology without favourable lensing studies.

William Lake, Smadar Naoz, Federico Marinacci, Blakesley Burkhart, Mark Vogelsberger, Claire E. Williams, Yeou S. Chiou, Gen Chiaki, Yurina Nakazato, Naoki Yoshida

11 pages, 5 figures

The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that Supersonically Induced Gas Objects (SIGOs), which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination, could be the progenitors of early globular clusters. In order to become GCs, SIGOs must form stars relatively efficiently despite forming outside of dark matter halos. We investigate the potential for star formation in SIGOs using cosmological hydrodynamic simulations, including the aforementioned relative motions of baryons and dark matter, molecular hydrogen cooling in primordial gas clouds, and including explicit star formation. We find that SIGOs do form stars and that the nascent star clusters formed through this process are accreted by dark matter halos on short timescales (a few hundreds of Myr). Thus, SIGOs may be found as intact substructures within these halos, analogous to many present-day GCs. From this result, we conclude that SIGOs are capable of forming star clusters with similar properties to globular clusters in the early Universe and we discuss their detectablity by upcoming JWST surveys.

Qingzheng Yu, Taotao Fang, Junfeng Wang, Jianfeng Wu

12 pages, 5 figures, accepted for publication in ApJ

We report the discovery of three HI absorbers toward low-power radio active galactic nuclei (AGNs) in a pilot HI absorption survey with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Compared to past studies, FAST observations have explored lower radio powers by $\sim$0.4 dex and detected these weakest absorbers at given redshifts. By comparing the gas properties and kinematics of sources along radio powers, we aim to explore the interplay between AGN and the surrounding interstellar medium (ISM). Compared to brighter sources at similar redshifts, our observations suggest a slightly lower detection rate of HI absorption lines ($\sim$$11.5\%$) in low-power radio AGNs with $\text{log}(P_{\text{1.4 GHz}}/\text{W Hz}^{-1})=21.8-23.7$. The low-power sources with $\text{log}(P_{\text{1.4 GHz}}/\text{W Hz}^{-1})<23$ have a lower detection rate of $\sim$$6.7\%$. Due to the incompleteness of the sample, these detection rates may represent the lower limits. The selection of more extended sources and dilution by HI emission at lower redshifts may contribute to the lower detection rate of HI absorption lines. These detected absorbers present relatively narrow line widths and comparable column densities consistent with previous observations. One absorber has a symmetric profile with a large velocity offset, while the other two show asymmetric profiles that can be decomposed into multiple components, suggesting various possibilities of gas origins and kinematics. These HI absorbers may have connections with rotating disks, gas outflows, galactic gas clouds, gas fueling of the AGN, and jet-ISM interactions, which will be further investigated with the upcoming systematic survey and spatially resolved observations.

Eirini Angeloudi, Jesús Falcón-Barroso, Marc Huertas-Company, Regina Sarmiento, Annalisa Pillepich, Daniel Walo-Martín, Lukas Eisert

23 pages, 15 figures. Accepted for publication in MNRAS

Quantifying the contribution of mergers to the stellar mass of galaxies is key for constraining the mechanisms of galaxy assembly across cosmic time. However, the mapping between observable galaxy properties and merger histories is not trivial: cosmological galaxy simulations are the only tools we have for calibration. We study the robustness of a simulation-based inference of the ex-situ stellar mass fraction of nearby galaxies to different observables -- integrated and spatially-resolved -- and to different galaxy formation models -- IllustrisTNG and EAGLE -- with Machine Learning. We find that at fixed simulation, the fraction of accreted stars can be inferred with very high accuracy, with an error $\sim5$ per cent (10 per cent) from 2D integral-field spectroscopic maps (integrated quantities) throughout the considered stellar mass range. A bias (> 5 per cent) and an increase in scatter by a factor of 2 are introduced when testing with a different simulation, revealing a lack of generalization to distinct galaxy-formation models. Interestingly, upon using only stellar mass and kinematics maps in the central galactic regions for training, we find that this bias is removed and the ex-situ stellar mass fraction can be recovered in both simulations with < 15 per cent scatter, independently of the training set's origin. This opens up the door to a potential robust inference of the accretion histories of galaxies from existing Integral Field Unit surveys, such as MaNGA, covering a similar field of view (FOV) and containing spatially-resolved spectra for tens of thousands of nearby galaxies.

We have examined the distribution of the position angle (PA) of the Galactic center filaments with lengths $L > 66''$ and $ < 66''$ as well as their length distribution as a function of PA. We find bimodal PA distributions of the filaments, long and short populations of radio filaments. Our PA study shows the evidence for a distinct population of short filaments with PA close to the Galactic plane. Mainly thermal short radio filaments ($<66''$) have PAs concentrated close to the Galactic plane within $60^\circ < \rm PA <120^\circ$. Remarkably, the short filament PAs are radial with respect to the Galactic center at $l <0^\circ$, and extend in the direction toward Sgr A*. On a smaller scale, the prominent Sgr E HII complex G358.7-0.0 provides a vivid example of the nearly radial distribution of short filaments. The bimodal PA distribution suggests different origin for two distinct filament populations. We argue that alignment of the short filament population results from the ram pressure of a degree-scale outflow from Sgr A* that exceeds the internal filament pressure, and aligns them along the Galactic plane. The ram pressure is estimated to be 2$\times10^6\,$ cm$^{-3}\,$ K at a distance of 300pc, requiring biconical mass outflow rate $10^{-4}$ \msol\, yr$^{-1}$ with an opening angle of $\sim40^\circ$. This outflow aligns not only the magnetized filaments along the Galactic plane but also accelerates thermal material associated with embedded or partially embedded clouds. This places an estimate of $\sim$6 Myr as the age of the outflow.

Christopher Immer, Philip Metzger

9 pages, 8 figures. Presented at Earth & Space 2010 conference

With NASA's planned return to the moon and possibly with lunar outposts being formed, repeated landings at the same site will be necessary. Understanding rocket plume interaction with lunar and Martian surfaces is of paramount importance in order to safely land and protect hardware surrounding the landing site. This work will report on results of three small experiments intended to explore plume impingement onto lunar and Martian surfaces: Handheld Observation of Scour Holes (HOOSH), Handheld Angle of Repose Measurements of Lunar Simulants (HARMLuS), and Mars Architecture Team study (MATS). The first two experiments were performed during two sorties of reduced gravity flights. HOOSH was designed to investigate crater formation as a function of gravitational level (lunar and Martian gravity). HARMLuS was designed to measure the Angle of Failure (related to the angle of repose) at lunar and Martian gravity. Both experiments have complex findings indicative of the hysteretic behavior of granular materials, especially resulting from reduced gravity. The MATS experiment was designed to investigate the effects of regolith compaction on the granular mechanics of crater formation. In general, the granular mechanics is a much stronger function of compaction than gravitation acceleration. Crater formation is greatly enhanced at reduced gravity (resulting in much larger craters). The angle of failure of the lunar simulants increases with decreasing gravitational acceleration, and occasionally becomes infinite for some compactions at lunar gravity. The angle of failure also increases with increasing compaction. While compaction does play a role in the time development of crater formation, the asymptotic behavior is largely unaffected.

Laila A. Rahmatian, Philip T. Metzger

15 pages, 14 figures. Presented at Earth & Space 2010 conference

We have studied several field geotechnical test instruments for their applicability to lunar soil simulants and analog soils. Their performance was evaluated in a series of tests in lunar simulants JSC-1A, NU-LHT-2M, and CHENOBI each prepared in carefully controlled states of compaction through vibration on a shake table with overburden. In general, none of the instruments is adequate for a low-cohesion, frictional soil, but we find that a modified version of a shear vane tester allows us to extract several of the important soil parameters. This modified instrument may be useful for use on the lunar surface by astronauts or a robotic lander. We have also found that JSC-1A does not behave mechanically like the other lunar soil simulants, probably because its particle shapes are more rounded. Furthermore we have studied a soil material, BP-1, identified as very lunar-like at a lunar analog location. We find this material has a natural particle size distribution similar to that of lunar soil and arguably better than JSC-1A. We find that BP-1 behaves very similarly to the high fidelity lunar simulants NU-LHT-2M and CHENOBI.

John Lane, Steven Trigwell, Paul Hintze, Philip Metzger

10 pages, 10 figures. Presented at Earth & Space 2012 conference

The Apollo 12 lunar module (LM) landing near the Surveyor III spacecraft at the end of 1969 has remained the primary experimental verification of the predicted physics of plume ejecta effects from a rocket engine interacting with the surface of the moon. This was made possible by the return of the Surveyor III camera housing by the Apollo 12 astronauts, allowing detailed analysis of the composition of dust deposited by the LM plume. It was soon realized after the initial analysis of the camera housing that the LM plume tended to remove more dust than it had deposited. In the present study, coupons from the camera housing have been reexamined. In addition, plume effects recorded in landing videos from each Apollo mission have been studied for possible clues. Several likely scenarios are proposed to explain the Surveyor III dust observations. These include electrostatic levitation of the dust from the surface of the Moon as a result of periodic passing of the day-night terminator; dust blown by the Apollo 12 LM flyby while on its descent trajectory; dust ejected from the lunar surface due to gas forced into the soil by the Surveyor III rocket nozzle, based on Darcy's law; and mechanical movement of dust during the Surveyor landing. Even though an absolute answer may not be possible based on available data and theory, various computational models are employed to estimate the feasibility of each of these proposed mechanisms. Scenarios are then discussed which combine multiple mechanisms to produce results consistent with observations.

Conor Larison, Saurabh W. Jha, Lindsey A. Kwok, Yssavo Camacho-Neves

Submitted to AAS journals

We present an analysis of 102 type Ia supernovae (SNe Ia) in nearby (z < 0.1), x-ray selected galaxy clusters. This is the largest such sample to date and is based on archival data primarily from ZTF and ATLAS. We divide our SNe Ia into an inner cluster sample projected within $r_{500}$ of the cluster center and an outer cluster sample projected between $r_{500}$ and $2\,r_{500}$. We compare these to field samples of SNe Ia at similar redshifts in both quiescent and star-forming host galaxies. Based on SALT3 fits to the light curves, we find that the inner cluster SNe Ia have a higher fraction of fast-evolving objects (SALT3 $x_1 < -1$) than the outer cluster or field quiescent samples. This implies an intrinsically different population of SNe Ia occurs in inner cluster environments, beyond known correlations based on host galaxy alone. Our cluster samples show a strongly bimodal $x_1$ distribution with a fast-evolving component that dominates the inner cluster objects ($\gtrsim$ 75%) but is just a small fraction of SNe Ia in field star-forming galaxies ($\lesssim$ 10%). We do not see strong evidence for variations in the color (SALT3 $c$) distributions among the samples and find only minor differences in SN Ia standardization parameters and Hubble residuals. We suggest that the age of the stellar population drives the observed distributions, with the oldest populations nearly exclusively producing fast-evolving SNe Ia.

Barak Rom, Itai Linial, Re'em Sari

6 pages, Accepted for publication in ApJ

We examine the effects of two-body interactions in a nuclear star cluster surrounding a supermassive black hole. We evaluate the energy flux, analogously to the particle flux calculation of Bahcall and Wolf (1976). We show that there are two types of power-law steady-state solutions: one with zero energy flux and constant particle flux and the other with constant energy flux and zero particle flux. We therefore prove that a zero particle flux solution, which corresponds to the case of an accreting supermassive black hole, can be obtained by requiring a constant energy flux. Consequently, this solution can be derived by simple dimensional analysis, bypassing the need for detailed calculation. Finally, we show that this characteristic, of zero particle flux for constant energy flux and vice versa, is not unique to the Keplerian potential of a supermassive black hole but holds for any central potential of the form $\phi\propto r^{-\beta}$.

Tawny Sit, Mansi M. Kasliwal, Anastasios Tzanidakis, Kishalay De, Christoffer Fremling, Jesper Sollerman, Avishay Gal-Yam, Adam A. Miller, Scott Adams, Robert Aloisi, Igor Andreoni, Matthew Chu, David Cook, Kaustav Kashyap Das, Alison Dugas, Steven L. Groom, Anna Y. Q. Ho, Viraj Karambelkar, James D. Neill, Frank J. Masci, Michael S. Medford, Josiah Purdum, Yashvi Sharma, Roger Smith, Robert Stein, Lin Yan, Yuhan Yao, Chaoran Zhang

32 pages, 17 figures, 5 tables. Submitted to ApJ

SN 1987A was an unusual hydrogen-rich core-collapse supernova originating from a blue supergiant star. Similar blue supergiant explosions remain a small family of events, and are broadly characterized by their long rises to peak. The Zwicky Transient Facility (ZTF) Census of the Local Universe (CLU) experiment aims to construct a spectroscopically complete sample of transients occurring in galaxies from the CLU galaxy catalog. We identify 13 long-rising (>40 days) Type II supernovae from the volume-limited CLU experiment during a 3.5 year period from June 2018 to December 2021, approximately doubling the previously known number of these events. We present photometric and spectroscopic data of these 13 events, finding peak r-band absolute magnitudes ranging from -15.6 to -17.5 mag and the tentative detection of Ba II lines in 9 events. Using our CLU sample of events, we derive a long-rising Type II supernova rate of $1.37^{+0.26}_{-0.30}\times10^{-6}$ Mpc$^{-3}$ yr$^{-1}$, $\approx$1.4% of the total core-collapse supernova rate. This is the first volumetric rate of these events estimated from a large, systematic, volume-limited experiment.

D. J. Matthews, R. Margutti, B. D. Metzger, D. Milisavljevic, G. Migliori, T. Laskar, D. Brethauer, E. Berger, R. Chornock, M. Drout, E. Ramirez-Ruiz

5 pages, 1 figure, 1 table. Submitted to RNAAS

We present the X-ray monitoring campaign of AT2022tsd in the time range $\delta t_{rest} = 23 - 116$ d rest-frame since discovery. With an initial 0.3-10 keV X-ray luminosity of $L_x \approx 10^{44}$ erg s$^{-1}$ at $\delta t_{rest}\approx$ 23 d, AT2022tsd is the most luminous FBOT to date and rivals even the most luminous GRBs. We find no statistical evidence for spectral evolution. The average X-ray spectrum is well described by an absorbed simple power-law spectral model with best-fitting photon index $\Gamma = 1.89 ^{+0.09}_{-0.08}$ and marginal evidence at the 3$\sigma$ confidence level for intrinsic absorption $NH_{int}\approx 4\times10^{19}$ cm$^{-2}$. The X-ray light-curve behavior can be either interpreted as a power-law decay $L_x\propto t^{\alpha}$ with $\alpha\approx -2$ and superimposed X-ray variability, or as a broken power-law with a steeper post-break decay as observed in other FBOTs such as AT2018cow. We briefly compare these results to accretion models of TDEs and GRB afterglow models.

Genaro Suárez (1, 2), Stanimir Metchev (2) ((1) American Museum of Natural History, (2) The University of Western Ontario)

Accepted for publication in MNRAS

Analysis of all archival 5--14 micron spectra of field ultracool dwarfs from the Infrared Spectrograph on the Spitzer Space Telescope has shown that absorption by silicates in the 8--11 micron region is seen in most L-type (1300 K to 2200 K) dwarfs. The absorption is caused by silicate-rich clouds in the atmospheres of L dwarfs and is strongest at L4--L6 spectral types. Herein we compare averages of the mid-infrared silicate absorption signatures of L3--L7 dwarfs that have low ($\lesssim$10$^{4.5}$ cm s$^{-2}$) vs.\ high ($\gtrsim$10$^5$ cm s$^{-2}$) surface gravity. We find that the silicate absorption feature is sensitive to surface gravity and indicates a difference in grain size and composition between dust condensates in young and old mid-L dwarfs. The mean silicate absorption profile of low-gravity mid-L dwarfs matches expectations for $\sim$1 micron-sized amorphous iron- and magnesium-bearing pyroxene (Mg$_x$Fe$_{1-x}$SiO$_3$) grains. High-gravity mid-L dwarfs have silicate absorption better represented by smaller ($\lesssim$0.1 $\mu$m) and more volatile amorphous enstatite (MgSiO$_3$) or SiO grains. This is the first direct spectroscopic evidence for gravity-dependent sedimentation of dust condensates in ultracool atmospheres. It confirms theoretical expectations for lower sedimentation efficiencies in low-gravity atmospheres and independently confirms their increased dustiness.

M. Bzowski, M.A. Kubiak, M. Strumik, I. Kowalska-Leszczynska, C. Porowski, E. Quemerais

Accepted for ApJ

Interstellar neutral hydrogen flows into the heliosphere as a mixture of the primary and secondary populations from two somewhat different directions due to splitting occurring in the magnetized outer heliosheath. The direction of inflow of interstellar neutral H observed in the inner heliosphere, confronted with that of the unperturbed flow of interstellar neutral helium, is important for understanding the geometry of the distortion of the heliosphere from axial symmetry. It is also needed for facilitating remote-sensing studies of the solar wind structure based on observations of the helioglow, such as those presently performed by SOHO/SWAN, and in a near future by IMAP/GLOWS. In the past, the only means to measure the flow direction of interstellar hydrogen were spectroscopic observations of the helioglow. Here, we propose a new method to determine this parameter based on a long series of photometric observations of the helioglow. The method is based on purely geometric considerations and does not depend on any model and absolute calibration of the measurements. We apply this method to sky maps of the helioglow available from the SOHO/SWAN experiment and derive the mean flow longitude of interstellar hydrogen. We obtain $253.1\degr \pm 2.8\degr$, which is in perfect agreement with the previously obtained results based on spectroscopic observations.

Rogelio Medina Covarrubias, Fabio De Colle, Gerardo Urrutia, Felipe Vargas

11 pages, 9 figures, accepted for publication in MNRAS

We compute the linear polarisation during the afterglow phase of gamma-ray bursts, for both on-axis and off-axis observers. We use numerical simulations of the deceleration of a relativistic jet, and compute the polarisation by post-processing the results of the numerical simulations. In our simulations, we consider a magnetic field that is chaotic in the plane of the shock, in addition to a magnetic field component that is parallel to the shock velocity. While the linear polarisation computed for on-axis observers is consistent with previous analytical estimates, we found that lateral expansion, which is accurately handled in our simulations, plays a crucial role in determining the linear polarisation for off-axis observers. Our results show that the off-axis linear polarisation, as seen by off-axis observers, exhibits a single peak, in contrast to the two peaks inferred by previous analytical studies. The maximum polarisation degree is 40\% at an observing angle $\theta_{\rm obs}=0.4$ rad, and it decreases as the observing angle increases, which is opposite to what predicted by analytical models, where polarisation increases with larger observing angles. From the upper limit of 12\% in the linear polarisation obtained at 244 days for the GRB 170817A, we also infer an anisotropy factor of $B_\parallel/B_\perp = 0.5-0.9$, consistent with the post-shock magnetic field being amplified by turbulence.

James Reeves, Valentina Braito, Delphine Porquet, Marco Laurenti, Andrew Lobban, Gabriele Matzeu

15 pages, 10 figures, accepted for publication in the Astrophysical Journal

PG 1448+273 is a luminous, nearby ($z=0.0645$), narrow line Seyfert 1 galaxy, which likely accretes close to the Eddington limit. XMM-Newton observations of PG 1448+273 in 2017 revealed the presence of an ultra fast outflow, as seen through its blueshifted iron K absorption profile, with an outflow velocity of about $0.1c$. Here, the first NuSTAR observation of PG 1448+273, performed in 2022 and coordinated with XMM-Newton is presented, which shows remarkable variability of its ultra fast outflow. The average count rate is a factor of 2 lower during the last 60 ks of the NuSTAR observation, where a much faster component of the ultra fast outflow was detected with a terminal velocity of $0.26\pm0.04c$. This is significantly faster than the outflow component which was initially detected in 2017, when overall PG 1448+273 was observed at a lower X-ray flux and which implies an order of magnitude increase in the wind kinetic power between the 2017 and 2022 epochs. Furthermore, the rapid variability of the ultra fast outflow in 2022, on timescales down to 10 ks, suggests we are viewing through a highly inhomogeneous disk wind in PG 1448+273, where the passage of a denser wind clump could account for the increase in obscuration in the last 60 ks of the NuSTAR observation.

We present mean horizontal branch absolute magnitudes and iron abundances for a sample of 39 globular clusters. These quantities were calculated in an unprecedented homogeneous fashion based on Fourier decomposition of ligt curves of RR Lyrae cluster members. Zero points for the luminosity calibrations are discussed. Our photometrically derived metallicities and distances compare very well with spectroscopic determinations of [Fe/H] and accurate distances obtained using {\sl Gaia} and {\sl Hubble Space Telescope} data. The need to distinguish between the results for RRab and RRc stars for a correct evaluation of the $M_V$--[Fe/H] relation is discussed. For RRab stars, the relation is non-linear, and the horizontal branch structure plays a significant role. For RRc stars, the relation remains linear and tight, and the slope is very shallow. Hence, the RRc stars seem better indicators of the parental cluster distances. Systematic time-series CCD imaging performed over the last 20 years enabled to discover and classify 330 variables in our sample of globular clusters.

Satoshi Itoh, Taro Matsuo, Shunsuke Ota, Kensuke Hara, Yuji Ikeda, Reiki Kojima, Toru Yamada, Takahiro Sumi

18 pages, 10 figures, accepted for the Publications of the Astronomical Society of the Pacific

We performed an experimental verification of a coronagraph. As a result, we confirmed that, at the focal region where the planetary point spread function exists, the coronagraph system mitigates the raw contrast of a star-planet system by at least $1\times10^{-5}$ even for the 1-$\lambda/D$ star-planet separation. In addition, the verified coronagraph keeps the shapes of the off-axis point spread functions when the setup has the source angular separation of 1$\lambda/D$. The low-order wavefront error and the non-zero extinction ratio of the linear polarizer may affect the currently confirmed contrast. The sharpness of the off-axis point spread function generated by the sub-$\lambda/D$ separated sources is promising for the fiber-based observation of exoplanets. The coupling efficiency with a single mode fiber exceeds 50% when the angular separation is greater than 3--4$\times 10^{-1}\lambda/D$. For sub-$\lambda/D$ separated sources, the peak positions (obtained with Gaussian fitting) of the output point spread functions are different from the angular positions of sources; the peak position moved from about $0.8\lambda/D$ to $1.0\lambda/D$ as the angular separation of the light source varies from $0.1\lambda/D$ to $1.0\lambda/D$. The off-axis throughput including the fiber-coupling efficiency (with respect to no focal plane mask) is about 40% for 1-$\lambda/D$ separated sources and 10% for 0.5-$\lambda/D$ separated ones (excluding the factor of the ratio of pupil aperture width and Lyot stop width), where we assumed a linear-polarized-light injection. In addition, because this coronagraph can remove point sources on a line in the sky, it has another promising application for high-contrast imaging of exoplanets in binary systems.

Recently different cosmological measurements have shown a tension in the value of the Hubble constant, $H_0$. Assuming the $\Lambda$CDM model, the Planck satellite mission has inferred the Hubble constant from the cosmic microwave background (CMB) anisotropies to be $H_0 = 67.4 \pm 0.5 \, \rm{km \, s^{-1} \, Mpc^{-1}}$. On the other hand, low redshift measurements such as those using Cepheid variables and supernovae Type Ia (SNIa) have obtained a significantly larger value. For instance, Riess et al. reported $H_0 = 73.04 \pm 1.04 \, \rm{km \, s^{-1} \, Mpc^{-1}}$, which is $5\sigma$ apart of the prediction from Planck observations. This tension is a major problem in cosmology nowadays, and it is not clear yet if it comes from systematic effects or new physics. The use of new methods to infer the Hubble constant is therefore essential to shed light on this matter. In this paper, we discuss using the age of the oldest astrophysical objects (OAO) to probe the Hubble tension. We show that, although this data can provide additional information, the method can also artificially introduce a tension.

W. C. Jing, J. L. Han, Tao Hong, Chen Wang, X. Y. Gao, L. G. Hou, D. J. Zhou, J. Xu, Z. L. Yang

7 pages, 7 figures. Under review in MNRAS

The pulsar B0458+46 was previously believed to have a distance of about 1.3$~$kpc and to be associated with a nearby supernova remnant, SNR HB9 (G160.9+2.6). We observe the neutral hydrogen (HI) absorption spectrum of PSR B0458+46 by using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), and detect two absorption lines at radial velocities of $V_{\rm LSR} = {-7.7}~{\rm km~s}^{-1}$ and $-28.1~{\rm km~s^{-1}}$. Based on the Galactic rotation curve with a modification factor correcting for the systematic stream in the anticenter region, we derive the kinematic distance of the farther absorption clouds, which are found to be located $2.7^{+0.9}_{-0.8}$ kpc away, just beyond the Perseus Arm. We also obtain a direct distance estimation of the absorption clouds, being $2.3_{-0.7}^{+1.1}$ kpc, based on a comparison of their velocity with the HI emission in the Perseus and Outer Arms that was well-defined by recently measured parallax tracers. As a result, we conclude that PSR B0458+46 should be located beyond the Perseus Arm, with a lower limit distance of 2.7 kpc, and therefore not associated with SNR HB9. The doubled distance indicates a deficiency of thermal electrons in the immediate outer Galaxy, with much less density than current models predict. Additionally, we detect a new high-velocity HI cloud in the direction of this pulsar.

Kai Wu, M.B.N. Kouwenhoven, Rainer Spurzem, Xiaoying Pang

17 pages, 20 figures, and 3 tables. Accepted for publication in MNRAS

Although debris disks may be common in exoplanet systems, only a few systems are known in which debris disks and planets coexist. Planets and the surrounding stellar population can have a significant impact on debris disk evolution. Here we study the dynamical evolution of debris structures around stars embedded in star clusters, aiming to determine how the presence of a planet affects the evolution of such structures. We combine NBODY6++GPU and REBOUND to carry out N-body simulations of planetary systems in star clusters (N=8000; Rh=0.78 pc) for a period of 100 Myr, in which 100 solar-type stars are assigned 200 test particles. Simulations are carried out with and without a Jupiter-mass planet at 50 au. We find that the planet destabilizes test particles and speeds up their evolution. The planet expels most particles in nearby and resonant orbits. Remaining test particles tend to retain small inclinations when the planet is present, and fewer test particles obtain retrograde orbits. Most escaping test particles with speeds smaller than the star cluster's escape speed originate from cold regions of the planetary system or from regions near the planet. We identify three regions within planetary systems in star clusters: (i) the private region of the planet, where few debris particles remain (40 - 60 au), (ii) the reach of the planet, in which particles are affected by the planet (0 - 400 au), and (iii) the territory of the planetary system, most particles outside which will eventually escape (0 - 700 au).

CN Cha is a slow symbiotic nova characterized by a three-years-long optical flat peak followed by a rapid decline. We present theoretical light curves for CN Cha, based on hydrostatic approximation, and estimate the white dwarf (WD) mass to be $\sim 0.6 ~M_\odot$ for a low metal abundance of Z = 0.004. This kind of flat peak novae are border objects between classical novae having a sharp optical peak and extremely slow novae, the evolutions of which are too slow to be recognized as a nova outburst in human timescale. Theoretically, there are two types of nova envelope solutions, static and optically-thick wind, in low mass WDs ($\lesssim 0.7 ~M_\odot$). Such a nova outburst begins first in a hydrostatic manner, and later it could change to an optically-thick wind evolution due to perturbation by the companion star in the nova envelope. Multiple peaks are a reflection of the relaxation process of transition. CN Cha supports our explanation on the difference between long-lasted flat peak novae like CN Cha and multiple peak novae like V723 Cas, because the companion star is located far outside, and does not perturb, the nova envelope in CN Cha.

Kunio Kaneta, Hye-Sung Lee, Jiheon Lee, Jaeok Yi

23 pages, 6 figures

A coherent field over the entire universe is an attractive picture in studying the dark sector of the universe. The misalignment mechanism, which relies on inflation to achieve homogeneousness of the field, is a popular mechanism for producing such a coherent dark matter. Nevertheless, unlike a scalar field case, a vector boson field suffers because its energy density is exponentially suppressed by the scale factor during the cosmic expansion. We show that if the vector field gets a mass from a scalar field, whose value increases by orders of magnitude, the suppression can be compensated, and the misalignment can produce the coherent vector boson that has a sizable amount of energy density in the present universe. Quintessence can be such a scalar field.

Li Xiao, Ming Zhu, Xiao-Hui Sun, Peng Jiang, Chun Sun

9 pages, 8 figures, accepted by ApJ

We have obtained the polarization data cube of the VRO 42.05.01 supernova remnant at 1240 MHz using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Three-dimensional Faraday Synthesis is applied to the FAST data to derive the Faraday depth spectrum. The peak Faraday depth map shows a large area of enhanced foreground RM of ~60 rad m-2 extending along the remnant's "wing" section, which coincides with a large-scale HI shell at -20 km/s. The two depolarization patches within the "wing" region with RM of 97 rad m-2 and 55 rad m-2 coincide with two HI structures in the HI shell. Faraday screen model fitting on the Canadian Galactic Plane Survey (CGPS) 1420 MHz full-scale polarization data reveals a distance of 0.7-0.8d_{SNR} in front of the SNR with enhanced regular magnetic field there. The highly piled-up magnetic field indicates that the HI shell at -20 km/s could originate from an old evolved SNR.

Aims. Assuming fast radio bursts (FRBs) are produced by matter travelling ultra-relativistically in a localised region of a smooth bundle of streamlines, we study the constraints applied by geometry to the morphology and polarisation of the burst in time and frequency independently of the intrinsic radiative process. Methods. We express the problem only in terms of the local properties of direction and curvature of a streamline. This allows us to cast the general results to any desired geometry. We illustrate by applying this framework to two geometries inspired by pulsar and magnetar magnetospheres, namely the dipolar polar-cap region and a magnetic dipole with an additional toroidal component. Results. Geometry constrains bursts to occur within an envelope in the frequency vs. time plane (dynamic spectrum). This envelope notably characterises spectral occupancy and frequency drifts (both burst-to-burst and within an individual burst). We illustrate how one can simulate bursts by specifying some basic properties of an intrinsic emission process. In particular we show that the typical properties of one-off bursts can be produced in polar-cap geometry by a star with spin period > 1s, while bursts from repeating sources are better accounted for with an additional strong toroidal component and a sub-second spin period. Conclusions. We propose that a relationship between burst morphologies and the properties of the source, such as its spin period and magnetospheric properties, can be established at least qualitatively based on geometrical considerations. Our results favour models where repeaters are younger and faster magnetars with highly twisted magnetospheres.

Yoshiaki Hagiwara, Willem A. Baan, Masatoshi Imanishi, Philip Diamond

11 pages, 9 figures, submitted to MNRAS on 30 March 2023

The results of high-resolution spectral-line observations of dense molecular gas are presented towards the nuclear region of the type 2 Seyfert galaxy NGC1068. MERLIN observations of the 22 GHz H$_2$O maser were made for imaging the known off-nuclear maser emission at radio jet component located about 0.3" north-east of the radio nucleus in the galaxy. High angular resolution ALMA observations have spatially resolved the molecular gas emissions of HCN and HCO$^{+}$ in this region. The off-nuclear maser spots are found to nearly overlap with a ring-like molecular gas structure and are tracing an evolving shock-like structure, which appears to be energized by interaction between the radio jet and circumnuclear medium. A dynamic jet-ISM interaction is further supported by a systematic shift of the centroid velocities of the off-nuclear maser features over a period of 35 years. The integrated flux ratios of the HCO$^{+}$ line emission features at component C suggest a kinetic temperature T$_{k}$ $\gtrsim$ 300K and an H$_2$ density of $\gtrsim$ 10$^6$ cm$^{-3}$, which are conditions where water masers may be formed. The diagnostics of the masering action in this jet-ISM interaction region is exemplary for galaxies hosting off-nuclear H$_2$O maser emission.

We present a new end-to-end pipeline for Mock Observations of X-ray Halos and Analysis (MOXHA) for hydrodynamic simulations of massive halos, and use it to investigate X-ray scaling relations and hydrostatic mass bias in the Simba cosmological hydrodynamic simulation for halos with $M_{500}\sim 10^{13-15}M_\odot$. MOXHA ties together existing yT-based software packages and adds new functionality to provide an end-to-end pipeline for generating mock X-ray halo data from large-scale or zoom simulation boxes. We compare MOXHA-derived halo properties in Simba to their emission-weighted counterparts, and forecast the systematic mass bias in mock Athena observations. Overall, we find inferred hydrostatic masses are biased low compared to true Simba values. For simple mass-weighting, we find $b_\text{MW} = 0.15^{+0.15}_{-0.14}$ ($16-84\%$ range), while emission-weighting increases this to $b_\text{LW}=0.30^{+0.19}_{-0.10}$. The larger bias versus mass-weighted values we attribute to the spectroscopic and emission-weighted temperatures being biased systematically lower than mass-weighted temperatures. The full MOXHA pipeline recovers the emission-weighted hydrostatic masses at $R_{500}$ reasonably well, yielding $b_\text{X}=0.33^{+0.28}_{-0.34}$. MOXHA-derived halo X-ray scalings are in very good agreement with observed scaling relations, with the inclusion of lower-mass groups significantly steepening the $L_\text{X}-M_{500}$, $M_{500}-T_\text{X}$, and $L_\text{X}-T_\text{X}$ relations. This indicates the strong effect the Simba feedback model has on low-mass halos, which strongly evacuates poor groups but still retains enough gas to reproduce observations. We find similar trends for analogous scaling relations measured at $R_{500}$, as expected for halo-wide gas evacuation.

C. Quintero Noda, E. Khomenko, M. Collados, B. Ruiz Cobo, R. Gafeira, N. Vitas, M. Rempel, R. J. Campbell, A. Pastor Yabar, H. Uitenbroek, D. Orozco Suárez

14 pages, 13 figures

In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magneto-hydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4m class telescopes DKIST and EST. We first study multiple inversion configurations using the SIR code and the Fe I transitions at 630 nm until we obtain minor differences between the input and the inferred atmosphere in a wide range of heights. Also, we examine how the inversion accuracy depends on the noise level of the Stokes profiles. The results indicate that when the majority of the inverted pixels come from strongly magnetised areas, there are almost no restrictions in terms of the noise, obtaining good results for noise amplitudes up to 1$\times10^{-3}$ of $I_c$. At the same time, the situation is different for observations where the dominant magnetic structures are weak, and noise restraints are more demanding. Moreover, we find that the accuracy of the fits is almost the same as that obtained without noise when the noise levels are on the order of 1$\times10^{-4}$of $I_c$. We, therefore, advise aiming for noise values on the order of or lower than 5$\times10^{-4}$ of $I_c$ if observers seek reliable interpretations of the results for the magnetic field vector reliably. We expect those noise levels to be achievable by next-generation 4m class telescopes thanks to an optimised polarisation calibration and the large collecting area of the primary mirror.

Baoqiang Lao, Sumit Jaiswal, Zhen Zhao, Leping Lin, Junyi Wang, Xiaohui Sun, Shengli Qin

Modern large radio continuum surveys have high sensitivity and resolution, and can resolve previously undetected extended and diffuse emissions, which brings great challenges for the detection and morphological classification of extended sources. We present HeTu-v2, a deep learning-based source detector that uses the combined networks of Mask Region-based Convolutional Neural Networks (Mask R-CNN) and a Transformer block to achieve high-quality radio sources segmentation and classification. The sources are classified into 5 categories: Compact or point-like sources (CS), Fanaroff-Riley Type I (FRI), Fanaroff-Riley Type II (FRII), Head-Tail (HT), and Core-Jet (CJ) sources. HeTu-v2 has been trained and validated with the data from the Faint Images of the Radio Sky at Twenty-one centimeters (FIRST). We found that HeTu-v2 has a high accuracy with a mean average precision ($AP_{\rm @50:5:95}$) of 77.8%, which is 15.6 points and 11.3 points higher than that of HeTu-v1 and the original Mask R-CNN respectively. We produced a FIRST morphological catalog (FIRST-HeTu) using HeTu-v2, which contains 835,435 sources and achieves 98.6% of completeness and up to 98.5% of accuracy compared to the latest 2014 data release of the FIRST survey. HeTu-v2 could also be employed for other astronomical tasks like building sky models, associating radio components, and classifying radio galaxies.

Varun, Gayathri Raman

11 pages, 6 figures, 7 tables. Accepted for publication in MNRAS

We present the results obtained from a comprehensive timing and spectral study of two high-mass X-ray binary sources using NuSTAR observations. These two sources, IGR J16320-4751 and IGR J16479-4514, were discovered by INTEGRAL and have been characterized for the first time in the hard X-ray band (beyond 10~keV) in this work. In these sources, we observe the occurrence of intense X-ray flares, with average luminosities exceeding 10$^{36}$~erg~s$^{-1}$. Our analysis reveals that these flares can be described consistently in the quasi-spherical accretion regime. The orbital phase of the first flare in NuSTAR observation of IGR J16479-4514 matches with the orbital phases of previous flares ($\phi=0.35$) in this source detected by other telescopes. We conclude that this flare occurs as a result of the periastron passage of the neutron star, rather than due to the presence of a corotating interaction region (CIR). Furthermore, from the energy-resolved pulse profile analysis of IGR J16320-4751, we find that the pulse fraction is lower in hard X-rays compared to soft X-rays. We present the hard X-ray spectral parameters of these two sources using several standard spectral model components. We do not detect a cyclotron absorption feature in either target. We provide estimates of the surface magnetic field strength of NS in IGR J16320-4751 using two indirect methods. Lastly, we observe spectral hardening during flaring segments compared to the off-flaring segments which indicates that comptonization is more effective during the flaring segments.

E.O. Angüner, G. Spengler, E. Amato, S. Casanova

Accepted for publication in MNRAS. 16 pages, 9 figures

The Pevatron Test Statistic (PTS) is applied to data from $\gamma$-ray observatories to test for the origin of Cosmic Rays (CRs) at energies around the knee of the CR spectrum. Several sources are analyzed within hadronic emission models. Previously derived results for RX J1713.7$-$3946, Vela Jr., and HESS J1745$-$290 are confirmed to demonstrate the concept, reliability, and advantages of the PTS. It is excluded with a significance more than $5\sigma$ that the sources RX J1713.7$-$3946 and Vela Jr. are Pevatrons, while strong indications exceeding $4\sigma$ are found for excluding HESS J1745$-$290 as a Pevatron. The importance to resolve source confusion with high angular resolution observations for Pevatrons searches is demonstrated using PTS for the region containing the SNR G106.3+2.7 and the Boomerang nebula. No statistically significant conclusion with respect to Pevatron associations could be drawn from this region, for the diffuse $\gamma$-ray emission around the Galactic Center, and the unidentified $\gamma$-ray sources LHAASO J2108$+$5157, HESS J1702$-$420A and MGRO J1908$+$06. Assuming the entire $\gamma$-ray emission from MGRO J1908+06 and the tail region of SNR G106.3+2.7 is hadronic, a statistical indication exceeding $3\sigma$ is found for the underlying proton spectrum to extend beyond 350$-$400 TeV as a power-law. This result can indicate that these sources are proton and helium Pevatrons, in which the accelerated particles contribute to the knee of proton and helium spectra observed at Earth.

Ralph P. Eatough, Gregory Desvignes, Kuo Liu, Robert S. Wharton, Aristedis Noutsos, Pablo Torne, Ramesh Karuppusamy, Lijing Shao, Michael Kramer, Heino Falcke, Luciano Rezzolla

4 pages, 3 figures, Proceedings of The Fifteenth Marcel Grossmann Meeting on General Relativity, University of Rome La Sapienza, Italy, 1 - 7 July 2018

BlackHoleCam is a project funded by a European Research Council Synergy Grant to build a complete astrophysical description of nearby supermassive black holes by using a combination of radio imaging, pulsar observations, stellar astrometry and general relativistic magneto-hydrodynamic models. BlackHoleCam scientists are active partners of the Event Horizon Telescope Consortium. In this talk I will discuss the use of pulsars orbiting Sagittarius A* for tests of General Relativity, the current difficulties in detecting such sources, recent results from the Galactic Centre magnetar PSR J1745-2900 and how BlackHoleCam aims to search for undiscovered pulsars in the Galactic Centre.

Ye-Peng Yan, Guo-Jian Wang, Si-Yu Li, Yang-Jie Yan, Jun-Qing Xia

12 pages, 8 figures, accepted by ApJ

The lensing effect of the cosmic microwave background (CMB) is a powerful tool for our study of the distribution of matter in the universe. Currently, the quadratic estimator (EQ) method, which is widely used to reconstruct lensing potential, has been known to be sub-optimal for the low-noise levels polarization data from next-generation CMB experiments. To improve the performance of the reconstruction, other methods, such as the maximum likelihood estimator and machine learning algorithms are developed. In this work, we present a deep convolutional neural network model named the Residual Dense Local Feature U-net (RDLFUnet) for reconstructing the CMB lensing convergence field. By simulating lensed CMB data with different noise levels to train and test network models, we find that for noise levels less than $5\mu$K-arcmin, RDLFUnet can recover the input gravitational potential with a higher signal-to-noise ratio than the previous deep learning and the traditional QE methods at almost the entire observation scales.

Dennis Zaritsky, Richard Donnerstein, Arjun Dey, Ananthan Karunakaran, Jennifer Kadowaki, Donghyeon J. Khim, Kristine Spekkens, Huanian Zhang

21 pages, accepted for publication in ApJS, full catalog available upon request

We present the completed catalog of ultra-diffuse galaxy (UDG) candidates (7070 objects) from our search of the DR9 Legacy Survey images, including distance and total mass estimates for 1529 and 1436 galaxies, respectively, that we provide and describe in detail. From the sample with estimated distances, we obtain a sample of 585 UDGs ($\mu_{0,g} \ge 24$ mag arcsec$^{-2}$ and $r_e \ge 1.5$ kpc) over 20,000 sq. deg of sky in various environments. We conclude that UDGs in our sample are limited to $10^{10} \lesssim$ M$_h$/M$_\odot \lesssim 10^{11.5}$ and are on average a factor of 1.5 to 7 deficient in stars relative to the general population of galaxies of the same total mass. That factor increases with increasing galaxy size and mass up to a factor of $\sim$10 when the total mass of the UDG increases beyond M$_h = 10^{11}$ M$_\odot$. We do not find evidence that this factor has a dependence on the UDG's large-scale environment.

Charlotte O. G. Waterfall, Silvia Dalla, Mike S. Marsh, Timo Laitinen, Adam Hutchinson

Accepted for publication in Space Weather 13 pages, 5 figures

The forecasting of solar energetic particles (SEPs) is a prominent area of space weather research. Numerous forecasting models exist that predict SEP event properties at proton energies <100MeV. One of these models is the SPARX system, a physics-based forecasting tool that calculates >10MeV and >60MeV flux profiles within minutes of a flare being detected. This work describes SPARX-H, the extension of SPARX to forecast SEP events above 300MeV . SPARX-H predicts fluxes in three high energy channels up to several hundred MeV. Correlations between SEP peak flux and peak intensity of the associated solar flare are seen to be weak at high energies, but improved when events are grouped based on the field polarity during the event. Initial results from this new high energy forecasting tool are presented here and the applications of high energy forecasts are discussed. Additionally, the new high energy version of SPARX is tested on a set of historic SEP events. We see that SPARX-H performs best when predicting peak fluxes from events with source locations in well-connected regions, where many large SEP events tend to originate.

Alexander Men'shchikov

20 pages, 17 figures, 5 tables, accepted by Astronomy & Astrophysics

A simple Gaussian size deconvolution method is routinely used to remove the blur of observed images caused by insufficient angular resolutions of existing telescopes, thereby to estimate the physical sizes of extracted sources and filaments. The size deconvolution method is expected to work when the structures, as well as the telescope beams, have Gaussian shapes. This study employs model images of the spherical and cylindrical objects with Gaussian and power-law shapes, representing the dense cores and filaments. The images are convolved to a wide range of angular resolutions to probe various degrees of resolvedness of the models. Simplified flat, convex, and concave backgrounds are added to the images, then planar backgrounds across the footprints of the structures are subtracted and sizes are measured and deconvolved. When background subtraction is inaccurate, the structures acquire profoundly non-Gaussian profiles. The deconvolved half maximum sizes can be strongly under- or overestimated, by factors of up to ~20 when the structures are unresolved or partially resolved. For resolved structures, the errors are within a factor of ~2, although for some power-law models show the factors of up to ~6. The size deconvolution method cannot be applied to unresolved structures, it can be used only for the Gaussian-like structures, including the critical Bonnor-Ebert spheres, when they are at least partially resolved. The method must be considered inapplicable for the power-law structures with shallow profiles. This work also reveals subtle properties of convolution for different geometries. When convolved with different kernels, spherical objects and cylindrical filaments with identical profiles obtain different widths and shapes. A filament, imaged by the telescope with a non-Gaussian PSF, could appear substantially shallower than the structure is in reality, even when it is resolved.

M. Van Cuyck, N. Ponthieu, G. Lagache, A. Beelen, M. Béthermin, A. Gkogkou, M. Aravena, A. Benoit, J. Bounmy, M. Calvo, A. Catalano, F.X. Désert, F.-X. Dupé, A. Fasano, A. Ferrara, J. Goupy, C. Hoarau, W. Hu, J.-C Lambert, J. F. Macías-Pérez, J. Marpaud, G. Mellema, A. Monfardini, A. Pallottini

15 pages, 12 figures, to be published in Astronomy & Astrophysics

CONCERTO is the first experiment to perform a [CII] line intensity mapping survey to target $z>5.2$. Measuring the [CII] power spectrum allows us to study the role of dusty star-forming galaxies in the star formation history during the Reionization and post-Reionization. The main obstacle to this measurement is the contamination by bright foregrounds. We evaluate our ability to retrieve the [CII] signal in mock observations using the Simulated Infrared Dusty Extragalactic Sky. We compared two methods for dealing with the dust continuum emission from galaxies: the standard PCA and the arPLS method. For line interlopers, the strategy relies on masking low-redshift galaxies using external catalogues. As we do not have observations of CO or classical CO proxies ,we relied on the COSMOS stellar mass catalogue. To measure the power spectrum of masked data, we adapted the P of K EstimatoR and discuss its use on LIM data. The arPLS method achieves a reduction of the continuum background to a sub-dominant level of the [CII] at z=7 by a factor of>70. When using PCA, this factor is only 0.7. The masking lowers the power amplitude of line contamination down to $2 \times 10^2 Jy^2/sr$ This residual level is dominated by faint undetected sources. For our [CII] model, this results in a detection at z = 5.2 with a power ratio [CII]/(residual interlopers) = $62 \pm 32$ for a 22 % area survey loss. However, at z = 7, [C II ] / (residual interlopers)$=2.0 \pm 1.4$. Thanks to the large area covered by SIDES-Uchuu, we show that the power amplitude of line residuals varies by 12-15% for z=5.2-7. We present an end-to-end simulation of the extragalactic foreground removal that we ran to detect the [CII] at high redshift via its power spectrum. We show that dust continuum emission are not a limiting foreground for [CII] LIM. Residual CO and [CI] limits our ability to measure the [CII] power spectrum at z>7.

Pooja Pareek, Akhilesh Nautiyal

14 pages, 5 figures

The large value of non-minimal coupling constant $\xi$ required to satisfy CMB observations in Higgs inflation violates unitarity. In this work we study Higgs-inflation with non-canonical kinetic term of DBI form to find whether $\xi$ can be reduced. To study the inflationary dynamics, we transform the action to the Einstein frame, in which the Higgs is minimally coupled to gravity with a non-canonical kinetic term and modified potential. We choose the Higgs self coupling constant $\lambda=0.14$ for our analysis. We find that the value of $\xi$ can be reduced from $10^{3}-10^{4}$ to $\mathcal{O}(10)$ to satisfy Planck constraints on amplitude of scalar power spectrum. However, this model produces a larger tensor-to-scalar ratio $r$, in comparison to the Higgs inflation with canonical kinetic term. We also find that, to satisfy joint constraints on scalar spectral index $n_s$ and tensor-to-scalar ratio $r$ from Planck-2018 and bounds on $r$ from Planck and BICEP3, the value of $\xi$ should be of the order of $10^4$.

Interacting dark energy (IDE) scenario is a natural and important extension to the standard $\Lambda$CDM cosmology. We develop a full numerical routine, called IDECAMB, as a patch to the public Einstein-Boltzmann solver CAMB, to solve the background and perturbation equations of the IDE models. The IDECAMB solver provides a unified interface for the widely studied IDE models by employing a parametrization model with five free functions. By configuring these five functions, one can easily map the coupled quintessence (CQ) and coupled fluid (CF) models into the parametrization. We handle the perturbation evolutions of the CF models with the parametrized post-Friedmann (PPF) approach to avoid the possible large-scale instability. Compared with the previous established PPF approach whose form depends on a specific IDE model, the PPF approach in this work are model-independent, making it easy to use. We constrain a specific CQ model with the IDECAMB package. The fitting results are consistent with those obtained by Planck Collaboration, which confirms the validity of the package.

Richard Monier, E. Niemczura, D.W. Kurtz, S. Rappaport, D.M. Bowman, Simon J. Murphy, Yveline Lebreton, Remko Stuik, Morgan Deal, Thibault Merle, Tolgahan Kılıçoğlu, Marwan Gebran, Ewen Le Ster

54 pages, accepted in The Astronomical Journal. arXiv admin note: substantial text overlap with arXiv:1908.05023

We report on a detailed abundance study of six bright, mostly southern, slowly rotating late B stars: HD~1279 (B8III), HD~99803 (B9V), HD~123445 (B9V), HD~147550 (B9V), HD~171961 (B8III) and HD~202671 (B5II/III), hitherto reported as normal stars. We compare them to the two classical HgMn stars $\mu$ Lep and $\beta$ Scl and to the superficially normal star, $\nu$ Cap. In the spectra of the six stars, the \ion{Hg}{2} line at 3984 \AA\ line is clearly seen and numerous lines of P, Ti, Mn, Fe, Ga, Sr, Y, and Zr appear to be strong absorbers. A comparison of newly acquired and archival spectra of these objects with a grid of synthetic spectra for selected unblended lines reveals large overabundances of P, Ti, Cr, Mn, Sr, Y, Zr, Ba, Pt and Hg and underabundances of He, Mg, Sc and Ni. The effective temperatures, surface gravities, low projected rotational velocities and the peculiar abundance patterns of the six investigated stars show that they are new chemically peculiar stars, mostly new HgMn stars, and are reclassified as such. The evolutionary status of these stars has been inferred and their ages and masses estimated. The two most massive objects, HD~1279 and HD~202671, might have evolved away from the main-sequence recently, the other stars are main-sequence objects. HD~99803A is a sharp lined HgMn star with grazing eclipses; from TESS and MASCARA photometry we determine an orbital period of $P_{\rm orb} = 26.12022 \pm 0.00004$\,d.

Rasmus Maråk, Emmanuel Blazquez, Pablo Gómez

To be presented at the ICATT 2023 conference

Proper modelling of the gravitational fields of irregularly shaped asteroids and comets is an essential yet challenging part of any spacecraft visit and flyby to these bodies. Accurate density representations provide crucial information for proximity missions which rely heavily on it to design safe and efficient trajectories. This work explores using a spacecraft swarm to maximise the measured gravitational signal in a hypothetical mission around the comet 67P/Churyumov-Gerasimenko. Spacecraft trajectories are simultaneously propagated with an evolutionary optimisation approach to maximise overall signal return. The propagation is based on an open-source polyhedral gravity model using a detailed mesh of 67P and considers the comet's sidereal rotation. We compare performance on a mission scenario using one and four spacecraft. The results show that the swarm achieved almost twice the single spacecraft coverage over a fixed mission duration. However, optimising for a single spacecraft results in a more effective trajectory. Overall, this work serves as a testbed for efficiently designing a set of trajectories in this complex gravitational environment balancing measured signals and risks in a swarm scenario. The codebase and results are publicly available at https://github.com/rasmusmarak/TOSS

Ye-Wei Mao, Jun-Yu Gong, Hua Gao, Si-Yue Yu

48 Pages, 38 Figures, 5 Tables; Accepted for Publication in The Astrophysical Journal Supplement Series

A panchromatic investigation of morphology for the early-type spiral galaxy M81 is presented in this paper. We perform bulge-disk decomposition in M81 images at totally 20 wavebands from FUV to NIR obtained with GALEX, Swift, SDSS, WIYN, 2MASS, WISE, and Spitzer. Morphological parameters such as Sersic index, effective radius, position angle, and axis ratio for the bulge and the disk are thus derived at all the wavebands, which enables quantifying the morphological K-correction for M81 and makes it possible to reproduce images for the bulge and the disk in the galaxy at any waveband. The morphology as a function of wavelength appears as a variable-slope trend of the Sersic index and the effective radius, in which the variations are steep at UV--optical and shallow at optical--NIR bands; the position angle and the axis ratio keep invariable at least at optical--NIR bands. It is worth noting that, the Sersic index for the bulge reaches to about 4--5 at optical and NIR bands, but drops to about 1 at UV bands. This difference brings forward a caveat that, a classical bulge is likely misidentified for a pseudo-bulge or no bulge at high redshifts where galaxies are observed through rest-frame UV channels with optical telescopes. The next work of this series is planned to study spatially resolved SEDs for the bulge and the disk, respectively, and thereby explore stellar population properties and star formation/quenching history for the the galaxy composed of the subsystems.

Rebecca Lin, Marten H. van Kerkwijk, Franz Kirsten, Ue-Li Pen, Adam T. Deller

16 pages, 8 figures, submitted to ApJ

We use four observations with the European VLBI network to measure the first precise radio parallax of the Crab Pulsar. We found two in-beam extragalactic sources just outside the Crab Nebula, with one bright enough to use as a background reference source in our data. We use the Crab Pulsar's giant pulses to determine fringe and bandpass calibration solutions, which greatly improved the sensitivity and reliability of our images and allowed us to determine precise positional offsets between the pulsar and the background source. From those offsets, we determine a parallax of $\pi=0.53\pm0.06\rm{\;mas}$ and proper motion of $(\mu_{\alpha},\mu_{\delta})=(-11.34\pm0.06,2.65\pm0.14)\rm{\;mas\;yr^{-1}}$, yielding a distance of $d=1.90^{+0.22}_{-0.18}\rm{\;kpc}$ and transverse velocity of $v_{\perp}=104^{+13}_{-11}\rm{\;km\;s^{-1}}$. These results are consistent with the Gaia 3 measurements, and open up the possibility of far more accurate astrometry with further VLBI observations.

R. Tripodi, J. Scholtz, R. Maiolino, S. Fujimoto, S. Carniani, J. D. Silverman, C. Feruglio, M. Ginolfi, L. Zappacosta, T. Costa, G. C. Jones, E. Piconcelli, M. Bischetti, F. Fiore

10 pages, 7 figures, 3 Tables. Submitted to A&A

We present ALMA deep observations of the [CII] 158 $\mu$m emission line and the continuum at 253 GHz and 99 GHz towards SDSS J0100+2802 at $z\simeq 6.3$, the most luminous QSO at z$>$6. It belongs to the HYPERION sample of luminous QSOs at $z\sim 6-7.5$. The observations (at 2.2$''$ resolution in band 3 and 0.9$''$ resolution in band 6) are optimized to detect extended emission around the QSO. We detect a merging, tidally disrupted companion both in [CII] and in continuum, stretching on scales up to 20 kpc from the quasar, with a knotty morphology. For the newly-detected companion we estimate a dust mass of $M_{\rm dust}=(0.6-4.3)\times 10^7\ \rm M_\odot$, an SFR in the range $[43-402]\ \rm M_\odot$, that is remarkably similar to the SFR of the QSO, and a neutral gas mass of $M_{\rm HI}=3.3\times 10^9\ \rm M_{\odot}$, suggesting that both the QSO and its companion are gas rich and that the major merging may be at the origin of the boosted star formation. This close merging companion is undetected by deep JWST imaging observations, supporting the effectiveness of ALMA in detecting dust obscured sources especially in the vicinity of optically bright quasars. We also detect a broad blueshifted component in the [CII] spectrum aligned with the radio jet of the QSO, suggesting that this may be the first detection of a radio jet - driven outflow at such high redshift. We estimate a mass outflow rate in the range $\dot{M}_{\rm out}=(115-269)\ \rm M_\odot\ yr^{-1}$. The outflow energetics is similar to that of ionized outflows found in other QSOs host at lower redshift, and the low momentum loading factor suggests that this outflow would not be very effective in removing the gas from the entire galaxy. These results highlight the importance of deep medium-resolution ALMA observations for the study of QSOs and their environment at the Epoch of Reionization.

A. Morbidelli, K. Batygin, E. Lega

In press in Astronomy & Astrophysics

Context: Radius and mass measurements of short-period giant planets reveal that many of these planets contain a large amount of heavy elements, in sharp contrast with the expectations of the conventional core-accretion model for the origin of giant planets. Aims: The proposed explanations for the heavy-element enrichment of giant planets fall short of explaining the most enriched planets. We look for additional processes that can explain the full envelope of inferred enrichments. Methods: We revisit the dynamics of pebbles and dust in the vicinity of giant planets using analytic estimates. Although our results are derived in the framework of a viscous alpha-disk we also discuss the case of disks driven by angular momentum removal in magnetized winds. Results: When giant planets are far from the star, dust and pebbles are confined in a pressure bump at the outer edge of the planet-induced gap. Instead, when the planets reach the inner part of the disk (r << 2 au), dust penetrates the gap together with the gas. The dust/gas ratio can be enhanced by more than an order of magnitude if radial drift of dust is not impeded farther out by other barriers. Thus, hot planets undergoing runaway gas accretion can swallow a large amount of dust. Conclusions: Whereas the gas accreted by giant planets in the outer disk is very dust-poor, that accreted by hot planets can be extremely dust-rich. Thus, provided that a large fraction of the atmosphere of hot-Jupiters is accreted in situ, a large amount of dust can be accreted as well. We draw a distinction between this process and pebble accretion, which is ineffective at small stellocentric radii, even for super-Earths. Giant planets farther out in the disk are extremely effective barriers against the flow of pebbles and dust across their gap.

Natalie O. Butterfield, David T. Chuss, Jordan A. Guerra, Mark R. Morris, Dylan Pare, Edward J. Wollack, C. Darren Dowell, Matthew J. Hankins, Javad Siah, Johannes Staguhn, Ellen Zweibel

Submitted to ApJ

We present the first data release (DR1) of the Far-Infrared Polarimetric Large Area CMZ Exploration (FIREPLACE) survey. The survey was taken using the 214-micron band of the HAWC+ instrument with the SOFIA telescope (19.6" resolution; 0.7 pc). In this first data release we present dust polarization observations covering a ~0.5degree region of the Galactic Center's Central Molecular Zone (CMZ), centered on the Sgr B2 complex. We detect ~25,000 Nyquist-sampled polarization pseudovectors, after applying the standard SOFIA cuts for minimum signal-to-noise in fractional polarization and total intensity of 3 and 200, respectively. Analysis of the magnetic field orientation suggests a bimodal distribution in the field direction. This bimodal distribution shows enhancements in the distribution of field directions for orientations parallel and perpendicular to the Galactic plane, which is suggestive of a CMZ magnetic field configuration with polodial and torodial components. Furthermore, a detailed analysis of individual clouds included in our survey (i.e., Sgr B2, Sgr B2-NW, Sgr B2-Halo, Sgr B1, and Clouds-E/F) shows these clouds have fractional polarization values of 1-10% at 214-micron, with most of the emission having values <5%. A few of these clouds (i.e., Sgr B2, Clouds-E/F) show relatively low fractional polarization values toward the cores of the cloud, with higher fractional polarization values toward the less dense periphery. We also observe higher fractional polarization towards compact HII regions which could indicate an enhancement in the grain alignment in the dust surrounding these sources.

Hiroki Kawashimo, Ryo Sawada, Yudai Suwa, Takashi J. Moriya, Ataru Tanikawa, Nozomu Tominaga

16 pages, 12 figures

Nuclear reactions are key to our understanding of stellar evolution, particularly the $^{12}\rm{C}\left({\alpha},{\gamma}\right)^{16}\!\rm{O}$ rate, which is known to significantly influence the lower and upper ends of the black hole (BH) mass distribution due to pair-instability supernovae (PISNe). However, these reaction rates have not been sufficiently determined. We use the $\texttt{MESA}$ stellar evolution code to explore the impact of uncertainty in the $^{12}\rm{C}\left({\alpha},{\gamma}\right)^{16}\!\rm{O}$ rate on PISN explosions, focusing on nucleosynthesis and explosion energy by considering the high resolution of the initial mass. Our findings show that the mass of synthesized radioactive nickel ($^{56}{\rm Ni}$) and the explosion energy increase with $^{12}\rm{C}\left({\alpha},{\gamma}\right)^{16}\!\rm{O}$ rate for the same initial mass, except in the high-mass edge region. With a high (about twice the $\texttt{STARLIB}$ standard value) rate, the maximum amount of nickel produced falls below 70 $M_\odot$, while with a low rate (about half of the standard value) it increases up to 83.7 $M_\odot$. These results highlight that carbon burning plays a crucial role in PISNe by determining when a star initiates expansion. The initiation of expansion competes with collapse caused by helium photodisintegration, and the maximum mass that can lead to an explosion depends on the $^{12}\rm{C}\left({\alpha},{\gamma}\right)^{16}\!\rm{O}$ reaction rate.

Keitaro Ishikawa, Tomomi Sunayama, Atsushi J. Nishizawa, Hironao Miyatake

10 pages, 11 figures, 3 tables

We investigate the robustness of baryon acoustic oscillations (BAO) measurements with a photometric galaxy sample using mock galaxy catalogs with various sizes of photometric redshift (photo-$z$) uncertainties. We first investigate the robustness of BAO measurements, assuming we have a perfect knowledge of photo-$z$ uncertainties. We find that the BAO shift parameter $\alpha$ can be constrained in an unbiased manner for various sizes of photometric redshift uncertainties at $z=0.251$, $0.617$, and $1.03$ as long as the number density of galaxies is high. A sparse galaxy sample causes additional noise in the covariance matrix calculation and it can bias the constraint on $\alpha$. Next, we investigate the scenario where incorrect photometric redshift uncertainties are assumed in the fitting model and find that underestimating the photo-$z$ uncertainty leads to a degradation in the constraining power on $\alpha$. In addition, we investigate BAO measurements with a cross-correlation signal between a spec-$z$ sample and a photo-$z$ sample. We find BAO constraints are unbiased and slightly tighter than the auto-correlation signal of a photo-$z$ sample. We also quantify the constraining power on $\Omega_{\rm m0}$ assuming the LSST-like covariance and find that the 95\% confidence level is $\sigma(\Omega_{\rm m0})\sim0.03$-$0.05$ corresponding to the photo-$z$ uncertainties of 1\% to 3\% respectively. Finally, we examine whether the skewness in the photometric redshift can bias the constraint on $\alpha$ and confirm that the constraint on $\alpha$ is unbiased even if we use a fitting model assuming a Gaussian photo-$z$ uncertainty.

J. García-Vázquez (1), William J. Henney (2), H. O. Castañeda (1) ((1) ESFM-IPN, Mexico, (2) IRyA-UNAM, Morelia, Mexico)

Accepted for publication in MNRAS. 28 pages with numerous figures, tables, and appendices. Deluxe addition with supplementary figures and tables included

Radial velocity fluctuations on the plane of the sky are a powerful tool for studying the turbulent dynamics of emission line regions. We conduct a systematic statistical analysis of the H alpha velocity field for a diverse sample of 9 H II regions, spanning two orders of magnitude in size and luminosity, located in the Milky Way and other Local Group galaxies. By fitting a simple model to the second-order spatial structure function of velocity fluctuations, we extract three fundamental parameters: the velocity dispersion, the correlation length, and the power law slope. We determine credibility limits for these parameters in each region, accounting for observational limitations of noise, atmospheric seeing, and the finite map size. The plane-of-sky velocity dispersion is found to be a better diagnostic of turbulent motions than the line width, especially for lower luminosity regions where the turbulence is subsonic. The correlation length of velocity fluctuations is found to be always roughly 2% of the H II region diameter, implying that turbulence is driven on relatively small scales. No evidence is found for any steepening of the structure function in the transition from subsonic to supersonic turbulence, possibly due to the countervailing effect of projection smoothing. Ionized density fluctuations are too large to be explained by the action of the turbulence in any but the highest luminosity sources. A variety of behaviors are seen on scales larger than the correlation length, with only a minority of sources showing evidence for homogeneity on the largest scales.

David Wittman (1), Rodrigo Stancioli (1), Kyle Finner (2), Faik Bouhrik (3), Reinout van Weeren (4), Andrea Botteon (5) ((1) University of California, Davis, (2) IPAC, California Institute of Technology, (3) California Northstate University, (4) Leiden Observatory, (5) INAF - IRA)

accepted to ApJ

We report the discovery of a binary galaxy cluster merger via a search of the redMaPPer optical cluster catalog, with a projected separation of 535 kpc between the BCGs. Archival XMM-Newton spectro-imaging reveals a gas peak between the BCGs, suggesting a recent pericenter passage. We conduct a galaxy redshift survey to quantify the line-of-sight velocity difference ($153\pm281$ km/s) between the two subclusters. We present weak lensing mass maps from archival HST/ACS imaging, revealing masses of $M_{200}=4.5\pm0.8\times10^{14}$ and $2.8\pm0.7\times10^{14}$ M$_\odot$ associated with the southern and northern galaxy subclusters respectively. We also present deep GMRT 650 MHz data revealing extended emission, 420 kpc long, which may be an AGN tail but is potentially also a candidate radio relic. We draw from cosmological n-body simulations to find analog systems, which imply that this system is observed fairly soon (60-271 Myr) after pericenter, and that the subcluster separation vector is within 22$^\circ$ of the plane of the sky, making it suitable for an estimate of the dark matter scattering cross section. We find $\sigma_{\rm DM}=1.1\pm0.6$ cm$^2$/g, suggesting that further study of this system could support interestingly tight constraints.

A. Aguilar-Arevalo, I. Arnquist, N. Avalos, L. Barak, D. Baxter, X. Bertou, I.M. Bloch, A.M. Botti, M. Cababie, G. Cancelo, N. Castelló-Mor, B.A. Cervantes-Vergara, A.E. Chavarria, J. Cortabitarte-Gutiérrez, M. Crisler, J. Cuevas-Zepeda, A. Dastgheibi-Fard, C. De Dominicis, O. Deligny, A. Drlica-Wagner, J. Duarte-Campderros, J.C. D'Olivo, R. Essig, E. Estrada, J. Estrada, E. Etzion, F. Favela-Perez, N. Gadola, R. Gaïor, S.E. Holland, T. Hossbach, L. Iddir, B. Kilminster, Y. Korn, A. Lantero-Barreda, I. Lawson, S. Lee, A. Letessier-Selvon, P. Loaiza, A. Lopez-Virto, S. Luoma, E. Marrufo-Villalpando, K.J. McGuire, G.F. Moroni, S. Munagavalasa, D. Norcini, A. Orly, G. Papadopoulos, S. Paul, S.E. Perez, A. Piers, P. Privitera, P. Robmann, D. Rodrigues, N.A. Saffold, S. Scorza, et al. (16 additional authors not shown)

8 pages, 4 figures

We present results from a 3.1 kg-day target exposure of two charge-coupled devices (CCDs), each with 24 megapixels and skipper readout, deployed in the DAMIC (DArk Matter In CCDs) setup at SNOLAB. With a reduction in pixel readout noise of a factor of 10 relative to the previous detector, we investigate the excess population of low-energy bulk events previously observed above expected backgrounds. We address the dominant systematic uncertainty of the previous analysis through a depth fiducialization designed to reject surface backgrounds on the CCDs. The measured bulk ionization spectrum confirms with higher significance the presence of an excess population of low-energy events in the CCD target with characteristic rate of ${\sim}7$ events per kg-day and electron-equivalent energies of ${\sim}80~$eV, whose origin remains unknown.

We comment on the methods and the conclusion of Ref. [1], "Does gravitational confinement sustain flat galactic rotation curves without dark matter?" The article employs two methods to investigate whether non-perturbative corrections from General Relativity are important for galactic rotation curves, and concludes that they are not. This contradicts a series of articles [2-4] that had determined that such corrections are large. We comment here that Ref. [1] use approximations known to exclude the specific mechanism studied in [2-4] and therefore is not testing the finding of Refs. [2-4].

Alessandro Lella, Pierluca Carenza, Giampaolo Co', Giuseppe Lucente, Maurizio Giannotti, Alessandro Mirizzi, Thomas Rauscher

12 pages, 6 figures

Axion-like particles (ALPs) coupled to nucleons might be copiously emitted from a supernova (SN) core. We extend existing bounds on free-streaming ALPs to the case in which these are so strongly-interacting with the nuclear matter to be trapped in the SN core. For strongly-interacting ALPs, we also extend the bound from the absence of an ALP-induced signal in Kamiokande-II neutrino detector at the time of SN 1987A. We find that combining the different arguments, SNe exclude values of axion-nucleon coupling $g_{aN}\gtrsim10^{-9}$ for ALP masses $m_{a}\lesssim1$ MeV. Remarkably, in the case of canonical QCD axion models, the SN bounds exclude all values of $m_{a}\gtrsim 10^{-2}$ eV. This result prevents the possibility for current and future cosmological surveys to detect any axion signal.

Paul A. Cassak, M. Hasan Barbhuiya, Haoming Liang, Matthew R. Argall

16 pages, 3 figures, includes both main paper and supplementary material

Weakly collisional and collisionless plasmas are typically far from local thermodynamic equilibrium (LTE), and understanding energy conversion in such systems is a forefront research problem. The standard approach is to investigate changes in internal (thermal) energy and density, but this omits energy conversion that changes any higher order moments of the phase space density. In this study, we calculate from first principles the energy conversion associated with all higher moments of the phase space density for systems not in LTE. Particle-in-cell simulations of collisionless magnetic reconnection reveal that energy conversion associated with higher order moments can be locally significant. The results may be useful in numerous plasma settings, such as reconnection, turbulence, shocks, and wave-particle interactions in heliospheric, planetary, and astrophysical plasmas.

Metric perturbations in General Relativity are usually separated into three distinct classes: scalar, vector, and tensor. In many cases these modes are separable, i.e. they satisfy independent equations of motion for each mode. However, in the present paper we argue that in many cases tensor and scalar modes are not separable, no matter what gauge conditions are chosen. The propagation of any of these mode depends on the other. A realistic example providing such mixing is presented.

Kohei Fujikura, Yuichiro Nakai, Ryosuke Sato, Yaoduo Wang

17 pages, 2 figures

We investigate cosmological phase transitions in various composite Higgs models consisting of four-dimensional asymptotically-free gauge field theories. Each model may lead to a confinement-deconfinement transition and a phase transition associated with the spontaneous breaking of a global symmetry that realizes the Standard Model Higgs field as a pseudo-Nambu-Goldstone boson. Based on the argument of universality, we discuss the order of the phase transition associated with the global symmetry breaking by studying the renormalization group flow of the corresponding linear sigma model at finite temperature, which is calculated by utilizing the $\epsilon$-expansion technique at the one-loop order. Our analysis indicates that some composite Higgs models accommodate phenomenologically interesting first-order phase transitions. We also explore the confinement-deconfinement transition in a UV-completed composite Higgs model based on a $Sp(2N_c)$ gauge theory. It is found that the first-order phase transition is favored when the number of degrees of freedom for the $Sp(2N_c)$ gauge field is much larger than that of matter fields in the fundamental representation of $Sp(2N_c)$. We comment on the gravitational wave signal generated by the confinement-deconfinement transition and its detectability at future observations. Our discussions motivate further studies on phase transitions in composite Higgs models with the use of lattice simulations.

This paper utilized the high temporal and spatial resolution temperature profile data observed by the multi-channel microwave radiometer at the Large High Altitude Air Shower Observatory (LHAASO) on the eastern slope of the Tibetan Plateau from February to May and August to November 2021, combined with the ERA5 reanalysis data products for the whole year of 2021, to study the daily, monthly, and seasonal variations of the atmospheric boundary layer height (ABLH). The results are as follows: (1) The ABLH on sunny days showed obvious fluctuations with peaks and valleys. The ABLH continued to rise with the increase of surface temperature after sunrise and usually reached its maximum value in the afternoon around 18:00, then rapidly decreased until sunset. (2) The average ABLH in April was the highest at about 1200 m, while it was only around 600 m in November. The ABLH fluctuated greatly during the day and was stable at around 400 m at night. The ABLH results obtained from ERA5 were slightly smaller overall but had a consistent trend of change with the microwave radiometer. (3) The maximum ABLH appeared in spring, followed by summer and autumn, and winter had the lowest value, with all peaks reached around 14:00-15:00. These results are of great significance for understanding the ABLH on the eastern slope of the Tibetan Plateau, and provide reference for the absolute calibration of photon numbers of the LHAASO telescope and the atmospheric monitoring plan, as well as for evaluating the authenticity and accuracy of existing reanalysis datasets.

We study evolution of quantum fluctuations of gravity around an inflationary solution in renormalizable quantum gravity that asymptotically indicates background freedom represented by a special conformal invariance. Inflation ignites at the Planck scale and continues until spacetime phase transition occurs at a dynamical scale about $10^{17}$GeV. We can show that during inflation, the initially large scale-invariant fluctuations reduce in amplitude to the appropriate magnitude suggested by tiny CMB anisotropies. The goal of this research is to derive spectra of fluctuations at the phase transition point, that is, the primordial spectra. A system of nonlinear evolution equations for the fluctuations is derived from the quantum gravity effective action. The running coupling constant is then expressed by a time-dependent average following the spirit of the mean field approximation. In this paper, we determine and examine various nonlinear terms, not treated in previous studies such as the exponential factor of the conformal mode, that contribute at the early stage of inflation with still large amplitude. Moreover, in order to verify their effects concretely, we numerically solve the evolution equation by making a simplification to extract the most contributing parts of the terms in comoving momentum space. The result indicates that they serve to maintain the initial scale invariance over a wide range beyond the comoving Planck scale. This is a challenge toward derivation of the precise primordial spectra, and we expect in the future that it will lead to the resolution of the tensions that have arisen in cosmology.

Avi Friedlander, Ningqiang Song, Aaron C. Vincent

22 pages, 9 figures

The evaporation of primordial black holes provides a promising dark matter production mechanism without relying on any non-gravitational interactions between the dark sector and the Standard Model. In theories of ``Large'' Extra Dimensions (LEDs), the true scale of quantum gravity, $M_*$, could be well below the Planck scale, thus allowing for energetic particle collisions to produce microscopic black holes in the primordial plasma at temperatures as low as $T \gtrsim 100$ GeV. Additionally, LEDs modify the relationship between black hole mass, radius, and temperature, allowing microscopic black holes to grow to macroscopic sizes in the early Universe. In this work we study three scenarios for the production of dark matter via LED black holes: 1) Delayed Evaporating Black Holes (DEBHs) which grow to macroscopic sizes before ultimately evaporating, 2) Instantly Evaporating Black Holes (IEBHs) which immediately evaporate, and 3) stable black hole relics with a mass $M_*$ known as Planckeons. For a given reheating temperature, $T_\mathrm{RH}$, we show that DEBHs produce significantly less dark matter than both IEBHs and Planckeons. IEBHs are able to produce the observed relic abundance of dark matter so long as the reheating scale is in the range $10^{-2} \leq T_\mathrm{RH}/M_* \leq 10^{-1}$. We calculate the average speed for the resulting dark matter and show that it would be sufficiently cold for all dark matter masses $m_{dm} \gtrsim 10^{-4}$ GeV. This mechanism is viable for any scale of quantum gravity in the range $10^4\,\mathrm{ GeV} \leq M_* \leq M_{Pl}$ and for any number of LEDs.

We discuss how string theory, and in particular the "fuzzball" paradigm, has already made and can make meaningful contributions to the phenomenology of strong gravity observations. We outline pertinent research directions for the near-future within this program, and emphasize the unique viewpoints that string theory and fuzzballs bring to phenomenology.