Abstract visibility
Abstract text size

Papers for Friday, Oct 14 2022

Naveen Dukiya, Kuntal Misra, Bikram Pradhan, Vibhore Negi, Bhavya Ailawadhi, Brajesh Kumar, Paul Hickson, Jean Surdej

Accepted for Publication in Journal of Astronomical Instrumentation; 12 pages + references; 12 figures, 2 tables

The International Liquid Mirror Telescope (ILMT) is a 4-meter class survey telescope that has recently achieved first light and is expected to swing into full operations by 1st January 2023. It scans the sky in a fixed 22' wide strip centered at the declination of $+29^o21'41''$ and works in Time Delay Integration (TDI) mode. We present a full catalog of sources in the ILMT strip that can serve as astrometric calibrators. The characteristics of the sources for astrometric calibration are extracted from Gaia EDR3 as it provides a very precise measurement of astrometric properties such as RA ($\alpha$), Dec ($\delta$), parallax ($\pi$), and proper motions ($\mu_{\alpha^{*}}$ & $\mu_{\delta}$). We have crossmatched the Gaia EDR3 with SDSS DR17 and PanSTARRS-1 (PS1) and supplemented the catalog with apparent magnitudes of these sources in g, r, and i filters. We also present a catalog of spectroscopically confirmed white dwarfs with SDSS magnitudes that may serve as photometric calibrators. The catalogs generated are stored in a SQLite database for query-based access. We also report the offsets in equatorial positions compared to Gaia for an astrometrically calibrated TDI frame observed with the ILMT.

B. J. Lynch, B. E. Wood, M. Jin, T. Török, X. Sun, E. Palmerio, R. A. Osten, A. A. Vidotto, O. Cohen, J. D. Alvarado-Gómez, J. J. Drake, V. S. Airapetian, Y. Notsu, A. Veronig, K. Namekata, R. M. Winslow, L. K. Jian, A. Vourlidas, N. Lugaz, N. Al-Haddad, W. B. Manchester, C. Scolini, C. J. Farrugia, E. E. Davies, T. Nieves-Chinchilla, F. Carcaboso, C. O. Lee, T. M. Salman

9 pages, 5 figures, white paper submitted to the Heliophysics 2024--2033 Decadal Survey

The aim of this white paper is to briefly summarize some of the outstanding gaps in the observations and modeling of stellar flares, CMEs, and exoplanetary space weather, and to discuss how the theoretical and computational tools and methods that have been developed in heliophysics can play a critical role in meeting these challenges. The maturity of data-inspired and data-constrained modeling of the Sun-to-Earth space weather chain provides a natural starting point for the development of new, multidisciplinary research and applications to other stars and their exoplanetary systems. Here we present recommendations for future solar CME research to further advance stellar flare and CME studies. These recommendations will require institutional and funding agency support for both fundamental research (e.g. theoretical considerations and idealized eruptive flare/CME numerical modeling) and applied research (e.g. data inspired/constrained modeling and estimating exoplanetary space weather impacts). In short, we recommend continued and expanded support for: (1.) Theoretical and numerical studies of CME initiation and low coronal evolution, including confinement of "failed" eruptions; (2.) Systematic analyses of Sun-as-a-star observations to develop and improve stellar CME detection techniques and alternatives; (3.) Improvements in data-inspired and data-constrained MHD modeling of solar CMEs and their application to stellar systems; and (4.) Encouraging comprehensive solar--stellar research collaborations and conferences through new interdisciplinary and multi-agency/division funding mechanisms.

Danial Langeroodi, Jens Hjorth

8 pages, 5 figures, Submitted to ApJL, Comments are Welcome!

The mid-infrared spectra of star-forming galaxies (SFGs) are characterized by characteristic broad PAH emission features at 3$-$20 $\mu$m. As these features are redshifted, they are predicted to dominate the flux at specific mid-infrared wavelengths, leading to substantial redshift-dependent color variations in broad-band photometry. The advent of JWST for the first time allows the study of this effect for normal SFGs. Based on spectral energy distribution templates, we here present tracks in mid-infrared (4.4, 7.7, 10, 15, and 18 $\mu$m) color-color diagrams describing the redshift dependence of SFG colors. In addition, we present simulated color-color diagrams by populating these tracks using the cosmic star-formation history and the star-formation rate function. Depending on redshift, we find that SFGs stand out in the color-color diagrams by several magnitudes. We provide the first observational demonstration of this effect for galaxies detected in the JWST Early Release Observations of the field towards the lensing cluster SMACS J0723.3$-$7327. While the distribution of detected galaxies is consistent with the simulations, the numbers are substantially boosted by lensing effects. The PAH emitter with the highest spectroscopic redshift, detected in all bands, is a multiply-imaged galaxy at $z=1.45$. There is also a substantial number of cluster members, which do not exhibit PAH emission, except for one Seyfert galaxy at $z=0.38$. Future wider-field observations will further populate mid-infrared color-color diagrams and provide insight into the evolution of typical SFGs.

Yimeng Tang, Bojun Tao, Hong-Xin Zhang, Guangwen Chen, Yulong Gao, Zesen Lin, Yao Yao, Yong Shi, Xu Kong

14 pages, 16 figures. Accepted for publication in A&A

NGC 2915 is a unique nearby galaxy that is classified as an isolated blue compact dwarf based on its optical appearance but has an extremely extended H I gas disk with prominent Sd-type spiral arms. To unveil the starburst-triggering mystery of NGC 2915, we perform a comprehensive analysis of deep VLT/MUSE integral field spectroscopic observations that cover the central kiloparsec star-forming region. We find that episodes of bursty star formation have recurred in different locations throughout the central region, and the most recent one peaked around 50 Myr ago. The bursty star formation has significantly disturbed kinematics of the ionized gas but not the neutral atomic gas, which implies that the two gas phases are largely spatially decoupled along the line of sight. No evidence for an active galactic nucleus is found based on the classical line-ratio diagnostic diagrams. The ionized gas metallicities have a positive radial gradient, which confirms the previous study based on several individual H II regions and may be attributed to both the stellar feedback-driven outflows and metal-poor gas inflow. Evidence for metal-poor gas infall/inflow includes discoveries of high-speed collisions between gas clouds of different metallicities, localized gas metallicity drops and unusually small metallicity differences between gas and stars. The central stellar disk appears to be counter-rotating with respect to the extended H I disk, implying that the recent episodes of bursty star formation have been sustained by externally accreted gas.

D. M. Rowan, T. Jayasinghe, K. Z. Stanek, C. S. Kochanek, Todd A. Thompson, B. J. Shappee, T. W.-S. Holoien, J. L. Prieto, W. Giles

19 pages, 14 figures. Submitted to MNRAS

Detached eclipsing binaries are the primary tool used to measure precise masses and radii of stars. In our previous paper estimating the parameters of more than 30,000 detached eclipsing binaries, we identified 766 eclipsing binaries with additional features in their All-Sky Automated Survey for Supernovae (ASAS-SN) and Transiting Exoplanet Survey Satellite (TESS) light curves. Here, we characterize these "extra-physics" systems, identifying eclipsing binaries with spotted stars, pulsating components, and candidate triple/quadruple systems. We use the Gaia, ATLAS, ZTF, and ASAS-SN variable star catalogs to consider possible blends. We use MIST isochrones and evolutionary tracks to identify systems with main sequence, subgiant, and giant primaries and highlight systems in sparsely populated regions of the color-magnitude diagram. We find that the orbital period distribution of spotted binaries is divided by evolutionary state and find 68 with X-ray detections. For the candidate triple/quadruples and pulsating systems, we calculate the extra orbital/pulsational period and identify systems with resonances. Finally, we highlight a number of exotic systems, including eclipsing CVs, subdwarfs, and binaries with disks.

Cole Johnston, Andrew Tkachenko, Timothy Van Reeth, Dominic M. Bowman, Kresimir Pavlovski, Hugues Sana, Sanjay Sekaran

Accepted for publication in Astronomy & Astrophysics. 13 pages, 10 figures, 3 tables, 1 appendix

Context. Unambiguous examples of the influence of tides on self-excited, free stellar pulsations have recently been observationally detected in space-based photometric data. Aims. We aim to investigate U Gru and contextualise it within the growing class of tidally influenced pulsators. Initial analysis of U Gru revealed frequencies spaced by the orbital frequency that are difficult to explain by currently proposed tidal mechanisms. Methods. We re-investigate the TESS photometry of U Gru alongside new uves spectroscopy. We analyse the uves spectroscopy with least-squares deconvolution and spectral disentangling techniques, and perform an atmospheric analysis. We remove the binary signature from the light curve using an effective model in order to investigate the pulsation signal in the residuals. We track the amplitudes and phases of the residual pulsations as a function of the orbital period to reveal their tidal influence. Results. We establish that U Gru is likely a hierarchical triple system. We identify a single p mode oscillation to exhibit amplitude and phase variation over the binary orbit. We propose a toy model to demonstrate that the series of frequencies separated by the orbital frequency can be reproduced by eclipse mapping. We find no evidence of modulation to the other independent oscillation modes. Conclusions. We demonstrate that U Gru hosts at least one tidally perturbed pulsation. Additionally we argue that eclipse mapping of the dominant, tidally perturbed mode can produce the series of frequencies separated by the observed orbital frequency. Our simulations show that the effects of eclipse mapping are mode dependent, and are not expected to produce an observable signature for all pulsation modes in an eclipse binary.

K. D. Neumann, T. W.-S. Holoien, C. S. Kochanek, K. Z. Stanek, P. J. Vallely, B. J. Shappee, J. L. Prieto, T. Pessi, T. Jayasinghe, J. Brimacombe, D. Bersier, E. Aydi, C. Basinger, J. F. Beacom, S. Bose, J. S. Brown, P. Chen, A. Clocchiatti, D. D. Desai, S. Dong, E. Falco, S. Holmbo, N. Morrell, J. V. Shields, K. Sokolovsky, J. Strader, M. D. Stritzinger, S. Swihart, T. A. Thompson, Z. Way, L. Aslan, D. W. Bishop, G. Bock, J. Bradshaw, P. Cacella, N. Castro, E. Conseil, R. Cornect, I. Cruz, R. G. Farfan, J. M. Fernandez, A. Gabuya, J.-L. Gonzalez-Carballo, M. R. Kendurkar, S. Kiyota, R. A. Koff, G. Krannich, P. Marples, G. Masi, L. A. G. Monard, J. A. Muñoz, B. Nicholls, R. S. Post, Z. Pujic, G. Stone, L. Tomasella, D. L. Trappett, W. S. Wiethoff

14 pages, 7 figures, 4 tables. MNRAS submission in progress. Tables containing the catalog data presented in this submission are included in machine-readable format as ancillary files

We catalog the 443 bright supernovae discovered by the All-Sky Automated Survey for Supernovae (ASAS-SN) in $2018-2020$ along with the 519 supernovae recovered by ASAS-SN and 516 additional $m_{peak}\leq18$ mag supernovae missed by ASAS-SN. Our statistical analysis focuses primarily on the 984 supernovae discovered or recovered in ASAS-SN $g$-band observations. The complete sample of 2427 ASAS-SN supernovae includes earlier $V$-band samples and unrecovered supernovae. For each supernova, we identify the host galaxy, its UV to mid-IR photometry, and the offset of the supernova from the center of the host. Updated light curves, redshifts, classifications, and host galaxy identifications supersede earlier results. With the increase of the limiting magnitude to $g\leq18$ mag, the ASAS-SN sample is roughly complete up to $m_{peak}=16.7$ mag and is $90\%$ complete for $m_{peak}\leq17.0$ mag. This is an increase from the $V$-band sample where it was roughly complete up to $m_{peak}=16.2$ mag and $70\%$ complete for $m_{peak}\leq17.0$ mag.

Brent Tan, S. Peng Oh, Max Gronke

23 pages, 20 figures, submitted to MNRAS

Understanding the survival, growth and dynamics of cold gas is fundamental to galaxy formation. While there has been a plethora of work on `wind tunnel' simulations that study such cold gas in winds, the infall of this gas under gravity is at least equally important, and fundamentally different since cold gas can never entrain. Instead, velocity shear increases and remains unrelenting. If these clouds are growing, they can experience a drag force due to the accretion of low momentum gas, which dominates over ram pressure drag. This leads to sub-virial terminal velocities, in line with observations. We develop simple analytic theory and predictions based on turbulent radiative mixing layers. We test these scalings in 3D hydrodynamic simulations, both for an artificial constant background, as well as a more realistic stratified background. We find that the survival criterion for infalling gas is more stringent than in a wind, requiring that clouds grow faster than they are destroyed ($t_{\rm grow} < 4\,t_{\rm cc} $). This can be translated to a critical pressure, which for Milky Way like conditions is $P \sim 3000 {\rm k}_B {\rm K}\,{\rm cm}^{-3}$ . Cold gas which forms via linear thermal instability ($t_{\rm cool}/t_{\rm ff} < 1$) in planar geometry meets the survival threshold. In stratified environments, larger clouds need only survive infall until cooling becomes effective. We discuss applications to high velocity clouds and filaments in galaxy clusters.

Elham Rahimi, Evan Vienneau, Nassim Bozorgnia, Andrew Robertson

39 pages, 12 figures, 15 tables

We study the effects of dark matter self-interactions on the local dark matter distribution in selected Milky Way-like galaxies in the EAGLE hydrodynamical simulations. The simulations were run with two different self-interacting dark matter models, a constant and velocity-dependent self-interaction cross-section. We find that the local dark matter velocity distribution of the Milky Way-like halos in the simulations with dark matter self-interactions and baryons are generally similar to those extracted from cold collisionless dark matter simulations with baryons. In both cases, the local dark matter speed distributions agree well with their best fit Maxwellian distributions. Including baryons in the simulations with or without dark matter self-interactions increases the local dark matter density and shifts the dark matter speed distributions to higher speeds. To study the implications for direct detection, we compute the dark matter halo integrals obtained directly from the simulations and compare them to those obtained from the best fit Maxwellian velocity distribution. We find that a Maxwellian distribution provides a good fit to the halo integrals of most halos, without any significant difference between the results of different dark matter self-interaction models.

Griffin Hosseinzadeh, David J. Sand, Jacob E. Jencson, Jennifer E. Andrews, Irene Shivaei, K. Azalee Bostroem, Stefano Valenti, Tamás Szalai, Jamison Burke, D. Andrew Howell, Curtis McCully, Megan Newsome, Estefania Padilla Gonzalez, Craig Pellegrino, Giacomo Terreran

submitted to ApJL

We present near- and mid-infrared (0.9-18 $\mu$m) photometry of supernova (SN) 2021afdx, which was imaged serendipitously with the James Webb Space Telescope (JWST) as part of its Early Release Observations of the Cartwheel Galaxy. Our ground-based optical observations show it is likely to be a Type IIb SN, the explosion of a yellow supergiant, and its infrared spectral energy distribution (SED) $\approx$200 days after explosion shows two distinct components, which we attribute to hot gas and warm dust in the SN ejecta. By fitting models of dust emission to the SED, we derive a lower limit on the dust mass of ${>}2.8 \times 10^{-3}\ M_\odot$, which is the highest yet observed in a Type IIb SN but consistent with other Type II SNe observed by the Spitzer Space Telescope. We also find that the radius of the dust is consistent with the radius of the ejecta, as derived from spectroscopic velocities during the photospheric phase, which could imply that the dust formed inside the ejecta. However, we cannot rule out an infrared echo off of pre-existing dust in the progenitor environment. Our results show the power of JWST to address questions of dust formation in SNe, and therefore the presence of dust in the early Universe, with much larger samples than have been previously possible.

Shifra Mandel, Michael M. Shara, David Zurek, Charlie Conroy, Pieter van Dokkum

Accurate determination of the rates of nova eruptions in different kinds of galaxies give us strong constraints on those galaxies' underlying white dwarf and binary populations, and those stars' spatial distributions. Until 2016, limitations inherent in ground-based surveys of external galaxies - and dust extinction in the Milky Way - significantly hampered the determination of those rates and how much they differ between different types of galaxies. Infrared Galactic surveys and dense cadence Hubble Space Telescope (HST)-based surveys are overcoming these limitations, leading to sharply increased nova-in-galaxy rates relative to those previously claimed. Here we present 14 nova candidates that were serendipitously observed during a year-long HST survey of the massive spiral galaxy M51 (the "Whirlpool Galaxy"). We use simulations based on observed nova light curves to model the incompleteness of the HST survey in unprecedented detail, determining a nova detection efficiency $\epsilon = 20.3$ percent. The survey's M51 area coverage, combined with $\epsilon$, indicates a conservative M51 nova rate of $172^{+46}_{-37}$ novae yr$^{-1}$, corresponding to a luminosity-specific nova rate (LSNR) of $\sim10.4^{+2.8}_{-2.2}$ novae yr$^{-1}$/$10^{10} L_{\odot,K}$. Both these rates are approximately an order of magnitude higher than those estimated by ground-based studies, contradicting claims of universal low nova rates in all types of galaxies determined by low cadence, ground-based surveys. They demonstrate that, contrary to theoretical models, the HST-determined LSNR in a giant elliptical galaxy (M87) and a giant spiral galaxy (M51) likely do not differ by an order of magnitude or more, and may in fact be quite similar.

A.C. Trapp, Steven R. Furlanetto, Frederick B. Davies

14 pages, 14 figures, submitted to MNRAS

Lyman-alpha emitters (LAEs) are excellent probes of the reionization process, as they must be surrounded by large ionized bubbles in order to be visible during the reionization era. Large ionized regions are thought to correspond to over-dense regions and may be protoclusters, making them interesting test-beds for early massive structures. Close associations containing several LAEs are often assumed to mark over-dense, ionized bubbles. Here, we develop the first framework to quantify the ionization and density fields of high-z galaxy associations. We explore the interplay between (i) the large-scale density of a survey field, (ii) Poisson noise due to the small number density of bright sources at high redshifts (z~7), and (iii) the effects of the ionized fraction on the observation of LAEs. We use Bayesian statistics, a simple model of reionization, and a Monte-Carlo simulation to construct a more comprehensive method for calculating the large-scale density of LAE regions than previous works. We find that Poisson noise has a strong effect on the inferred density of a region and show how the ionized fraction can be inferred. We then apply our framework to the strongest association yet identified: Hu et al. (2021) found 14 LAEs in a volume of ~50,000 cMpc^3 inside the COSMOS field at z~7. We show that this is most likely a 2.5-sigma over-density inside of an ionized or nearly ionized bubble. We also show that this LAE association implies that the global ionized fraction is Q = 0.60 (+0.08,-0.09), within the context of a simple reionization model.

Jose M. Diego, Ashish K. Meena, Nathan J. Adams, Tom Broadhurst, Liang Dai, Dan Coe, Brenda Frye, Patrick Kelly, Anton M. Koekemoer, Massimo Pascale, S. P. Willner, Erik Zackrisson, Adi Zitrin, Rogier A. Windhorst, Seth H. Cohen, Rolf A. Jansen, Jake Summers, Scott Tompkins, Christopher J. Conselice, Simon P. Driver, Haojing Yan, Norman Grogin, Madeline A. Marshall, Nor Pirzkal, Aaron Robotham, Russell E. Ryan, Jr., Christopher N. A. Willmer, Larry D. Bradley, Gabriel Caminha, Karina Caputi

27 pages, 21 figures

The first JWST data on the massive colliding cluster El Gordo confirm 23 known families of multiply lensed images and identify 8 new members of these families. Based on these families, which have been confirmed spectroscopically by MUSE, we derived an initial lens model. This model guided the identification of 37 additional families of multiply lensed galaxies, among which 28 are entirely new systems, and 9 were previously known. The initial lens model determined geometric redshifts for the 37 new systems. The geometric redshifts agree reasonably well with spectroscopic or photometric redshifts when those are available. The geometric redshifts enable two additional models that include all 60 families of multiply lensed galaxies spanning a redshift range $2<z<6$. The derived dark-matter distribution confirms the double-peak configuration of mass found by earlier work with the southern and northern clumps having similar masses. We confirm that El Gordo is the most massive known cluster at $z>0.8$ and has an estimated virial mass close the maximum mass allowed by standard cosmological models. The JWST images also reveal the presence of small-mass perturbers that produce small lensing distortions. The smallest of these is consistent with being a dwarf galaxy at $z=0.87$ and has an estimated mass of $3.8\times10^9$~\Msol, making it the smallest substructure found at $z>0.5$. The JWST images also show several candidate caustic-crossing events. One of them is detected at high significance at the expected position of the critical curve and is likely a red supergiant star at $z=2.1878$. This would be the first red supergiant found at cosmological distances. The cluster lensing should magnify background objects at $z>6$, making more of them visible than in blank fields of similar size, but there appears to be a deficiency of such objects.

Fengwu Sun, Eiichi Egami

6 pages, 4 figures, accepted for publication in MNRAS Letter

We analysed the high-resolution (up to $\sim$0.2") ALMA CO (2-1) and 1.3 mm dust continuum data of eight gas-rich post-starburst galaxies (PSBs) in the local Universe, six of which had been studied by Smercina et al. (2022). In contrast to this study reporting the detections of extraordinarily compact (i.e., unresolved) reservoirs of molecular gas in the six PSBs, our visibility-plane analysis resolves the CO (2-1) emission in all eight PSBs with effective radii ($R_\mathrm{e,CO}$) of $0.8_{-0.4}^{+0.9}$ kpc, typically consisting of gaseous components at both circumnuclear and extended disc scales. With this new analysis, we find that the CO sizes of gas-rich PSBs are compact with respect to their stellar sizes (median ratio $=0.43_{-0.21}^{+0.27}$), but comparable to the sizes of the gas discs seen in local luminous infrared galaxies (LIRGs) and early-type galaxies. We also find that the CO-to-stellar size ratio of gas-rich PSBs is potentially correlated with the gas depletion time scale, placing them as transitional objects between LIRGs and early-type galaxies from an evolutionary perspective. Finally, the star formation efficiency of the observed PSBs appear consistent with those of star-forming galaxies on the Kennicutt-Schmidt relation, showing no sign of suppressed star formation from turbulent heating.

A. J. Cooper, A. Rowlinson, R. A. M. J. Wijers, C. Bassa, K. Gourdji, J. Hessels, A. J. van der Horst, V. Kondratiev, Z. Pleunis, T. Shimwell, S. ter Veen

Accepted for publication in MNRAS

We present LOFAR imaging observations from the April/May 2020 active episode of magnetar SGR 1935+2154. We place the earliest radio limits on persistent emission following the low-luminosity fast radio burst FRB 200428 from the magnetar. We also perform an image-plane search for transient emission and find no radio flares during our observations. We examine post-FRB radio upper limits in the literature and find that all are consistent with the multi-wavelength afterglow predicted by the synchrotron maser shock model interpretation of FRB 200428. However, early optical observations appear to rule out the simple versions of the afterglow model with constant-density circumburst media. We show that these constraints may be mitigated by adapting the model for a wind-like environment, but only for a limited parameter range. In addition, we suggest that late-time non-thermal particle acceleration occurs within the afterglow model when the shock is no longer relativistic, which may prove vital for detecting afterglows from other Galactic FRBs. We also discuss future observing strategies for verifying either magnetospheric or maser shock FRB models via rapid radio observations of Galactic magnetars and nearby FRBs.

Yinuo Han, Peter G. Tuthill, Ryan M. Lau, Anthony Soulain

Published in Nature

The Wolf-Rayet (WR) binary system WR140 is a close (0.9-16.7 mas) binary star consisting of an O5 primary and WC7 companion and is known as the archetype of episodic dust-producing WRs. Dust in WR binaries is known to form in a confined stream originating from the collision of the two stellar winds, with orbital motion of the binary sculpting the large-scale dust structure into arcs as dust is swept radially outwards. It is understood that sensitive conditions required for dust production in WR140 are only met around periastron when the two stars are sufficiently close. Here we present multiepoch imagery of the circumstellar dust shell of WR140. We constructed geometric models that closely trace the expansion of the intricately structured dust plume, showing that complex effects induced by orbital modulation may result in a 'Goldilocks zone' for dust production. We find that the expansion of the dust plume cannot be reproduced under the assumption of a simple uniform-speed outflow, finding instead the dust to be accelerating. This constitutes a direct kinematic record of dust motion under acceleration by radiation pressure and further highlights the complexity of the physical conditions in colliding-wind binaries.

Paola Pinilla

Invited contribution as an article in the European Physics Journal + (accepted for publication)

Brown dwarfs and very low mass stars are a significant fraction of stars in our galaxy, and they are interesting laboratories to investigate planet formation in extreme conditions of low temperature and densities. In addition, the dust radial drift of particles is expected to be a more difficult barrier to overcome during the first steps of planet formation in these disks. ALMA high-angular resolution observations of few protoplanetary disks around BDs and VLMS have shown substructures as in the disks around Sun-like stars. Such observations suggest that giant planets embedded in the disks are the most likely origin of the observed substructures. However, this type of planets represent less than 2% of the confirmed exoplanets so far around all stars, and they are difficult to form by different core accretion models (either pebble or planetesimal accretion). Dedicated deep observations of disks around BDs and VLMS with ALMA and JWST will provide significant progress on understanding the main properties of these objects (e.g., disk size and mass), which is crucial for determining the physical mechanisms that rule the evolution of these disks and the effect on the potential planets that may form in these environments.

Andreas Seifahrt, Jacob L. Bean, David Kasper, Julian Stürmer, Madison Brady, Robert Liu, Mathias Zechmeister, Gudmundur Stefansson, Ben Montet, John White, Eduardo Tapia, Teo Mocnik, Siyi Xu, Christian Schwab

MAROON-X is a fiber-fed, optical EPRV spectrograph at the 8-m Gemini North Telescope on Mauna Kea, Hawai'i. MAROON-X was commissioned as a visiting instrument in December 2019 and is in regular use since May 2020. Originally designed for RV observations of M-dwarfs, the instrument is used for a broad range of exoplanet and stellar science cases and has transitioned to be the second-most requested instrument on Gemini North over a number of semesters. We report here on the first two years of operations and radial velocity observations. MAROON-X regularly achieves sub-m/s RV performance on sky with a short-term instrumental noise floor at the 30 cm/s level. We will discuss various technical aspects in achieving this level of precision and how to further improve long-term performance

Kunal P. Mooley (Caltech, NRAO), Jay Anderson (STSci), Wenbin Lu (Caltech, Princeton, Berkeley)

Published in Nature main journal (full-text SharedIt link: this https URL ). Jet proper motion data available at: this https URL and this http URL

The afterglow of the binary neutron star merger GW170817 gave evidence for a structured relativistic jet and a link between such mergers and short gamma-ray bursts. Superluminal motion, found using radio very long baseline interferometry (VLBI), together with the afterglow light curve provided constraints on the viewing angle (14-28 degrees), the opening angle of the jet core (less than about 5 degrees), and a modest limit on the initial Lorentz factor of the jet core (more than 4). Here we report on another superluminal motion measurement, at seven times the speed of light, leveraging Hubble Space Telescope precision astrometry and previous radio VLBI data of GW170817. We thereby obtain a unique measurement of the Lorentz factor of the wing of the structured jet, as well as substantially improved constraints on the viewing angle (19-25 degrees) and the initial Lorentz factor of the jet core (more than 40).

Viktória Kecskeméthy, Csaba Kiss, Róbert Szakáts, András Pál, Gyula M. Szabó, László Molnár, Krisztián Sárneczky, József Vinkó, Róbert Szabó, Gábor Marton, Anikó Farkas-Takács, Csilla E. Kalup, László L. Kiss

Accepted for publication in ApJS

The K2 mission of the Kepler Space Telescope allowed the observations of light curves of small solar system bodies throughout the whole Solar system. In this paper we present the results of a collection of K2 transneptunian object observations, between Campaigns C03 (November 2014 -- February 2015) to C19 (August -- September, 2018), which includes 66 targets. Due to the faintness of our targets the detectability rate of a light curve period is $\sim$56%, notably lower than in the case of other small body populations, like Hildas or Jovian trojans. We managed to obtain light curve periods with an acceptable confidence for 37 targets; the majority of these cases are new identifications. We were able to give light curve amplitude upper limits for the other 29 targets. Several of the newly detected light curve periods are longer than $\sim$24 h, in many cases close to $\sim$80 h, i.e., these targets are slow rotators. This relative abundance of slowly rotating objects is similar to that observed among Hildas, Jovian trojans and Centaurs in the K2 mission, and also among main belt asteroids measured with the TESS Space Telescope. Transneptunian objects show notably higher light curve amplitudes at large (D $\gtrsim$ 300 km) sizes than that found among large main belt asteroids, in contrast to the general expectation that due to their lower compressive strength they reach hydrostatic equlibrium at smaller sizes than their inner solar system counterparts.

The intermediate period gap, discovered by Kepler, is an observed dearth of stellar rotation periods in the temperature-period diagram at $\sim$ 20 days for G dwarfs and up to $\sim$ 30 days for early-M dwarfs. However, because Kepler mainly targeted solar-like stars, there is a lack of measured periods for M dwarfs, especially those at the fully convective limit. Therefore it is unclear if the intermediate period gap exists for mid- to late-M dwarfs. Here, we present a period catalog containing 40,553 rotation periods (9,535 periods $>$ 10 days), measured using the Zwicky Transient Facility (ZTF). To measure these periods, we developed a simple pipeline that improves directly on the ZTF archival light curves and reduces the photometric scatter by 26%, on average. This new catalog spans a range of stellar temperatures that connect samples from Kepler with MEarth, a ground-based time domain survey of bright M-dwarfs, and reveals that the intermediate period gap closes at the theoretically predicted location of the fully convective boundary ($G_{\rm BP} - G_{\rm RP} \sim 2.45$ mag). This result supports the hypothesis that the gap is caused by core-envelope interactions. Using gyro-kinematic ages, we also find a potential rapid spin-down of stars across this period gap.

Lizhong Zhang, Omer Blaes, Yan-Fei Jiang

18 pages, 17 figures, and 3 table, submitted to MNRAS

We perform 2D axisymmetric radiative relativistic MHD simulations of radiation pressure supported neutron star accretion columns in split-monopole magnetic fields. The accretion columns exhibit quasi-periodic oscillations, which manifest in the luminosity power spectrum as 2-10 kHz peaks, together with broader extensions to somewhat higher frequencies. The peak frequency decreases for wider columns or higher mass accretion rates. In contrast to the case of shorter columns in uniform magnetic fields, pdV work contributes substantially to maintaining the radiation pressure inside the column against sideways radiative cooling. This is in part due to the compression associated with accretion along the converging magnetic field lines towards the stellar surface. Propagating entropy waves which are associated with the slow-diffusion photon bubble instability form in all our simulations. Radial advection of radiation from the oscillation itself as well as the entropy waves is also important in maintaining radiation pressure inside the column. The time-averaged profile of our fiducial simulation accretion is approximately consistent with the classical 1D stationary model provided one incorporates the correct column shape. We also quantify the porosity in all our accretion column simulations so that this may also in principle be used to improve 1D models.

Luciano Gottardi, Stephen Smith

45 pages, 18 figures. This Chapter will appear in the Section "Detectors for X-ray Astrophysics" (Section Editors: J-W. den Harder, M. Feroci, N. Meidinger) of the "Handbook of X-ray and Gamma-ray Astrophysics" (Editors in chief: C. Bambi and A. Santangelo)

Large arrays of superconducting transition-edge sensor (TES) microcalorimeters are becoming the key technology for future space-based X-ray observatories and ground-based experiments in the fields of astrophysics, laboratory astrophysics, plasma physics, particle physics and material analysis. Thanks to their sharp superconducting-to-normal transition, TESs can achieve very high sensitivity in detecting small temperature changes at very low temperature. TES based X-ray detectors are non-dispersive spectrometers bringing together high resolving power, imaging capability and high-quantum efficiency simultaneously. In this chapter, we highlight the basic principles behind the operation and design of TESs, and their fundamental noise limits. We will further elaborate on the key fundamental physics processes that guide the design and optimization of the detector. We will then describe pulse-processing and important calibration considerations for space flight instruments, before introducing novel multi-pixel TES designs and discussing applications in future X-ray space missions over the coming decades.

Rafael C. R. de Lima, Jonas P. Pereira, Jaziel G. Coelho, Rafael C. Nunes, Paulo E. F. Stecchini, Manuel Castro, Pierre Gomes, Rodrigo R. da Silva, Claudia V. Rodrigues, José C. N. de Araujo, Michał Bejger, Paweł Haensel, J. Leszek Zdunik

23 pages, 10 figures, 2 tables. Comments welcome

3XMM J185246.6+003317 is a slowly rotating soft-gamma repeater (neutron star) in the vicinity of the supernova remnant Kes\,79. So far, observations have only set upper limits to its surface magnetic field and spindown, and there is no estimate for its mass and radius. Using ray-tracing modelling and Bayesian inference for the analysis of several light curves spanning a period of around three weeks, we have found that it may be one of the most massive neutron stars to date. In addition, our analysis suggests a multipolar magnetic field structure with a subcritical field strength and a carbon atmosphere composition. Due to the time-resolution limitation of the available light curves, we estimate the surface magnetic field and the mass to be $\log_{10} (B/{\rm G}) = 11.89^{+0.19}_{-0.93}$ and $M=2.09^{+0.16}_{-0.09}$~$M_{\odot}$ at $1\sigma$ confidence level, while the radius is estimated to be $R=12.02^{+1.44}_{-1.42}$ km at $2\sigma$ confidence level. The robustness of these estimates was verified by simulations, i.e., data injections with known model parameters, and their subsequent recovery. The best-fit model has three small hot spots, two of them in the southern hemisphere. We interpret the above results as due to accretion of supernova layers/interstellar medium onto 3XMM J185246.6+003317 leading to burying and a subsequent re-emergence of the magnetic field, and a carbon atmosphere being formed possibly due to hydrogen/helium diffusive nuclear burning. Finally, we briefly discuss some consequences of our findings for superdense matter constraints.

L. R. Cullinane, A. D. Mackey, G. S. Da Costa, S. E. Koposov, D. Erkal

5 pages, 4 figures; Accepted by MNRAS

We use data from the Magellanic Edges Survey (MagES) in combination with Gaia EDR3 to study the extreme southern outskirts of the Small Magellanic Cloud (SMC), focussing on a field at the eastern end of a long arm-like structure which wraps around the southern periphery of the Large Magellanic Cloud (LMC). Unlike the remainder of this structure, which is thought to be comprised of perturbed LMC disk material, the aggregate properties of the field indicate a clear connection with the SMC. We find evidence for two stellar populations in the field: one having properties consistent with the outskirts of the main SMC body, and the other significantly perturbed. The perturbed population is on average ~0.2 dex more metal-rich, and is located ~7 kpc in front of the dominant population with a total space velocity relative to the SMC centre of ~230 km/s broadly in the direction of the LMC. We speculate on possible origins for this perturbed population, the most plausible of which is that it comprises debris from the inner SMC that has been recently tidally stripped by interactions with the LMC.

Casey Brinkman, Lauren M. Weiss, Fei Dai, Daniel Huber, Edwin S. Kite, Diana Valencia, Jacob L. Bean, Corey Beard, Aida Behmard, Sarah Blunt, Madison Brady, Benjamin Fulton, Steven Giacalone, Andrew W. Howard, Howard Isaacson, David Kasper, Jack Lubin, Mason MacDougall, Joseph M. Akana Murphy, Mykhalo Plotnykov, Alex S. Polanski, Malena Rice, Andreas Seifahrt, Gudmundur Stefansson, Julian Sturmer

Submitted to AJ

TOI-561 is a galactic thick disk star hosting an ultra-short period (0.45 day orbit) planet with a radius of 1.37 R$_{\oplus}$, making it one of the most metal-poor ([Fe/H] = -0.41) and oldest ($\sim$10 Gyr) sites where an Earth-sized planet has been found. We present new simultaneous radial velocity measurements (RVs) from Gemini-N/MAROON-X and Keck/HIRES, which we combined with literature RVs to derive a mass of M$_{b}$=2.24 $\pm$ 0.20 M$_{\oplus}$. We also used two new Sectors of TESS photometry to improve the radius determination, finding R$_{b}$=$1.37 \pm 0.04 R_\oplus$, and confirming that TOI-561 b is one of the lowest-density super-Earths measured to date ($\rho_b$= 4.8 $\pm$ 0.5 g/cm$^{3}$). This density is consistent with an iron-poor rocky composition reflective of the host star's iron and rock-building element abundances; however, it is also consistent with a low-density planet with a volatile envelope. The equilibrium temperature of the planet ($\sim$2300 K) suggests that this envelope would likely be composed of high mean molecular weight species, such as water vapor, carbon dioxide, or silicate vapor, and is likely not primordial. We also demonstrate that the composition determination is sensitive to the choice of stellar parameters, and that further measurements are needed to determine if TOI-561 b is a bare rocky planet, a rocky planet with an optically thin atmosphere, or a rare example of a non-primordial envelope on a planet with a radius smaller than 1.5 R$_{\oplus}$.

L. Aldoroty, L. Wang, P. Hoeflich, J. Yang, N. Suntzeff, G. Aldering, P. Antilogus, C. Aragon, S. Bailey, C. Baltay, S. Bongard, K. Boone, C. Buton, Y. Copin, S. Dixon, D. Fouchez, E. Gangler, R. Gupta, B. Hayden, Mitchell Karmen, A. G. Kim, M. Kowalski, D. Küsters, P.-F. Léget, F. Mondon, J. Nordin, R. Pain, E. Pecontal, R. Pereira, S. Perlmutter, K. A. Ponder, D. Rabinowitz, M. Rigault, D. Rubin, K. Runge, C. Saunders, G. Smadja, N. Suzuki, C. Tao, R. C. Thomas, M. Vincenzi

19 pages, 9 figures

We apply the color-magnitude intercept calibration method (CMAGIC) to the Nearby Supernova Factory SNe Ia spectrophotometric dataset. The currently existing CMAGIC parameters are the slope and intercept of a straight line fit to the first linear region in the color-magnitude diagram, which occurs over a span of approximately 30 days after maximum brightness. We define a new parameter, $\omega_{XY}$, the size of the ``bump'' feature near maximum brightness for arbitrary filters $X$ and $Y$. We find a significant correlation between the slope of the first linear region, $\beta_{XY, 1}$, in the CMAGIC diagram and $\omega_{XY}$. These results may be used to our advantage, as they are less affected by extinction than parameters defined as a function of time. Additionally, $\omega_{XY}$ is computed independently of templates. We find that current empirical templates are successful at reproducing the features described in this work, particularly SALT3, which correctly exhibits the negative correlation between slope and bump size seen in our data. In 1-D simulations, we show that the correlation between the size of the bump feature and $\beta_{XY, 1}$ can be understood as a result of chemical mixing due to large-scale Rayleigh-Taylor instabilities.

Baitian Tang, Jiajun Zhang, Zhiqiang Yan, Zhiyu Zhang, Leticia Carigi, José G. Fernández-Trincado

Resubmitted to A&A after addressing the second referee report

Owing to the recent identification of major substructures in our Milky Way (MW), the astronomical community start to reevaluate the importance of dissolved and existing dwarf galaxies. In this work, we investigate up to 13 elements in 43 giant stars of Sculptor dwarf galaxy (Scl) using high signal-to-noise ratio near-infrared APOGEE spectra. Thanks to the strong feature lines in the near-infrared, high resolution O, Si, and Al abundances are determined at such large group of sample stars for the first time in Scl. By comparing [$\alpha$/Fe] (i.e., O, Mg, Si, Ca, and Ti) of stars from Scl, Sagittarius, and MW, we confirm the general trend that less massive galaxy tends to show lower [$\alpha$/Fe]. The low [Al/Fe] ($\sim -0.5$) in Scl proves it as a trustworthy discriminator to identify stars born in dwarf galaxies (from MW field stars). Chemical evolution model suggests Scl have a top-light initial mass function (IMF), with a high-mass IMF power index $\sim -2.7$, and a minimum SN Ia delay time of $\sim 100$ Myr. Furthermore, linear regression analysis indicates negative radial metallicity gradient and positive radial gradients for [Mg/Fe] and [Ca/Fe], qualitatively agree with the outside-in formation scenario.

Tianyu Zhang, Xiangkun Liu, Chengliang Wei, Guoliang Li, Yu Luo, Xi Kang, Zuhui Fan

17 pages, 17 figures, accepted for publication in ApJ

The galaxy intrinsic alignment (IA) is a dominant source of systematics in weak lensing (WL) studies. In this paper, by employing large simulations with semi-analytical galaxy formation, we investigate the IA effects on WL peak statistics. Different simulated source galaxy samples of different redshift distributions are constructed, where both WL shear and IA signals are included. Convergence reconstruction and peak statistics are then performed for these samples. Our results show that the IA effects on peak abundances mainly consist of two aspects. One is the additional contribution from IA to the shape noise. The other is from the satellite IA that can affect the peak signals from their host clusters significantly. The latter depends on the level of inclusion in a shear sample of the satellite galaxies of the clusters that contribute to WL peaks, and thus is sensitive to the redshift distribution of source galaxies. We pay particular attention to satellite IA and adjust it artificially in the simulations to analyze the dependence of the satellite IA impacts on its strength. This information can potentially be incorporated into the modeling of WL peak abundances, especially for high peaks physically originated from massive clusters of galaxies, and thus to mitigate the IA systematics on the cosmological constraints derived from WL peaks.

Florian Kirchschlager, Franziska D. Schmidt, Mike J. Barlow, Ilse De Looze, Nina S. Sartorio

21 pages, 18 Figures

Dust grains form in the clumpy ejecta of core-collapse supernovae where they are subject to the reverse shock, which is able to disrupt the clumps and destroy the grains. Important dust destruction processes include thermal and kinetic sputtering as well as fragmentation and grain vaporization. In the present study, we focus on the effect of magnetic fields on the destruction processes. We have performed magneto-hydrodynamical simulations using AstroBEAR to model a shock wave interacting with an ejecta clump. The dust transport and destruction fractions are computed using our post-processing code Paperboats in which the acceleration of grains due to the magnetic field and a procedure that allows partial grain vaporization have been newly implemented. For the oxygen-rich supernova remnant Cassiopeia A we found a significantly lower dust survival rate when magnetic fields are aligned perpendicular to the shock direction compared to the non-magnetic case. For a parallel field alignment, the destruction is also enhanced but at a lower level. The survival fractions depend sensitively on the gas density contrast between the clump and the ambient medium and on the grain sizes. For a low-density contrast of $100$, e.g., $5\,$nm silicate grains are completely destroyed while the survival fraction of $1\,\mu$m grains is $86\,$per cent. For a high-density contrast of $1000$, $95\,$per cent of the $5\,$nm grains survive while the survival fraction of $1\,\mu$m grains is $26\,$per cent. Alternative clump sizes or dust materials (carbon) have non-negligible effects on the survival rate but have a lower impact compared to density contrast, magnetic field strength, and grain size.

Solar extreme-ultraviolet (EUV) waves generally refer to large-scale disturbances propagating outward from sites of solar eruptions in EUV imaging observations. Using the recent observations from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO), we report a quasi-periodic wave train propagating outward at an average speed of $\sim$308 km s$^{-1}$. At least five wavefronts can be clearly identified with the period being $\sim$120 s. These wavefronts originate from the coronal loop expansion, which propagates with an apparent speed of $\sim$95 km s$^{-1}$, about 3 times slower than the wave train. In the absence of a strong lateral expansion, these observational results might be explained by the theoretical model of Chen et al. (2002), which predicted that EUV waves may have two components: a faster component that is a fast-mode magnetoacoustic wave or shock wave and a slower apparent front formed as a result of successive stretching of closed magnetic field lines. In this scenario, the wave train and the successive loop expansion we observed likely correspond to the fast and slow components in the model, respectively.

Zhihuan Zhou, Yuhao Mu, Gang Liu, Lixin Xu, Jianbo Lu

12 pages, 7 figures

The Hubble horizon at matter-radiation equality ($k^{-1}_{\rm{eq}}$) and the sound horizon at the last scattering surface ($r_s(z_*)$) provides interesting consistency check for the $\Lambda$CDM model and its extensions. It is well known that the reduction of $r_s$ can be compensated by the increase of $H_0$, while the same is true for the standard rulers $k_{\rm{eq}}$. Adding extra radiational component to the early universe can reduce $k_{\rm{eq}}$. The addition of early dark energy (EDE), however, tends to increase $k_{\rm{eq}}$. We perform $k_{\rm{eq}}$- and $r_s$-based analyses in both the EDE model and the Wess-Zumino Dark Radiation (WZDR) model. In the latter case we find $\Delta H_0 = 0.4$ between the $r_s$- and $k_{\rm{eq}}$-based datasets, while in the former case we find $\Delta H_0 = 1.2$. This result suggests that the dark radiation scenario is more consistent in the fit of the two standard rulers ($k_{\rm{eq}}$ and $r_s$). As a forecast analyses, we fit the two models with a mock $k_{\rm{eq}}$ prior derived from \emph{Planck} best-fit $\Lambda$CDM model. Compared with the best-fit $H_0$ in baseline $\Lambda$CDM model, we find $\Delta H_0 = 1.1$ for WZDR model and $\Delta H_0 = - 2.4$ for EDE model.

T. Azevedo Silva (1 and 2), O. D. S. Demangeon (1 and 2), N. C. Santos (1 and 2), R. Allart (3 and 4), F. Borsa (5), E. Cristo (1 and 2), E. Esparza-Borges (6 and 7), J. V. Seidel (8), E. Palle (9), S. G. Sousa (1), H. M. Tabernero (10), M. R. Zapatero Osorio (10), S. Cristiani (11), F. Pepe (4), R. Rebolo (6 and 7), V. Adibekyan (1 and 2), Y. Alibert (12), S. C. C. Barros (1 and 2), F. Bouchy (4), V. Bourrier (4), G. Lo Curto (8), P. Di Marcantonio (11), V. D'Odorico (11 and 13 and 14), D. Ehrenreich (4 and 15), P. Figueira (4 and 1), J. I. González Hernández (6 and 16), C. Lovis (4), C. J. A. P. Martins (1 and 17), A. Mehner (8), G. Micela (18), P. Molaro (11 and 19), D. Mounzer (4), N. J. Nunes (20 and 21), A. Sozzetti (22), A. Suárez Mascareño (6 and 7), S. Udry (4) ((1) Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP Porto, Portugal, (2) Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto, Porto, Portugal, (3) Department of Physics, and Institute for Research on Exoplanets, Université de Montréal, Montréal, Canada, (4) Observatoire astronomique de l'Université de Genève, Versoix, Switzerland, (5) INAF - Osservatorio Astronomico di Brera, Merate, Italy, (6) Instituto de Astrofísica de Canarias (IAC), Tenerife, Spain, (7) Departamento de Astrofísica, Universidad de La Laguna, Tenerife, Spain, (8) European Southern Observatory, Vitacura, Región Metropolitana, Chile, (9) Instituto de Astrofísica de Canarias, Láctea, Spain, (10) Centro de Astrobiología (CSIC-INTA), Madrid, Spain, (11) INAF - Osservatorio Astronomico di Trieste, Trieste, Italy, (12) Physikalisches Institut, University of Bern, Bern, Switzerland, (13) Scuola Normale Superiore, Pisa, Italy, (14) Institute for Fundamental Physics of the Universe, Trieste, Italy, (15) Centre Vie dans l'Univers, Faculté des sciences de l'Université de Genève, Geneva, Switzerland, (16) Universidad de La Laguna (ULL), Departamento de Astrofísica, Tenerife, Spain, (17) Centro de Astrofísica da Universidade do Porto, Porto, Portugal, (18) INAF - Osservatorio Astronomico di Palermo, Palermo, Italy, (19) Institute for Fundamental Physics of the Universe, Miramare, Trieste, Italy, (20) Instituto de Astrofísica e Ciências do Espaço, Faculdade de Ciências da Universidade de Lisboa, (21) Departamento de Física, Faculdade de Ciências da Universidade de Lisboa, Lisboa, Portugal, (22) INAF - Osservatorio Astrofisico di Torino, Pino Torinese (TO), Italy)

21 pages, 16 figures. Published in A&A Volume 666, October 2022 - Letter to the Editor

High-resolution spectroscopy studies of ultra-hot Jupiters have been key in our understanding of exoplanet atmospheres. Observing into the atmospheres of these giant planets allows for direct constraints on their atmospheric compositions and dynamics while laying the groundwork for new research regarding their formation and evolution environments. Two of the most well-studied ultra-hot Jupiters are WASP-76b and WASP-121b, with multiple detected chemical species and strong signatures of their atmospheric dynamics. We take a new look at these two exceptional ultra-hot Jupiters by reanalyzing the transit observations taken with ESPRESSO at the Very Large Telescope and attempt to detect additional species. To extract the planetary spectra of the two targets, we corrected for the telluric absorption and removed the stellar spectrum contributions. We then exploited new synthetic templates that were specifically designed for ultra-hot Jupiters in combination with the cross-correlation technique to unveil species that remained undetected by previous analyses. We add a novel detection of Ba+ to the known atmospheric compositions of WASP-76b and WASP-121b, the heaviest species detected to date in any exoplanetary atmosphere, with additional new detections of Co and Sr+ and a tentative detection of Ti+ for WASP-121b. We also confirm the presence of Ca+, Cr, Fe, H, Li, Mg, Mn, Na, and V on both WASP-76b and WASP-121b, with the addition of Ca, Fe+, and Ni for the latter. Finally, we also confirm the clear asymmetric absorption feature of Ca+ on WASP-121b, with an excess absorption at the bluer wavelengths and an effective planet radius beyond the Roche lobe. This indicates that the signal may arise from the escape of planetary atmosphere.

Masanobu Kunitomo, Tristan Guillot, Gaël Buldgen

Accepted for publication in A&A. 10 pages, 7 figures. Supplemental materials are available at this https URL

Solar evolutionary models are thus far unable to reproduce spectroscopic, helioseismic, and neutrino constraints consistently, resulting in the so-called solar modeling problem. In parallel, planet formation models predict that the evolving composition of the protosolar disk and, thus, of the gas accreted by the proto-Sun must have been variable. We show that solar evolutionary models that include a realistic planet formation scenario lead to an increased core metallicity of up to 5%, implying that accurate neutrino flux measurements are sensitive to the initial stages of the formation of the Solar System. Models with homogeneous accretion match neutrino constraints to no better than 2.7$\sigma$. In contrast, accretion with a variable composition due to planet formation processes, leading to metal-poor accretion of the last $\sim$4% of the young Sun's total mass, yields solar models within 1.3$\sigma$ of all neutrino constraints. We thus demonstrate that in addition to increased opacities at the base of the convective envelope, the formation history of the Solar System constitutes a key element in resolving the current crisis of solar models.

J. Gomes da Silva (IA), A. Bensabat (IA), T. Monteiro (IA), N. C. Santos (IA)

25 pages, 14 figures, 3 tables, accepted for publication in A&A

In this work we investigate the effect of the H$\alpha$ bandpass width in the correlation between the CaII H&K and H$\alpha$ indices with the aim of improving the H$\alpha$ index to better identify and model the signals coming from activity variability. We used a sample of 152 FGK dwarfs observed with HARPS for more than 13 years with enough cadence to be able to detect rotational modulations and cycles in activity proxies. We calculated the CaII H&K and H$\alpha$ activity indices using a range of bandwidths for H$\alpha$ between 0.1 and 2.0 Ang. We studied the correlation between the indices time series at long and short timescales and analysed the impact of stellar parameters, activity level and variability on the correlations. The correlation between CaII H&K and H$\alpha$ both at short and long timespans is maximised when using narrow H$\alpha$ bandwidths, with a maximum at 0.6 Ang. For some inactive stars, as the activity level increases, the flux in the H$\alpha$ line core increases while the flux in the line wings decreases as the line becomes shallower and broader. The balance between these fluxes can cause stars to show the negative correlations observed in the literature when using a wide bandwidth on H$\alpha$. These anti-correlations may become positive correlations if using the 0.6 Ang bandwidth. Calculating the H$\alpha$ index using a bandpass of 0.6 Ang maximises the correlation between CaII H&K and H$\alpha$ both at short and long timescales. On the other hand, the use of the broader 1.6 Ang, generally used in exoplanet detection to identify stellar activity signals, degrades the signal by including the flux in the line wings. In face to these results we strongly recommend the use of a 0.6 Ang bandwidth when computing the H$\alpha$ index for the identification of activity rotational modulation and magnetic cycle signals in solar-type stars. (Abridged)

Chen Jiang

13 pages, 11 figures, Accepted for publication by Astronomische Nachrichten

The advent of ultra-precise photometry space missions enable the possibility of investigating stellar interior with mixed modes. The structural variations induced by the discontinuity of the chemical composition left behind during the first dredge--up is an important feature in the stellar mid-layers located between the hydrogen-burning shell and the base of the convective zone of red giants, as the mixed-mode properties can be significantly affected by these variations. In this paper, the contributing factors to variations of the mixed-mode coupling factor, $q$, are discussed with stellar models. In general, the structural variations give rise to a subtle displacement in the Lamb frequency and a sharp change in the buoyancy frequency, which lead to variations in the value of $q$ computed using the asymptotic formalisms that assuming a smooth background free of structural variations. The impact of these two factors can be felt in detectable mixed modes in low-luminosity red giants. Furthermore, the different nature of variations of the two characteristic frequencies with radius near the base of the convective zone, produces a sudden increase in $q$ in evolved red giants. This is followed by a quick drop in $q$ as the star evolves further along the red giant branch.

Luciano Gottardi, Kenichiro Nagayoshi (NWO-I/SRON Netherlands Institute for Space Research, The Netherlands)

48 pages, 14 figures, published in Appl. Sci. 2021, 11, 3793

The state-of-the-art technology of X-ray microcalorimeters based on superconducting transition edge sensors (TESs), for applications in astrophysics and particle physics, is reviewed. We will show the advance in understanding the detector physics and describe the recent breakthroughs in the TES design that are opening the way towards the fabrication and the read-out of very large arrays of pixels with unprecedented energy resolution. The most challenging low temperature instruments for space- and ground-base experiments will be described.

Bibhuti Kumar Jha, Manjunath Hegde, Aditya Priyadarshi, Sudip Mandal, B Ravindra, Dipankar Banerjee

13 Pages, 7 Figures, Accepted for Publication in Frontiers in Astronomy and Space Sciences Stellar and Solar Physics

Kodaikanal Solar Observatory (KoSO) possesses one of world's longest and homogeneous records of sunspot observations that span more than a century (1904-2017). Interestingly, these observations (originally recorded in photographic plates/films) were taken with the same setup over this entire time period which makes this data unique and best suitable for long-term solar variability studies. A large part of this data, between 1921-2011, were digitized earlier and a catalog containing the detected sunspot parameters (e.g., area and location) was published in Mandal et al.(2017). In this article, we extend the earlier catalog by including new sets of data between 1904-1921 and 2011-2017. To this end, we digitize and calibrate these new datasets which include resolving the issue of random image orientation. We fix this by comparing the KoSO images with co-temporal data from Royal Greenwich Observatory. Following that, a semi-automated sunspot detection and automated umbra detection algorithm are implemented onto these calibrated images to detect sunspots and umbra. Additionally, during this catalog updation, we also filled data gaps in the existing KoSO sunspot catalog (1921-2011) by virtue of re-calibrating the 'rouge' plates. This updated sunspot area series covering nearly 115 years (1904-2017) are being made available to the community and will be a unique source to study the long term variability of the Sun

Sahana Kumar, Eric Y. Hsiao, Chris Ashall, Mark M. Phillips, Nidia Morrell, Peter Hoeflich, Chris R. Burns, Lluis Galbany, Eddie Baron, Carlos Contreras, Scott Davis, Tiara Diamond, Francisco Forster, Melissa L. Graham, Emir Karamehmetoglu, Robert P. Kirshner, Baerbel Koribalski, Kevin Krisciunas, Jing Lu, G. H. Marion, Priscila J. Pessi, Anthony L. Piro, Melissa Shahbandeh, Maximillian D Stritzinger, Nicholas B. Suntzeff, Syed A. Uddin

submitted to ApJ on Oct 4, 2022

We present multi-wavelength time-series spectroscopy of SN 2013aa and SN 2017cbv, two Type Ia supernovae (SNe Ia) on the outskirts of the same host galaxy, NGC 5643. This work utilizes new nebular-phase near-infrared (NIR) spectra obtained by the Carnegie Supernova Project-II, in addition to previously published optical and NIR spectra. By measuring nebular-phase [Fe II] lines in both the optical and NIR, we examine the explosion kinematics and test the efficacy of several emission line fitting techniques commonly used in the literature. The NIR [Fe II] 1.644 $\mu$m line provides the most robust velocity measurements against variations due to the choice of the fit method and line blending. The resulting effects on velocity measurements due to choosing different fit methods, initial fit parameters, continuum and line profile functions, and fit region boundaries were also investigated. The NIR [Fe II] velocities yield the same radial shift direction as velocities measured using the optical [Fe II] 7155 A line, but the sizes of the shifts are consistently and substantially lower, pointing to a potential issue in optical studies. The NIR [Fe II] 1.644 $\mu$m emission profile shows a lack of significant asymmetry in both SNe Ia, and the observed low velocities elevate the importance for correcting for any radial velocity contribution from the host galaxy's rotation. The low [Fe II] velocities measured in the NIR at nebular phases disfavors most progenitor scenarios in close double-degenerate systems for both SN 2013aa and SN 2017cbv. The time evolution of the NIR [Fe II] 1.644 $\mu$m line also indicates moderately high progenitor white dwarf central density and potentially high magnetic fields. These sibling SNe Ia were well observed at both early and late times, providing an excellent opportunity to study the intrinsic diversity of SNe Ia.

Marcelo M. Miller Bertolami

Accepted for Publication in The Astrophysical Journal. 26 Pages. 15 Figures. 5 appendixes

In spite of the spectacular progress accomplished by stellar evolution theory some simple questions remain unanswered. One of these questions is ``Why do stars become Red Giants?''. Here we present a relatively simple analytical answer to this question. We validate our analysis by constructing a quantitative toy-model of a red giant and comparing its predictions to full stellar evolutionar models. We find that the envelope forces the value of $\nabla=d \ln T/d \ln P$ at, and above, the burning shell into a very narrow range of possible values. Together with the fact that the stellar material at the burning shell both provides and transports most of the stellar luminosity, this leads to tight relations between the thermodynamic variables at the burning shell and the mass and radius of the core -- $T_s(M_c,R_s)$, $P_s(M_c,R_s)$ and $\rho_s(M_c,R_s)$. When complemented by typical mass-radius relations of the helium cores, this implies that for all stellar masses the evolution of the core dictates the values of $T_s$, $P_s$ and $\rho_s$. We show that for all stellar masses evolution leads to an increase in the pressure and density contrasts between the shell and the core, forcing a huge expansion of the layers on top of the burning shell. Besides explaining why stars become red giants our analysis also offers a mathematical demonstration of the so-called shell homology relations, and provides simple quantitative answers to some properties of low-mass red giants.

In this paper, we present a spherical Fast Multipole Method (sFMM) for ray tracing simulation of gravitational lensing (GL) on a curved sky. The sFMM is a non-trivial extension of the Fast Multiple Method (FMM) to sphere $\mathbb S^2$, and it can accurately solve the Poisson equation with time complexity of $O(N)\log(N)$, where $N$ is the number of particles. It is found that the time complexity of the sFMM is near $O(N)$ and the computational accuracy can reach $10^{-10}$ in our test. In addition, comparing with the Fast Spherical Harmonic Transform (FSHT), the sFMM is not only faster but more accurate, as it has the ability to reserve high frequency component of the density field. These merits make the sFMM an optimum method to simulation the gravitational lensing on a curved sky, which is the case for upcoming large-area sky surveys, such as the Vera Rubin Observatory and the China Space Station Telescope.

Matteo Braglia, Xingang Chen, Dhiraj Kumar Hazra, Lucas Pinol

36 pages, 15 figures

Future Cosmic Microwave Background (CMB) experiments will deliver extremely accurate measurements of the E-modes pattern of the CMB polarization field. Given the sharpness of the E-modes transfer functions, such surveys make for a powerful detector of high-frequency signals from primordial features that may be lurking in current data sets. With a handful of toy models that increase the fit to the latest Planck data, but are of marginal statistical significance, we use a state-of-the-art forecast pipeline to illustrate the promising prospects to test primordial features in the next decade. Not only will future experiments allow us to detect such features in data, but they will also be able to discriminate between models and narrow down the physical mechanism originating them with high statistical significance. On the other hand, if the anomalies in the currently measured CMB spectra are just statistical fluctuations, all the current feature best fit candidates will be ruled out. Either way, our results show that primordial features are a clear target of forthcoming CMB surveys beyond the detection of tensor modes.

Misao Sasaki, Valeri Vardanyan, Vicharit Yingcharoenrat

10 pages, 6 figures

We propose a novel mechanism for significantly enhancing the amplitude of primordial electromagnetic fields during inflation. Similar to existing proposals, our idea is based on parametric resonance effects due to conformal-symmetry-breaking coupling of the gauge field and the inflaton. Our proposed scenario, however, significantly differs from previously studied models, and avoids their shortcomings. We, particularly, construct a viable system where the gauge field is exponentially amplified on super-horizon scales, therefore evading the no-go theorem formulated on the basis of widely encountered drastic back-reaction of the magnetic field energy on the inflationary background. We compute the spectrum of the produced magnetic fields and demonstrate the compatibility with current observational constraints.

Akash Biswas, Bidya Binay Karak, Robert Cameron

Accepted in Phys.Rev.Lett. (PRL)

A striking feature of the solar cycle is that at the beginning, sunspots appear around mid-latitudes, and over time the latitudes of emergences migrate towards the equator.The maximum level of activity (e.g., sunspot number) varies from cycle to cycle.For strong cycles, the activity begins early and at higher latitudes with wider sunspot distributions than for weak cycles. The activity and the width of sunspot belts increase rapidly and begin to decline when the belts are still at high latitudes. Surprisingly, it has been reported that in the late stages of the cycle the level of activity (sunspot number) as well as the widths and centers of the butterfly wings all have the same statistical properties independent of how strong the cycle was during its rise and maximum phases.We have modeled these features using a Babcock--Leighton type dynamo model and show that the flux loss through magnetic buoyancy is an essential nonlinearity in the solar dynamo.Our study shows that the nonlinearity is effective if the flux emergence becomes efficient at the mean-field strength of the order of $10^4$~G in the lower part of the convection zone.

Maciej Kierkla, Alexandros Karam, Bogumila Swiezewska

We present an updated analysis of the first-order phase transition associated with symmetry breaking in the early Universe in a classically scale-invariant model extended with a new SU(2) gauge group. Including recent developments in understanding supercooled phase transitions, we compute all of its characteristics and significantly constrain the parameter space. We then predict gravitational wave spectra generated during this phase transition and by computing the signal-to-noise ratio we conclude that this model is well testable (and falsifiable) with LISA. We also provide predictions for the relic dark matter abundance. It is consistent with observations in a rather narrow part of the parameter space, since we exclude the so-called supercool dark matter scenario based on an improved description of percolation and reheating after the phase transition as well as inclusion of the running of couplings. Finally, we devote attention to renormalisation-scale dependence of the results. Even though our main results are obtained with the use of renormalisation-group improved effective potential, we also perform a fixed-scale analysis which proves that the dependence on the scale is not only qualitative but also quantitative.

Martin Lopez-Corredoira

accepted to be published in MNRAS

Statistical analyses of the measurements of the Hubble-Lema\^itre constant $H_0$ (163 measurements between 1976 and 2019) show that the statistical error bars associated with the observed parameter measurements have been underestimated -- or the systematic errors were not properly taken into account -- in at least 15-20\% of the measurements. The fact that the underestimation of error bars for $H_0$ is so common might explain the apparent discrepancy of values, which is formally known today as the Hubble tension. Here we have carried out a recalibration of the probabilities with this sample of measurements. We find that $x\sigma $ deviation is indeed equivalent in a normal distribution to $x_{\rm eq.}\sigma $s deviation in the frequency of values, where $x_{\rm eq.}=0.83x^{0.62}$. Hence, a tension of 4.4$\sigma $, estimated between the local Cepheid-supernova distance ladder and cosmic microwave background (CMB) data, is indeed a 2.1$\sigma $ tension in equivalent terms of a normal distribution of frequencies, with an associated probability $P(>x_{\rm eq.})=0.036$ (1 in 28). This can be increased up to a equivalent tension of 2.5$\sigma $ in the worst of the cases of claimed 6$\sigma $ tension, which may anyway happen as a random statistical fluctuation.

Danny C. Price, Cheuk-Yu Edward Tong, Adrian T. Sutinjo, Lincoln J. Greenhill, Nipanjana Patra

IEEE Transactions on Microwave Theory and Techniques

Noise parameters are a set of four measurable quantities which determine the noise performance of a radio-frequency device under test. The noise parameters of a 2-port device can be extracted by connecting a set of 4 or more source impedances at the device's input, measuring the noise power of the device with each source connected, and then solving a matrix equation. However, sources with high reflection coefficients cannot be used due to a singularity that arises in entries of the matrix. Here, we detail a new method of noise parameter extraction using a singularity-free matrix that is compatible with high-reflection sources. We show that open, short, load and an open cable ("OSLC") can be used to extract noise parameters, and we detail a practical measurement approach. The OSLC approach is particularly well-suited for low-noise amplifier measurement at frequencies below 1 GHz, where alternative methods require physically large apparatus.

Searching for extrasolar biosignatures is important to understand life on Earth and its origin. Astronomical observations of exoplanets may find such signatures, but it is difficult and may be impossible to claim unambiguous detection of life by remote sensing of exoplanet atmospheres. Here, another approach is considered: collecting grains ejected by asteroid impacts from habitable exoplanets in the Milky Way and then traveling to the Solar System. The optimal grain size for this purpose is around 1 $\mu$m, and about $10^5$ such grains are expected to be accreting on Earth every year, which may contain biosignatures of life that existed on their home planets. These grains may be collected by detectors placed in space, or extracted from Antarctic ice or deep-sea sediments, depending on future technological developments. In the foreseeable future, this is probably the only approach for humankind to search for extrasolar biosignatures by directly sampling biological materials.

Björn Eichmann, Michael Kachelrieß

The observational data on ultrahigh energy cosmic rays (UHECR), in particular their mass composition, show strong indications for extremely hard spectra of individual mass groups of CR nuclei at Earth. In this work, we show that such hard spectra can be the result of the finite life-time of UHECR sources, if a few individual sources dominate the UHECR flux at the highest energies. In this case, time delays induced by deflections in the turbulent extragalactic magnetic field as well as from the diffusive or advective escape from the source environment can suppress low-energy CRs, leading to a steepening of the observed spectrum. Considering radio galaxies as the main source of UHECRs, we discuss the necessary conditions that few individual sources dominate over the total contribution from the bulk of sources that have been active in the past. We provide two proof-of-principle scenarios showing that for a turbulent extragalactic magnetic field with a strength $B$ and a coherence length $l_{\rm coh}$, the life-time of a source at a distance $d_{\rm src}$ should satisfy ${t_{\rm act} \sim \left( B/1\,\text{nG} \right)^2\,\left( d_{\rm src}/10\,\text{Mpc} \right)^2\,\left( l_{\rm coh}/1\,\text{Mpc} \right)\,\text{Myr}}$ to obtain the necessary hardening of the CR spectrum at Earth.

Toshiya Ueta, Masaaki Otsuka

5 pages (main text with 2 tables) plus 5 supplemental pages (5 tables in Appendix); to be published in Astronomy and Astrophysics Letter to the Editor

Plasma diagnostics are the bases of investigation into the physical and chemical properties of line-emitting gaseous systems. To perform plasma diagnostics properly, it is essential to correct the input spectrum for extinction properly. This is simply because determining the degree of extinction is dependent on the physical properties of the line-emitting gas. Hence, both extinction correction and plasma diagnostics have to be performed simultaneously and self-consistently. By comparing the results of analyses performed for a sample of nine bright planetary nebulae in M 31 with and without the proper extinction correction and plasma diagnostics, we demonstrate how initial assumptions for the physical conditions of the line-emitting gas in extinction correction would compromise the results of the entire analyses. While the electron density/temperature are relatively immune to the imposed inconsistent assumptions, the compromised extinction would cause systematic offsets in the extinction-corrected line strengths, which consequently would impose adverse effects on the resulting ionic and elemental abundances, and other inferences made from the incorrect results. We find that this M 31 PN sample simply represents those around the high-mass end of the mass range for low-mass planetary nebula progenitor stars as expected from the existing theoretical models. It appears that the suspicion raised in the previous study - these PNe being anomalously nitrogen overabundant - is simply caused by the apparent underestimate in extinction that originates from the imposed inconsistent assumptions in extinction correction. In a larger context, the results of plasma diagnostics in the literature without seeking simultaneous self-consistency with extinction correction have to be handled cautiously. Ideally, such previous results should be re-evaluated by seeking simultaneous self-consistency.

Daniel C. M Palumbo, Zachary Gelles, Paul Tiede, Dominic O. Chang, Dominic W. Pesce, Andrew Chael, Michael D. Johnson

Accepted to ApJ, 16 pages, 10 figures

The Event Horizon Telescope (EHT) has produced images of two supermassive black holes, Messier~87* (M 87*) and Sagittarius~A* (Sgr A*). The EHT collaboration used these images to indirectly constrain black hole parameters by calibrating measurements of the sky-plane emission morphology to images of general relativistic magnetohydrodynamic (GRMHD) simulations. Here, we develop a model for directly constraining the black hole mass, spin, and inclination through signatures of lensing, redshift, and frame dragging, while simultaneously marginalizing over the unknown accretion and emission properties. By assuming optically thin, axisymmetric, equatorial emission near the black hole, our model gains orders of magnitude in speed over similar approaches that require radiative transfer. Using 2017 EHT M 87* baseline coverage, we use fits of the model to itself to show that the data are insufficient to demonstrate existence of the photon ring. We then survey time-averaged GRMHD simulations fitting EHT-like data, and find that our model is best-suited to fitting magnetically arrested disks, which are the favored class of simulations for both M 87* and Sgr A*. For these simulations, the best-fit model parameters are within ${\sim}10\%$ of the true mass and within ${\sim}10^\circ$ for inclination. With 2017 EHT coverage and 1\% fractional uncertainty on amplitudes, spin is unconstrained. Accurate inference of spin axis position angle depends strongly on spin and electron temperature. Our results show the promise of directly constraining black hole spacetimes with interferometric data, but they also show that nearly identical images permit large differences in black hole properties, highlighting degeneracies between the plasma properties, spacetime, and most crucially, the unknown emission geometry when studying lensed accretion flow images at a single frequency.

Maria Lomaeva, Ilse De Looze, Amélie Saintonge, Marjorie Decleir

15 pages, 9 figures, accepted to MNRAS

Star formation histories (SFHs) are integral to our understanding of galaxy evolution. We can study recent SFHs by comparing the star formation rate (SFR) calculated using different tracers, as each probes a different timescale. We aim to calibrate a proxy for the present-day rate of change in SFR, dSFR/dt, which does not require full spectral energy distribution (SED) modeling and depends on as few observables as possible, to guarantee its broad applicability. To achieve this, we create a set of models in CIGALE and define an SFR change diagnostic as the ratio of the SFR averaged over the past 5 and 200 Myr, <SFR5>/<SFR200>, probed by the H$\alpha$-FUV colour. We apply <SFR5>/<SFR200> to the nearby spiral NGC 628 and find that its star formation activity has overall been declining in the recent past, with the spiral arms, however, maintaining a higher level of activity. The impact of the spiral arm structure is observed to be stronger on <SFR5>/<SFR200> than on the star formation efficiency (SFE$_\text{H$_2$}$). In addition, increasing disk pressure tends to increase recent star formation, and consequently <SFR5>/<SFR200>. We conclude that <SFR5>/<SFR200> is sensitive to the molecular gas content, spiral arm structure, and disk pressure. The <SFR5>/<SFR200> indicator is general and can be used to reconstruct the recent SFH of any star-forming galaxy for which H$\alpha$, FUV, and either mid- or far-IR photometry is available, without the need of detailed modeling.

G. Banyard, L. Mahy, H. Sana, J. Bodensteiner, J. I. Villaseńor, K. Sen, N. Langer, S. de Mink, A. Picco, T. Shenar

Submitted to A&A

Recent evolutionary computations predict that a few percent of massive OB stars in binary systems should have a dormant BH companion. Despite several reported X-ray quiet OB+BH systems over the last couple of years, finding them with certainty remains challenging. These have great importance as they can be gravitational wave (GW) source progenitors, and are landmark systems in constraining supernova kick physics. This work aims to characterise the hidden companions to the single-lined spectroscopic binaries (SB1s) in the B star population of the young open Galactic cluster NGC 6231 to find candidate systems for harbouring compact object companions. With the orbital solutions for each SB1 previously constrained, we applied Fourier spectral disentangling to multi-epoch optical VLT/FLAMES spectra of each target to extract a potential signature of a faint companion, and to identify newly disentangled double-lined spectroscopic binaries (SB2s). For targets where the disentangling does not reveal any signature of a stellar companion, we performed atmospheric and evolutionary modelling on the primary to obtain constraints on the unseen companion. Seven newly classified SB2 systems with mass ratios down to near 0.1 were identified. From the remaining targets, for which no faint companion could be extracted from the spectra, four are found to have companion masses in the predicted mass ranges of neutron stars (NSes) and BHes. Two of these have companion masses between 1 and 3.5 $M_{\odot}$, making them potential hosts of NSes (or lower mass main sequence stars). The other two are between 2.5 to 8 $M_{\odot}$ and 1.6 and 26 $M_{\odot}$, respectively, and so are identified as candidates for harbouring BH companions. However, unambiguous identification of these systems as X-ray quiet compact object harbouring binaries requires follow up observations.

Kyle Boone, Matthew McQuinn

34 pages in preprint format; 3 figures; to be submitted to JCAP; comments welcome!

The light from a source at a distance d will arrive at detectors separated by 100 AU at times that differ by as much as 120 (d/100 Mpc)^{-1} nanoseconds because of the curvature of the wavefront. At gigahertz frequencies, the arrival time difference can be determined to better than a nanosecond with interferometry. If the space-time positions of the detectors are known to a few centimeters, comparable to the accuracy to which very long baseline interferometry baselines and global navigation satellite systems (GNSS) geolocations are constrained, nanosecond timing would allow competitive cosmological constraints. We show that a four-detector constellation at Solar radii of >10 AU could measure distances to individual sources with sub-percent precision and, hence, cosmological parameters such as the Hubble constant to this precision. FRBs are the only known bright extragalactic radio source that are sufficiently point-like. Galactic scattering limits the timing precision at <5 GHz, whereas at higher frequencies the precision is set by removing dispersion. Furthermore, for baselines greater than 100 AU, Shapiro time delays limit the precision, but their effect can be cleaned with two additional detectors. Accelerations that result in ~1 cm uncertainty in detector positions (from variations in the Sun's irradiance, dust collisions and gaseous drag) could be corrected for with weekly GNSS-like trilaterations. Gravitational accelerations from asteroids occur over longer timescales, and so a setup with a precise accelerometer and calibrating the space-time detector positions off of distant FRBs may also be sufficient. The proposed interferometer would also resolve the radio emission region of Galactic pulsars, constrain the mass distribution in the outer Solar System, and reach interesting sensitivities to ~0.01-100 micro-Hz gravitational waves.

Ze-Wei Zhao, Ling-Feng Wang, Ji-Guo Zhang, Jing-Fei Zhang, Xin Zhang

16 pages, 6 figures

Interacting dark energy (IDE) scenario assumes that there exists a direct interaction between dark energy and cold dark matter, but this interaction is hard to be tightly constrained by the current data. Finding new cosmological probes to precisely measure this interaction could deepen our understanding of dark energy and dark matter. Fast radio bursts (FRBs) will be seen in large numbers by future radio telescopes, and thus they have potential to become a promising low-redshift cosmological probe. In this work, we investigate the capability of future FRBs of constraining the dimensionless coupling parameter $\beta$ in four phenomenological IDE models. We find that in the IDE models with the interaction proportional to the energy density of dark energy, about $10^5$ FRB data can give constraint on $\beta$ tighter than the current cosmic microwave background data. In all the IDE models, about $10^6$ FRB data can constrain the absolute errors of $\beta$ to less than $0.10$, providing a way to precisely measure $\beta$ by only one cosmological probe. The reconstruction of the interaction term also shows that the FRB data could help constrain the redshift evolution of interaction.

Shun-Sheng Li, Konrad Kuijken, Henk Hoekstra, Lance Miller, Catherine Heymans, Hendrik Hildebrandt, Jan Luca van den Busch, Angus H. Wright, Mijin Yoon, Maciej Bilicki, Matías Bravo, Claudia del P. Lagos

28 pages, 31 figures, 2 tables, submitted to A&A

We present SKiLLS, a suite of multi-band image simulations for the KiDS-Legacy analysis, the weak lensing analysis of the complete Kilo-Degree Survey. The resulting catalogues enable joint shear and redshift calibration, enhancing the realism and hence accuracy over previous efforts. To create a large volume of simulated galaxies with faithful properties and to a sufficient depth, we integrate cosmological simulations with high-quality imaging observations. We also improve the realism of simulated images by allowing the point spread function (PSF) to differ between CCD images, including stellar density variations and varying noise levels between pointings. Using realistic variable shear fields, we account for the impact of blended systems at different redshifts. Although the overall correction is minor, we find a clear redshift-bias correlation in the blending-only variable shear simulations, indicating the non-trivial impact of this higher-order blending effect. We also explore the impact of the PSF modelling errors and find a small but noticeable effect on the shear bias. Finally, we conduct a series of sensitivity tests, including changing the input galaxy properties. We conclude that our fiducial shape measurement algorithm, lensfit, is robust within the requirements of lensing analyses with KiDS. As for future weak lensing surveys with tighter requirements, we suggest further investments in understanding the impact of blends at different redshifts, improving the PSF modelling algorithm and developing the shape measurement method to be less sensitive to the galaxy properties.

Pierre Salati

Contribution given at the 14th International Conference on Identification of Dark Matter, Vienna, Austria, 18-22 July 2022

Dark matter particles could be the major component of the haloes of galaxies. Their mutual annihilations or decays would produce an indirect signature under the form of high-energy cosmic-rays. The focus of this presentation is on antimatter species, a component so rare that any excess over the background should be easily detected. After a recap on Galactic propagation, I will discuss positrons, antiprotons and anti-nuclei. For each of these species, anomalies have been reported. The antiproton excess, for instance, is currently a hot topic. Alas, it does not resist a correct treatment of theoretical and data errors.

Subhajeet Karmakar, Sachindra Naik, Jeewan C. Pandey, Igor S. Savanov

23 pages, 12 figures, 5 tables, accepted for publication in MNRAS

We present an in-depth study of a large and long duration ($>$1.3 days) X-ray flare observed on an RS CVn type eclipsing binary system SZ Psc using observations from Swift observatory. In the 0.35$-$10 keV energy band, the peak luminosity is estimated to be 4.2$\times$10$^{33}$ erg s$^{-1}$. The quiescent corona of SZ Psc was observed $\sim$5.67 d after the flare using Swift observatory, and also $\sim$1.4 yr after the flare using the XMM-Newton satellite. The quiescent corona is found to consist of three temperature plasma: 4, 13, and 48 MK. High-resolution X-ray spectral analysis of the quiescent corona of SZ Psc suggests that the high first ionization potential (FIP) elements are more abundant than the low-FIP elements. The time-resolved X-ray spectroscopy of the flare shows a significant variation in the flare temperature, emission measure, and abundance. The peak values of temperature, emission measure, and abundances during the flare are estimated to be 199$\pm$11 MK, 2.13$\pm$0.05 $\times 10^{56}$ cm$^{-3}$, 0.66$\pm$0.09 Z$_{\odot}$, respectively. Using the hydrodynamic loop modeling, we derive the loop length of the flare as 6.3$\pm$0.5 $\times 10^{11}$ cm, whereas the loop pressure and density at the flare peak are derived to be 3.5$\pm$0.7 $\times 10^{3}$ dyne cm$^{-2}$ and 8$\pm$2 $\times 10^{10}$ cm$^{-3}$, respectively. The total magnetic field to produce the flare is estimated to be 490$\pm$60 G. The large magnetic field at the coronal height is supposed to be due to the presence of an extended convection zone of the sub-giant and the high orbital velocity.

Ali Baktash, Dieter Horns, Manuel Meyer

The nearby GRB22109A at redshift $z=0.1505$ has been observed up to a maximum energy of 18 TeV with the LHAASO air shower array. The expected optical depth for a photon with energy $E_\gamma=18$~TeV varies between 9.2 and 27.1 according to existing models of the extra-galactic background light in the relevant mid infra-red range. The resulting suppression of the flux makes it unlikely that this photon could be observed. If the photon event and its energy are however confirmed and possibly even more photons above 10 TeV have been observed, the photon-pairproduction process would have to be suppressed by mechanisms predicted in extensions of the standard model of particle physics. We consider the possibilities of photon mixing with a light pseudo-scalar (e.g., axion-like particle: ALP) in the magnetic field of the host galaxy and the Milky Way and Lorentz-invariance violation (LIV). In the case of photon-ALP mixing, the boost factor would reach values $\sim10^6$ for photon couplings not ruled out by the CAST experiment. In the case of LIV, required boost factors are achievable for a LIV breaking energy scale $\lesssim 10^{28}$ eV for the linear modification of the dispersion relation and $\lesssim 10^{20}$ eV for the quadratic modification. A more simple explanation would be a misidentification of a charged cosmic-ray air shower.

Dustin D. Nguyen, Todd A. Thompson, Evan E. Schneider, Sebastian Lopez, Laura A. Lopez

submitted (13 October 2022)

The galactic wind model of Chevalier1985 (CC85) assumes $\textit{uniform}$ energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are inherently non-uniform. We generalize to cases with non-uniform energy/mass injection that scale as $r^{-\Delta}$ within the starburst volume $R$, providing solutions for $\Delta = 0$, 1/2, 1, 3/2, and 2. In marked contrast with the CC85 model ($\Delta=0$), which predicts zero velocity at the center, for a singular isothermal sphere profile ($\Delta=2$), we find that the flow maintains a $\textit{constant}$ Mach number of $\mathcal{M}=\sqrt{3/5} \simeq 0.77$ throughout the volume. The fast interior flow can be as $v_{r < R} = (\dot{E}_T/3\dot{M}_T)^{1/2} \simeq 0.41 \, v_\infty$, where $v_\infty$ is the asymptotic velocity, and $\dot{E}_T$ and $\dot{M}_T$ are the total energy and mass injection rates. For $v_\infty \simeq 2000 \, \mathrm{km \, s^{-1}}$, $v_{r<R} \simeq 820 \, \mathrm{km\, s^{-1}}$ throughout the wind-driving region. The temperature and density profiles of the non-uniform models may be important for interpreting spatially-resolved maps of starburst nuclei. We compute velocity resolved spectra to contrast the CC85 and $\Delta=2$ models. Next generation X-ray space telescopes such as XRISM may assess these kinematic predictions.

Natalie B. Hogg, Pierre Fleury, Julien Larena, Matteo Martinelli

16 + 3 pages, 14 figures, prepared for submission to MNRAS

Line-of-sight effects in strong gravitational lensing have long been treated as a nuisance. However, it was recently proposed that the line-of-sight shear could be a cosmological observable in its own right, if it is not degenerate with lens model parameters. Using the formalism introduced by Fleury et al. (2021a), we firstly demonstrate that the line-of-sight shear can be accurately measured from a simple simulated strong lensing image with percent precision. We then extend our analysis to more complex simulated images and stress test the recovery of the line-of-sight shear when using deficient fitting models, finding that it escapes from degeneracies with lens model parameters, albeit at the expense of the precision. Lastly, we check the validity of the tidal approximation by simulating and fitting an image generated in the presence of many line-of-sight dark matter haloes, finding that an explicit violation of the tidal approximation does not necessarily prevent one from measuring the line-of-sight shear.

Tariq Yasin, Harry Desmond, Julien Devriendt, Adrianne Slyz

13 pages, 7 figures

Rotation curves are the key observational manifestation of the dark matter distribution around late-type galaxies. In a halo model context, the precision of constraints on halo parameters is a complex function of the properties of the measurements as well as properties of the galaxy itself. Forthcoming surveys will resolve rotation curves to varying degrees of precision, or measure their integrated effect in the HI linewidth. To ascertain the relative significance of the relevant quantities for constraining halo properties, we study the information on halo mass and concentration as quantified by the Kullback-Leibler divergence of the kinematics-informed posterior from the uninformative prior. We calculate this divergence as a function of the different types of spectroscopic observation, properties of the measurement, galaxy properties and auxiliary observational data on the baryonic components. Using the SPARC sample, we find that fits to the full rotation curve exhibit a large variation in information gain between galaxies, ranging from 1 to 11 bits. This range is equivalent to the difference between the prior shrinking to a flat posterior that is only a factor of 2 smaller compared to a factor of 2000 smaller. It is predominantly caused by the vast differences in the number of data points and the size of velocity uncertainties between the SPARC galaxies. Out of the galaxy properties, only the degree of dark matter dominance is important. We also study the relative importance of the minimum HI surface density probed and the size of velocity uncertainties on the constraining power on the inner halo slope, finding the latter to be significantly more important. We spell out the implications of these results for optimising the strategy of galaxy surveys aiming to constrain galaxies' dark matter distributions, highlighting spectroscopic precision as the most important factor.

We derive the fixed-$\Lambda$ and unimodular propagators using the path integral formalism as applied to the Einstein-Cartan action. The simplicity of the action (which is linear in the lapse function) allows for an exact integration starting from the lapse function and the enforcement of the Hamiltonian constraint, leading to a product of Chern-Simons states if the connection is fixed at the endpoints. No saddle point approximation is needed. Should the metric be fixed at the endpoints, then, depending on the contour chosen for the connection, Hartle-Hawking or Vilenkin propagators are obtained. Thus, in this approach one trades a choice of contour in the lapse function for one in the connection, where appropriate. The unimodular propagators are also trivial to obtain via the path integral, and the previously derived expressions are recovered.

Ali Zafari, Atefeh Khoshkhahtinat, Piyush M. Mehta, Nasser M. Nasrabadi, Barbara J. Thompson, Daniel da Silva, Michael S. F. Kirk

Accepted to ICMLA 2022 (Oral Presentation)

NASA's Solar Dynamics Observatory (SDO) mission gathers 1.4 terabytes of data each day from its geosynchronous orbit in space. SDO data includes images of the Sun captured at different wavelengths, with the primary scientific goal of understanding the dynamic processes governing the Sun. Recently, end-to-end optimized artificial neural networks (ANN) have shown great potential in performing image compression. ANN-based compression schemes have outperformed conventional hand-engineered algorithms for lossy and lossless image compression. We have designed an ad-hoc ANN-based image compression scheme to reduce the amount of data needed to be stored and retrieved on space missions studying solar dynamics. In this work, we propose an attention module to make use of both local and non-local attention mechanisms in an adversarially trained neural image compression network. We have also demonstrated the superior perceptual quality of this neural image compressor. Our proposed algorithm for compressing images downloaded from the SDO spacecraft performs better in rate-distortion trade-off than the popular currently-in-use image compression codecs such as JPEG and JPEG2000. In addition we have shown that the proposed method outperforms state-of-the art lossy transform coding compression codec, i.e., BPG.

Jeff A. Dror, Stefania Gori, Jacob M. Leedom, Nicholas L. Rodd

15 pages, 3 figures

A leading direction in the hunt for axion dark matter is to search for its influence on nuclear spins. The detection scheme involves polarizing a sample of nuclei within a strong static magnetic field and then looking for a spin precession induced by the oscillating axion field. We study the axion signal and background contributions that arise in such experiments (a prominent example being CASPEr), finding key differences with the existing literature. Most importantly, in the limit where the transverse spin-relaxation time of the material is the largest timescale of the problem, we show that the induced signal continues to grow even beyond the coherence time of the axion field. As a result, we find that spin-precession instruments are much more sensitive than what has been previously estimated in a sizable range of axion masses, with sensitivity improvement of up to a factor of 100 at an axion mass of 100 neV using a Xenon-129 sample. This improves the detection prospects for the QCD axion, and we estimate the experimental requirements to reach this motivated target. Our results apply to both the axion electric and magnetic dipole moment operators.

Gordan Krnjaic, Duncan Rocha, Anastasia Sokolenko

arXiv admin note: text overlap with arXiv:1309.4091 by other authors

We study a simple model of vector dark matter that couples to Standard Model particles via magnetic dipole interactions. In this scenario, the cosmological abundance arises through the freeze-in mechanism and depends on the dipole coupling, the vector mass, and the reheat temperature. To ensure cosmological metastability, the vector must be lighter than the fermions to which it couples, but rare decays can still produce observable 3$\gamma$ final states; two-body decays can also occur at one-loop with additional weak suppression, but are subdominant if the vector couples mainly to light fermions. For sufficiently heavy vectors, induced kinetic mixing with the photon can also yield additional two body decays to lighter fermions and predict indirect detection signals through final state radiation. We explore the implications of couplings to various flavors of visible particles and emphasize leptophilic dipoles involving electrons, muons, and taus, which offer the most promising indirect detection signatures through 3$\gamma$, $e^+ e^- \gamma$, and $\mu^+ \mu^- \gamma$ decay channels. We also present constraints from current and past telescopes, and sensitivity projections for future missions including e-ASTROGAM and AMEGO.

Kim V. Berghaus, Angelo Esposito, Rouven Essig, Mukul Sholapurkar

30 pages, 8 figures

Dark matter scattering off a nucleus has a small probability of inducing an observable ionization through the inelastic excitation of an electron, called the Migdal effect. We use an effective field theory to extend the computation of the Migdal effect in semiconductors to regions of small momentum transfer to the nucleus, where the final state of the nucleus is no longer well described by a plane wave. Our analytical result can be fully quantified by the measurable dynamic structure factor of the semiconductor, which accounts for the vibrational degrees of freedom (phonons) in a crystal. We show that, due to the sum rules obeyed by the structure factor, the inclusive Migdal rate and the shape of the electron recoil spectrum is well captured by approximating the nuclei in the crystal as free ions; however, the exclusive differential rate with respect to energy depositions to the crystal depends on the phonon dynamics encoded in the dynamic structure function of the specific material. Our results now allow the Migdal effect in semiconductors to be evaluated even for the lightest dark matter candidates ($m_\chi \gtrsim 1$ MeV) that can kinematically excite electrons.

We propose a method to determine the mass and spin of primordial black holes (PBHs) in the mass range $5\times 10^7-10^{12}$ kg (Hawking temperatures $\sim10$ MeV $-200$ GeV), based on measuring the energy of specific features in the photon Hawking emission spectrum, including both primary and secondary components. This is motivated by scenarios where PBHs in this mass range spin up as they evaporate, namely the string axiverse, where dimensionless spin parameters $\tilde{a}\sim 0.1-0.5$ can be achieved through the Hawking emission of hundreds or even thousands of light axion-like particles. Measuring the present PBH mass-spin distribution may thus be an important probe of physics beyond the Standard Model. Since the proposed method relies on the energy of the photons emitted by a given PBH, rather than on the associated flux, it is independent of the PBH-Earth distance and, as a byproduct, can also be used to infer the latter.

We investigate the Fermi acceleration of charged particles in 2D MHD anti-parallel plasmoid reconnection, finding a drastic enhancement in energization rate $\dot{\varepsilon}$ over a standard Fermi model of $\dot{\varepsilon} \sim \varepsilon$. The shrinking particle orbit width around a magnetic island due to $\vec{E}\times\vec{B}$ drift produces a $\dot{\varepsilon}_\parallel \sim \varepsilon_\parallel^{1+1/2\chi}$ power law with $\chi \sim 0.75$. The increase in the maximum possible energy gain of a particle within a plasmoid due to the enhanced efficiency increases with the plasmoid size, and is by multiple factors of 10 in the case of solar flares and much more for larger plasmas. Including effects of the non-constant $\vec{E}\times\vec{B}$ drift rates leads to further variation of power law indices from $\gtrsim 2$ to $\lesssim 1$, decreasing with plasmoid size at the time of injection.

Giovanni Otalora, Emmanuel N. Saridakis

7 pages, 4 figures

We investigate cosmological scenarios with spin-gravity coupling. In particular, due to the spin of the baryonic and dark matter particles and its coupling to gravity, they probe an effective spin-dependent metric, which can be calculated semi-classically in the Mathisson-Papapetrou-Tulczyjew-Dixon formalism. Hence, the usual field equations give rise to modified Friedmann equations, in which the extra terms can be identified as an effective dark-energy sector. Additionally, we obtain an effective interaction between the matter and dark-energy sectors. In the case where the spin-gravity coupling switches off, we recover standard $\Lambda$CDM cosmology. We perform a dynamical system analysis and we find a matter-dominated point that can describe the matter era, and a stable late-time solution corresponding to acceleration and dark-energy domination. For small values of the spin coupling parameter, deviations from $\Lambda$CDM concordance scenario are small, however for larger values they can be brought to the desired amount, leading to different dark-energy equation-of-state parameter behavior, as well as to different transition redshift from acceleration to deceleration. Finally, we confront the model predictions with Hubble function data.

Numerical simulation of atmospheric turbulence is one of the biggest bottlenecks in developing computational techniques for solving the inverse problem in long-range imaging. The classical split-step method is based upon numerical wave propagation which splits the propagation path into many segments and propagates every pixel in each segment individually via the Fresnel integral. This repeated evaluation becomes increasingly time-consuming for larger images. As a result, the split-step simulation is often done only on a sparse grid of points followed by an interpolation to the other pixels. Even so, the computation is expensive for real-time applications. In this paper, we present a new simulation method that enables \emph{real-time} processing over a \emph{dense} grid of points. Building upon the recently developed multi-aperture model and the phase-to-space transform, we overcome the memory bottleneck in drawing random samples from the Zernike correlation tensor. We show that the cross-correlation of the Zernike modes has an insignificant contribution to the statistics of the random samples. By approximating these cross-correlation blocks in the Zernike tensor, we restore the homogeneity of the tensor which then enables Fourier-based random sampling. On a $512\times512$ image, the new simulator achieves 0.025 seconds per frame over a dense field. On a $3840 \times 2160$ image which would have taken 13 hours to simulate using the split-step method, the new simulator can run at approximately 60 seconds per frame.

Heliudson Bernardo, Benjamin Bose, Guilherme Franzmann, Steffen Hagstotz, Yutong He, Aliki Litsa, Florian Niedermann

54 pages, 2 figures, 1 table. Invited review for Universe, Special Issue "Cosmological Constant" ( this https URL )

The cosmological constant and its phenomenology remain among the greatest puzzles in theoretical physics. We review how modifications of Einstein's general relativity could alleviate the different problems associated with it that result from the interplay of classical gravity and quantum field theory. We introduce a modern and concise language to describe the problems associated with its phenomenology, and inspect no-go theorems and their loopholes to motivate the approaches discussed here. Constrained gravity approaches exploit minimal departures from general relativity; massive gravity introduces mass to the graviton; Horndeski theories lead to the breaking of translational invariance of the vacuum; and models with extra dimensions change the symmetries of the vacuum. We also review screening mechanisms that have to be present in some of these theories if they aim to recover the success of general relativity on small scales as well. Finally, we summarise the statuses of these models in their attempt to solve the different cosmological constant problems while being able to account for current astrophysical and cosmological observations.

Tanmay Kumar Poddar

6 pages, 3 figures, 3 tables, proceedings for 41st International Conference on High Energy Physics - ICHEP2022

The Standard Model of particle physics is a $SU(3)_c\times SU(2)_L\times U(1)_Y$ gauge theory that can explain the strong, weak, and electromagnetic interactions between the particles. The gravitational interaction is described by Einstein's General Relativity theory which is a classical theory of gravity. These theories can explain all the four fundamental forces of nature with great level of accuracy. However, there are several theoretical and experimental motivations of studying physics beyond the Standard Model of particle physics and Einstein's General Relativity theory. Probing these new physics scenarios with ultralight particles has its own importance as they can be a promising candidates for dark matter that can evade the constraints from dark matter direct detection experiments and solve the small scale structure problems of the universe. In this paper, we have considered axions and gauge bosons as light particles and their possible searches through astrophysical observations. In particular, we obtain constraints on ultralight axions from orbital period loss of compact binary systems, gravitational light bending, and Shapiro time delay. We also derive constraints on ultralight gauge bosons from indirect evidence of gravitational waves, and perihelion precession of planets. Such type of observations can also constrain several particle physics models and are discussed.

Mohsen Khodadi, Reza Pourkhodabakhshi

13 pages (two column), 6 figures, 1 table

Does circumventing the curvature singularity of the Kerr black hole affects the timescale of the scalar cloud formation around it? By definition, the scalar cloud, forms a gravitational atom with hydrogen-like bound states, lying on the threshold of a massive scalar field's superradiant instability regime (time-growing quasi-bound states) and beyond (time-decaying quasi-bound states). By taking a novel type of rotating hollow regular black hole proposed by Simpson and Visser which unlike its standard rivals has an asymptotically Minkowski core, we address this question. The metric has a minimal extension relative to the standard Kerr, originating from a single regularization parameter $\ell$, with length dimension. We show with the inclusion of the regularization length scale $\ell$ into the Kerr spacetime, without affecting the standard superradiant instability regime, the timescale of scalar cloud formation gets shorter. Since the scalar cloud after its formation, via energy dissipation, can play the role of a continuum source for gravitational waves, such a reduction in the instability growth time improves the phenomenological detection prospects of new physics because the shorter the time, the more astrophysically important.

L. Gottardi (1), M. de Wit (1), E. Taralli (1), K. Nagayashi (1), A. Kozorezov (2) ((1) NWO-I/SRON Netherlands Institute for Space Research, Leiden, Netherlands (2) Department of Physics, Lancaster University, United Kingdom)

9 pages, 7 figures

We present a detailed analysis of the fundamental noise sources in superconducting transition-edge sensors (TESs), ac voltage biased at MHz frequencies and treated as superconducting weak links. We have studied the noise in the resistive transition as a function of bath temperature of several detectors with different normal resistances and geometries. We show that the excess noise, typically observed in the TES electrical bandwidth, can be explained by the equilibrium Johnson noise of the quasiparticles generated within the weak link. The fluctuations at the Josephson frequency and higher harmonics contribute significantly to the measured voltage noise at the detector bandwidth through the nonlinear response of the weak link with a sinusoidal current-phase relation.