Abstract visibility
Abstract text size

Papers for Monday, Apr 17 2023

Papers with local authors

Lizvette Villafaña, Peter R. Williams, Tommaso Treu, Brendon J. Brewer, Aaron J. Barth, Vivian U, Vardha N. Bennert, Hengxiao Guo, Misty C. Bentz, Gabriela Canalizo, Alexei V. Filippenko, Elinor Gates, Michael D. Joner, Matthew A. Malkan, Jong-Hak Woo, Bela Abolfathi, Thomas Bohn, K. Azalee Bostroem, Andrew Brandel, Thomas G. Brink, Sanyum Channa, Maren Cosens, Edward Donohue, Goni Halevi, Carol E. Hood, J. Chuck Horst, Maxime de Kouchkovsky, Benjamin Kuhn, Douglas C. Leonard, Raul Michel, Melanie Kae B. Olaes, Daeseong Park, Jordan N. Runco, Remington O. Sexton, Isaac Shivvers, Chance L. Spencer, Benjamin E. Stahl, Samantha Stegman, Jonelle L. Walsh, WeiKang Zheng

17 pages, 10 figures

0 votes
Paper 9 — arXiv:2304.06764
0 votes
Paper 9 — arXiv:2304.06764

We combine our dynamical modeling black hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors, f, used in traditional reverberation mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have utilized our methods. Aiming to improve the precision of black hole mass estimates, as well as uncover any regularities in the behavior of the broad-line region (BLR), we search for correlations between f and other AGN/BLR parameters. We find (i) evidence for a correlation between the virial coefficient log10(fmean,{\sigma}) and black hole mass, (ii) marginal evidence for a similar correlation between log10(frms,{\sigma}) and black hole mass, (iii) marginal evidence for an anti-correlation of BLR disk thickness with log10(fmean,FWHM)and log10(frms,FWHM), and (iv) marginal evidence for an anti-correlation of inclination angle with log10(fmean,FWHM), log10(frms,{\sigma}), and log10(fmean,{\sigma}). Lastly, we find marginal evidence for a correlation between line-profile shape, when using the root-meansquare spectrum, log10(FWHM/{\sigma})rms, and the virial coefficient, log10(frms,{\sigma}), and investigate how BLR properties might be related to line-profile shape using caramel models.

All other papers

Ziming Wang, Dicong Liang, Junjie Zhao, Chang Liu, Lijing Shao

29 pages, 13 figures

In future gravitational-wave (GW) detections, a large number of overlapping GW signals will appear in the data stream of detectors. When extracting information from one signal, the presence of other signals can cause large parameter estimation biases. Using the Fisher matrix (FM), we develop a bias analysis procedure to investigate how each parameter of other signals affects the inference biases. Taking two-signal overlapping as an example, we show detailedly and quantitatively that the biases essentially originate from the overlapping of the frequency evolution. Furthermore, we find that the behaviors of the correlation coefficients between the parameters of the two signals are similar to the biases. Both of them can be used as characterization of the influence between signals. We also corroborate the bias results of the FM method with full Bayesian analysis. Our results provide powerful guidance for parameter estimation, and the analysis methodology is easy to generalize.

Daniel R. Rybarczyk, Snezana Stanimirovic, Antoine Gusdorf

19 pages, 6 figures. Accepted to ApJ

Interstellar shocks, a key element of stellar feedback processes, shape the structure of the interstellar medium (ISM) and are essential for the chemistry, thermodynamics, and kinematics of interstellar gas. Powerful, high-velocity shocks are driven by stellar winds, young supernova explosions, more evolved supernova remnants, cloud-cloud collisions, and protostellar outflows, whereas the existence and origin of much-lower-velocity shocks ($\lesssim $ 10 km$~$s$^{-1}$) are not understood. Direct observational evidence for interstellar shocks in diffuse and translucent ISM environments have been especially lacking. We present the most sensitive survey to date of SiO -- often considered an unambiguous tracer of interstellar shocks -- in absorption, obtained with the Northern Extended Millimeter Array interferometer. We detect SiO in 5/8 directions probing diffuse and translucent environments without ongoing star formation. Our results demonstrate that SiO formation in the diffuse ISM (i.e., in the absence of significant star formation and stellar feedback) is more widespread and effective than previously reported. The observed SiO linewidths are all $\lesssim$ 4 km$~$s$^{-1}$, excluding high-velocity shocks as a formation mechanism. Yet, the SiO abundances we detect are mostly 1 to 2 orders of magnitude higher than typically assumed in quiescent environments and are often accompanied with other molecular transitions whose column densities cannot be explained with UV-dominated chemical models. Our results challenge the traditional view of SiO production via stellar-feedback sources and emphasize the need for observational constraints on the distribution of Si in the gas phase and grain mantles, which are crucial for understanding the physics of grain processing and diffuse interstellar chemistry.

Xuejian Shen, Josh Borrow, Mark Vogelsberger, Enrico Garaldi, Aaron Smith, Rahul Kannan, Sandro Tacchella, Jesús Zavala, Lars Hernquist, Jessica Y.-C. Yeh, Chunyuan Zheng

To be submitted to MNRAS 22 pages, 12 Figures

Using high-resolution cosmological radiation-hydrodynamic (RHD) simulations (THESAN-HR), we explore the impact of alternative dark matter (altDM) models on galaxies during the Epoch of Reionization. The simulations adopt the IllustrisTNG galaxy formation model. We focus on altDM models that exhibit small-scale suppression of the matter power spectrum, namely warm dark matter (WDM), fuzzy dark matter (FDM), and interacting dark matter (IDM) with strong dark acoustic oscillations (sDAO). In altDM scenarios, both the halo mass functions and the UV luminosity functions at $z\gtrsim 6$ are suppressed at the low-mass/faint end, leading to delayed global star formation and reionization histories. However, strong non-linear effects enable altDM models to "catch up" with cold dark matter (CDM) in terms of star formation and reionization. The specific star formation rates are enhanced in halos below the half-power mass in altDM models. This enhancement coincides with increased gas abundance, reduced gas depletion times, more compact galaxy sizes, and steeper metallicity gradients at the outskirts of the galaxies. These changes in galaxy properties can help disentangle altDM signatures from a range of astrophysical uncertainties. Meanwhile, it is the first time that altDM models have been studied in RHD simulations of galaxy formation. We uncover significant systematic uncertainties in reionization assumptions on the faint-end luminosity function. This underscores the necessity of accurately modeling the small-scale morphology of reionization in making predictions for the low-mass galaxy population. Upcoming James Webb Space Telescope (JWST) imaging surveys of deep, lensed fields hold potential for uncovering the faint, low-mass galaxy population, which could provide constraints on altDM models.

Yossef Zenati, C. Albertus, M. Ángeles Pérez-García, Joseph Silk

comments are welcome

We calculate the neutrino luminosity in an astrophysical scenario where dark matter is captured by a neutron star which eventually implodes to form a low mass black hole. The Trojan horse scenario involves the collapse of a neutron star (NS) due to the accumulation of a critical amount of dark matter (DM) during its lifetime. As a result, a central disk forms out of the ejected material with a finite radial extension, density, temperature, and lepton fraction, producing fainter neutrino luminosities and colder associated spectra than found in a regular core-collapse supernova. The emitted gravitational wave (GW) signal from the imploding NS should be detectable at ultra-high $\gtrsim 0.1$ GHz frequencies.

Anirban Roy, Dariannette Valentín-Martínez, Kailai Wang, Nicholas Battaglia, Alexander van Engelen

19 pages, 10 figures, comments are welcome

Mapping of multiple lines such as the fine-structure emission from [CII] (157.7 $\mu \text{m}$), [OIII] (52 \& 88.4 $\mu \text{m}$), and rotational emission lines from CO are of particular interest for upcoming line intensity mapping (LIM) experiments at millimetre wavelengths, due to their brightness features. Several upcoming experiments aim to cover a broad range of scientific goals, from detecting signatures of the epoch of reionization to the physics of star formation and its role in galaxy evolution. In this paper, we develop a semi-analytic approach to modelling line strengths as functions of the star formation rate (SFR) or infrared (IR) luminosity based on observations of local and high-z galaxies. This package, $\texttt{LIMpy}$ (Line Intensity Mapping in Python), estimates the intensity and power spectra of [CII], [OIII], and CO rotational transition lines up to the $J$-levels (1-0) to (13-12) based both on analytic formalism and on simulations. We develop a relation among halo mass, SFR, and multi-line intensities that permits us to construct a generic formula for the evolution of several line strengths up to $z \sim 10$. We implement a variety of star formation models and multi-line luminosity relations to estimate the astrophysical uncertainties on the intensity power spectrum of these lines. As a demonstration, we predict the signal-to-noise ratio of [CII] detection for an EoR-Spec-like instrument on the Fred Young Submillimeter Telescope (FYST). Furthermore, the ability to use any halo catalogue allows the $\texttt{LIMpy}$ code to be easily integrated into existing simulation pipelines, providing a flexible tool to study intensity mapping in the context of complex galaxy formation physics.

Keith Inight, Boris Gänsicke, Elmé Breedt, Henry Israel, Stuart Littlefair, Christopher Manser, Thomas Marsh, Timothy Mulvany, Anna Pala, John Thorstensen

Submitted to MNRAS. Supplementary information available at this https URL

We present a catalogue of 507 cataclysmic variables (CVs) observed in SDSS I to IV including 70 new classifications collated from multiple archival data sets. This represents the largest sample of CVs with high-quality and homogeneous optical spectroscopy. We have used this sample to derive unbiased space densities and period distributions for the major sub-types of CVs. We also report on some peculiar CVs, period bouncers and also CVs exhibiting large changes in accretion rates. We report 70 new CVs, 59 new periods, 178 unpublished spectra and 262 new or updated classifications. From the SDSS spectroscopy, we also identified 18 systems incorrectly identified as CVs in the literature. We discuss the observed properties of 13 peculiar CVS, and we identify a small set of eight CVs that defy the standard classification scheme. We use this sample to investigate the distribution of different CV sub-types, and we estimate their individual space densities, as well as that of the entire CV population. The SDSS I to IV sample includes 14 period bounce CVs or candidates. We discuss the variability of CVs across the Hertzsprung-Russell diagram, highlighting selection biases of variability-based CV detection. Finally, we searched for, and found eight tertiary companions to the SDSS CVs. We anticipate that this catalogue and the extensive material included in the Supplementary Data will be useful for a range of observational population studies of CVs.

Keith Doore, Erik B. Monson, Rafael T. Eufrasio, Bret D. Lehmer, Kristen Garofali, Antara Basu-Zych

34 pages, 17 figures. Accepted to ApJ

We present an updated version of Lightning, a galaxy spectral energy distribution (SED) fitting code that can model X-ray to submillimeter observations. The models in Lightning include the options to contain contributions from stellar populations, dust attenuation and emission, and active galactic nuclei (AGN). X-ray emission, when utilized, can be modeled as originating from stellar compact binary populations with the option to include emission from AGN. We have also included a variety of algorithms to fit the models to observations and sample parameter posteriors; these include an adaptive Markov-Chain Monte-Carlo (MCMC), affine-invariant MCMC, and Levenberg-Marquardt gradient decent (MPFIT) algorithms. To demonstrate some of the capabilities of Lightning, we present several examples using a variety of observational data. These examples include (1) deriving the spatially resolved stellar properties of the nearby galaxy M81, (2) demonstrating how X-ray emission can provide constrains on the properties of the supermassive black hole of a distant AGN, (3) exploring how to rectify the attenuation effects of inclination on the derived the star formation rate of the edge-on galaxy NGC 4631, (4) comparing the performance of Lightning to similar Bayesian SED fitting codes when deriving physical properties of the star-forming galaxy NGC 628, and (5) comparing the derived X-ray and UV-to-IR AGN properties from Lightning and CIGALE for a distant AGN. Lightning is an open-source application developed in the Interactive Data Language (IDL) and is available at https://github.com/rafaeleufrasio/lightning.

M. Perna, S. Arribas, M. Marshall, F. D'Eugenio, H. Übler, A. Bunker, S. Charlot, S. Carniani, P. Jakobsen, R. Maiolino, B. Rodríguez Del Pino, C. J. Willott, T. Böker, C. Circosta, G. Cresci, M. Curti, B. Husemann, N. Kumari, I. Lamperti, P. G. Pérez-González, J. Scholtz

23 pages, 23 figures; Submitted to A&A; Comments welcome!

LBQS 0302-0019 is a blue quasar (QSO) at z $\sim$ 3.3, hosting powerful outflows, and residing in a complex environment consisting of an obscured AGN candidate, and multiple companions, all within 30 kpc in projection. We use JWST NIRSpec IFS observations to characterise the ionized gas in this complex system. We develop a procedure to correct for the spurious oscillations (or 'wiggles') in NIRSpec single-spaxel spectra, due to the spatial under-sampling of the point spread function. We perform a quasar-host decomposition with the QDeblend3D tools, and use multi-component kinematic decomposition of the optical emission line profiles to infer the physical properties of the emitting gas. The quasar-host decomposition allows us to identify i) a low-velocity component possibly tracing a warm rotating disk, with a dynamical mass Mdyn $\sim 10^{11}$ Msun and a rotation-to-random motion ratio $v_{rot}$/$\sigma_0 \sim 2$; ii) a spatially unresolved ionised outflow, with a velocity of $\sim$ 1000 km/s and an outflow mass rate of $\sim 10^4$ Msun/yr. We also detect eight interacting companion objects close to LBQS 0302-0019. Optical line ratios confirm the presence of a second, obscured AGN at $\sim 20$ kpc of the primary QSO; the dual AGN dominates the ionization state of the gas in the entire NIRSpec field-of-view. This work has unveiled with unprecedented detail the complex environment of this dual AGN, which includes nine interacting companions (five of which were previously unknown), all within 30 kpc of the QSO. Our results support a scenario where mergers can trigger dual AGN, and can be important drivers for rapid early SMBH growth.

Christopher J. Lindsay, J. M. Joel Ong, Sarbani Basu

14 pages, 5 figures, Accepted for publication in ApJ: April 12, 2023

Asteroseismology has been used extensively in recent years to study the interior structure and physical processes of main sequence stars. We consider prospects for using pressure modes (p-modes) near the frequency of maximum oscillation power to probe the structure of the near-core layers of main sequence stars with convective cores by constructing stellar model tracks. Within our mass range of interest, the inner turning point of p modes as determined by the JWKB approximation evolves in two distinct phases during the main sequence, implying a sudden loss of near-core sensitivity during the discontinuous transition between the two phases. However, we also employ non-JWKB asymptotic analysis to derive a contrasting set of expressions for the effects that these structural properties will have on the mode frequencies, which do not encode any such transition. We show analytically that a sufficiently near-core perturbation to the stellar structure results in non-oscillatory, degree-dependent perturbations to the star's oscillation mode frequencies, contrasting with the case of an outer glitch. We also demonstrate numerically that these near-core acoustic glitches exhibit strong angular degree dependence, even at low degree, agreeing with the non-JWKB analysis, rather than the degree-independent oscillations which emerge from JWKB analyses. These properties have important implications for using p-modes to study near-core mixing processes for intermediate-mass stars on the main sequence, as well as for the interpretation of near-center acoustic glitches in other astrophysical configurations, such as red giants.

Paulo E. Stecchini, Marcos P. Diaz, Flavio D'Amico, Francisco Jablonski

9 pages, 8 figures. To be published in MNRAS

In this study we revisit public data on the supersoft X-ray source CAL 83 in the Large Magellanic Cloud. A significant part of our analysis is focused on XMM-Newton X-ray observations, in which updated data reduction procedures and quality assessment were applied. We report on the capability of publicly available hot atmosphere models in describing the source's soft X-ray spectrum. By gathering historical flux measurements in multiple wavelengths and comparing them with the fluxes derived from the X-ray analysis, we find that a $\sim$ 360 kK phenomenological blackbody model describes the spectral energy distribution of CAL 83 fairly well. We also retrieve data from the XMM-Newton UV/optical camera, which is co-alligned with the X-ray instruments and provides strictly simultaneous measurements. These observations demonstrate that the X-ray emission is definitely anti-correlated with emission at longer wavelengths in a time-scale of days to weeks. A closer look at simultaneous X-ray and UV count rates in single light curves reveals that the anti-correlated behaviour is actually present in time scales as short as minutes, suggesting that the origin of variable emission in the system is not unique.

J. A. Guerra, E. Lopez-Rodriguez, D. T. Chuss, N. O. Butterfield, J. T. Schmelz

Submitted to ApJ; Under review

Recent high-resolution 53-$\mu$m polarimetric observations from SOFIA/HAWC+ revealed the inferred plane-of-the-sky magnetic field (B-field) orientation in the Galactic center's Circum-Nuclear Disk (CND). The B-field is mostly aligned with the steamers of ionized material falling onto Sgr A* at large, differential velocities (shear). In such conditions, estimating the B-field strength with the ``classical" Davis-Chandrasekhar-Fermi (DCF) method does not provide accurate results. We derive a ``modified'' DCF method by solving the ideal MHD equations from first principles considering the effects of a large-scale, shear flow on the propagation of a fast magnetosonic wave. In the context of the DCF approximation, both the value of the shear and its Laplacian affect the inferred B-field strength. Using synthetic polarization data from MHD simulations for a medium dominated by shear flows, we find that the ``classical'' DCF determines B-field strengths only within $>50$\% of the true value where the ``modified" DCF results are improved significantly ($\sim$3-22\%). Applying our ``modified'' DCF method to the CND revealed B-field strengths of $\sim1 - 16$ mG in the northern arm, $\sim1 - 13$ mG in the eastern arm, and $\sim3 - 27$ mG in the western arm at spatial scales $\lesssim1$ pc. The balance between turbulent gas energy (kinetic + hydrostatic) and turbulent magnetic energy densities suggest that, along the magnetic-field-flow direction, magnetic effects become less dominant as the shear flow increases and weakens the B-field via magnetic convection. Our results indicate that the transition from magnetically to gravitationally dominated accretion of material onto Sgr A* starts at distances $\sim$ 1 pc.

Everett Schlawin, Ryan Challener, Megan Mansfield, Emily Rauscher, Arthur D. Adams, Jacob Lustig-Yaeger

AJ, accepted, 22 pages

High precision lightcurves combined with eclipse mapping techniques can reveal the horizontal and vertical structure of a planet's thermal emission and the dynamics of hot Jupiters. Someday, they even may reveal the surface maps of rocky planets. However, inverting lightcurves into maps requires an understanding of the planet, star and instrumental trends because they can resemble the gradual flux variations as the planet rotates (ie. partial phase curves). In this work, we simulate lightcurves with baseline trends and assess the impact on planet maps. Baseline trends can be erroneously modeled by incorrect astrophysical planet map features, but there are clues to avoid this pitfall in both the residuals of the lightcurve during eclipse and sharp features at the terminator of the planet. Models that use a Gaussian process or polynomial to account for a baseline trend successfully recover the input map even in the presence of systematics but with worse precision for the m=1 spherical harmonic terms. This is also confirmed with the ThERESA eigencurve method where fewer lightcurve terms can model the planet without correlations between the components. These conclusions help aid the decision on how to schedule observations to improve map precision. If the m=1 components are critical, such as measuring the East/West hotspot shift on a hot Jupiter, better characterization of baseline trends can improve the m=1 terms' precision. For latitudinal North/South information from the remaining mapping terms, it is preferable to obtain high signal-to-noise at ingress/egress with more eclipses.

We investigate the two-halo galactic conformity effect for central galaxies, which is the spatial correlation of the star formation activities for central galaxies to several Mpcs, by studying the dependence of the star formation activities of central galaxies on their large-scale structure in our local Universe using the SDSS data. Here we adopt a novel environment metric using only central galaxies quantified by the distance to the $n$-th nearest central galaxy. This metric measures the environment within an aperture from $\sim$ 1 Mpc to $\gtrsim$ 10 Mpc, with a median value of $\sim$ 4 Mpc. We found that two kinds of conformity effects in our local Universe. The first one is that low-mass central galaxies are more quenched in high-density regions, and we found that this effect mainly comes from low-mass centrals that are close to a more massive halo. A similar trend is also found in the IllustrisTNG simulation, which can be entirely explained by backsplash galaxies. The second conformity effect is that massive central galaxies in low-density regions are more star-forming. This population of galaxies also possesses a higher fraction of spiral morphology and lower central stellar velocity dispersion, suggesting that their low quiescent fraction is due to less-frequent major merger events experienced in the low-density regions, and as a consequence, less-massive bulges and central black holes.

Christopher Sneden, Ann Merchant Boesgaard, John J. Cowan, Ian U. Roederer, Elizabeth A. Den Hartog, James E. Lawler

ApJ, in press

We have derived new detailed abundances of Mg, Ca, and the Fe-group elements Sc through Zn (Z = 21-30) for 37 main sequence turnoff very metal-poor stars ([Fe/H] <= -2.1). We analyzed Keck HIRES optical and near-UV high signal-to-noise spectra originally gathered for a beryllium abundance survey. Using typically about 400 Fe-group lines with accurate laboratory transition probabilities for each star, we have determined accurate LTE metallicities and abundance ratios for neutral and ionized species of the 10 Fe-group elements as well as alpha elements Mg and Ca. We find good neutral/ion abundance agreement for the 6 elements that have detectable transitions of both species in our stars in the 3100-5800 Angstrom range. Earlier reports of correlated Sc-Ti-V relative overabundances are confirmed, and appear to slowly increase with decreasing metallicity. To this element trio we add Zn; it also appears to be increasingly overabundant in the lowest metallicity regimes. Co appears to mimic the behavior of Zn, but issues surrounding its abundance reliability cloud its interpretation.

Andrew B. Pace (CMU), Sergey E. Koposov, Matthew G. Walker, Nelson Caldwell, Mario Mateo, Edward W. Olszewski, Ian U. Roederer, John I. Bailey III, Vasily Belokurov, Kyler Kuehn, Ting S. Li, Daniel B. Zucker

19 pages, 15 figures, submitted to MNRAS, associated data products available at this https URL

We present Magellan/M2FS spectroscopy of four recently discovered Milky Way star clusters (Gran 3, Gran 4, Garro 01, LP 866) and two newly discovered open clusters (Gaia 9, Gaia 10) at low Galactic latitudes. We measure line-of-sight velocities and stellar parameters ([Fe/H], $\log{g}$, $T_{\rm eff}$, [Mg/Fe]) from high resolution spectroscopy centered on the Mg triplet and identify 20-80 members per star cluster. We determine the kinematics and chemical properties of each cluster and measure the systemic proper motion and orbital properties by utilizing Gaia astrometry. We find Gran 3 to be an old, metal-poor (mean metallicity of [Fe/H]=-1.84) globular cluster located in the Galactic bulge on a retrograde orbit. Gran 4 is an old, metal-poor ([Fe/H]}=-1.84) globular cluster with a halo-like orbit that happens to be passing through the Galactic plane. The orbital properties of Gran 4 are consistent with the proposed LMS-1/Wukong and/or Helmi streams merger events. Garro 01 is an old, metal-rich ([Fe/H]=-0.30) globular cluster on a near circular orbit in the outer disk. Gaia 9 and Gaia 10 are among the most distant known open clusters at $R_{GC}\sim 18, 21.2~kpc$ and most metal-poor with [Fe/H]~-0.50,-0.46 for Gaia 9 and Gaia 10, respectively. LP 866 is a nearby, metal-rich open cluster ([Fe/H]$=+0.1$). The discovery and confirmation of multiple star clusters in the Galactic plane shows the power of {\it Gaia} astrometry and the star cluster census remains incomplete.

Federico Spada

MATLAB code will be shared upon reasonable request to the author. Comments are welcome!

I present a re-analysis of the available observational constraints on the trajectory of 'Oumuamua, the first confirmed interstellar object discovered in the solar system. 'Oumuamua passed through the inner solar system on a hyperbolic (i.e., unbound) trajectory. Its discovery occurred after perihelion passage, and near the time of its closest approach to Earth. After being observable for approximately four months, the object became too faint and was lost at a heliocentric distance of around 3 au. Intriguingly, analysis of the trajectory of 'Oumuamua revealed that a dynamical model including only gravitational accelerations does not provide a satisfactory fit of the data, and a non-gravitational term must be included. The detected non-gravitational acceleration is compatible with either solar radiation pressure or recoil due to outgassing. It has, however, proved challenging to reconcile either interpretation with the existing quantitative models of such effects without postulating unusual physical properties for 'Oumuamua (such as extremely low density and/or unusual geometry, non-standard chemistry). My analysis independently confirms the detection of the non-gravitational acceleration. After comparing several possible parametrizations for this effects, I find a strong preference for a radially directed non-gravitational acceleration, pointing away from the Sun, and a moderate preference for a power-law scaling with the heliocentric distance, with an exponent between 1 and 2. These results provide valuable constraints on the physical mechanism behind the effect; a conclusive identification, however, is probably not possible on the basis of dynamical arguments alone.

Frédéric Marin, Eugene Churazov, Ildar Khabibullin, Riccardo Ferrazzoli, Laura Di Gesu, Thibault Barnouin, Alessandro Di Marco, Riccardo Middei, Alexey Vikhlinin, Enrico Costa, Paolo Soffitta, Fabio Muleri, Rashid Sunyaev, William Forman, Ralph Kraft, Stefano Bianchi, Immacolata Donnarumma, Pierre-Olivier Petrucci, Teruaki Enoto, Iván Agudo, Lucio A. Antonelli, Matteo Bachetti, Luca Baldini, Wayne H. Baumgartner, Ronaldo Bellazzini, Stephen D. Bongiorno, Raffaella Bonino, Alessandro Brez, Niccolò Bucciantini, Fiamma Capitanio, Simone Castellano, Elisabetta Cavazzuti, Chien-Ting Chen, Stefano Ciprini, Alessandra De Rosa, Ettore Del Monte, Niccolò Di Lalla, Victor Doroshenko, Michal Dovciak, Steven R. Ehlert, Yuri Evangelista, Sergio Fabiani, Javier A. Garcia, Shuichi Gunji, et al. (60 additional authors not shown)

24 pages, 5 figures, 4 tables, author's version of the paper accepted for publication in Nature

The center of the Milky Way Galaxy hosts a $\sim$4 million solar mass black hole (Sgr A$^*$) that is currently very quiescent with a luminosity many orders of magnitude below those of active galactic nuclei. Reflection of X-rays from Sgr A$^*$ by dense gas in the Galactic Center region offers a means to study its past flaring activity on times scales of hundreds and thousands of years. The shape of the X-ray continuum and the strong fluorescent iron line observed from giant molecular clouds in the vicinity of Sgr A$^*$ are consistent with the reflection scenario. If this interpretation is correct, the reflected continuum emission should be polarized. Here we report observations of polarized X-ray emission in the direction of the Galactic center molecular clouds using the Imaging X-ray Polarimetry Explorer (IXPE). We measure a polarization degree of 31\% $\pm$ 11\%, and a polarization angle of $-$48$^\circ$ $\pm$ 11$^\circ$. The polarization angle is consistent with Sgr A$^*$ being the primary source of the emission, while the polarization degree implies that some 200 years ago the X-ray luminosity of Sgr A$^*$ was briefly comparable to a Seyfert galaxy.

Meng Zhang, Andrea Ferrara, Bin Yue

12 pages, 8 figures

ALMA observations have detected extended ($\simeq 10$ kpc) [C II] halos around high-redshift ($z \gtrsim 5$) star-forming galaxies. If such extended structures are common, they may have an impact on the line intensity mapping (LIM) signal. We compute the LIM power spectrum including both the central galaxy and the [C II] halo, and study the detectability of such signal in an ALMA LIM survey. We model the central galaxy and the [C II] halo brightness with a S\'ersic+exponential profile. The model has two free parameters: the effective radius ratio $f_{R_e}$, and the central surface brightness ratio, $f_{\Sigma}$, between the two components. [C II] halos can significantly boost the LIM power spectrum signal. For example, for relatively compact [C II] halos ($f_\Sigma=0.4$, $f_{R_{\rm e}}=2.0$), the signal is boosted by $\simeq 20$ times; for more extended and diffuse halos ($f_\Sigma=0.1, f_{R_{\rm e}}=6.0$), the signal is boosted by $\simeq 100$ times. For the ALMA ASPECS survey (resolution $\theta_{\rm beam} = 1.13''$, survey area $\Omega_{\rm survey}=2.9\,\rm arcmin^{2}$), the [C II] power spectrum is detectable only if the deL14d [C II] - SFR relation holds. However, with an optimized survey ($\theta_{\rm beam} = 0.232''$, $\Omega_{\rm survey}=2.0\,\rm deg^{2}$), the power spectrum is detectable for almost all the [C II] - SFR relations considered in this paper. Such a survey can constrain $f_\Sigma$ ($f_{R_{\rm e}}$) with a relative uncertainty of $60\%$ ($20\%$). A successful LIM experiment will provide unique constraints on the nature, origin, and frequency of extended [C II] halos, and the [C II] - SFR relation at early times.

Roger Dufresne, Giulio Del Zanna, Helen Mason

Accepted by MNRAS, 16 pages, 3 figures

The use of the coronal approximation to model line emission from the solar transition region has led to discrepancies with observations over many years, particularly for Li- and Na-like ions. Studies have shown that a number of atomic processes are required to improve the modelling for this region, including the effects of high densities, solar radiation and charge transfer on ion formation. Other non-equilibrium processes, such as time dependent ionisation and radiative transfer, are also expected to play a role. A set of models which include the three relevant atomic processes listed above in ionisation equilibrium has recently been built. These new results cover the main elements observed in the transition region. To assess the effectiveness of the results, the present work predicts spectral line intensities using differential emission measure modelling. Although limited in some respects, this differential emission measure modelling does give a good indication of the impact of the new atomic calculations. The results are compared to predictions of the coronal approximation and to observations of the average, quiet Sun from published literature. Significant improvements are seen for the line emission from Li- and Na-like ions, inter-combination lines and many other lines. From this study, an assessment is made of how far down into the solar atmosphere the coronal approximation can be applied, and the range over which the new atomic models are valid.

M Atemkeng, S Perkins, E Seck, S Makhathini, O Smirnov, L Bester, B Hugo

This work proposes to reduce visibility data volume using a baseline-dependent lossy compression technique that preserves smearing at the edges of the field-of-view. We exploit the relation of the rank of a matrix and the fact that a low-rank approximation can describe the raw visibility data as a sum of basic components where each basic component corresponds to a specific Fourier component of the sky distribution. As such, the entire visibility data is represented as a collection of data matrices from baselines, instead of a single tensor. The proposed methods are formulated as follows: provided a large dataset of the entire visibility data; the first algorithm, named $simple~SVD$ projects the data into a regular sampling space of rank$-r$ data matrices. In this space, the data for all the baselines has the same rank, which makes the compression factor equal across all baselines. The second algorithm, named $BDSVD$ projects the data into an irregular sampling space of rank$-r_{pq}$ data matrices. The subscript $pq$ indicates that the rank of the data matrix varies across baselines $pq$, which makes the compression factor baseline-dependent. MeerKAT and the European Very Long Baseline Interferometry Network are used as reference telescopes to evaluate and compare the performance of the proposed methods against traditional methods, such as traditional averaging and baseline-dependent averaging (BDA). For the same spatial resolution threshold, both $simple~SVD$ and $BDSVD$ show effective compression by two-orders of magnitude higher than traditional averaging and BDA. At the same space-saving rate, there is no decrease in spatial resolution and there is a reduction in the noise variance in the data which improves the S/N to over $1.5$ dB at the edges of the field-of-view.

H$\alpha$ heliographs are imaging instruments designed to produce monochromatic images of the solar chromosphere at fast cadence (60 s or less). They are designed to monitor efficiently dynamic phenomena of solar activity, such as flares or material ejections. Meudon and Haute Provence observatories started systematic observations in the frame of the International Geophysical Year (1957) with Lyot filters. This technology evolved several times until 1985 with tunable filters allowing to observe alternatively the line wings and core (variable wavelength). More than 6 million images were produced during 50 years, mostly on 35 mm films (catalogs are available on-line). We present in this paper the optical characteristics and the capabilities of the successive versions of the H$\alpha$ heliographs in operation between 1954 and 2004, and describe briefly the new heliograph (MeteoSpace) which will be commissioned in 2023 at Calern observatory.

Clàudia Soriano-Guerrero, Daniele Viganò, Rosalba Perna, Taner Akgün, Carlos Palenzuela

submitted for publication in MNRAS

While magnetism in exoplanets remains largely unknown, Hot Jupiters have been considered as natural candidates to harbour intense magnetic fields, both due to their large masses and their high energy budgets coming from irradiation as a consequence of their vicinity to their host stars. In this work we perform MHD simulations of a narrow day-side atmospheric column of ultra-hot Jupiters, suitable for very high local temperatures (T > 3000 K). Since the conductivity in this regime is very high, the dominant effect is winding due to the intense zonal winds. By including a forcing that mimics the wind profiles obtained in global circulation models, the shear layer induces a strong toroidal magnetic field (locally reaching hundreds of gauss), supported by meridional currents. Such fields and the sustaining currents don$'$t depend on the internally generated field, but are all confined in the thin (less than a scale-height) shear layer around 1 bar. Additionally, we add random perturbations that induce turbulent motions, which lead to further (but much smaller) magnetic field generation to a broader range of depths. These results allow an evaluation of the currents induced by the atmospheric dynamo. Although here we use ideal MHD and the only resistivity comes from the numerical scheme, we estimate a-posteriori the amount of Ohmic heat deposited in the outer layers, which could be employed in evolutionary models for Hot Jupiters' inflated radii.

Yuta Yagi, Ryohei Konno, Tasuku Hayash, Keita Tanaka, Noriko Y. Yamasaki, Kazuhisa Mitsuda, Rumi Sato, Mikiko Saito, Takayuki Homma, Yoshiki Nishida, Shohei Mori, Naoko Iyomoto, Toru Hara

10 pages, 6 figures, published in Journal of Low Temperature Physics on 4 February 2023

A 57Fe nucleus in the solar core could emit a 14.4-keV monochromatic axion through the M1 transition if a hypothetical elementary particle, axion, exists to solve the strong CP problem. Transition edge sensor (TES) X-ray microcalorimeters can detect such axions very efficiently if they are again converted into photons by a 57Fe absorber. We have designed and produced a dedicated TES array with 57Fe absorbers for the solar axion search. The iron absorber is set next to the TES, keeping a certain distance to reduce the iron-magnetization effect on the spectroscopic performance. A gold thermal transfer strap connects them. A sample pixel irradiated from a 55Fe source detected 698 pulses. In contrast to thermal simulations, we consider that the pulses include either events produced in an iron absorber or gold strap at a fraction dependent on the absorption rate of each material. Furthermore, photons deposited on the iron absorber are detected through the strap as intended. The identification of all events still needs to be completed. However, we successfully operated the TES with the unique design under iron magnetization for the first time.

Avery Meiksin

Published in the Research Notes of the American Astronomical Society; 3 pages, 1 figure

The scattering of Ly$\alpha$ photons from the first radiating sources in the Universe plays a pivotal role in 21-cm radio detections of Cosmic Dawn and the Epoch of Reionization through the Wouthuysen-Field effect. New data from JWST show the Ly$\alpha$ photon scattering rate exceeds that required to decouple the intergalactic hydrogen spin temperature from that of the Cosmic Microwave Background up to $z\sim14$ and render the neutral hydrogen visible over the main redshift range expected for the Epoch of Reionization.

V.P. Grinin, L.V. Tambovtseva, O.Yu. Barsunova, D.N. Shakhovskoy

7 pages, 4 figures, accepted by Astrophysics

The star CQ Tau belongs to the family of UX Ori type stars. It has very complex photometric behavior and complex structure of the circumstellar environment. In our paper we constructed the historical 125 years light curve of this star basing on the published photometric observations. It follows that besides a random component characteristic of UX Ori type stars, the large amplitude periodic component with the 10 year period is also present. Its existence was suspected earlier in [11]. New observations confirm its reality. It points to an existence of the second component close to the star. The density waves and matter flows caused by the companion motion lead to periodic changes in the circumstellar extinction and brightness of the star. This result is discussed in context of the recent observations of CQ Tau with high angular resolution.

Martin M. Roth

9 pages, 10 figures. Presented at SDW2022, accepted for publication in Special Issue of Astronomische Nachrichten

Image sensors, most notably the Charge Coupled Device (CCD), have revolutionized observational astronomy as perhaps the most important innovation after photography. Since the 50th anniversary of the invention of the CCD has passed in 2019, it is time to review the development of detectors for the visible wavelength range, starting with the discovery of the photoelectric effect and first experiments to utilize it for the photometry of stars at Sternwarte Babelsberg in 1913, over the invention of the CCD, its development at the Jet Propulsion Laboratory, to the high performance CCD and CMOS imagers that are available off-the-shelf today.

Saniya Khan, Andrea Miglio, Emma Willett, Benoît Mosser, Yvonne P. Elsworth, Richard I. Anderson, Leo Girardi, Kévin Belkacem, Anthony G. A. Brown, Tristan Cantat-Gaudin, Luca Casagrande, Gisella Clementini, Antonella Vallenari

11 pages, 8 figures, Accepted for publication in A&A

Gaia EDR3 has provided unprecedented data that generate a lot of interest in the astrophysical community, despite the fact that systematics affect the reported parallaxes at the level of ~ 10 muas. Independent distance measurements are available from asteroseismology of red-giant stars with measurable parallaxes, whose magnitude and colour ranges more closely reflect those of other stars of interest. In this paper, we determine distances to nearly 12,500 red-giant branch and red clump stars observed by Kepler, K2, and TESS. This is done via a grid-based modelling method, where global asteroseismic observables, constraints on the photospheric chemical composition, and on the unreddened photometry are used as observational inputs. This large catalogue of asteroseismic distances allows us to provide a first comparison with Gaia EDR3 parallaxes. Offset values estimated with asteroseismology show no clear trend with ecliptic latitude or magnitude, and the trend whereby they increase (in absolute terms) as we move towards redder colours is dominated by the brightest stars. The correction model proposed by Lindegren et al. (2021) is not suitable for all the fields considered in this study. We find a good agreement between asteroseismic results and model predictions of the red clump magnitude. We discuss possible trends with the Gaia scan law statistics, and show that two magnitude regimes exist where either asteroseismology or Gaia provides the best precision in parallax.

Hao-Jui Kuan, Arthur G. Suvorov, Kostas D. Kokkotas

11 pages, 6 figures, 2 tables, with an appendix containing 1 figure

Gamma-ray bursts resulting from binary neutron-star mergers are sometimes preceded by precursor flares. These harbingers may be ignited by quasi-normal modes, excited by orbital resonances, shattering the stellar crust of one of the inspiralling stars up to $\gtrsim10$ seconds before coalescence. In the rare case that a system displays two precursors, successive overtones of either interface- or $g$-modes may be responsible for the overstrainings. Since the free-mode frequencies of these overtones have an almost constant ratio, and the inertial-frame frequencies for rotating stars are shifted relative to static ones, the spin frequency of the flaring component can be constrained as a function of the equation of state, the binary mass ratio, the mode quantum numbers, and the spin-orbit misalignment angle. As a demonstration of the method, we find that the precursors of GRB090510 hint at a spin frequency range of $2 \lesssim \nu_{\star}/\text{Hz} \lesssim 20$ for the shattering star if we allow for an arbitrary misalignment angle, assuming $\ell=2$ $g$-modes account for the events.

Micah Bowles, Hongming Tang, Eleni Vardoulaki, Emma L. Alexander, Yan Luo, Lawrence Rudnick, Mike Walmsley, Fiona Porter, Anna M. M. Scaife, Inigo Val Slijepcevic, Elizabeth A. K. Adams, Alexander Drabent, Thomas Dugdale, Gülay Gürkan, Andrew M. Hopkins, Eric F. Jimenez-Andrade, Denis A. Leahy, Ray P. Norris, Syed Faisal ur Rahman, Xichang Ouyang, Gary Segal, Stanislav S. Shabala, O. Ivy Wong

17 pages, 11 Figures, Accepted at MNRAS

We present a novel natural language processing (NLP) approach to deriving plain English descriptors for science cases otherwise restricted by obfuscating technical terminology. We address the limitations of common radio galaxy morphology classifications by applying this approach. We experimentally derive a set of semantic tags for the Radio Galaxy Zoo EMU (Evolutionary Map of the Universe) project and the wider astronomical community. We collect 8,486 plain English annotations of radio galaxy morphology, from which we derive a taxonomy of tags. The tags are plain English. The result is an extensible framework which is more flexible, more easily communicated, and more sensitive to rare feature combinations which are indescribable using the current framework of radio astronomy classifications.

Christian Glaser (for the ARIANNA collaboration)

Proceedings of the 9th ARENA workshop 2022

The ARIANNA in-ice radio detector explores the detection of UHE neutrinos with shallow detector stations on the Ross Ice Shelf and the South Pole. Here, we present recent results that lay the foundation for future large-scale experiments. We show a limit on the UHE neutrino flux derived from ARIANNA data, measurements of the more abundant air showers, results from in-situ measurement campaigns, a study of a potential background from internal reflection layers, and give an outlook of future detector improvements.

We investigate the effect of conductive heating of the gas surrounding a geometrically thick accretion disk on the growth of a black hole at high redshift. If a black hole is accreting the surrounding gas at a super-Eddington rate, the X-ray radiation from the vicinity of the black hole would be highly anisotropic due to the self-shielding of a geometrically thick accretion disk, and then the radiative feedback on the surrounding medium would be suppressed in the equatorial region, within which super-Eddington accretion can continue. However, if this region is sufficiently heated via thermal conduction from the adjacent region that is not shielded and heated by the X-ray irradiation, the surrounding gas becomes isotropically hot and the Bondi accretion rate would be suppressed and become sub-Eddington. We evaluate the condition under which such isotropic heating is realized, and derive new criteria required for super-Eddington accretion.

Kai Wang, Yangyao Chen, Qingyang Li, Xiaohu Yang

14 pages, 7 + 5 figures, MNRAS accepted

The star formation and quenching of central galaxies are regulated by the assembly histories of their host halos. In this work, we use the central stellar mass to halo mass ratio as a proxy of halo formation time, and we devise three different models, from the physical hydrodynamical simulation to the empirical statistical model, to demonstrate its robustness. With this proxy, we inferred the dependence of the central galaxy properties on the formation time of their host halos using the SDSS main galaxy sample, where central galaxies are identified with the halo-based group finder. We found that central galaxies living in late-formed halos have higher quiescent fractions and lower spiral fractions than their early-formed counterparts by $\lesssim$ 8%. Finally, we demonstrate that the group finding algorithm has a negligible impact on our results.

Anna L. H. Hughes, Irene Gonzalez-Hernandez, Sean G. McManus, Kiran Jain, Sushanta C. Tripathy

36 pages, 32 figures

In helioseismic studies, an observational parameter of primary concern is the P-angle, the angle along which lies the solar axis of rotation for a given image. For the six observing sites employed by The Global Oscillation Network Group (GONG), this angle acts additionally as a marker of relative image orientation, allowing concurrent images to be precisely aligned and merged to provide the highest possible quality data. In this report, we present and investigate two methods of determining the P-angle via the rotational signature embedded in solar Dopplergram images by examining the large-scale structure of the observed velocity field. As with other studies, we find that the Dopplergram produces a time-varying 'P-angle' signature according to the presentation of various physical phenomena across the solar surface, but with the potential for sub-degree identification of the axis of rotation. However, close agreement between separate P-angle-finding techniques also reveals current limitations to P-angle determination that are imposed by the calibration state of the GONG-site Dopplergrams, leaving these P-angle-finding methods for GONG with errors on the scale of less than a degree between two site.

Guoqing Zhen, Guoliang Lv, Helei Liu, Akira Dohi, Bobuya Nishimura, Chunhua Zhu, Liyu Song, Weiyang Wang, Renxin Xu

14 pages, 10 figures, accepted for publication in ApJ

Type I X-ray bursts are rapidly brightening phenomena triggered by thermonuclear burning on accreting layer of a neutron star (NS). The light curves represent the physical properties of NSs and the nuclear reactions on the proton-rich nuclei. The numerical treatments of the accreting NS and physics of the NS interior are not established, which shows uncertainty in modelling for observed X-ray light curves. In this study, we investigate theoretical X-ray-burst models, compared with burst light curves with GS~1826-24 observations. We focus on the impacts of the NS mass, the NS radius, and base-heating on the NS surface using the MESA code. We find a monotonic correlation between the NS mass and the parameters of the light curve. The higher the mass, the longer the recurrence time and the greater the peak luminosity. While the larger the radius, the longer the recurrence time, the peak luminosity remains nearly constant. In the case of increasing base heating, both the recurrence time and peak luminosity decrease. We also examine the above results using with a different numerical code, HERES, based on general relativity and consider the central NS. We find that the burst rate, burst energy and burst strength are almost same in two X-ray burst codes by adjusting the base-heat parameter in MESA (the relative errors $\lesssim5\%$), while the duration time and the rise time are significantly different between (the relative error is possibly $\sim50\%$). The peak luminosity and the e-folding time are ragged between two codes for different accretion rates.

T. Gessey-Jones, A. Fialkov, E. de Lera Acedo, W. J. Handley, R. Barkana

21 pages, 10 figures. Submitted to MNRAS

Cosmic rays generated by supernovae carry away a significant portion of the lifetime energy emission of their parent star, making them a plausible mechanism for heating the early universe intergalactic medium (IGM). Following a review of the existing literature on cosmic ray heating, we develop a flexible model of this heating mechanism for use in semi-numerical 21-cm signal simulations and conduct the first investigations of the signatures it imprints on the 21-cm power spectrum and tomographic maps. We find that cosmic ray heating of the IGM is short-ranged, leading to heating clustered around star-forming sites, and a sharp contrast between heated regions of 21-cm emission and unheated regions of absorption. This contrast results in greater small-scale power for cosmic ray heated scenarios compared to what is found for X-ray heating, thus suggesting a way to test the nature of IGM heating with future 21-cm observations. Finally, we find an unexpectedly rich thermal history in models where cosmic rays can only escape efficiently from low-mass halos, such as in scenarios where these energetic particles originate from population III star supernovae remnants. The interplay of heating and the Lyman-Werner feedback in these models can produce a local peak in the IGM kinetic temperature and, for a limited parameter range, a flattened absorption trough in the global 21-cm signal.

N. Hatamkhani, R. C. Kraan-Korteweg, S. L. Blyth, K. Said, A. Elagali

18 pages, 10 Figures, 13 Tables. Accepted for publication in MNRAS

We present six deep Near-InfraRed (JHK_s) photometric catalogues of galaxies identified in six cluster candidates (VC02, VC04, VC05, VC08, VC10, VC11) within the Vela Supercluster (VSCL) as part of our efforts to learn more about this large supercluster which extends across the zone of avoidance (l=272.5 \pm 20 deg, b= \pm 10 deg, at cz~ 18000 km/s). The observations were conducted with the InfraRed Survey Facility (IRSF), a 1.4m telescope situated at the South African Astronomical Observatory (SAAO) in Sutherland. The images in each cluster cover ~ 80% of their respective Abell radii. We identified a total number of 1715 galaxies distributed over the six cluster candidates, of which only ~ 15% were previously known. We study the structures and richnesses of the six clusters out to the cluster-centric completeness radius of r_c<1.5 Mpc and magnitude completeness limit of K_s^0<15.5 mag, using their iso-density contour maps and radial density profiles. The analysis shows VC04 to be the richest of the six. It is a massive cluster comparable to the Coma and Norma clusters, although its velocity dispersion, sigma=455 km/s, seems rather low for a rich cluster. VC02 and VC05 are found to be relatively rich clusters while VC08 is rather poor. Also, VC05 has the highest central number density among the six. VC11 is an intermediate cluster that contains two major subclusters while VC10 has a filament-like structure and is likely not to be a cluster after all.

Matthew E. Caplan

3 pages, 1 figure, submitted to RNAAS

The diffusion coefficients of neutron rich nuclei in crystallizing white dwarf (WD) stars are essential microphysics input for modeling the evolution of the composition profile. Recently, molecular dynamics simulations have been used to compute diffusion coefficients for realistic mixtures of C-O and O-Ne WDs with many trace nuclides that could be important sedimentary heat sources such as $^{22}$Ne, $^{23}$Na, $^{25}$Mg, and $^{27}$Mg. In this brief note, I repeat these simulations but now include $^{56}$Fe. I find that for the large charge ratios involved in these mixtures the empirical law developed in our earlier work tends to under-predict diffusion coefficients in the moderately coupled regime by 30 to 40 percent. As this formalism is presently implemented in the stellar evolution code MESA, it is important for authors studying mixtures containing heavy nuclides like $^{56}$Fe to be aware of these systematics. However, the impact on astrophysics is expected to be small.

Caroline Gieser, Henrik Beuther, Dmitry Semenov, Aida Ahmadi, Thomas Henning, Molly Wells

21 pages, 10 figures, submitted to A&A

In this study, we investigate how physical properties, such as the density and temperature profiles, evolve on core scales through the evolutionary sequence during high-mass star formation ranging from protostars in cold infrared dark clouds to evolved UCHII regions. We observed 11 high-mass star-forming regions with ALMA at 3 mm wavelengths. Based on the 3 mm continuum morphology and recombination line emission, tracing locations with free-free (ff) emission, the fragmented cores analyzed in this study are classified into either dust or dust+ff cores. In addition, we resolve three cometary UCHII regions with extended 3 mm emission that is dominated by free-free emission. The temperature structure and radial profiles (T~r^-q ) are determined by modeling molecular emission of CH3CN and CH313CN with XCLASS and by using the HCN-to- HNC intensity ratio as probes for the gas kinetic temperature. The density profiles (n~r^-p ) are estimated from the 3 mm continuum visibility profiles. The masses M and H2 column densities N(H2) are then calculated from the 3 mm dust continuum emission. Results. We find a large spread in mass and peak H2 column density in the detected sources ranging from 0.1-150 Msun and 10^23 - 10^26 cm-2 , respectively. Including the results of the CORE and CORE-extension studies (Gieser et al. 2021, 2022) to increase the sample size, we find evolutionary trends on core scales for the temperature power-~law index q increasing from 0.1 to 0.7 from infrared dark clouds to UCHII regions, while for the the density power-law index p on core scales, we do not find strong evidence for an evolutionary trend. However, we find that on the larger clump scales throughout these evolutionary phases the density profile flattens from p = 2.2 to p = 1.2. (abridged)

Rakhi R, Geethika Santhosh, Prajwel Joseph, Koshy George, Smitha Subramanian, Indulekha Kavila, J. Postma, Pierre-Alain Duc, Patrick Côté, Luca Cortese, S. K. Ghosh, Annapurni Subramaniam, Shyam Tandon, John Hutchings, P Samuel Wesley, Aditya Bharadwaj, Neeran Niroula

Accepted for publication in MNRAS. Appendix table can be found in published version

NGC 5291, an early-type galaxy surrounded by a giant HI ring, is believed to be formed from collision with another galaxy. Several star forming complexes and tidal dwarf galaxies are distributed along the collisional ring which are sites of star formation in environments where extreme dynamical effects are involved. Dynamical effects can affect the star formation properties and the spatial distribution of star forming complexes along the tidal features. To study and quantify the star formation activity in the main body and in the ring structure of the NGC 5291 system, we use high spatial resolution FUV and NUV imaging observations from the Ultraviolet Imaging Telescope onboard AstroSat. A total of 57 star-forming knots are identified to be part of this interacting system out of which 12 are new detections (star forming complexes that lie inside the HI contour) compared to the previous measurements from lower resolution UV imaging. We estimate the attenuation in UV for each of the resolved star-forming knots using the UV spectral slope $\beta$, derived from the FUV-NUV colour. Using the extinction corrected UV fluxes, we derive the star formation rate of the resolved star forming complexes. The extinction corrected total star formation rate of this system is estimated as 1.75 $\pm$ 0.04 $M_{\odot}/yr$. The comparison with dwarf galaxy populations (BCD, Sm and dIm galaxies) in the nearby Universe shows that many of the knots in the NGC 5291 system have SFR values comparable to the SFR of BCD galaxies.

Gisela N. Ortiz-Leon, Sergio A. Dzib, Laurent Loinard, Yan Gong, Thushara Pillai Adele Plunkett

Accepted to A&A Letters

In this Letter, we report Very Long Baseline Array observations of 22 GHz water masers toward the protostar CARMA-6, located at the center of the Serpens South young cluster. From the astrometric fits to maser spots, we derive a distance of 440.7+/-3.5 pc for the protostar (1% error). This represents the best direct distance determination obtained so far for an object this young and deeply embedded in this highly obscured region. Taking into account depth effects, we obtain a distance to the cluster of 440.7+/-4.6 pc. Stars visible in the optical that have astrometric solutions in the Gaia Data Release 3 are, on the other hand, all located in the periphery of the cluster. Their mean distance of 437 (+51, -41) pc is consistent within 1-sigma with the value derived from maser astrometry. As the maser source is just at the center of Serpens South, we finally solve the ambiguity of the distance to this region that has prevailed over the years.

Chris Hamilton (IAS), Tobias Heinemann (NBI)

5 pages, 2 figures

The linear response of a stellar system's gravitational potential to a perturbing mass comprises two distinct contributions. Most famously, the system will respond by forming a polarization `wake' around the perturber. At the same time, the perturber may also excite one or more `normal modes', i.e. coherent oscillations of the entire stellar system which are either stable or unstable depending on the system parameters. The amplitude of the first (wake) contribution is known to diverge as a system approaches marginal stability. In this paper we consider the linear response of a homogeneous stellar system to a point mass moving on a straight line orbit. We prove analytically that the divergence of the wake response is in fact cancelled by a corresponding divergence in the normal mode response, rendering the total response finite. We demonstrate this cancellation explicitly for a box of stars with Maxwellian velocity distribution. Our results imply that polarization wakes may be much less efficient drivers of secular evolution than previously thought. More generally, any prior calculation that accounted for wakes but ignored modes may need to be revised.

Yuewei Wen, Nhat-Minh Nguyen, Dragan Huterer

23 pages, 6 figures

We propose and numerically validate a new fitting formula that is sufficiently accurate to model the growth of structure in Horndeski theories of modified gravity for upcoming Stage IV and V large-scale structure surveys. Based on an analysis of more than 18,000 Horndeski models and adopting the popular parameterization of the growth rate $f(z) = \Omega_{M}(z)^{\gamma}$, we generalize the constant growth index $\gamma$ to a two-parameter redshift-dependent quantity, $\gamma(z)$, that more accurately fits these models. We demonstrate that the functional form $\gamma(z)=\gamma_0+\gamma_1z^2 / (1+z)$ improves the median $\chi^2$ of the fit to viable Horndeski models by a factor of $\sim40$ relative to that of a constant $\gamma$, and is sufficient to obtain unbiased results even for precise measurements expected in Stage IV and V surveys. Finally, we constrain the parameters of the new fitting formula using current cosmological data.

We explore the effect magnetic fields have on self-gravitating accretion disks around spinning black holes via numerical evolutions in full dynamical magnetohydrodynamic spacetimes. The configurations we study are unstable to the Papaloizou-Pringle Instability (PPI). PPI-saturated accretion tori have been shown to produce gravitational waves, detectable to cosmological distances by third-generation gravitational wave (GW) observatories. While the PPI operates strongly for purely hydrodynamic disks, the situation can be different for disks hosting initially small magnetic fields. Evolutions of disks without self-gravity in fixed BH spacetimes have shown that small seed fields can initiate the rapid growth of the magneto-rotational instability (MRI), which then strongly suppresses the PPI. Since realistic astrophysical disks are expected to be magnetized, PPI-generated GW signals may be suppressed as well. However, it is unclear what happens when the disk self-gravity is restored. Here, we study the impact of magnetic fields on the PPI-saturated state of a self-gravitating accretion disk around a spinning BH ($\chi = 0.7$) aligned with the disk angular momentum, as well as one around a non-spinning BH. We find the MRI is effective at reducing the amplitude of PPI modes and their associated GWs, but the systems still generate GWs. Estimating the detectability of these systems accross a wide range of masses, we show that magnetic fields reduce the maximum detection distance by Cosmic Explorer from 300Mpc (in the pure hydrodynamic case) to 45Mpc for a $10 M_{\odot}$ system, by LISA from 11500Mpc to 2700Mpc for a $2 \times 10^{5} M_{\odot}$ system, and by DECIGO from $z \approx 5$ down to $z \approx 2$ for a $1000 M_{\odot}$ system.

Jaemyoung Lee, J. Richard Bond, Pavel Motloch, Alexander van Engelen, George Stein

16 pages, 18 figures

Gravitational lensing deflects the paths of photons, altering the statistics of cosmic backgrounds and distorting their information content. We take the Cosmic Infrared Background (CIB), which provides a wealth of information about galaxy formation and evolution, as an example to probe the effect of gravitational lensing on non-Gaussian statistics. Using the Websky simulations, we first quantify the non-Gaussianity of the CIB, revealing additional detail on top of its well-measured power spectrum. To achieve this, we use needlet-like multipole-band-filters to calculate the variance and higher-point correlations. We show the 3-point and 4-point spectra, and compare our calculated bispectra to Planck values. We then lens the CIB, shell-by-shell with corresponding convergence maps, to capture the broad redshift extent of both the CIB and its lensing convergence. Using our simulations, we show that the lensed CIB power spectrum and bispectrum agree with observations: the lensing of the CIB changes the 3-point and 4-point functions by a few tens of percent at large scales, unlike with the power spectrum, which changes by less than two percent. We expand our analyses to encompass the full intensity probability distribution functions (PDFs) involving all n-point correlations as a function of scale. In particular we use the relative entropy between lensed and unlensed PDFs to create a spectrum of templates that can allow estimation of lensing. The underlying CIB model has uncertainties, in particular missing the important role of star-bursting, which has a larger effect on higher point correlations than on the variance. We test this by adding a stochastic log-normal term to the intensity distributions. The novel aspects of our filtering and lensing pipeline should prove useful for not just CIB applications, but for any radiant background, including line intensity maps.

Jose Arita-Escalante, Tyler Parsotan, S. Bradley Cenko

12 pages, 10 figures, submitted to ApJ, for calculations, see: this https URL

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydrodynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

We report tidal-induced latency variations on a transpacific subsea cable. Week-long recordings with a precision phase meter suggest length changes in the sub-meter range caused by the Poisson effect. The described method adds to the toolbox for the new field >>optical oceanic seismology<<.

Incredible progress on the theoretical uncertainty of the spatial correlations of the stochastic gravitational wave (GW) background were recently made. However, it remains to realize the impact of this theoretical uncertainty on PTA cross correlations analysis. This paper pushes forward in this direction, as a proof--of--principle: showing the potential role that theoretical uncertainty has on unburying the stochastic GW background signal in noisy PTA cross correlation measurements. We consider both a mock data set and the noise--marginalized 12.5 years NANOGrav spatial correlation measurements, and find optimistic conclusions regardless of the physical content of the GW background and the nature of the noise in the data. Very briefly, we show through various cases a modest, but profound result that looking out for a stochastic signal is better when two of its moments are utilized. Or, in terms of GWs, we show that the theoretical uncertainty can play a substantial role in the hunt for the stochastic GW background.

The parity violating model based on teleparallel gravity is a competitive scheme for parity violating gravity, which has been preliminary studied in the literature. To further investigate the parity violating model in teleparallel gravity, in this paper, we construct all independent parity-odd terms that are quadratic in torsion tensor and coupled to a scalar field in a way without higher-order derivatives. Using these parity-odd terms, we formulate a general parity violating scalar-tensor model in teleparallel gravity and obtain its equations of motion. To explore potentially viable models within the general model, we investigate the cosmological application of a submodel of the general model in which terms above the second power of torsion are eliminated. We focus on analyzing cosmological perturbations and identify the conditions that preserve the parity violating signal of gravitational waves at linear order while avoiding the ghost instability.

Simone Biondini, Nora Brambilla, Gramos Qerimi, Antonio Vairo

Discussion about hierarchies of scales improved with respect to the published version in the conference proceedings. 4 pages, 2 figures, conference proceeding of "NuDM-2022"

In this conference paper, we consider effective field theories of non-relativistic dark matter particles interacting with a light force mediator in the early expanding universe. We present a general framework, where to account in a systematic way for the relevant processes that may affect the dynamics during thermal freeze-out. In the temperature regime where near-threshold effects, most notably the formation of bound states and Sommerfeld enhancement, have a large impact on the dark matter relic density, we scrutinize possible contributions from higher excited states and radiative corrections in the annihilations and decays of dark-matter pairs.

Sergey V. Uchaikin, Boris I. Ivanov, Jinmyeong Kim, Çağlar Kutlu, Arjan F. Van Loo, Yasunobu Nakamura, Seonjeong OH, Violeta Gkika, Andrei Matlashov, Woohyun Chung, Yannis K. Semertzidis

6 pages, 7 figures, 29th International Conference on Low Temperature Physics, August 18-24, 2022, Sapporo, Japan

The axion is expected to solve the strong CP problem of quantum chromodynamics and is one of the leading candidates for dark matter. CAPP in South Korea has several axion search experiments based on cavity haloscopes in the frequency range of 1-6 GHz. The main effort focuses on operation of the experiments with the highest possible sensitivity. It requires maintenance of the haloscopes at the lowest physical temperature in the range of mK and usage of low noise components to amplify the weak axion signal. We report development and operation of low noise amplifiers for 5 haloscope experiments targeting at different frequency ranges. The amplifiers show noise temperatures approaching the quantum limit.