Abstract visibility
Abstract text size

Papers for Monday, Mar 18 2024

Papers with local authors

Ahmad Nemer, ChangHoon Hahn, Jiaxuan Li, Peter Melchior, Jeremy Goodman
1 vote
Paper 39 — arXiv:2403.10243
1 vote
Paper 39 — arXiv:2403.10243

Protoplanetary disks are the sites of vigorous hydrodynamic processes, such as accretion and outflows, and ultimately establish the conditions for the formation of planets. The properties of disk outflows are often inferred through analysis of forbidden emission lines. These lines contain multiple overlapping components, tracing different emission regions with different processes that excite them: a high-velocity component (tracing a jet), a broad low-velocity component (tracing inner disk wind), and a narrow low-velocity component (tracing outer disk wind). They are also heavily contaminated by background spectral features. All of these challenges call into question the traditional approach of fitting Gaussian components to the line profiles, and cloud the physical interpretation of those components. We introduce a novel statistical technique to analyze emission lines in protoplanetary disks. Simulation-Based Inference is a computationally efficient machine learning technique that produces posterior distributions of the parameters (e.g. magnetic field, radiation sources, geometry) of a representative wind model when given a spectrum, without any prior assumption about line shapes (e.g.symmetry). In this pathfinder study, we demonstrate that this technique indeed accurately recovers the parameters from simulated spectra without noise and background. A following work will deal with the analysis of observed spectra.

Danny Horta, Yuxi, Lu, Melissa K. Ness, Mariangela Lisanti, Adrian M. Price-Whelan

Submitted to ApJ. 15 pages, 13 figures, 1 Table

0 votes
Paper 7 — arXiv:2403.09777
0 votes
Paper 7 — arXiv:2403.09777

Reconstructing the mass assembly history of the Milky Way relies on obtaining detailed measurements of the properties of many stars in the Galaxy, especially in the stellar halo. One of the most constraining quantities is stellar age, as it can shed light on the accretion time and quenching of star formation in merging satellites. However, obtaining reliable age estimates for large samples of halo stars is difficult. We report published ages of 120 subgiant halo stars with highly-radial orbits that likely belong to the debris of the $Gaia-Enceladus/Sausage$ (GES) galaxy. The majority of these halo stars are old, with an age distribution characterized by a median of 11.6~Gyr and 16$^{\rm th}$(84$^{\rm th}$) percentile of 10.5~(12.7)~Gyr. However, the distribution is skewed, with a tail of younger stars that span ages down to $\sim6$-$9$ Gyr. All highly-radial halo stars have chemical and kinematic/orbital quantities that associate them with the GES debris. Initial results suggest that these intermediate-age stars are not a product of mass transfer and/or stellar mergers, which can bias their age determination low. If this conclusion is upheld by upcoming spectro-photometric studies, then the presence of these stars will pose an important challenge for constraining the properties of the GES merger and the accretion history of the Galaxy.

Hsun-Yeong Yang, Barry T. Chiang, Guan-Ming Su, Hsi-Yu Schive, Tzihong Chiueh, Jeremiah P. Ostriker

20 pages, 26 figures; Accepted for publication in MNRAS

0 votes
Paper 16 — arXiv:2403.09845
0 votes
Paper 16 — arXiv:2403.09845

Fuzzy dark matter (FDM), an attractive dark matter candidate comprising ultralight bosons (axions) with a particle mass $m_a\sim10^{-22}$ eV, is motivated by the small-scale challenges of cold dark matter and features a kpc-size de Broglie wavelength. Quantum wave interference inside an FDM halo gives rise to stochastically fluctuating density granulation; the resulting gravitational perturbations could drive significant disc thickening, providing a natural explanation for galactic thick discs. Here we present the first self-consistent simulations of FDM haloes and stellar discs, exploring $m_a=0.2-1.2\times10^{-22}$ eV and halo masses $M_\text{h} = 0.7-2.8\times10^{11}$ M$_\odot$. Disc thickening is observed in all simulated systems. The disc heating rates are approximately constant in time and increase substantially with decreasing $m_a$, reaching $dh/dt \simeq 0.04$ ($0.4$) kpc Gyr$^{-1}$ and $d\sigma_z^2/dt \simeq4$ ($150$) km$^2$s$^{-2}$Gyr$^{-1}$ for $m_a=1.2$ ($0.2$) $\times10^{-22}$ eV and $M_\text{h} =7\times10^{10} \text{M}_\odot$, where $h$ is the disc scale height and $\sigma_z$ is the vertical velocity dispersion. These simulated heating rates agree within a factor of two with the theoretical estimates of Chiang et al., confirming that the rough estimate of Church et al. overpredicts the granulation-driven disc heating rate by two orders of magnitude. However, the simulation-inferred heating rates scale less steeply than the theoretically predicted relation $d\sigma^2_z/dt \propto m_a^{-3}$. Finally, we examine the applicability of the Fokker-Planck approximation in FDM granulation modelling and the robustness of the $m_a$ exclusion bound derived from the Galactic disc kinematics.

The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y.X. Bai, Y.W. Bao, D. Bastieri, X.J. Bi, Y.J. Bi, J.T. Cai, Q. Cao, W.Y. Cao, Zhe Cao, J. Chang, J.F. Chang, A.M. Chen, E.S. Chen, Liang Chen, Lin Chen, Long Chen, M.J. Chen, M.L. Chen, Q.H. Chen, S.H. Chen, S.Z. Chen, T.L. Chen, Y. Chen, N. Cheng, Y.D. Cheng, M.Y. Cui, S.W. Cui, X.H. Cui, Y.D. Cui, B.Z. Dai, H.L. Dai, Z.G. Dai, Danzengluobu, D. della Volpe, X.Q. Dong, K.K. Duan, J.H. Fan, Y.Z. Fan, J. Fang, K. Fang, C.F. Feng, L. Feng, S.H. Feng, X.T. Feng, Y.L. Feng, S. Gabici, B. Gao, C.D. Gao, L.Q. Gao, Q. Gao, W. Gao, W.K. Gao, M.M. Ge, L.S. Geng, G. Giacinti, G.H. Gong, Q.B. Gou, M.H. Gu, F.L. Guo, X.L. Guo, Y.Q. Guo, Y.Y. Guo, Y.A. Han, H.H. He, H.N. He, J.Y. He, X.B. He, Y. He, M. Heller, Y.K. Hor, B.W. Hou, C. Hou, et al.

8 pages, 3 figures

0 votes
Paper 28 — arXiv:2403.10010
0 votes
Paper 28 — arXiv:2403.10010

We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components.

All other papers

Edwin A. Bergin, Arthur Bosman, Richard Teague, Jenny Calahan, Karen Willacy, L. Ilsedore Cleeves, Kamber Schwarz, Ke Zhang, Simon Bruderer

16 pages, 8 figures, accepted by the Astrophysical Journal

We present the first detection of 13CCH in a protoplanetary disk (TW Hya). Using observations of C2H we measure CCH/13CCH = 65 +/- 20 in gas with a CO isotopic ratio of 12CO/13CO = 21 +/- 5 (Yoshida et al. 2022a). The TW Hya disk exhibits a gas phase C/O that exceeds unity and C2H is the tracer of this excess carbon. We confirm that the TW Hya gaseous disk exhibits two separate carbon isotopic reservoirs as noted previously (Yoshida et al. 2022a). We explore two theoretical solutions for the development of this dichotomy. One model represents TW Hya today with a protoplanetary disk exposed to a cosmic ray ionization rate that is below interstellar as consistent with current estimates. We find that this model does not have sufficient ionization in cold (T < 40 K) layers to activate carbon isotopic fractionation. The second model investigates a younger TW Hya protostellar disk exposed to an interstellar cosmic ray ionization rate. We find that the younger model has sources of ionization deeper in a colder disk that generates two independent isotopic reservoirs. One reservoir is 12C-enriched carried by methane/hydrocarbon ices and the other is 13C-enriched carried by gaseous CO. The former potentially provides a source of methane/hydrocarbon ices to power the chemistry that generates the anomalously strong C$_2$H emission in this (and other) disk systems in later stages. The latter provides a source of gaseous 13C-rich material to generate isotopic enrichments in forming giant planets as recently detected in the super-Jupiter TYC 8998-760-1 b by Zhang et al. (2021).

A. Mishra, G. Mamatsashvili, M. Seilmayer, F. Stefani

11 pages, 5 figures, submitted to Journal of Fluid Mechanics

The effects of thermal convection on turbulence in accretion discs, and particularly its interaction with the magnetorotational instability (MRI), are of significant astrophysical interest. Despite extensive theoretical and numerical studies, such interactions have not been explored experimentally. We conduct linear analysis of the azimuthal version of MRI (AMRI) in the presence of convection and compare our findings with the experimental data. We show that AMRI sets in at lower critical Hartmann numbers ($Ha$) in the presence of convection. Importantly, convection breaks symmetry between the $m = \pm 1$ modes ($m$ is the azimuthal wavenumber). This preference for one mode over the other makes the AMRI-wave look like a "one-winged butterfly" observed in the experiments.

Colin Orion Chandler, Chadwick A. Trujillo, William J. Oldroyd, Jay K. Kueny, William A. Burris, Henry H. Hsieh, Jarod A. DeSpain, Nima Sedaghat, Scott S. Sheppard, Kennedy A. Farrell, David E. Trilling, Annika Gustafsson, Mark Jesus Mendoza Magbanua, Michele T. Mazzucato, Milton K. D. Bosch, Tiffany Shaw-Diaz, Virgilio Gonano, Al Lamperti, José A. da Silva Campos, Brian L. Goodwin, Ivan A. Terentev, Charles J. A. Dukes, Sam Deen

35 pages, 5 figures, 3 tables

We present the Citizen Science program Active Asteroids and describe discoveries stemming from our ongoing project. Our NASA Partner program is hosted on the Zooniverse online platform and launched on 2021 August 31, with the goal of engaging the community in the search for active asteroids -- asteroids with comet-like tails or comae. We also set out to identify other unusual active solar system objects, such as active Centaurs, active quasi-Hilda asteroids, and Jupiter-family comets (JFCs). Active objects are rare in large part because they are difficult to identify, so we ask volunteers to assist us in searching for active bodies in our collection of millions of images of known minor planets. We produced these cutout images with our project pipeline that makes use of publicly available Dark Energy Camera (DECam) data. Since the project launch, roughly 8,300 volunteers have scrutinized some 430,000 images to great effect, which we describe in this work. In total we have identified previously unknown activity on 15 asteroids, plus one Centaur, that were thought to be asteroidal (i.e., inactive). Of the asteroids, we classify four as active quasi-Hilda asteroids, seven as JFCs, and four as active asteroids, consisting of one Main-belt comet (MBC) and three MBC candidates. We also include our findings concerning known active objects that our program facilitated, an unanticipated avenue of scientific discovery. These include discovering activity occurring during an orbital epoch for which objects were not known to be active, and the reclassification of objects based on our dynamical analyses.

Manuel R. Izquierdo, Miguel Bezares, Steven Liebling, Carlos Palenzuela

16 pages, 11 figures. Comments are welcome!

The LIGO-Virgo-Kagra collaboration has observed gravitational waves consistent with the mergers of a black hole and a neutron star, namely GW200105 and GW200115, providing evidence for such cataclysmic events. Although no electromagnetic counterpart was reported for either of these two events, under certain conditions black hole--neutron star mergers are expected to form a significant accretion disk and to produce both a short gamma ray burst and a kilonova, much as observed in the binary neutron star merger GW170817. Here, we extend our publicly available code \mhduet to study numerically the merger of a magnetized neutron star with a black hole. \mhduet employs Large Eddy Simulation techniques to help capture the magnetic field amplification resulting from turbulence and other sub-grid scale dynamics in the post-merger stage. In particular, we simulate a merger with parameters favorable to producing an accretion disk, focusing on the formation and dynamics of the turbulent disk and the resulting magnetic field amplification. Following the tidal disruption and during the formation of the accretion disk, the magnetic field undergoes significant amplification driven by the Kelvin-Helmholtz instability, reaching strengths of more than $10^{14}\,\rm{G}$ from a realistic initial strength of $10^{11}\,\rm{G}$ in short timescales of approximately $20\,\rm{ms}$. At later times in our simulations the magnetic field growth occurs at larger scales and is dominated mainly by magnetic winding.

O. Springer, E. O. Ofek, B. Zackay, R. Konno, A. Sharon, G. Nir, A. Rubin, A. Haddad, J. Friedman, L. Schein Lubomirsky, I. Aizenberg, A. Krassilchtchikov, A. Gal-Yam

13 pages, 11 figures, submitted to AJ

Detection of moving sources over complicated background is important for several reasons. First is measuring the astrophysical motion of the source. Second is that such motion resulting from atmospheric scintillation, color refraction, or astrophysical reasons is a major source of false alarms for image subtraction methods. We extend the Zackay, Ofek, and Gal-Yam image subtraction formalism to deal with moving sources. The new method, named translient (translational transient) detector, applies hypothesis testing between the hypothesis that the source is stationary and that the source is moving. It can be used to detect source motion or to distinguish between stellar variability and motion. For moving source detection, we show the superiority of translient over the proper image subtraction, using the improvement in the receiver-operating characteristic curve. We show that in the small translation limit, Translient is an optimal detector of point source motion in any direction. Furthermore, it is numerically stable, fast to calculate, and presented in a closed form. Efficient transient detection requires both the proper image subtraction statistics and the translient statistics: when the translient statistic is higher, then the subtraction artifact is likely due to motion. We test our algorithm both on simulated data and on real images obtained by the Large Array Survey Telescope (LAST). We demonstrate the ability of translient to distinguish between motion and variability, which has the potential to reduce the number of false alarms in transients detection. We provide the translient implementation in Python and MATLAB.

Haruki Katoh, Chikako Yasui, Yuji Ikeda, Naoto Kobayashi, Noriyuki Matsunaga, Sohei Kondo, Hiroaki Sameshima, Satoshi Hamano, Misaki Mizumoto, Hideyo Kawakita, Kei Fukue, Shogo Otsubo, Keiichi Takenaka

Accepted for publication in ApJ

In this study, we performed high-resolution near-infrared (NIR) spectroscopy (R = 28,000; ${\lambda} = 0.90$-1.35 ${\mu}$m) with a high signal-to-noise ratio on HD 200775, a very young (${\sim}$0.1 Myr old) and massive intermediate-mass star (a binary star with a mass of about 10 $M_{\odot}$ each) with a protoplanetary disk. The obtained spectra show eight forbidden lines of three elements: two of [S II] (10289 and 10323 {\AA}), two of [N I] (10400 and 10410 {\AA}), and four of [Fe II] (12570, 12946, 12981, and 13209 {\AA}). This is the first time that the [N I] lines are detected in a young stellar object with a doublet deblended. Gaussian fitting of the spectra indicates that all line profiles have low-velocity components and exhibit blueshifted features, suggesting that all lines originate from the disk winds (magnetohydrodynamic disk wind and/or photoevaporative wind). Based on the fit, the [N I] and [Fe II] lines are categorized into narrow components, while the [S II] lines are at the boundary between broad and narrow components. These forbidden lines are suggested to be very promising disk wind tracers among the existing ones because they are in the NIR-wavelength range, which can be observed from early stages with high sensitivities. Among these lines, [N I] lines would be a rather powerful probe for deriving the basic physical parameters of disk wind gases. However, the study of these lines herein is limited to one object; thus, further studies are needed to examine their properties.

The formation of circumstellar discs is a critical step in the formation of stars and planets. Magnetic fields can strongly affect the evolution of angular momentum during prestellar core collapse, potentially leading to the failure of protostellar disc formation. This phenomenon, known as the magnetic braking catastrophe, has been observed in ideal-MHD simulations. In this work, we present results from ideal-MHD simulations of circumstellar disc formation from realistic initial conditions of strongly magnetised, massive cores with masses between $30 ~{\rm M}_\odot$ and $300 ~{\rm M}_\odot$ resolved by zooming into Giant Molecular Clouds with masses $\sim 10^4 \ {\rm M}_\odot$ and initial mass-to-flux ratios $0.6 \le \mu_0 \le 3$. Due to the large turbulence caused by the non-axisymmetric gravitational collapse of the gas, the dominant vertical support of discs is turbulent motion, while magnetic and turbulent pressures contribute equally in the outer toroid. The magnetic field topology is extremely turbulent and incoherent, reducing the effect of magnetic braking by roughly one order of magnitude and leading to the formation of large Keplerian discs even in magnetically critical or near-critical cores. Only cores in GMCs with $\mu_0 < 1$ fail to form discs. Instead, they collapse into a sheet-like structure and produce numerous low-mass stars. We also discuss a universal $B-\rho$ relation valid over a large range of scales from the GMC to massive cores, irrespective of the GMC magnetisation. This study differs from the vast literature on this topic which typically focus on smaller mass discs with idealised initial and boundary conditions, therefore providing insights into the initial conditions of massive prestellar core collapse and disc formation.

Andrew D. Sellek, Naman S. Bajaj, Ilaria Pascucci, Cathie J. Clarke, Richard Alexander, Chengyan Xie, Giulia Ballabio, Dingshan Deng, Uma Gorti, Andras Gaspar, Jane Morrison

32 pages, 16 figures, Accepted 14/03/24 to the Astronomical Journal. Complementary modeling to Bajaj et al. 2024 ( arXiv:2403.01060 )

[Ne II] 12.81 $\mu\mathrm{m}$ emission is a well-used tracer of protoplanetary disk winds due to its blueshifted line profile. MIRI-MRS recently observed T Cha, detecting this line along with lines of [Ne III], [Ar II] and [Ar III], with the [Ne II] and [Ne III] lines found to be extended while the [Ar II] was not. In this complementary work, we use these lines to address long-debated questions about protoplanetary disk winds regarding their mass-loss rate, the origin of their ionization, and the role of magnetically-driven winds as opposed to photoevaporation. To this end, we perform photoionization radiative transfer on simple hydrodynamic wind models to map the line emission. We compare the integrated model luminosities to those observed with MIRI-MRS to identify which models most closely reproduce the data and produce synthetic images from these to understand what information is captured by measurements of the line extents. Along with the low degree of ionization implied by the line ratios, the relative compactness of [Ar II] compared to [Ne II] is particularly constraining. This requires Ne II production by hard X-rays and Ar II production by soft X-rays (and/or EUV) in an extended ($\gtrsim 10$ au) wind that is shielded from soft X-rays - necessitating a dense wind with material launched on scales down to ~1 au. Such conditions could be produced by photoevaporation, whereas an extended MHD wind producing equal shielding would likely underpredict the line fluxes. However, a tenuous inner MHD wind may still contribute to shielding the extended wind. This picture is consistent with constraints from spectrally-resolved line profiles.

Atila Poro, Mahya Hedayatjoo, Maryam Nastaran, Mahshid Nourmohammad, Hossein Azarara, Sepideh AlipourSoudmand, Fatemeh AzarinBarzandig, Razieh Aliakbari, Sadegh Nasirian, Nazanin Kahali Poor

Accepted by the New Astronomy Journal

The accuracy of absolute parameters' estimation in contact binary systems is important for investigating their evolution and solving some challenges. The Gaia DR3 parallax is one of the methods used for estimating the absolute parameters, in cases where photometric data is the only one that is available. The use of this method includes advantages and limitations that we have described and examined in this study. We selected 48 contact binary systems whose mass ratios were mostly obtained by spectroscopic data, in addition to a number of photometric studies. The target systems were suitable for A_V and the Re-normalised Unit Weight Error (RUWE), and their absolute parameters were calculated based on Gaia DR3 parallax, observational information, orbital period, and light cure solution from the literature and catalogs. The outcomes of OO Aql differed significantly from those reported in the literature. Upon analyzing the system's light curve with TESS data, we concluded that the stars' temperatures were the reason for this difference, and utilizing Gaia DR3 parallax provided reasonable results. We displayed the target systems on the Hertzsprung-Russell (HR), q-L_ratio, P-M_V, and logM_tot-logJ_0 diagrams, and the systems are in good agreement with the theoretical fits. We showed that the estimation of absolute parameters with this method might be acceptable if Delta a(R_Sun) is less than about 0.1(R_Sun). There are open questions regarding the existence of l_3 in the light curve analysis and its effect on the estimation of absolute parameters with this method.

N. Biava, A. Bonafede, F. Gastaldello, A. Botteon, M. Brienza, T. W. Shimwell, G. Brunetti, L. Bruno, K. Rajpurohit, C. J. Riseley, R.J. van Weeren, M. Rossetti, R. Cassano, F. De Gasperin, A. Drabent, H.J.A. Rottgering, A.C. Edge, C. Tasse

34 pages, 30 figures. Accepted for publication in A&A

Radio observations of a few cool-core galaxy clusters have revealed the presence of diffuse emission on cluster scales, similar to what was found in merging clusters in the form of radio halos. These sources might suggest that a minor merger, while not sufficiently energetic to disrupt the cool core, could still trigger particle acceleration in the intracluster medium on scales of hundreds of kpc. We observed with LOFAR at 144 MHz a sample of twelve cool-core galaxy clusters presenting some level of dynamical disturbances, according to X-ray data. We also performed a systematic search of cold fronts in these clusters, re-analysing archival Chandra data. The clusters PSZ1G139.61+24, A1068 (new detection), MS 1455.0+2232, and RX J1720.1+2638 present diffuse radio emission on a cluster scale. This emission is characterised by a double component: a central mini-halo confined by cold fronts and diffuse emission on larger scales, whose radio power at 144 MHz is comparable to that of radio halos detected in merging systems. The cold fronts in A1068 are a new detection. We also found a candidate plasma depletion layer in this cluster. No sloshing features are found in the other eight clusters. Two of them present a mini-halo, with diffuse radio emission confined to the cluster core. We also found a new candidate mini-halo. Whereas, for the remaining five clusters, we did not detect halo-like emission. For clusters without cluster-scale halos, we derived upper limits to the radio halo power. We found that cluster-scale diffuse radio emission is not present in all cool-core clusters when observed at a low frequency, but it is correlated to the presence of cold fronts. This morphology requires a specific configuration of the merger and so it puts some constraints on the turbulence, which deserves to be investigated in the future with theoretical works.

E. Bulbul, X. Zhang, M.Kluge, M. Brueggen, B. Koribalski, A. Liu, E. Artis, Y. E. Bahar, F. Balzer, C. Garrel, V. Ghirardini, N. Malavasi, A. Merloni, K. Nandra, M. E. Ramos-Ceja, J. S. Sanders, S. Zelmer

Submitted to A&A Letters. 8 pages, 6 figures, 2 tables

Odd radio circles (ORC) are a newly discovered class of extended faint radio sources of unknown origin. We report the first detection of diffuse X-ray gas at the location of a low-redshift ORC (z=0.046), known as Cloverleaf ORC. This observation was performed with the XMM-Newton X-ray telescope. The physical extent of the diffuse X-ray emission corresponds to approximately a 230 kpc by 160 kpc region, lying perpendicular to the radio emission detected by ASKAP. The X-ray spectrum shows characteristics of thermal multi-phase gas with temperatures of $1.10\pm0.08$ keV and $0.22\pm0.01$ keV and a central density of $(4.9\pm0.6)\times10^{-4}$ cm$^{-3}$, indicating that the Cloverleaf ORC resides in a low-mass galaxy group. Using velocity dispersion measurements in the optical, the halo mass of the galaxy group has been reported as $\sim1.6\times10^{13}$ M$_{sun}$. The presence of a high-velocity subgroup identified in optical data, the orientation of the BCG, the disturbed morphologies of galaxies towards the east of the Cloverleaf ORC, and the irregular morphology of the X-ray emission suggest that this system is undergoing a galaxy group merger. The radio power of the ORC could be explained by the shock reacceleration of fossil cosmic rays generated by a previous episode of black hole activity in the central active galactic nucleus.

Fabio Rigamonti, Luca Cortese, Francesco Bollati, Stefano Covino, Massimo Dotti, A. Fraser-McKelvie, Francesco Haardt

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

We investigate the relevance of kinematically identified bulges, discs and their role relative to galaxy quenching. We utilize an analysis of the SDSS-MaNGA survey conducted with the GPU-based code BANG which simultaneously models galaxy photometry and kinematics to decompose galaxies into their structural components. Below M~1011 Msun, galaxies exhibit a wide range of dynamical properties, determined by the relative prominence of a dispersion-supported inner region and a rotationally-supported disc. Our analysis reveals a natural separation between these classes, with only a minor fraction of stellar mass retained by structures exhibiting intermediate dynamical support. When examining galaxies in terms of their star formation activity, an apparent decrease in rotational support is observed as they move below the star-forming main sequence. This behaviour is evident with luminosity-weighted tracers of kinematics, while it almost vanishes with mass-weighted tracers. Luminosity-weighted quantities not only capture differences in kinematics but also in the stellar population, potentially leading to biased interpretations of galaxy dynamical properties and quenching. Our findings suggest that quenching does not imply almost any structural transformation in galaxies below M~10^11 Msun. Processes as disc fading more likely account for observed differences in mass-weighted and luminosity-weighted galaxy properties; when the galactic disc ceases star formation, its mass-to-light ratio grows without any significant morphological transformation. The picture is remarkably different above M~10^11 Msun. Regardless of the tracer used, a substantial increase in galaxy dispersion support is observed along with a significant structural change. A different quenching mechanism, most likely associated with mergers, dominates. Notably, this mechanism is confined to a very limited range of high masses.

Kartick C. Sarkar

Invited article at The Astronomy and Astrophysics Review. This is the author's version. Comments are welcome

Galactic outflows are ubiquitous in galaxies containing active star formation or supermassive black hole activity. The presence of a large-scale outflow from the center of our own Galaxy was confirmed after the discovery of two large ($\sim 8-10$ kpc) \gamma-ray bubbles using the \textit{Fermi-LAT} telescope. These bubbles, known as the Fermi Bubbles, are highly symmetric about the Galactic disk as well as around the Galactic rotation axis and appear to emanate from the center of our Galaxy. The sharp edges of these bubbles suggest that they are related to the Galactic outflow. These bubbles are surrounded by two even bigger ($\sim 12-14$ kpc) X-ray structures, known as the eROSITA bubbles. Together, they represent the characteristics of an outflow from the Galaxy into the circumgalactic medium. Multi-wavelength observations such as in radio, microwave, and UV toward the Fermi Bubbles have provided us with much information in the last decade. However, the origin and the nature of these bubbles remain elusive. In this review, I summarize the observations related to the Fermi/eROSITA Bubbles at different scales and wavelengths, and give a brief overview of our current understanding of them.

E. Goffo, P. Chaturvedi, F. Murgas, G. Morello, J. Orell-Miquel, L. Acuña, L. Peña-Moñino, E. Pallé, A. P. Hatzes, S. Geraldía-González, F. J. Pozuelos, A. F. Lanza, D. Gandolfi, J. A. Caballero, M. Schlecker, M. Pérez-Torres, N. Lodieu, A. Schweitzer, C. Hellier, S. V. Jeffers, C. Duque-Arribas, C. Cifuentes, V. J. S. Béjar, M. Daspute, F. Dubois, S. Dufoer, E. Esparza-Borges, A. Fukui, Y. Hayashi, E. Herrero, M. Mori, N. Narita, H. Parviainen, L. Tal-Or, S. Vanaverbeke, I. Hermelo, P. J. Amado, S. Dreizler, Th. Henning, J. Lillo-Box, R. Luque, M. Mallorquín, E. Nagel, A. Quirrenbach, S. Reffert, A. Reiners, I. Ribas, P. Schöfer, H. M. Tabernero, M. Zechmeister

Accepted for publication on Astronomy & Astrophysics

We report the confirmation and mass determination of a mini-Neptune transiting the M3.5 V star TOI-4438 (G 182-34) every 7.44 days. A transit signal was detected with NASA's TESS space mission in the sectors 40, 52, and 53. In order to validate the planet TOI-4438 b and to determine the system properties, we combined TESS data with high-precision radial velocity measurements from the CARMENES spectrograph, spanning almost one year, and ground-based transit photometry. We found that TOI-4438 b has a radius of Rb = 2.52 +/- 0.13 R_Earth (5% precision), which together with a mass of Mb=5.4 +/- 1.1 M_Earth (20% precision), results in a bulk density of rho = 1.85+0.51-0.44 g cm-3 (28% precision), aligning the discovery with a volatile-rich planet. Our interior structure retrieval with a pure water envelope yields a minimum water mass fraction of 46% (1-sigma). TOI-4438 b is a volatile-rich mini-Neptune with likely H/He mixed with molecules, such as water, CO_2, and CH_4. The primary star has a J-band magnitude of 9.7, and the planet has a high transmission spectroscopy metric (TSM) of 136 +/- 13. Taking into account the relatively warm equilibrium temperature of T_eq = 435 +/- 15 K, and the low activity level of its host star, TOI-4438 b is one of the most promising mini-Neptunes around an M dwarf for transmission spectroscopy studies.

Simon R. Dicker, Karen Perez Sarmiento, Brian Mason, Tanay Bhandarkar, Mark J. Devlin, Luca Di Mascolo, Saianeesh Haridas, Matt Hilton, Mathew Madhavacheril, Emily Moravec, Tony Mroczkowski, John Orlowski-Scherer, Charles Romero, Craig L. Sarazin, Jonathan Sievers

11 pages, 7 figures and 2 tables Extended version of figure 4 is included and the full data for table 1 can be found as the auxiliary file tab1.txt Submitted to ApJ

In this paper we present the results of a blind survey for compact sources in 243 Galaxy clusters that were identified using the thermal Sunyaev-Zeldovich effect (tSZ). The survey was carried out at 90 GHz using MUSTANG2 on the Green Bank telescope and achieved a $5\sigma$ detection limit of 1 mJy in the center of each cluster. We detected 24 discrete sources. The majority (18) of these correspond to known radio sources, and of these, 5 show signs of significant variability, either with time or in spectral index. The remaining sources have no clear counterparts at other wavelengths. Searches for galaxy clusters via the tSZ effect strongly rely on observations at 90 GHz, and the sources found have the potential to bias mass estimates of clusters. We compare our results to the simulation Websky that can be used to estimate the source contamination in galaxy cluster catalogs. While the simulation showed a good match to our observations at the clusters' centers, it does not match our source distribution further out. Sources over 104" from a cluster's center bias the tSZ signal high, for some of our sources, by over 50%. When averaged over the whole cluster population the effect is smaller but still at a level of 1 to 2%. We also discovered that unlike previous measurements and simulations we see an enhancement of source counts in the outer regions of the clusters and fewer sources than expected in the centers of this tSZ selected sample.

Sheila Sagear, Adrian M. Price-Whelan, Sarah Ballard, Yuxi (Lucy)Lu, Ruth Angus, David W. Hogg

23 pages, 12 figures, 2 tables. Submitted to AAS Journals

Stellar age measurements are fundamental to understanding a wide range of astronomical processes, including galactic dynamics, stellar evolution, and planetary system formation. However, extracting age information from Main Sequence stars is complicated, with techniques often relying on age proxies in the absence of direct measurements. The Gaia data releases have enabled detailed studies of the dynamical properties of stars within the Milky Way, offering new opportunities to understand the relationship between stellar age and dynamics. In this study, we leverage high-precision astrometric data from Gaia DR3 to construct a stellar age prediction model based only on stellar dynamical properties; namely, the vertical action. We calibrate two distinct, hierarchical stellar age--vertical action relations, first employing asteroseismic ages for red giant branch stars, then isochrone ages for main-sequence turn-off stars. We describe a framework called "zoomies" based on this calibration, by which we can infer ages for any star given its vertical action. This tool is open-source and intended for community use. We compare dynamical age estimates from "zoomies" with ages derived from other techniques for a sample of open clusters and main-sequence stars with asteroseismic age measurements. We also compare dynamical age estimates for stellar samples from the Kepler, K2, and TESS exoplanet transit surveys. While dynamical age relations are associated with large uncertainty, they are generally mass-independent and depend on homogeneously measured astrometric data. These age predictions are uniquely useful for large-scale demographic investigations, especially in disentangling the relationship between planet occurrence, metallicity, and age for low-mass stars.

Europa's surface features many regions of complex topography termed 'chaos terrains'. One set of hypotheses for chaos terrain formation require upward migration of liquid water from perched water bodies formed by convection and tidal heating within the icy shell. However, consideration of the behaviour of terrestrial ice sheets suggests the rapid downwards, not upwards, movement of water from perched water bodies, initiated through hydrofracture and without a propagation limit provided there is a sufficient volume of water to fill the resulting vertical fracture. I suggest that drainage of perched water bodies to the subsurface ocean through hydrofracture is a likely phenomenon at Europa and other icy moons. This may at first appear to rule out perched water bodies in the formation of chaos areas, but the violence of such drainage events and loss of mechanical support could lead to the collapse of the surface shell and a temporary surface-ocean connection. Alternatively, hydrofracture provides a mechanism for entire ice-shell fracture if a surface impact melts enough near-surface ice.

B.P. Blakley, Will M. Grundy, Jordan K. Steckloff, Sugata P. Tan, Jennifer Hanley, Anna E. Engle, Stephen C. Tegler, Gerrick E. Lindberg, Shae M. Raposa, Kendall J. Koga, Cecilia L. Thieberger

29 pages, 9 figures, to be published in Planetary and Space Science

Few laboratory studies have investigated the vapor pressures of the volatiles that may be present as ices in the outer solar system; even fewer studies have investigated these species at the temperatures and pressures suitable to the surfaces of icy bodies in the Saturnian and Uranian systems ($\lt$100 K, $\lt10^{-9}$ bar). This study adds to the work of Grundy et al. (2024) in extending the known equilibrium vapor pressures of outer solar system ices through laboratory investigations at very low temperatures. Our experiments with ammonia and oxygen ices provide new thermodynamic models for these species' respective enthalpies of sublimation. We find that ammonia ice, and to a lesser degree oxygen ice, are stable at higher temperatures than extrapolations in previous literature have predicted. Our results show that these ices should be retained over longer periods of time than previous extrapolations would predict, and a greater amount of these solids is required to support observation in exospheres of airless bodies in the outer solar system.

Sungsoon Lim, Eric W. Peng, Patrick Côté, Laura Ferrarese, Joel C. Roediger, Chengze Liu, Chelsea Spengler, Elisabeth Sola, Pierre-Alain Duc, Laura V. Sales, John P. Blakeslee, Jean-Charles Cuillandre, Patrick R. Durrell, Eric Emsellem, Stephen D. J. Gwyn, Ariane Lançon, Francine R. Marleau, J. Christopher Mihos, Oliver Müller, Thomas H. Puzia, Rubén Sánchez-Janssen

28 pages, 18 Figures, 3 tables, accepted for publication in ApJ

We study the size and structure of globular clusters (GC) systems of 118 early-type galaxies from the NGVS, MATLAS, and ACSVCS surveys. Fitting S\'ersic profiles, we investigate the relationship between effective radii of GC systems ($R_{e, \rm gc}$) and galaxy properties. GC systems are 2--4 times more extended than host galaxies across the entire stellar mass range of our sample ($10^{8.3} < M_* < 10^{11.6}~M_{\odot}$). The relationship between $R_{e, \rm gc}$ and galaxy stellar mass exhibits a characteristic "knee" at a stellar mass of $M_p \simeq 10^{10.8}$, similar to galaxy $R_e$--stellar mass relationship. We present a new characterization of the traditional blue and red GC color sub-populations, describing them with respect to host galaxy $(g'-i')$ color ($\Delta_{gi}$): GCs with similar colors to their hosts have a "red" $\Delta_{gi}$, and those significantly bluer GCs have a "blue" $\Delta_{gi}$. The GC populations with red $\Delta_{gi}$, even in dwarf galaxies, are twice as extended as the stars, suggesting that formation or survival mechanisms favor the outer regions. We find a tight correlation between $R_{e, \rm gc}$ and the total number of GCs, with intrinsic scatter $\lesssim 0.1$ dex spanning two and three orders of magnitude in size and number, respectively. This holds for both red and blue subpopulations, albeit with different slopes. Assuming that $N_{GC, Total}$ correlates with $M_{200}$, we find that the red GC systems have effective radii of roughly 1-5\% $R_{\rm 200}$, while the blue GC systems in massive galaxies can have sizes as large as $\sim$10\% $R_{\rm 200}$. Environmental dependence on $R_{e, \rm gc}$ is also found, with lower density environments exhibiting more extended GC systems at fixed mass.

Jack M. Jenkins, Christopher M. J. Osborne, Ye Qiu, Rony Keppens, Chuan Li

9 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters (ApJL)

Solar prominences observed close to the limb commonly include a bright feature that, from the perspective of the observer, runs along the interface between itself and the underlying chromosphere. Despite several idealised models being proposed to explain the underlying physics, a more general approach remains outstanding. In this manuscript we demonstrate as a proof-of-concept the first steps in applying the Lightweaver radiative transfer framework's 2.5D extension to a `toy' model prominence + VAL3C chromosphere, inspired by recent 1.5D experiments that demonstrated a significant radiative chromosphere--prominence interaction. We find the radiative connection to be significant enough to enhance both the electron number density within the chromosphere, as well as its emergent intensity across a range of spectral lines in the vicinity of the filament absorption signature. Inclining the viewing angle from the vertical, we find these enhancements to become increasingly asymmetric and merge with a larger secondary enhancement sourced directly from the prominence underside. In wavelength, the enhancements are then found to be the largest in both magnitude and horizontal extent for the spectral line cores, decreasing into the line wings. Similar behaviour is found within new Chinese H$\alpha$ Solar Explorer (CHASE)/H$\alpha$ Imaging Spectrograph (HIS) observations, opening the door for subsequent statistical confirmations of the theoretical basis we develop here.

Andrew Sullivan, Chris Power, Connor Bottrell, Aaron Robotham, Stas Shabala

This paper has been accepted for publication with PASA

We investigate the connection between a cluster's structural configuration and observable measures of its gas emission that can be obtained in X-ray and Sunyaev-Zeldovich (SZ) surveys. We present an analytic model for the intracluster gas density profile: parameterised by the dark matter halo's inner logarithmic density slope, $\alpha$, the concentration, $c$, the gas profile's inner logarithmic density slope, $\varepsilon$, the dilution, $d$, and the gas fraction, $\eta$, normalised to cosmological content. We predict four probes of the gas emission: the emission-weighted, $T_\mathrm{X}$, and mean gas mass-weighted, $T_\mathrm{m_g}$, temperatures, and the spherically, $Y_\mathrm{sph}$, and cylindrically, $Y_\mathrm{cyl}$, integrated Compton parameters. Over a parameter space of clusters, we constrain the X-ray temperature scaling relations, $M_{200} - T_\mathrm{X}$ and $M_{500} - T_\mathrm{X}$, within $57.3\%$ and $41.6\%$, and $M_{200} - T_\mathrm{m_g}$ and $M_{500} - T_\mathrm{m_g}$, within $25.7\%$ and $7.0\%$, all respectively. When excising the cluster's core, the $M_{200} - T_\mathrm{X}$ and $M_{500} - T_\mathrm{X}$ relations are further constrained, to within $31.3\%$ and $17.1\%$, respectively. Similarly, we constrain the SZ scaling relations, $M_{200} - Y_\mathrm{sph}$ and $M_{500} - Y_\mathrm{sph}$, within $31.1\%$ and $17.7\%$, and $M_{200} - Y_\mathrm{cyl}$ and $M_{500} - Y_\mathrm{cyl}$, within $25.2\%$ and $22.0\%$, all respectively. The temperature observable $T_\mathrm{m_g}$ places the strongest constraint on the halo mass, whilst $T_\mathrm{X}$ is more sensitive to the parameter space. The SZ constraints are sensitive to the gas fraction, whilst insensitive to the form of the gas profile itself. In all cases, the halo mass is recovered with an uncertainty that suggests the cluster's structural profiles only contribute a minor uncertainty in its scaling relations.

Shuai Feng, Shi-Yin Shen, Fang-Ting Yuan, Wen-Xin Zhong, Wen-Yuan Cui, Lin-Lin Li

10 pages, 4 figures, accepted for publication in ApJ

We investigate the suppression of star formation in galaxy pairs based on the isolated galaxy pair sample derived from the SDSS survey. By comparing the star formation rate between late-type galaxies in galaxy pairs and those in the isolated environment, we detect the signal of star formation suppression in galaxy pairs at $d_p < 100$kpc and $200$kpc$ < d_p < 350$kpc. The occurrence of star formation suppression in these late-type galaxies requires their companion galaxies to have an early-type morphology ($n_s > 2.5$). Star formation suppression in wide galaxy pairs with $200$kpc$ < d_p < 350$kpc mainly occurs in massive late-type galaxies, while in close galaxy pairs with $d_p < 100$kpc, it only appears in late-type galaxies with a massive companion ( $\log M_\star > 11.0$), nearly independent of their own stellar mass. Based on these findings, we infer that star formation suppression in wide galaxy pairs is actually a result of galaxy conformity, while in close galaxy pairs, it stems from the influence of hot circum-galactic medium surrounding companion galaxies.

Rui-Jing Lu, Wen-Hao Chen, Wen-Qiang Liang, Cheng-Feng Peng

Comments: 8 pages, 8 figures,Accepted for publication in MNRAS

The properties of the progenitors of gamma-ray bursts (GRBs) and of their environment are encoded in their luminosity function and cosmic formation rate. They are usually recovered from a flux-limited sample based on Lynden-Bell's $c^{-}$ method. However, this method is based on the assumption that the luminosity is independent of the redshift. Observationally, if correlated, people use nonparametric $\tau$ statistical method to remove this correlation through the transformation, $L^{\prime}=L/g(z)$, where $z$ is the burst redshift, and $g(z)=(1+z)^{k}$ parameterizes the underlying luminosity evolution. However, the application of this method to different observations could result in very different luminosity functions. By the means of Monte Carlo simulation, in this paper, we demonstrate that the origin of an observed correlation, measured by the $\tau$ statistical method, is a complex combination of multiple factors when the underlying data are correlated. Thus, in this case, it is difficult to unbiasedly reconstruct the underlying population distribution from a truncated sample, unless the detailed information of the intrinsic correlation is accurately known in advance. In addition, we argue that an intrinsic correlation between luminosity function and formation rate is unlikely eliminated by a misconfigured transformation, and the $g(z)$, derived from a truncated sample with the $\tau$ statistical method, does not necessarily represent its underlying luminosity evolution.

Eric Gaidos, Thanawuth Thanathibodee, Andrew Hoffman, Joel Ong, Jason Hinkle, Benjamin J. Shappee, Andrea Banzatti

Accepted to The Astrophysical Journal

Transition disks, with inner regions depleted in dust and gas, could represent later stages of protoplanetary disk evolution when newly-formed planets are emerging. The PDS 70 system has attracted particular interest because of the presence of two giant planets at tens of au orbits within the inner disk cavity, at least one of which is itself accreting. However, the region around PDS 70 most relevant to understanding the planet populations revealed by exoplanet surveys of middle-aged stars is the inner disk, which is the dominant source of the system's excess infrared emission but only marginally resolved by ALMA. Here we present and analyze time-series optical and infrared photometry and spectroscopy that reveal the inner disk to be dynamic on timescales of days to years, with occultation of sub-micron dust dimming the star at optical wavelengths and 3-5 $\mu$m emission varying due to changes in disk structure. Remarkably, the infrared emission from the innermost region (nearly) disappears for ~1 year. We model the spectral energy distribution of the system and its time variation with a flattened warm (T <~ 600K) disk and a hotter (1200K) dust that could represent an inner rim or wall. The high dust-to-gas ratio of the inner disk relative to material accreting from the outer disk, means that the former could be a chimera consisting of depleted disk gas that is subsequently enriched with dust and volatiles produced by collisions and evaporation of planetesimals in the inner zone.

Maksym Mohorian (1 and 2), Devika Kamath (1, 2 and 3), Meghna Menon (1 and 2), Paolo Ventura (3 and 4), Hans Van Winckel (5), D. A. García-Hernández (6 and 7), Thomas Masseron (6 and 7) ((1) School of Mathematical and Physical Sciences, Macquarie University, Sydney, NSW, Australia, (2) Astrophysics and Space Technologies Research Centre, Macquarie University, Sydney, NSW, Australia, (3) INAF, Osservatorio Astronomico di Roma, Monte Porzio Catone, Italy, (4) Osservatorio Astronomico di Roma, Monte Porzio Catone, Lazio, Italy, (5) Institute of Astronomy, KU Leuven, Leuven, Belgium, (6) Instituto de Astrofísica de Canarias (IAC), La Laguna, Tenerife, Spain, (7) Departamento de Astrofísica, Universidad de La Laguna (ULL), La Laguna, Tenerife, Spain)

Dusty post-red giant branch (post-RGB) stars are low- and intermediate-mass stars where the RGB evolution was prematurely terminated by a poorly understood binary interaction. These binary stars are considered to be low-luminosity analogues of post-asymptotic giant branch (post-AGB) binary stars. In this study, we investigated the chemical composition of two dusty post-RGB binary stars, SZ Mon and DF Cyg, using multi-wavelength spectroscopic data from HERMES/Mercator (optical) and the APOGEE survey (near-infrared). Owing to challenges posed by existing spectral analysis tools for the study of evolved stars with complex atmospheres, we developed E-iSpec: a dedicated spectral analysis tool for evolved stars, to consistently determine atmospheric parameters, elemental abundances, and carbon isotopic ratios. Our abundance analysis revealed that observed depletion patterns and estimated depletion efficiencies resemble those found in post-AGB binary stars. However, the onset of chemical depletion in post-RGB targets occurs at higher condensation temperatures ($T_{\rm turn-off, post-RGB}\approx1400$ K), than in most post-AGB stars ($T_{\rm turn-off, post-AGB}\approx1100$ K). Additionally, our study resulted in the first estimates of carbon isotopic ratios for post-RGB stars ($^{12}$C/$^{13}$C$_{\rm SZ Mon}=8\pm4$, $^{12}$C/$^{13}$C$_{\rm DF Cyg}=12\pm3$). We found that the observationally derived CNO abundances and the carbon isotopic ratios of our post-RGB binary targets are in good agreement with theoretical predictions from the ATON single star evolutionary models involving first dredge-up and moderately-deep extra mixing. This agreement emphasises that in post-RGB binary targets, the observed CNO abundances reflect the chemical composition expected from single star nucleosynthesis (i.e., convective and non-convective mixing processes) occurring during the RGB phase before it is terminated.

Gregory P. Sercel, Pradip R. Gatkine, Nemanja Jovanovic, Jeffrey B. Jewell, Luis Pereira da Costa, J. Kent Wallace, Dimitri P. Mawet

High-resolution astronomical spectroscopy carried out with a photonic Fourier transform spectrograph (FTS) requires long asymmetrical optical delay lines that can be dynamically tuned. For example, to achieve a spectral resolution of R = 30,000, a delay line as long as 1.5 cm would be required. Such delays are inherently prone to phase errors caused by temperature fluctuations. This is due to the relatively large thermo-optic coefficient and long lengths of the waveguides, in this case composed of SiN, resulting in thermally dependent changes to the optical path length. To minimize phase error to the order of 0.05 radians, thermal stability of the order of 0.05{\deg} C is necessary. A thermal control system capable of stability such as this would require a fast thermal response and minimal overshoot/undershoot. With a PID temperature control loop driven by a Peltier cooler and thermistor, we minimized interference fringe phase error to +/- 0.025 radians and achieved temperature stability on the order of 0.05{\deg} C. We present a practical system for precision temperature control of a foundry-fabricated and packaged FTS device on a SiN platform with delay lines ranging from 0.5 to 1.5 cm in length using inexpensive off-the-shelf components, including design details, control loop optimization, and considerations for thermal control of integrated photonics.

Wei Zhu, Xu-Rong Chen, Yu-Chen Tang

27 pages and 6 figures

We use a special gluon distribution in the nucleon, which was predicted by a QCD evolution equation to consistently explain several intriguing phenomena associated with gamma-ray bursts. They are the GeV-TeV spectra of GRB 221009A, the remarkably symmetrical explosion cloud in kilonova AT2017gfo, and the absence of a very high-energy gamma-ray signature in GRB 170817A. We find that these occurrences can be attributed to the gluon condensation within nucleons, i.e., a significant number of soft gluons within nucleons are condensed at a critical momentum, resulting in the emergence of a steep and high peak in the gluon distributions. Through this profound connection between microscopic and macroscopic phenomena, we have not only expanded the applications of the hadronic scenario in cosmic gamma-ray emissions but also presented new evidence for the existence of gluon condensation.

Many dozens of protoplanetary discs show signatures of sculpting by planets. To help find these protoplanets by direct imaging, we compute their broadband spectral energy distributions, which overlap with the JWST (James Webb Space Telescope) and ALMA (Atacama Large Millimeter Array) passbands. We consider circumplanetary spherical envelopes, and combinations of circumplanetary discs and envelopes, heated by protoplanet accretion and external irradiation. Searches with JWST's NIRCam (Near-Infrared Camera) and the blue portion of MIRI (Mid-Infrared Instrument) are most promising since $\sim$300--1000 K protoplanets outshine their $\sim$20--50 K circumstellar environs at wavelengths of $\sim$2--20 $\mu$m. Detection is easier if circumplanetary dust settles into discs (more likely for more massive planets) or is less abundant per unit mass gas (because of grain growth or filtration). At wavelengths longer than 20 $\mu$m, circumplanetary material is difficult to see against the circumstellar disc's surface layers that directly absorb starlight and reprocess it to the far-infrared. Such contaminating circumstellar emission can be serious even within the evacuated gaps observed by ALMA. Only in strongly depleted regions, like the cavity of the transitional disc PDS 70 where two protoplanets have been confirmed, may long-wavelength windows open for protoplanet study. We compile a list of candidate protoplanets and identify those with potentially the highest accretion luminosities, all peaking in the near-infrared.

Alma Maria Sebastian, Emma Ryan-Weber, Rebecca L.Davies, George D.Becker, Laura C. Keating, Valentina D'Odorico, Romain A.Meyer, Sarah E.I. Bosman, Guido Cupani, Girish Kulkarni, Martin G. Haehnelt, Samuel Lai, Anna-Christina Eilers, Manuela Bischetti, Simona Gallerani

20 pages, 21 figures

Intervening metal absorbers in quasar spectra at $z > 6$ can be used as probes to study the chemical enrichment of the Universe during the Epoch of Reionization (EoR). This work presents the comoving line densities ($dn/dX$) of low ionisation absorbers, namely, Mg II (2796\r{A}), C II (1334\r{A}) and O I (1302\r{A}) across $2 <z < 6$ using the E-XQR-30 metal absorber catalog prepared from 42 XSHOOTER quasar spectra at $5.8 < z < 6.6$. Here, we analyse 280 Mg II ($1.9 < z < 6.4$), 22 C II ($5.2 < z < 6.4$) and 10 O I ($5.3 < z < 6.4$) intervening absorbers, thereby building up on previous studies with improved sensitivity of 50% completeness at an equivalent width of $W > 0.03$\r{A}. For the first time, we present the comoving line densities of 131 weak ($W < 0.3$\r{A}) intervening Mg II absorbers at $1.9 < z < 6.4$ which exhibit constant evolution with redshift similar to medium ($0.3 < W < 1.0$\r{A}) absorbers. However, the cosmic mass density of Mg II - dominated by strong Mg II systems - traces the evolution of global star formation history from redshift 1.9 to 5.5. E-XQR-30 also increases the absorption path length by a factor of 50% for C II and O I whose line densities show a rising trend towards $z > 5$, in agreement with previous works. In the context of a decline in metal enrichment of the Universe at $z > 5$, the overall evolution in the incidence rates of absorption systems can be explained by a weak - possibly soft fluctuating - UV background. Our results, thereby, provide evidence for a late reionization continuing to occur in metal-enriched and therefore, biased regions in the Universe.

Y. Omiya, K. Nakazawa, T. Tamura, H. Akamatsu, K. Matsushita, N. Okabe, K. Sato, Y. Fujita, L. Gu, A. Simionescu, Y. Ichinohe, C.J. Riseley, T. Akahori, D. Ito, K. Sakai, K. Kurahara

Abell 3667 is a near-by (z=0.056) merging cluster with the most prominent cold front and a pair of two bright radio relics. Assuming a face-to-face merger scenario, the origin of the cold front is often considered to be a remnant core of the cluster stripped of its surrounding ICM. Sarazin et al. (2016) proposed an offset merger scenario in which the sub-cluster cores rotate after the first core-crossing. To distinguish between these scenarios, we revisited the ICM distribution and measured the line-of-sight bulk ICM velocity using a calibration technique proposed by Sanders et al. (2020). In the unsharp masked image, we identified several ICM features, some of which we detected for the first time. There is an enhancement of the X-ray surface brightness extending from the 1st BCG to the cold front, which is named the "BCG-E tail". The notable feature is the "RG1 vortex", which is a clockwise vortex-like enhancement with a radius of about 250 kpc connecting the 1st BCG to the radio galaxy (RG1). It is particularly enhanced near the north of the 1st BCG, which is named the "BCG-N tail". The thermodynamic map shows that the ICM in the RG1 vortex has a relatively high abundance of 0.5-0.6 solar compared to the surrounding regions. The ICM of the BCG-E tail also has high abundance, and low pseudo-entropy, and can be interpreted as the remnant of the ICM of the cluster core. Including its arc-like shape, the RG1 vortex supports the idea that the ICM around the cluster center is rotating, which is natural in an offset merger scenario. The results of the line-of-sight bulk ICM velocity measurement show that the ICM around the BCG-N tail is redshifted with a velocity difference of 940$\pm$440 km s$^{-1}$ compared to the optical redshift of the 1st BCG. Other symptoms of diversity in the line-of-sight velocity of the ICM were also obtained and discussed in the context of the offset merger.

Francisco Oyarzún, Vincent Chambouleyron, Benoit Neichel, Thierry Fusco, Andrés Guesalaga

The use of artificial Laser Guide Stars (LGS) is planned for the new generation of giant segmented mirror telescopes, to extend the sky coverage of their adaptive optics systems. The LGS, being a 3D object at a finite distance will have a large elongation that will affect its use with the Shack-Hartmann (SH) wavefront sensor. In this paper, we compute the expected performance for a Pyramid WaveFront Sensor (PWFS) using a LGS for a 40 m telescope affected by photon noise, and also extend the analysis to a flat 2D object as reference. We developed a new way to discretize the LGS, and a new, faster method of propagating the light for any Fourier Filtering wavefront sensors (FFWFS) when using extended objects. We present the use of a sensitivity model to predict the performance of a closed-loop adaptive optic system. We optimized a point source calibrated interaction matrix to accommodate the signal of an extended object, by means of computing optical gains using a convolutional model. We found that the sensitivity drop, given the size of the extended laser source, is large enough to make the system operate in a low-performance regime given the expected return flux of the LGS. The width of the laser beam, rather than the thickness of the sodium layer was identified as the limiting factor. Even an ideal, flat LGS will have a drop in performance due to the flux of the LGS, and small variations in the return flux will result in large variations in performance. We conclude that knife-edge-like wavefront sensors, such as the PWFS, are not recommended for their use with LGS for a 40 m telescope, as they will operate in a low-performance regime, given the size of the extended object.

Peng Jia, Chao Lv, Yushan Li, Yongyang Sun, Shu Niu, Zhuoxiao Wang

Accepted by the AJ, comments are welcome. The complete code could be downloaded from: DOI: 10.12149/101387

In recent years, there has been a gradual increase in the performance of Complementary Metal Oxide Semiconductor (CMOS) cameras. These cameras have gained popularity as a viable alternative to charge-coupled device (CCD) cameras in a wide range of applications. One particular application is the CMOS camera installed in small space telescopes. However, the limited power and spatial resources available on satellites present challenges in maintaining ideal observation conditions, including temperature and radiation environment. Consequently, images captured by CMOS cameras are susceptible to issues such as dark current noise and defective pixels. In this paper, we introduce a data-driven framework for mitigating dark current noise and bad pixels for CMOS cameras. Our approach involves two key steps: pixel clustering and function fitting. During pixel clustering step, we identify and group pixels exhibiting similar dark current noise properties. Subsequently, in the function fitting step, we formulate functions that capture the relationship between dark current and temperature, as dictated by the Arrhenius law. Our framework leverages ground-based test data to establish distinct temperature-dark current relations for pixels within different clusters. The cluster results could then be utilized to estimate the dark current noise level and detect bad pixels from real observational data. To assess the effectiveness of our approach, we have conducted tests using real observation data obtained from the Yangwang-1 satellite, equipped with a near-ultraviolet telescope and an optical telescope. The results show a considerable improvement in the detection efficiency of space-based telescopes.

Qinchun Ma, Yuhan Wen, Xue-Bing Wu, Huapeng Gu, Yuming Fu

23 pages, 34 figures, accepted for publication in the Astrophysical Journal

In our previous work on broadband photometric reverberation mapping (PRM), we proposed the ICCF-Cut process to obtain the time lags of H$\alpha$ emission line from two broadband lightcurves via subtracting the continuum emission from the line band. Extending the work, we enlarge our sample to the Zwicky Transient Facility (ZTF) database. We adopt two criteria to select 123 type 1 AGNs with sufficient variability and smooth lightcurves from 3537 AGNs at $z<0.09$ with more than 100 epoch observations in the $g$ and $r$ bands from the ZTF database. We calculate the H$\alpha$ time lags for 23 of them which have previous spectroscopic reverberation mapping (SRM) results using ICCF-Cut, JAVELIN and $\chi ^2$ methods. Our obtained H$\alpha$ time lags are slightly larger than the H$\beta$ time lags, which is consistent with the previous SRM results and the theoretical model of the AGN broad line region. The comparisons between SRM and PRM lag distributions and between the subtracted emission line lightcurves indicate that after selecting AGNs with the two criteria, combining the ICCF-Cut, JAVELIN and $\chi^2$ methods provides an efficient way to get the reliable H$\alpha$ lags from the broadband PRM. Such techniques can be used to estimate the black hole masses of a large sample of AGNs in the large multi-epoch photometric sky surveys such as the Legacy Survey of Space and Time (LSST) and the survey from the Wide Field Survey Telescope (WFST) in the near future.

Swayamtrupta Panda, Szymon Kozłowski, Mariusz Gromadzki, Marcin Wrona, Patryk Iwanek, Andrzej Udalski, Michał K. Szymański, Igor Soszyński, Paweł Pietrukowicz, Krzysztof Ulaczyk, Jan Skowron, Radosław Poleski, Przemek Mróz, Dorota M. Skowron, Krzysztof Rybicki, Mateusz Mróz

15 pages, 8 figures, 8 tables, accepted for publication in The Astrophysical Journal Supplement Series

We use the spectroscopic data collected by the Magellanic Quasars Survey (MQS) as well as the photometric V- and I-band data from the Optical Gravitational Lensing Experiment (OGLE) to measure the physical parameters for active galactic nuclei (AGNs) located behind the Magellanic Clouds. The flux-uncalibrated MQS spectra were obtained with the 4-m Anglo-Australian Telescope and the AAOmega spectroscope (R=1300) in a typical ~1.5-hour visit. They span a spectral range of 3700-8500 Angstroms and have S/N ratios in a range of 3-300. We report the discovery and observational properties of 161 AGNs in this footprint, which expands the total number of spectroscopically confirmed AGNs by MQS to 919. After converting the OGLE mean magnitudes to the monochromatic luminosities at 5100 Angstroms, 3000 Angstroms, and 1350 Angstroms, we reliably measured the black hole masses for 165 out of 919 AGNs. The remaining physical parameters we provide are the bolometric luminosities and the Eddington ratios. A fraction of these AGNs have been observed by the OGLE survey since 1997 (all of them since 2001), enabling studies of correlations between their variability and physical parameters.

Marco Castellano, Lorenzo Napolitano, Adriano Fontana, Guido Roberts-Borsani, Tommaso Treu, Eros Vanzella, Jorge A. Zavala, Pablo Arrabal Haro, Antonello Calabrò, Mario Llerena, Sara Mascia, Emiliano Merlin, Diego Paris, Laura Pentericci, Paola Santini, Tom J. L. C. Bakx, Pietro Bergamini, Guido Cupani, Mark Dickinson, Alexei V. Filippenko, Karl Glazebrook, Claudio Grillo, Patrick L. Kelly, Matthew A. Malkan, Charlotte A. Mason, Takahiro Morishita, Themiya Nanayakkara, Piero Rosati, Eleonora Sani, Xin Wang, Ilsang Yoon

20 pages, 7 figures, 2 tables; Submitted to The Astrophysical Journal

We spectroscopically confirm the $M_{\rm UV} = -20.5$ mag galaxy GHZ2/GLASS-z12 to be at redshift $z=12.34$. The source was selected via NIRCam photometry in GLASS-JWST Early Release Science data, providing the first evidence of a surprising abundance of bright galaxies at $z \gtrsim 10$. The NIRSpec PRISM spectrum is remarkable and unlike any local analog. It shows significant detections of N IV, C IV, He II, O III, C III, O II, and Ne III lines, and the first detection in a high-redshift object of the O III Bowen fluorescence line at 3133 {\AA} rest-frame. The prominent C IV line with rest-frame equivalent width (EW) $\sim 46$ {\AA} puts GHZ2 in the category of extreme C IV emitters characterised by hard radiation fields. GHZ2 displays UV lines with EWs that are only found in active galactic nuclei (AGNs) or composite objects at low/intermediate redshifts, and UV line-intensity ratios that are compatible both with AGNs and star formation in a low-metallicity environment. The nondetection of the very high-ionization lines [Ne IV] and [Ne V], and the remarkable similarity between GHZ2 and other known C IV emitters, favors a scenario in which the high ionizing output is due to very low metallicity, massive stars forming in a dense environment. We estimate a metallicity $\lesssim 0.1 Z/{\rm Z}_{\odot}$, a high ionization parameter logU $>$ -2, a N/O abundance 4--5 times the solar value, and a subsolar C/O ratio similar to the recently discovered class of nitrogen-enhanced objects at high redshift. Considering its abundance patterns and the high stellar mass density ($10^4$ M$_{\odot}$ pc$^{-2}$), GHZ2 is an ideal formation site for the progenitors of today's globular clusters. The remarkable brightness of GHZ2 makes it a "Rosetta stone" for understanding the physics of galaxy formation within just 360 Myr after the Big Bang.

F. Tinaut-Ruano, J. de León, E. Tatsumi, D. Morate, M. Mahlke, P. Tanga, J. Licandro

14 pages 13 figures, accepted for publication in A&A

In the context of charge-coupled devices (CCDs), the ultraviolet (UV) region has mostly remained unexplored after the 1990s. Gaia DR3 offers the community a unique opportunity to explore tens of thousands of asteroids in the near-UV as a proxy of the UV absorption. This absorption has been proposed in previous works as a diagnostic of hydration, organics, and space weathering. Aims. In this work, we aim to explore the potential of the NUV as a diagnostic region for primitive asteroids using Gaia DR3. We used a corrective factor over the blue part of Gaia spectra to erase the solar analog selection effect. We identified an artificial relation between the band noise and slope and applied a signal-to-noise ratio (S/N) threshold for Gaia bands. Meeting the quality standards, we employed a Markov chain Monte Carlo (MCMC) algorithm to compute the albedo threshold, maximizing primitive asteroid inclusion. Utilizing one- and two-dimensional (1D and 2D) projections, along with dimensionality-reduction methods (such as PCA and UMAP), we identified primitive asteroid populations. We uncovered: (a) the first observational evidence linking UV absorption to the 0.7 {\mu}m band, tied to hydrated iron-rich phyllosilicates; and (b) a 2D space revealing a split in C-type asteroids based on spectral features, including UV absorption. The computed average depth (3.5 +- 1.0 %) and center (0.70 +- 0.03 {\mu}m) of the 0.7 {\mu}m absorption band for primitive asteroids observed with Gaia is in agreement with the literature values. In this paper, we shed light on the importance of the UV absorption feature to discriminate among different mineralogies (i.e., iron-rich phyllosilicates vs. iron-poor) or to identify taxonomies that are conflated in the visible (i.e., F-types vs. B-types). We have shown that this is a promising region for diagnostic studies of the composition of primitive asteroids.

Yingjie Zhu, Shadia R. Habbal, Adalbert Ding, Bryan Yamashiro, Enrico Landi, Benjamin Boe, Sage Constantinou, Michael Nassir

29 pages, 19 figures, accepted for publication in The Astrophysical Journal. Jupyter notebooks and IDL codes are available at this https URL and this https URL

The spectroscopic observations presented here were acquired during the 2017 August 21 total solar eclipse with a three-channel partially multiplexed imaging spectrometer (3PAMIS) operating at extremely high orders ($>$ 50). The 4 $R_\odot$ extent of the slit in the North-South direction scanned the corona starting from the central meridian out to approximately 1.0 $R_\odot$ off the east limb throughout totality. The line widths and Doppler shifts of the Fe X (637.4 nm) and Fe XIV (530.3 nm) emission lines, characteristic of $1.1 \times 10^6$ K and $1.8 \times 10^6$ K electron temperatures respectively, varied across the different coronal structures intercepted by the slit. Fe XIV was the dominant emission in the closed fields of an active region and the base of a streamer, with relatively constant 20 - 30 km s$^{-1}$ line widths independent of the height. In contrast, Fe X emission exhibited broader ($>40 $km s$^{-1}$) line widths in open fields which increased with height, in particular in the polar coronal hole. Inferences of line widths and Doppler shifts were consistent with extreme ultraviolet (EUV) observations from Hinode/EIS, as well as with the near-infrared Fe XIII 1074 nm line observed by CoMP. The differences in the spectral line widths between distinct coronal structures are interpreted as an indication of the predominance of wave heating in open structures versus localized heating in closed structures. This study underscores the unparalleled advantages and the enormous potential of TSE spectroscopy in measuring line widths simultaneously in open and closed fields at high altitudes, with minimal exposure times, stray light levels, and instrumental widths.

Some short gamma-ray bursts (SGRBs) exhibit a short duration and spectral hard emission (referred to as a "hard spike") followed by a slightly longer soft emission (known as a "soft tail"). We identified nine SGRBs with known redshift in the \textit{Swift}/BAT gamma-ray burst catalog by specifically searching for the soft tail. We found that spectra of these SGRBs can be described as a cutoff power-law model for both hard spike and soft tail, and both show time variation keeping the $E_{\rm peak}$--$L_{\rm iso}$ correlation. This suggests that the emission mechanism of both phenomena is identical. Furthermore, we found a trend of luminosity evolution as a function of redshift. This phenomenon suggests that these bursts originate from sources that have intrinsically bright and/or energy density concentrated within a narrower jet at higher redshift. We demonstrate that the average jet opening angle, derived from the jet break, can be explained by considering a model based on a strongly redshift-dependent jet opening angle.

The most recent SH0ES measurement of the Hubble constant, based on type Ia supernovae from the Pantheon+ compilation, employs corrections of supernova peak magnitudes which effectively accounts for extinction in the supernova host galaxies. These corrections are estimated using a probabilistic model which is trained on Hubble flow (z>0.03) supernovae and extrapolated to the calibration galaxies (those with observed Cepheid distances). By comparing the corrected peak magnitudes to distance moduli from Cepheids, we show that this standard approach underestimates the brightness of reddened supernovae in the high stellar-mass ($M_{\star}>10^{10}M_{\odot}$) calibration galaxies. This can be traced back to the fact that for these galaxies, a low total-to-selective extinction coefficient (R_B~3) is assumed, while for the low stellar-mass analogues a more standard R_B~4 is assumed. We propose a minimalistic modification of the Pantheon+ extinction model in order to alleviate this systematic effect. The modification is twofold and it involves: (i) the same, Milky Way-like distribution of R_B in all calibration galaxies (with mean R_B of 4.3 -- consistent with the extinction curve used for colour corrections of the Cepheids -- and scatter 0.4) and (ii) a modified shape of the E(B-V) reddening distribution while keeping the same effective slope of the supernova peak magnitude-colour relation and the same mean E(B-V) reddening as measured for supernovae in the Hubble flow. We show that this new approach yields a significantly better fit ($\Delta$BIC=-11) to the calibration data and results in a lower value of the derived Hubble constant through a stronger extinction correction of supernovae in the calibration galaxies. Our result is $H_{0}=70.5\pm1$ km/s/Mpc implying a reduction of the tension with the Planck $H_{0}$ measurement assuming a flat LCDM cosmology from $5.2\sigma$ to $2.8\sigma$.

Kartika Sangal, A.K. Srivastava, P. Kayshap, Ding Yuan, E. Scullion

18 pages, 10 figures, accepted for publication in ApJ

In the current analysis, we use spectroscopic observations of the quiet-Sun made by IRIS instrument, and investigate wave propagation. We analyze various spectral lines formed in different atmospheric layers such as the photosphere, chromosphere, and transition region. We examine Doppler velocity time-series at various locations in the quiet-Sun to determine the dominant oscillation periods. Our results executing statistical analysis resemble those of the classical physical scenario, indicating that the photosphere is mainly characterized by the dominant 5-minute period, while the chromosphere is primarily associated with the 3-minute oscillation period. In the transition region, we observe a variety of oscillation periods, with dominant periods of 3, 8, and 12 minutes. We estimate the cut-off frequency by deducing phase difference between two Doppler velocity time-series obtained from spectral line pairs in different atmospheric layers formed at different temperatures. It reveals a significant correlation between 3-minute periods in TR and photospheric oscillations, suggesting that these oscillations in the TR might propagate from the photosphere. Additionally, we analyze the phase difference between chromospheric oscillations and photospheric oscillations, demonstrating that only the 3-minute oscillations propagate upwards. Based on the statistical analyses, we suggest the presence of magnetoacoustic waves in the solar atmosphere in which some are propagating from the lower solar atmosphere upward, while some others are propagating downward. TR carries both long-period oscillations generated in situ, and some photospheric oscillations which are also able to reach there from below.

A. Reguitti, G. Pignata, A. Pastorello, R. Dastidar, D. E. Reichart, J. B. Haislip, V. V. Kouprianov

14 pages, 13 figures, 2 tables, accepted for publication on A&A

We conducted a search for luminous outbursts prior to the explosion of Type IIn Supernovae (SNe IIn). We built a sample of 27 objects spectroscopically classified as SNe IIn, all located at $z<0.015$. Using deep archival SN fields images taken up to nearly 20 years prior from transient surveys (PTF, ZTF, DES, CHASE) and major astronomical observatories (ESO and NOAO), we found at least one outburst years to months before the explosion of seven SNe IIn, the earliest precursor being 10 years prior to the explosion of SN 2019bxq. The maximum absolute magnitudes of the outbursts range between -11.5 mag and -15 mag, and the eruptive phases last for a few weeks to a few years. The $g-r$ colour measured for three objects during their outburst is relatively red, with $g-r$ ranging between 0.5 and 1.0 mag. This is similar to the colour expected during the eruptions of Luminous Blue Variables. We noticed that the SNe with pre-SN outbursts have light curves with faster decline rates than those that do not show pre-SN outbursts. SN 2011fh is remarkable, as it is still visible 12 years after the luminous SN-like event, indicating that the progenitor possibly survived, or that the interaction is still on-going. We detect precursor activity in 29% of bona-fide SNe~IIn in our sample. However, a quantitative assessment of the observational biases affecting the sample suggests that this fraction underestimates the intrinsic precursor occurrence rate.

Naoki Ogino, Makoto Arimoto, Tatsuya Sawano, Daisuke Yonetoku, Hsien-chieh Shen, Takanori Sakamoto, Junko S. Hiraga, Yoichi Yatsu, Tatehiro Mihara

We developed an FPGA-based high-speed readout system for a complementary metal-oxide-semiconductor (CMOS) image sensor to observe soft X-ray transients in future satellite missions, such as HiZ-GUNDAM. Our previous research revealed that the CMOS image sensor has low-energy X-ray detection capability (0.4-4 keV) and strong radiation tolerance, which satisfies the requirements of the HiZ-GUNDAM mission. However, CMOS sensors typically have small pixel sizes (e.g., $\sim$10 ${\rm \mu m}$), resulting in large volumes of image data. GSENSE400BSI has 2048$\times$2048 pixels, producing 6 Mbyte per frame. These large volumes of observed raw image data cannot be stored in a satellite bus system with a limited storage size. Therefore, only X-ray photon events must be extracted from the raw image data. Furthermore, the readout time of CMOS image sensors is approximately ten times faster than that of typical X-ray CCDs, requiring faster event extraction on a timescale of $\sim$0.1 s. To address these issues, we have developed an FPGA-based image signal processing system capable of high-speed X-ray event extraction onboard without storing raw image data. The developed compact system enabled mounting on a CubeSat mission, facilitating early in-orbit operation demonstration. Here, we present the design and results of the performance evaluation tests of the proposed FPGA-based readout system. Utilizing X-ray irradiation experiments, the results of the X-ray event extraction with the onboard and offline processing methods were consistent, validating the functionality of the proposed system.

Theoretical arguments as well as observations of young stellar objects (YSO) support the presence of a diversified circumstellar environment. A stellar jet is thought to account for most of the stellar spin down and disk wind outflow for the observed high mass loss rate, thus playing a major role in the launching of powerful jets. RY Tau, for instance, is an extensively studied intermediate mass pre-main sequence star. Observational data reveal a small scale jet called microjet. Nevertheless, it is not clear how the microjet shapes the jet observed at a large scale. The goal is to investigate the spatial stability and structure of the central jet at a large scale by mixing the stellar and disk components. We mix two existing analytical self-similar models for the disk and the stellar winds to build the initial set-ups. Instead of using a polytropic equation of state, we map from the analytical solutions, the heating and cooling sources. The heating exchange rate is controlled by two parameters, its spatial extent and its intensity. The central jet and the surrounding disk are strongly affected by these two parameters. We separate the results in three categories, which show different emissivity, temperature, and velocity maps. We reached this categorization by looking at the opening angle of the stellar solution. For cylindrically, well collimated jets, we have opening angles as low as 10 degrees between 8 and 10 au, and for the wider jets, we can reach 30 degrees with a morphology closer to radial solar winds. Our parametric study shows that the less heated the outflow is, the more collimated it appears. We also show that recollimation shocks appear consistently with UV observations in terms of temperature but not density.

Roberto Bruno (1), Rossana DeMarco (1), Raffaella D Amicis (1), Denise Perrone (2), Maria Federica Marcucci (1), Daniele Telloni (3), Raffaele Marino (4), Luca Sorriso Valvo (5), Vito Fortunato (6), Gennaro Mele (7), Francesco Monti (8), Andrei Fedorov (9), Philippe Louarn (9), Chris Owen (10), Stefano Livi (11) ((1) INAF-IAPS, (2) ASI, (3) INAF-OATO, (4) Univ Lyon, CNRS, (5) CNR-ISTP, (6) Planetek, (8) Leonardo, (9) TSD, (10) UCL-MSSL, (11) SwRI)

paper submitted to APJ, 20 pages, 17 figures

The problems of heating and acceleration of solar wind particles are of significant and enduring interest in astrophysics. The interactions between waves and particles are crucial in determining the distributions of proton and alpha particles, resulting in non-Maxwellian characteristics including temperature anisotropies and particle beams. These processes can be better understood as long as the beam can be separated from the core for the two major components of the solar wind. We utilized an alternative numerical approach that leverages the clustering technique employed in Machine Learning to differentiate the primary populations within the velocity distribution, rather than employing the conventional biMaxwellian fitting method. Separation of the core and beam revealed new features for protons and alphas. We estimated that the total temperature of the two beams was slightly higher than that of their respective cores, and the temperature anisotropy for the cores and beams was larger than 1. We concluded that the temperature ratio between alphas and protons largely over 4 is due to the presence of a massive alpha beam, which is approximately 50% of the alpha core. We provided evidence that the alpha core and beam populations are sensitive to Alfv\'enic fluctuations and the surfing effect found in the literature can be recovered only when considering the core and beam as a single population. Several similarities between proton and alpha beams would suggest a common and local generation mechanism not shared with the alpha core, which may not have necessarily been accelerated and heated locally.

We revisit the claimed detection of a new cosmic microwave background (CMB) foreground based on the correlation between low-redshift 2MASS Redshift Survey (2MRS) galaxies and CMB temperature maps from the Planck and WMAP missions. We reproduce the reported measurements but argue that the original analysis significantly underestimated the uncertainties. We cross-correlate the 2MRS galaxy positions with simulated CMB maps and show that the correlation measured with the real data for late-type spiral galaxies at angular scales $\theta\geq0.1^{\circ}$ and redshift $cz<4500$ km s$^{-1}$ is consistent with zero at the $1.7\sigma$ level or less, depending on the exact CMB map and simulation construction. This was the sample that formed the basis for the original detection claim. For smaller angular separations the results are not robust to galaxy type or CMB cleaning method, and we are unable to draw firm conclusions. The original analysis did not propose a specific, falsifiable physical correlation mechanism, and it is impossible to rule out any contribution from an underlying physical effect. However, given our calculations, the lack of signal from expanding the redshift range, and the lack of corroboration from other galaxy surveys, we do not find the evidence for a new CMB foreground signal compelling.

Jorge A. Zavala, Marco Castellano, Hollis B. Akins, Tom J. L. C. Bakx, Denis Burgarella, Caitlin M. Casey, Óscar A. Chávez Ortiz, Mark Dickinson, Steven L. Finkelstein, Ikki Mitsuhashi, Kimihiko Nakajima, Pablo G. Pérez-González, Pablo Arrabal Haro, Veronique Buat, Bren Backhaus, Antonello Calabrò, Nikko J. Cleri, David Fernández-Arenas, Adriano Fontana, Maximilien Franco, Mauro Giavalisco, Norman A. Grogin, Nimish Hathi, Michaela Hirschmann, Ryota Ikeda, Intae Jung, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Rebeca L. Larson, Jed McKinney, Casey Papovich, Toshiki Saito, Paola Santini, Roberto Terlevich, Elena Terlevich, Tommaso Treu, L. Y. Aaron Yung

14 pages, 6 figures. Submitted for publication. Comments are most welcome

The James Webb Space Telescope (JWST) has discovered a surprising population of bright galaxies in the very early universe (< 500 Myrs after the Big Bang) that is hard to explain with conventional galaxy formation models and whose physical properties remain to be fully understood. Insight into the internal physics of galaxies is captured best via observations of excited-state atomic transitions of ionized gas, but beyond z~7-9, the brightest spectral signatures are redshifted into the mid-infrared regime, where observations are increasingly more difficult. Here, we present the first detection of a hydrogen recombination line (H{\alpha}) and doubly-ionized oxygen ([OIII]4959,5007{\AA}) at z>10 using the JWST Mid-Infrared Instrument, MIRI. These detections place the bright galaxy GHZ2/GLASS-z12 at z=12.33+/-0.02, making it the most distant astronomical object with direct spectroscopic detection of these lines and the brightest confirmed object at this epoch. These observations provide key insights into the conditions of this primeval galaxy, which shows hard ionizing conditions rarely seen in the local Universe and likely driven by compact, young (<30 Myr) star formation. Its oxygen-to-hydrogen abundance is close to a tenth of the solar value, indicating a rapid metal enrichment during the earliest phases of galaxy formation. This study confirms the unique conditions of the brightest and most distant galaxies recently discovered by JWST and the huge potential of mid-IR observations to characterize these systems, opening a range of new possibilities in the study of the very early Universe.

Pedro Labraña, Juan Ortiz

22 pages, 4 figures, changes in abstract and introduction, improvements in the text and comments added in conclusions. References added

The study of Emergent Universe models is based on the assumption that the universe emerged from a past eternal Einstein Static (ES) state towards an inflationary phase and then evolves into a hot big bang era. These models are appealing since they provide specific examples of nonsingular (geodesically complete) inflationary universes. However, it has been pointed out by Mithani-Vilenkin that certain Emergent Universe scenarios which have a classically stable ES state could experience a semiclassical instability and collapse. In this paper, we investigate the classical and quantum stability of the ES regime of Emergent Universes within the framework of Jordan-Brans-Dicke theory. We demonstrate that when considering these models, it is possible to have both, classical and semiclassical stability of the ES state without addressing the instability highlighted by Mithani-Vilenkin.

N.Galikyan, Sh.Khlghatyan, A.A.Kocharyan, V.G.Gurzadyan

6 pages, 3 figs; Eur Phys J Plus (in press)

Physics-informed neural network (PINN) analysis of the dynamics of S-stars in the vicinity of the supermassive black hole in the Galactic center is performed within General Relativity treatment. The aim is to reveal the role of possible extended mass (dark matter) configuration in the dynamics of the S-stars, in addition to the dominating central black hole's mass. The PINN training fails to detect the extended mass perturbation in the observational data for S2 star within the existing data accuracy, and the precession constraint indicates no signature of extended mass up to 0.01% of the central mass inside the apocenter of S2. Neural networks analysis thus confirm its efficiency in the analysis of the S-star dynamics.

Recently, the Event Horizon Telescope (EHT) achieved the realization of an image of the supermassive black hole $\textrm{Sgr~A}^\star$ showing an angular shadow diameter $\mathcal{D}= 48.7 \pm 7\mu as$ and the fractional deviation $\mathbf{\delta} = -0.08^{+0.09}_{-0.09}~\text{(VLTI)},-0.04^{+0.09}_{-0.10}~\text{(Keck)}$, alongside the earlier image of $\textrm{M87}^\star$ with angular diameter $ \mathcal{D}=42 \pm 3 \mu as$, deviation $\mathbf{\delta}=-0.01^{+0.17}_{-0.17}$ and deviations from circularity estimated to be $\Delta \mathcal{C}\lesssim 10\%$. In addition, the shadow radii are assessed within the ranges $3.38 \le \frac{r_{\text{s}}}{M} \le 6.91$ for $\textrm{M87}^\star$ and $3.85 \le \frac{r_{\text{s}}}{M} \le 5.72$ as well as $3.95 \le \frac{r_{\text{s}}}{M} \le 5.92$ for $\textrm{Sgr~A}^\star$ using the Very Large Telescope Interferometer (VLTI) and Keck observatories, respectively. These values are provided with $1$-$\sigma$ and $2$-$\sigma$ measurements. Such realizations can unveil a better comprehension of gravitational physics at the horizon scale. In this paper, we use the EHT observational results for $\textrm{M87}^\star$ and $\textrm{Sgr~A}^\star$ to elaborate the constraints on parameters of accelerating black holes with a cosmological constant. Concretely, we utilize the mass and distance of both black holes to derive the observables associated with the accelerating black hole shadow. First, we compare our findings with observed quantities such as angular diameter, circularity, shadow radius, and the fractional deviation from the $\textrm{M87}^\star$ data. This comparison reveals constraints within the acceleration parameter and the cosmological constant... Lastly, one cannot rule out the possibility of the negative values for the cosmological constant on the emergence of accelerated black hole solutions within the context of minimal gauged supergravity...

Sophia Yi, Adrien Kuntz, Enrico Barausse, Emanuele Berti, Mark Ho-Yeuk Cheung, Konstantinos Kritos, Andrea Maselli

20 pages, 15 figures, 2 tables

In the aftermath of a binary black hole merger event, the gravitational wave signal emitted by the remnant black hole is modeled as a superposition of damped sinusoids known as quasinormal modes. While the dominant quasinormal modes originating from linear black hole perturbation theory have been studied extensively in this post-merger "ringdown" phase, more accurate models of ringdown radiation include the nonlinear modes arising from higher-order perturbations of the remnant black hole spacetime. We explore the detectability of quadratic quasinormal modes with both ground- and space-based next-generation detectors. We quantify how predictions of the quadratic mode detectability depend on the quasinormal mode starting times. We then calculate the signal-to-noise ratio of quadratic modes for several detectors and binary black hole populations, focusing on the ($220\times220$) mode - i.e., on the quadratic term sourced by the square of the linear $(220)$ mode. For the events with the loudest quadratic mode signal-to-noise ratios, we additionally compute statistical errors on the mode parameters in order to further ascertain the distinguishability of the quadratic mode from the linear quasinormal modes. The astrophysical models used in this paper suggest that while the quadratic mode may be detectable in at most a few events with ground-based detectors, the prospects for detection with the Laser Interferometer Space Antenna (LISA) are more optimistic.

In this work we study the power spectrum of the Stochastic Gravitational Wave Background produced by standard and biased domain wall networks, using the Velocity-dependent One-Scale model to compute the cosmological evolution of their characteristic scale and root-mean-squared velocity. We consider a standard radiation + $\Lambda \rm CDM$ background and assume that a constant fraction of the energy of collapsing domain walls is emitted in the form of gravitational waves. We show that, in an expanding background, the total energy density in gravitational radiation decreases with cosmic time (after a short initial period of quick growth). We also propose a two parameter model for the scale-dependence of the frequency distribution of the gravitational waves emitted by collapsing domain walls. We determine the corresponding power spectrum of the Stochastic Gravitational Wave Background generated by domain walls, showing that it is a monotonic decreasing function of the frequency for frequencies larger than that of the peak generated by the walls that have decayed most recently. We also develop an analytical approximation to this spectrum, assuming perfect linear scaling during both the radiation and matter eras, in order to characterize the dependence of the amplitude, peak frequency and slope of the power spectrum on the model parameters.

Sabrina Guastavino, Katsiaryna Bahamazava, Emma Perracchione, Fabiana Camattari, Gianluca Audone, Daniele Telloni, Roberto Susino, Gianalfredo Nicolini, Silvano Fineschi, Michele Piana, Anna Maria Massone

This study addresses the prediction of geomagnetic disturbances by exploiting machine learning techniques. Specifically, the Long-Short Term Memory recurrent neural network, which is particularly suited for application over long time series, is employed in the analysis of in-situ measurements of solar wind plasma and magnetic field acquired over more than one solar cycle, from $2005$ to $2019$, at the Lagrangian point L$1$. The problem is approached as a binary classification aiming to predict one hour in advance a decrease in the SYM-H geomagnetic activity index below the threshold of $-50$ nT, which is generally regarded as indicative of magnetospheric perturbations. The strong class imbalance issue is tackled by using an appropriate loss function tailored to optimize appropriate skill scores in the training phase of the neural network. Beside classical skill scores, value-weighted skill scores are then employed to evaluate predictions, suitable in the study of problems, such as the one faced here, characterized by strong temporal variability. For the first time, the content of magnetic helicity and energy carried by solar transients, associated with their detection and likelihood of geo-effectiveness, were considered as input features of the network architecture. Their predictive capabilities are demonstrated through a correlation-driven feature selection method to rank the most relevant characteristics involved in the neural network prediction model. The optimal performance of the adopted neural network in properly forecasting the onset of geomagnetic storms, which is a crucial point for giving real warnings in an operational setting, is finally showed.

Alnadhief H. A. Alfedeel, M. Koussour, N. Myrzakulov

Astronomy and Computing accepted version

In this paper, we investigate the cosmological implications and constraints of Weyl-type $f(Q, T)$ gravity. This theory introduces a coupling between the non-metricity $Q$ and the trace $T$ of the energy-momentum tensor, using the principles of proper Weyl geometry. In this geometry, the scalar non-metricity $Q$, which characterizes the deviations from Riemannian geometry, is expressed in its standard Weyl form $\nabla _{\mu }g_{\alpha \beta }=-w_{\mu }g_{\alpha \beta }$ and is determined by a vector field $w_{\mu }$. To study the implications of this theory, we propose a deceleration parameter with a single unknown parameter $\chi$, which we constrain by using the latest cosmological data. By solving the field equations derived from Weyl-type $f(Q, T)$ gravity, we aim to understand the behavior of the energy conditions within this framework. In the present work, we consider two well-motivated forms of the function $f(Q, T)$: (i) the linear model represented by $f(Q, T) = \alpha Q + \frac{\beta}{6\kappa^2} T$, and (ii) the coupling model represented by $f(Q, T) = \frac{\gamma}{6H_0^2 \kappa^2} QT$, where $\alpha$, $\beta$, and $\gamma$ are free parameters. Here, $\kappa^2 = \frac{1}{16\pi G}$ represents the gravitational coupling constant. In both of the models considered, the strong energy condition is violated, indicating consistency with the present accelerated expansion. However, the null, weak, and dominant energy conditions are satisfied in these models.

We study the linear stability of spontaneously scalarized black holes (BHs) induced by a scalar field $\phi$ coupled to a Gauss-Bonnet (GB) invariant $R_{\rm GB}^2$. For the scalar-GB coupling $\xi(\phi)=(\eta/8) (\phi^2+\alpha \phi^4)$, where $\eta$ and $\alpha$ are constants, we first show that there are no angular Laplacian instabilities of even-parity perturbations far away from the horizon for large multipoles $l \gg 1$. The deviation of angular propagation speeds from the speed of light is largest on the horizon, whose property can be used to put constraints on the model parameters. For $\alpha \gtrsim -1$, the region in which the scalarized BH is subject to angular Laplacian instabilities can emerge. Provided that $\alpha \lesssim -1$ and $-1/2<\alpha \phi_0^2<-0.1155$, where $\phi_0$ is the field value on the horizon with a unit of the reduced Planck mass $M_{\rm Pl}=1$, there are scalarized BH solutions satisfying all the linear stability conditions throughout the horizon exterior. We also study the stability of spontaneously scalarized BHs in scalar-GB theories with a nonminimal coupling $-\beta \phi^2 R/16$, where $\beta$ is a positive constant and $R$ is a Ricci scalar. As the amplitude of the field on the horizon approaches an upper limit $|\phi_0|=4/\sqrt{\beta}$, one of the squared angular propagation speeds $c_{\Omega-}^2$ enters the instability region $c_{\Omega-}^2<0$. So long as $|\phi_0|$ is smaller than a maximum value determined for each $\beta$ in the range $\beta>5$, however, the scalarized BHs are linearly stable in both angular and radial directions.

Scalar-tensor theories have taken on a key role in attempts to confront the growing open questions in standard cosmology. It is important to understand entirely their dynamics at perturbative level including any possible spatial dependence in their growth of large scale structures. In this work, we investigate the spatial dependence of the growth rate of scalar-tensor theories through the M\'{e}sz\'{a}ros equation. We confirm that at subhorizon level this dependence does not play a major role for viable models. However, we establish conditions on which this criterion is met which may be important for developing new models. In our work, we consider three specific models that exhibit spatial dependence of the growth rate at superhorizon modes, which may also be important for early Universe models.

Waleska P. F. de Medeiros, Matheus J. Lazo, Daniel Müller, Dinalva A. Sales

In this work, non perfect fluid, tilted source solutions in both Einstein-Hilbert General Relativity (GR) and Quadratic Gravity (QG) for the anisotropic Bianchi V model are addressed. Since the excellent CMBR match of Starobinsky's inflation with Planck's team measurements data, QG has acquired a prominent status in the effective sense, for sufficiently strong gravity fields. The main interest is in the numeric time evolution to the past towards the singularity and the behavior of the kinematic variables, vorticity, acceleration, and the expansion of this source substance. In QG we found that for universes with higher and smaller matter densities fall into the Kasner or isotropic singularity attractors to the past, respectively. We also found that the Kasner singularity attractor to the past has always zero vorticity, for both GR and QG theories. While for QG the isotropic singularity attractor may have divergent vorticity. For the set of assumptions and conditions supposed in this work, the isotropic singularity attractor, favors QG as compared to GR. Only in QG we were able to find a geometric singularity with divergences in all of the kinematic variables of the substance, decreasing to finite values to the future, upon time reversing. That is, we obtained an initial kinematic singularity substance, that approaches a perfect fluid source.

Praveen Kumar, Thomas Dent (IGFAE, Universidade de Santiago de Compostela)

22 pages, 6 figures, 2 tables

Maximizing the number of detections in matched filter searches for compact binary coalescence (CBC) gravitational wave (GW) signals requires a model of the source population distribution. In previous searches using the PyCBC framework, sensitivity to the population of binary black hole (BBH) mergers was improved by restricting the range of filter template mass ratios and use of a simple one-dimensional population model. However, this approach does not make use of our full knowledge of the population and cannot be extended to a full parameter space search. Here, we introduce a new ranking method, based on kernel density estimation (KDE) with adaptive bandwidth, to accurately model the probability distributions of binary source parameters over a template bank, both for signals and for noise events. We demonstrate this ranking method by conducting a search over LIGO-Virgo O3 data for BBH with unrestricted mass ratio, using a signal model derived from previous significant detected events. We achieve over 10% increase in sensitive volume for a simple power-law simulated signal population, compared to the previous BBH search. Correspondingly, with the new ranking, 8 additional candidate events above an inverse false alarm rate (IFAR) threshold 0.5 yr are identified.

Arhum Ansari, Pinak Banerjee, Prateksh Dhivakar, Sachin Jain, Nilay Kundu

40 pages + appendices; 2 figures

We study the conformal invariance of inflationary non-Gaussianities associated with scalar fluctuations in a non-Bunch-Davies initial state, known as the $\alpha$-vacuum, in single-field slow-roll inflation. The $\alpha$-vacuum is a one-parameter family of states, including the Bunch-Davies one, that preserves the conformal symmetry of inflationary dynamics in a nearly de-Sitter space-time. Working within the leading slow-roll approximation, we compute the four-point scalar correlator (the trispectrum) in $\alpha$-vacuum using the in-in formalism. We check that the conformal Ward identities are met between the three and four-point scalar $\alpha$-vacua correlators. Surprisingly, this contrasts the previously reported negative result of the Ward identities being violated between the two and the three-point correlators. We have also extended the wave-functional method, previously used for correlators with Bunch-Davies initial condition, to compute the three and four-point scalar correlators in $\alpha$-vacua. The results obtained from the wave-function method match the corresponding in-in results, adding further justification to our check of Ward identities with $\alpha$-vacua correlators.