Abstract visibility
Abstract text size

Papers for Monday, Dec 12 2022

Papers with local authors

Yaoxin Chen, Yingzhong Xu, Xi Kang

Accepted for publication in MNRAS. 14 pages (including references), 11 figures

0 votes
Paper 40 — arXiv:2212.05007
0 votes
Paper 40 — arXiv:2212.05007

Galaxies display several well-behaved scaling relations between their properties, such as the star formation rate-stellar mass relation (the main sequence) and the stellar mass-halo mass relation (SHMR). In principle, these scaling relations could imply different star formation histories (SFHs) of galaxies and different constraints on galaxy formation physics. In this paper, we derive the SFHs of galaxies by assuming that they always follow the SHMRs at different redshifts and use an empirical model to constrain key processes in their baryon cycle. It is found that, besides cold accretion due to halo growth, outflow of gas produced by stellar feedback has to be recycled to sustain the derived SFHs of galaxies. The recycled fraction is strongly affected by the baryon fraction in accreted low-mass haloes and the mass loading factor which quantifies the ratio between the galactic outflow rate and star formation rate. Our fiducial model predicts that around $20-60\%$ of outflow is recycled in $\sim0.5-4Gyrs$, while simulations predict a slightly higher recycle fraction and a lower recycle time. We argue that strong constraints on the baryon cycle process can be obtained from future observation of the circum-galactic medium (CGM) of galaxies, such as the gas cooling rate of CGM. We also find that the implied SFHs from the SHMRs indicate that galaxies stay on the main sequences only for part of their lifetimes. Our model reproduces the evolution of the mass-metallicity relation as well.

All other papers

Shinya Wanajo, Sho Fujibayashi, Kota Hayashi, Kenta Kiuchi, Yuichiro Sekiguchi, Masaru Shibata

5 pages, 4 figures

We examine nucleosynthesis in the ejecta of black hole-neutron star mergers based on the results of self-consistent, long-term neutrino-radiation-magnetohydrodynamics simulations for the first time. We find that the combination of dynamical and post-merger ejecta reproduces a solar-like r-process pattern. Moreover, the enhancement level of actinides is highly sensitive to the distribution of both electron fraction and the velocity of the dynamical ejecta. Our result implies that the mean electron fraction of dynamical ejecta should be ~ 0.05-0.08 in order to reconcile the nucleosynthetic abundances with those in r-process-enhanced, actinide-boosted stars. This result provides an important constraint for nuclear equations of state.

Jacob B. Simon, Jürgen Blum, Til Birnstiel, David Nesvorný

Planetesimal Formation Review accepted for publication in Comets III

In this chapter, we review the processes involved in the formation of planetesimals and comets. We will start with a description of the physics of dust grain growth and how this is mediated by gas-dust interactions in planet-forming disks. We will then delve into the various models of planetesimal formation, describing how these planetesimals form as well as their resulting structure. In doing so, we focus on and compare two paradigms for planetesimal formation: the gravitational collapse of particle over-densities (which can be produced by a variety of mechanisms) and the growth of particles into planetesimals via collisional and gravitational coagulation. Finally, we compare the predictions from these models with data collected by the Rosetta and New Horizons missions and that obtained via observations of distant Kuiper Belt Objects.

Sophie E. van Mierlo, Karina I. Caputi, Vasily Kokorev

Submitted to ApJL

Recent controversy regarding the existence of massive ($\log(M_*/M_\odot) \gtrsim 11$) galaxies at $z>6$ is posing a challenge for galaxy formation theories. Hence, it is of critical importance to understand the effects of SED fitting methods on stellar mass estimates of Epoch of Re-ionisation galaxies. In this work, we perform a case study on the AGN-host galaxy candidate COS-87259 with spectroscopic redshift $z_{\rm spec}=6.853$, that is claimed to have an extremely high stellar mass of $\log(M_*/M_\odot) \sim 11.2$. We test a suite of different SED fitting algorithms and stellar population models on our independently measured photometry in 17 broad bands for this source. Between five different code set-ups, the stellar mass estimates for COS-87259 span $\log(M_*/M_\odot) = 10.24$--11.00, whilst the reduced $\chi^2$ values of the fits are all close to unity within $\Delta\chi^2_\nu=0.9$, so that the quality of the SED fits is basically indistinguishable. Only the Bayesian inference code Prospector using a non-parametric star formation history model yields a stellar mass exceeding $\log(M_*/M_\odot)=11$. As this SED fitting prescription is becoming increasingly popular for James Webb Space Telescope high-redshift science, we stress the absolute importance to test various SED fitting routines particularly on apparently very massive galaxies at such high redshifts. Ultimately, we conclude that the extremely high stellar mass estimate for COS-87259 is not necessary, deriving equally good fits with stellar masses $\sim 1$ dex lower.

Nayantara Mudur, Core Francisco Park, Douglas P Finkbeiner

21 pages, 10 figures

Cosmological surveys must correct their observations for the reddening of extragalactic objects by Galactic dust. Existing dust maps, however, have been found to have spatial correlations with the large-scale structure of the Universe. Errors in extinction maps can propagate systematic biases into samples of dereddened extragalactic objects and into cosmological measurements such as correlation functions between foreground lenses and background objects and the primordial non-gaussianity parameter $f_{NL}$. Emission-based maps are contaminated by the cosmic infrared background, while maps inferred from stellar-reddenings suffer from imperfect removal of quasars and galaxies from stellar catalogs. Thus, stellar-reddening based maps using catalogs without extragalactic objects offer a promising path to making dust maps with minimal correlations with large-scale structure. We present two high-latitude integrated extinction maps based on stellar reddenings, with a point spread function of full-width half-maximum 6.1' and 15'. We employ a strict selection of catalog objects to filter out galaxies and quasars and measure the spatial correlation of our extinction maps with extragalactic structure. Our galactic extinction maps have reduced spatial correlation with large scale structure relative to most existing stellar-reddening based and emission-based extinction maps.

Yuxi (Lucy)Lu, Ivan Minchev, Tobias Buck, Sergey Khoperskov, Matthias Steinmetz, Noam Libeskind, Gabriele Cescutti, Ken C. Freeman

Under review at Nature Letters, submitted Oct. 19, 2022

Stars move away from their birth places over time via a process known as radial migration, which blurs chemo-kinematic relations used for reconstructing the Milky Way formation history. One of the ultimate goals of Galactic Archaeology, therefore, is to find stars' birth aggregates in the disk via chemical tagging. Here we show that stellar birth radii can be derived directly from the data with minimum prior assumptions on the Galactic enrichment history. We recover the time evolution of the stellar birth metallicity gradient, $d$[Fe/H]($R$, $\tau$)/$dR$, through its inverse relation to the metallicity range as a function of age today, allowing us to place any star with age and metallicity measurements back to its birthplace, $R_b$. Applying our method to a high-precision large data set of Milky Way disk subgiant stars, we find a steepening of the birth metallicity gradient from 11 to 8 Gyr ago, which coincides with the time of the last major merger, Gaia-Sausage-Enceladus (GSE). This transition appears to play a major role in shaping both the age-metallicity relation and the bimodality in the [$\alpha$/Fe]-[Fe/H] plane. By dissecting the disk into mono-$R_b$ populations, clumps in the low-[$\alpha$/Fe] sequence appear, which are not seen in the total sample and coincide in time with known star-formation bursts. We estimated that the Sun was born at $4.5 \pm 0.4$ kpc from the Galactic center. Our $R_b$ estimates provide the missing piece needed to recover the Milky Way formation history, while the by-product,[Fe/H]$(R$, $\tau)$, can be used as the thus-far missing prior for chemical evolution modeling.

Bogumił Pilecki, Ian B. Thompson, Felipe Espinoza-Arancibia, Richard I. Anderson, Wolfgang Gieren, Weronika Narloch, Javier Minniti, Grzegorz Pietrzyński, Mónica Taormina, Giuseppe Bono, Gergely Hajdu

6 pages, 4 figures, 1 table, published in The Astrophysical Journal Letters

We report the discovery of a surprising binary configuration of the double-mode Cepheid OGLE-LMC-CEP-1347 pulsating in the first (P_1=0.690d) and second overtone (P_2=0.556d) modes. The orbital period (P_orb=59d) of the system is five times shorter than the shortest known to date (310d) for a binary Cepheid. The Cepheid itself is also the shortest-period one ever found in a binary system and the first double-mode Cepheid in a spectroscopically double-lined binary. OGLE-LMC-CEP-1347 is most probably on its first crossing through the instability strip, as inferred from both its short period and fast period increase, consistent with evolutionary models, and from the short orbital period (not expected for binary Cepheids whose components have passed through the red giant phase). Our evolutionary analysis yielded a first-crossing Cepheid with a mass in a range of 2.9-3.4 Msun (lower than any measured Cepheid mass), consistent with observations. The companion is a stable star, at least two times fainter and less massive than the Cepheid (preliminary mass ratio q=0.55), while also redder and thus at the subgiant or more advanced evolutionary stage. To match these characteristics, the Cepheid has to be a product of binary interaction, most likely a merger of two less massive stars, which makes it the second known classical Cepheid of binary origin. Moreover, further evolution of the components may lead to another binary interaction.

Kiyam Lin, Maximilian von Wietersheim-Kramsta, Benjamin Joachimi, Stephen Feeney

14 pages, 8 figures; submitted to MNRAS

The standard approach to inference from cosmic large-scale structure data employs summary statistics that are compared to analytic models in a Gaussian likelihood with pre-computed covariance. To overcome the idealising assumptions about the form of the likelihood and the complexity of the data inherent to the standard approach, we investigate simulation-based inference (SBI), which learns the likelihood as a probability density parameterised by a neural network. We construct suites of simulated, exactly Gaussian-distributed data vectors for the most recent Kilo-Degree Survey (KiDS) weak gravitational lensing analysis and demonstrate that SBI recovers the full 12-dimensional KiDS posterior distribution with just under $10^4$ simulations. We optimise the simulation strategy by initially covering the parameter space by a hypercube, followed by batches of actively learnt additional points. The data compression in our SBI implementation is robust to suboptimal choices of fiducial parameter values and of data covariance. Together with a fast simulator, SBI is therefore a competitive and more versatile alternative to standard inference.

Samuel Brieden, Héctor Gil-Marín, Licia Verde

42 pages, 12 figures, 3 tables, to be submitted to JCAP. Comments welcome

We use the large-scale structure galaxy data (LSS) from the BOSS and eBOSS surveys, in combination with abundances information from Big Bang Nucleosynthesis (BBN) to measure two values of the Hubble expansion rate, $H_0=100h\,[{\rm km}\, {\rm s}^{-1}\,{\rm Mpc}^{-1}]$, each of them based on very different physical processes. One is a (traditional) late-time-background measurement based on determining the BAO scale and using BBN abundances on baryons for calibrating its absolute size (BAO+BBN). This method anchors $H_0$ to the (standard) physics of the sound horizon scale at pre-recombination times. The other is a newer early-time based measurement associated with the broadband shape of the power spectrum. This second method anchors $H_0$ to the physics of the matter-radiation equality scale, which also needs BBN information for determining the suppression of baryons in the power spectrum shape (shape+BBN). Within the $\Lambda$CDM model, we find very good consistency among these two $H_0$'s: BAO+BBN (+growth) delivers $H_0=67.42_{-0.94}^{+0.88}$ $(67.37_{-0.95}^{+0.86})$ km s$^{-1}$Mpc$^{-1}$ , whereas the shape+BBN (+growth) delivers $H_0 = 70.1_{-2.1}^{+2.1}$ $(70.1_{-2.1}^{+1.9})$ km s$^{-1}$ Mpc$^{-1}$, where "growth" stands for information from the late-time-perturbations captured by the growth of structure parameter. These are the tightest sound-horizon free $H_0$ constraints from LSS data to date. As a consequence to be viable, any $\Lambda$CDM extension proposed to address the so-called "Hubble tension" needs to modify consistently not only the sound horizon scale physics, but also the matter-radiation equality scale, in such a way that both late- and early-based $H_0$'s return results mutually consistent and consistent with the high $H_0$ value recovered by the standard cosmic distance ladder (distance-redshift relation) determinations.

Alan E. E. Rogers, John P. Barrett, Judd D. Bowman, Rigel Cappallo, Colin J. Lonsdale, Nivedita Mahesh, Raul A. Monsalve, Steven G. Murray, Peter H. Sims

Scattering from objects near an antenna produce correlated signals from strong compact radio sources in a manner similar to those used by the Sea Interferometer to measure the radio source positions using the fine frequency structure in the total power spectrum of a single antenna. These fringes or ripples due to correlated signal interference are present at a low level in the spectrum of any single antenna and are a major source of systematics in systems used to measure the global redshifted 21-cm signal from the early universe. In the Sea Interferometer a single antenna on a cliff above the sea is used to add the signal from the direct path to the signal from the path reflected from the sea thereby forming an interferometer. This was used for mapping radio sources with a single antenna by Bolton and Slee in the 1950s. In this paper we derive analytic expressions to determine the level of these ripples and compare these results in a few simple cases with electromagnetic modeling software to verify that the analytic calculations are sufficient to obtain the magnitude of the scattering effects on the measurements of the global 21-cm signal. These analytic calculations are needed to evaluate the magnitude of the effects in cases that are either too complex or take too much time to be modeled using software.

Erick Pastén, Víctor H. Cárdenas

We revisited the problem of describing, on average, a fractal distribution of matter using a Lemaitre-Tolman-Bondi (LTB) solution. Here we study the fractal structure of our local universe having a fractal dimension and a scale transition. We test our model with the latest type Ia supernova data, the Pantheon compilation, and discuss problems and possible improvements for it, concluding that a fractal transition in LTB cosmology cannot be used to explain the effects of dark energy but it can be useful to study structures at low scales.

Kielan K. W. Hoch, Quinn M. Konopacky, Christopher A. Theissen, Jean-Baptiste Ruffio, Travis S. Barman, Marshall D. Perrin, Bruce Macintosh, Christian Marois

20 pages, 12 figures

The carbon-to-oxygen (C/O) ratio in an exoplanet atmosphere has been suggested as a potential diagnostic of planet formation. Now that a number of exoplanets have measured C/O ratios, it is possible to examine this diagnostic at a population level. Here, we present an analysis of currently measured C/O ratios of directly imaged and transit/eclipse planets. First, we derive atmospheric parameters for the substellar companion HD 284149 b using data that were taken with the OSIRIS integral field spectrograph at the W.M. Keck Observatory and report two non-detections from our ongoing imaging spectroscopy survey of exoplanetary atmospheres with Keck/OSIRIS. We find an effective temperature of $T_\mathrm{eff} = 2502$~K, with a range of 2291--2624~K, $\log g=4.52$, with a range of 4.38--4.91, and [M/H] = 0.37, with a range of 0.10--0.55. These values are in agreement with previous studies done by Bonavita et al. (2014, 2017). We derive a C/O of 0.589$^{+0.148}_{-0.295}$ for HD 284149 b. We then add this measurement to the growing list of C/O ratios for directly imaged planets from the literature, and compare them with those available from a sample of transit/eclipse planets. There is a trend in C/O ratio with companion mass (M$_{\mathrm{Jup}}$), with a break seen around 4 M$_{\mathrm{Jup}}$. We run a Kolmogorov-Smirnov and an Anderson-Darling test on planets above and below this mass boundary, and find that they are two distinct populations. This could be additional evidence of two distinct populations possibly having two different formation pathways, with companion mass as a primary indicator of most likely formation scenario.

We model Pop III star formation in different FUV and X-ray backgrounds, including radiation feedback from protostars. We confirm previous results that a moderate X-ray background increases the number of Pop III systems per unit cosmological volume, but masses and multiplicities of the system are reduced. The stellar mass function also agrees with previous results, and we confirm the outward migration of the stars within the protostellar discs. We find that nearly all Pop III star systems are hierarchical, i.e., binaries of binaries. Typically, two equal-mass stars form near the centre of the protostellar disc and migrate outward. Around these stars, mini-discs fragment forming binaries that also migrate outward. Stars may also form at Lagrange points L4/L5 of the system. Afterward, star formation becomes more stochastic due to the large multiplicity, and zero-metallicity low-mass stars can form when rapidly ejected from the disc. Stars in the disc often have eccentric orbits, leading to a periodic modulation of their accretion rates and luminosities. At the pericenter, due to strong accretion, the star can enter a red-supergiant phase reaching nearly Eddington luminosity in the optical bands ($m_{\rm AB} \sim 34$ for a $100~M_{odot}$ star at $z=6$). During this phase, the star, rather than its nebular lines, can be observed directly by JWST, if sufficiently magnified by a gravitational lens. The $\sim 10,000$ AU separations and high eccentricities of many Pop III star binaries in our simulations are favorable parameters for IMBH mergers - and gravitational waves emission - through orbital excitation by field stars.

Emma Curtis-Lake, Stefano Carniani, Alex Cameron, Stephane Charlot, Peter Jakobsen, Roberto Maiolino, Andrew Bunker, Joris Witstok, Renske Smit, Jacopo Chevallard, Chris Willott, Pierre Ferruit, Santiago Arribas, Nina Bonaventura, Mirko Curti, Francesco D'Eugenio, Marijn Franx, Giovanna Giardino, Tobias J. Looser, Nora Lützgendorf, Michael V. Maseda, Tim Rawle, Hans-Walter Rix, Bruno Rodriguez del Pino, Hannah Übler, Marko Sirianni, Alan Dressler, Eiichi Egami, Daniel J. Eisenstein, Ryan Endsley, Kevin Hainline, Ryan Hausen, Benjamin D. Johnson, Marcia Rieke, Brant Robertson, Irene Shivaei, Daniel P. Stark, Sandro Tacchella, Christina C. Williams, Christopher N. A. Willmer, Rachana Bhatawdekar, Rebecca Bowler, Kristan Boyett, Zuyi Chen, Anna de Graaff, Jakob M. Helton, Raphael E. Hviding, et al. (12 additional authors not shown)

32 pages, 9 figures, Submitted

Finding and characterising the first galaxies that illuminated the early Universe at cosmic dawn is pivotal to understand the physical conditions and the processes that led to the formation of the first stars. In the first few months of operations, imaging from the James Webb Space Telescope (JWST) have been used to identify tens of candidates of galaxies at redshift (z) greater than 10, less than 450 million years after the Big Bang. However, none of these candidates has yet been confirmed spectroscopically, leaving open the possibility that they are actually low-redshift interlopers. Here we present spectroscopic confirmation and analysis of four galaxies unambiguously detected at redshift 10.3<z<13.2, previously selected from NIRCam imaging. The spectra reveal that these primeval galaxies are extremely metal poor, have masses between 10^7 and a few times 10^8 solar masses, and young ages. The damping wings that shape the continuum close to the Lyman edge are consistent with a fully neutral intergalactic medium at this epoch. These findings demonstrate the rapid emergence of the first generations of galaxies at cosmic dawn.

Hyunjin Shim, Ho Seong Hwang, Woong-Seob Jeong, Yoshiki Toba, Minjin Kim, Dohyeong Kim, Hyunmi Song, Tetsuya Hashimoto, Takao Nakagawa, Ambra Nanni, William J. Pearson, Toshinobu Takagi

19 pages, 9 figures, AJ, in press

We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ($z<0.4$) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as $q_\mathrm{PAH}$, which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous mid-infrared photometric data points covering 7-24$\mu$m and far-infrared flux densities, $q_\mathrm{PAH}$ values are derived through spectral energy distribution fitting. The correlation between $q_\mathrm{PAH}$ and other physical properties of galaxies, i.e., gas-phase metallicity ($12+\mathrm{log(O/H)}$), stellar mass, and specific star-formation rate (sSFR) are explored. As in previous studies, $q_\mathrm{PAH}$ values of galaxies with high metallicity are found to be higher than those with low metallicity. The strength of PAH emission is also positively correlated with the stellar mass and negatively correlated with the sSFR. The correlation between $q_\mathrm{PAH}$ and each parameter still exists even after the other two parameters are fixed. In addition to the PAH strength, the application of metallicity-dependent gas-to-dust mass ratio appears to work well to estimate gas mass that matches the observed relationship between molecular gas and physical parameters. The result obtained will be used to calibrate the observed PAH luminosity-total infrared luminosity relation, based on the variation of MIR-FIR SED, which is used in the estimation of hidden star formation.

Kun Wang, Anbing Ren, Mads Fredslund Andersen, Frank Grundahl, Tao Chen, Pere L. Palle

15+4 pages, 8 figures, Submitted to AJ,

We present a detailed analysis of an eclipsing double-lined binary FX UMa based on TESS photometry and newly acquired spectroscopic observations. The radial velocities and atmospheric parameters for each component star are obtained from the SONG high-resolution spectra. Combined with the radial-velocity measurements, our light-curve modeling yields absolute masses and radii of the two components. The Fourier amplitude spectrum of the residual light curve reveals a total of 103 frequencies with signal-to-noise ratio (S/N) > 4, including 12 independent frequencies, 17 multiples of the orbital frequency (Nforb), and 74 combination frequencies. Ten Nforb peaks with S/N > 10 have very high amplitudes and are likely due to tidally excited oscillations (TEOs). The remaining Nforb peaks (4 < S/N < 10) may be originated from the imperfect removal, or they are actually real TEOs. Four anharmonic frequencies can pair up and sum to give exact harmonics of the orbital frequency, suggesting the existence of non-linear tidal processes in the eccentric binary system FX UMa. Eight independent frequencies in the range of 20 to 32 day$^{-1}$ are typical low-order pressure modes of delta Scuti pulsators.

Matthew T. Miles, Ryan M. Shannon, Matthew Bailes, Daniel J. Reardon, Michael J. Keith, Andrew D. Cameron, Aditya Parthasarathy, Mohsen Shamohammadi, Renee Spiewak, Willem van Straten, Sarah Buchner, Fernando Camilo, Marisa Geyer, Aris Karastergiou, Michael Kramer, Maciej Serylak, Gilles Theureau, Vivek Venkatraman Krishnan

18 pages, 6 figures

We present the first 2.5 years of data from the MeerKAT Pulsar Timing Array (MPTA), part of MeerTime, a MeerKAT Large Survey Project. The MPTA aims to precisely measure pulse arrival times from an ensemble of 88 pulsars visible from the Southern Hemisphere, with the goal of contributing to the search, detection and study of nanohertz-frequency gravitational waves as part of the International Pulsar Timing Array. This project makes use of the MeerKAT telescope, and operates with a typical observing cadence of two weeks using the L-band receiver that records data from 856-1712 MHz. We provide a comprehensive description of the observing system, software, and pipelines used and developed for the MeerTime project. The data products made available as part of this data release are from the 78 pulsars that had at least $30$ observations between the start of the MeerTime programme in February 2019 and October 2021. These include both sub-banded and band-averaged arrival times, as well as the initial timing ephemerides, noise models, and the frequency-dependent standard templates (portraits) used to derive pulse arrival times. After accounting for detected noise processes in the data, the frequency-averaged residuals of $67$ of the pulsars achieved a root-mean-square residual precision of $< 1 \mu \rm{s}$. We also present a novel recovery of the clock correction waveform solely from pulsar timing residuals, and an exploration into preliminary findings of interest to the international pulsar timing community. The arrival times, standards and full Stokes parameter calibrated pulsar timing archives are publicly available.

Ye Feng, James F. Steiner, Santiago Ubach Ramirez, Lijun Gou

13 pages, 10 figures, submitted to MNRAS

MAXI J1305-704 is a transient X-ray binary with a black hole primary. It was discovered on April 9, 2012, during its only known outburst. MAXI J1305-704 is also a high inclination low-mass X-ray binary with prominent dip features in its light curves, so we check the full catalog of 92 \emph{Swift}/XRT continuous observations of MAXI J1305-704, focusing only on the stable spectra. We select 13 ``gold" spectra for which the root mean square RMS <0.075 and the coronal scattered fraction $f_{\mathrm{sc}} \lesssim 25 \%$. These ``gold" data are optimal thermal-state observations for continuum-fitting modeling, in which the disk extends to the innermost-stable circular orbit and is geometrically thin. The black hole spin was unknown for this object before. By utilizing the X-ray continuum fitting method with the relativistic thin disk model \texttt{kerrbb2} and supplying the known dynamical binary system parameters, we find MAXI J1305-704 has a moderate spin ($a_{*}=0.87_{-0.13}^{+0.07}$) at a 68.3\% confidence level. This is the first determination of MAXI J1305-704's spin.

Cheng-Min Zhang, Xiang-Han Cui, Di Li, De-Hua Wang, Shuang-Qiang Wang, Na Wang, Jian-Wei Zhang, Bo Peng, Wei-Wei Zhu, Yi-Yan Yang, Yuan-Yue Pan

21pages, 3 figures, 2 tables, and 5 appendices, published in Universe

The evolutions of a neutron star's rotation and magnetic field (B-field) have remained unsolved puzzles for over half a century. We ascribe the rotational braking torques of pulsar to both components, the standard magnetic dipole radiation (MDR) and particle wind flow ( MDR + Wind, hereafter named MDRW), which we apply to the Crab pulsar (B0531 + 21), the only source with a known age and long-term continuous monitoring by radio telescope. Based on the above presumed simple spin-down torques, we obtain the exact analytic solution on the rotation evolution of the Crab pulsar, together with the related outcomes as described below: (1) unlike the constant characteristic B-field suggested by the MDR model, this value for the Crab pulsar increases by a hundred times in 50~kyr while its real B-field has no change; (2) the rotational braking index evolves from $\sim$3 to 1 in the long-term, however, it drops from 2.51 to 2.50 in $\sim$45 years at the present stage, while the particle flow contributes approximately 25% of the total rotational energy loss rate; (3) strikingly, the characteristic age has the maximum limit of $\sim$10 kyr, meaning that it is not always a good indicator of real age. Furthermore, we discussed the evolutionary path of the Crab pulsar from the MDR to the wind domination by comparing it with the possible wind braking candidate pulsar PSR J1734-3333.

Urmila Chadayammuri, Akos Bogdan, Angelo Ricarte, Priyamvada Natarajan

Submitted to AAS Journals

Dwarf galaxies are considered to be potential ideal test-beds for constraining models of the seeding and tracing of the growth of supermassive and intermediate mass black holes (MBH) via their black hole occupation fraction (BHOF). Disentangling seeding from the confounding effects of mass assembly is, however, challenging. In this work, we use semi-analytical models (SAMs) to probe how various surveys perform at teasing apart different seed and growth scenarios. We check for differences in the measured BHOF given various cuts to black hole mass and AGN luminosity and develop a scheme to robustly compare SAMs, with their intrinsic uncertainties, to X-ray observations. We demonstrate that to tell seeding models apart, we need to detect or model all AGN brighter than $10^{37}\ \rm{erg \ s^{-1}}$ in galaxies of $M_* \sim 10^{8-10} \ \rm{M_{\odot}}$ Shallower surveys, like eRASS, cannot distinguish between seed models even with the compensation of a much larger survey volume. We show that the AMUSE survey strongly favours heavy seed models, growing with empirically motivated growth models either a power-law Eddington Ratio Distribution Function (ERDF) or one in which black hole accretion is tagged to the star-formation rate (AGN-MS). These two growth channels in turn can then be distinguished by the AGN luminosity function at $< 10^{44}\ \rm{erg \ s^{-1}}$. The different models also predict different radio scaling relations, which we quantify using the fundamental plane of black hole activity. We close with recommendations for the design of upcoming multi-wavelength campaigns that can optimally detect MBHs in dwarf galaxies.

Pedro Alonso, Wenting Wang, Jun Zhang, Hekun Li, Shi Shao, Qi Guo, Yanqin He, Cai-Na Hao, Rui Shi

28 pages, 12 figures

Based on independent shear measurements using the DECaLS/DR8 imaging data, we measure the weak lensing signals around isolated central galaxies (ICGs) from SDSS/DR7 at $z\sim0.1$. The projected stellar mass density profiles of surrounding satellite galaxies are further deduced, using photometric sources from the Hyper Suprime-Cam (HSC) survey (pDR3). The signals of ICGs $+$ their extended stellar halos are taken from Wang et al.(2021). All measurements are compared with predictions by the Illustris-TNG300-1 simulation. We find, overall, a good agreement between observation and TNG300. In particular, a correction to the stellar mass of massive observed ICGs is applied based on the calibration of He et al.(2013), which brings a much better agreement with TNG300 predicted lensing signals at $\log_{10}M_\ast/M_\odot>11.1$. In real observation, red ICGs are hosted by more massive dark matter halos, have more satellites and more extended stellar halos than blue ICGs at fixed stellar mass. However, in TNG300 there are more satellites around blue ICGs at fixed stellar mass, and the outer stellar halos of red and blue ICGs are similar. The stellar halos of TNG galaxies are more extended compared with real observed galaxies, especially for blue ICGs with $\log_{10}M_\ast/M_\odot>10.8$. We find the same trend for TNG100 galaxies and for true halo central galaxies. The tensions between TNG and real galaxies might indicate that satellite disruptions are stronger in TNG. In both TNG300 and observation, satellites approximately trace the underlying dark matter distribution beyond $0.1R_{200}$, but the fraction of total stellar mass in TNG300 does not show the same radial distribution as real galaxies.

Marta Volonteri, Melanie Habouzit, Monica Colpi

Submitted to MNRAS

JWST is discovering star forming `candidate' galaxies with photometric redshifts $z>9$ and little attenuation. We model presumptive massive black holes (MBHs) in such galaxies and find that their unobscured emission is fainter than the galaxy starlight in JWST filters, and difficult to be detected via color-color selection, and X-ray and radio observations. Only MBHs overmassive relative to expected galaxy scaling relations, accreting at high Eddington rates, would be detectable. Their discovery would point to the presence of heavy MBH seeds, but care is needed to exclude the existence of lighter seeds as only overmassive MBHs are detectable in this type of galaxies. Conversely, if no overmassive MBHs are hosted in these galaxies, either there are no heavy seeds or they are rare. The most massive/highest redshift candidate galaxies can attain stellar masses in excess of 5e10 Msun by z~6 if they grow along the SFR-mass sequence, and can nurse a MBH growing from ~1e5 Msun up to >3e7 Msun by z~6, to become hosts of some z>6 quasars. Candidate galaxies of log(M_{gal}/Msun)~8 can not grow their putative seeds fast, unless seeds are >1e6 Msun. The number density of the JWST candidate galaxies far outnumbers that of the highest-z quasar hosts and this allows for about only 1 high-redshift quasar every 1000 of these galaxies.

Barbara Catinella, Luca Cortese, Alfred L. Tiley, Steven Janowiecki, Adam B. Watts, Julia J. Bryant, Scott M. Croom, Francesco d'Eugenio, Jesse van de Sande, Joss Bland-Hawthorn, Amelia Fraser-McKelvie, Samuel N. Richards, Sarah M. Sweet, Daniel J. Pisano, Nickolas Pingel, Rebecca A. Koopmann, Dillion Cottrill, Meghan Hill

18 pages, 12 figures. Accepted for publication in MNRAS

We present SAMI-HI, a survey of the atomic hydrogen content of 296 galaxies with integral field spectroscopy available from the SAMI Galaxy Survey. The sample spans nearly 4 dex in stellar mass ($M_\star = 10^{7.4}-10^{11.1}~ \rm M_\odot$), redshift $z<0.06$, and includes new Arecibo observations of 153 galaxies, for which we release catalogues and HI spectra. We use these data to compare the rotational velocities obtained from optical and radio observations and to show how systematic differences affect the slope and scatter of the stellar-mass and baryonic Tully-Fisher relations. Specifically, we show that H$\alpha$ rotational velocities measured in the inner parts of galaxies (1.3 effective radii in this work) systematically underestimate HI global measurements, with HI/H$\alpha$ velocity ratios that increase at low stellar masses, where rotation curves are typically still rising and H$\alpha$ measurements do not reach their plateau. As a result, the H$\alpha$ stellar mass Tully-Fisher relation is steeper (when $M_\star$ is the independent variable) and has larger scatter than its HI counterpart. Interestingly, we confirm the presence of a small fraction of low-mass outliers of the H$\alpha$ relation that are not present when HI velocity widths are used and are not explained by "aperture effects". These appear to be highly disturbed systems for which H$\alpha$ widths do not provide a reliable estimate of the rotational velocity. Our analysis reaffirms the importance of taking into account differences in velocity definitions as well as tracers used when interpreting offsets from the Tully-Fisher relation, at both low and high redshifts and when comparing with simulations.

Adam B. Watts, Luca Cortese, Barbara Catinella, Chris Power, Amelia Fraser-McKelvie, Julia J. Bryant, Scott M. Croom, Jesse van de Sande, Brent Groves

11 pages, 6 figures, 1 appendix

Observations of the neutral atomic hydrogen (HI) gas in galaxies are predominantly spatially unresolved, in the form of a global HI spectral line. There has been substantial work on quantifying asymmetry in global HI spectra (`global HI asymmetry'), but due to being spatially unresolved, it remains unknown what physical regions of galaxies the asymmetry traces, and whether the other gas phases are affected. Using optical integral field spectrograph (IFS) observations from the Sydney AAO Multi-object IFS (SAMI) survey for which global HI spectra are also available (SAMI-HI), we study the connection between asymmetry in galaxies' ionised and neutral gas reservoirs to test if and how they can help us better understand the origin of global HI asymmetry. We reconstruct the global H$\alpha$ spectral line from the IFS observations and find that, while some global H$\alpha$ asymmetries can arise from disturbed ionised gas kinematics, the majority of asymmetric cases are driven by the distribution of H$\alpha$-emitting gas. When compared to the HI, we find no evidence for a relationship between the global H$\alpha$ and HI asymmetry. Further, a visual inspection reveals that cases where galaxies have qualitatively similar H$\alpha$ and HI spectral profiles can be spurious, with the similarity originating from an irregular 2D H$\alpha$ flux distribution. Our results highlight that comparisons between global H$\alpha$ and HI asymmetry are not straightforward, and that many global HI asymmetries trace disturbances that do not significantly impact the central regions of galaxies.

Zhuyun Zhuang, Nicha Leethochawalit, Evan N. Kirby, J. W. Nightingale, Charles C. Steidel, Karl Glazebrook, Tania M. Barone, Hannah Skobe, Sarah M. Sweet, Themiya Nanayakkara, Rebecca J. Allen, Keerthi Vasan G. C., Tucker Jones, Glenn G. Kacprzak, Kim-Vy H. Tran, Colin Jacobs

18 pages, 11 figures, submitted to ApJ; comments welcome

Gravitational lenses can magnify distant galaxies, allowing us to discover and characterize the stellar populations of intrinsically faint, quiescent galaxies that are otherwise extremely difficult to directly observe at high redshift from ground-based telescopes. Here, we present the spectral analysis of two lensed, quiescent galaxies at $z\gtrsim 1$ discovered by the ASTRO 3D Galaxy Evolution with Lenses survey: AGEL1323 ($M_*\sim 10^{11.1}M_{\odot}$, $z=1.016$, $\mu \sim 14.6$) and AGEL0014 ($M_*\sim 10^{11.3}M_{\odot}$, $z=1.374$, $\mu \sim 4.3$). We measured the age, [Fe/H], and [Mg/Fe] of the two lensed galaxies using deep, rest-frame-optical spectra (S/N $\gtrsim$ 40\AA$^{-1}$) obtained on the Keck I telescope. The ages of AGEL1323 and AGEL0014 are $5.6^{+0.8}_{-0.8}$ Gyr and $3.1^{+0.8}_{-0.3}$ Gyr, respectively, indicating that most of the stars in the galaxies were formed less than 2 Gyr after the Big Bang. Compared to nearby quiescent galaxies of similar masses, the lensed galaxies have lower [Fe/H] and [Mg/H]. Surprisingly, the two galaxies have comparable [Mg/Fe] to similar-mass galaxies at lower redshifts, despite their old ages. Using a simple analytic chemical evolution model connecting the instantaneously recycled element Mg with the mass-loading factors of outflows averaged over the entire star formation history, we found that the lensed galaxies may have experienced enhanced outflows during their star formation compared to lower-redshift galaxies, which may explain why they quenched early.

Active regions (ARs) are typical magnetic structures found in the solar atmosphere. We calculate several magnetohydrostatic (MHS) equilibrium models that include the effect of a finite plasma-$\beta$ and gravity and that are representative of these structures in three dimensions. The construction of the models is based on the use of two Euler potentials, $\alpha$ and $\beta$, that represent the magnetic field as ${\bf B}=\nabla \alpha \times \nabla \beta$. The ideal MHS nonlinear partial differential equations are solved numerically using finite elements in a fixed 3D rectangular domain. The boundary conditions are initially chosen to correspond to a potential magnetic field (current-free) with known analytical expressions for the corresponding Euler potentials. The distinctive feature is that we incorporate the effect of shear by progressively deforming the initial potential magnetic field. This procedure is quite generic and allows us to generate a vast variety of MHS models. The thermal structure of the ARs is incorporated through the dependence of gas pressure and temperature on the Euler potentials. Using this method we achieve the characteristic hot and over-dense plasma found in ARs, but we demonstrate that the method can also be applied to study configurations with open magnetic field lines. Furthermore, we investigate basic topologies that include neutral lines. Our focus is on the force balance of the structures and we do not consider the energy balance in the constructed models. In addition, we address the difficult question of the stability of the calculated 3D models. We find that if the plasma is convectively stable, then the system is not prone in general to develop magnetic Rayleigh-Taylor instabilities.

David Dahiya, Deepak Jain

16 pages, 4 Figures, 3 Tables

We revisit the epoch of cosmic speed-up characterized by the redshift of transition from a decelerated to an accelerated phase. This redshift is termed the transition redshift ($z_t$). We use the spatially Flat and Non-Flat variants of the most common $\Lambda$CDM and XCDM models to put constraints on the transition redshift along with the other model parameters. The data for this analysis comes from the recent and updated Pantheon+ Supernova dataset and the Hubble parameter measurements obtained from Cosmic Chronometers. We consider both datasets with their respective covariance matrices incorporating all kinds of statistical and systematic uncertainties. We observe that using the combined datasets of H(z) and SNe, the best fit value of transition redshift lies in the range $0.61 < z_t < 0.82$ for all four dark energy models. Incidentally, we observe a positive curvature for the Non-Flat models and correlations between several model parameters.

Yuan-Ming Hsu, Hiroyuki Hirashita, Yen-Hsing Lin, Peter Camps, Maarten Baes

11 pages, 6 figures, accepted for publication in MNRAS

Dust emission from high-redshift galaxies gives us a clue to the origin and evolution of dust in the early Universe. Previous studies have shown that different sources of dust (stellar dust production and dust growth in dense clouds) predict different ultraviolet (UV) extinction curves for galaxies at $z\sim 7$ but that the observed attenuation curves depend strongly on the geometry of dust and star distributions. Thus, we perform radiative transfer calculations under the dust-stars geometries computed by a cosmological hydrodynamic simulation (IllustrisTNG). This serves to investigate the dust attenuation curves predicted from `realistic' geometries. We choose objects with stellar mass and star formation rate appropriate for Lyman break galaxies at $z\sim 7$. We find that the attenuation curves are very different from the original extinction curves in most of the galaxies. This makes it difficult to constrain the dominant dust sources from the observed attenuation curves. We further include infrared dust emission in the analysis and plot the infrared excess (IRX)-UV spectral slope ($\beta$) diagram. We find that different sources of dust cause different IRX-$\beta$ relations for the simulated galaxies. In particular, if dust growth is the main source of dust, a variation of dust-to-metal ratio causes a more extended sequence with smaller IRX in the IRX-$\beta$ diagram. Thus, the comprehensive analysis of the abundances of dust and metals, the UV slope, and the dust emission could provide a clue to the dominant dust sources in the Universe.

Blaz Bortolato, Jernej F. Kamenik, Michele Tammaro

28 pages + 16 page Appendix, 27 figures

We apply statistical inference on the Pierre Auger Open Data to discern for the first time the full mass composition of cosmic rays at different energies. Working with longitudinal electromagnetic profiles of cosmic ray showers, in particular their peaking depths $X_{\rm max}$, we employ central moments of the $X_{\rm max}$ distributions as features to discriminate between different shower compositions. We find that already the first few moments entail the most relevant information to infer the primary cosmic ray mass spectrum. Our approach, based on an unbinned likelihood, allows us to consistently account for sources of statistical uncertainties due to finite datasets, both measured and simulated, as well as systematic effects. Finally, we provide a quantitative comparison of different high energy hadronic interaction models available in the atmospheric shower simulation codes.

Supriyo Ghosh (TIFR, Mumbai), Joe P. Ninan (TIFR, Mumbai), Devendra K. Ojha (TIFR, Mumbai), Saurabh Sharma (ARIES, Nainital)

10 pages, 6 figures, accepted for publication in the Special Issue of Journal of Astrophysics & Astronomy, 2022, Star formation studies in context of NIR instruments on 3.6m DOT, held at ARIES, Nainital during 4-7, May, 2022

The TIFR-ARIES Near Infrared Spectrometer (TANSPEC) instrument provides simultaneous wavelength coverage from 0.55 to 2.5 micron, mounted on India's largest ground-based telescope, 3.6-m Devasthal Optical Telescope at Nainital, India. The TANSPEC offers three modes of observations, imaging with various filters, spectroscopy in the low-resolution prism mode with derived R~ 100-400 and the high-resolution cross-dispersed mode (XD-mode) with derived median R~ 2750 for a slit of width 0.5 arcsec. In the XD-mode, ten cross-dispersed orders are packed in the 2048 x 2048 pixels detector to cover the full wavelength regime. As the XD-mode is most utilized as well as for consistent data reduction for all orders and to reduce data reduction time, a dedicated pipeline is at the need. In this paper, we present the code for the TANSPEC XD-mode data reduction, its workflow, input/output files, and a showcase of its implementation on a particular dataset. This publicly available pipeline pyTANSPEC is fully developed in Python and includes nominal human intervention only for the quality assurance of the reduced data. Two customized configuration files are used to guide the data reduction. The pipeline creates a log file for all the fits files in a given data directory from its header, identifies correct frames (science, continuum and calibration lamps) based on the user input, and offers an option to the user for eyeballing and accepting/removing of the frames, does the cleaning of raw science frames and yields final wavelength calibrated spectra of all orders simultaneously.

R. A. Main, S. Bethapudi, V. R. Marthi, M. L. Bause, D. Z. Li, H.-H. Lin, L. G. Spitler, R. S. Wharton

5 pages, 5 Figures, submitted to MNRAS Letters

Compact radio sources exhibit scintillation, an interference pattern arising from propagation through inhomogeneous plasma, where scintillation patterns encode the relative distances and velocities of the source, scattering material, and Earth. In Main et al. 2022, we showed that the scintillation velocity of the repeating fast radio burst FRB20201124A can be measured by correlating pairs of burst spectra, and suggested that the scattering was nearby the Earth at $\sim0.4\,$kpc from the low values of the scintillation velocity and scattering timescale. In this work, we have measured the scintillation velocity at 10 epochs spanning a year, observing an annual variation which strongly implies the screen is within the Milky Way. Modelling the annual variation with a 1D anisotropic or 2D isotropic screen results in a screen distance $d_{l} = 0.24\pm0.04\,$pc or $d_{l} = 0.37\pm0.07\,$pc from Earth respectively, possibly associated with the Local Bubble or the edge of the Orion-Eridanus Superbubble. Continued monitoring, and using measurements from other telescopes particularly at times of low effective velocity will help probe changes in screen properties, and distinguish between screen models. Where scintillation of an FRB originates in its host galaxy or local environment, these techniques could be used to detect orbital motion, and probe the FRB's local ionized environment.

Richard Qiu, Angelo Ricarte, Ramesh Narayan, George N. Wong, Andrew Chael, Daniel Palumbo

24 pages, 27 figures

We introduce a new library of 535,194 model images of the supermassive black holes and Event Horizon Telescope (EHT) targets Sgr A* and M87*, computed by performing general relativistic radiative transfer calculations on general relativistic magnetohydrodynamics simulations. Then, to infer underlying black hole and accretion flow parameters (spin, inclination, ion-to-electron temperature ratio, and magnetic field polarity), we train a random forest machine learning model on various hand-picked polarimetric observables computed from each image. Our random forest is capable of making meaningful predictions of spin, inclination, and the ion-to-electron temperature ratio, but has more difficulty inferring magnetic field polarity. To disentangle how physical parameters are encoded in different observables, we apply two different metrics to rank the importance of each observable at inferring each physical parameter. Details of the spatially resolved linear polarization morphology stand out as important discriminators between models. Bearing in mind the theoretical limitations and incompleteness of our image library, for the real M87* data, our machinery favours high-spin retrograde models with large ion-to-electron temperature ratios. Due to the time-variable nature of these targets, repeated polarimetric imaging will further improve model inference as the EHT and next-generation (EHT) continue to develop and monitor their targets.

Luca Caloni, Serena Giardiello, Margherita Lembo, Martina Gerbino, Giulia Gubitosi, Massimiliano Lattanzi, Luca Pagano

We perform a comprehensive study of the signatures of Lorentz violation in electrodynamics on the Cosmic Microwave Background (CMB) anisotropies. In the framework of the minimal Standard Model Extension (SME), we consider effects generated by renormalizable operators, both CPT-odd and CPT-even. These operators are responsible for sourcing, respectively, cosmic birefringence and circular polarization. We propagate jointly the effects of all the relevant Lorentz-violating parameters to CMB observables and provide constraints with the most recent CMB datasets. We bound the CPT-even coefficient to $k_{F,E+B} < 2.31 \times 10^{-31}$ at 95\% CL. This improves previous CMB bounds by one order of magnitude. The limits we obtain on the CPT-odd coefficients, i.e. $|k_{(V)00}^{(3)}| < 1.54 \times 10^{-44} \; {\rm GeV}$ and $|\mathbf{k_{AF}}| < 0.74 \times 10^{-44} \; {\rm GeV}$ at 95\% CL, are respectively one and two orders of magnitude stronger than previous CMB-based limits, superseding also bounds from non-CMB searches. This analysis provides the strongest constraints to date on CPT-violating coefficients in the minimal SME from CMB searches.

Caroline Guandalin, Jade Piat, Chris Clarkson, Roy Maartens

13 pages, 3 figures, 2 tables

The Cosmological Principle is part of the foundation that underpins the standard model of the Universe. In the era of precision cosmology, when stress tests of the standard model are uncovering various tensions and possible anomalies, it is critical to check the viability of this Principle. A key test is the consistency between the kinematic dipoles of the cosmic microwave background and of the large-scale matter distribution. Results using radio continuum and quasar samples indicate a rough agreement in the directions of the two dipoles, but a larger than expected amplitude of the matter dipole. The resulting tension with the radiation dipole has been estimated at $\sim 5\sigma$ for some cases, suggesting a potential new cosmological tension and a possible violation of the Cosmological Principle. However, the standard formalism for predicting the dipole in the 2-dimensional projection of sources overlooks possible evolution effects in the luminosity function. In fact, radial information from the luminosity function is necessary for a correct projection of the 3-dimensional source distribution. Using a variety of current models of the quasar luminosity function, we show that neglecting redshift evolution can significantly overestimate the relative velocity amplitude. While the models we investigate are consistent with each other and with current data, the dipole derived from these, which depends on derivatives of the luminosity function, can disagree by more than $3\sigma$. This theoretical systematic bias needs to be resolved before robust conclusions can be made about a new cosmic tension.

M. Alhallak (Damascus Univ.), N. Chamoun (HIAST), M. S. Eldaher (Damascus Univ.)

pdflatex, 10 pages (epjc style), 5 tables, 6 figures

Power-Law inflation with scale factor $a \propto t^m$ is investigated in the context of warm inflation. The treatment is performed in the weak and strong dissipation limits. In addition, we discuss the common three cases for the thermal dissipation coefficient $\Gamma(T)$. We compare the theoretical results of the Power-Law model within warm inflation with the observational constraints from Planck $2018$ and BICEP/Keck 2018, as presented by the tensor-to-scalar ratio $r$ and spectral index $n_s$. The model results agree largely with the observations for most of the $\Gamma(T)$ cases. Furthermore, in order to addresses the problem of exiting the inflationary epoch, we suggest a perturbed modification to the power-law definition so that it becomes affine, and find that this small change is sufficient to provide a successful exit scenario with suitable e-foldings number. Finally, we examine this perturbation ansatz within the context of cold inflation with exponential potential, and we find that it can accommodate the observational data with sufficient e-foldings. Thus, the study successfully rescues the Power-Law inflation and the exponential potential in both warm and cold inflation contexts.

J. A. Hsu, J. C. Jiang, H. Xu, K. J. Lee, R. X. Xu

10 pages, 9 figures, accepted by MNRAS

More than 100 rotating radio transients (RRATs) have been discovered since 2006. However, it is unclear whether RRATs radiate in the nulling states. PSR J0628+0909 has been classified as an RRAT. In this paper, we study the single pulses and integrated pulse profile of PSR J0628+0909 to check whether we can detect pulsed radio emission in the nulling states. We also aim to study the polarization of the RRAT and its relationship to the general pulsar population. We used the Five-hundred-meter Aperture Spherical radio Telescope (FAST) to observe PSR J0628+0909 in the frequency range from 1.0 to 1.5 GHz. We searched for strong single pulses and looked for pulsed emission in the RRAT nulling states. Polarisation profiles, the single-pulse energy distribution, and waiting-time statistics were measured. The Faraday rotation measure and dispersion measure values are updated with the current observation. The single-pulse polarisation behaviours show great diversity, similar to the case of pulsars. Based on the integrated pulse profile and single-pulse energy statistics, we argue that continuous pulsar-like emission exists in addition to the transient-like burst emission for PSR J0628+0909. We find that the pulse waiting time is not correlated with the pulse energy and conclude that the strong transient emission of RRAT is not generated by the energy store-release mechanism.

Cristiano Longarini, Giuseppe Lodato, Giuseppe Bertin, Philip J. Armitage

14 pages, 8 figures, 1 appendix, accepted for publication on MNRAS

Recent observations show that planet formation is already underway in young systems, when the protostar is still embedded into the molecular cloud and the accretion disc is massive. In such environments, the role of self gravity (SG) and gravitational instability (GI) is crucial in determining the dynamical evolution of the disc. In this work, we study the dynamical role of drag force in self-gravitating discs as a way to form planetesimals in early protoplanetary stages. We obtain the dispersion relation for density-wave perturbations on a fluid composed of two phases (gas and dust) interacting through the common gravitation field and the mutual drag force, and we find that the stability threshold is determined by three parameters: the local dust-to-gas density ratio, the dust relative temperature and the relevant Stokes number. In a region of parameters space, where young protoplanetary discs are likely to be found, the instability can be dust driven, occurring at small wavelengths. In this regime, the Jeans mass is much smaller than the one predicted by the standard gravitational instability model. This mechanism can be a viable way to form planetary cores in protostellar discs, since their predicted mass is about 10 Earth masses.

Vicki L. Herde, Phillip C. Chamberlin, Don Schmit, Souvik Bose, Adrian Daw, Ryan O. Milligan, Vanessa Polito

17 pages, 9 figures, presented at the AGU Fall 2022 conference, Submitted to AAS Journal

We have developed an algorithm to identify solar spicules in the first-ever systematic survey of on-disk spicules using exclusively Mg II spectral observations. Using this algorithm we identify 2219 events in three IRIS datasets with unique solar feature targets spanning a total of 300 minutes: 1) an active region, 2) decayed active region/active network, and 3) a coronal hole. We present event statistics and relate occurrence rates to underlying photospheric magnetic field strength. This method identifies spicule event densities and occurrence rates similar to previous studies performed using H{\alpha} and Ca II observations of active regions. Additionally, this study identifies spicule-like events at very low rates at magnetic field intensities below 20 Gauss and increasing significantly between 100-200 Gauss in active regions and above 20 Gauss in coronal holes, which can be used to inform future observation campaigns. This information can be be used to help characterize spicules over their full lifetime, and compliments existing H-{\alpha} spectral capabilities and upcoming Ly-{\alpha} spectral observations on the SNIFS Sounding Rocket. In total, this study presents a method for detecting solar spicules using exclusively Mg II spectra, and provides statistics for spicule occurrence in Mg II wavelengths with respect to magnetic field strength for the purpose of predicting spicule occurrences.

We study evolution of hydrogenated amorphous carbon (HAC) grains under harsh UV radiation in photo-dissociation regions (PDRs) near young massive stars. Our aim is to evaluate the impact of the HAC grains on formation of observed small hydrocarbons: C$_2$H, C$_2$H$_2$, C$_3$H$^+$, C$_3$H, C$_3$H$_2$, C$_4$H, in PDRs. We developed a microscopic model of the HAC grains based on available experimental results. The model includes processes of photo- and thermodesorption, accretion of hydrogen and carbon atoms and subsequent formation of carbonaceous mantle on dust surface. H$_2$, CH$_4$, C$_2$H$_2$, C$_2$H$_4$, C$_2$H$_6$, C$_3$H$_4$, C$_3$H$_6$, C$_3$H$_8$ are considered as the main fragments of the HAC photo-destruction. We simulated evolution of the HAC grains under the physical conditions of two PDRs, the Orion Bar and the Horsehead nebula. We estimated the production rates of the HAC fragments in gas phase chemical reactions and compared them with the production rates of fragments due to the HAC destruction. The latter rates may dominate under some conditions, namely, at A$_V$=0.1 in both PDRs. We coupled our model with the gas-grain chemical model MONACO and calculated abundances of observed small hydrocarbons. We conclude that the contribution of the HAC destruction fragments to chemistry is not enough to match the observed abundances, although it increases the abundances by several orders of magnitude in the Orion Bar at A$_V$=0.1. Additionally, we found that the process of carbonaceous mantle formation on dust surface can be an inhibitor for the formation of observed small hydrocarbons in PDRs.

Dragan Huterer

Invited review for The Astronomy and Astrophysics Review, aimed at a non-expert; 29 pages + references

We review one of the most fruitful areas in cosmology today that bridge theory and data - the temporal growth of large-scale structure. We go over the growth's physical foundations, and derive its behavior in simple cosmological models. While doing so, we explain how measurements of growth can be used to understand theory. We then review how some of the most mature cosmological probes - galaxy clustering, gravitational lensing, the abundance of clusters of galaxies, cosmic velocities, and cosmic microwave background - can be used to probe the growth of structure. We report the current constraints on growth, which are summarized as measurements of the parameter combination $f\sigma_8$ as a function of redshift, or else as the mass fluctuation amplitude parameter $S_8$. We finally illustrate several statistical approaches, ranging from the "growth index" parameterization to more general comparisons of growth and geometry, that can sharply test the standard cosmological model and indicate the presence of modifications to general relativity.

Dionysios Gakis, Konstantinos N. Gourgouliatos

11 pages, 9 figures, revised afer reviewer's comments

Recent observations by TESS revealed the existence of circumbinary planets in the systems of TOI-1338 and TIC-172900988. The purpose of this work is to model the planetary orbits in these two systems and study them under the perspective of previous theoretical models. Each planet's distance from the barycenter through time is simulated using n-body integrations and is compared with outcomes from a semi-analytic, a geometric and a Keplerian-based approach. Furthermore, we infer the most prominent frequencies of both planets' orbits induced by the central binaries. We confirm that both systems appear to be stable. Lastly, we examine the implications of an additional candidate planet in TOI-1338 system finding that an extra, 48 M$_\oplus$ planet that has been hinted from observations could be located at 0.8 AU without generating any radical changes to the orbits of the other members of the system.

A. Arellano Ferro, F. C. Rojas Galindo, I. H. Bustos Fierro, S. Muneer, M.A. Yepez, Sunetra Giridhar

16 pages, 5 Tables, 8 Figures. Accepted by MNRAS November 8, 2022

A five year CCD photometric \emph{VI} time-series of NGC 7006 is employed to perform a detailed analysis of the known population of variable stars. In the process we have corrected inconsistent classifications, sky coordinates and found ten new cluster member variables. An independent reddening estimate with a value $E(B-V)=0.08\pm0.05$ is made. Using Fourier decompositions of RR Lyrae light curves and well established calibrations, the cluster mean metallicity and distance [Fe/H]$_{\rm ZW}= -1.53\pm0.15$ and $41.2\pm1.4$ kpc are estimated based on an extended sample of cluster member RRab stars. Using the $Gaia$-DR3 data we performed an extensive membership analysis which leads to a clean Colour-Magnitude diagram, and hence to the identification of variables that are likely field stars, and to considerations on the variables distribution in the Horizontal Branch (HB). A double mode RR Lyrae and three CW star are discussed. The origin of CW stars from precursors in the blue tail of the HB with very thin ($\sim 0.06 \pm 0.01 M_{\odot}$) envelopes is argued. Our models indicate that the Main Sequence predecessor of RR Lyrae stars had a mass of 0.82-0.85 $M_{\odot}$ and lost about 25-35\% of its mass during the red giant branch events before settling in the HB some 12-13.5 Gyrs later.

Sudip Mandal, Hardi Peter, Lakshmi Pradeep Chitta, Regina A. Cuadrado, Udo Schühle, Luca Teriaca, Sami K. Solanki, Louise Harra, David Berghmans, Frédéric Auchère, Susanna Parenti, Andrei N. Zhukov, Éric Buchlin, Cis Verbeeck, Emil Kraaikamp, Luciano Rodriguez, David M. Long, Conrad Schwanitz, Krzysztof Barczynski, Gabriel Pelouze, Philip J. Smith, Wei Liu, Mark C. Cheung

Accepted for publication in A&A Letters. Event movie can be downloaded from this https URL

The solar chromosphere hosts a wide variety of transients, including dynamic fibrils (DFs) that are characterised as elongated, jet-like features seen in active regions, often through H$\alpha$ diagnostics. So far, these features have been difficult to identify in coronal images primarily due to their small size and the lower spatial resolution of the current EUV imagers. Here we present the first unambiguous signatures of DFs in coronal EUV data using high-resolution images from the Extreme Ultraviolet Imager (EUI) on board Solar Orbiter. Using the data acquired with the 174~{\AA} High Resolution Imager (HRI$_{EUV}$) of EUI, we find many bright dot-like features (of size 0.3-0.5 Mm) that move up and down (often repeatedly) in the core of an active region. In a space-time map, these features produce parabolic tracks akin to the chromospheric observations of DFs. Properties such as their speeds (14 km~s$^{-1}$), lifetime (332~s), deceleration (82 m~s$^{-2}$) and lengths (1293~km) are also reminiscent of the chromospheric DFs. The EUI data strongly suggest that these EUV bright dots are basically the hot tips (of the cooler chromospheric DFs) that could not be identified unambiguously before because of a lack of spatial resolution.

Yoav Rotman, Thaddeus D. Komacek, Geronimo L. Villanueva, Thomas J. Fauchez, Erin M. May

11 pages, 4 figures, 1 table, accepted for publication in ApJL

Terrestrial exoplanets such as TRAPPIST-1e will be observed in a new capacity with JWST/NIRSpec, which is expected to be able to detect CO$_2$, CH$_4$, and O$_2$ signals, if present, with multiple co-added transit observations. The CO$_2$-CH$_4$ pair in particular is theorized to be a potential biosignature when inferred to be in chemical disequilibrium. Here, we simulate TRAPPIST-1e's atmosphere using the ExoCAM General Circulation Model (GCM), assuming an optimistic haze-free, tidally locked planet with an aquaplanet surface, with varying atmospheric compositions from $10^{-4}$ bar to 1 bar of partial CO$_2$ pressure with 1 bar of background N$_2$. We investigate cases both with and without a modern Earth-like CH$_4$ mixing ratio to examine the effect of CO$_2$ and CH$_4$ on the transmission spectrum and climate state of the planet. We demonstrate that in the optimistic haze-free cloudy case, H$_2$O, CO$_2$, and CH$_4$ could all be detectable in less than 50 transits within an atmosphere of 1 bar N$_2$ and 10 mbar CO$_2$ during JWST's lifespan with NIRSpec as long as the noise floor is $\lesssim$ 10 ppm. We find that in these optimistic cases, JWST may be able to detect potential biosignature pairs such as CO$_2$-CH$_4$ in TRAPPIST-1e's atmosphere across a variety of atmospheric CO$_2$ content, and that temporal climate variability does not significantly affect spectral feature variability for NIRSpec PRISM.

Ke Gao, Lei-Hua Liu, Mian Zhu

10 pages, 11 figures, references updated, typos corrected

In this paper, we investigate the microlensing effects of blackhole-like wormholes. We evaluate the deflection angle upon the second order under weak field approximation with Gauss-Bonnet theorem. We elaborate on the deflection angle of the Ellis-Bronnikov wormhole as an example. Following the same procedure, we study the magnification of three typical wormholes (WH): Schwarzschild WH, Kerr-like WH, and RN WH, as well as their blackhole correspondence. We find that the prograde case of Kerr-like metric will lead to the multi-peaks of magnification as the mass part is compatible with the charge part. Moreover, the first two gentle peaks of Kerr blackhole are larger than the wormhole case by one order of magnitude, while the main peak of Kerr blackholes and wormholes are of the same order. For other cases, the behavior of magnification wormholes and their corresponding blackholes is similar. Our result may shed new light on exploring compact objects through the microlensing effect.

Peizhi Du, Daniel Egaña-Ugrinovic, Rouven Essig, Mukul Sholapurkar

12 pages, 7 figures

Dopant atoms in semiconductors can be ionized with $\sim10$ meV energy depositions, allowing for the design of low-threshold detectors. We propose using doped semiconductor targets to search for sub-MeV dark matter scattering or sub-eV dark matter absorption on electrons. Currently unconstrained cross sections could be tested with a 1 g-day exposure in a doped detector with backgrounds at the level of existing pure semiconductor detectors, but improvements would be needed to probe the freeze-in target. We discuss the corresponding technological requirements and lay out a possible detector design.

Debasish Borah, Suruj Jyoti Das, Nobuchika Okada

16 pages, 3 captioned figures

We propose a mechanism for cogenesis of baryon and dark matter (DM) in the universe via the Affleck-Dine (AD) route. An AD field which breaks the lepton number symmetry, leads to the generation of lepton asymmetry by virtue of its cosmic evolution, which then gets transferred into lepton and dark sectors. While the lepton asymmetry gets converted into baryon asymmetry via sphalerons, the dark sector asymmetry leads to the final DM abundance with the symmetric part being annihilated away due to resonantly enhanced annihilation, which we choose to be provided by a gauged $B-L$ portal. Stringent constraints from DM direct detection forces DM and $B-L$ gauge boson masses to be light, in the few GeV ballpark. While some part of the model parameter space is already ruled out, the remaining parameter is within sensitivity of laboratory as well as cosmology based experiments. The AD field also plays the role of inflaton with the required dynamics by virtue of its non-minimal coupling to gravity, consistent with observations.

The population of stellar origin black hole binaries (SOBHBs) detected by existing ground-based gravitational wave detectors is an exciting target for the future space-based Laser Interferometer Space Antenna (LISA). LISA is sensitive to signals at significantly lower frequencies than ground-based detectors. SOBHB signals will thus be detected much earlier in their evolution, years to decades before they merge. The mergers will then occur in the frequency band covered by ground-based detectors. Observing SOBHBs years before merger can help distinguish between progenitor models for these systems. We present a new Bayesian parameter estimation algorithm for LISA observations of SOBHBs that uses a time-frequency (wavelet) based likelihood function. Our technique accelerates the analysis by several orders of magnitude compared to the standard frequency domain approach and allows for an efficient treatment of non-stationary noise.

Arpan Kar, Tanmoy Kumar, Sourov Roy, Jure Zupan

24 pages, 12 figures

Relativistic axions produced in decays of ${\mathcal O}(10^{-7}-10^{-2}$ $\text{eV})$ dark matter (DM) partially convert to photons after traversing the galactic magnetic field, giving rise to a signal observable by the Square Kilometer Array (SKA) radio telescope. We show that for axions lighter than a few $\times$ $10^{-13}$ eV a 100\,h SKA observation of the local dwarf galaxy Seg I would probe parameter space not constrained by stellar cooling and cosmological observations, with sensitivity several orders of magnitude better than the planned dedicated axion dark matter search experiments. We quantify the uncertainties in the SKA sensitivity projections due to two effects that enhance the photon flux: the presence of turbulent magnetic fields inside the galaxy, and the Bose enhancement of the DM decays to axions, where the latter, in particular, warrants further study.

We revisit the problem of inertial r-modes in stratified stars, drawing on a more precise description of the composition stratification in a mature neutron star. The results highlight issues with the traditional approach to the problem, leading us to rethink the computational strategy for r-modes of non-barotropic neutron stars. We outline two strategies for dealing with the problem. For moderate to slowly rotating neutron stars the only viable alternative may be to approach the problem numerically from the outset, while a meaningful slow-rotation calculation can be carried out for the fastest known spinning stars (which may be close to being driven unstable by the emission of gravitational waves). We demonstrate that the latter approach leads to a problem close, but not identical, to that for barotropic inertial modes. We also suggest that these reformulations of the problem likely resolve the long-standing problem of singular behaviour associated with a co-rotation point in rotating relativistic neutron stars. This issue needs to be resolved in order to guide future gravitational-wave searches.

The only r-modes that exist in a globally barotropic, rotating, Newtonian star are the fundamental $l = |m|$ solutions, where $l$ and $m$ are the indices of the spherical harmonic $Y_l^m$ that describe the mode's angular dependence. This is in stark contrast to a stellar model that is non-barotropic throughout its interior, which hosts all the $l \geq |m|$ perturbations including radial overtones. In reality, neutron stars are stratified with locally barotropic regions. Therefore, we explore how stratification alters a star's ability to support r-modes. We consider the globally stratified case and examine the behaviour of the modes as the star gets close to barotropicity. In this limit, we find that all but the fundamental $l = |m|$ perturbations change character and become generic inertial modes. Restricting the analysis to $l = |m|$ perturbations, we develop the r-mode equations in order to consider stellar models that exhibit local barotropicity. Our results for such models show that the r-mode overtones diverge and join the inertial modes. This suggests that neutron stars can only support the fundamental $l = |m|$ r-modes.

Hao Hao, Qing-Yuan Zhao, Yang-Hui Huang, Jie Deng, Hui Wang, Jia-Wei Guo, Shi Chen, Sai-Ying Ru, Zhen Liu, Yi-Jin Zhou, Shun-Hua Wang, Chao Wan, Hao Liu, Zhi-Jian Li, Hua-bing Wang, Xue-Cou Tu, La-Bao Zhang, Xiao-Qing Jia, Jian Chen, Lin Kang, Pei-Heng Wu

Deep space explorations require transferring huge amounts of data quickly from very distant targets. Laser communication is a promising technology that can offer a data rate of magnitude faster than conventional microwave communication due to the fundamentally narrow divergence of light. This study demonstrated a photon-sensitive receiver prototype with over Gigabit data rate, immunity to strong background photon noise, and simultaneous tracking ability. The advantages are inherited from a joint-optimized superconducting nanowire single-photon detector (SNSPD) array, designed into a four-quadrant structure with each quadrant capable of resolving six photons. Installed in a free-space coupled and low-vibration cryostat, the system detection efficiency reached 72.7%, the detector efficiency was 97.5%, and the total photon counting rate was 1.6 Gcps. Additionally, communication performance was tested for pulse position modulation (PPM) format. A series of signal processing methods were introduced to maximize the performance of the forward error correction (FEC) code. Consequently, the receiver exhibits a faster data rate and better sensitivity by about twofold (1.76 photons/bit at 800 Mbps and 3.40 photons/bit at 1.2 Gbps) compared to previously reported results (3.18 photon/bit at 622 Mbps for the Lunar Laser Communication Demonstration). Furthermore, communications in strong background noise and with simultaneous tracking ability were demonstrated aimed at the challenges of daylight operation and accurate tracking of dim beacon light in deep space scenarios.